Activity-by-contact model of enhancer-promoter regulation from thousands of CRISPR perturbations

Context-aware single-cell multiomics approach identifies cell-type-specific lung cancer susceptibility genes

Genome-wide association studies (GWAS) identified over fifty loci associated with lung cancer risk. However, underlying mechanisms and target genes are largely unknown, as most risk-associated variants might regulate gene expression in a context-specific manner. Here, we generate a barcode-shared transcriptome and chromatin accessibility map of 117,911 human lung cells from age/sex-matched ever- and never-smokers to profile context-specific gene regulation. Identified candidate cis-regulatory elements (cCREs) are largely cell type-specific, with 37% detected in one cell type. Colocalization of lung cancer candidate causal variants (CCVs) with these cCREs combined with transcription factor footprinting prioritize the variants for 68% of the GWAS loci. CCV-colocalization and trait relevance score indicate that epithelial and immune cell categories, including rare cell types, contribute to lung cancer susceptibility the most. A multi-level cCRE-gene linking system identifies candidate susceptibility genes from 57% of the loci, where most loci display cell-category-specific target genes, suggesting context-specific susceptibility gene function.

Genetic and functional analysis of Raynaud's syndrome implicates loci in vasculature and immunity

Raynaud's syndrome is a dysautonomia where exposure to cold causes vasoconstriction and hypoxia, particularly in the extremities. We performed meta-analysis in four cohorts and discovered eight loci (ADRA2A, IRX1, NOS3, ACVR2A, TMEM51, PCDH10-DT, HLA, and RAB6C) where ADRA2A, ACVR2A, NOS3, TMEM51, and IRX1 co-localized with expression quantitative trait loci (eQTLs), particularly in distal arteries. CRISPR gene editing further showed that ADRA2A and NOS3 loci modified gene expression and in situ RNAscope clarified the specificity of ADRA2A in small vessels and IRX1 around small capillaries in the skin. A functional contraction assay in the cold showed lower contraction in ADRA2A-deficient and higher contraction in ADRA2A-overexpressing smooth muscle cells. Overall, our study highlights the power of genome-wide association testing with functional follow-up as a method to understand complex diseases. The results indicate temperature-dependent adrenergic signaling through ADRA2A, effects at the microvasculature by IRX1, endothelial signaling by NOS3, and immune mechanisms by the HLA locus in Raynaud's syndrome.

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

Extensive co-regulation of neighboring genes complicates the use of eQTLs in target gene prioritization

Identifying causal genes underlying genome-wide association studies (GWASs) is a fundamental problem in human genetics. Although colocalization with gene expression quantitative trait loci (eQTLs) is often used to prioritize GWAS target genes, systematic benchmarking has been limited due to unavailability of large ground truth datasets. Here, we re-analyzed plasma protein QTL data from 3,301 individuals of the INTERVAL cohort together with 131 eQTL Catalog datasets. Focusing on variants located within or close to the affected protein identified 793 proteins with at least one cis-pQTL where we could assume that the most likely causal gene was the gene coding for the protein. We then benchmarked the ability of cis-eQTLs to recover these causal genes by comparing three Bayesian colocalization methods (coloc.susie, coloc.abf, and CLPP) and five Mendelian randomization (MR) approaches (three varieties of inverse-variance weighted MR, MR-RAPS, and MRLocus). We found that assigning fine-mapped pQTLs to their closest protein coding genes outperformed all colocalization methods regarding both precision (71.9%) and recall (76.9%). Furthermore, the colocalization method with the highest recall (coloc.susie - 46.3%) also had the lowest precision (45.1%). Combining evidence from multiple conditionally distinct colocalizing QTLs with MR increased precision to 81%, but this was accompanied by a large reduction in recall to 7.1%. Furthermore, the choice of the MR method greatly affected performance, with the standard inverse-variance-weighted MR often producing many false positives. Our results highlight that linking GWAS variants to target genes remains challenging with eQTL evidence alone, and prioritizing novel targets requires triangulation of evidence from multiple sources.

Super-Enhancer Driven LIF/LIFR-STAT3-SOX2 Regulatory Feedback Loop Promotes Cancer Stemness in Head and Neck Squamous Cell Carcinoma

Super-enhancers (SEs) have been recognized as key epigenetic regulators underlying cancer stemness and malignant traits by aberrant transcriptional control and promising therapeutic targets against human cancers. However, the SE landscape and their roles during head and neck squamous cell carcinoma (HNSCC) development especially in cancer stem cells (CSCs) maintenance remain underexplored yet. Here, we identify leukemia inhibitory factor (LIF)-SE as a representative oncogenic SE to activate LIF transcription in HNSCC. LIF secreted from cancer cells and cancer-associated fibroblasts promotes cancer stemness by driving SOX2 transcription in an autocrine/paracrine manner, respectively. Mechanistically, enhancer elements E1, 2, 4 within LIF-SE recruit SOX2/SMAD3/BRD4/EP300 to facilitate LIF transcription; LIF activates downstream LIFR-STAT3 signaling to drive SOX2 transcription, thus forming a previously unknown regulatory feedback loop (LIF-SE-LIF/LIFR-STAT3-SOX2) to maintain LIF overexpression and CSCs stemness. Clinically, increased LIF abundance in clinical samples correlate with malignant clinicopathological features and patient prognosis; higher LIF concentrations in presurgical plasma dramatically diminish following cancer eradication. Therapeutically, pharmacological targeting LIF-SE-LIF/LIFR-STAT3 significantly impairs tumor growth and reduces CSC subpopulations in xenograft and PDX models. Our findings reveal a hitherto uncharacterized LIF-SE-mediated auto-regulatory loop in regulating HNSCC stemness and highlight LIF as a novel noninvasive biomarker and potential therapeutic target for HNSCC.

Biological Insights from Schizophrenia-associated Loci in Ancestral Populations

Large-scale genome-wide association studies of schizophrenia have uncovered hundreds of associated loci but with extremely limited representation of African diaspora populations. We surveyed electronic health records of 200,000 individuals of African ancestry in the Million Veteran and All of Us Research Programs, and, coupled with genotype-level data from four case-control studies, realized a combined sample size of 13,012 affected and 54,266 unaffected persons. Three genome-wide significant signals - near PLXNA4 , PMAIP1 , and TRPA1 - are the first to be independently identified in populations of predominantly African ancestry. Joint analyses of African, European, and East Asian ancestries across 86,981 cases and 303,771 controls, yielded 376 distinct autosomal loci, which were refined to 708 putatively causal variants via multi-ancestry fine-mapping. Utilizing single-cell functional genomic data from human brain tissue and two complementary approaches, transcriptome-wide association studies and enhancer-promoter contact mapping, we identified a consensus set of 94 genes across ancestries and pinpointed the specific cell types in which they act. We identified reproducible associations of schizophrenia polygenic risk scores with schizophrenia diagnoses and a range of other mental and physical health problems. Our study addresses a longstanding gap in the generalizability of research findings for schizophrenia across ancestral populations, underlining shared biological underpinnings of schizophrenia across global populations in the presence of broadly divergent risk allele frequencies.

The Genetic Determinants and Genomic Consequences of Non-Leukemogenic Somatic Point Mutations

Clonal hematopoiesis (CH) is defined by the expansion of a lineage of genetically identical cells in blood. Genetic lesions that confer a fitness advantage, such as point mutations or mosaic chromosomal alterations (mCAs) in genes associated with hematologic malignancy, are frequent mediators of CH. However, recent analyses of both single cell-derived colonies of hematopoietic cells and population sequencing cohorts have revealed CH frequently occurs in the absence of known driver genetic lesions. To characterize CH without known driver genetic lesions, we used 51,399 deeply sequenced whole genomes from the NHLBI TOPMed sequencing initiative to perform simultaneous germline and somatic mutation analyses among individuals without leukemogenic point mutations (LPM), which we term CH-LPMneg. We quantified CH by estimating the total mutation burden. Because estimating somatic mutation burden without a paired-tissue sample is challenging, we developed a novel statistical method, the Genomic and Epigenomic informed Mutation (GEM) rate, that uses external genomic and epigenomic data sources to distinguish artifactual signals from true somatic mutations. We performed a genome-wide association study of GEM to discover the germline determinants of CH-LPMneg. After fine-mapping and variant-to-gene analyses, we identified seven genes associated with CH-LPMneg ( TCL1A, TERT, SMC4, NRIP1, PRDM16 , MSRA , SCARB1 ), and one locus associated with a sex-associated mutation pathway ( SRGAP2C) . We performed a secondary analysis excluding individuals with mCAs, finding that the genetic architecture was largely unaffected by their inclusion. Functional analyses of SMC4 and NRIP1 implicated altered HSC self-renewal and proliferation as the primary mediator of mutation burden in blood. We then performed comprehensive multi-tissue transcriptomic analyses, finding that the expression levels of 404 genes are associated with GEM. Finally, we performed phenotypic association meta-analyses across four cohorts, finding that GEM is associated with increased white blood cell count and increased risk for incident peripheral artery disease, but is not significantly associated with incident stroke or coronary disease events. Overall, we develop GEM for quantifying mutation burden from WGS without a paired-tissue sample and use GEM to discover the genetic, genomic, and phenotypic correlates of CH-LPMneg.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

Toward a foundation model of causal cell and tissue biology with a Perturbation Cell and Tissue Atlas

Comprehensively charting the biologically causal circuits that govern the phenotypic space of human cells has often been viewed as an insurmountable challenge. However, in the last decade, a suite of interleaved experimental and computational technologies has arisen that is making this fundamental goal increasingly tractable. Pooled CRISPR-based perturbation screens with high-content molecular and/or image-based readouts are now enabling researchers to probe, map, and decipher genetically causal circuits at increasing scale. This scale is now eminently suitable for the deployment of artificial intelligence and machine learning (AI/ML) to both direct further experiments and to predict or generate information that was not-and sometimes cannot-be gathered experimentally. By combining and iterating those through experiments that are designed for inference, we now envision a Perturbation Cell Atlas as a generative causal foundation model to unify human cell biology.

Evaluating Methods for the Prediction of Cell Type-Specific Enhancers in the Mammalian Cortex

Identifying cell type-specific enhancers in the brain is critical to building genetic tools for investigating the mammalian brain. Computational methods for functional enhancer prediction have been proposed and validated in the fruit fly and not yet the mammalian brain. We organized the 'Brain Initiative Cell Census Network (BICCN) Challenge: Predicting Functional Cell Type-Specific Enhancers from Cross-Species Multi-Omics' to assess machine learning and feature-based methods designed to nominate enhancer DNA sequences to target cell types in the mouse cortex. Methods were evaluated based on in vivo validation data from hundreds of cortical cell type-specific enhancers that were previously packaged into individual AAV vectors and retro-orbitally injected into mice. We find that open chromatin was a key predictor of functional enhancers, and sequence models improved prediction of non-functional enhancers that can be deprioritized as opposed to pursued for in vivo testing. Sequence models also identified cell type-specific transcription factor codes that can guide designs of in silico enhancers. This community challenge establishes a benchmark for enhancer prioritization algorithms and reveals computational approaches and molecular information that are crucial for the identification of functional enhancers for mammalian cortical cell types. The results of this challenge bring us closer to understanding the complex gene regulatory landscape of the mammalian brain and help us design more efficient genetic tools and potential gene therapies for human neurological diseases.

IVEA: an integrative variational Bayesian inference method for predicting enhancer-gene regulatory interactions

Motivation

Enhancers play critical roles in cell-type-specific transcriptional control. Despite the identification of thousands of candidate enhancers, unravelling their regulatory relationships with their target genes remains challenging. Therefore, computational approaches are needed to accurately infer enhancer-gene regulatory relationships.

Results

In this study, we propose a new method, IVEA, that predicts enhancer-gene regulatory interactions by estimating promoter and enhancer activities. Its statistical model is based on the gene regulatory mechanism of transcriptional bursting, which is characterized by burst size and frequency controlled by promoters and enhancers, respectively. Using transcriptional readouts, chromatin accessibility, and chromatin contact data as inputs, promoter and enhancer activities were estimated using variational Bayesian inference, and the contribution of each enhancer-promoter pair to target gene transcription was calculated. Our analysis demonstrates that the proposed method can achieve high prediction accuracy and provide biologically relevant enhancer-gene regulatory interactions.

Availability and implementation

The IVEA code is available on GitHub at https://github.com/yasumasak/ivea. The publicly available datasets used in this study are described in Supplementary Table S4.

Hi-C-guided many-polymer model to decipher 3D genome organization

We propose a high-throughput chromosome conformation capture data-based many-polymer model that allows us to generate an ensemble of multi-scale genome structures. We demonstrate the efficacy of our model by validating the generated structures against experimental measurements and employ them to address key questions regarding genome organization. Our model first confirms a significant correlation between chromosome size and nuclear positioning. Specifically, smaller chromosomes are distributed at the core region, whereas larger chromosomes are at the periphery, interacting with the nuclear envelope. The spatial distribution of A- and B-type compartments, which is nontrivial to infer from the corresponding high-throughput chromosome conformation capture maps alone, can also be elucidated using our model, accounting for an issue such as the effect of chromatin-lamina interaction on the compartmentalization of conventional and inverted nuclei. In accordance with imaging data, the overall shape of the 3D genome structures generated from our model displays significant variation. As a case study, we apply our method to the yellow fever mosquito genome, finding that the predicted morphology displays, on average, a more globular shape than the previously suggested spindle-like organization and that our prediction better aligns with the fluorescence in situ hybridization data. Our model has great potential to be extended to investigate many outstanding issues concerning 3D genome organization.

Single-cell analysis of the amphioxus hepatic caecum and vertebrate liver reveals genetic mechanisms of vertebrate liver evolution

The evolution of the vertebrate liver is a prime example of the evolution of complex organs, yet the driving genetic factors behind it remain unknown. Here we study the evolutionary genetics of liver by comparing the amphioxus hepatic caecum and the vertebrate liver, as well as examining the functional transition within vertebrates. Using in vivo and in vitro experiments, single-cell/nucleus RNA-seq data and gene knockout experiments, we confirm that the amphioxus hepatic caecum and vertebrate liver are homologous organs and show that the emergence of ohnologues from two rounds of whole-genome duplications greatly contributed to the functional complexity of the vertebrate liver. Two ohnologues, kdr and flt4, play an important role in the development of liver sinusoidal endothelial cells. In addition, we found that liver-related functions such as coagulation and bile production evolved in a step-by-step manner, with gene duplicates playing a crucial role. We reconstructed the genetic footprint of the transfer of haem detoxification from the liver to the spleen during vertebrate evolution. Together, these findings challenge the previous hypothesis that organ evolution is primarily driven by regulatory elements, underscoring the importance of gene duplicates in the emergence and diversification of a complex organ.

Single-nucleus multiome analysis of human cerebellum in Alzheimer's disease-related dementia

Although human cerebellum is known to be neuropathologically impaired in Alzheimer's disease (AD) and AD-related dementias (ADRD), the cell type-specific transcriptional and epigenomic changes that contribute to this pathology are not well understood. Here, we report single-nucleus multiome (snRNA-seq and snATAC-seq) analysis of 103,861 nuclei isolated from cerebellum from 9 human cases of AD/ADRD and 8 controls, and with frontal cortex of 6 AD donors for additional comparison. Using peak-to-gene linkage analysis, we identified 431,834 significant linkages between gene expression and cell subtype-specific chromatin accessibility regions enriched for candidate cis-regulatory elements (cCREs). These cCREs were associated with AD/ADRD-specific transcriptomic changes and disease-related gene regulatory networks, especially for RAR Related Orphan Receptor A (RORA) and E74 Like ETS Transcription Factor 1 (ELF1) in cerebellar Purkinje cells and granule cells, respectively. Trajectory analysis of granule cell populations further identified disease-relevant transcription factors, such as RORA, and their regulatory targets. Finally, we prioritized two likely causal genes, including Seizure Related 6 Homolog Like 2 (SEZ6L2) in Purkinje cells and KAT8 Regulatory NSL Complex Subunit 1 (KANSL1) in granule cells, through integrative analysis of cCREs derived from snATAC-seq, genome-wide AD/ADRD loci, and Hi-C looping data. This first cell subtype-specific regulatory landscape in the human cerebellum identified here offer novel genomic and epigenomic insights into the neuropathology and pathobiology of AD/ADRD and other neurological disorders if broadly applied.

Genome-wide large-scale multi-trait analysis characterizes global patterns of pleiotropy and unique trait-specific variants

Genome-wide association studies (GWAS) have found widespread evidence of pleiotropy, but characterization of global patterns of pleiotropy remain highly incomplete due to insufficient power of current approaches. We develop fastASSET, a method that allows efficient detection of variant-level pleiotropic association across many traits. We analyze GWAS summary statistics of 116 complex traits of diverse types collected from the GRASP repository and large GWAS Consortia. We identify 2293 independent loci and find that the lead variants in nearly all these loci (~99%) to be associated with ≥ 2 traits (median = 6). We observe that degree of pleiotropy estimated from our study predicts that observed in the UK Biobank for a much larger number of traits (K = 4114) (correlation = 0.43, p-value < 2.2 × 10 - 16 ). Follow-up analyzes of 21 trait-specific variants indicate their link to the expression in trait-related tissues for a small number of genes involved in relevant biological processes. Our findings provide deeper insight into the nature of pleiotropy and leads to identification of highly trait-specific susceptibility variants.

Comprehensive network modeling approaches unravel dynamic enhancer-promoter interactions across neural differentiation

BACKGROUND

Increasing evidence suggests that a substantial proportion of disease-associated mutations occur in enhancers, regions of non-coding DNA essential to gene regulation. Understanding the structures and mechanisms of the regulatory programs this variation affects can shed light on the apparatuses of human diseases.

RESULTS

We collect epigenetic and gene expression datasets from seven early time points during neural differentiation. Focusing on this model system, we construct networks of enhancer-promoter interactions, each at an individual stage of neural induction. These networks serve as the base for a rich series of analyses, through which we demonstrate their temporal dynamics and enrichment for various disease-associated variants. We apply the Girvan-Newman clustering algorithm to these networks to reveal biologically relevant substructures of regulation. Additionally, we demonstrate methods to validate predicted enhancer-promoter interactions using transcription factor overexpression and massively parallel reporter assays.

CONCLUSIONS

Our findings suggest a generalizable framework for exploring gene regulatory programs and their dynamics across developmental processes; this includes a comprehensive approach to studying the effects of disease-associated variation on transcriptional networks. The techniques applied to our networks have been published alongside our findings as a computational tool, E-P-INAnalyzer. Our procedure can be utilized across different cellular contexts and disorders.

Genetic determinants and phenotypic consequences of blood T-cell proportions in 207,000 diverse individuals

T-cells play a critical role in multiple aspects of human health and disease. However, to date the genetic determinants of human T-cell abundance have not been studied at scale because assays quantifying T-cell abundance are not widely used in clinical or research settings. The complete blood count clinical assay quantifies lymphocyte abundance which includes T-cells, B-cells, and NK-cells. To address this gap, we directly estimate T-cell fractions from whole genome sequencing data in over 200,000 individuals from the multi-ethnic TOPMed and All of Us studies. We identified 27 loci associated with T-cell fraction. Interrogating electronic health records identified clinical phenotypes associated with T-cell fraction, including notable changes in T-cell proportions that were highly dynamic over the course of pregnancy. In summary, by estimating T-cell fraction, we obtained new insights into the genetic regulation of T-cells and identified disease consequences of T-cell fractions across the human phenome.

Single cell multiome profiling of pancreatic islets reveals physiological changes in cell type-specific regulation associated with diabetes risk

Physiological variability in pancreatic cell type gene regulation and the impact on diabetes risk is poorly understood. In this study we mapped gene regulation in pancreatic cell types using single cell multiomic (joint RNA-seq and ATAC-seq) profiling in 28 non-diabetic donors in combination with single cell data from 35 non-diabetic donors in the Human Pancreas Analysis Program. We identified widespread associations with age, sex, BMI, and HbA1c, where gene regulatory responses were highly cell type- and phenotype-specific. In beta cells, donor age associated with hypoxia, apoptosis, unfolded protein response, and external signal-dependent transcriptional regulators, while HbA1c associated with inflammatory responses and gender with chromatin organization. We identified 10.8K loci where genetic variants were QTLs for cis regulatory element (cRE) accessibility, including 20% with lineage- or cell type-specific effects which disrupted distinct transcription factor motifs. Type 2 diabetes and glycemic trait associated variants were enriched in both phenotype- and QTL-associated beta cell cREs, whereas type 1 diabetes showed limited enrichment. Variants at 226 diabetes and glycemic trait loci were QTLs in beta and other cell types, including 40 that were statistically colocalized, and annotating target genes of colocalized QTLs revealed genes with putatively novel roles in disease. Our findings reveal diverse responses of pancreatic cell types to phenotype and genotype in physiology, and identify pathways, networks, and genes through which physiology impacts diabetes risk.

Identification, Design, and Application of Noncoding Cis-Regulatory Elements

Cis-regulatory elements (CREs) play a pivotal role in orchestrating interactions with trans-regulatory factors such as transcription factors, RNA-binding proteins, and noncoding RNAs. These interactions are fundamental to the molecular architecture underpinning complex and diverse biological functions in living organisms, facilitating a myriad of sophisticated and dynamic processes. The rapid advancement in the identification and characterization of these regulatory elements has been marked by initiatives such as the Encyclopedia of DNA Elements (ENCODE) project, which represents a significant milestone in the field. Concurrently, the development of CRE detection technologies, exemplified by massively parallel reporter assays, has progressed at an impressive pace, providing powerful tools for CRE discovery. The exponential growth of multimodal functional genomic data has necessitated the application of advanced analytical methods. Deep learning algorithms, particularly large language models, have emerged as invaluable tools for deconstructing the intricate nucleotide sequences governing CRE function. These advancements facilitate precise predictions of CRE activity and enable the de novo design of CREs. A deeper understanding of CRE operational dynamics is crucial for harnessing their versatile regulatory properties. Such insights are instrumental in refining gene therapy techniques, enhancing the efficacy of selective breeding programs, pushing the boundaries of genetic innovation, and opening new possibilities in microbial synthetic biology.

EPInformer: a scalable deep learning framework for gene expression prediction by integrating promoter-enhancer sequences with multimodal epigenomic data

Transcriptional regulation, critical for cellular differentiation and adaptation to environmental changes, involves coordinated interactions among DNA sequences, regulatory proteins, and chromatin architecture. Despite extensive data from consortia like ENCODE, understanding the dynamics of cis-regulatory elements (CREs) in gene expression remains challenging. Deep learning is a powerful tool for learning gene expression and epigenomic signals from DNA sequences, exhibiting superior performance compared to conventional machine learning approaches. However, even the most advanced deep learning-based methods may fall short in capturing the regulatory effects of distal elements such as enhancers, limiting their predictive accuracy. In addition, these methods may require significant resources to train or to adapt to newly generated data. To address these challenges, we present EPInformer, a scalable deep-learning framework for predicting gene expression by integrating promoter-enhancer interactions with their sequences, epigenomic signals, and chromatin contacts. Our model outperforms existing gene expression prediction models in rigorous cross-chromosome validation, accurately recapitulates enhancer-gene interactions validated by CRISPR perturbation experiments, and identifies crucial transcription factor motifs within regulatory sequences. EPInformer is available as open-source software at https://github.com/pinellolab/EPInformer.

Epigenome editing technologies for discovery and medicine

Epigenome editing has rapidly evolved in recent years, with diverse applications that include elucidating gene regulation mechanisms, annotating coding and noncoding genome functions and programming cell state and lineage specification. Importantly, given the ubiquitous role of epigenetics in complex phenotypes, epigenome editing has unique potential to impact a broad spectrum of diseases. By leveraging powerful DNA-targeting technologies, such as CRISPR, epigenome editing exploits the heritable and reversible mechanisms of epigenetics to alter gene expression without introducing DNA breaks, inducing DNA damage or relying on DNA repair pathways.

Scalable genetic screening for regulatory circuits using compressed Perturb-seq

Pooled CRISPR screens with single-cell RNA sequencing readout (Perturb-seq) have emerged as a key technique in functional genomics, but they are limited in scale by cost and combinatorial complexity. In this study, we modified the design of Perturb-seq by incorporating algorithms applied to random, low-dimensional observations. Compressed Perturb-seq measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq with an order of magnitude cost reduction and greater power to learn genetic interactions. We identified known and novel regulators of immune responses and uncovered evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing genome-wide association studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Simplifying causal gene identification in GWAS loci

Genome-wide association studies (GWAS) help to identify disease-linked genetic variants, but pinpointing the most likely causal genes in GWAS loci remains challenging. Existing GWAS gene prioritization tools are powerful, but often use complex black box models trained on datasets containing unaddressed biases. Here we present CALDERA, a gene prioritization tool that achieves similar or better performance than state-of-the-art methods, but uses just 12 features and a simple logistic regression model with L1 regularization. We use a data-driven approach to construct a truth set of causal genes in 406 GWAS loci and correct for potential confounders. We demonstrate that CALDERA is well-calibrated in external datasets and prioritizes genes with expected properties, such as being mutation-intolerant (OR = 1.751 for pLI > 90%, P = 8.45×10-3). CALDERA facilitates the prioritization of potentially causal genes in GWAS loci and may help identify novel genetics-driven drug targets.

CatLearning: highly accurate gene expression prediction from histone mark

Histone modifications, known as histone marks, are pivotal in regulating gene expression within cells. The vast array of potential combinations of histone marks presents a considerable challenge in decoding the regulatory mechanisms solely through biological experimental approaches. To overcome this challenge, we have developed a method called CatLearning. It utilizes a modified convolutional neural network architecture with a specialized adaptation Residual Network to quantitatively interpret histone marks and predict gene expression. This architecture integrates long-range histone information up to 500Kb and learns chromatin interaction features without 3D information. By using only one histone mark, CatLearning achieves a high level of accuracy. Furthermore, CatLearning predicts gene expression by simulating changes in histone modifications at enhancers and throughout the genome. These findings help comprehend the architecture of histone marks and develop diagnostic and therapeutic targets for diseases with epigenetic changes.

SnapHiC-G: identifying long-range enhancer-promoter interactions from single-cell Hi-C data via a global background model

Harnessing the power of single-cell genomics technologies, single-cell Hi-C (scHi-C) and its derived technologies provide powerful tools to measure spatial proximity between regulatory elements and their target genes in individual cells. Using a global background model, we propose SnapHiC-G, a computational method, to identify long-range enhancer-promoter interactions from scHi-C data. We applied SnapHiC-G to scHi-C datasets generated from mouse embryonic stem cells and human brain cortical cells. SnapHiC-G achieved high sensitivity in identifying long-range enhancer-promoter interactions. Moreover, SnapHiC-G can identify putative target genes for noncoding genome-wide association study (GWAS) variants, and the genetic heritability of neuropsychiatric diseases is enriched for single-nucleotide polymorphisms (SNPs) within SnapHiC-G-identified interactions in a cell-type-specific manner. In sum, SnapHiC-G is a powerful tool for characterizing cell-type-specific enhancer-promoter interactions from complex tissues and can facilitate the discovery of chromatin interactions important for gene regulation in biologically relevant cell types.

Epigenetic control of metabolic identity across cell types

Constraint-based network modelling is a powerful tool for analysing cellular metabolism at genomic scale. Here, we conducted an integrative analysis of metabolic networks reconstructed from RNA-seq data with paired epigenomic data from the EpiATLAS resource of the International Human Epigenome Consortium (IHEC). Applying a state-of-the-art contextualisation algorithm, we reconstructed metabolic networks across 1,555 samples corresponding to 58 tissues and cell types. Analysis of these networks revealed the distribution of metabolic functionalities across human cell types and provides a compendium of human metabolic activity. This integrative approach allowed us to define, across tissues and cell types, i) reactions that fulfil the basic metabolic processes (core metabolism), and ii) cell type-specific functions (unique metabolism), that shape the metabolic identity of a cell or a tissue. Integration with EpiATLAS-derived cell type-specific gene-level chromatin states and enhancer-gene interactions identified enhancers, transcription factors, and key nodes controlling core and unique metabolism. Transport and first reactions of pathways were enriched for high expression, active chromatin state, and Polycomb-mediated repression in cell types where pathways are inactive, suggesting that key nodes are targets of repression. This integrative analysis forms the basis for identifying regulation points that control metabolic identity in human cells.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.

Diversity and scale: Genetic architecture of 2068 traits in the VA Million Veteran Program

One of the justifiable criticisms of human genetic studies is the underrepresentation of participants from diverse populations. Lack of inclusion must be addressed at-scale to identify causal disease factors and understand the genetic causes of health disparities. We present genome-wide associations for 2068 traits from 635,969 participants in the Department of Veterans Affairs Million Veteran Program, a longitudinal study of diverse United States Veterans. Systematic analysis revealed 13,672 genomic risk loci; 1608 were only significant after including non-European populations. Fine-mapping identified causal variants at 6318 signals across 613 traits. One-third (n = 2069) were identified in participants from non-European populations. This reveals a broadly similar genetic architecture across populations, highlights genetic insights gained from underrepresented groups, and presents an extensive atlas of genetic associations.

Long range regulation of transcription scales with genomic distance in a gene specific manner

While critical for tuning the timing and level of transcription, enhancer communication with distal promoters is not well understood. Here we bypass the need for sequence-specific transcription factors and recruit activators directly using CARGO-VPR, an approach for targeting dCas9-VPR using a multiplexed array of RNA guides. We show that this approach achieves effective activator recruitment to arbitrary genomic sites, even those inaccessible by single dCas9. We utilize CARGO-VPR across the Prdm8-Fgf5 locus in mESCs, where neither gene is expressed. We demonstrate that while activator recruitment to any tested region results in transcriptional induction of at least one gene, the expression level strongly depends on the genomic distance between the promoter and activator recruitment site. However, the expression-distance relationship for each gene scales distinctly in a manner not attributable to differences in 3D contact frequency, promoter DNA sequence or presence of the repressive chromatin marks at the locus.

Endogenous retroviruses mediate transcriptional rewiring in response to oncogenic signaling in colorectal cancer

Cancer cells exhibit rewired transcriptional regulatory networks that promote tumor growth and survival. However, the mechanisms underlying the formation of these pathological networks remain poorly understood. Through a pan-cancer epigenomic analysis, we found that primate-specific endogenous retroviruses (ERVs) are a rich source of enhancers displaying cancer-specific activity. In colorectal cancer and other epithelial tumors, oncogenic MAPK/AP1 signaling drives the activation of enhancers derived from the primate-specific ERV family LTR10. Functional studies in colorectal cancer cells revealed that LTR10 elements regulate tumor-specific expression of multiple genes associated with tumorigenesis, such as ATG12 and XRCC4. Within the human population, individual LTR10 elements exhibit germline and somatic structural variation resulting from a highly mutable internal tandem repeat region, which affects AP1 binding activity. Our findings reveal that ERV-derived enhancers contribute to transcriptional dysregulation in response to oncogenic signaling and shape the evolution of cancer-specific regulatory networks.

Identifying dysregulated regions in amyotrophic lateral sclerosis through chromatin accessibility outliers

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

Rhabdomyosarcoma fusion oncoprotein initially pioneers a neural signature in vivo

Fusion-positive rhabdomyosarcoma is an aggressive pediatric cancer molecularly characterized by arrested myogenesis. The defining genetic driver, PAX3::FOXO1, functions as a chimeric gain-of-function transcription factor. An incomplete understanding of PAX3::FOXO1's in vivo epigenetic mechanisms has hindered therapeutic development. Here, we establish a PAX3::FOXO1 zebrafish injection model and semi-automated ChIP-seq normalization strategy to evaluate how PAX3::FOXO1 initially interfaces with chromatin in a developmental context. We investigated PAX3::FOXO1's recognition of chromatin and subsequent transcriptional consequences. We find that PAX3::FOXO1 interacts with inaccessible chromatin through partial/homeobox motif recognition consistent with pioneering activity. However, PAX3::FOXO1-genome binding through a composite paired-box/homeobox motif alters chromatin accessibility and redistributes H3K27ac to activate neural transcriptional programs. We uncover neural signatures that are highly representative of clinical rhabdomyosarcoma gene expression programs that are enriched following chemotherapy. Overall, we identify partial/homeobox motif recognition as a new mode for PAX3::FOXO1 pioneer function and identify neural signatures as a potentially critical PAX3::FOXO1 tumor initiation event.

CRISPR-Cas9-mediated deletion enhancer of MECOM play a tumor suppressor role in ovarian cancer

MDS1 and EVI1 complex locus (MECOM), a transcription factor encoding several variants, has been implicated in progression of ovarian cancer. The function of regulatory regions in regulating MECOM expression in ovarian cancer is not fully understood. In this study, MECOM expression was evaluated in ovarian cancer cell lines treated with bromodomain and extraterminal (BET) inhibitor JQ-1. Oncogenic phenotypes were assayed using assays of CCK-8, colony formation, wound-healing and transwell. Oncogenic phenotypes were estimated in stable sgRNA-transfected OVCAR3 cell lines. Xenograft mouse model was assayed via subcutaneous injection of enhancer-deleted OVCAR3 cell lines. The results displayed that expression of MECOM is downregulated in cell lines treated with JQ-1. Data from published ChIP-sequencing (H3K27Ac) in 3 ovarian cancer cell lines displayed a potential enhancer around the first exon. mRNA and protein expression were downregulated in OVCAR3 cells after deletion of the MECOM enhancer. Similarly, oncogenic phenotypes both in cells and in the xenograft mouse model were significantly attenuated. This study demonstrates that JQ-1 can inhibit the expression of MECOM and tumorigenesis. Deletion of the enhancer activity of MECOM has an indispensable role in inhibiting ovarian cancer progress, which sheds light on a promising opportunity for ovarian cancer treatment through the application of this non-coding DNA deletion.

High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility.

Single cell transcriptomes and multiscale networks from persons with and without Alzheimer's disease

The emergence of single nucleus RNA sequencing (snRNA-seq) offers to revolutionize the study of Alzheimer's disease (AD). Integration with complementary multiomics data such as genetics, proteomics and clinical data provides powerful opportunities to link cell subpopulations and molecular networks with a broader disease-relevant context. We report snRNA-seq profiles from superior frontal gyrus samples from 101 well characterized subjects from the Banner Brain and Body Donation Program in combination with whole genome sequences. We report findings that link common AD risk variants with CR1 expression in oligodendrocytes as well as alterations in hematological parameters. We observed an AD-associated CD83(+) microglial subtype with unique molecular networks and which is associated with immunoglobulin IgG4 production in the transverse colon. Our major observations were replicated in two additional, independent snRNA-seq data sets. These findings illustrate the power of multi-tissue molecular profiling to contextualize snRNA-seq brain transcriptomics and reveal disease biology.

Dynamic stress- and inflammatory-based regulation of psychiatric risk loci in human neurons

The prenatal environment can alter neurodevelopmental and clinical trajectories, markedly increasing risk for psychiatric disorders in childhood and adolescence. To understand if and how fetal exposures to stress and inflammation exacerbate manifestation of genetic risk for complex brain disorders, we report a large-scale context-dependent massively parallel reporter assay (MPRA) in human neurons designed to catalogue genotype x environment (GxE) interactions. Across 240 genome-wide association study (GWAS) loci linked to ten brain traits/disorders, the impact of hydrocortisone, interleukin 6, and interferon alpha on transcriptional activity is empirically evaluated in human induced pluripotent stem cell (hiPSC)-derived glutamatergic neurons. Of ~3,500 candidate regulatory risk elements (CREs), 11% of variants are active at baseline, whereas cue-specific CRE regulatory activity range from a high of 23% (hydrocortisone) to a low of 6% (IL-6). Cue-specific regulatory activity is driven, at least in part, by differences in transcription factor binding activity, the gene targets of which show unique enrichments for brain disorders as well as co-morbid metabolic and immune syndromes. The dynamic nature of genetic regulation informs the influence of environmental factors, reveals a mechanism underlying pleiotropy and variable penetrance, and identifies specific risk variants that confer greater disorder susceptibility after exposure to stress or inflammation. Understanding neurodevelopmental GxE interactions will inform mental health trajectories and uncover novel targets for therapeutic intervention.

HyperCas12a enables highly-multiplexed epigenome editing screens

Interactions between multiple genes or cis-regulatory elements (CREs) underlie a wide range of biological processes in both health and disease. High-throughput screens using dCas9 fused to epigenome editing domains have allowed researchers to assess the impact of activation or repression of both coding and non-coding genomic regions on a phenotype of interest, but assessment of genetic interactions between those elements has been limited to pairs. Here, we combine a hyper-efficient version of Lachnospiraceae bacterium dCas12a (dHyperLbCas12a) with RNA Polymerase II expression of long CRISPR RNA (crRNA) arrays to enable efficient highly-multiplexed epigenome editing. We demonstrate that this system is compatible with several activation and repression domains, including the P300 histone acetyltransferase domain and SIN3A interacting domain (SID). We also show that the dCas12a platform can perform simultaneous activation and repression using a single crRNA array via co-expression of multiple dCas12a orthologues. Lastly, demonstrate that the dCas12a system is highly effective for high-throughput screens. We use dHyperLbCas12a-KRAB and a ~19,000-member barcoded library of crRNA arrays containing six crRNAs each to dissect the independent and combinatorial contributions of CREs to the dose-dependent control of gene expression at a glucocorticoid-responsive locus. The tools and methods introduced here create new possibilities for highly multiplexed control of gene expression in a wide variety of biological systems.

The contribution of silencer variants to human diseases

BACKGROUND

Although disease-causal genetic variants have been found within silencer sequences, we still lack a comprehensive analysis of the association of silencers with diseases. Here, we profiled GWAS variants in 2.8 million candidate silencers across 97 human samples derived from a diverse panel of tissues and developmental time points, using deep learning models.

RESULTS

We show that candidate silencers exhibit strong enrichment in disease-associated variants, and several diseases display a much stronger association with silencer variants than enhancer variants. Close to 52% of candidate silencers cluster, forming silencer-rich loci, and, in the loci of Parkinson's-disease-hallmark genes TRIM31 and MAL, the associated SNPs densely populate clustered candidate silencers rather than enhancers displaying an overall twofold enrichment in silencers versus enhancers. The disruption of apoptosis in neuronal cells is associated with both schizophrenia and bipolar disorder and can largely be attributed to variants within candidate silencers. Our model permits a mechanistic explanation of causative SNP effects by identifying altered binding of tissue-specific repressors and activators, validated with a 70% of directional concordance using SNP-SELEX. Narrowing the focus of the analysis to individual silencer variants, experimental data confirms the role of the rs62055708 SNP in Parkinson's disease, rs2535629 in schizophrenia, and rs6207121 in type 1 diabetes.

CONCLUSIONS

In summary, our results indicate that advances in deep learning models for the discovery of disease-causal variants within candidate silencers effectively "double" the number of functionally characterized GWAS variants. This provides a basis for explaining mechanisms of action and designing novel diagnostics and therapeutics.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

GET: a foundation model of transcription across human cell types

Transcriptional regulation, involving the complex interplay between regulatory sequences and proteins, directs all biological processes. Computational models of transcription lack generalizability to accurately extrapolate in unseen cell types and conditions. Here, we introduce GET, an interpretable foundation model designed to uncover regulatory grammars across 213 human fetal and adult cell types. Relying exclusively on chromatin accessibility data and sequence information, GET achieves experimental-level accuracy in predicting gene expression even in previously unseen cell types. GET showcases remarkable adaptability across new sequencing platforms and assays, enabling regulatory inference across a broad range of cell types and conditions, and uncovering universal and cell type specific transcription factor interaction networks. We evaluated its performance on prediction of regulatory activity, inference of regulatory elements and regulators, and identification of physical interactions between transcription factors. Specifically, we show GET outperforms current models in predicting lentivirus-based massive parallel reporter assay readout with reduced input data. In fetal erythroblasts, we identify distal (>1Mbp) regulatory regions that were missed by previous models. In B cells, we identified a lymphocyte-specific transcription factor-transcription factor interaction that explains the functional significance of a leukemia-risk predisposing germline mutation. In sum, we provide a generalizable and accurate model for transcription together with catalogs of gene regulation and transcription factor interactions, all with cell type specificity.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remains unclear. Here, we systematically identified enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We found that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further showed that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

Genome-wide Cas9-mediated screening of essential non-coding regulatory elements via libraries of paired single-guide RNAs

The functions of non-coding regulatory elements (NCREs), which constitute a major fraction of the human genome, have not been systematically studied. Here we report a method involving libraries of paired single-guide RNAs targeting both ends of an NCRE as a screening system for the Cas9-mediated deletion of thousands of NCREs genome-wide to study their functions in distinct biological contexts. By using K562 and 293T cell lines and human embryonic stem cells, we show that NCREs can have redundant functions, and that many ultra-conserved elements have silencer activity and play essential roles in cell growth and in cellular responses to drugs (notably, the ultra-conserved element PAX6_Tarzan may be critical for heart development, as removing it from human embryonic stem cells led to defects in cardiomyocyte differentiation). The high-throughput screen, which is compatible with single-cell sequencing, may allow for the identification of druggable NCREs.

Enhancer selectivity in space and time: from enhancer-promoter interactions to promoter activation

The primary regulators of metazoan gene expression are enhancers, originally functionally defined as DNA sequences that can activate transcription at promoters in an orientation-independent and distance-independent manner. Despite being crucial for gene regulation in animals, what mechanisms underlie enhancer selectivity for promoters, and more fundamentally, how enhancers interact with promoters and activate transcription, remain poorly understood. In this Review, we first discuss current models of enhancer-promoter interactions in space and time and how enhancers affect transcription activation. Next, we discuss different mechanisms that mediate enhancer selectivity, including repression, biochemical compatibility and regulation of 3D genome structure. Through 3D polymer simulations, we illustrate how the ability of 3D genome folding mechanisms to mediate enhancer selectivity strongly varies for different enhancer-promoter interaction mechanisms. Finally, we discuss how recent technical advances may provide new insights into mechanisms of enhancer-promoter interactions and how technical biases in methods such as Hi-C and Micro-C and imaging techniques may affect their interpretation.

Discovering mechanisms of human genetic variation and controlling cell states at scale

Population-scale sequencing efforts have catalogued substantial genetic variation in humans such that variant discovery dramatically outpaces interpretation. We discuss how single-cell sequencing is poised to reveal genetic mechanisms at a rate that may soon approach that of variant discovery. The functional genomics toolkit is sufficiently modular to systematically profile almost any type of variation within increasingly diverse contexts and with molecularly comprehensive and unbiased readouts. As a result, we can construct deep phenotypic atlases of variant effects that span the entire regulatory cascade. The same conceptual approach to interpreting genetic variation should be applied to engineering therapeutic cell states. In this way, variant mechanism discovery and cell state engineering will become reciprocating and iterative processes towards genomic medicine.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

Five autism-associated transcriptional regulators target shared loci proximal to brain-expressed genes

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

Fine-mapping causal tissues and genes at disease-associated loci

Heritable diseases often manifest in a highly tissue-specific manner, with different disease loci mediated by genes in distinct tissues or cell types. We propose Tissue-Gene Fine-Mapping (TGFM), a fine-mapping method that infers the posterior probability (PIP) for each gene-tissue pair to mediate a disease locus by analyzing GWAS summary statistics (and in-sample LD) and leveraging eQTL data from diverse tissues to build cis-predicted expression models; TGFM also assigns PIPs to causal variants that are not mediated by gene expression in assayed genes and tissues. TGFM accounts for both co-regulation across genes and tissues and LD between SNPs (generalizing existing fine-mapping methods), and incorporates genome-wide estimates of each tissue's contribution to disease as tissue-level priors. TGFM was well-calibrated and moderately well-powered in simulations; unlike previous methods, TGFM was able to attain correct calibration by modeling uncertainty in cis-predicted expression models. We applied TGFM to 45 UK Biobank diseases/traits (average N = 316K) using eQTL data from 38 GTEx tissues. TGFM identified an average of 147 PIP > 0.5 causal genetic elements per disease/trait, of which 11% were gene-tissue pairs. Implicated gene-tissue pairs were concentrated in known disease-critical tissues, and causal genes were strongly enriched in disease-relevant gene sets. Causal gene-tissue pairs identified by TGFM recapitulated known biology (e.g., TPO-thyroid for Hypothyroidism), but also included biologically plausible novel findings (e.g., SLC20A2-artery aorta for Diastolic blood pressure). Further application of TGFM to single-cell eQTL data from 9 cell types in peripheral blood mononuclear cells (PBMC), analyzed jointly with GTEx tissues, identified 30 additional causal gene-PBMC cell type pairs at PIP > 0.5-primarily for autoimmune disease and blood cell traits, including the biologically plausible example of CD52 in classical monocyte cells for Monocyte count. In conclusion, TGFM is a robust and powerful method for fine-mapping causal tissues and genes at disease-associated loci.

Genome-wide enhancer RNA profiling adds molecular links between genetic variation and human cancers

BACKGROUND

Dysregulation of enhancer transcription occurs in multiple cancers. Enhancer RNAs (eRNAs) are transcribed products from enhancers that play critical roles in transcriptional control. Characterizing the genetic basis of eRNA expression may elucidate the molecular mechanisms underlying cancers.

METHODS

Initially, a comprehensive analysis of eRNA quantitative trait loci (eRNAQTLs) was performed in The Cancer Genome Atlas (TCGA), and functional features were characterized using multi-omics data. To establish the first eRNAQTL profiles for colorectal cancer (CRC) in China, epigenomic data were used to define active enhancers, which were subsequently integrated with transcription and genotyping data from 154 paired CRC samples. Finally, large-scale case-control studies (34,585 cases and 69,544 controls) were conducted along with multipronged experiments to investigate the potential mechanisms by which candidate eRNAQTLs affect CRC risk.

RESULTS

A total of 300,112 eRNAQTLs were identified across 30 different cancer types, which exert their influence on eRNA transcription by modulating chromatin status, binding affinity to transcription factors and RNA-binding proteins. These eRNAQTLs were found to be significantly enriched in cancer risk loci, explaining a substantial proportion of cancer heritability. Additionally, tumor-specific eRNAQTLs exhibited high responsiveness to the development of cancer. Moreover, the target genes of these eRNAs were associated with dysregulated signaling pathways and immune cell infiltration in cancer, highlighting their potential as therapeutic targets. Furthermore, multiple ethnic population studies have confirmed that an eRNAQTL rs3094296-T variant decreases the risk of CRC in populations from China (OR = 0.91, 95%CI 0.88-0.95, P = 2.92 × 10-7) and Europe (OR = 0.92, 95%CI 0.88-0.95, P = 4.61 × 10-6). Mechanistically, rs3094296 had an allele-specific effect on the transcription of the eRNA ENSR00000155786, which functioned as a transcriptional activator promoting the expression of its target gene SENP7. These two genes synergistically suppressed tumor cell proliferation. Our curated list of variants, genes, and drugs has been made available in CancereRNAQTL ( http://canernaqtl.whu.edu.cn/#/ ) to serve as an informative resource for advancing this field.

CONCLUSION

Our findings underscore the significance of eRNAQTLs in transcriptional regulation and disease heritability, pinpointing the potential of eRNA-based therapeutic strategies in cancers.