Functional annotation of the animal genomes: An integrated annotation resource for the horse.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Human Cancer Genomes Display Recurrent Repeat Expansions.

Recurrent repeat expansions were identified near known regulatory elements across seven cancer types.

Recurrent repeat expansions in human cancer genomes.

Expansion of a single repetitive DNA sequence, termed a tandem repeat (TR), is known to cause more than 50 diseases1,2. However, repeat expansions are often not explored beyond neurological and neurodegenerative disorders. In some cancers, mutations accumulate in short tracts of TRs, a phenomenon termed microsatellite instability; however, larger repeat expansions have not been systematically analysed in cancer3-8. Here we identified TR expansions in 2,622 cancer genomes spanning 29 cancer types. In seven cancer types, we found 160 recurrent repeat expansions (rREs), most of which (155/160) were subtype specific. We found that rREs were non-uniformly distributed in the genome with enrichment near candidate cis-regulatory elements, suggesting a potential role in gene regulation. One rRE, a GAAA-repeat expansion, located near a regulatory element in the first intron of UGT2B7 was detected in 34% of renal cell carcinoma samples and was validated by long-read DNA sequencing. Moreover, in preliminary experiments, treating cells that harbour this rRE with a GAAA-targeting molecule led to a dose-dependent decrease in cell proliferation. Overall, our results suggest that rREs may be an important but unexplored source of genetic variation in human cancer, and we provide a comprehensive catalogue for further study.

Mapping Transcription Regulation with Run-on and Sequencing Data Using the Web-Based tfTarget Gateway.

Run-on and sequencing assays like GRO-seq, PRO-seq, and ChRO-seq allow for joint profiling of transcription activity of transcriptional regulatory elements (TREs), i.e., promoters and active enhancers, and target genes. Variation in biological conditions, such as treated vs. control, results in changes in the activity of transcription factors (TFs), which induces concerted changes in TREs and target genes. By modeling the differences between two biological conditions, we developed the computational pipeline known as tfTarget that predicts a set of putative TREs and target genes responding to each TF under the biological condition of interest. In this chapter, we demonstrate the use of the new web-based tfTarget in mapping transcription regulation using run-on sequencing data.

Develop an efficient and specific AAV-based labeling system for Muller glia in mice.

Reprogramming Müller glia (MG) into functional cells is considered a promising therapeutic strategy to treat ocular diseases and vision loss. However, current AAV-based system for MG-tracing was reported to have high leakage in recent studies. Here, we focused on reducing the leakage of AAV-based labeling systems and found that different AAV serotypes showed a range of efficiency and specificity in labeling MG, leading us to optimize a human GFAP-Cre reporter system packaged in the AAV9 serotype with the woodchuck hepatitis virus post-transcriptional regulatory element (WPRE) removed. The leakage ratio of the AAV9-hGFAP-Cre-ΔWPRE decreased by an approximate 40-fold compared with the AAV9-hGFAP-Cre-WPRE labeling system. In addition, we validated the specificity of the AAV-ΔWPRE system for tracing MG reprogramming under Ptbp1-suppression and observed strict non-MG-conversion, similar to previous studies using genetic lineage tracking mouse models. Thus, the AAV9-hGFAP-Cre-ΔWPRE system showed high efficiency and specificity for MG labeling, providing a promising tool for tracing cell fate in vivo.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome.

BACKGROUND: Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome. RESULTS: Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively. CONCLUSIONS: Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

Notch-dependent DNA cis-regulatory elements and their dose-dependent control of C. elegans stem cell self-renewal.

A long-standing biological question is how DNA cis-regulatory elements shape transcriptional patterns during metazoan development. Reporter constructs, cell culture assays and computational modeling have made major contributions to answering this question, but analysis of elements in their natural context is an important complement. Here, we mutate Notch-dependent LAG-1 binding sites (LBSs) in the endogenous Caenorhabditis elegans sygl-1 gene, which encodes a key stem cell regulator, and analyze the consequences on sygl-1 expression (nascent transcripts, mRNA, protein) and stem cell maintenance. Mutation of one LBS in a three-element cluster approximately halved both expression and stem cell pool size, whereas mutation of two LBSs essentially abolished them. Heterozygous LBS mutant clusters provided intermediate values. Our results lead to two major conclusions. First, both LBS number and configuration impact cluster activity: LBSs act additively in trans and synergistically in cis. Second, the SYGL-1 gradient promotes self-renewal above its functional threshold and triggers differentiation below the threshold. Our approach of coupling CRISPR/Cas9 LBS mutations with effects on both molecular and biological readouts establishes a powerful model for in vivo analyses of DNA cis-regulatory elements.

PCRMS: a database of predicted cis-regulatory modules and constituent transcription factor binding sites in genomes.

More accurate and more complete predictions of cis-regulatory modules (CRMs) and constituent transcription factor (TF) binding sites (TFBSs) in genomes can facilitate characterizing functions of regulatory sequences. Here, we developed a database predicted cis-regulatory modules (PCRMS) (https://cci-bioinfo.uncc.edu) that stores highly accurate and unprecedentedly complete maps of predicted CRMs and TFBSs in the human and mouse genomes. The web interface allows the user to browse CRMs and TFBSs in an organism, find the closest CRMs to a gene, search CRMs around a gene and find all TFBSs of a TF. PCRMS can be a useful resource for the research community to characterize regulatory genomes. Database URL: https://cci-bioinfo.uncc.edu/.

Conservation of dynamic characteristics of transcriptional regulatory elements in periodic biological processes.

BACKGROUND: Cell and circadian cycles control a large fraction of cell and organismal physiology by regulating large periodic transcriptional programs that encompass anywhere from 15 to 80% of the genome despite performing distinct functions. In each case, these large periodic transcriptional programs are controlled by gene regulatory networks (GRNs), and it has been shown through genetics and chromosome mapping approaches in model systems that at the core of these GRNs are small sets of genes that drive the transcript dynamics of the GRNs. However, it is unlikely that we have identified all of these core genes, even in model organisms. Moreover, large periodic transcriptional programs controlling a variety of processes certainly exist in important non-model organisms where genetic approaches to identifying networks are expensive, time-consuming, or intractable. Ideally, the core network components could be identified using data-driven approaches on the transcriptome dynamics data already available. RESULTS: This study shows that a unified set of quantified dynamic features of high-throughput time series gene expression data are more prominent in the core transcriptional regulators of cell and circadian cycles than in their outputs, in multiple organism, even in the presence of external periodic stimuli. Additionally, we observe that the power to discriminate between core and non-core genes is largely insensitive to the particular choice of quantification of these features. CONCLUSIONS: There are practical applications of the approach presented in this study for network inference, since the result is a ranking of genes that is enriched for core regulatory elements driving a periodic phenotype. In this way, the method provides a prioritization of follow-up genetic experiments. Furthermore, these findings reveal something unexpected-that there are shared dynamic features of the transcript abundance of core components of unrelated GRNs that control disparate periodic phenotypes.

Impact of Androgen Receptor Activity on Prostate-Specific Membrane Antigen Expression in Prostate Cancer Cells.

Prostate-specific membrane antigen (PSMA) is an essential molecular regulator of prostate cancer (PCa) progression coded by the FOLH1 gene. The PSMA protein has become an important factor in metastatic PCa diagnosis and radioligand therapy. However, low PSMA expression is suggested to be a resistance mechanism to PSMA-based imaging and therapy. Clinical studies revealed that androgen receptor (AR) inhibition increases PSMA expression. The mechanism has not yet been elucidated. Therefore, this study investigated the effect of activation and inhibition of androgen signaling on PSMA expression levels in vitro and compared these findings with PSMA levels in PCa patients receiving systemic therapy. To this end, LAPC4, LNCaP, and C4-2 PCa cells were treated with various concentrations of the synthetic androgen R1881 and antiandrogens. Changes in FOLH1 mRNA were determined using qPCR. Open access databases were used for ChIP-Seq and tissue expression analysis. Changes in PSMA protein were determined using western blot. For PSMA staining in patients' specimens, immunohistochemistry (IHC) was performed. Results revealed that treatment with the synthetic androgen R1881 led to decreased FOLH1 mRNA and PSMA protein. This effect was partially reversed by antiandrogen treatment. However, AR ChIP-Seq analysis revealed no canonical AR binding sites in the regulatory elements of the FOLH1 gene. IHC analysis indicated that androgen deprivation only resulted in increased PSMA expression in patients with low PSMA levels. The data demonstrate that AR activation and inhibition affects PSMA protein levels via a possible non-canonical mechanism. Moreover, analysis of PCa tissue reveals that low PSMA expression rates may be mandatory to increase PSMA by androgen deprivation.

Three classes of epigenomic regulators converge to hyperactivate the essential maternal gene deadhead within a heterochromatin mini-domain.

The formation of a diploid zygote is a highly complex cellular process that is entirely controlled by maternal gene products stored in the egg cytoplasm. This highly specialized transcriptional program is tightly controlled at the chromatin level in the female germline. As an extreme case in point, the massive and specific ovarian expression of the essential thioredoxin Deadhead (DHD) is critically regulated in Drosophila by the histone demethylase Lid and its partner, the histone deacetylase complex Sin3A/Rpd3, via yet unknown mechanisms. Here, we identified Snr1 and Mod(mdg4) as essential for dhd expression and investigated how these epigenomic effectors act with Lid and Sin3A to hyperactivate dhd. Using Cut&Run chromatin profiling with a dedicated data analysis procedure, we found that dhd is intriguingly embedded in an H3K27me3/H3K9me3-enriched mini-domain flanked by DNA regulatory elements, including a dhd promoter-proximal element essential for its expression. Surprisingly, Lid, Sin3a, Snr1 and Mod(mdg4) impact H3K27me3 and this regulatory element in distinct manners. However, we show that these effectors activate dhd independently of H3K27me3/H3K9me3, and that dhd remains silent in the absence of these marks. Together, our study demonstrates an atypical and critical role for chromatin regulators Lid, Sin3A, Snr1 and Mod(mdg4) to trigger tissue-specific hyperactivation within a unique heterochromatin mini-domain.

ReMap 2022: a database of Human, Mouse, Drosophila and Arabidopsis regulatory regions from an integrative analysis of DNA-binding sequencing experiments.

ReMap (https://remap.univ-amu.fr) aims to provide manually curated, high-quality catalogs of regulatory regions resulting from a large-scale integrative analysis of DNA-binding experiments in Human, Mouse, Fly and Arabidopsis thaliana for hundreds of transcription factors and regulators. In this 2022 update, we have uniformly processed >11 000 DNA-binding sequencing datasets from public sources across four species. The updated Human regulatory atlas includes 8103 datasets covering a total of 1210 transcriptional regulators (TRs) with a catalog of 182 million (M) peaks, while the updated Arabidopsis atlas reaches 4.8M peaks, 423 TRs across 694 datasets. Also, this ReMap release is enriched by two new regulatory catalogs for Mus musculus and Drosophila melanogaster. First, the Mouse regulatory catalog consists of 123M peaks across 648 TRs as a result of the integration and validation of 5503 ChIP-seq datasets. Second, the Drosophila melanogaster catalog contains 16.6M peaks across 550 TRs from the integration of 1205 datasets. The four regulatory catalogs are browsable through track hubs at UCSC, Ensembl and NCBI genome browsers. Finally, ReMap 2022 comes with a new Cis Regulatory Module identification method, improved quality controls, faster search results, and better user experience with an interactive tour and video tutorials on browsing and filtering ReMap catalogs.

A cis-regulatory-directed pipeline for the identification of genes involved in cardiac development and disease.

BACKGROUND: Congenital heart diseases are the major cause of death in newborns, but the genetic etiology of this developmental disorder is not fully known. The conventional approach to identify the disease-causing genes focuses on screening genes that display heart-specific expression during development. However, this approach would have discounted genes that are expressed widely in other tissues but may play critical roles in heart development. RESULTS: We report an efficient pipeline of genome-wide gene discovery based on the identification of a cardiac-specific cis-regulatory element signature that points to candidate genes involved in heart development and congenital heart disease. With this pipeline, we retrieve 76% of the known cardiac developmental genes and predict 35 novel genes that previously had no known connectivity to heart development. Functional validation of these novel cardiac genes by RNAi-mediated knockdown of the conserved orthologs in Drosophila cardiac tissue reveals that disrupting the activity of 71% of these genes leads to adult mortality. Among these genes, RpL14, RpS24, and Rpn8 are associated with heart phenotypes. CONCLUSIONS: Our pipeline has enabled the discovery of novel genes with roles in heart development. This workflow, which relies on screening for non-coding cis-regulatory signatures, is amenable for identifying developmental and disease genes for an organ without constraining to genes that are expressed exclusively in the organ of interest.

Transcriptional network orchestrating regional patterning of cortical progenitors.

We uncovered a transcription factor (TF) network that regulates cortical regional patterning in radial glial stem cells. Screening the expression of hundreds of TFs in the developing mouse cortex identified 38 TFs that are expressed in gradients in the ventricular zone (VZ). We tested whether their cortical expression was altered in mutant mice with known patterning defects (Emx2, Nr2f1, and Pax6), which enabled us to define a cortical regionalization TF network (CRTFN). To identify genomic programming underlying this network, we performed TF ChIP-seq and chromatin-looping conformation to identify enhancer-gene interactions. To map enhancers involved in regional patterning of cortical progenitors, we performed assays for epigenomic marks and DNA accessibility in VZ cells purified from wild-type and patterning mutant mice. This integrated approach has identified a CRTFN and VZ enhancers involved in cortical regional patterning in the mouse.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits.

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

Structure of the human SAGA coactivator complex.

The SAGA complex is a regulatory hub involved in gene regulation, chromatin modification, DNA damage repair and signaling. While structures of yeast SAGA (ySAGA) have been reported, there are noteworthy functional and compositional differences for this complex in metazoans. Here we present the cryogenic-electron microscopy (cryo-EM) structure of human SAGA (hSAGA) and show how the arrangement of distinct structural elements results in a globally divergent organization from that of yeast, with a different interface tethering the core module to the TRRAP subunit, resulting in a dramatically altered geometry of functional elements and with the integration of a metazoan-specific splicing module. Our hSAGA structure reveals the presence of an inositol hexakisphosphate (InsP6) binding site in TRRAP and an unusual property of its pseudo-(Ψ)PIKK. Finally, we map human disease mutations, thus providing the needed framework for structure-guided drug design of this important therapeutic target for human developmental diseases and cancer.

Quantitative spatial and temporal assessment of regulatory element activity in zebrafish.

Mutations or genetic variation in noncoding regions of the genome harbouring cis-regulatory elements (CREs), or enhancers, have been widely implicated in human disease and disease risk. However, our ability to assay the impact of these DNA sequence changes on enhancer activity is currently very limited because of the need to assay these elements in an appropriate biological context. Here, we describe a method for simultaneous quantitative assessment of the spatial and temporal activity of wild-type and disease-associated mutant human CRE alleles using live imaging in zebrafish embryonic development. We generated transgenic lines harbouring a dual-CRE dual-reporter cassette in a pre-defined neutral docking site in the zebrafish genome. The activity of each CRE allele is reported via expression of a specific fluorescent reporter, allowing simultaneous visualisation of where and when in development the wild-type allele is active and how this activity is altered by mutation.

A multi-omic single-cell landscape of human gynecologic malignancies.

Deconvolution of regulatory mechanisms that drive transcriptional programs in cancer cells is key to understanding tumor biology. Herein, we present matched transcriptome (scRNA-seq) and chromatin accessibility (scATAC-seq) profiles at single-cell resolution from human ovarian and endometrial tumors processed immediately following surgical resection. This dataset reveals the complex cellular heterogeneity of these tumors and enabled us to quantitatively link variation in chromatin accessibility to gene expression. We show that malignant cells acquire previously unannotated regulatory elements to drive hallmark cancer pathways. Moreover, malignant cells from within the same patients show substantial variation in chromatin accessibility linked to transcriptional output, highlighting the importance of intratumoral heterogeneity. Finally, we infer the malignant cell type-specific activity of transcription factors. By defining the regulatory logic of cancer cells, this work reveals an important reliance on oncogenic regulatory elements and highlights the ability of matched scRNA-seq/scATAC-seq to uncover clinically relevant mechanisms of tumorigenesis in gynecologic cancers.

The Role of Long Non-Coding RNAs in Endometriosis.

Endometriosis is a chronic gynecological disorder affecting the quality of life and fertility of many women around the world. Heterogeneous and non-specific symptoms may lead to a delay in diagnosis, with treatment options limited to surgery and hormonal therapy. Hence, there is a need to better understand the pathogenesis of the disease to improve diagnosis and treatment. Long non-coding RNAs (lncRNAs) have been increasingly shown to be involved in gene regulation but remain relatively under investigated in endometriosis. Mutational and transcriptomic studies have implicated lncRNAs in the pathogenesis of endometriosis. Single-nucleotide polymorphisms (SNPs) in lncRNAs or their regulatory regions have been associated with endometriosis. Genome-wide transcriptomic studies have identified lncRNAs that show deregulated expression in endometriosis, some of which have been subjected to further experiments, which support a role in endometriosis. Mechanistic studies indicate that lncRNAs may regulate genes involved in endometriosis by acting as a molecular sponge for miRNAs, by directly targeting regulatory elements via interactions with chromatin or transcription factors or by affecting signaling pathways. Future studies should concentrate on determining the role of uncharacterized lncRNAs revealed by endometriosis transcriptome studies and the relevance of lncRNAs implicated in the disease by in vitro and animal model studies.