Enhancer Variants Synergistically Drive Dysfunction of a Gene Regulatory Network In Hirschsprung Disease.

Common sequence variants in cis-regulatory elements (CREs) are suspected etiological causes of complex disorders. We previously identified an intronic enhancer variant in the RET gene disrupting SOX10 binding and increasing Hirschsprung disease (HSCR) risk 4-fold. We now show that two other functionally independent CRE variants, one binding Gata2 and the other binding Rarb, also reduce Ret expression and increase risk 2- and 1.7-fold. By studying human and mouse fetal gut tissues and cell lines, we demonstrate that reduced RET expression propagates throughout its gene regulatory network, exerting effects on both its positive and negative feedback components. We also provide evidence that the presence of a combination of CRE variants synergistically reduces RET expression and its effects throughout the GRN. These studies show how the effects of functionally independent non-coding variants in a coordinated gene regulatory network amplify their individually small effects, providing a model for complex disorders.

Distilling pathophysiology from complex disease genetics.

Technologies for genome-wide sequence interrogation have dramatically improved our ability to identify loci associated with complex human disease. However, a chasm remains between correlations and causality that stems, in part, from a limiting theoretical framework derived from Mendelian genetics and an incomplete understanding of disease physiology. Here we propose a set of criteria, akin to Koch's postulates for infectious disease, for assigning causality between genetic variants and human disease phenotypes.

The road ahead in genetics and genomics.

In celebration of the 20th anniversary of Nature Reviews Genetics, we asked 12 leading researchers to reflect on the key challenges and opportunities faced by the field of genetics and genomics. Keeping their particular research area in mind, they take stock of the current state of play and emphasize the work that remains to be done over the next few years so that, ultimately, the benefits of genetic and genomic research can be felt by everyone.

RET enhancer haplotype-dependent remodeling of the human fetal gut development program.

Hirschsprung disease (HSCR) is associated with deficiency of the receptor tyrosine kinase RET, resulting in loss of cells of the enteric nervous system (ENS) during fetal gut development. The major contribution to HSCR risk is from common sequence variants in RET enhancers with additional risk from rare coding variants in many genes. Here, we demonstrate that these RET enhancer variants specifically alter the human fetal gut development program through significant decreases in gene expression of RET, members of the RET-EDNRB gene regulatory network (GRN), other HSCR genes, with an altered transcriptome of 2,382 differentially expressed genes across diverse neuronal and mesenchymal functions. A parsimonious hypothesis for these results is that beyond RET's direct effect on its GRN, it also has a major role in enteric neural crest-derived cell (ENCDC) precursor proliferation, its deficiency reducing ENCDCs with relative expansion of non-ENCDC cells. Thus, genes reducing RET proliferative activity can potentially cause HSCR. One such class is the 23 RET-dependent transcription factors enriched in early gut development. We show that their knockdown in human neuroblastoma SK-N-SH cells reduces RET and/or EDNRB gene expression, expanding the RET-EDNRB GRN. The human embryos we studied had major remodeling of the gut transcriptome but were unlikely to have had HSCR: thus, genetic or epigenetic changes in addition to those in RET are required for aganglionosis.

Ret deficiency decreases neural crest progenitor proliferation and restricts fate potential during enteric nervous system development.

The receptor tyrosine kinase RET plays a critical role in the fate specification of enteric neural crest-derived cells (ENCDCs) during enteric nervous system (ENS) development. RET loss of function (LoF) is associated with Hirschsprung disease (HSCR), which is marked by aganglionosis of the gastrointestinal (GI) tract. Although the major phenotypic consequences and the underlying transcriptional changes from Ret LoF in the developing ENS have been described, cell type- and state-specific effects are unknown. We performed single-cell RNA sequencing on an enriched population of ENCDCs from the developing GI tract of Ret null heterozygous and homozygous mice at embryonic day (E)12.5 and E14.5. We demonstrate four significant findings: 1) Ret-expressing ENCDCs are a heterogeneous population comprising ENS progenitors as well as glial- and neuronal-committed cells; 2) neurons committed to a predominantly inhibitory motor neuron developmental trajectory are not produced under Ret LoF, leaving behind a mostly excitatory motor neuron developmental program; 3) expression patterns of HSCR-associated and Ret gene regulatory network genes are impacted by Ret LoF; and 4) Ret deficiency leads to precocious differentiation and reduction in the number of proliferating ENS precursors. Our results support a model in which Ret contributes to multiple distinct cellular phenotypes during development of the ENS, including the specification of inhibitory neuron subtypes, cell cycle dynamics of ENS progenitors, and the developmental timing of neuronal and glial commitment.

Quality assessment and refinement of chromatin accessibility data using a sequence-based predictive model.

Chromatin accessibility assays are central to the genome-wide identification of gene regulatory elements associated with transcriptional regulation. However, the data have highly variable quality arising from several biological and technical factors. To surmount this problem, we developed a sequence-based machine learning method to evaluate and refine chromatin accessibility data. Our framework, gapped k-mer SVM quality check (gkmQC), provides the quality metrics for a sample based on the prediction accuracy of the trained models. We tested 886 DNase-seq samples from the ENCODE/Roadmap projects to demonstrate that gkmQC can effectively identify "high-quality" (HQ) samples with low conventional quality scores owing to marginal read depths. Peaks identified in HQ samples are more accurately aligned at functional regulatory elements, show greater enrichment of regulatory elements harboring functional variants, and explain greater heritability of phenotypes from their relevant tissues. Moreover, gkmQC can optimize the peak-calling threshold to identify additional peaks, especially for rare cell types in single-cell chromatin accessibility data.

Interleukin-1 Receptor Antagonist Gene (IL1RN) Variants Modulate the Cytokine Release Syndrome and Mortality of COVID-19.

BACKGROUND: We examined effects of single-nucleotide variants (SNVs) of IL1RN, the gene encoding the anti-inflammatory interleukin 1 receptor antagonist (IL-1Ra), on the cytokine release syndrome (CRS) and mortality in patients with acute severe respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. METHODS: IL1RN CTA haplotypes formed from 3 SNVs (rs419598, rs315952, rs9005) and the individual SNVs were assessed for association with laboratory markers of inflammation and mortality. We studied 2589 patients hospitalized with SARS-CoV-2 between March 2020 and March 2021. RESULTS: Mortality was 15.3% and lower in women than men (13.1% vs 17.3%, P = .0003). Carriers of the CTA-1/2 IL1RN haplotypes exhibited decreased inflammatory markers and increased plasma IL-1Ra. Evaluation of the individual SNVs of the IL1RN, carriers of the rs419598 C/C SNV exhibited significantly reduced inflammatory biomarker levels and numerically lower mortality compared to the C/T-T/T genotype (10.0% vs 17.8%, P = .052) in men, with the most pronounced association observed in male patients ≤74 years old, whose mortality was reduced by 80% (3.1% vs 14.0%, P = .030). CONCLUSIONS: The IL1RN haplotype CTA and C/C variant of rs419598 are associated with attenuation of the CRS and decreased mortality in men with acute SARS-CoV-2 infection. The data suggest that the IL1RN pathway modulates the severity of coronavirus disease 2019 (COVID-19) via endogenous anti-inflammatory mechanisms.

Sequence-based correction of barcode bias in massively parallel reporter assays.

Massively parallel reporter assays (MPRAs) are a high-throughput method for evaluating in vitro activities of thousands of candidate cis-regulatory elements (CREs). In these assays, candidate sequences are cloned upstream or downstream from a reporter gene tagged by unique DNA sequences. However, tag sequences may themselves affect reporter gene expression and lead to major potential biases in the measured cis-regulatory activity. Here, we present a sequence-based method for correcting tag-sequence-specific effects and show that our method can significantly reduce this source of variation and improve the identification of functional regulatory variants by MPRAs. We also show that our model captures sequence features associated with post-transcriptional regulation of mRNA. Thus, this new method helps not only to improve detection of regulatory signals in MPRA experiments but also to design better MPRA protocols.

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.

Genome-wide pleiotropy analysis identifies novel blood pressure variants and improves its polygenic risk scores.

Systolic and diastolic blood pressure (S/DBP) are highly correlated modifiable risk factors for cardiovascular disease (CVD). We report here a bidirectional Mendelian Randomization (MR) and horizontal pleiotropy analysis of S/DBP summary statistics from the UK Biobank (UKB)-International Consortium for Blood Pressure (ICBP) (UKB-ICBP) BP genome-wide association study and construct a composite genetic risk score (GRS) by including pleiotropic variants. The composite GRS captures greater (1.11-3.26 fold) heritability for BP traits and increases (1.09- and 2.01-fold) Nagelkerke's R2 for hypertension and CVD. We replicated 118 novel BP horizontal pleiotropic variants including 18 novel BP loci using summary statistics from the Million Veteran Program (MVP) study. An additional 219 novel BP signals and 40 novel loci were identified after a meta-analysis of the UKB-ICBP and MVP summary statistics but without further independent replication. Our study provides further insight into BP regulation and provides a novel way to construct a GRS by including pleiotropic variants for other complex diseases.

Pathogenic alleles in microtubule, secretory granule and extracellular matrix-related genes in familial keratoconus.

Keratoconus is a common corneal defect with a complex genetic basis. By whole exome sequencing of affected members from 11 multiplex families of European ancestry, we identified 23 rare, heterozygous, potentially pathogenic variants in 8 genes. These include nonsynonymous single amino acid substitutions in HSPG2, EML6 and CENPF in two families each, and in NBEAL2, LRP1B, PIK3CG and MRGPRD in three families each; ITGAX had nonsynonymous single amino acid substitutions in two families and an indel with a base substitution producing a nonsense allele in the third family. Only HSPG2, EML6 and CENPF have been associated with ocular phenotypes previously. With the exception of MRGPRD and ITGAX, we detected the transcript and encoded protein of the remaining genes in the cornea and corneal cell cultures. Cultured stromal cells showed cytoplasmic punctate staining of NBEAL2, staining of the fibrillar cytoskeletal network by EML6, while CENPF localized to the basal body of primary cilia. We inhibited the expression of HSPG2, EML6, NBEAL2 and CENPF in stromal cell cultures and assayed for the expression of COL1A1 as a readout of corneal matrix production. An upregulation in COL1A1 after siRNA inhibition indicated their functional link to stromal cell biology. For ITGAX, encoding a leukocyte integrin, we assayed its level in the sera of 3 affected families compared with 10 unrelated controls to detect an increase in all affecteds. Our study identified genes that regulate the cytoskeleton, protein trafficking and secretion, barrier tissue function and response to injury and inflammation, as being relevant to keratoconus.

Analysis of putative cis-regulatory elements regulating blood pressure variation.

Hundreds of loci have been associated with blood pressure (BP) traits from many genome-wide association studies. We identified an enrichment of these loci in aorta and tibial artery expression quantitative trait loci in our previous work in ~100 000 Genetic Epidemiology Research on Aging study participants. In the present study, we sought to fine-map known loci and identify novel genes by determining putative regulatory regions for these and other tissues relevant to BP. We constructed maps of putative cis-regulatory elements (CREs) using publicly available open chromatin data for the heart, aorta and tibial arteries, and multiple kidney cell types. Variants within these regions may be evaluated quantitatively for their tissue- or cell-type-specific regulatory impact using deltaSVM functional scores, as described in our previous work. We aggregate variants within these putative CREs within 50 Kb of the start or end of 'expressed' genes in these tissues or cell types using public expression data and use deltaSVM scores as weights in the group-wise sequence kernel association test to identify candidates. We test for association with both BP traits and expression within these tissues or cell types of interest and identify the candidates MTHFR, C10orf32, CSK, NOV, ULK4, SDCCAG8, SCAMP5, RPP25, HDGFRP3, VPS37B and PPCDC. Additionally, we examined two known QT interval genes, SCN5A and NOS1AP, in the Atherosclerosis Risk in Communities Study, as a positive control, and observed the expected heart-specific effect. Thus, our method identifies variants and genes for further functional testing using tissue- or cell-type-specific putative regulatory information.

Molecular Genetic Anatomy and Risk Profile of Hirschsprung's Disease.

BACKGROUND: Hirschsprung's disease, or congenital aganglionosis, is a developmental disorder of the enteric nervous system and is the most common cause of intestinal obstruction in neonates and infants. The disease has more than 80% heritability, including significant associations with rare and common sequence variants in genes related to the enteric nervous system, as well as with monogenic and chromosomal syndromes. METHODS: We genotyped and exome-sequenced samples from 190 patients with Hirschsprung's disease to quantify the genetic burden in patients with this condition. DNA sequence variants, large copy-number variants, and karyotype variants in probands were considered to be pathogenic when they were significantly associated with Hirschsprung's disease or another neurodevelopmental disorder. Novel genes were confirmed by functional studies in the mouse and human embryonic gut and in zebrafish embryos. RESULTS: The presence of five or more variants in four noncoding elements defined a widespread risk of Hirschsprung's disease (48.4% of patients and 17.1% of controls; odds ratio, 4.54; 95% confidence interval [CI], 3.19 to 6.46). Rare coding variants in 24 genes that play roles in enteric neural-crest cell fate, 7 of which were novel, were also common (34.7% of patients and 5.0% of controls) and conferred a much greater risk than noncoding variants (odds ratio, 10.02; 95% CI, 6.45 to 15.58). Large copy-number variants, which were present in fewer patients (11.4%, as compared with 0.2% of controls), conferred the highest risk (odds ratio, 63.07; 95% CI, 36.75 to 108.25). At least one identifiable genetic risk factor was found in 72.1% of the patients, and at least 48.4% of patients had a structural or regulatory deficiency in the gene encoding receptor tyrosine kinase (RET). For individual patients, the estimated risk of Hirschsprung's disease ranged from 5.33 cases per 100,000 live births (approximately 1 per 18,800) to 8.38 per 1000 live births (approximately 1 per 120). CONCLUSIONS: Among the patients in our study, Hirschsprung's disease arose from common noncoding variants, rare coding variants, and copy-number variants affecting genes involved in enteric neural-crest cell fate that exacerbate the widespread genetic susceptibility associated with RET. For individual patients, the genotype-specific odds ratios varied by a factor of approximately 67, which provides a basis for risk stratification and genetic counseling. (Funded by the National Institutes of Health.).

Identifying Needs, Challenges, and Benefits Among Adults and Parents of Children With Hirschsprung Disease.

ABSTRACT: Adults with Hirschsprung disease (HD), unaffected parents of children with HD, and affected adults with an affected child completed a cross-sectional survey with open-ended questions about greatest needs at diagnosis and at current time, greatest challenges encountered, and any benefits of having HD or having a child with HD. In the 297 respondents, information and good medical care were common needs at diagnosis and at the time of survey, but the information needed evolved with time. Managing ongoing symptoms was a frequently cited need and challenge, along with managing medical care and the social and emotional impact of HD. Perceived benefits included empathy for others and new perspectives on life. The needs and challenges identified in this study can guide healthcare providers in discussions with families. Provision of information, recommendations, and referrals based on each individual family's needs can support families with HD throughout the lifecycle and facilitate adaptation.

Gene- and tissue-level interactions in normal gastrointestinal development and Hirschsprung disease.

The development of the gut from endodermal tissue to an organ with multiple distinct structures and functions occurs over a prolonged time during embryonic days E10.5-E14.5 in the mouse. During this process, one major event is innervation of the gut by enteric neural crest cells (ENCCs) to establish the enteric nervous system (ENS). To understand the molecular processes underpinning gut and ENS development, we generated RNA-sequencing profiles from wild-type mouse guts at E10.5, E12.5, and E14.5 from both sexes. We also generated these profiles from homozygous Ret null embryos, a model for Hirschsprung disease (HSCR), in which the ENS is absent. These data reveal 4 major features: 1) between E10.5 and E14.5 the developmental genetic programs change from expression of major transcription factors and its modifiers to genes controlling tissue (epithelium, muscle, endothelium) specialization; 2) the major effect of Ret is not only on ENCC differentiation to enteric neurons but also on the enteric mesenchyme and epithelium; 3) a muscle genetic program exerts significant effects on ENS development; and 4) sex differences in gut development profiles are minor. The genetic programs identified, and their changes across development, suggest that both cell autonomous and nonautonomous factors, and interactions between the different developing gut tissues, are important for normal ENS development and its disorders.

Multiple SCN5A variant enhancers modulate its cardiac gene expression and the QT interval.

The rationale for genome-wide association study (GWAS) results is sequence variation in cis-regulatory elements (CREs) modulating a target gene's expression as the major cause of trait variation. To understand the complete molecular landscape of one of these GWAS loci, we performed in vitro reporter screens in cardiomyocyte cell lines for CREs overlapping nearly all common variants associated with any of five independent QT interval (QTi)-associated GWAS hits at the SCN5A-SCN10A locus. We identified 13 causal CRE variants using allelic reporter activity, cardiomyocyte nuclear extract-based binding assays, overlap with human cardiac tissue DNaseI hypersensitive regions, and predicted impact of sequence variants on DNaseI sensitivity. Our analyses identified at least one high-confidence causal CRE variant for each of the five sentinel hits that could collectively predict SCN5A cardiac gene expression and QTi association. Although all 13 variants could explain SCN5A gene expression, the highest statistical significance was obtained with seven variants (inclusive of the five above). Thus, multiple, causal, mutually associated CRE variants can underlie GWAS signals.

A gene regulatory network explains RET-EDNRB epistasis in Hirschsprung disease.

Disruptions in gene regulatory networks (GRNs), driven by multiple deleterious variants, potentially underlie complex traits and diseases. Hirschsprung disease (HSCR), a multifactorial disorder of enteric nervous system (ENS) development, is associated with at least 24 genes and seven chromosomal loci, with RET and EDNRB as its major genes. We previously demonstrated that RET transcription in the ENS is controlled by an extensive GRN involving the transcription factors (TFs) RARB, GATA2 and SOX10 and other HSCR genes. We now demonstrate, using human and mouse cellular and animal models, that EDNRB is transcriptionally regulated in the ENS by GATA2, SOX10 and NKX2.5 TFs. Significantly, RET and EDNRB expression is regulated by their shared use of GATA2 and SOX10, and in turn, these TFs are controlled by EDNRB and RET in a dose-dependent manner. This study expands the ENS development GRN to include both RET and EDNRB, uncovers the mechanistic basis for RET-EDNRB epistasis and emphasizes how functionally different genes associated with a complex disorder can be united through a common GRN.

Human cardiac cis-regulatory elements, their cognate transcription factors, and regulatory DNA sequence variants.

Cis-regulatory elements (CRE), short DNA sequences through which transcription factors (TFs) exert regulatory control on gene expression, are postulated to be the major sites of causal sequence variation underlying the genetics of complex traits and diseases. We present integrative analyses, combining high-throughput genomic and epigenomic data with sequence-based computations, to identify the causal transcriptional components in a given tissue. We use data on adult human hearts to demonstrate that (1) sequence-based predictions detect numerous, active, tissue-specific CREs missed by experimental observations, (2) learned sequence features identify the cognate TFs, (3) CRE variants are specifically associated with cardiac gene expression, and (4) a significant fraction of the heritability of exemplar cardiac traits (QT interval, blood pressure, pulse rate) is attributable to these variants. This general systems approach can thus identify candidate causal variants and the components of gene regulatory networks (GRN) to enable understanding of the mechanisms of complex disorders on a tissue- or cell-type basis.

Magnitude of Mendelian versus complex inheritance of rare disorders.

In medical genetics, the vast majority of patients with a currently known genetic basis harbor a rare deleterious allele explaining its Mendelian inheritance. Increasingly, for these phenotypes, we recognize exceptions to Mendelian expectations from non-penetrance of clinical disease to significant inter-individual variation in clinical manifestations, likely reflecting the actions of additional modifier genes. Despite recent progress, we still remain ignorant about the molecular basis for many rare disorders presumed to be Mendelian. The molecular evidence increasingly suggests a role for multiple genes in some of these cases, but for how many? In this article, I discuss why equating a phenotype as Mendelian or complex may be short-sighted or even erroneous. As we learn more about the functions of the human genome with its genes in networks, we should view the phenotype of an individual patient as arising from his or her total genomic deleterious burden in a set of functionally inter-related genes affecting that phenotype. This can sometimes arise from deleterious allele(s) at a single gene (Mendelian inheritance) creating a specific biochemical deficiency (or excess) but could just as frequently arise from the cumulative effects of multiple disease alleles (complex inheritance) leading to the same biochemical deficiency (or excess).

Loss of δ-catenin function in severe autism.

Autism is a multifactorial neurodevelopmental disorder affecting more males than females; consequently, under a multifactorial genetic hypothesis, females are affected only when they cross a higher biological threshold. We hypothesize that deleterious variants at conserved residues are enriched in severely affected patients arising from female-enriched multiplex families with severe disease, enhancing the detection of key autism genes in modest numbers of cases. Here we show the use of this strategy by identifying missense and dosage sequence variants in the gene encoding the adhesive junction-associated δ-catenin protein (CTNND2) in female-enriched multiplex families and demonstrating their loss-of-function effect by functional analyses in zebrafish embryos and cultured hippocampal neurons from wild-type and Ctnnd2 null mouse embryos. Finally, through gene expression and network analyses, we highlight a critical role for CTNND2 in neuronal development and an intimate connection to chromatin biology. Our data contribute to the understanding of the genetic architecture of autism and suggest that genetic analyses of phenotypic extremes, such as female-enriched multiplex families, are of innate value in multifactorial disorders.