Discovery of transcription factors and regulatory regions driving in vivo tumor development by ATAC-seq and FAIRE-seq open chromatin profiling

A narrative review of epigenetic marker in H3K27ac and its emerging potential as a therapeutic target in cancer

Histone acetylation, particularly H3 K27 acetylation (H3K27ac), is a critical post-translational modification that regulates chromatin structure and gene expression, which plays a significant role in various cancers, including breast, colon, lung, hepatocellular, and prostate cancer. However, the mechanisms of H3K27ac in tumorigenesis are not yet comprehensive, especially its epigenetic mechanisms. This review endeavors to discuss findings on the involvement of H3K27ac in carcinogenesis within the past 5 years through a literature search using academic databases such as Web of Science. Firstly, we provide an overview of the diverse landscape of histone modifications, emphasizing the distinctive characteristics and critical significance of H3K27ac. Secondly, we summarize and compare advanced high-throughput sequencing technologies that have been utilized in the construction of the H3K27ac epigenetic map. Thirdly, we elucidate the role of H3K27ac in mediating gene transcription. Fourthly, we venture into the potential molecular mechanism of H3K27ac in cancer development. Finally, we engage in discussing future therapeutic approaches in oncology, with a spotlight on strategies that harness the potential of H3K27 modifications. In conclusion, this review comprehensively summarizes the characteristics of H3K27ac and underscores its pivotal role in cancer, providing valuable insights into its potential as a therapeutic target for cancer intervention.

Histone H4 lysine 20 monomethylation is not a mark of transcriptional silencers

Transcriptional silencers are cis-regulatory elements that downregulate the expression of target genes. Although thousands of silencers have been identified experimentally, a predictive chromatin signature of silencers has not been found. H4K20me1 previously was reported to be highly enriched among human silencers, but our reanalysis of those data using an appropriate background revealed that the enrichment is only marginal. We generated H4K20me1 ChIP-seq profiles in Drosophila S2 cells, which similarly showed that H4K20me1 does not mark Drosophila silencers and instead is associated with active transcription. Silencers remain a poorly annotated, difficult to predict class of cis-regulatory elements whose specific chromatin features remain to be identified.

Is the Cis-Element CACCC-Box a Master Regulatory Element during Cardiovascular Disease? A Bioinformatics Approach from the Perspective of the Krüppel-like Family of Transcription Factors

The CACCC-box motif emerges as a pivotal cis-regulatory element implicated in diverse developmental processes and diseases, particularly cardiovascular diseases (CVDs). This study centers on the intricate interplay between the CACCC-box and its binding proteins such as: the Krüppel-Like Family (KLF) of transcription factors as primary effectors in the context of CVDs. Our analysis was through a bioinformatics approach, which revealed significant transcriptional activity among KLF subgroup 2, exhibiting the highest number of interactions focusing on the established roles: pluripotency, cancer, and cardiovascular development and diseases. Our analysis reveals KLF's interactions with GATA4, MEF2C, NKX2.5 and other ~90 potential genes that participate in the regulation of the hypertrophic environment (or CVDs' Environment). Also, the GO analysis showed that genes containing the motif CACCC were enriched for multiple CVDs; in combination with STRING analysis, these results pointed to a link between KLFs and these diseases. The analysis further identifies other potential CACCC-box binding factors, such as SP family members, WT1, VEZF1, and -SALL4, which are implicated in cardiac contraction, remodeling, and inflammation processes.

DNA methylation and chromatin accessibility predict age in the domestic dog

Across mammals, the epigenome is highly predictive of chronological age. These "epigenetic clocks," most of which have been built using DNA methylation (DNAm) profiles, have gained traction as biomarkers of aging and organismal health. While the ability of DNAm to predict chronological age has been repeatedly demonstrated, the ability of other epigenetic features to predict age remains unclear. Here, we use two types of epigenetic information-DNAm, and chromatin accessibility as measured by ATAC-seq-to develop age predictors in peripheral blood mononuclear cells sampled from a population of domesticated dogs. We measured DNAm and ATAC-seq profiles for 71 dogs, building separate predictive clocks from each, as well as the combined dataset. We also use fluorescence-assisted cell sorting to quantify major lymphoid populations for each sample. We found that chromatin accessibility can accurately predict chronological age (R2 ATAC = 26%), though less accurately than the DNAm clock (R2 DNAm = 33%), and the clock built from the combined datasets was comparable to both (R2 combined = 29%). We also observed various populations of CD62L+ T cells significantly correlated with dog age. Finally, we found that all three clocks selected features that were in or near at least two protein-coding genes: BAIAP2 and SCARF2, both previously implicated in processes related to cognitive or neurological impairment. Taken together, these results highlight the potential of chromatin accessibility as a complementary epigenetic resource for modeling and investigating biologic age.

Differential Regulation of Male-Hormones-Related Enhancers Revealed by Chromatin Accessibility and Transcriptional Profiles in Pig Liver

Surgical castration can effectively avoid boar taint and improve pork quality by removing the synthesis of androstenone in the testis, thereby reducing its deposition in adipose tissue. The expression of genes involved in testis-derived hormone metabolism was altered following surgical castration, but the upstream regulatory factors and underlying mechanism remain unclear. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics in liver tissue of castrated and intact full-sibling Yorkshire pigs. First, we identified 897 differentially expressed genes and 6864 differential accessible regions (DARs) using RNA- and ATAC-seq. By integrating the RNA- and ATAC-seq results, 227 genes were identified, and a significant positive correlation was revealed between differential gene expression and the ATAC-seq signal. We constructed a transcription factor regulatory network after motif analysis of DARs and identified a candidate transcription factor (TF) SP1 that targeted the HSD3B1 gene, which was responsible for the metabolism of androstenone. Subsequently, we annotated DARs by incorporating H3K27ac ChIP-seq data, marking 2234 typical enhancers and 245 super enhancers involved in the regulation of all testis-derived hormones. Among these, four typical enhancers associated with HSD3B1 were identified. Furthermore, an in-depth investigation was conducted on the androstenone-related enhancers, and an androstenone-related mutation was identified in a newfound candidatetypical enhancer (andEN) with dual-luciferase assays. These findings provide further insights into how enhancers function as links between phenotypic and non-coding area variations. The discovery of upstream TF and enhancers of HSD3B1 contributes to understanding the regulatory networks of androstenone metabolism and provides an important foundation for improving pork quality.

Cell-fate conversion of intestinal cells in adult Drosophila midgut by depleting a single transcription factor

The manipulation of cell identity by reprograming holds immense potential in regenerative medicine, but is often limited by the inefficient acquisition of fully functional cells. This problem can potentially be resolved by better understanding the reprogramming process using in vivo genetic models, which are currently scarce. Here we report that both enterocytes (ECs) and enteroendocrine cells (EEs) in adult Drosophila midgut show a surprising degree of cell plasticity. Depleting the transcription factor Tramtrack in the differentiated ECs can initiate Prospero-mediated cell transdifferentiation, leading to EE-like cells. On the other hand, depletion of Prospero in the differentiated EEs can lead to the loss of EE-specific transcription programs and the gain of intestinal progenitor cell identity, allowing cell cycle re-entry or differentiation into ECs. We find that intestinal progenitor cells, ECs, and EEs have a similar chromatin accessibility profile, supporting the concept that cell plasticity is enabled by pre-existing chromatin accessibility with switchable transcription programs. Further genetic analysis with this system reveals that the NuRD chromatin remodeling complex, cell lineage confliction, and age act as barriers to EC-to-EE transdifferentiation. The establishment of this genetically tractable in vivo model should facilitate mechanistic investigation of cell plasticity at the molecular and genetic level.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

Advanced Omics Techniques for Understanding Cochlear Genome, Epigenome, and Transcriptome in Health and Disease

Advanced genomics, transcriptomics, and epigenomics techniques are providing unprecedented insights into the understanding of the molecular underpinnings of the central nervous system, including the neuro-sensory cochlea of the inner ear. Here, we report for the first time a comprehensive and updated overview of the most advanced omics techniques for the study of nucleic acids and their applications in cochlear research. We describe the available in vitro and in vivo models for hearing research and the principles of genomics, transcriptomics, and epigenomics, alongside their most advanced technologies (like single-cell omics and spatial omics), which allow for the investigation of the molecular events that occur at a single-cell resolution while retaining the spatial information.

Prediction accuracy of regulatory elements from sequence varies by functional sequencing technique

Introduction

Various sequencing based approaches are used to identify and characterize the activities of cis-regulatory elements in a genome-wide fashion. Some of these techniques rely on indirect markers such as histone modifications (ChIP-seq with histone antibodies) or chromatin accessibility (ATAC-seq, DNase-seq, FAIRE-seq), while other techniques use direct measures such as episomal assays measuring the enhancer properties of DNA sequences (STARR-seq) and direct measurement of the binding of transcription factors (ChIP-seq with transcription factor-specific antibodies). The activities of cis-regulatory elements such as enhancers, promoters, and repressors are determined by their sequence and secondary processes such as chromatin accessibility, DNA methylation, and bound histone markers.

Methods

Here, machine learning models are employed to evaluate the accuracy with which cis-regulatory elements identified by various commonly used sequencing techniques can be predicted by their underlying sequence alone to distinguish between cis-regulatory activity that is reflective of sequence content versus secondary processes.

Results and discussion

Models trained and evaluated on D. melanogaster sequences identified through DNase-seq and STARR-seq are significantly more accurate than models trained on sequences identified by H3K4me1, H3K4me3, and H3K27ac ChIP-seq, FAIRE-seq, and ATAC-seq. These results suggest that the activity detected by DNase-seq and STARR-seq can be largely explained by underlying DNA sequence, independent of secondary processes. Experimentally, a subset of DNase-seq and H3K4me1 ChIP-seq sequences were tested for enhancer activity using luciferase assays and compared with previous tests performed on STARR-seq sequences. The experimental data indicated that STARR-seq sequences are substantially enriched for enhancer-specific activity, while the DNase-seq and H3K4me1 ChIP-seq sequences are not. Taken together, these results indicate that the DNase-seq approach identifies a broad class of regulatory elements of which enhancers are a subset and the associated data are appropriate for training models for detecting regulatory activity from sequence alone, STARR-seq data are best for training enhancer-specific sequence models, and H3K4me1 ChIP-seq data are not well suited for training and evaluating sequence-based models for cis-regulatory element prediction.

Genetic variation in chromatin state across multiple tissues in Drosophila melanogaster

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

Evaluation of MYBL1 as the master regulator for pachytene spermatocyte genes dysregulated in interspecific hybrid dzo

Misregulation of spermatogenesis transcription factors (TF) in hybrids can lead to misexpression, which is a mechanism for hybrid male sterility (HMS). We used dzo (male offspring of Bos taurus ♂ × Bos grunniens ♀) in bovines to investigate the relationship of the key TF with HMS via RNA sequencing and assay for transposase-accessible chromatin with high-throughput sequencing analyses. RNA sequencing showed that the widespread misexpression in dzo was associated with spermatogenesis-related genes and somatic or progenitor genes. The transition from leptotene or zygotene spermatocytes to pachytene spermatocytes may be the key stage for meiosis arrest in dzo. The analysis of TF-binding motif enrichment revealed that the male meiosis-specific master TF MYB proto-oncogene like 1 (MYBL1, known as A-MYB) motif was enriched on the promoters of downregulated pachytene spermatocyte genes in dzo. Assay for transposase-accessible chromatin with high-throughput sequencing revealed that TF-binding sites for MYBL1, nuclear transcription factor Y, and regulatory factor X were enriched in the low-chromatin accessibility region of dzo. The target genes of the MYBL1-binding motif were associated with meiosis-specific genes and significantly downregulated in dzo testis. The transcription factor MYBL1 may be the candidate master regulator for pachytene spermatocyte genes dysregulated in interspecific HMS dzo. This study reported that a few upstream TF regulation changes might exert a cascading effect downstream in a regulatory network as a mechanism for HMS.

Universal NicE-Seq: A Simple and Quick Method for Accessible Chromatin Detection in Fixed Cells

Genome-wide accessible chromatin sequencing and identification has enabled deciphering the epigenetic information encoded in chromatin, revealing accessible promoters, enhancers, nucleosome positioning, transcription factor occupancy, and other chromosomal protein binding. The starting biological materials are often fixed using formaldehyde crosslinking. Here, we describe accessible chromatin library preparation from low numbers of formaldehyde-crosslinked cells using a modified nick translation method, where a nicking enzyme nicks one strand of DNA and DNA polymerase incorporates biotin-conjugated dATP, dCTP, and methyl-dCTP. Once the DNA is labeled, it can be isolated for NGS library preparation. We termed this method as universal NicE-seq (nicking enzyme-assisted sequencing). We also demonstrate a single tube method that enables direct NGS library preparation from low cell numbers without DNA purification. Furthermore, we demonstrated universal NicE-seq on FFPE tissue section sample.

FAIRE-MS reveals mitotic retention of transcriptional regulators on a proteome-wide scale

Mitosis entails global and dramatic alterations, such as higher-order chromatin organization disruption, concomitant with global transcription downregulation. Cells reliably re-establishing gene expression patterns upon mitotic exit and maintaining cellular identities remain poorly understood. Previous studies indicated that certain transcription factors (TFs) remain associated with individual loci during mitosis and serve as mitotic bookmarkers. However, it is unclear which regulatory factors remain bound to the compacted mitotic chromosomes. We developed formaldehyde-assisted isolation of regulatory elements-coupled mass spectrometry (FAIRE-MS) that combines FAIRE-based open chromatin-associated protein pull-down and mass spectrometry (MS) to quantify the open chromatin-associated proteome during the interphase and mitosis. We identified 189 interphase and mitosis maintained (IM) regulatory factors using FAIRE-MS and found intrinsically disordered proteins and regions (IDP(R)s) are highly enriched, which plays a crucial role in liquid-liquid phase separation (LLPS) and chromatin organization during the cell cycle. Notably, in these IDP(R)s, we identified mitotic bookmarkers, such as CEBPB, HMGB1, and TFAP2A, and several factors, including MAX, HMGB3, hnRNP A2/B1, FUS, hnRNP D, and TIAL1, which are at least partially bound to the mitotic chromosome. Furthermore, it will be essential to study whether these IDP(R)s through LLPS helps cells transit from mitosis to the G1 phase during the cell cycle.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

Optimized assay for transposase-accessible chromatin by sequencing (ATAC-seq) library preparation from adult Drosophila melanogaster neurons

Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is rapidly becoming the assay of choice to investigate chromatin-mediated gene regulation, largely because of low input requirements, a fast workflow, and the ability to interrogate the entire genome in an untargeted manner. Many studies using ATAC-seq use mammalian or human-derived tissues, and established protocols work well in these systems. However, ATAC-seq is not yet widely used in Drosophila. Vinegar flies present several advantages over mammalian systems that make them an excellent model for ATAC-seq studies, including abundant genetic tools that allow straightforward targeting, transgene expression, and genetic manipulation that are not available in mammalian models. Because current ATAC-seq protocols are not optimized to use flies, we developed an optimized workflow that accounts for several complicating factors present in Drosophila. We examined parameters affecting nuclei isolation, including input size, freezing time, washing, and possible confounds from retinal pigments. Then, we optimized the enzymatic steps of library construction to account for the smaller Drosophila genome size. Finally, we used our optimized protocol to generate ATAC-seq libraries that meet ENCODE quality metrics. Our optimized protocol enables extensive ATAC-seq experiments in Drosophila, thereby leveraging the advantages of this powerful model system to understand chromatin-mediated gene regulation.

ATAC-Seq Reveals the Landscape of Open Chromatin and cis-Regulatory Elements in the Phytophthora sojae Genome

Nucleosome-free open chromatin often harbors transcription factor (TF)-binding sites that are associated with active cis-regulatory elements. However, analysis of open chromatin regions has rarely been applied to oomycete or fungal plant pathogens. In this study, we performed the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to identify open chromatin and cis-regulatory elements in Phytophthora sojae at the mycelial stage. We identified 10,389 peaks representing nucleosome-free regions (NFRs). The peaks were enriched in gene-promoter regions and associated with 40% of P. sojae genes; transcription levels were higher for genes with multiple peaks than genes with a single peak and were higher for genes with a single peak than genes without peak. Chromatin accessibility was positively correlated with gene transcription level. Through motif discovery based on NFR peaks in core promoter regions, 25 candidate cis-regulatory motifs with evidence of TF-binding footprints were identified. These motifs exhibited various preferences for location in the promoter region and associations with the transcription level of their target genes, which included some putative pathogenicity-related genes. As the first study revealing the landscape of open chromatin and the correlation between chromatin accessibility and gene transcription level in oomycetes, the results provide a technical reference and data resources for future studies on the regulatory mechanisms of gene transcription.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A shared ancient enhancer element differentially regulates the bric-a-brac tandem gene duplicates in the developing Drosophila leg

Gene duplications and transcriptional enhancer emergence/modifications are thought having greatly contributed to phenotypic innovations during animal evolution. Nevertheless, little is known about how enhancers evolve after gene duplication and how regulatory information is rewired between duplicated genes. The Drosophila melanogaster bric-a-brac (bab) complex, comprising the tandem paralogous genes bab1 and bab2, provides a paradigm to address these issues. We previously characterized an intergenic enhancer (named LAE) regulating bab2 expression in the developing legs. We show here that bab2 regulators binding directly the LAE also govern bab1 expression in tarsal cells. LAE excision by CRISPR/Cas9-mediated genome editing reveals that this enhancer appears involved but not strictly required for bab1 and bab2 co-expression in leg tissues. Instead, the LAE enhancer is critical for paralog-specific bab2 expression along the proximo-distal leg axis. Chromatin features and phenotypic rescue experiments indicate that LAE functions partly redundantly with leg-specific regulatory information overlapping the bab1 transcription unit. Phylogenomics analyses indicate that (i) the bab complex originates from duplication of an ancestral singleton gene early on within the Cyclorrhapha dipteran sublineage, and (ii) LAE sequences have been evolutionarily-fixed early on within the Brachycera suborder thus predating the gene duplication event. This work provides new insights on enhancers, particularly about their emergence, maintenance and functional diversification during evolution.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.

Mechanisms underlying the cooperation between loss of epithelial polarity and Notch signaling during neoplastic growth in Drosophila

Aggressive neoplastic growth can be initiated by a limited number of genetic alterations, such as the well-established cooperation between loss of cell architecture and hyperactive signaling pathways. However, our understanding of how these different alterations interact and influence each other remains very incomplete. Using Drosophila paradigms of imaginal wing disc epithelial growth, we have monitored the changes in Notch pathway activity according to the polarity status of cells (scrib mutant). We show that the scrib mutation impacts the direct transcriptional output of the Notch pathway, without altering the global distribution of Su(H), the Notch-dedicated transcription factor. The Notch-dependent neoplasms require, however, the action of a group of transcription factors, similar to those previously identified for Ras/scrib neoplasm (namely AP-1, Stat92E, Ftz-F1 and basic leucine zipper factors), further suggesting the importance of this transcription factor network during neoplastic growth. Finally, our work highlights some Notch/scrib specificities, in particular the role of the PAR domain-containing basic leucine zipper transcription factor and Notch direct target Pdp1 for neoplastic growth.

The Clock:Cycle complex is a major transcriptional regulator of Drosophila photoreceptors that protects the eye from retinal degeneration and oxidative stress

The aging eye experiences physiological changes that include decreased visual function and increased risk of retinal degeneration. Although there are transcriptomic signatures in the aging retina that correlate with these physiological changes, the gene regulatory mechanisms that contribute to cellular homeostasis during aging remain to be determined. Here, we integrated ATAC-seq and RNA-seq data to identify 57 transcription factors that showed differential activity in aging Drosophila photoreceptors. These 57 age-regulated transcription factors include two circadian regulators, Clock and Cycle, that showed sustained increased activity during aging. When we disrupted the Clock:Cycle complex by expressing a dominant negative version of Clock (ClkDN) in adult photoreceptors, we observed changes in expression of 15-20% of genes including key components of the phototransduction machinery and many eye-specific transcription factors. Using ATAC-seq, we showed that expression of ClkDN in photoreceptors leads to changes in activity of 37 transcription factors and causes a progressive decrease in global levels of chromatin accessibility in photoreceptors. Supporting a key role for Clock-dependent transcription in the eye, expression of ClkDN in photoreceptors also induced light-dependent retinal degeneration and increased oxidative stress, independent of light exposure. Together, our data suggests that the circadian regulators Clock and Cycle act as neuroprotective factors in the aging eye by directing gene regulatory networks that maintain expression of the phototransduction machinery and counteract oxidative stress.

Comprehensive epigenomic profiling of human alveolar epithelial differentiation identifies key epigenetic states and transcription factor co-regulatory networks for maintenance of distal lung identity

BACKGROUND

Disruption of alveolar epithelial cell (AEC) differentiation is implicated in distal lung diseases such as chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, and lung adenocarcinoma that impact morbidity and mortality worldwide. Elucidating underlying disease pathogenesis requires a mechanistic molecular understanding of AEC differentiation. Previous studies have focused on changes of individual transcription factors, and to date no study has comprehensively characterized the dynamic, global epigenomic alterations that facilitate this critical differentiation process in humans.

RESULTS

We comprehensively profiled the epigenomic states of human AECs during type 2 to type 1-like cell differentiation, including the methylome and chromatin functional domains, and integrated this with transcriptome-wide RNA expression data. Enhancer regions were drastically altered during AEC differentiation. Transcription factor binding analysis within enhancer regions revealed diverse interactive networks with enrichment for many transcription factors, including NKX2-1 and FOXA family members, as well as transcription factors with less well characterized roles in AEC differentiation, such as members of the MEF2, TEAD, and AP1 families. Additionally, associations among transcription factors changed during differentiation, implicating a complex network of heterotrimeric complex switching in driving differentiation. Integration of AEC enhancer states with the catalog of enhancer elements in the Roadmap Epigenomics Mapping Consortium and Encyclopedia of DNA Elements (ENCODE) revealed that AECs have similar epigenomic structures to other profiled epithelial cell types, including human mammary epithelial cells (HMECs), with NKX2-1 serving as a distinguishing feature of distal lung differentiation.

CONCLUSIONS

Enhancer regions are hotspots of epigenomic alteration that regulate AEC differentiation. Furthermore, the differentiation process is regulated by dynamic networks of transcription factors acting in concert, rather than individually. These findings provide a roadmap for understanding the relationship between disruption of the epigenetic state during AEC differentiation and development of lung diseases that may be therapeutically amenable.

RasV12; scrib-/- Tumors: A Cooperative Oncogenesis Model Fueled by Tumor/Host Interactions

The phenomenon of how oncogenes and tumor-suppressor mutations can synergize to promote tumor fitness and cancer progression can be studied in relatively simple animal model systems such as Drosophila melanogaster. Almost two decades after the landmark discovery of cooperative oncogenesis between oncogenic Ras^V12 and the loss of the tumor suppressor scribble in flies, this and other tumor models have provided new concepts and findings in cancer biology that has remarkable parallels and relevance to human cancer. Here we review findings using the Ras^V12; scrib^-/- tumor model and how it has contributed to our understanding of how these initial simple genetic insults cooperate within the tumor cell to set in motion the malignant transformation program leading to tumor growth through cell growth, cell survival and proliferation, dismantling of cell-cell interactions, degradation of basement membrane and spreading to other organs. Recent findings have demonstrated that cooperativity goes beyond cell intrinsic mechanisms as the tumor interacts with the immediate cells of the microenvironment, the immune system and systemic organs to eventually facilitate malignant progression.

Computational identification and experimental characterization of preferred downstream positions in human core promoters

Metazoan core promoters, which direct the initiation of transcription by RNA polymerase II (Pol II), may contain short sequence motifs termed core promoter elements/motifs (e.g. the TATA box, initiator (Inr) and downstream core promoter element (DPE)), which recruit Pol II via the general transcription machinery. The DPE was discovered and extensively characterized in Drosophila, where it is strictly dependent on both the presence of an Inr and the precise spacing from it. Since the Drosophila DPE is recognized by the human transcription machinery, it is most likely that some human promoters contain a downstream element that is similar, though not necessarily identical, to the Drosophila DPE. However, only a couple of human promoters were shown to contain a functional DPE, and attempts to computationally detect human DPE-containing promoters have mostly been unsuccessful. Using a newly-designed motif discovery strategy based on Expectation-Maximization probabilistic partitioning algorithms, we discovered preferred downstream positions (PDP) in human promoters that resemble the Drosophila DPE. Available chromatin accessibility footprints revealed that Drosophila and human Inr+DPE promoter classes are not only highly structured, but also similar to each other, particularly in the proximal downstream region. Clustering of the corresponding sequence motifs using a neighbor-joining algorithm strongly suggests that canonical Inr+DPE promoters could be common to metazoan species. Using reporter assays we demonstrate the contribution of the identified downstream positions to the function of multiple human promoters. Furthermore, we show that alteration of the spacing between the Inr and PDP by two nucleotides results in reduced promoter activity, suggesting a spacing dependency of the newly discovered human PDP on the Inr. Taken together, our strategy identified novel functional downstream positions within human core promoters, supporting the existence of DPE-like motifs in human promoters.

Transcription of genes by the RNA polymerase II enzyme initiates at a genomic region termed the core promoter. The core promoter is a regulatory region that may contain diverse short DNA sequence motifs/elements that confer specific properties to it. Interestingly, core promoter motifs can be located both upstream and downstream of the transcription start site. Variable compositions of core promoter elements were identified. The initiator (Inr) motif and the downstream core promoter element (DPE) is a combination of elements that has been identified and extensively characterized in fruit flies. Although a few Inr+DPE -containing human promoters were identified, the presence of transcriptionally important downstream core promoter positions within human promoters has been a matter of controversy in the literature. Here, using a newly-designed motif discovery strategy, we discovered preferred downstream positions in human promoters that resemble fruit fly DPE. Clustering of the corresponding sequence motifs in eight additional species indicated that such promoters could be common to multicellular non-plant organisms. Importantly, functional characterization of the newly discovered preferred downstream positions supports the existence of Inr+DPE-containing promoters in human genes.

Dpp and Hedgehog promote the glial response to neuronal apoptosis in the developing Drosophila visual system

Damage in the nervous system induces a stereotypical response that is mediated by glial cells. Here, we use the eye disc of Drosophila melanogaster as a model to explore the mechanisms involved in promoting glial cell response after neuronal cell death induction. We demonstrate that these cells rapidly respond to neuronal apoptosis by increasing in number and undergoing morphological changes, which will ultimately grant them phagocytic abilities. We found that this glial response is controlled by the activity of Decapentaplegic (Dpp) and Hedgehog (Hh) signalling pathways. These pathways are activated after cell death induction, and their functions are necessary to induce glial cell proliferation and migration to the eye discs. The latter of these 2 processes depend on the function of the c-Jun N-terminal kinase (JNK) pathway, which is activated by Dpp signalling. We also present evidence that a similar mechanism controls glial response upon apoptosis induction in the leg discs, suggesting that our results uncover a mechanism that might be involved in controlling glial cells response to neuronal cell death in different regions of the peripheral nervous system (PNS).

In reaction to neuronal damage, glial cells proliferate, change their morphology and alter their behaviour; this response is associated with glial cells’ regenerative function and is necessary to preserve the structural integrity and function of the nervous system. This study identifies a role for the Decapentaplegic and Hedgehog pathways in controlling the glial response.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Deregulation of Transcriptional Enhancers in Cancer

Simple Summary

One of the major challenges in cancer treatments is the dynamic adaptation of tumor cells to cancer therapies. In this regard, tumor cells can modify their response to environmental cues without altering their DNA sequence. This cell plasticity enables cells to undergo morphological and functional changes, for example, during the process of tumour metastasis or when acquiring resistance to cancer therapies. Central to cell plasticity, are the dynamic changes in gene expression that are controlled by a set of molecular switches called enhancers. Enhancers are DNA elements that determine when, where and to what extent genes should be switched on and off. Thus, defects in enhancer function can disrupt the gene expression program and can lead to tumour formation. Here, we review how enhancers control the activity of cancer-associated genes and how defects in these regulatory elements contribute to cell plasticity in cancer. Understanding enhancer (de)regulation can provide new strategies for modulating cell plasticity in tumour cells and can open new research avenues for cancer therapy.

Abstract

Epigenetic regulations can shape a cell’s identity by reversible modifications of the chromatin that ultimately control gene expression in response to internal and external cues. In this review, we first discuss the concept of cell plasticity in cancer, a process that is directly controlled by epigenetic mechanisms, with a particular focus on transcriptional enhancers as the cornerstone of epigenetic regulation. In the second part, we discuss mechanisms of enhancer deregulation in adult stem cells and epithelial-to-mesenchymal transition (EMT), as two paradigms of cell plasticity that are dependent on epigenetic regulation and serve as major sources of tumour heterogeneity. Finally, we review how genetic variations at enhancers and their epigenetic modifiers contribute to tumourigenesis, and we highlight examples of cancer drugs that target epigenetic modifications at enhancers.

Multi-omics analysis of intertumoral heterogeneity within medulloblastoma uncharted-pathway subtypes

Medulloblastoma is a common pediatric malignant brain tumor. There were four consensus molecular subgroups (WNT, SHH, Group3 and Group4). Group 3 and Group 4 tumors exhibited a great degree of transcriptional overlap, and were neither derived from exact pathway aberration. We investigated transcriptional and chromatin accessibility of medulloblastoma by multi-omics single-cell analysis. Our work identified inter- and intra-tumoral heterogeneity within the Group 3, Group 4 and Group 3/4 intermediate subgroups. Unsupervised cluster of each tumor identified 9 cell clusters with transcriptional profiles and 6 cell clusters with chromatin accessibility profiles. OTX2 had the highest activity and expression level across the clusters in a special cluster based on open chromatin single-cell profilings. We identified multiple genes as a significant targeted gene within the OTX2 target genes, which made sense in prognosis. We analyzed the copy-number-variations which presented with expected subgroup distribution from transcriptional and chromatin accessibility profiles. Collectively, these data provide novel insights into Group 3 and Group 4 medulloblastoma and provide a potential therapeutic target.

In vivo tissue-specific chromatin profiling in Drosophila melanogaster using GFP-tagged nuclei

The chromatin landscape defines cellular identity in multicellular organisms with unique patterns of DNA accessibility and histone marks decorating the genome of each cell type. Thus, profiling the chromatin state of different cell types in an intact organism under disease or physiological conditions can provide insight into how chromatin regulates cell homeostasis in vivo. To overcome the many challenges associated with characterizing chromatin state in specific cell types, we developed an improved approach to isolate Drosophila melanogaster nuclei tagged with a GFPKASH protein. The perinuclear space-localized KASH domain anchors GFP to the outer nuclear membrane, and expression of UAS-GFPKASH can be controlled by tissue-specific Gal4 drivers. Using this protocol, we profiled chromatin accessibility using an improved version of Assay for Transposable Accessible Chromatin followed by sequencing (ATAC-seq), called Omni-ATAC. In addition, we examined the distribution of histone marks using Chromatin immunoprecipitation followed by sequencing (ChIP-seq) and Cleavage Under Targets and Tagmentation (CUT&Tag) in adult photoreceptor neurons. We show that the chromatin landscape of photoreceptors reflects the transcriptional state of these cells, demonstrating the quality and reproducibility of our approach for profiling the transcriptome and epigenome of specific cell types in Drosophila.

Single Nucleotide Resolution Analysis Reveals Pervasive, Long-Lasting DNA Methylation Changes by Developmental Exposure to a Mitochondrial Toxicant

Mitochondrial-driven alterations of the epigenome have been reported, but whether they are relevant at the organismal level remains unknown. The viable yellow agouti mouse (A^vy) is a powerful epigenetic biosensor model that reports on the DNA methylation status of the A^vy locus, which is established prior to the three-germ-layer separation, through the coat color of the animals. Here we show that maternal exposure to rotenone, a potent mitochondrial complex I inhibitor, not only changes the DNA methylation status of the A^vy locus in the skin but broadly affects the liver DNA methylome of the offspring. These effects are accompanied by altered gene expression programs that persist throughout life, and which associate with impairment of antioxidant activity and mitochondrial function in aged animals. These pervasive and lasting genomic effects suggest a putative role for mitochondria in regulating life-long gene expression programs through developmental nuclear epigenetic remodeling.

Lozoya et al. provide in vivo evidence of the epigenetic effects of mitochondrial dysfunction. Developmental-only exposure to rotenone through the mother’s diet inhibits mitochondrial complex I in the dams and results in lifelong nuclear DNA methylation and gene expression changes in the offspring. Aged offspring also show functional outcomes.

Evidence of widespread, independent sequence signature for transcription factor cobinding

Transcription factors (TFs) are the vocabulary that genomes use to regulate gene expression and phenotypes. The interactions among TFs enrich this vocabulary and orchestrate diverse biological processes. Although simple models identify open chromatin and the presence of TF motifs as the two major contributors to TF binding patterns, it remains elusive what contributes to the in vivo TF cobinding landscape. In this study, we developed a machine learning algorithm to explore the contributors of the cobinding patterns. The algorithm substantially outperforms the state-of-the-field models for TF cobinding prediction. Game theory-based feature importance analysis reveals that, for most of the TF pairs we studied, independent motif sequences contribute one or more of the two TFs under investigation to their cobinding patterns. Such independent motif sequences include, but are not limited to, transcription initiation-related proteins and known TF complexes. We found the motif sequence signatures and the TFs are rarely mutual, corroborating a hierarchical and directional organization of the regulatory network and refuting the possibility of artifacts caused by shared sequence similarity with the TFs under investigation. We modeled such regulatory language with directed graphs, which reveal shared, global factors that are related to many binding and cobinding patterns.

Chromatin accessibility profiling methods

Chromatin accessibility, or the physical access to chromatinized DNA, is a widely studied characteristic of the eukaryotic genome. As active regulatory DNA elements are generally 'accessible', the genome-wide profiling of chromatin accessibility can be used to identify candidate regulatory genomic regions in a tissue or cell type. Multiple biochemical methods have been developed to profile chromatin accessibility, both in bulk and at the single-cell level. Depending on the method, enzymatic cleavage, transposition or DNA methyltransferases are used, followed by high-throughput sequencing, providing a view of genome-wide chromatin accessibility. In this Primer, we discuss these biochemical methods, as well as bioinformatics tools for analysing and interpreting the generated data, and insights into the key regulators underlying developmental, evolutionary and disease processes. We outline standards for data quality, reproducibility and deposition used by the genomics community. Although chromatin accessibility profiling is invaluable to study gene regulation, alone it provides only a partial view of this complex process. Orthogonal assays facilitate the interpretation of accessible regions with respect to enhancer-promoter proximity, functional transcription factor binding and regulatory function. We envision that technological improvements including single-molecule, multi-omics and spatial methods will bring further insight into the secrets of genome regulation.

The regulatory genome of the malaria vector Anopheles gambiae: integrating chromatin accessibility and gene expression

Anopheles gambiae mosquitoes are primary human malaria vectors, but we know very little about their mechanisms of transcriptional regulation. We profiled chromatin accessibility by the assay for transposase-accessible chromatin by sequencing (ATAC-seq) in laboratory-reared A. gambiae mosquitoes experimentally infected with the human malaria parasite Plasmodium falciparum. By integrating ATAC-seq, RNA-seq and ChIP-seq data, we showed a positive correlation between accessibility at promoters and introns, gene expression and active histone marks. By comparing expression and chromatin structure patterns in different tissues, we were able to infer cis-regulatory elements controlling tissue-specific gene expression and to predict the in vivo binding sites of relevant transcription factors. The ATAC-seq assay also allowed the precise mapping of active regulatory regions, including novel transcription start sites and enhancers that were annotated to mosquito immune-related genes. Not only is this study important for advancing our understanding of mechanisms of transcriptional regulation in the mosquito vector of human malaria, but the information we produced also has great potential for developing new mosquito-control and anti-malaria strategies.

Genomic retargeting of p53 and CTCF is associated with transcriptional changes during oncogenic HRas-induced transformation

Gene transcription is regulated by distant regulatory elements via combinatorial binding of transcription factors. It is increasingly recognized that alterations in chromatin state and transcription factor binding in these distant regulatory elements may have key roles in cancer development. Here we focused on the first stages of oncogene-induced carcinogenic transformation, and characterized the regulatory network underlying transcriptional changes associated with this process. Using Hi-C data, we observe spatial coupling between differentially expressed genes and their differentially accessible regulatory elements and reveal two candidate transcription factors, p53 and CTCF, as determinants of transcriptional alterations at the early stages of oncogenic HRas-induced transformation in human mammary epithelial cells. Strikingly, the malignant transcriptional reprograming is promoted by redistribution of chromatin binding of these factors without major variation in their expression level. Our results demonstrate that alterations in the regulatory landscape have a major role in driving oncogene-induced transcriptional reprogramming.

Lineage-Resolved Enhancer and Promoter Usage during a Time Course of Embryogenesis

Enhancers are essential drivers of cell states, yet the relationship between accessibility, regulatory activity, and in vivo lineage commitment during embryogenesis remains poorly understood. Here, we measure chromatin accessibility in isolated neural and mesodermal lineages across a time course of Drosophila embryogenesis. Promoters, including tissue-specific genes, are often constitutively open, even in contexts where the gene is not expressed. In contrast, the majority of distal elements have dynamic, tissue-specific accessibility. Enhancer priming appears rarely within a lineage, perhaps reflecting the speed of Drosophila embryogenesis. However, many tissue-specific enhancers are accessible in other lineages early on and become progressively closed as embryogenesis proceeds. We demonstrate the usefulness of this tissue- and time-resolved resource to definitively identify single-cell clusters, to uncover predictive motifs, and to identify many regulators of tissue development. For one such predicted neural regulator, l(3)neo38, we generate a loss-of-function mutant and uncover an essential role for neuromuscular junction and brain development.

Repurposing of KLF5 activates a cell cycle signature during the progression from a precursor state to oesophageal adenocarcinoma

Oesophageal adenocarcinoma (OAC) is one of the most common causes of cancer deaths. Barrett's oesophagus (BO) is the only known precancerous precursor to OAC, but our understanding about the molecular events leading to OAC development is limited. Here, we have integrated gene expression and chromatin accessibility profiles of human biopsies and identified a strong cell cycle gene expression signature in OAC compared to BO. Through analysing associated chromatin accessibility changes, we have implicated the transcription factor KLF5 in the transition from BO to OAC. Importantly, we show that KLF5 expression is unchanged during this transition, but instead, KLF5 is redistributed across chromatin to directly regulate cell cycle genes specifically in OAC cells. This new KLF5 target gene programme has potential prognostic significance as high levels correlate with poorer patient survival. Thus, the repurposing of KLF5 for novel regulatory activity in OAC provides new insights into the mechanisms behind disease progression.

Chromatin-modifying elements for recombinant protein production in mammalian cell systems

Mammalian cells are the preferred choice system for the production of complex molecules, such as recombinant therapeutic proteins. Although the technology for increasing the yield of proteins has improved rapidly, the process of selecting, identifying as well as maintaining high-yield cell clones is still troublesome, time-consuming and usually uncertain. Optimization of expression vectors is one of the most effective methods for enhancing protein expression levels. Several commonly used chromatin-modifying elements, including the matrix attachment region, ubiquitous chromatin opening elements, insulators, stabilizing anti-repressor elements can be used to increase the expression level and stability of recombinant proteins. In this review, these chromatin-modifying elements used for the expression vector optimization in mammalian cells are summarized, and future strategies for the utilization of expression cassettes are also discussed.

Total Functional Score of Enhancer Elements Identifies Lineage-Specific Enhancers That Drive Differentiation of Pancreatic Cells

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers.

Rounding up the Usual Suspects: Assessing Yorkie, AP-1, and Stat Coactivation in Tumorigenesis

Can hyperactivation of a few key signaling effectors be the underlying reason for the majority of epithelial cancers despite different driver mutations? Here, to address this question, we use the Drosophila model, which allows analysis of gene expression from tumors with known initiating mutations. Furthermore, its simplified signaling pathways have numerous well characterized targets we can use as pathway readouts. In Drosophila tumor models, changes in the activities of three pathways, Jun N-terminal Kinase (JNK), Janus Kinase / Signal Transducer and Activator of Transcription (JAK/STAT), and Hippo, mediated by AP-1 factors, Stat92E, and Yorkie, are reported frequently. We hypothesized this may indicate that these three pathways are commonly deregulated in tumors. To assess this, we mined the available transcriptomic data and evaluated the activity levels of eight pathways in various tumor models. Indeed, at least two out of our three suspects contribute to tumor development in all Drosophila cancer models assessed, despite different initiating mutations or tissues of origin. Surprisingly, we found that Notch signaling is also globally activated in all models examined. We propose that these four pathways, JNK, JAK/STAT, Hippo, and Notch, are paid special attention and assayed for systematically in existing and newly developed models.

BICORN: An R package for integrative inference of de novo cis-regulatory modules

Genome-wide transcription factor (TF) binding signal analyses reveal co-localization of TF binding sites based on inferred cis-regulatory modules (CRMs). CRMs play a key role in understanding the cooperation of multiple TFs under specific conditions. However, the functions of CRMs and their effects on nearby gene transcription are highly dynamic and context-specific and therefore are challenging to characterize. BICORN (Bayesian Inference of COoperative Regulatory Network) builds a hierarchical Bayesian model and infers context-specific CRMs based on TF-gene binding events and gene expression data for a particular cell type. BICORN automatically searches for a list of candidate CRMs based on the input TF bindings at regulatory regions associated with genes of interest. Applying Gibbs sampling, BICORN iteratively estimates model parameters of CRMs, TF activities, and corresponding regulation on gene transcription, which it models as a sparse network of functional CRMs regulating target genes. The BICORN package is implemented in R (version 3.4 or later) and is publicly available on the CRAN server at https://cran.r-project.org/web/packages/BICORN/index.html.

Identification of genomic enhancers through spatial integration of single-cell transcriptomics and epigenomics

Single-cell technologies allow measuring chromatin accessibility and gene expression in each cell, but jointly utilizing both layers to map bona fide gene regulatory networks and enhancers remains challenging. Here, we generate independent single-cell RNA-seq and single-cell ATAC-seq atlases of the Drosophila eye-antennal disc and spatially integrate the data into a virtual latent space that mimics the organization of the 2D tissue using ScoMAP (Single-Cell Omics Mapping into spatial Axes using Pseudotime ordering). To validate spatially predicted enhancers, we use a large collection of enhancer-reporter lines and identify ~ 85% of enhancers in which chromatin accessibility and enhancer activity are coupled. Next, we infer enhancer-to-gene relationships in the virtual space, finding that genes are mostly regulated by multiple, often redundant, enhancers. Exploiting cell type-specific enhancers, we deconvolute cell type-specific effects of bulk-derived chromatin accessibility QTLs. Finally, we discover that Prospero drives neuronal differentiation through the binding of a GGG motif. In summary, we provide a comprehensive spatial characterization of gene regulation in a 2D tissue.

BAMscale: quantification of next-generation sequencing peaks and generation of scaled coverage tracks

Background

Next-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches.

Results

We have developed BAMscale, a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale, we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~ 100 Gb size) in minutes, outperforming currently available tools.

Conclusions

BAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github (https://github.com/ncbi/BAMscale).

Transposable elements contribute to the genomic response to insecticides in Drosophila melanogaster

Most of the genotype–phenotype analyses to date have largely centred attention on single nucleotide polymorphisms. However, transposable element (TE) insertions have arisen as a plausible addition to the study of the genotypic–phenotypic link because of to their role in genome function and evolution. In this work, we investigate the contribution of TE insertions to the regulation of gene expression in response to insecticides. We exposed four Drosophila melanogaster strains to malathion, a commonly used organophosphate insecticide. By combining information from different approaches, including RNA-seq and ATAC-seq, we found that TEs can contribute to the regulation of gene expression under insecticide exposure by rewiring cis-regulatory networks.

This article is part of a discussion meeting issue ‘Crossroads between transposons and gene regulation’.

Molecular evolution of gland cell types and chemical interactions in animals

Across the Metazoa, the emergence of new ecological interactions has been enabled by the repeated evolution of exocrine glands. Specialized glands have arisen recurrently and with great frequency, even in single genera or species, transforming how animals interact with their environment through trophic resource exploitation, pheromonal communication, chemical defense and parental care. The widespread convergent evolution of animal glands implies that exocrine secretory cells are a hotspot of metazoan cell type innovation. Each evolutionary origin of a novel gland involves a process of 'gland cell type assembly': the stitching together of unique biosynthesis pathways; coordinated changes in secretory systems to enable efficient chemical release; and transcriptional deployment of these machineries into cells constituting the gland. This molecular evolutionary process influences what types of compound a given species is capable of secreting, and, consequently, the kinds of ecological interactions that species can display. Here, we discuss what is known about the evolutionary assembly of gland cell types and propose a framework for how it may happen. We posit the existence of 'terminal selector' transcription factors that program gland function via regulatory recruitment of biosynthetic enzymes and secretory proteins. We suggest ancestral enzymes are initially co-opted into the novel gland, fostering pleiotropic conflict that drives enzyme duplication. This process has yielded the observed pattern of modular, gland-specific biosynthesis pathways optimized for manufacturing specific secretions. We anticipate that single-cell technologies and gene editing methods applicable in diverse species will transform the study of animal chemical interactions, revealing how gland cell types are assembled and functionally configured at a molecular level.

Diverse families of transposable elements affect the transcriptional regulation of stress-response genes in Drosophila melanogaster

Although transposable elements are an important source of regulatory variation, their genome-wide contribution to the transcriptional regulation of stress-response genes has not been studied yet. Stress is a major aspect of natural selection in the wild, leading to changes in the transcriptional regulation of a variety of genes that are often triggered by one or a few transcription factors. In this work, we take advantage of the wealth of information available for Drosophila melanogaster and humans to analyze the role of transposable elements in six stress regulatory networks: immune, hypoxia, oxidative, xenobiotic, heat shock, and heavy metal. We found that transposable elements were enriched for caudal, dorsal, HSF, and tango binding sites in D. melanogaster and for NFE2L2 binding sites in humans. Taking into account the D. melanogaster population frequencies of transposable elements with predicted binding motifs and/or binding sites, we showed that those containing three or more binding motifs/sites are more likely to be functional. For a representative subset of these TEs, we performed in vivo transgenic reporter assays in different stress conditions. Overall, our results showed that TEs are relevant contributors to the transcriptional regulation of stress-response genes.

Widespread activation of developmental gene expression characterized by PRC1-dependent chromatin looping

Polycomb repressive complexes 1 and 2 have been historically described as transcriptional repressors, but recent reports suggest that PRC1 might also support activation, although the underlying mechanisms remain elusive. Here, we show that stage-specific PRC1 binding at a subset of active promoters and enhancers during Drosophila development coincides with the formation of three-dimensional (3D) loops, an increase in expression during development and repression in PRC1 mutants. Dissection of the dachshund locus indicates that PRC1-anchored loops are versatile architectural platforms that persist when surrounding genes are transcriptionally active and fine-tune their expression. The analysis of RING1B binding profiles and 3D contacts during neural differentiation in mice suggests that this role is conserved in mammals.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

Drug-induced chromatin accessibility changes associate with sensitivity to liver tumor promotion

Liver cancer susceptibility varies amongst humans and between experimental animal models because of multiple genetic and epigenetic factors. The molecular characterization of such susceptibilities has the potential to enhance cancer risk assessment of xenobiotic exposures and disease prevention strategies. Here, using DNase I hypersensitivity mapping coupled with transcriptomic profiling, we investigate perturbations in cis-acting gene regulatory elements associated with the early stages of phenobarbital (PB)-mediated liver tumor promotion in susceptible versus resistant mouse strains (B6C3F1 versus C57BL/6J). Integrated computational analyses of strain-selective changes in liver chromatin accessibility underlying PB response reveal differential epigenetic regulation of molecular pathways associated with PB-mediated tumor promotion, including Wnt/β-catenin signaling. Complementary transcription factor motif analyses reveal mouse strain-selective gene regulatory networks and a novel role for Stat, Smad, and Fox transcription factors in the early stages of PB-mediated tumor promotion. Mapping perturbations in cis-acting gene regulatory elements provides novel insights into the molecular basis for susceptibility to xenobiotic-induced rodent liver tumor promotion and has the potential to enhance mechanism-based cancer risk assessments of xenobiotic exposures.

α-Ketoglutarate links p53 to cell fate during tumour suppression

The tumour suppressor TP53 is mutated in the majority of human cancers, and in over 70% of pancreatic ductal adenocarcinoma (PDAC)^1,2. Wild-type p53 accumulates in response to cellular stress, and regulates gene expression to alter cell fate and prevent tumour development². Wild-type p53 is also known to modulate cellular metabolic pathways³, although p53-dependent metabolic alterations that constrain cancer progression remain poorly understood. Here we find that p53 remodels cancer-cell metabolism to enforce changes in chromatin and gene expression that favour a premalignant cell fate. Restoring p53 function in cancer cells derived from KRAS-mutant mouse models of PDAC leads to the accumulation of α-ketoglutarate (αKG, also known as 2-oxoglutarate), a metabolite that also serves as an obligate substrate for a subset of chromatin-modifying enzymes. p53 induces transcriptional programs that are characteristic of premalignant differentiation, and this effect can be partially recapitulated by the addition of cell-permeable αKG. Increased levels of the αKG-dependent chromatin modification 5-hydroxymethylcytosine (5hmC) accompany the tumour-cell differentiation that is triggered by p53, whereas decreased 5hmC characterizes the transition from premalignant to de-differentiated malignant lesions that is associated with mutations in Trp53. Enforcing the accumulation of αKG in p53-deficient PDAC cells through the inhibition of oxoglutarate dehydrogenase-an enzyme of the tricarboxylic acid cycle-specifically results in increased 5hmC, tumour-cell differentiation and decreased tumour-cell fitness. Conversely, increasing the intracellular levels of succinate (a competitive inhibitor of αKG-dependent dioxygenases) blunts p53-driven tumour suppression. These data suggest that αKG is an effector of p53-mediated tumour suppression, and that the accumulation of αKG in p53-deficient tumours can drive tumour-cell differentiation and antagonize malignant progression.

Condition-Specific Modeling of Biophysical Parameters Advances Inference of Regulatory Networks

Large-scale inference of eukaryotic transcription-regulatory networks remains challenging. One underlying reason is that existing algorithms typically ignore crucial regulatory mechanisms, such as RNA degradation and post-transcriptional processing. Here, we describe InfereCLaDR, which incorporates such elements and advances prediction in Saccharomyces cerevisiae. First, InfereCLaDR employs a high-quality Gold Standard dataset that we use separately as prior information and for model validation. Second, InfereCLaDR explicitly models transcription factor activity and RNA half-lives. Third, it introduces expression subspaces to derive condition-responsive regulatory networks for every gene. InfereCLaDR’s final network is validated by known data and trends and results in multiple insights. For example, it predicts long half-lives for transcripts of the nucleic acid metabolism genes and members of the cytosolic chaperonin complex as targets of the proteasome regulator Rpn4p. InfereCLaDR demonstrates that more biophysically realistic modeling of regulatory networks advances prediction accuracy both in eukaryotes and prokaryotes.