Genome-wide Profiling Methods Research Articles

Ligand-induced native G-quadruplex stabilization impairs transcription initiation.

G-quadruplexes (G4s) are noncanonical DNA secondary structures formed through the self-association of guanines, and G4s are distributed widely across the genome. G4 participates in multiple biological processes including gene transcription, and G4-targeted ligands serve as potential therapeutic agents for DNA-targeted therapies. However, genome-wide studies of the exact roles of G4s in transcriptional regulation are still lacking. Here, we establish a sensitive G4-CUT&Tag method for genome-wide profiling of native G4s with high resolution and specificity. We find that native G4 signals are cell type–specific and are associated with transcriptional regulatory elements carrying active epigenetic modifications. Drug-induced promoter-proximal RNA polymerase II pausing promotes nearby G4 formation. In contrast, G4 stabilization by G4-targeted ligands globally reduces RNA polymerase II occupancy at gene promoters as well as nascent RNA synthesis. Moreover, ligand-induced G4 stabilization modulates chromatin states and impedes transcription initiation via inhibition of general transcription factors loading to promoters. Together, our study reveals a reciprocal genome-wide regulation between native G4 dynamics and gene transcription, which will deepen our understanding of G4 biology toward therapeutically targeting G4s in human diseases.

Rapid and Low-Input Profiling of Histone Marks in Plants Using Nucleus CUT&Tag.

Characterizing genome-wide histone posttranscriptional modifications and transcriptional factor occupancy is crucial for deciphering their biological functions. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a powerful method for genome-wide profiling of histone modifications and transcriptional factor-binding sites. However, the current ChIP-seq experimental procedure in plants requires significant material and several days for completion. CUT&Tag is an alternative method of ChIP-seq for low-sample and single-cell epigenomic profiling using protein A-Tn5 transposase fusion proteins (PAT). In this study, we developed a nucleus CUT&Tag (nCUT&Tag) protocol based on the live-cell CUT&Tag technology. Our results indicate that nCUT&Tag could be used for histone modifications profiling in both monocot rice and dicot rapeseed using crosslinked or fresh tissues. In addition, both active and repressive histone marks such as H3K4me3 and H3K9me2 can be identified using our nCUT&Tag. More importantly, all the steps in nCUT&Tag can be finished in only 1 day, and the assay can be performed with as little as 0.01 g of plant tissue as starting materials. Therefore, our results demonstrate that nCUT&Tag is an efficient alternative strategy for plant epigenomic studies.

Different Plant Species Have Common Sequence Features Related to mRNA Degradation Intermediates.

mRNA degradation is an important cellular mechanism involved in the control of gene expression. Several genome-wide profiling methods have been developed for detecting mRNA degradation in plants and animals. However, because many of these techniques use poly (A) mRNA for library preparation, degradation intermediates are often only detected near the 3'-ends of transcripts. Previously, we developed the Truncated RNA End Sequencing (TREseq) method using Arabidopsis thaliana, and demonstrated that this method ameliorates 3'-end bias. In analyses using TREseq, we observed G-rich sequences near the 5'-ends of degradation intermediates. However, this finding remained to be confirmed in other plant species. Hence, in this study, we conducted TREseq analyses in Lactuca sativa (lettuce), Oryza sativa (rice) and Rosa hybrida (rose). These species including A. thaliana were selected to encompass a diverse range in the angiosperm phylogeny. The results revealed similar sequence features near the 5'-ends of degradation intermediates, and involvement of translation process in all four species. In addition, homologous genes have similar efficiencies of mRNA degradation in different plants, suggesting that similar mechanisms of mRNA degradation are conserved across plant species. These strong sequence features were not observed in previous degradome analyses among different species in plants.

Endogenous DNA Double-Strand Breaks during DNA Transactions: Emerging Insights and Methods for Genome-Wide Profiling.

DNA double-strand breaks (DSBs) jeopardize genome integrity and can—when repaired unfaithfully—give rise to structural rearrangements associated with cancer. Exogenous agents such as ionizing radiation or chemotherapy can invoke DSBs, but a vast amount of breakage arises during vital endogenous DNA transactions, such as replication and transcription. Additionally, chromatin looping involved in 3D genome organization and gene regulation is increasingly recognized as a possible contributor to DSB events. In this review, we first discuss insights into the mechanisms of endogenous DSB formation, showcasing the trade-off between essential DNA transactions and the intrinsic challenges that these processes impose on genomic integrity. In the second part, we highlight emerging methods for genome-wide profiling of DSBs, and discuss future directions of research that will help advance our understanding of genome-wide DSB formation and repair.

Alternative polyadenylation analysis in animals and plants: newly developed strategies for profiling, processing and validation.

Alternative polyadenylation is an essential RNA processing event that contributes significantly to regulation of transcriptome diversity and functional dynamics in both animals and plants. Here we review newly developed next generation sequencing methods for genome-wide profiling of alternative polyadenylation (APA) sites, bioinformatics pipelines for data processing and both wet and dry laboratory approaches for APA validation. The library construction methods LITE-Seq (Low-Input 3'-Terminal sequencing) and PAC-seq (PolyA Click sequencing) tag polyA+ cDNA, while BAT-seq (BArcoded, three-prime specific sequencing) and PAPERCLIP (Poly(A) binding Protein-mediated mRNA 3′End Retrieval by CrossLinking ImmunoPrecipitation) enrich polyA+ RNA. Interestingly, only WTTS-seq (Whole Transcriptome Termini Site sequencing) targets both polyA+ RNA and polyA+ cDNA. Varieties of bioinformatics pipelines are well established to pursue read quality control, mapping, clustering, characterization and pathway analysis. The RHAPA (RNase H alternative polyadenylation assay) and 3'RACE-seq (3' rapid amplification of cDNA end sequencing) methods directly validate APA sites, while WTSS-seq (whole transcriptome start site sequencing), RNA-seq (RNA sequencing) and public APA databases can serve as indirect validation methods. We hope that these tools, pipelines and resources trigger huge waves of interest in the research community to investigate APA events underlying physiological, pathological and psychological changes and thus understand the information transfer events from genome to phenome relevant to economically important traits in both animals and plants.

Enzymatic methods for genome-wide profiling of protein binding sites

Genome-wide mapping of protein-DNA interactions is a staple approach in many areas of modern molecular biology. Genome-wide profiles of protein-binding sites are most commonly generated by chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq). Although ChIP-seq has played a central role in studying genome-wide protein binding, recent work has highlighted systematic biases in the technique that warrant technical and interpretive caution and underscore the need for orthogonal techniques to both confirm the results of ChIP-seq studies and uncover new insights not accessible to ChIP. Several such techniques, based on genetic or immunological targeting of enzymatic activity to specific genomic loci, have been developed. Here, we review the development, applications and future prospects of these methods as complements to ChIP-based approaches and as powerful techniques in their own right.

DNA Methylation in Embryo Development: Epigenetic Impact of ART (Assisted Reproductive Technologies).

DNA methylation can be considered a component of epigenetic memory with a critical role during embryo development, and which undergoes dramatic reprogramming after fertilization. Though it has been a focus of research for many years, the reprogramming mechanism is still not fully understood. Recent results suggest that absence of maintenance at DNA replication is a major factor, and that there is an unexpected role for TET3-mediated oxidation of 5mC to 5hmC in guarding against de novo methylation. Base-resolution and genome-wide profiling methods are enabling more comprehensive assessments of the extent to which ART might impair DNA methylation reprogramming, and which sequence elements are most vulnerable. Indeed, as we also review here, studies showing the effect of culture media, ovarian stimulation or embryo transfer on the methylation pattern of embryos emphasize the need to face ART-associated defects and search for strategies to mitigate adverse effects on the health of ART-derived children.

A novel method for genome-wide profiling of dynamic host-pathogen interactions using 3\u2032 end enriched RNA-seq

Marek’s disease is a contagious lymphoproliferative disease of chickens and typical model of viral oncogenesis. Mapping changes or different states over the course of infection for both host and pathogen would provide important insights into dynamic host-pathogen interactions. Here we introduced 3′ end enriched RNA-seq as a novel method to study host-pathogen interactions in chicken embryo fibroblasts cells challenged with Marek’s disease virus. The method allowed accurate profiling of gene expression and alternative polyadenylation sites for host and pathogen simultaneously. We totally identified 476 differentially expressed genes and 437 APA switching genes in host, including switching in tandem 3′ UTRs and switching between coding region and 3′ UTR. Most of these genes were related to innate immunity, apoptosis and metabolism, but two sets of genes overlapped a little, suggesting two complementary mechanisms in gene regulation during MDV infection. In summary, our results provided a relatively comprehensive insight into dynamic host-pathogen interactions in regulation of gene transcription during infection of Marek’s disease virus and suggested that 3′ end enriched RNA-seq was a promising method to investigate global host-pathogen interactions.

BLISS is a versatile and quantitative method for genome-wide profiling of DNA double-strand breaks

Precisely measuring the location and frequency of DNA double-strand breaks (DSBs) along the genome is instrumental to understanding genomic fragility, but current methods are limited in versatility, sensitivity or practicality. Here we present Breaks Labeling In Situ and Sequencing (BLISS), featuring the following: (1) direct labelling of DSBs in fixed cells or tissue sections on a solid surface; (2) low-input requirement by linear amplification of tagged DSBs by in vitro transcription; (3) quantification of DSBs through unique molecular identifiers; and (4) easy scalability and multiplexing. We apply BLISS to profile endogenous and exogenous DSBs in low-input samples of cancer cells, embryonic stem cells and liver tissue. We demonstrate the sensitivity of BLISS by assessing the genome-wide off-target activity of two CRISPR-associated RNA-guided endonucleases, Cas9 and Cpf1, observing that Cpf1 has higher specificity than Cas9. Our results establish BLISS as a versatile, sensitive and efficient method for genome-wide DSB mapping in many applications.

Size no longer matters.

A method for genome-wide profiling of small RNAs, such as microRNAs, from single cells is developed and applied to human cancer and stem cells.

Genome-wide footprinting: ready for prime time?

High-throughput sequencing technologies have allowed many gene locus-level molecular biology assays to become genome-wide profiling methods. DNA-cleaving enzymes such as DNase I have been used to probe accessible chromatin. The accessible regions contain functional regulatory sites, including promoters, insulators and enhancers. Deep sequencing of DNase-seq libraries and computational analysis of the cut profiles have been used to infer protein occupancy in the genome at the nucleotide level, a method introduced as 'digital genomic footprinting'. The approach has been proposed as an attractive alternative to the analysis of transcription factors (TFs) by chromatin immunoprecipitation followed by sequencing (ChIP-seq), and in theory it should overcome antibody issues, poor resolution and batch effects. Recent reports point to limitations of the DNase-based genomic footprinting approach and call into question the scope of detectable protein occupancy, especially for TFs with short-lived chromatin binding. The genomics community is grappling with issues concerning the utility of genomic footprinting and is reassessing the proposed approaches in terms of robust deliverables. Here we summarize the consensus as well as different views emerging from recent reports, and we describe the remaining issues and hurdles for genomic footprinting.

Abstract 374: A novel sequencing method for genome-wide profiling of 5-hydroxymethylcytosine with single-base resolution

Abstract 5-hydroxymethylcytosine (5hmC) is an epigenetic mark abundant in embryo stem cells and brain tissues. The exact biological functions of 5hmC are still under close investigation although several lines of evidence have indicated it could be involved in active DNA demethylation. Meanwhile, extensive studies have been carried out to determine its genomic distribution. A number of approaches have been developed using either affinity based enrichment, such as hMeDIP, that rely on antibody and other specific binding proteins to target 5hmC, or modified bisulfite sequencing, namely oxidative bisulfite sequencing (OxBS) and TET assisted bisulfite sequencing (TAB-sequencing). However, all those methods have limitations which hamper their application. For example, affinity based methods lack single base resolution while modified bisulfite sequencing methods require efficient chemical or enzymatic oxidation which cannot be easily achieved or guaranteed. As an alternative, we have developed a novel genome-wide sequencing method that utilizes an enzyme based modification approach coupled with bisulfite-sequencing for detecting 5hmC. This methodology allows quantification of 5hmC levels with single CpG resolution and can also be employed for locus-specific assays. Using this method, we were able to map and quantify 5hmC sites at the genomic scale for several different biological samples. This novel method can determine the exact location and abundance of 5hmC, which will facilitate our understanding of 5hmC in regulating gene expression in different biological contexts. Citation Format: Xueguang Sun, Tzu-Hung Chung, Yap Ching Chew, Darany Tan, Xi-Yu Jia. A novel sequencing method for genome-wide profiling of 5-hydroxymethylcytosine with single-base resolution. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 374. doi:10.1158/1538-7445.AM2014-374

JBP1-seq: A fast and efficient method for genome-wide profiling of 5hmC

We developed a novel approach, J-binding protein 1 sequencing (JBP1-seq), that combines the benefits of an improved recombinant JBP1 protein, Nextera-based library construction, and next-generation sequencing (NGS) for genome-wide profiling of 5-hydroxymethylcytosine (5hmC). Compared with the original JBP1, this new recombinant JBP1 was biotinylated in vivo and conjugated to magnetic beads via biotin–streptavidin interactions. These modifications allowed a more efficient and consistent pull-down of β-glucosyl-5-hydroxymethylcytosine (β-glu-5hmC), and sequence-ready libraries can be generated within 4.5h from DNA inputs as low as 50ng. 5hmC enrichment of human brain DNA using the new JBP1 resulted in over 25,000 peaks called, which is significantly higher than the 4003 peaks enriched using the old JBP1. Comparison of the technical duplicates and validations with other platforms indicated the results are reproducible and reliable. Thus, JBP1-seq provides a fast, efficient, and cost-effective method for accurate 5hmC genome-wide profiling.

A Comprehensive Profile of ChIP-Seq-Based STAT1 Target Genes Suggests the Complexity of STAT1-Mediated Gene Regulatory Mechanisms

Interferon-gamma (IFNγ) plays a key role in macrophage activation, T helper and regulatory cell differentiation, defense against intracellular pathogens, tissue remodeling, and tumor surveillance. The diverse biological functions of IFNγ are mediated by direct activation of signal transducer and activator of transcription 1 (STAT1) as well as numerous downstream effector genes. Because a perturbation in STAT1 target gene networks is closely associated with development of autoimmune diseases and cancers, it is important to characterize the global picture of these networks. Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) provides a highly efficient method for genome-wide profiling of DNA-binding proteins. We analyzed the STAT1 ChIP-Seq dataset of IFNγ-stimulated HeLa S3 cells derived from the ENCODE project, along with transcriptome analysis on microarray. We identified 1,441 stringent ChIP-Seq peaks of protein-coding genes. They were located in the promoter (21.5%) and more often in intronic regions (72.2%) with an existence of IFNγ-activated site (GAS) elements. Among the 1,441 STAT1 target genes, 212 genes are known IFN-regulated genes (IRGs) and 194 genes (13.5%) are actually upregulated in response to IFNγ by transcriptome analysis. The panel of upregulated genes constituted IFN-signaling molecular networks pivotal for host defense against infections, where interferon-regulatory factor (IRF) and STAT transcription factors serve as a hub on which biologically important molecular connections concentrate. The genes with the peak location in intronic regions showed significantly lower expression levels in response to IFNγ. These results indicate that the binding of STAT1 to GAS is not sufficient to fully activate target genes, suggesting the high complexity of STAT1-mediated gene regulatory mechanisms.

Abstract 5098: Identification of copy number alterations associated with the progression and clinical outcomes of oral cavity squamous cell carcinoma

Abstract The molecular mechanisms contributing to the progression of oral cavity squamous cell carcinoma (OSCC) remain poorly understood. The primary treatment for advanced-stage OSCC is radical surgery with post-operative chemoradiation. However, the prognoses are usually discouraging. If we can better understand the characteristics of OSCCs, this may ultimately help clinicians to provide OSCC patients with more appropriate treatment. The advent of genome-wide profiling methods has opened up new possibility to catalogue chromosomal abnormal regions at high resolution. In the present study, 72 advanced-stage OSCCs were analyzed using two microarray platforms (26 cases with Affymetrix SNP 500K and 46 cases with Affymetrix SNP 6.0). Genomic copy number was determined by calculating the intensity ratio between the OSCC and a set of 270 normal DNAs from HapMap project. Segmented copy numbers for each tumor were inferred with the hidden markov model (HMM) using the Partek GS software. To identify regions for copy-number alterations (CNAs), GISTIC (genomic identification of significant targets in cancer) was applied to determine the significance. With this approach, 46 (14 gain and 32 loss) and 208 (48 gain and 160 loss) significant regions (FDR q-value < 0.25) were identified from SNP 500K and SNP 6.0 platform, respectively. Among them, 15 focal CNAs (8 gains (3q26.1, 7p11.2, 11q13.3, 13q22.1, 15q11.2, 17p13.3, 18p11.23-18p11.31, and 20p11.22) and 7 losses (1q21.1, 6p21.33, 7q31.3, 12q14.1, 2 focal regions in 14q11.2, 15q11.2)) are identical in these 2 platforms. Amplification of 11q13.3 containing a known oncogene CCND1 was confirmed by fluorescence in situ hybridization (FISH) with Vysis CCND1/CEP11 FISH probe kit. The clinicopathologic significances of these 15 CNAs were examined. Gain of 3q26.1 and 18p11.23-18p11.31, and loss of 15q11.2 were correlated significantly with lymph node metastasis. Furthermore, loss of two focal regions in 14q11.2 was significantly associated with disease-free and overall survival. In conclusion, CNAs associated with the progression and clinical outcomes of Taiwanese OSCC were identified by genomic approach using high-density SNP array. These results provide useful clues for future characterization of important oncogenes or tumor suppressor genes associated with the behaviors of the disease. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5098. doi:1538-7445.AM2012-5098

Hypermethylation of the gene LARP2 for noninvasive prenatal diagnosis of β-thalassemia based on DNA methylation profile

In order to identify epigenetic markers of β-thalassemia, a genome-wide profiling method named differential methylation hybridization was used to search these differentially methylated genes. Unsupervised hierarchical clustering and molecular annotation system were used to analyze the data, and methylation-specific PCR and real-time PCR were used to confirm the differentially methylated genes. This system was validated by detecting 13 cases, 10 of which were homo-zygous β-thalassaemia. Totally 113 genes were identified as methlyation-enriched genes (ratio ≥ 2.0, P < 0.05) and 96 genes were identified as hypomethylated genes in both groups (ratio ≤ 0.5, P < 0.05). The promoter of the gene of La ribonucleoprotein domain family (LARP2) was significantly hypermethylated in β-thalassemia, and the expression of LARP2 was significantly lower in β-thalassemia. Hypermethylation of the LARP2 promoter was correlated with its lower expression in β-thalassemia and our chip-based DNA methylation detection system can provide earlier diagnosis of β-thalassemia using this epigenetic marker.

Transcriptome Analysis of Neural Progenitor Cells by a Genetic Dual Reporter Strategy

Global analysis of stem/progenitor cells promises new insight into mechanisms that govern self-renewal and cellular potential, an unresolved question of stem/progenitor cell biology. Despite rapid advance of genome-wide profiling methods, the difficulty in cell purification remains a major challenge for global analysis of somatic stem/progenitor cells. Genetic tagging with a reporter provides a powerful tool for identification and isolation of a specific mature cell type; however, for stem/progenitor cells, reporter retention by progeny may be a concern for impurity. Here, we describe a genetic system combining a progenitor cell specific label with a second tag for marking differentiation. We present evidence that differential labeling of neural progenitor cells and their progeny enables prospective purification of these two cell types, whereas isolation based on a single marker compromises the purity of the intended progenitor population. Comparative expression profiling between the purified progenitors and progeny documents a neural progenitor cell transcriptome and uncovers an important role of Tyro3/Axl/Mer receptor tyrosine kinases in the maintenance of neural progenitor cells. This study establishes a general strategy for isolation of somatic stem/progenitor cells and provides a transcriptome database of neural progenitor cells useful for identification of causal factors of neural progenitor cell state, global dissection of epigenetic control of cellular potential, as well as for developing biomarkers or targets of brain cancer stem/initiating cells for therapeutic interventions.

Mapping Allele-Specific DNA Methylation: A New Tool for Maximizing Information from GWAS

In this issue of The Journal, an article by Schalkwyk et al.(1) shows the landscape of allele-specific DNA methylation (ASM) in the human genome. ASM has long been studied as a hallmark of imprinted genes, and a chromosome-wide version of this phenomenon occurs, in a random fashion, during X chromosome inactivation in female cells. But the type of ASM motivating the study by Schalkwyk et al. is different. They used a high-resolution, methylation-sensitive SNP array (MSNP) method for genome-wide profiling of ASM in total peripheral-blood leukocytes (PBL) and buccal cells from a series of monozygotic twin pairs. Their data bring a new level of detail to our knowledge of a newly recognized phenomenon-nonimprinted, sequence-dependent ASM. They document the widespread occurrence of this phenomenon among human genes and discuss its basic implications for gene regulation and genetic-epigenetic interactions. But this paper and recent work from other laboratories(2,3) raises the possibility of a more immediate and practical application for ASM mapping, namely to help extract maximum information from genome-wide association studies.

Genome-wide profiling of promoter methylation in human

DNA methylation in the promoter region of a gene is associated with a loss of that gene's expression and plays an important role in gene silencing. The inactivation of tumor-suppressor genes by aberrant methylation in the promoter region is well recognized in carcinogenesis. However, there has been little study in this area when it comes to genome-wide profiling of the promoter methylation. Here, we developed a genome-wide profiling method called Microarray-based Integrated Analysis of Methylation by Isoschizomers to analyse the DNA methylation of promoter regions of 8091 human genes. With this method, resistance to both the methylation-sensitive restriction enzyme HpaII and the methylation-insensitive isoschizomer MspI was compared between samples by using a microarray with promoter regions of the 8091 genes. The reliability of the difference in HpaII resistance was judged using the difference in MspI resistance. We demonstrated the utility of this method by finding epigenetic mutations in cancer. Aberrant hypermethylation is known to inactivate tumour suppressor genes. Using this method, we found that frequency of the aberrant promoter hypermethylation in cancer is higher than previously hypothesized. Aberrant hypomethylation is known to induce activation of oncogenes in cancer. Genome-wide analysis of hypomethylated promoter sequences in cancer demonstrated low CG/GC ratio of these sequences, suggesting that CpG-poor genes are sensitive to demethylation activity in cancer.