Gene Regulatory Elements, Major Drivers of Human Disease

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.

Unexplored power of CRISPR-Cas9 in neuroscience, a multi-OMICs review

The neuroscience and neurobiology of gene editing to enhance learning and memory is of paramount interest to the scientific community. The advancements of CRISPR system have created avenues to treat neurological disorders by means of versatile modalities varying from expression to suppression of genes and proteins. Neurodegenerative disorders have also been attributed to non-canonical DNA secondary structures by affecting neuron activity through controlling gene expression, nucleosome shape, transcription, translation, replication, and recombination. Changing DNA regulatory elements which could contribute to the fate and function of neurons are thoroughly discussed in this review. This study presents the ability of CRISPR system to boost learning power and memory, treat or cure genetically-based neurological disorders, and alleviate psychiatric diseases by altering the activity and the irritability of the neurons at the synaptic cleft through DNA manipulation, and also, epigenetic modifications using Cas9. We explore and examine how each different OMIC techniques can come useful when altering DNA sequences. Such insight into the underlying relationship between OMICs and cellular behaviors leads us to better neurological and psychiatric therapeutics by intelligently designing and utilizing the CRISPR/Cas9 technology.

Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A genome-wide association study identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we use mouse enhancer assays, three mouse enhancer knockouts, and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterize a sequence, PEC7, which overlaps the AIS-associated variant, and find it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, or deletion of both sequences lead to a kinky tail phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicates several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Adaptive Selection of Cis-regulatory Elements in the Han Chinese

Cis-regulatory elements have an important role in human adaptation to the living environment. However, the lag in population genomic cohort studies and epigenomic studies, hinders the research in the adaptive analysis of cis-regulatory elements in human populations. In this study, we collected 4,013 unrelated individuals and performed a comprehensive analysis of adaptive selection of genome-wide cis-regulatory elements in the Han Chinese. In total, 12.34% of genomic regions are under the influence of adaptive selection, where 1.00% of enhancers and 2.06% of promoters are under positive selection, and 0.06% of enhancers and 0.02% of promoters are under balancing selection. Gene ontology enrichment analysis of these cis-regulatory elements under adaptive selection reveals that many positive selections in the Han Chinese occur in pathways involved in cell-cell adhesion processes, and many balancing selections are related to immune processes. Two classes of adaptive cis-regulatory elements related to cell adhesion were in-depth analyzed, one is the adaptive enhancers derived from neanderthal introgression, leads to lower hyaluronidase level in skin, and brings better performance on UV-radiation resistance to the Han Chinese. Another one is the cis-regulatory elements regulating wound healing, and the results suggest the positive selection inhibits coagulation and promotes angiogenesis and wound healing in the Han Chinese. Finally, we found that many pathogenic alleles, such as risky alleles of type 2 diabetes or schizophrenia, remain in the population due to the hitchhiking effect of positive selections. Our findings will help deepen our understanding of the adaptive evolution of genome regulation in the Han Chinese.

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.

Aberrant gene expression yet undiminished retinal ganglion cell genesis in iPSC-derived models of optic nerve hypoplasia

BACKGROUND

Optic nerve hypoplasia (ONH), the leading congenital cause of permanent blindness, is characterized by a retinal ganglion cell (RGC) deficit at birth. Multifactorial developmental events are hypothesized to underlie ONH and its frequently associated neurologic and endocrine abnormalities; however, environmental influences are unclear and genetic underpinnings are unexplored. This work investigates the genetic contribution to ONH RGC production and gene expression using patient induced pluripotent stem cell (iPSC)-derived retinal organoids (ROs).

MATERIALS AND METHODS

iPSCs produced from ONH patients and controls were differentiated to ROs. RGC genesis was assessed using immunofluorescence and flow cytometry. Flow-sorted BRN3+ cells were collected for RNA extraction for RNA-Sequencing. Differential gene expression was assessed using DESeq2 and edgeR. PANTHER was employed to identify statistically over-represented ontologies among the differentially expressed genes (DEGs). DEGs of high interest to ONH were distinguished by assessing function, mutational constraint, and prior identification in ONH, autism and neurodevelopmental disorder (NDD) studies.

RESULTS

RGC genesis and survival were similar in ONH and control ROs. Differential expression of 70 genes was identified in both DESeq2 and edgeR analyses, representing a ~ 4-fold higher percentage of DEGs than in randomized study participants. DEGs showed trends towards over-representation of validated NDD genes and ONH exome variant genes. Among the DEGs, RAPGEF4 and DMD had the greatest number of disease-relevant features.

CONCLUSIONS

ONH genetic background was not associated with impaired RGC genesis but was associated with DEGs exhibiting disease contribution potential. This constitutes some of the first evidence of a genetic contribution to ONH.

CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions

Machine Learning-based scoring and classification of genetic variants aids the assessment of clinical findings and is employed to prioritize variants in diverse genetic studies and analyses. Combined Annotation-Dependent Depletion (CADD) is one of the first methods for the genome-wide prioritization of variants across different molecular functions and has been continuously developed and improved since its original publication. Here, we present our most recent release, CADD v1.7. We explored and integrated new annotation features, among them state-of-the-art protein language model scores (Meta ESM-1v), regulatory variant effect predictions (from sequence-based convolutional neural networks) and sequence conservation scores (Zoonomia). We evaluated the new version on data sets derived from ClinVar, ExAC/gnomAD and 1000 Genomes variants. For coding effects, we tested CADD on 31 Deep Mutational Scanning (DMS) data sets from ProteinGym and, for regulatory effect prediction, we used saturation mutagenesis reporter assay data of promoter and enhancer sequences. The inclusion of new features further improved the overall performance of CADD. As with previous releases, all data sets, genome-wide CADD v1.7 scores, scripts for on-site scoring and an easy-to-use webserver are readily provided via https://cadd.bihealth.org/ or https://cadd.gs.washington.edu/ to the community.

Minimum entropy framework identifies a novel class of genomic functional elements and reveals regulatory mechanisms at human disease loci

We introduce CoRE-BED, a framework trained using 19 epigenomic features in 33 major cell and tissue types to predict cell-type-specific regulatory function. CoRE-BED identifies nine functional classes de-novo, capturing both known and new regulatory categories. Notably, we describe a previously undercharacterized class that we term Development Associated Elements (DAEs), which are highly enriched in cell types with elevated regenerative potential and distinguished by the dual presence of either H3K4me2 and H3K9ac (an epigenetic signature associated with kinetochore assembly) or H3K79me3 and H4K20me1 (a signature associated with transcriptional pause release). Unlike bivalent promoters, which represent a transitory state between active and silenced promoters, DAEs transition directly to or from a non-functional state during stem cell differentiation and are proximal to highly expressed genes. CoRE-BED's interpretability facilitates causal inference and functional prioritization. Across 70 complex traits, distal insulators account for the largest mean proportion of SNP heritability (~49%) captured by the GWAS. Collectively, our results demonstrate the value of exploring non-conventional ways of regulatory classification that enrich for trait heritability, to complement existing approaches for cis-regulatory prediction.

SilenceREIN: seeking silencers on anchors of chromatin loops by deep graph neural networks

Silencers are repressive cis-regulatory elements that play crucial roles in transcriptional regulation. Experimental methods for identifying silencers are always costly and time-consuming. Computational methods, which relies on genomic sequence features, have been introduced as alternative approaches. However, silencers do not have significant epigenomic signature. Therefore, we explore a new way to computationally identify silencers, by incorporating chromatin structural information. We propose the SilenceREIN method, which focuses on finding silencers on anchors of chromatin loops. By using graph neural networks, we extracted chromatin structural information from a regulatory element interaction network. SilenceREIN integrated the chromatin structural information with linear genomic signatures to find silencers. The predictive performance of SilenceREIN is comparable or better than other states-of-the-art methods. We performed a genome-wide scanning to systematically find silencers in human genome. Results suggest that silencers are widespread on anchors of chromatin loops. In addition, enrichment analysis of transcription factor binding motif support our prediction results. As far as we can tell, this is the first attempt to incorporate chromatin structural information in finding silencers. All datasets and source codes of SilenceREIN have been deposited in a GitHub repository (https://github.com/JianHPan/SilenceREIN).

Genetic Variations in TrkB.T1 Isoform and Their Association with Somatic and Psychological Symptoms in Individuals with IBS

Irritable bowel syndrome (IBS), a disorder of gut-brain interaction, is often comorbid with somatic pain and psychological disorders. Dysregulated signaling of brain-derived neurotrophic factor (BDNF) and its receptor, tropomyosin-related kinase B (TrkB), has been implicated in somatic-psychological symptoms in individuals with IBS. Thus, we investigated the association of 10 single nucleotide polymorphisms (SNPs) in the regulatory 3' untranslated region (UTR) of NTRK2 (TrkB) kinase domain-deficient truncated isoform (TrkB.T1) and the BDNF Val66Met SNP with somatic and psychological symptoms and quality of life in a U.S. cohort (IBS n=464; healthy controls n=156). We found that the homozygous recessive genotype (G/G) of rs2013566 in individuals with IBS is associated with worsened somatic symptoms, including headache, back pain, joint pain, muscle pain, and somatization as well as diminished sleep quality, energy level and overall quality of life. Validation using U.K. BioBank (UKBB) data confirmed the association of rs2013566 with increased likelihood of headache. Several SNPs (rs1627784, rs1624327, rs1147198) showed significant associations with muscle pain in our U.S. cohort. Notably, these SNPs are predominantly located in H3K4Me1-enriched regions, suggesting their enhancer and/or transcription regulation potential. Together, our findings suggest that genetic variation within the 3'UTR region of the TrkB.T1 isoform may contribute to comorbid conditions in individuals with IBS, resulting in a spectrum of somatic and psychological symptoms that may influence their quality of life. These findings advance our understanding of the genetic interaction between BDNF/TrkB pathways and somatic-psychological symptoms in IBS, highlighting the importance of further exploring this interaction for potential clinical applications.

Human genomic data have different statistical properties than the data of randomised controlled trials

Madole & Harden argue that the Mendelian reshuffling of genes and genomes is analogous to randomised controlled trials. We are not convinced by their arguments. First, their recipe for meeting the demands on randomised experiments is inherently inconsistent. Second, disequilibrium across chromosomes conflicts with their assumption of statistical independence. Third, the genome-wide association study (GWAS) method has many pitfalls, including low repeatability.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

SCING: Inference of robust, interpretable gene regulatory networks from single cell and spatial transcriptomics

Gene regulatory network (GRN) inference is an integral part of understanding physiology and disease. Single cell/nuclei RNA-seq (scRNA-seq/snRNA-seq) data has been used to elucidate cell-type GRNs; however, the accuracy and speed of current scRNAseq-based GRN approaches are suboptimal. Here, we present Single Cell INtegrative Gene regulatory network inference (SCING), a gradient boosting and mutual information-based approach for identifying robust GRNs from scRNA-seq, snRNA-seq, and spatial transcriptomics data. Performance evaluation using Perturb-seq datasets, held-out data, and the mouse cell atlas combined with the DisGeNET database demonstrates the improved accuracy and biological interpretability of SCING compared to existing methods. We applied SCING to the entire mouse single cell atlas, human Alzheimer's disease (AD), and mouse AD spatial transcriptomics. SCING GRNs reveal unique disease subnetwork modeling capabilities, have intrinsic capacity to correct for batch effects, retrieve disease relevant genes and pathways, and are informative on spatial specificity of disease pathogenesis.

Integrative epigenomic and functional characterization assay based annotation of regulatory activity across diverse human cell types

We introduce ChromActivity, a computational framework for predicting and annotating regulatory activity across the genome through integration of multiple epigenomic maps and various functional characterization datasets. ChromActivity generates genomewide predictions of regulatory activity associated with each functional characterization dataset across many cell types based on available epigenomic data. It then for each cell type produces (1) ChromScoreHMM genome annotations based on the combinatorial and spatial patterns within these predictions and (2) ChromScore tracks of overall predicted regulatory activity. ChromActivity provides a resource for analyzing and interpreting the human regulatory genome across diverse cell types.

Forkhead transcription factor FKH-8 cooperates with RFX in the direct regulation of sensory cilia in Caenorhabditis elegans

Cilia, either motile or non-motile (a.k.a primary or sensory), are complex evolutionarily conserved eukaryotic structures composed of hundreds of proteins required for their assembly, structure and function that are collectively known as the ciliome. Ciliome gene mutations underlie a group of pleiotropic genetic diseases known as ciliopathies. Proper cilium function requires the tight coregulation of ciliome gene transcription, which is only fragmentarily understood. RFX transcription factors (TF) have an evolutionarily conserved role in the direct activation of ciliome genes both in motile and non-motile cilia cell-types. In vertebrates, FoxJ1 and FoxN4 Forkhead (FKH) TFs work with RFX in the direct activation of ciliome genes, exclusively in motile cilia cell-types. No additional TFs have been described to act together with RFX in primary cilia cell-types in any organism. Here we describe FKH-8, a FKH TF, as a direct regulator of the sensory ciliome genes in Caenorhabditis elegans. FKH-8 is expressed in all ciliated neurons in C. elegans, binds the regulatory regions of ciliome genes, regulates ciliome gene expression, cilium morphology and a wide range of behaviors mediated by sensory ciliated neurons. FKH-8 and DAF-19 (C. elegans RFX) physically interact and synergistically regulate ciliome gene expression. C. elegans FKH-8 function can be replaced by mouse FOXJ1 and FOXN4 but not by other members of other mouse FKH subfamilies. In conclusion, RFX and FKH TF families act jointly as direct regulators of ciliome genes also in sensory ciliated cell types suggesting that this regulatory logic could be an ancient trait predating functional cilia sub-specialization.

Defining the fine structure of promoter activity on a genome-wide scale with CISSECTOR

Classic promoter mutagenesis strategies can be used to study how proximal promoter regions regulate the expression of particular genes of interest. This is a laborious process, in which the smallest sub-region of the promoter still capable of recapitulating expression in an ectopic setting is first identified, followed by targeted mutation of putative transcription factor binding sites. Massively parallel reporter assays such as survey of regulatory elements (SuRE) provide an alternative way to study millions of promoter fragments in parallel. Here we show how a generalized linear model (GLM) can be used to transform genome-scale SuRE data into a high-resolution genomic track that quantifies the contribution of local sequence to promoter activity. This coefficient track helps identify regulatory elements and can be used to predict promoter activity of any sub-region in the genome. It thus allows in silico dissection of any promoter in the human genome to be performed. We developed a web application, available at cissector.nki.nl, that lets researchers easily perform this analysis as a starting point for their research into any promoter of interest.

Calling Cards: a customizable platform to longitudinally record protein-DNA interactions over time in cells and tissues

Calling Cards is a platform technology to record a cumulative history of transient protein-DNA interactions in the genome of genetically targeted cell types. The record of these interactions is recovered by next generation sequencing. Compared to other genomic assays, whose readout provides a snapshot at the time of harvest, Calling Cards enables correlation of historical molecular states to eventual outcomes or phenotypes. To achieve this, Calling Cards uses the piggyBac transposase to insert self-reporting transposon (SRT) "Calling Cards" into the genome, leaving permanent marks at interaction sites. Calling Cards can be deployed in a variety of in vitro and in vivo biological systems to study gene regulatory networks involved in development, aging, and disease. Out of the box, it assesses enhancer usage but can be adapted to profile specific transcription factor binding with custom transcription factor (TF)-piggyBac fusion proteins. The Calling Cards workflow has five main stages: delivery of Calling Card reagents, sample preparation, library preparation, sequencing, and data analysis. Here, we first present a comprehensive guide for experimental design, reagent selection, and optional customization of the platform to study additional TFs. Then, we provide an updated protocol for the five steps, using reagents that improve throughput and decrease costs, including an overview of a newly deployed computational pipeline. This protocol is designed for users with basic molecular biology experience to process samples into sequencing libraries in 1-2 days. Familiarity with bioinformatic analysis and command line tools is required to set up the pipeline in a high-performance computing environment and to conduct downstream analyses. Basic Protocol 1: Preparation and delivery of Calling Cards reagentsBasic Protocol 2: Sample preparationBasic Protocol 3: Sequencing library preparationBasic Protocol 4: Library pooling and sequencingBasic Protocol 5: Data analysis.

Short tandem repeats are important contributors to silencer elements in T cells

The action of cis-regulatory elements with either activation or repression functions underpins the precise regulation of gene expression during normal development and cell differentiation. Gene activation by the combined activities of promoters and distal enhancers has been extensively studied in normal and pathological contexts. In sharp contrast, gene repression by cis-acting silencers, defined as genetic elements that negatively regulate gene transcription in a position-independent fashion, is less well understood. Here, we repurpose the STARR-seq approach as a novel high-throughput reporter strategy to quantitatively assess silencer activity in mammals. We assessed silencer activity from DNase hypersensitive I sites in a mouse T cell line. Identified silencers were associated with either repressive or active chromatin marks and enriched for binding motifs of known transcriptional repressors. CRISPR-mediated genomic deletions validated the repressive function of distinct silencers involved in the repression of non-T cell genes and genes regulated during T cell differentiation. Finally, we unravel an association of silencer activity with short tandem repeats, highlighting the role of repetitive elements in silencer activity. Our results provide a general strategy for genome-wide identification and characterization of silencer elements.

A Variant in the IRF6 Promoter Associated with the Risk for Orofacial Clefting

The single-nucleotide polymorphism (SNP) rs2235371 (IRF6 V274I) is associated with nonsyndromic cleft lip with or without cleft palate (NSCL/P) in Han Chinese and other populations but appears to be without a functional effect. To find the common etiologic variant or variants within the haplotype tagged by rs2235371, we carried out targeted sequencing of an interval containing IRF6 in 159 Han Chinese with NSCL/P. This study revealed that the SNP rs12403599, within the IRF6 promoter, is associated with all phenotypes of NSCL/P, especially nonsyndromic cleft lip (NSCLO) and a subphenotype of it, microform cleft lip (MCL). This association was replicated in 2 additional much larger cohorts of cases and controls from the Han Chinese. Conditional logistic analysis indicated that association of rs2235371 with NSCL/P was lost if rs12403599 was excluded. rs12403599 contributes the most risk to MCL: its G allele is responsible for 38.47% of the genetic contribution to MCL, and the odds ratios of G/C and G/G genotypes were 2.91 and 6.58, respectively, for MCL. To test if rs12403599 is functional, we carried out reporter assays in a fetal oral epithelium cells (GMSM-K). Unexpectedly, the risk allele G yielded higher promoter activity in GMSM-K. Consistent with the reporter studies, expression of IRF6 in lip tissues from NSCLO and MCL patients with the G/G phenotype was higher than in those from patients with the C/C phenotype. These results indicate that rs12403599 is tagging the risk haplotype for NSCL/P better than rs2235371 in Han Chinese and supports investigation of the mechanisms by which the allele of rs12403599 affects IRF6 expression and tests of this association in different populations.

Transcription factor binding site orientation and order are major drivers of gene regulatory activity

The gene regulatory code and grammar remain largely unknown, precluding our ability to link phenotype to genotype in regulatory sequences. Here, using a massively parallel reporter assay (MPRA) of 209,440 sequences, we examine all possible pair and triplet combinations, permutations and orientations of eighteen liver-associated transcription factor binding sites (TFBS). We find that TFBS orientation and order have a major effect on gene regulatory activity. Corroborating these results with genomic analyses, we find clear human promoter TFBS orientation biases and similar TFBS orientation and order transcriptional effects in an MPRA that tested 164,307 liver candidate regulatory elements. Additionally, by adding TFBS orientation to a model that predicts expression from sequence we improve performance by 7.7%. Collectively, our results show that TFBS orientation and order have a significant effect on gene regulatory activity and need to be considered when analyzing the functional effect of variants on the activity of these sequences.

Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A GWAS identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we used mouse enhancer assays, three mouse enhancer knockouts and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterized a sequence, PEC7, that overlaps the AIS-associated variant, and found it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, and deletion of both sequences led to a kinky phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicated several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Massively parallel characterization of transcriptional regulatory elements in three diverse human cell types

The human genome contains millions of candidate cis-regulatory elements (CREs) with cell-type-specific activities that shape both health and myriad disease states. However, we lack a functional understanding of the sequence features that control the activity and cell-type-specific features of these CREs. Here, we used lentivirus-based massively parallel reporter assays (lentiMPRAs) to test the regulatory activity of over 680,000 sequences, representing a nearly comprehensive set of all annotated CREs among three cell types (HepG2, K562, and WTC11), finding 41.7% to be functional. By testing sequences in both orientations, we find promoters to have significant strand orientation effects. We also observe that their 200 nucleotide cores function as non-cell-type-specific 'on switches' providing similar expression levels to their associated gene. In contrast, enhancers have weaker orientation effects, but increased tissue-specific characteristics. Utilizing our lentiMPRA data, we develop sequence-based models to predict CRE function with high accuracy and delineate regulatory motifs. Testing an additional lentiMPRA library encompassing 60,000 CREs in all three cell types, we further identified factors that determine cell-type specificity. Collectively, our work provides an exhaustive catalog of functional CREs in three widely used cell lines, and showcases how large-scale functional measurements can be used to dissect regulatory grammar.

Massively parallel characterization of psychiatric disorder-associated and cell-type-specific regulatory elements in the developing human cortex

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and disease. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 sequences, including differentially accessible cell-type specific regions in the developing cortex and single-nucleotide variants associated with psychiatric disorders. In primary cells, we identified 46,802 active enhancer sequences and 164 disorder-associated variants that significantly alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning, we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

Hnf1b renal expression directed by a distal enhancer responsive to Pax8

Xenopus provides a simple and efficient model system to study nephrogenesis and explore the mechanisms causing renal developmental defects in human. Hnf1b (hepatocyte nuclear factor 1 homeobox b), a gene whose mutations are the most commonly identified genetic cause of developmental kidney disease, is required for the acquisition of a proximo-intermediate nephron segment in Xenopus as well as in mouse. Genetic networks involved in Hnf1b expression during kidney development remain poorly understood. We decided to explore the transcriptional regulation of Hnf1b in the developing Xenopus pronephros and mammalian renal cells. Using phylogenetic footprinting, we identified an evolutionary conserved sequence (CNS1) located several kilobases (kb) upstream the Hnf1b transcription start and harboring epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. By means of functional expression assays in Xenopus and mammalian renal cell lines we showed that CNS1 displays enhancer activity in renal tissue. Using CRISPR/cas9 editing in Xenopus tropicalis, we demonstrated the in vivo functional relevance of CNS1 in driving hnf1b expression in the pronephros. We further showed the importance of Pax8-CNS1 interaction for CNS1 enhancer activity allowing us to conclude that Hnf1b is a direct target of Pax8. Our work identified for the first time a Hnf1b renal specific enhancer and may open important perspectives into the diagnosis for congenital kidney anomalies in human, as well as modeling HNF1B-related diseases.

HiChIPdb: a comprehensive database of HiChIP regulatory interactions

Elucidating the role of 3D architecture of DNA in gene regulation is crucial for understanding cell differentiation, tissue homeostasis and disease development. Among various chromatin conformation capture methods, HiChIP has received increasing attention for its significant improvement over other methods in profiling of regulatory (e.g. H3K27ac) and structural (e.g. cohesin) interactions. To facilitate the studies of 3D regulatory interactions, we developed a HiChIP interactions database, HiChIPdb (http://health.tsinghua.edu.cn/hichipdb/). The current version of HiChIPdb contains ∼262M annotated HiChIP interactions from 200 high-throughput HiChIP samples across 108 cell types. The functionalities of HiChIPdb include: (i) standardized categorization of HiChIP interactions in a hierarchical structure based on organ, tissue and cell line and (ii) comprehensive annotations of HiChIP interactions with regulatory genes and GWAS Catalog SNPs. To the best of our knowledge, HiChIPdb is the first comprehensive database that utilizes a unified pipeline to map the functional interactions across diverse cell types and tissues in different resolutions. We believe this database has the potential to advance cutting-edge research in regulatory mechanisms in development and disease by removing the barrier in data aggregation, preprocessing, and analysis.

RegEl corpus: identifying DNA regulatory elements in the scientific literature

High-throughput technologies led to the generation of a wealth of data on regulatory DNA elements in the human genome. However, results from disease-driven studies are primarily shared in textual form as scientific articles. Information extraction (IE) algorithms allow this information to be (semi-)automatically accessed. Their development, however, is dependent on the availability of annotated corpora. Therefore, we introduce RegEl (Regulatory Elements), the first freely available corpus annotated with regulatory DNA elements comprising 305 PubMed abstracts for a total of 2690 sentences. We focus on enhancers, promoters and transcription factor binding sites. Three annotators worked in two stages, achieving an overall 0.73 F1 inter-annotator agreement and 0.46 for regulatory elements. Depending on the entity type, IE baselines reach F1-scores of 0.48–0.91 for entity detection and 0.71–0.88 for entity normalization. Next, we apply our entity detection models to the entire PubMed collection and extract co-occurrences of genes or diseases with regulatory elements. This generates large collections of regulatory elements associated with 137 870 unique genes and 7420 diseases, which we make openly available.

Database URL: https://zenodo.org/record/6418451#.YqcLHvexVqg

Rewiring of 3D Chromatin Topology Orchestrates Transcriptional Reprogramming and the Development of Human Dilated Cardiomyopathy

BACKGROUND

Transcriptional reconfiguration is central to heart failure, the most common cause of which is dilated cardiomyopathy (DCM). The effect of 3-dimensional chromatin topology on transcriptional dysregulation and pathogenesis in human DCM remains elusive.

METHODS

We generated a compendium of 3-dimensional epigenome and transcriptome maps from 101 biobanked human DCM and nonfailing heart tissues through highly integrative chromatin immunoprecipitation (H3K27ac [acetylation of lysine 27 on histone H3]), in situ high-throughput chromosome conformation capture, chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin using sequencing, and RNA sequencing. We used human induced pluripotent stem cell-derived cardiomyocytes and mouse models to interrogate the key transcription factor implicated in 3-dimensional chromatin organization and transcriptional regulation in DCM pathogenesis.

RESULTS

We discovered that the active regulatory elements (H3K27ac peaks) and their connectome (H3K27ac loops) were extensively reprogrammed in DCM hearts and contributed to transcriptional dysregulation implicated in DCM development. For example, we identified that nontranscribing NPPA-AS1 (natriuretic peptide A antisense RNA 1) promoter functions as an enhancer and physically interacts with the NPPA (natriuretic peptide A) and NPPB (natriuretic peptide B) promoters, leading to the cotranscription of NPPA and NPPB in DCM hearts. We revealed that DCM-enriched H3K27ac loops largely resided in conserved high-order chromatin architectures (compartments, topologically associating domains) and their anchors unexpectedly had equivalent chromatin accessibility. We discovered that the DCM-enriched H3K27ac loop anchors exhibited a strong enrichment for HAND1 (heart and neural crest derivatives expressed 1), a key transcription factor involved in early cardiogenesis. In line with this, its protein expression was upregulated in human DCM and mouse failing hearts. To further validate whether HAND1 is a causal driver for the reprogramming of enhancer-promoter connectome in DCM hearts, we performed comprehensive 3-dimensional epigenome mappings in human induced pluripotent stem cell-derived cardiomyocytes. We found that forced overexpression of HAND1 in human induced pluripotent stem cell-derived cardiomyocytes induced a distinct gain of enhancer-promoter connectivity and correspondingly increased the expression of their connected genes implicated in DCM pathogenesis, thus recapitulating the transcriptional signature in human DCM hearts. Electrophysiology analysis demonstrated that forced overexpression of HAND1 in human induced pluripotent stem cell-derived cardiomyocytes induced abnormal calcium handling. Furthermore, cardiomyocyte-specific overexpression of Hand1 in the mouse hearts resulted in dilated cardiac remodeling with impaired contractility/Ca2+ handling in cardiomyocytes, increased ratio of heart weight/body weight, and compromised cardiac function, which were ascribed to recapitulation of transcriptional reprogramming in DCM.

CONCLUSIONS

This study provided novel chromatin topology insights into DCM pathogenesis and illustrated a model whereby a single transcription factor (HAND1) reprograms the genome-wide enhancer-promoter connectome to drive DCM pathogenesis.

Clinical impact of variants in non-coding regions of SHOX - Current knowledge

The short stature homeobox-containing (SHOX) is the most frequently analysed gene in patients classified as short stature patients (ISS) or diagnosed with Leri-Weill dyschondrosteosis (LWD), Langer mesomelic dysplasia (LMD), or Madelung deformity (MD). However, clinical testing of this gene focuses primarily on single nucleotide variants (SNV) in its coding sequences and copy number variants (CNV) overlapping SHOX gene. This review summarizes the clinical impact of variants in noncoding regions of SHOX. RECENT FINDINGS: CNV extending exclusively into the regulatory elements (i.e., not interrupting the coding sequence) are found more frequently in downstream regulatory elements of SHOX. Further, duplications are more frequent than deletions. Interestingly, downstream duplications are more common than deletions in patients with ISS or LWD but no such differences exist for upstream CNV. Moreover, the presence of specific CNVs in the patient population suggests the involvement of additional unknown factors. Some of its intronic variants, notably NM_000451.3(SHOX):c.-9delG and c.-65C>A in the 5'UTR, have unclear clinical roles. However, these intronic SNV may increase the probability that other CNV will arise de novo in the SHOX gene based on homologous recombination or incorrect splicing of mRNA. SUMMARY: This review highlights the clinical impact of noncoding changes in the SHOX gene and the need to apply new technologies and genotype-phenotype correlation in their analysis.

Best practices for the interpretation and reporting of clinical whole genome sequencing

Whole genome sequencing (WGS) shows promise as a first-tier diagnostic test for patients with rare genetic disorders. However, standards addressing the definition and deployment practice of a best-in-class test are lacking. To address these gaps, the Medical Genome Initiative, a consortium of leading health care and research organizations in the US and Canada, was formed to expand access to high quality clinical WGS by convening experts and publishing best practices. Here, we present best practice recommendations for the interpretation and reporting of clinical diagnostic WGS, including discussion of challenges and emerging approaches that will be critical to harness the full potential of this comprehensive test.

Computational and experimental methods for classifying variants of unknown clinical significance

The increase in sequencing capacity, reduction in costs, and national and international coordinated efforts have led to the widespread introduction of next-generation sequencing (NGS) technologies in patient care. More generally, human genetics and genomic medicine are gaining importance for more and more patients. Some communities are already discussing the prospect of sequencing each individual's genome at time of birth. Together with digital health records, this shall enable individualized treatments and preventive measures, so-called precision medicine. A central step in this process is the identification of disease causal mutations or variant combinations that make us more susceptible for diseases. Although various technological advances have improved the identification of genetic alterations, the interpretation and ranking of the identified variants remains a major challenge. Based on our knowledge of molecular processes or previously identified disease variants, we can identify potentially functional genetic variants and, using different lines of evidence, we are sometimes able to demonstrate their pathogenicity directly. However, the vast majority of variants are classified as variants of uncertain clinical significance (VUSs) with not enough experimental evidence to determine their pathogenicity. In these cases, computational methods may be used to improve the prioritization and an increasing toolbox of experimental methods is emerging that can be used to assay the molecular effects of VUSs. Here, we discuss how computational and experimental methods can be used to create catalogs of variant effects for a variety of molecular and cellular phenotypes. We discuss the prospects of integrating large-scale functional data with machine learning and clinical knowledge for the development of accurate pathogenicity predictions for clinical applications.

Evaluation of chromatin mesoscale organization

Chromatin organization in the nucleus represents an important aspect of transcription regulation. Most of the studies so far focused on the chromatin structure in cultured cells or in fixed tissue preparations. Here, we discuss the various approaches for deciphering chromatin 3D organization with an emphasis on the advantages of live imaging approaches.

Developmental disorders caused by haploinsufficiency of transcriptional regulators: a perspective based on cell fate determination

Many human birth defects and neurodevelopmental disorders are caused by loss-of-function mutations in a single copy of transcription factor (TF) and chromatin regulator genes. Although this dosage sensitivity has long been known, how and why haploinsufficiency (HI) of transcriptional regulators leads to developmental disorders (DDs) is unclear. Here I propose the hypothesis that such DDs result from defects in cell fate determination that are based on disrupted bistability in the underlying gene regulatory network (GRN). Bistability, a crucial systems biology concept to model binary choices such as cell fate decisions, requires both positive feedback and ultrasensitivity, the latter often achieved through TF cooperativity. The hypothesis explains why dosage sensitivity of transcriptional regulators is an inherent property of fate decisions, and why disruption of either positive feedback or cooperativity in the underlying GRN is sufficient to cause disease. I present empirical and theoretical evidence in support of this hypothesis and discuss several issues for which it increases our understanding of disease, such as incomplete penetrance. The proposed framework provides a mechanistic, systems-level explanation of HI of transcriptional regulators, thus unifying existing theories, and offers new insights into outstanding issues of human disease. This article has an associated Future Leader to Watch interview with the author of the paper.

Fine mapping with epigenetic information and 3D structure

Since 2005, thousands of genome-wide association studies (GWAS) have been published, identifying hundreds of thousands of genetic variants that increase risk of complex traits such as autoimmune diseases. This wealth of data has the potential to improve patient care, through personalized medicine and the identification of novel drug targets. However, the potential of GWAS for clinical translation has not been fully achieved yet, due to the fact that the functional interpretation of risk variants and the identification of causal variants and genes are challenging. The past decade has seen the development of great advances that are facilitating the overcoming of these limitations, by utilizing a plethora of genomics and epigenomics tools to map and characterize regulatory elements and chromatin interactions, which can be used to fine map GWAS loci, and advance our understanding of the biological mechanisms that cause disease.

Controllable SiOx Nanorod Memristive Neuron for Probabilistic Bayesian Inference

Modern artificial neural network technology using a deterministic computing framework is faced with a critical challenge in dealing with massive data that are largely unstructured and ambiguous. This challenge demands the advances of an elementary physical device for tackling these uncertainties. Here, we designed and fabricated a SiO_x nanorod memristive device by employing the glancing angle deposition (GLAD) technique, suggesting a controllable stochastic artificial neuron that can mimic the fundamental integrate-and-fire signaling and stochastic dynamics of a biological neuron. The nanorod structure provides the random distribution of multiple nanopores all across the active area, capable of forming a multitude of Si filaments at many SiO_x nanorod edges after the electromigration process, leading to a stochastic switching event with very high dynamic range (≈5.15 × 10¹⁰ ) and low energy (≈4.06 pJ). Different probabilistic activation (ProbAct) functions in a sigmoid form are implemented, showing its controllability with low variation by manufacturing and electrical programming schemes. Furthermore, as an application prospect, based on the suggested memristive neuron, we demonstrated the self-resting neural operation with the local circuit configuration and revealed probabilistic Bayesian inferences for genetic regulatory networks with low normalized mean squared errors (≈2.41 × 10^-2 ) and its robustness to the ProbAct variation.

Methods of massive parallel reporter assays for investigation of enhancers

The correct deployment of genetic programs for development and differentiation relies on finely coordinated regulation of specific gene sets. Genomic regulatory elements play an exceptional role in this process. There are few types of gene regulatory elements, including promoters, enhancers, insulators and silencers. Alterations of gene regulatory elements may cause various pathologies, including cancer, congenital disorders and autoimmune diseases. The development of high-throughput genomic assays has made it possible to significantly accelerate the accumulation of information about the characteristic epigenetic properties of regulatory elements. In combination with high-throughput studies focused on the genome-wide distribution of epigenetic marks, regulatory proteins and the spatial structure of chromatin, this significantly expands the understanding of the principles of epigenetic regulation of genes and allows potential regulatory elements to be searched for in silico. However, common experimental approaches used to study the local characteristics of chromatin have a number of technical limitations that may reduce the reliability of computational identification of genomic regulatory sequences. Taking into account the variability of the functions of epigenetic determinants and complex multicomponent regulation of genomic elements activity, their functional verification is often required. A plethora of methods have been developed to study the functional role of regulatory elements on the genome scale. Common experimental approaches for in silico identification of regulatory elements and their inherent technical limitations will be described. The present review is focused on original high-throughput methods of enhancer activity reporter analysis that are currently used to validate predicted regulatory elements and to perform de novo searches. The methods described allow assessing the functional role of the nucleotide sequence of a regulatory element, to determine its exact boundaries and to assess the influence of the local state of chromatin on the activity of enhancers and gene expression. These approaches have contributed substantially to the understanding of the fundamental principles of gene regulation.

RefSeq Functional Elements as experimentally assayed nongenic reference standards and functional interactions in human and mouse

Eukaryotic genomes contain many nongenic elements that function in gene regulation, chromosome organization, recombination, repair, or replication, and mutation of those elements can affect genome function and cause disease. Although numerous epigenomic studies provide high coverage of gene regulatory regions, those data are not usually exposed in traditional genome annotation and can be difficult to access and interpret without field-specific expertise. The National Center for Biotechnology Information (NCBI) therefore provides RefSeq Functional Elements (RefSeqFEs), which represent experimentally validated human and mouse nongenic elements derived from the literature. The curated data set is comprised of richly annotated sequence records, descriptive records in the NCBI Gene database, reference genome feature annotation, and activity-based interactions between nongenic regions, target genes, and each other. The data set provides succinct functional details and transparent experimental evidence, leverages data from multiple experimental sources, is readily accessible and adaptable, and uses a flexible data model. The data have multiple uses for basic functional discovery, bioinformatics studies, genetic variant interpretation; as known positive controls for epigenomic data evaluation; and as reference standards for functional interactions. Comparisons to other gene regulatory data sets show that the RefSeqFE data set includes a wider range of feature types representing more areas of biology, but it is comparatively smaller and subject to data selection biases. RefSeqFEs thus provide an alternative and complementary resource for experimentally assayed functional elements, with future data set growth expected.

Polymorphic Variants in the GRK5 Gene Promoter Are Associated With Diastolic Dysfunction in Coronary Artery Bypass Graft Surgery Patients

BACKGROUND

The G-protein-coupled receptor kinase 5 (GRK5) is a mediator of cardiovascular homeostasis and participates in inflammation and cardiac fibrosis, both being involved in the development of diastolic dysfunction (DD). While mechanisms of transcriptional regulation of the GRK5 promoter are unclear, we tested the hypotheses, that (1) GRK5 expression varies depending on functional single nucleotide polymorphisms (SNPs) in the GRK5 promoter and (2) this is associated with DD in patients undergoing coronary artery bypass graft (CABG) surgery.

METHODS

We amplified and sequenced the GRK5 promoter followed by cloning, reporter assays, and electrophoretic mobility shift assays (EMSA). GRK5 messenger ribonucleic acid (mRNA) expression was determined in right atrial tissue sampled from 50 patients undergoing CABG surgery. In another prospective study, GRK5 genotypes were associated with determinants of diastolic function using transesophageal echocardiography in 255 patients with CABG with normal systolic left ventricular (LV) function. Specifically, we measured ejection fraction (EF), transmitral Doppler early filling velocity (E), tissue Doppler early diastolic lateral mitral annular velocity (E' lateral), and calculated E/E', E' norm and the difference of E' lateral and E' norm to account for age-related changes in diastolic function.

RESULTS

We identified 6 SNPs creating 3 novel haplotypes with the greatest promoter activation in haplotype tagging (ht) SNP T(-678)C T-allele constructs (P < .001). EMSAs showed allele-specific transcription factor binding proving functional activity. GRK5 mRNA expression was greatest in TT genotypes (TT: 131 fg/µg [95% CI, 108-154]; CT: 109 [95% confidence interval {CI}, 93-124]; CC: 83 [95% CI, 54-112]; P = .012). Moreover, GRK5 genotypes were significantly associated with determinants of diastolic function. Grading of DD revealed more grade 3 patients in TT compared to CT and CC genotypes (58% vs 38% vs 4%; P = .023). E´ lateral was lowest in TT genotypes (P = .007) and corresponding E/E' measurements showed 1.27-fold increased values in TT versus CC genotypes (P = .01), respectively. While E' norm values were not different between genotypes (P = .182), the difference between E' lateral and E' norm was significantly higher in TT genotypes compared to CC and CT genotypes (-1.2 [interquartile range {IQR}, 2.7], -0.5 [IQR, 3.4], and -0.4 [IQR, 4.2; P = .035], respectively).

CONCLUSIONS

A functional GRK5 SNP results in allele-dependent differences in GRK5 promoter activity and mRNA expression. This is associated with altered echocardiographic determinants of diastolic function. Thus, SNPs in the GRK5 promoter are associated with altered perioperative diastolic cardiac function. In the future, preoperative testing for these and other SNPs might allow to initiate more specific diagnostic and perioperative pathways to benefit patients at risk.

Mechanistic insights into genetic susceptibility to prostate cancer

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

A multi-enhancer RET regulatory code is disrupted in Hirschsprung disease

The major genetic risk factors for Hirschsprung disease (HSCR) are three common polymorphisms within cis-regulatory elements (CREs) of the receptor tyrosine kinase gene RET, which reduce its expression during enteric nervous system (ENS) development. These risk variants attenuate binding of the transcription factors RARB, GATA2, and SOX10 to their cognate CREs, reduce RET gene expression, and dysregulate other ENS and HSCR genes in the RET-EDNRB gene regulatory network (GRN). Here, we use siRNA, ChIP, and CRISPR-Cas9 deletion analyses in the SK-N-SH cell line to ask how many additional HSCR-associated risk variants reside in RET CREs that affect its gene expression. We identify 22 HSCR-associated variants in candidate RET CREs, of which seven have differential allele-specific in vitro enhancer activity, and four of these seven affect RET gene expression; of these, two enhancers are bound by the transcription factor PAX3. We also show that deleting multiple variant-containing enhancers leads to synergistic effects on RET gene expression. These, coupled with our prior results, show that common sequence variants in at least 10 RET enhancers affect HSCR risk, seven with experimental evidence of affecting RET gene expression, extending the known RET-EDNRB GRN to reveal an extensive regulatory code modulating disease risk at a single gene.

Quantitative spatial and temporal assessment of regulatory element activity in zebrafish

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

Epromoters function as a hub to recruit key transcription factors required for the inflammatory response

Gene expression is controlled by the involvement of gene-proximal (promoters) and distal (enhancers) regulatory elements. Our previous results demonstrated that a subset of gene promoters, termed Epromoters, work as bona fide enhancers and regulate distal gene expression. Here, we hypothesized that Epromoters play a key role in the coordination of rapid gene induction during the inflammatory response. Using a high-throughput reporter assay we explored the function of Epromoters in response to type I interferon. We find that clusters of IFNa-induced genes are frequently associated with Epromoters and that these regulatory elements preferentially recruit the STAT1/2 and IRF transcription factors and distally regulate the activation of interferon-response genes. Consistently, we identified and validated the involvement of Epromoter-containing clusters in the regulation of LPS-stimulated macrophages. Our findings suggest that Epromoters function as a local hub recruiting the key TFs required for coordinated regulation of gene clusters during the inflammatory response.

Building a Mammalian Retina: An Eye on Chromatin Structure

Regulation of gene expression by chromatin structure has been under intensive investigation, establishing nuclear organization and genome architecture as a potent and effective means of regulating developmental processes. The substantial growth in our knowledge of the molecular mechanisms underlying retinogenesis has been powered by several genome-wide based tools that mapped chromatin organization at multiple cellular and biochemical levels. Studies profiling the retinal epigenome and transcriptome have allowed the systematic annotation of putative cis-regulatory elements associated with transcriptional programs that drive retinal neural differentiation, laying the groundwork to understand spatiotemporal retinal gene regulation at a mechanistic level. In this review, we outline recent advances in our understanding of the chromatin architecture in the mammalian retina during development and disease. We focus on the emerging roles of non-coding regulatory elements in controlling retinal cell-type specific transcriptional programs, and discuss potential implications in untangling the etiology of eye-related disorders.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance

BACKGROUND

Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging.

METHODS

Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain.

RESULTS

Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism.

CONCLUSION

This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

Targeted chromatin conformation analysis identifies novel distal neural enhancers of ZEB2 in pluripotent stem cell differentiation

The transcription factor zinc finger E-box binding protein 2 (ZEB2) controls embryonic and adult cell fate decisions and cellular maturation in many stem/progenitor cell types. Defects in these processes in specific cell types underlie several aspects of Mowat–Wilson syndrome (MOWS), which is caused by ZEB2 haplo-insufficiency. Human ZEB2, like mouse Zeb2, is located on chromosome 2 downstream of a ±3.5 Mb-long gene-desert, lacking any protein-coding gene. Using temporal targeted chromatin capture (T2C), we show major chromatin structural changes based on mapping in-cis proximities between the ZEB2 promoter and this gene desert during neural differentiation of human-induced pluripotent stem cells, including at early neuroprogenitor cell (NPC)/rosette state, where ZEB2 mRNA levels increase significantly. Combining T2C with histone-3 acetylation mapping, we identified three novel candidate enhancers about 500 kb upstream of the ZEB2 transcription start site. Functional luciferase-based assays in heterologous cells and NPCs reveal co-operation between these three enhancers. This study is the first to document in-cis Regulatory Elements located in ZEB2’s gene desert. The results further show the usability of T2C for future studies of ZEB2 REs in differentiation and maturation of multiple cell types and the molecular characterization of newly identified MOWS patients that lack mutations in ZEB2 protein-coding exons.

Study on potential differentially expressed genes in stroke by bioinformatics analysis

Stroke is a sudden cerebrovascular circulatory disorder with high morbidity, disability, mortality, and recurrence rate, but its pathogenesis and key genes are still unclear. In this study, bioinformatics was used to deeply analyze the pathogenesis of stroke and related key genes, so as to study the potential pathogenesis of stroke and provide guidance for clinical treatment. Gene Expression profiles of GSE58294 and GSE16561 were obtained from Gene Expression Omnibus (GEO), the differentially expressed genes (DEGs) were identified between IS and normal control group. The different expression genes (DEGs) between IS and normal control group were screened with the GEO2R online tool. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were performed. Using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and gene set enrichment analysis (GSEA), the function and pathway enrichment analysis of DEGS were performed. Then, a protein–protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. Cytoscape with CytoHubba were used to identify the hub genes. Finally, NetworkAnalyst was used to construct the targeted microRNAs (miRNAs) of the hub genes. A total of 85 DEGs were screened out in this study, including 65 upward genes and 20 downward genes. In addition, 3 KEGG pathways, cytokine − cytokine receptor interaction, hematopoietic cell lineage, B cell receptor signaling pathway, were significantly enriched using a database for labeling, visualization, and synthetic discovery. In combination with the results of the PPI network and CytoHubba, 10 hub genes including CEACAM8, CD19, MMP9, ARG1, CKAP4, CCR7, MGAM, CD79A, CD79B, and CLEC4D were selected. Combined with DEG-miRNAs visualization, 5 miRNAs, including hsa-mir-146a-5p, hsa-mir-7-5p, hsa-mir-335-5p, and hsa-mir-27a- 3p, were predicted as possibly the key miRNAs. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of ischemic stroke, and provide a new strategy for clinical therapy.

Single-nucleotide-level mapping of DNA regulatory elements that control fetal hemoglobin expression

Pinpointing functional noncoding DNA sequences and defining their contributions to health-related traits is a major challenge for modern genetics. We developed a high-throughput framework to map noncoding DNA functions with single-nucleotide resolution in four loci that control erythroid fetal hemoglobin (HbF) expression, a genetically determined trait that modifies sickle cell disease (SCD) phenotypes. Specifically, we used the adenine base editor ABEmax to introduce 10,156 separate A•T to G•C conversions in 307 predicted regulatory elements and quantified the effects on erythroid HbF expression. We identified numerous regulatory elements, defined their epigenomic structures and linked them to low-frequency variants associated with HbF expression in an SCD cohort. Targeting a newly discovered γ-globin gene repressor element in SCD donor CD34+ hematopoietic progenitors raised HbF levels in the erythroid progeny, inhibiting hypoxia-induced sickling. Our findings reveal previously unappreciated genetic complexities of HbF regulation and provide potentially therapeutic insights into SCD.

Genome-Wide Role of HSF1 in Transcriptional Regulation of Desiccation Tolerance in the Anhydrobiotic Cell Line, Pv11

The Pv11, an insect cell line established from the midge Polypedilum vanderplanki, is capable of extreme hypometabolic desiccation tolerance, so-called anhydrobiosis. We previously discovered that heat shock factor 1 (HSF1) contributes to the acquisition of desiccation tolerance by Pv11 cells, but the mechanistic details have yet to be elucidated. Here, by analyzing the gene expression profiles of newly established HSF1-knockout and -rescue cell lines, we show that HSF1 has a genome-wide effect on gene regulation in Pv11. The HSF1-knockout cells exhibit a reduced desiccation survival rate, but this is completely restored in HSF1-rescue cells. By comparing mRNA profiles of the two cell lines, we reveal that HSF1 induces anhydrobiosis-related genes, especially genes encoding late embryogenesis abundant proteins and thioredoxins, but represses a group of genes involved in basal cellular processes, thus promoting an extreme hypometabolism state in the cell. In addition, HSF1 binding motifs are enriched in the promoters of anhydrobiosis-related genes and we demonstrate binding of HSF1 to these promoters by ChIP-qPCR. Thus, HSF1 directly regulates the transcription of anhydrobiosis-related genes and consequently plays a pivotal role in the induction of anhydrobiotic ability in Pv11 cells.

Panomics: New Databases for Advancing Cardiology

The multifactorial nature of cardiology makes it challenging to separate noisy signals from confounders and real markers or drivers of disease. Panomics, the combination of various omic methods, provides the deepest insights into the underlying biological mechanisms to develop tools for personalized medicine under a systems biology approach. Questions remain about current findings and anticipated developments of omics. Here, we search for omic databases, investigate the types of data they provide, and give some examples of panomic applications in health care. We identified 104 omic databases, of which 72 met the inclusion criteria: genomic and clinical measurements on a subset of the database population plus one or more omic datasets. Of those, 65 were methylomic, 59 transcriptomic, 41 proteomic, 42 metabolomic, and 22 microbiomic databases. Larger database sample sizes and longer follow-up are often better suited for panomic analyses due to statistical power calculations. They are often more complete, which is important when dealing with large biological variability. Thus, the UK BioBank rises as the most comprehensive panomic resource, at present, but certain study designs may benefit from other databases.

Comprehensive co-expression analysis reveals TMC8 as a prognostic immune-associated gene in head and neck squamous cancer

The occurrence and prognosis of head and neck squamous cell cancer (HNSC) is closely associated with human papillomavirus (HPV) infection. Transmembrane channel-like 8 (TMC8) is a key gene affecting the susceptibility of HPV and that plays an important role in T cell regulation. However, the mechanism by which TMC8 affects T cells and whether it further affects the prognosis of patients with HNSC remains unclear. In the present study, oral cancer cell lines and independent tumor specimens were used to detect TMC8 expression in HNSC. Differential expression of TMC8, methylation status, function and associated signaling pathways were further analyzed. Then, multiple databases were cross-analyzed for the relationship of TMC8 with immune cell infiltration and its impact on the prognosis of numerous types of cancer. The results showed that TMC8 was upregulated in HNSC and high expression was predictive of an improved prognosis. Furthermore, TMC8 was concentrated in multiple immune-associated signaling pathways and the expression of TMC8 was associated with the infiltration of CD4⁺ T cells and their subsets, including CD8⁺ T cells, B cells and macrophages, suggesting that TMC8 may play an anti-HPV role by regulating CD4⁺ T cells. Thus, TMC8 plays an anti-HPV role by regulating the infiltration level of CD4⁺ T cells, and could therefore be used as a potential prognostic marker for patients with HNSC.