Species-aware DNA language models capture regulatory elements and their evolution.

BACKGROUND: The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution. RESULTS: Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery. CONCLUSIONS: Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

p53 rapidly restructures 3D chromatin organization to trigger a transcriptional response.

Activation of the p53 tumor suppressor triggers a transcriptional program to control cellular response to stress. However, the molecular mechanisms by which p53 controls gene transcription are not completely understood. Here, we uncover the critical role of spatio-temporal genome architecture in this process. We demonstrate that p53 drives direct and indirect changes in genome compartments, topologically associating domains, and DNA loops prior to one hour of its activation, which escort the p53 transcriptional program. Focusing on p53-bound enhancers, we report 340 genes directly regulated by p53 over a median distance of 116 kb, with 74% of these genes not previously identified. Finally, we showcase that p53 controls transcription of distal genes through newly formed and pre-existing enhancer-promoter loops in a cohesin dependent manner. Collectively, our findings demonstrate a previously unappreciated architectural role of p53 as regulator at distinct topological layers and provide a reliable set of new p53 direct target genes that may help designs of cancer therapies.

TP63-TRIM29 axis regulates enhancer methylation and chromosomal instability in prostate cancer.

BACKGROUND: Prostate adenocarcinoma (PRAD) is the second leading cause of cancer-related deaths in men. High variability in DNA methylation and a high rate of large genomic rearrangements are often observed in PRAD. RESULTS: To investigate the reasons for such high variance, we integrated DNA methylation, RNA-seq, and copy number alterations datasets from The Cancer Genome Atlas (TCGA), focusing on PRAD, and employed weighted gene co-expression network analysis (WGCNA). Our results show that only single cluster of co-expressed genes is associated with genomic and epigenomic instability. Within this cluster, TP63 and TRIM29 are key transcription regulators and are downregulated in PRAD. We discovered that TP63 regulates the level of enhancer methylation in prostate basal epithelial cells. TRIM29 forms a complex with TP63 and together regulates the expression of genes specific to the prostate basal epithelium. In addition, TRIM29 binds DNA repair proteins and prevents the formation of the TMPRSS2:ERG gene fusion typically observed in PRAD. CONCLUSION: Our study demonstrates that TRIM29 and TP63 are important regulators in maintaining the identity of the basal epithelium under physiological conditions. Furthermore, we uncover the role of TRIM29 in PRAD development.

Dynamic enhancer landscapes in human craniofacial development.

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory.

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Decoding Non-coding Variants: Recent Approaches to Studying Their Role in Gene Regulation and Human Diseases.

Genome-wide association studies (GWAS) have mapped over 90% of disease- and quantitative-trait-associated variants within the non-coding genome. Non-coding regulatory DNA (e.g., promoters and enhancers) and RNA (e.g., 5' and 3' UTRs and splice sites) are essential in regulating temporal and tissue-specific gene expressions. Non-coding variants can potentially impact the phenotype of an organism by altering the molecular recognition of the cis-regulatory elements, leading to gene dysregulation. However, determining causality between non-coding variants, gene regulation, and human disease has remained challenging. Experimental and computational methods have been developed to understand the molecular mechanism involved in non-coding variant interference at the transcriptional and post-transcriptional levels. This review discusses recent approaches to evaluating disease-associated single-nucleotide variants (SNVs) and determines their impact on transcription factor (TF) binding, gene expression, chromatin conformation, post-transcriptional regulation, and translation.

txci-ATAC-seq: a massive-scale single-cell technique to profile chromatin accessibility.

We develop a large-scale single-cell ATAC-seq method by combining Tn5-based pre-indexing with 10× Genomics barcoding, enabling the indexing of up to 200,000 nuclei across multiple samples in a single reaction. We profile 449,953 nuclei across diverse tissues, including the human cortex, mouse brain, human lung, mouse lung, mouse liver, and lung tissue from a club cell secretory protein knockout (CC16-/-) model. Our study of CC16-/- nuclei uncovers previously underappreciated technical artifacts derived from remnant 129 mouse strain genetic material, which cause profound cell-type-specific changes in regulatory elements near many genes, thereby confounding the interpretation of this commonly referenced mouse model.

Exploring transient global transcriptional changes induced by ascorbic acid revealed via atKAS-seq profiling.

Transcription initiates the formation of single-stranded DNA (ssDNA) regions within the genome, delineating transcription bubbles, a highly dynamic genomic process. Kethoxal-assisted single-stranded DNA sequencing (KAS-seq) utilizing N3-kethoxal has emerged as a potent tool for mapping specific guanine positions in ssDNA on a genome-wide scale. However, the original KAS-seq method required the costly Accel-NGS Methyl-seq DNA library kit. This study introduces an optimized iteration of the KAS-seq technique, referred to as adapter-tagged KAS-seq (atKAS-seq), incorporating an adapter tagging strategy. This modification involves integrating sequencing adapters via complementary strand synthesis using random N9 tagging. Additionally, by harnessing the potential of ascorbic acid (ASC), recognized for inducing global epigenetic changes, we employed the atKAS-seq methodology to elucidate critical pathways influenced by short-term, high-dose ASC treatment. Our findings underscore that atKAS-seq enables rapid and precise analyses of transcription dynamics and enhancer activities concurrently. This method offers a streamlined, cost-efficient, and low-input approach, affirming its utility in probing intricate genomic regulatory mechanisms.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution.

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits.

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

The immunoglobulin heavy chain super enhancer controls class switch recombination in developing B cells.

Class switch recombination (CSR) plays an important role in adaptive immune response by enabling mature B cells to replace the initial IgM by another antibody class (IgG, IgE or IgA). CSR is preceded by transcription of the IgH constant genes and is controlled by the super-enhancer 3' regulatory region (3'RR) in an activation-specific manner. The 3'RR is composed of four enhancers (hs3a, hs1-2, hs3b and hs4). In mature B cells, 3'RR activity correlates with transcription of its enhancers. CSR can also occur in primary developing B cells though at low frequency, but in contrast to mature B cells, the transcriptional elements that regulate the process in developing B cells are ill-known. In particular, the role of the 3'RR in the control of constant genes' transcription and CSR has not been addressed. Here, by using a mouse line devoid of the 3'RR and a culture system that highly enriches in pro-B cells, we show that the 3'RR activity is indeed required for switch transcription and CSR, though its effect varies in an isotype-specific manner and correlates with transcription of hs4 enhancer only.

Genome-Wide Investigation of the PLD Gene Family in Tomato: Identification, Analysis, and Expression.

Phospholipase Ds (PLDs) are important phospholipid hydrolases in plants that play crucial roles in the regulation of plant growth, development, and stress tolerance. In this study, 14 PLD genes were identified in the tomato genome and were localized on eight chromosomes, and one tandem-duplicated gene pair was identified. According to a phylogenetic analysis, the genes were categorized into four subtypes: SlPLDα, ß, and δ belonged to the C2-PLD subfamily, while SlPLDζ belonged to the PXPH-PLD subfamily. The gene structure and protein physicochemical properties were highly conserved within the same subtype. The promoter of all the SlPLD genes contained hormone-, light-, and stress-responsive cis-acting regulatory elements, but no significant correlation between the number, distribution, and type of cis-acting elements was observed among the members of the same subtype. Transcriptome data showed that the expression of the SlPLD genes was different in multiple tissues. A quantitative RT-PCR analysis revealed that the SlPLD genes responded positively to cold, salt, drought, and abscisic acid treatments, particularly to salt stress. Different expression patterns were observed for different genes under the same stress, and for the same gene under different stresses. The results provide important insights into the functions of SlPLD genes and lay a foundation for further studies of the response of SlPLD genes to abiotic stresses.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator.

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera).

BACKGROUND: Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm. RESULTS: We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses. CONCLUSION: The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

A SWI/SNF-dependent transcriptional regulation mediated by POU2AF2/C11orf53 at enhancer.

Recent studies have identified a previously uncharacterized protein C11orf53 (now named POU2AF2/OCA-T1), which functions as a robust co-activator of POU2F3, the master transcription factor which is critical for both normal and neoplastic tuft cell identity and viability. Here, we demonstrate that POU2AF2 dictates opposing transcriptional regulation at distal enhance elements. Loss of POU2AF2 leads to an inhibition of active enhancer nearby genes, such as tuft cell identity genes, and a derepression of Polycomb-dependent poised enhancer nearby genes, which are critical for cell viability and differentiation. Mechanistically, depletion of POU2AF2 results in a global redistribution of the chromatin occupancy of the SWI/SNF complex, leading to a significant 3D genome structure change and a subsequent transcriptional reprogramming. Our genome-wide CRISPR screen further demonstrates that POU2AF2 depletion or SWI/SNF inhibition leads to a PTEN-dependent cell growth defect, highlighting a potential role of POU2AF2-SWI/SNF axis in small cell lung cancer (SCLC) pathogenesis. Additionally, pharmacological inhibition of SWI/SNF phenocopies POU2AF2 depletion in terms of gene expression alteration and cell viability decrease in SCLC-P subtype cells. Therefore, impeding POU2AF2-mediated transcriptional regulation represents a potential therapeutic approach for human SCLC therapy.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination.

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

Pathogenic GATA2 genetic variants utilize an obligate enhancer mechanism to distort a multilineage differentiation program.

Mutations in genes encoding transcription factors inactivate or generate ectopic activities to instigate pathogenesis. By disrupting hematopoietic stem/progenitor cells, GATA2 germline variants create a bone marrow failure and leukemia predisposition, GATA2 deficiency syndrome, yet mechanisms underlying the complex phenotypic constellation are unresolved. We used a GATA2-deficient progenitor rescue system to analyze how genetic variation influences GATA2 functions. Pathogenic variants impaired, without abrogating, GATA2-dependent transcriptional regulation. Variants promoted eosinophil and repressed monocytic differentiation without regulating mast cell and erythroid differentiation. While GATA2 and T354M required the DNA-binding C-terminal zinc finger, T354M disproportionately required the N-terminal finger and N terminus. GATA2 and T354M activated a CCAAT/Enhancer Binding Protein-ε (C/EBPε) enhancer, creating a feedforward loop operating with the T-cell Acute Lymphocyte Leukemia-1 (TAL1) transcription factor. Elevating C/EBPε partially normalized hematopoietic defects of GATA2-deficient progenitors. Thus, pathogenic germline variation discriminatively spares or compromises transcription factor attributes, and retaining an obligate enhancer mechanism distorts a multilineage differentiation program.

A 76-base pair duplication within the enhancer region of the HMX1 gene causes sheep microtia.

Sheep congenital microtia is characterized by underdeveloped ears and provides an ideal basis for studying human microtia. This study identified the causal mutation and regulatory mechanisms underlying this disorder. Whole-genome association analysis was conducted using 23 ear tissue samples from sheep with microtia and 28 samples from normal-eared sheep. A significant correlation was found between microtia and a 76-base pair duplication in the enhancer region of the HMX1 gene. Further analysis of offspring phenotypes confirmed an autosomal dominant inheritance pattern. Genotypic analysis showed that individuals that are homozygous for this duplication were earless, heterozygous individuals exhibited shortened ears, and wild-type individuals had normal ears. Moreover, luciferase assays confirmed that this duplication increased HMX1 gene expression, and duplication knock-in mice also exhibited shorter and narrower external ears compared to wild-type mice. Transcriptomic analysis further demonstrated that this duplication enhanced HMX1 gene expression in animal models. This study characterized the causal regulatory mutation underlying sheep microtia.

Cis-regulatory control of mammalian Trps1 gene expression.

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1. To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers. This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Novel Enhanced Mammalian Cell Transient Expression Vector via Promoter Combination.

During the emergence of infectious diseases, evaluating the efficacy of newly developed vaccines requires antigen proteins. Available methods enhance antigen protein productivity; however, structural modifications may occur. Therefore, we aimed to construct a novel transient overexpression vector capable of rapidly producing large quantities of antigenic proteins in mammalian cell lines. This involved expanding beyond the exclusive use of the human cytomegalovirus (CMV) promoter, and was achieved by incorporating a transcriptional enhancer (CMV enhancer), a translational enhancer (woodchuck hepatitis virus post-transcriptional regulatory element), and a promoter based on the CMV promoter. Twenty novel transient expression vectors were constructed, with the vector containing the human elongation factor 1-alpha (EF-1a) promoter showing the highest efficiency in expressing foreign proteins. This vector exhibited an approximately 27-fold higher expression of enhanced green fluorescent protein than the control vector containing only the CMV promoter. It also expressed the highest level of severe acute respiratory syndrome coronavirus 2 receptor-binding domain protein. These observations possibly result from the simultaneous enhancement of the transcriptional activity of the CMV promoter and the human EF-1a promoter by the CMV enhancer. Additionally, the synergistic effect between the CMV and human EF-1a promoters likely contributed to the further enhancement of protein expression.