Noncoding deletions reveal a gene that is critical for intestinal function

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.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. In this study, we examine the positioning and functional significance of human genes, grouped by their evolutionary age, within the 3D organization of the genome. We reveal that genes of different evolutionary origin have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, such associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young genes to ancient genes. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Nonsense mutation in the novel PERCC1 gene as a genetic cause of congenital diarrhea and enteropathy

Congenital diarrheas and enteropathies (CODEs) constitute a heterogeneous group of individually rare disorders manifesting with infantile-onset chronic diarrhea. Genomic deletions in chromosome 16, encompassing a sequence termed the 'intestine-critical region (ICR)', were recently identified as the cause of an autosomal recessive congenital enteropathy. The regulatory sequence within the ICR is flanked by an unannotated open reading frame termed PERCC1, which plays a role in enteroendocrine cell (EEC) function. We investigated two unrelated children with idiopathic congenital diarrhea requiring home parenteral nutrition attending the Irish Intestinal Failure Program. Currently 12 and 19-years old, these Irish male patients presented with watery diarrhea and hypernatremic dehydration in infancy. Probands were phenotyped by comprehensive clinical investigations, including endoscopic biopsies and serum gastrin level measurements. Following negative exome sequencing, PCR and Sanger sequencing of the entire coding region and intron boundaries of PERCC1 were performed for each proband and their parents. In both patients, serum gastrin levels were low and failed to increase following a meal challenge. While no deletions involving the ICR were detected, targeted sequencing of the PERCC1 gene revealed a shared homozygous c.390C > G stop gain variant. We report clinical and molecular findings in two unrelated patients harboring a shared homozygous variant in PERCC1, comprising the first description of a point mutation in this gene in association with CODE. That both parenteral nutrition dependent children with unexplained diarrhea at our institution harbored a PERCC1 mutation underscores the importance of its inclusion in exome sequencing interpretation.

The genetics of monogenic intestinal epithelial disorders

Monogenic intestinal epithelial disorders, also known as congenital diarrheas and enteropathies (CoDEs), are a group of rare diseases that result from mutations in genes that primarily affect intestinal epithelial cell function. Patients with CoDE disorders generally present with infantile-onset diarrhea and poor growth, and often require intensive fluid and nutritional management. CoDE disorders can be classified into several categories that relate to broad areas of epithelial function, structure, and development. The advent of accessible and low-cost genetic sequencing has accelerated discovery in the field with over 45 different genes now associated with CoDE disorders. Despite this increasing knowledge in the causal genetics of disease, the underlying cellular pathophysiology remains incompletely understood for many disorders. Consequently, clinical management options for CoDE disorders are currently limited and there is an urgent need for new and disorder-specific therapies. In this review, we provide a general overview of CoDE disorders, including a historical perspective of the field and relationship to other monogenic disorders of the intestine. We describe the genetics, clinical presentation, and known pathophysiology for specific disorders. Lastly, we describe the major challenges relating to CoDE disorders, briefly outline key areas that need further study, and provide a perspective on the future genetic and therapeutic landscape.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Identification of key genes for hypertrophic cardiomyopathy using integrated network analysis of differential lncRNA and gene expression

Objective

Hypertrophic cardiomyopathy (HCM) is a complex heterogeneous heart disease. Recent reports found that long non-coding RNAs (lncRNAs) play an important role in the progression of cardiovascular diseases. The present study aimed to identify the novel lncRNAs and messenger RNAs (mRNAs) and determine the key pathways involved in HCM.

Methods

The lncRNA and mRNA sequencing datasets of GSE68316 and GSE130036 were downloaded from the Gene Expression Omnibus (GEO) database. An integrated co-expression network analysis was conducted to identify differentially expressed lncRNAs and differentially expressed mRNAs in patients with HCM. Then, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were explored to identify the biological functions and signaling pathways of the co-expression network. Protein–protein interaction (PPI) and hub gene networks were constructed by using Cytoscape software. Plasma samples of patients with HCM and the GSE89714 dataset were used to validate the bioinformatics results.

Results

A total of 1,426 differentially expressed long non-coding RNAs (lncRNAs) and 1,715 differentially expressed mRNAs were obtained from GSE68316, of which 965 lncRNAs and 896 mRNAs were upregulated and 461 lncRNAs and 819 mRNAs were downregulated. A total of 469 differentially expressed lncRNAs and 2,407 differentially expressed mRNAs were screened from GSE130036, of which 183 lncRNAs and 1,283 mRNAs were upregulated and 286 lncRNAs and 1,124 mRNAs were downregulated. A co-expression network was constructed and contained 30 differentially expressed lncRNAs and 63 differentially expressed mRNAs, which were primarily involved in ‘G-protein beta/gamma-subunit complex binding,' ‘polyubiquitin modification-dependent protein binding,' ‘Apelin signaling pathway,' and ‘Wnt signaling pathway.' The 10 hub genes in the upregulated network [G Protein Subunit Alpha I2 (GNAI2), G Protein Subunit Alpha I1 (GNAI1), G Protein Subunit Alpha I3 (GNAI3), G Protein Subunit Gamma 2 (GNG2), G Protein Subunit Beta 1 (GNB1), G Protein Subunit Gamma 13 (GNG13), G Protein Subunit Gamma Transducin 1 (GNGT1), G Protein Subunit Gamma 12 (GNG12), AKT Serine/Threonine Kinase 1 (AKT1) and GNAS Complex Locus (GNAS)] and the 10 hub genes in the downregulated network [Nucleotide-Binding Oligomerization Domain Containing Protein 2 (NOD2), Receptor-Interacting Serine/Threonine Kinase 2 (RIPK2), Nucleotide-Binding Oligomerization Domain Containing Protein 1 (NOD1), Mitochondrial Antiviral Signaling Protein (MAVS), Autophagy Related 16-Like 1 (ATG16L1), Interferon Induced With Helicase C Domain 1 (IFIH1), Autophagy Related 5 (ATG5), TANK-Binding Kinase 1 (TBK1), Caspase Recruitment Domain Family Member 9 (CARD9), and von Willebrand factor (VWF)] were screened using cytoHubba. The expression of LA16c-312E8.2 and RP5-1160K1.3 in the plasma of patients with HCM was elevated, and the expression of the MIR22 host gene (MIR22HG) was decreased, which was consistent with our analysis, while the expression of LINC00324 and Small Nucleolar RNA Host Gene 12 (SNHG12) was not significantly different between the two groups. Verification analyses performed on GSE89714 showed the upregulated mRNAs of Chloride Voltage-Gated Channel 7 (CLCN7), N-Acetylglucosamine-1-Phosphate Transferase Subunit Gamma (GNPTG), Unk Like Zinc Finger (UNKL), Adenosine Monophosphate Deaminase 2 (AMPD2), GNAI3, WD Repeat Domain 81 (WDR81), and Serpin Family F Member 1 (SERPINF1) and downregulated mRNAs of TATA-Box Binding Protein Associated Factor 12 (TAF12) co-expressed with five crucial lncRNAs. Moreover, GNAI2, GNAI3, GNG12, and vWF were upregulated and GNAS was downregulated in the top 10 hub genes of upregulated and downregulated PPI networks.

Conclusion

These findings from integrative biological analysis of lncRNA-mRNA co-expression networks explored the key genes and pathways and provide new insights into the understanding of the mechanism and discovering new therapeutic targets for HCM. Three differentially expressed pivotal lncRNAs (LA16c-312E8.2, RP5-1160K1.3, and MIR22HG) in the co-expression network may serve as biomarkers and intervention targets for the diagnosis and treatment of HCM.

PERCC1, a new member of the Yap/TAZ/FAM181 transcriptional co-regulator family

Motivation

Disrupted PERCC1 gene expression causes an intractable congenital diarrhoea in infants. However, this gene's molecular mechanism is unknown and no homologous proteins have been reported.

Results

Our detailed evolutionary analysis of PERCC1 sequence reveals it to be a previously unappreciated member of the YAP/TAZ/FAM181 family of homologous transcriptional regulators. Like YAP and TAZ, PERCC1 likely interacts with DNA via binding to TEA/ATTS domain transcription factors (TEADs) using its conserved interface-2 and -3 sequences. We compare the expression patterns of PERCC1 with those of YAP, TAZ, TEADs. Our report provides the identification and first in-depth bioinformatic analysis of a YAP/TAZ homologue, and a likely new regulator of the YAP/TAZ-TEAD transcriptional complex.

Availability and implementation

The data underlying this article are available in UniProt Database.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Congenital enteropathies involving defects in enterocyte structure or differentiation

Congenital enteropathies (CE) are a group of rare inherited diseases with a typical onset early in life. They involve defects in enterocyte structure or differentiation. They can cause a severe condition of intestinal failure (IF). The diagnostic approach is based first on clinical presentation (consanguinity, prenatal expression, polyhydramnios, early neonatal onset, aspect of stools, persistence at bowel rest, associated extra-digestive manifestations….) and histo-pathological analyses. These rare intestinal diseases cause protracted diarrhea that might resolve, for a few, with a dietetic approach. However, protracted or permanent IF may require long term parenteral nutrition and, in limited cases, intestinal transplantation. With the progresses in both clinical nutrition and genetics, many of these CE are nowadays associated with recognized gene mutations. It improved our knowledge and the understanding in the patho-physiology of these diseases, thus, leading potentially to therapeutic perspectives. These review cover most of the early onset CE and excludes the immune related diarrhea.

Host genetic control of gut microbiome composition

The gut microbiome plays a significant role in health and disease, and there is mounting evidence indicating that the microbial composition is regulated in part by host genetics. Heritability estimates for microbial abundance in mice and humans range from (0.05-0.45), indicating that 5-45% of inter-individual variation can be explained by genetics. Through twin studies, genetic association studies, systems genetics, and genome-wide association studies (GWAS), hundreds of specific host genetic loci have been shown to associate with the abundance of discrete gut microbes. Using genetically engineered knock-out mice, at least 30 specific genes have now been validated as having specific effects on the microbiome. The relationships among of host genetics, microbiome composition, and abundance, and disease is now beginning to be unraveled through experiments designed to test causality. The genetic control of disease and its relationship to the microbiome can manifest in multiple ways. First, a genetic variant may directly cause the disease phenotype, resulting in an altered microbiome as a consequence of the disease phenotype. Second, a genetic variant may alter gene expression in the host, which in turn alters the microbiome, producing the disease phenotype. Finally, the genetic variant may alter the microbiome directly, which can result in the disease phenotype. In order to understand the processes that underlie the onset and progression of certain diseases, future research must take into account the relationship among host genetics, microbiome, and disease phenotype, and the resources needed to study these relationships.

Identification and characterization of conserved noncoding cis-regulatory elements that impact Mecp2 expression and neurological functions

In this study, Shao et al. investigated the transcriptional regulation of MeCP2, and identified six putative noncoding regulatory elements of Mecp2, two of which are conserved in humans. Their findings provide insight into transcriptional regulation of Mecp2/MECP2 and highlight genomic sites that could serve as diagnostic and therapeutic targets in Rett syndrome (RTT) and MECP2 duplication syndrome (MDS).

While changes in MeCP2 dosage cause Rett syndrome (RTT) and MECP2 duplication syndrome (MDS), its transcriptional regulation is poorly understood. Here, we identified six putative noncoding regulatory elements of Mecp2, two of which are conserved in humans. Upon deletion in mice and human iPSC-derived neurons, these elements altered RNA and protein levels in opposite directions and resulted in a subset of RTT- and MDS-like behavioral deficits in mice. Our discovery provides insight into transcriptional regulation of Mecp2/MECP2 and highlights genomic sites that could serve as diagnostic and therapeutic targets in RTT or MDS.

Chromatin regulatory dynamics of early human small intestinal development using a directed differentiation model

The establishment of the small intestinal (SI) lineage during human embryogenesis ensures functional integrity of the intestine after birth. The chromatin dynamics that drive SI lineage formation and regional patterning in humans are essentially unknown. To fill this knowledge void, we apply a cutting-edge genomic technology to a state-of-the-art human model of early SI development. Specifically, we leverage chromatin run-on sequencing (ChRO-seq) to define the landscape of active promoters, enhancers and gene bodies across distinct stages of directed differentiation of human pluripotent stem cells into SI spheroids with regional specification. Through comprehensive ChRO-seq analysis we identify candidate stage-specific chromatin activity states, novel markers and enhancer hotspots during the directed differentiation. Moreover, we propose a detailed transcriptional network associated with SI lineage formation or regional patterning. Our ChRO-seq analyses uncover a previously undescribed pattern of enhancer activity and transcription at HOX gene loci underlying SI regional patterning. We also validated this unique HOX dynamics by the analysis of single cell RNA-seq data from human fetal SI. Overall, the results lead to a new proposed working model for the regulatory underpinnings of human SI development, thereby adding a novel dimension to the literature that has relied almost exclusively on non-human models.

Translational genomics and beyond in bipolar disorder

Genome-wide association studies (GWAS) have revealed multiple genomic loci conferring risk of bipolar disorder (BD), providing hints for its underlying pathobiology. However, there are still remaining questions to answer. For example, discordance exists between BD heritability estimated with earlier epidemiological evidence and that calculated based on common GWAS variations. Where is the "missing heritability"? How can we explain the biology of the disease based on genetic findings? In this review, we summarize the accomplishments and limitations of current BD GWAS, and discuss potential reasons for the "missing heritability." In addition, progresses of research for the biological mechanisms underlying BD genetic risk using brain tissues, reprogrammed cells, and model animals are reviewed. While our knowledge of BD genetic basis is significantly promoted by these efforts, the complexities of gene regulation in the genome, the spatial-temporal heterogeneity during brain development, and the limitations of different experimental models should always be considered. Notably, several genes have been widely studied given their relatively well-characterized involvement in BD (e.g., CACAN1C and ANK3), and findings of these genes are summarized to both outline possible biological mechanisms of BD and describe examples of translating GWAS discoveries into the pathophysiology.

New Diagnostic Approaches for Undiagnosed Rare Genetic Diseases

Accurate diagnosis is the cornerstone of medicine; it is essential for informed care and promoting patient and family well-being. However, families with a rare genetic disease (RGD) often spend more than five years on a diagnostic odyssey of specialist visits and invasive testing that is lengthy, costly, and often futile, as 50% of patients do not receive a molecular diagnosis. The current diagnostic paradigm is not well designed for RGDs, especially for patients who remain undiagnosed after the initial set of investigations, and thus requires an expansion of approaches in the clinic. Leveraging opportunities to participate in research programs that utilize new technologies to understand RGDs is an important path forward for patients seeking a diagnosis. Given recent advancements in such technologies and international initiatives, the prospect of identifying a molecular diagnosis for all patients with RGDs has never been so attainable, but achieving this goal will require global cooperation at an unprecedented scale.

Apcdd1 is a dual BMP/Wnt inhibitor in the developing nervous system and skin

Animal development and homeostasis depend on precise temporal and spatial intercellular signaling. Components shared between signaling pathways, generally thought to decrease specificity, paradoxically can also provide a solution to pathway coordination. Here we show that the Bone Morphogenetic Protein (BMP) and Wnt signaling pathways share Apcdd1 as a common inhibitor and that Apcdd1 is a taxon-restricted gene with novel domains and signaling functions. Previously, we showed that Apcdd1 inhibits Wnt signaling (Shimomura et al., 2010), here we find that Apcdd1 potently inhibits BMP signaling in body axis formation and neural differentiation in chicken, frog, zebrafish. Furthermore, we find that Apcdd1 has an evolutionarily novel protein domain. Our results from experiments and modeling suggest that Apcdd1 may coordinate the outputs of two signaling pathways that are central to animal development and human disease.

Analysis of transcript-deleterious variants in Mendelian disorders: implications for RNA-based diagnostics

Background

At least 50% of patients with suspected Mendelian disorders remain undiagnosed after whole-exome sequencing (WES), and the extent to which non-coding variants that are not captured by WES contribute to this fraction is unclear. Whole transcriptome sequencing is a promising supplement to WES, although empirical data on the contribution of RNA analysis to the diagnosis of Mendelian diseases on a large scale are scarce.

Results

Here, we describe our experience with transcript-deleterious variants (TDVs) based on a cohort of 5647 families with suspected Mendelian diseases. We first interrogate all families for which the respective Mendelian phenotype could be mapped to a single locus to obtain an unbiased estimate of the contribution of TDVs at 18.9%. We examine the entire cohort and find that TDVs account for 15% of all “solved” cases. We compare the results of RT-PCR to in silico prediction. Definitive results from RT-PCR are obtained from blood-derived RNA for the overwhelming majority of variants (84.1%), and only a small minority (2.6%) fail analysis on all available RNA sources (blood-, skin fibroblast-, and urine renal epithelial cells-derived), which has important implications for the clinical application of RNA-seq. We also show that RNA analysis can establish the diagnosis in 13.5% of 155 patients who had received “negative” clinical WES reports. Finally, our data suggest a role for TDVs in modulating penetrance even in otherwise highly penetrant Mendelian disorders.

Conclusions

Our results provide much needed empirical data for the impending implementation of diagnostic RNA-seq in conjunction with genome sequencing.

Loss-of-function tolerance of enhancers in the human genome

Previous studies have surveyed the potential impact of loss-of-function (LoF) variants and identified LoF-tolerant protein-coding genes. However, the tolerance of human genomes to losing enhancers has not yet been evaluated. Here we present the catalog of LoF-tolerant enhancers using structural variants from whole-genome sequences. Using a conservative approach, we estimate that individual human genomes possess at least 28 LoF-tolerant enhancers on average. We assessed the properties of LoF-tolerant enhancers in a unified regulatory network constructed by integrating tissue-specific enhancers and gene-gene interactions. We find that LoF-tolerant enhancers tend to be more tissue-specific and regulate fewer and more dispensable genes relative to other enhancers. They are enriched in immune-related cells while enhancers with low LoF-tolerance are enriched in kidney and brain/neuronal stem cells. We developed a supervised learning approach to predict the LoF-tolerance of all enhancers, which achieved an area under the receiver operating characteristics curve (AUROC) of 98%. We predict 3,519 more enhancers would be likely tolerant to LoF and 129 enhancers that would have low LoF-tolerance. Our predictions are supported by a known set of disease enhancers and novel deletions from PacBio sequencing. The LoF-tolerance scores provided here will serve as an important reference for disease studies.

Enhancers are elements where transcription factors bind and regulate the expression of protein-coding genes. Although multiple previous studies have focused on which genes can tolerate loss-of-function (LoF), none has systematically evaluated the tolerance of all enhancers in the human genome to LoF. Individual studies have shown a broad range of phenotypic effects of enhancer LoF. The phenotypic effects of enhancer LoF likely fall into a spectrum where deletion of LoF-tolerant enhancers would not elicit substantial phenotypic impact, while some enhancers are likely to cause fitness defects when deleted. Here we report a systematic computational approach that uses machine learning and properties of enhancers in a unified human regulatory network with tissue-specific annotations to predict the LoF-tolerance of all enhancers identified in the human genome. The LoF-tolerance scores of enhancers provided in this study can significantly facilitate the interpretation and prioritization of non-coding sequence variants for disease and functional studies.

Congenital Diarrheal Syndromes

Congenital diarrheal disorders are heterogeneous conditions characterized by diarrhea with onset in the first years of life. They range from simple temporary conditions, such as cow's milk protein intolerance to irreversible complications, such as microvillous inclusion disease with significant morbidity and mortality. Advances in genomic medicine have improved our understanding of these disorders, leading to an ever-increasing list of identified causative genes. The diagnostic approach to these conditions consists of establishing the presence of diarrhea by detailed review of the history, followed by characterizing the composition of the diarrhea, the response to fasting, and with further specialized testing.