Name: Author Spotlight: Investigating the Role of Repetitive DNA Misregulation in Cancer Initiation and Immunotherapy Resistance
Uploaded: 2024-12-13T14:00:00
Description: The CUT&amp;RUN technique facilitates detection of protein-DNA interactions across the genome. Typical applications of CUT&amp;RUN include profiling ...

All illustrated figures were created with BioRender.com. CAI acknowledges support provided through an Ovarian Cancer Research Alliance Early Career Investigator Award, a Forbeck Foundation Accelerator Grant, and the Minnestoa Ovarian Cancer Alliance National Early Detection Research Award.

NOTE: Information for CUT&#38;RUN fastq files in GSE126612 are available in Table 1. Information related to the software applications used in this study are listed in the Table of Materials.
1. Downloading Easy-Shells_CUTnRUN pipeline from its Github page
<ol>
	<li>Open terminal from operating system. 
	NOTE: If user is not sure how to open terminal in macOS and Windows, review this webpage (https://discovery.cs.illinois.edu/guides/Syst...

Beginner-Friendly Single-Language Pipeline for CUT&RUN Data Analysis

<ol>
	<li>Hainer, S. J., Fazzio, T. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=High-resolution+chromatin+profiling+using+CUT&amp;RUN.">High-resolution chromatin profiling using CUT&amp;RUN.</a> Curr Protoc Mol Biol. 126 (1), e85 (2019).</li><li>Zhang, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Model-based+analysis+of+ChiP-Seq+(MACS).">Model-based analysis of ChiP-Seq (MACS).</a> Genome Biology. 9 (9), R137 (2008).</li><li>Xu, S., Grullon, S., Ge, K., Peng, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Stem cell transcriptional networks: Methods and Protocols. , (2014).</li><li>Meers, M. P., Tenenbaum, D., Henikoff, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Peak+calling+by+sparse+enrichment+analysis+for+cut&amp;run+chromatin+profiling.">Peak calling by sparse enrichment analysis for cut&amp;run chromatin profiling.</a> Epigenetics Chromatin. 12 (1), 42 (2019).</li><li>Ashburner, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Gene+ontology:+Tool+for+the+unification+of+biology.+The+gene+ontology+consortium.">Gene ontology: Tool for the unification of biology. The gene ontology consortium.</a> Nat Genet. 25 (1), 25-29 (2000).</li><li>Harris, M. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+gene+ontology+(GO)+database+and+informatics+resource.">The gene ontology (GO) database and informatics resource.</a> Nucleic Acids Res. 32 (Database issue), D258-D261 (2004).</li><li>The Gene Ontology Consortium. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+gene+ontology+resource:+20+years+and+still+going+strong.">The gene ontology resource: 20 years and still going strong.</a> Nucleic Acids Res. 47 (D1), D330-D338 (2019).</li><li>Conesa, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Blast2go:+A+universal+tool+for+annotation,+visualization+and+analysis+in+functional+genomics+research.">Blast2go: A universal tool for annotation, visualization and analysis in functional genomics research.</a> Bioinformatics. 21 (18), 3674-3676 (2005).</li><li>Carbon, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=AmiGO:+Online+access+to+ontology+and+annotation+data.">AmiGO: Online access to ontology and annotation data.</a> Bioinformatics. 25 (2), 288-289 (2009).</li><li>Eden, E., Navon, R., Steinfeld, I., Lipson, D., Yakhini, Z. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Gorilla:+A+tool+for+discovery+and+visualization+of+enriched+go+terms+in+ranked+gene+lists.">Gorilla: A tool for discovery and visualization of enriched go terms in ranked gene lists.</a> BMC Bioinformatics. 10, 48 (2009).</li><li>Huang Da, W., Sherman, B. T., Lempicki, R. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Bioinformatics+enrichment+tools:+Paths+toward+the+comprehensive+functional+analysis+of+large+gene+lists.">Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists.</a> Nucleic Acids Res. 37 (1), 1-13 (2009).</li><li>Huang Da, W., Sherman, B. T., Lempicki, R. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Systematic+and+integrative+analysis+of+large+gene+lists+using+david+bioinformatics+resources.">Systematic and integrative analysis of large gene lists using david bioinformatics resources.</a> Nat Protoc. 4 (1), 44-57 (2009).</li><li>Ge, S. X., Jung, D., Yao, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=ShinyGO:+A+graphical+gene-set+enrichment+tool+for+animals+and+plants.">ShinyGO: A graphical gene-set enrichment tool for animals and plants.</a> Bioinformatics. 36 (8), 2628-2629 (2020).</li><li>Tang, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=SRplot:+A+free+online+platform+for+data+visualization+and+graphing.">SRplot: A free online platform for data visualization and graphing.</a> PLoS One. 18 (11), e0294236 (2023).</li><li>Ram&#237;rez, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Deeptools2:+A+next+generation+web+server+for+deep-sequencing+data+analysis.">Deeptools2: A next generation web server for deep-sequencing data analysis.</a> Nucleic Acids Res. 44 (W1), W160-W165 (2016).</li><li>Robinson, J. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Integrative+genomics+viewer.">Integrative genomics viewer.</a> Nat Biotechnol. 29 (1), 24-26 (2011).</li><li>Kent, W. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+human+genome+browser+at+ucsc.">The human genome browser at ucsc.</a> Genome Res. 12 (6), 996-1006 (2002).</li><li>Yu, F., Sankaran, V. G., Yuan, G. -. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=CUT&amp;RUNTools+2.0:+A+pipeline+for+single-cell+and+bulk-level+CUT&amp;RUN+and+CUT&amp;Tag+data+analysis.">CUT&amp;RUNTools 2.0: A pipeline for single-cell and bulk-level CUT&amp;RUN and CUT&amp;Tag data analysis.</a> Bioinformatics. 38 (1), 252-254 (2021).</li><li>Zhu, Q., Liu, N., Orkin, S. H., Yuan, G. -. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=CUT&amp;RUNTools:+A+flexible+pipeline+for+CUT&amp;RUN+processing+and+footprint+analysis.">CUT&amp;RUNTools: A flexible pipeline for CUT&amp;RUN processing and footprint analysis.</a> Genome Biol. 20 (1), 192 (2019).</li><li>. Nf-core/cutandrun: Nf-core/cutandrun v3.2.2 iridium ibis Available from: <a target="_blank" href="https://github.com/nf-core/cutandrun/tree/3.2.2">https://github.com/nf-core/cutandrun/tree/3.2.2</a> (2024)</li><li>Kong, N. R., Chai, L., Tenen, D. G., Bassal, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+modified+CUT&amp;RUN+protocol+and+analysis+pipeline+to+identify+transcription+factor+binding+sites+in+human+cell+lines.">A modified CUT&amp;RUN protocol and analysis pipeline to identify transcription factor binding sites in human cell lines.</a> STAR Protoc. 2 (3), 100750 (2021).</li><li>Meers, M. P., Bryson, T. D., Henikoff, J. G., Henikoff, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Improved+CUT&amp;RUN+chromatin+profiling+tools.">Improved CUT&amp;RUN chromatin profiling tools.</a> eLife. 8, e46314 (2019).</li><li>Amemiya, H. M., Kundaje, A., Boyle, A. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+encode+blacklist:+Identification+of+problematic+regions+of+the+genome.">The encode blacklist: Identification of problematic regions of the genome.</a> Sci Rep. 9 (1), 9354 (2019).</li><li>Deberardine, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=BRgenomics+for+analyzing+high-resolution+genomics+data+in+R.">BRgenomics for analyzing high-resolution genomics data in R.</a> Bioinformatics. 39 (6), btad331 (2023).</li><li>Deberardine, M., Booth, G. T., Versluis, P. P., Lis, J. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+nelf+pausing+checkpoint+mediates+the+functional+divergence+of+cdk9.">The nelf pausing checkpoint mediates the functional divergence of cdk9.</a> Nat Commun. 14 (1), 2762 (2023).</li><li>Krueger, F., James, F. O., Ewels, P. A., Afyounian, E., Schuster-Boeckler, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> FelixKrueger/TrimGalore: v0.6.7 - DOI via Zenodo. , (2021).</li><li>. Easy bam downsampling Available from: <a target="_blank" href="https://davemcg.github.io/post/easy-bam-downsampling/">https://davemcg.github.io/post/easy-bam-downsampling/</a> (2018)</li><li>. Positionbaseddownsamplesam (picard) Available from: <a target="_blank" href="https://gatk.broadinstitute.org/hc/en-us/articles/360041850311-PositionBasedDownsampleSam-Picard">https://gatk.broadinstitute.org/hc/en-us/articles/360041850311-PositionBasedDownsampleSam-Picard</a> (2020)</li></ol>

The ability to map protein occupancy on chromatin is fundamental to conducting mechanistic studies in the field of chromatin biology. As laboratories adopt new wet lab techniques to profile chromatin, the ability to analyze sequencing data from those wet lab experiments becomes a common bottleneck for wet lab scientists. Therefore, we describe an introductory step-by-step protocol to enable bioinformatics beginners to overcome the analysis bottleneck, and initiate analysis and quality control checks of their own CUT&amp;...

The ability to measure interactions between proteins and genomic DNA is fundamental to understanding the biology of chromatin regulation. Effective assays that measure chromatin occupancy for a given protein provide at least two key pieces of information: i) genomic localization and ii) protein abundance at a given genomic region. Tracking the recruitment and localization changes of a protein of interest in chromatin can reveal direct target loci of the protein and reveal mechanistic roles of that protein in chromatin-based biological processes such as regulation of transcription, DNA repair, or DNA replication. The techniques available today to profile protein-DNA interactions are enabling researchers to explore regulation at unprecedented resolution. Such technical advances have been enabled through the introduction of new chromatin profiling techniques that include the development of Cleavage Under Targets and Release Using Nuclease (CUT&#38;RUN) by the Henikoff laboratory. CUT&#38;RUN offers several technical advantages over conventional chromatin immunoprecipitation (ChIP) that include lower cell input requirements, lower sequencing depth requirements, and increased sensitivity with reduced background signal due to a lack of cross-linking agents that otherwise mask antibody epitopes. Adopting this technique to study chromatin regulation requires a thorough understanding of the principle underlying the technique, and an understanding of how to analyze, validate, and interpret CUT&#38;RUN data.
The CUT&#38;RUN procedure begins with binding of cells to Concanavalin A conjugated to magnetic beads to enable manipulation of low cell numbers throughout the procedure. Isolated cells are permeabilized using a mild detergent to facilitate introduction of an antibody that targets the protein of interest. Micrococcal nuclease (MNase) is then recruited to the bound antibody using a Protein A or Protein A/G tag tethered to the enzyme. Calcium is introduced to initiate enzymatic activity. MNase digestion results in mono-nucleosomal DNA-protein complexes. Calcium is subsequently chelated to end the digestion reaction, and short DNA fragments from the MNase digestion are released from nuclei, then subjected to DNA purification, library preparation, and high-throughput sequencing1 (Figure 1).
In silico approaches to map and quantify protein occupancy across the genome have developed in parallel with the wet lab approaches used to enrich those DNA-protein interactions. Identification of regions of enriched signals (peaks) is one of the most critical steps in the bioinformatics analysis. Initial ChIP-seq analysis methods used algorithms such as MACS2 and SICER3, which employed statistical models to distinguish bona fide protein-DNA binding sites from background noise. However, the lower background noise and higher resolution of CUT&#38;RUN data render some peak calling programs employed in ChIP-seq analysis unsuitable for CUT&#38;RUN analysis4. This challenge highlights the need for new tools better suited for the analysis of CUT&#38;RUN data. SEACR4 represents one such tool recently developed to enable peak calling from CUT&#38;RUN data while overcoming limitations associated with tools typically employed toward ChIP-seq analysis.
Biological interpretations from CUT&#38;RUN sequencing data are drawn from the outputs downstream of peak calling in the analysis pipeline. Several functional annotation programs can be implemented to predict the potential biological relevance of the called peaks from CUT&#38;RUN data. For example, the Gene Ontology (GO) project provides well-established functional identification of genes of interest5,6,7. Various software tools and resources facilitate GO analysis to reveal genes and gene sets enriched amongst CUT&#38;RUN peaks8,9,10,11,12,13,14. Furthermore, visualization software such as Deeptools15, Integrative genomics viewer (IGV)16, and UCSC Genome Browser17 enable visualization of signal distribution and patterns at regions of interest across the genome.
The ability to draw biological interpretations from CUT&#38;RUN data depends critically upon validation of data quality. Critical components to validate include the assessment of: i) CUT&#38;RUN library sequencing quality, ii) replicate similarity, and iii) signal distribution at peak centers. Completing the validation of all three components is crucial to ensure the reliability of CUT&#38;RUN library samples and downstream analysis results. Therefore, it is essential to establish introductory CUT&#38;RUN analysis guides to enable bioinformatics beginners and wet lab researchers to conduct such validation steps as part of their standard CUT&#38;RUN analysis pipelines.
Alongside the development of wet lab CUT&#38;RUN experiment, various in silico CUT&#38;RUN analysis pipelines, such as CUT&#38;RUNTools 2.018,19, nf-core/cutandrun20, and CnRAP21, have been developed to support CUT&#38;RUN data analysis. These tools provide powerful approaches to analyzing single-cell and bulk CUT&#38;RUN and CUT&#38;Tag datasets. However, the relatively complex modular program structure and required familiarity with multiple programming languages to conduct these analysis pipelines may hinder adoption by bioinformatics beginners who seek to thoroughly understand the CUT&#38;RUN analysis steps and customize their own pipelines. Circumvention of this barrier requires a new introductory CUT&#38;RUN analysis pipeline that is provided in simple step-by-step scripts encoded using a simple single programming language.
In this article, we describe a simple single-language CUT&#38;RUN analysis pipeline protocol that provides step-by-step scripts supported with detailed descriptions to enable new and novice users to conduct CUT&#38;RUN sequencing analysis. Programs used in this pipeline are publicly available by the original developer groups. Major steps described in this protocol include read alignment, peak calling, functional analysis and, most critically, validation steps to assess sample quality to determine data suitability and reliability for biological interpretation (Figure 2). Furthermore, this pipeline provides users the opportunity to cross-reference analysis results against publicly available CUT&#38;RUN datasets. Ultimately, this CUT&#38;RUN analysis pipeline protocol serves as an introductory guide and reference for bioinformatic analysis beginners and wet lab researchers.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>bedGraphToBigWig</td>
			<td>ENCODE</td>
			<td>https://hgdownload.soe.ucsc.edu/admin/exe/</td>
			<td>Software to compress and convert readcounts bedGraph to bigWig</td>
		</tr>
		<tr>
			<td>bedtools-2.31.1</td>
			<td>The Quinlan Lab @ the U. of Utah</td>
			<td>https://bedtools.readthedocs.io/en/latest/index.html</td>
			<td>Software to process bam/bed/bedGraph files</td>
		</tr>
		<tr>
			<td>bowtie2 2.5.4</td>
			<td>Johns Hopkins University</td>
			<td>https://bowtie-bio.sourceforge.net/bowtie2/index.shtml</td>
			<td>Software to build bowtie index and perform alignment</td>
		</tr>
		<tr>
			<td>CollectInsertSizeMetrics (Picard)</td>
			<td>Broad institute</td>
			<td>https://github.com/broadinstitute/picard</td>
			<td>Software to perform insert size distribution analysis</td>
		</tr>
		<tr>
			<td>Cutadapt</td>
			<td>NBIS</td>
			<td>https://cutadapt.readthedocs.io/en/stable/index.html</td>
			<td>Software to perform adapter trimming</td>
		</tr>
		<tr>
			<td>Deeptoolsv3.5.1</td>
			<td>Max Planck Institute</td>
			<td>https://deeptools.readthedocs.io/en/develop/index.html</td>
			<td>Software to perform Pearson coefficient correlation analysis, Principal component analysis, and Heatmap/average plot analysis</td>
		</tr>
		<tr>
			<td>FastQC Version 0.12.0</td>
			<td>Babraham Bioinformatics</td>
			<td>https://github.com/s-andrews/FastQC</td>
			<td>Software to check quality of fastq file</td>
		</tr>
		<tr>
			<td>Intervenev0.6.1</td>
			<td>Computational Biology &#38; Gene regulation - Mathelier group</td>
			<td>https://intervene.readthedocs.io/en/latest/index.html</td>
			<td>Software to perform venn diagram analysis using peak files</td>
		</tr>
		<tr>
			<td>MACSv2.2.9.1</td>
			<td>Chan Zuckerberg initiative</td>
			<td>https://github.com/macs3-project/MACS/tree/macs_v2</td>
			<td>Software to call peaks</td>
		</tr>
		<tr>
			<td>MACSv3.0.2</td>
			<td>Chan Zuckerberg initiative</td>
			<td>https://github.com/macs3-project/MACS/tree/master</td>
			<td>Software to call peaks</td>
		</tr>
		<tr>
			<td>Samtools-1.21</td>
			<td>Wellcome Sanger Institute</td>
			<td>https://github.com/samtools/samtools</td>
			<td>Software to process sam/bam files</td>
		</tr>
		<tr>
			<td>SEACRv1.3</td>
			<td>Howard Hughes Medial institute</td>
			<td>https://github.com/FredHutch/SEACR</td>
			<td>Software to call peaks</td>
		</tr>
		<tr>
			<td>SRA Toolkit Release 3.1.1</td>
			<td>NCBI</td>
			<td>https://github.com/ncbi/sra-tools</td>
			<td>Software to download SRR from GEO</td>
		</tr>
		<tr>
			<td>Trim_Galore v0.6.10</td>
			<td>Babraham Bioinformatics</td>
			<td>https://github.com/FelixKrueger/TrimGalore</td>
			<td>Software to perform quality and atapter trimming</td>
		</tr>
	</tbody>
</table>

introductory analysis and validation of cut&#38;run sequencing data

The CUT&amp;RUN technique facilitates detection of protein-DNA interactions across the genome. Typical applications of CUT&amp;RUN include profiling changes in histone tail modifications or mapping transcription factor chromatin occupancy. Widespread adoption of CUT&amp;RUN is driven, in part, by technical advantages over conventional ChIP-seq that include lower cell input requirements, lower sequencing depth requirements, and increased sensitivity with reduced background signal due to a lack of cross-linking agents that otherwise mask antibody epitopes. Widespread adoption of CUT&amp;RUN has also been achieved through the generous sharing of reagents by the Henikoff lab and the development of commercial kits to accelerate adoption for beginners. As technical adoption of CUT&amp;RUN increases, CUT&amp;RUN sequencing analysis and validation become critical bottlenecks that must be surmounted to enable complete adoption by predominantly wet lab teams. CUT&amp;RUN analysis typically begins with quality control checks on raw sequencing reads to assess sequencing depth, read quality, and potential biases. Reads are then aligned to a reference genome sequence assembly, and several bioinformatics tools are subsequently employed to annotate genomic regions of protein enrichment, confirm data interpretability, and draw biological conclusions. Although multiple in silico analysis pipelines have been developed to support CUT&amp;RUN data analysis, their complex multi-module structure and usage of multiple programming languages render the platforms difficult for bioinformatics beginners who may lack familiarity with multiple programming languages but wish to understand the CUT&amp;RUN analysis procedure and customize their analysis pipelines. Here, we provide a single-language step-by-step CUT&amp;RUN analysis pipeline protocol designed for users with any level of bioinformatics experience. This protocol includes completing critical quality checks to validate that the sequencing data is suitable for biological interpretation. We expect that following the introductory protocol provided in this article combined with downstream peak annotation will allow users to draw biological insights from their own CUT&amp;RUN datasets.

Author Spotlight: Investigating the Role of Repetitive DNA Misregulation in Cancer Initiation and Immunotherapy Resistance

Introductory Analysis and Validation of CUT&#38;RUN Sequencing Data

Our research explores the role of repetitive DNA elements in cancer initiation and progression. We ask whether these elements are directly involved in mechanisms of cancer initiation, and, if so, can we use this information to improve patient outcomes. Many next-generation sequencing technology and analytic pipeline has been developed to capture the feature and changes of the epigenetic factors such as transcriptional or the post-translational chromatin changes and the binding at the specific localization on the chromatin and sequences of DNA or RNA which are captured by the defensive protein in both bulk and single-cell resolution.The technologies used to advance research in our field are those that can map the genomic positions of modified histones or transcription factors on chromatin. For our purposes, the most relevant technique within this umbrella of technologies is a technique called CUT&RUN. Many publicly available CUT&RUN analysis pipeline are relatively complex for the bioinformatic beginner who want to understand the flow of the analysis and to develop their own analysis pipeline for their project.But our CUT&RUN analysis pipeline are developed in step-by-step manner in single program language, which are much easier to understand the flow and easier to modify to get their publication-quality data and to build up the experience for the bioinformatics. We are key focusing on to find the answer about how chronic double-stranded RNA expression affect the biomimicry evasion and the immunotherapy rates in cancer development, and how TP53 protein can regulate the chronic double-stranded RNA expression and the biomimicry evasion. And at the end, we want to understand how we can cure the immunotherapy resistant cancer by understanding and targeting the chronic double-stranded RNA expression and the biomimicry evasion process.

Quality and adapter trimming retains reads with high sequencing quality 
High-throughput sequencing techniques are prone to generating sequencing errors such as sequence &#39;mutations&#39; in reads. Furthermore, sequencing adapter dimers can be enriched in sequencing datasets due to poor adapter removal during library preparation. Excessive sequencing errors, such as read mutations, generation of reads shorter than required for proper mapping, and enrichment of adapter dimers, can increase read map...

This protocol guides bioinformatics beginners through an introductory CUT&amp;RUN analysis pipeline that enables users to complete an initial analysis and validation of CUT&amp;RUN sequencing data. Completing the analysis steps described here, combined with downstream peak annotation, will allow users to draw mechanistic insights into chromatin regulation.

The CUT&amp;RUN technique facilitates detection of protein-DNA interactions across the genome. Typical applications of CUT&amp;RUN include profiling ...

introductory-analysis-and-validation-of-cutrun-sequencing-data

Research

JoVE Journal

Genetics

Detecting Protein-DNA Interactions Using CUT&RUN Analysis

To begin, download the compressed analysis pipeline from GitHub by typing the command in the terminal. After downloading, enter the command to decompress the zip file. Once the decompression is complete, remove the zip file and change the decompressed folder name.Then, set executable permissions for all the shell scripts in the working directory. From the available options, locate the appropriate installation shell script. Look for a script named Script_01_Installation.sh.Open the terminal and type echo SHELL to check the default shell used in the active terminal. If bash is not the default shell, set it by typing chsh s in the terminal. Type the fast queue download script in the terminal, or drag the script file into the terminal and press enter.After running the script, check the log file for error messages. Now, enter the trimming script in the terminal and execute it. Next, type the bowtie2 index creation script in the terminal and execute it.Finally, run the bowtie2 mapping script and press enter. Type the bam filtering and sorting script in the terminal and execute it. Next, convert bam to BEDPE, BED, and bedGraph formats by typing the appropriate command in the terminal and pressing enter.Perform normalization using the SRPMC method, then run the insert size analysis script. Perform peak calling using the MACS command in the terminal. Afterward, type the command for peak calling using SEACR and press enter.Create peak BED files by using the appropriate command in the terminal, and then filter the peaks. Next, enter the command to merge MACS peaks. Type the command in the terminal to merge SEACR peaks.Generate a correlation plot, and run the command for principle component analysis in the terminal. Then, create a Venn diagram comparing different methods using the appropriate command. Enter the command to generate a Venn diagram for replicates.Finally, generate a focused peak center heat map. Then, generate a whole peak center heat map.