Proteome-scale prediction of molecular mechanisms underlying dominant genetic diseases

Many dominant genetic disorders result from protein-altering mutations, acting primarily through dominant-negative (DN), gain-of-function (GOF), and loss-of-function (LOF) mechanisms. Deciphering the mechanisms by which dominant diseases exert their effects is often experimentally challenging and resource intensive, but is essential for developing appropriate therapeutic approaches. Diseases that arise via a LOF mechanism are more amenable to be treated by conventional gene therapy, whereas DN and GOF mechanisms may require gene editing or targeting by small molecules. Moreover, pathogenic missense mutations that act via DN and GOF mechanisms are more difficult to identify than those that act via LOF using nearly all currently available variant effect predictors. Here, we introduce a tripartite statistical model made up of support vector machine binary classifiers trained to predict whether human protein coding genes are likely to be associated with DN, GOF, or LOF molecular disease mechanisms. We test the utility of the predictions by examining biologically and clinically meaningful properties known to be associated with the mechanisms. Our results strongly support that the models are able to generalise on unseen data and offer insight into the functional attributes of proteins associated with different mechanisms. We hope that our predictions will serve as a springboard for researchers studying novel variants and those of uncertain clinical significance, guiding variant interpretation strategies and experimental characterisation. Predictions for the human UniProt reference proteome are available at https://osf.io/z4dcp/.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

BACKGROUND

Cells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and conditions. How the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear.

RESULTS

We analyze existing bulk and single-cell transcriptomic datasets to uncover the prevalence of transcriptional adaptation in mammalian systems across diverse contexts and cell types. We perform regulon gene expression analyses of transcription factor target sets in both bulk and pooled single-cell genetic perturbation datasets. Our results reveal greater robustness in expression of regulons of transcription factors exhibiting transcriptional adaptation compared to those of transcription factors that do not. Stochastic mathematical modeling of minimal compensatory gene networks qualitatively recapitulates several aspects of transcriptional adaptation, including paralog upregulation and robustness to mutation. Combined with machine learning analysis of network features of interest, our framework offers potential explanations for which regulatory steps are most important for transcriptional adaptation.

CONCLUSIONS

Our integrative approach identifies several putative hits-genes demonstrating possible transcriptional adaptation-to follow-up on experimentally and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Inverse and Proportional Trans Modulation of Gene Expression in Human Aneuploidies

Genomic imbalance in aneuploidy is often detrimental to organisms. To gain insight into the molecular basis of aneuploidies in humans, we analyzed transcriptome data from several autosomal and sex chromosome aneuploidies. The results showed that in human aneuploid cells, genes located on unvaried chromosomes are inversely or proportionally trans-modulated, while a subset of genes on the varied chromosomes are compensated. Less genome-wide modulation is found for sex chromosome aneuploidy compared with autosomal aneuploidy due to X inactivation and the retention of dosage sensitive regulators on both sex chromosomes to limit the effective dosage change. We also found that lncRNA and mRNA can have different responses to aneuploidy. Furthermore, we analyzed the relationship between dosage-sensitive transcription factors and their targets, which illustrated the modulations and indicates genomic imbalance is related to stoichiometric changes in components of gene regulatory complexes.In summary, this study demonstrates the existence of trans-acting effects and compensation mechanisms in human aneuploidies and contributes to our understanding of gene expression regulation in unbalanced genomes and disease states.

Activation of the YAP/KLF5 transcriptional cascade in renal tubular cells aggravates kidney injury

The Hippo/YAP pathway plays a critical role in tissue homeostasis. Our previous work demonstrated that renal tubular YAP activation induced by double knockout (dKO) of the upstream Hippo kinases Mst1 and Mst2 promotes tubular injury and renal inflammation under basal conditions. However, the importance of tubular YAP activation remains to be established in injured kidneys in which many other injurious pathways are simultaneously activated. Here, we show that tubular YAP was already activated 6 h after unilateral ureteral obstruction (UUO). Tubular YAP deficiency greatly attenuated tubular cell overproliferation, tubular injury, and renal inflammation induced by UUO or cisplatin. YAP promoted the transcription of the transcription factor KLF5. Consistent with this, the elevated expression of KLF5 and its target genes in Mst1/2 dKO or UUO kidneys was blocked by ablation of Yap in tubular cells. Inhibition of KLF5 prevented tubular cell overproliferation, tubular injury, and renal inflammation in Mst1/2 dKO kidneys. Therefore, our results demonstrate that tubular YAP is a key player in kidney injury. YAP and KLF5 form a transcriptional cascade, where tubular YAP activation induced by kidney injury promotes KLF5 transcription. Activation of this cascade induces tubular cell overproliferation, tubular injury, and renal inflammation.

PRPF31-retinitis pigmentosa: Challenges and opportunities for clinical translation

Mutations in pre-mRNA processing factor 31 cause autosomal dominant retinitis pigmentosa (PRPF31-RP), for which there is currently no efficient treatment, making this disease a prime target for the development of novel therapeutic strategies. PRPF31-RP exhibits incomplete penetrance due to haploinsufficiency, in which reduced levels of gene expression from the mutated allele result in disease. A variety of model systems have been used in the investigation of disease etiology and therapy development. In this review, we discuss recent advances in both in vivo and in vitro model systems, evaluating their advantages and limitations in the context of therapy development for PRPF31-RP. Additionally, we describe the latest approaches for treatment, including AAV-mediated gene augmentation, genome editing, and late-stage therapies such as optogenetics, cell transplantation, and retinal prostheses.

An unsupervised deep learning framework for predicting human essential genes from population and functional genomic data

BACKGROUND

The ability to accurately predict essential genes intolerant to loss-of-function (LOF) mutations can dramatically improve the identification of disease-associated genes. Recently, there have been numerous computational methods developed to predict human essential genes from population genomic data. While the existing methods are highly predictive of essential genes of long length, they have limited power in pinpointing short essential genes due to the sparsity of polymorphisms in the human genome.

RESULTS

Motivated by the premise that population and functional genomic data may provide complementary evidence for gene essentiality, here we present an evolution-based deep learning model, DeepLOF, to predict essential genes in an unsupervised manner. Unlike previous population genetic methods, DeepLOF utilizes a novel deep learning framework to integrate both population and functional genomic data, allowing us to pinpoint short essential genes that can hardly be predicted from population genomic data alone. Compared with previous methods, DeepLOF shows unmatched performance in predicting ClinGen haploinsufficient genes, mouse essential genes, and essential genes in human cell lines. Notably, at a false positive rate of 5%, DeepLOF detects 50% more ClinGen haploinsufficient genes than previous methods. Furthermore, DeepLOF discovers 109 novel essential genes that are too short to be identified by previous methods.

CONCLUSION

The predictive power of DeepLOF shows that it is a compelling computational method to aid in the discovery of essential genes.

Buffering of genetic dominance by allele-specific protein complex assembly

Protein complex assembly often occurs while subunits are being translated, resulting in complexes whose subunits were translated from the same mRNA in an allele-specific manner. It has thus been hypothesized that such cotranslational assembly may counter the assembly-mediated dominant-negative effect, whereby co-assembly of mutant and wild-type subunits "poisons" complex activity. Here, we show that cotranslationally assembling subunits are much less likely to be associated with autosomal dominant relative to recessive disorders, and that subunits with dominant-negative disease mutations are significantly depleted in cotranslational assembly compared to those associated with loss-of-function mutations. We also find that complexes with known dominant-negative effects tend to expose their interfaces late during translation, lessening the likelihood of cotranslational assembly. Finally, by combining complex properties with other features, we trained a computational model for predicting proteins likely to be associated with non-loss-of-function disease mechanisms, which we believe will be of considerable utility for protein variant interpretation.

Precise modulation of transcription factor levels identifies features underlying dosage sensitivity

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

The spectrum of GATA2 deficiency syndrome

Inherited or de novo germ line heterozygous mutations in the gene encoding the transcription factor GATA2 lead to its deficiency. This results in a constellation of clinical features including nontuberculous mycobacterial, bacterial, fungal, and human papillomavirus infections, lymphedema, pulmonary alveolar proteinosis, and myelodysplasia. The onset, or even the presence, of disease is highly variable, even in kindreds with the identical mutation in GATA2. The clinical manifestations result from the loss of a multilineage progenitor that gives rise to B lymphocytes, monocytes, natural killer cells, and dendritic cells, leading to cytopenias of these lineages and subsequent infections. The bone marrow failure is typically characterized by hypocellularity. Dysplasia may either be absent or subtle but typically evolves into multilineage dysplasia with prominent dysmegakaryopoiesis, followed in some instances by progression to myeloid malignancies, specifically myelodysplastic syndrome, acute myelogenous leukemia, and chronic myelomonocytic leukemia. The latter 3 malignancies often occur in the setting of monosomy 7, trisomy 8, and acquired mutations in ASXL1 or in STAG2. Importantly, myeloid malignancy may represent the primary presentation of disease without recognition of other syndromic features. Allogeneic hematopoietic stem cell transplantation (HSCT) results in reversal of the phenotype. There remain important unanswered questions in GATA2 deficiency, including the following: (1) Why do some family members remain asymptomatic despite harboring deleterious mutations in GATA2? (2) What are the genetic changes that lead to myeloid progression? (3) What causes the apparent genetic anticipation? (4) What is the role of preemptive HSCT?

A Critical Review of Zebrafish Neurological Disease Models-1. The Premise: Neuroanatomical, Cellular and Genetic Homology and Experimental Tractability

The last decade has seen a dramatic rise in the number of genes linked to neurological disorders, necessitating new models to explore underlying mechanisms and to test potential therapies. Over a similar period, many laboratories adopted zebrafish as a tractable model for studying brain development, defining neural circuits and performing chemical screens. Here we discuss strengths and limitations of using the zebrafish system to model neurological disorders. The underlying premise for many disease models is the high degree of homology between human and zebrafish genes, coupled with the conserved vertebrate Bauplan and repertoire of neurochemical signaling molecules. Yet, we caution that important evolutionary divergences often limit the extent to which human symptoms can be modeled meaningfully in zebrafish. We outline advances in genetic technologies that allow human mutations to be reproduced faithfully in zebrafish. Together with methods that visualize the development and function of neuronal pathways at the single cell level, there is now an unprecedented opportunity to understand how disease-associated genetic changes disrupt neural circuits, a level of analysis that is ideally suited to uncovering pathogenic changes in human brain disorders.

A Novel Loss-of-function Mutation in MYBPC3 Causes Familial Hypertrophic Cardiomyopathy with Extreme Intrafamilial Phenotypic Heterogeneity

Cardiomyopathies are a heterogeneous group of diseases predominantly affecting the heart muscle and often lead to progressive heart failure-related disability or cardiovascular death. Hypertrophic cardiomyopathy (HCM) is a cardiac muscle disorder mostly caused by the mutations in genes encoding cardiac sarcomere. Germ-line mutations in MYBPC3 causes hypertrophic cardiomyopathy (HCM). However, most of the HCM associated MYBPC3 mutations were truncating mutations. Extreme phenotypic heterogeneity was observed among HCM patients with MYBPC3 mutations. In this study, we investigated a Chinese man who presented with HCM. Whole exome sequencing identified a novel heterozygous deletion (c.3781_3785delGAGGC) in exon 33 of the MYBPC3 in the proband. This heterozygous variant causes frameshift (p.Glu1261Thrfs*3), which predicted to form a truncated MYBPC3 protein. The proband's father also carries this variant in a heterozygous state while the proband's mother did not harbor this variant. Here, we report on a novel deletion in the MYBPC3 gene associated with HCM. We also highlight the importance of whole exome sequencing for molecular diagnosis for the patients with familial HCM.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Consequences of Chromosome Loss: Why Do Cells Need Each Chromosome Twice?

Aneuploidy is a cellular state with an unbalanced chromosome number that deviates from the usual euploid status. During evolution, elaborate cellular mechanisms have evolved to maintain the correct chromosome content over generations. The rare errors often lead to cell death, cell cycle arrest, or impaired proliferation. At the same time, aneuploidy can provide a growth advantage under selective conditions in a stressful, frequently changing environment. This is likely why aneuploidy is commonly found in cancer cells, where it correlates with malignancy, drug resistance, and poor prognosis. To understand this "aneuploidy paradox", model systems have been established and analyzed to investigate the consequences of aneuploidy. Most of the evidence to date has been based on models with chromosomes gains, but chromosome losses and recurrent monosomies can also be found in cancer. We summarize the current models of chromosome loss and our understanding of its consequences, particularly in comparison to chromosome gains.

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.

Loss of proximal tubular transcription factor Krüppel-like factor 15 exacerbates kidney injury through loss of fatty acid oxidation

Loss of fatty acid β-oxidation (FAO) in the proximal tubule is a critical mediator of acute kidney injury and eventual fibrosis. However, transcriptional mediators of FAO in proximal tubule injury remain understudied. Krüppel-like factor 15 (KLF15), a highly enriched zinc-finger transcription factor in the proximal tubule, was significantly reduced in proximal tubule cells after aristolochic acid I (AAI) treatment, a proximal tubule-specific injury model. Proximal tubule specific knockout of Klf15 exacerbated proximal tubule injury and kidney function decline compared to control mice during the active phase of AAI treatment, and after ischemia-reperfusion injury. Furthermore, along with worsening proximal tubule injury and kidney function decline, knockout mice exhibited increased kidney fibrosis as compared to control mice during the remodeling phase after AAI treatment. RNA-sequencing of kidney cortex demonstrated increased transcripts involved in immune system and integrin signaling pathways and decreased transcripts encompassing metabolic pathways, specifically FAO, and PPARα signaling, in knockout versus control mice after AAI treatment. In silico and experimental chromatin immunoprecipitation studies collectively demonstrated that KLF15 occupied the promoter region of key FAO genes, CPT1A and ACAA2, in close proximity to transcription factor PPARα binding sites. While the loss of Klf15 reduced the expression of Cpt1a and Acaa2 and led to compromised FAO, induction of KLF15 partially rescued loss of FAO in AAI-treated cells. Klf15, Ppara, Cpt1a, and Acaa2 expression was also decreased in other mouse kidney injury models. Tubulointerstitial KLF15 independently correlated with eGFR, PPARA and CPT1A appearance in expression arrays from human kidney biopsies. Thus, proximal tubule-specific loss of Klf15 exacerbates acute kidney injury and fibrosis, likely due to loss of interaction with PPARα leading to loss of FAO gene transcription.

Two Novel C-Terminus RUNX2 Mutations in Two Cleidocranial Dysplasia (CCD) Patients Impairing p53 Expression

Cleidocranial dysplasia (CCD), a dominantly inherited skeletal disease, is characterized by a variable phenotype ranging from dental alterations to severe skeletal defects. Either de novo or inherited mutations in the RUNX2 gene have been identified in most CCD patients. Transcription factor RUNX2, the osteogenic master gene, plays a central role in the commitment of mesenchymal stem cells to osteoblast lineage. With the aim to analyse the effects of RUNX2 mutations in CCD patients, we investigated RUNX2 gene expression and the osteogenic potential of two CCD patients' cells. In addition, with the aim to better understand how RUNX2 mutations interfere with osteogenic differentiation, we performed string analyses to identify proteins interacting with RUNX2 and analysed p53 expression levels. Our findings demonstrated for the first time that, in addition to the alteration of downstream gene expression, RUNX2 mutations impair p53 expression affecting osteogenic maturation. In conclusion, the present work provides new insights into the role of RUNX2 mutations in CCD patients and suggests that an in-depth analysis of the RUNX2-associated gene network may contribute to better understand the complex molecular and phenotypic alterations in mutant subjects.

Human genetic analyses of organelles highlight the nucleus in age-related trait heritability

Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a 'hallmark' of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be 'haplosufficient,' while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration.

Does Developmental Variability in the Number of Midbrain Dopamine Neurons Affect Individual Risk for Sporadic Parkinson's Disease?

Parkinson’s disease (PD) is a slowly progressing neurodegenerative disorder that is coupled to both widespread protein aggregation and to loss of substantia nigra dopamine (DA) neurons, resulting in a wide variety of motor and non-motor signs and symptoms. Recent findings suggest that the PD process is triggered several years before there is sufficient degeneration of DA neurons to cause onset of overt motor symptoms. According to this concept, the number of DA neurons present in the substantia nigra at birth could influence the time from the molecular triggering event until the clinical diagnosis with lower number of neurons at birth increasing the risk to develop the disease. Conversely, the risk for diagnosis would be reduced if the number of DA neurons is high at birth. In this commentary, we discuss the genetic and epigenetic factors that might influence the number of nigral DA neurons that each individual is born with and how these may be linked to PD risk.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

Human GATA2 mutations and hematologic disease: how many paths to pathogenesis?

The surge of human genetic information, enabled by increasingly facile and economically feasible genomic technologies, has accelerated discoveries on the relationship of germline genetic variation to hematologic diseases. For example, germline variation in GATA2, encoding a vital transcriptional regulator of multilineage hematopoiesis, creates a predisposition to bone marrow failure and acute myeloid leukemia termed GATA2 deficiency syndrome. More than 300 GATA2 variants representing missense, truncating, and noncoding enhancer mutations have been documented. Although these variants can diminish GATA2 expression and/or function, the functional ramifications of many variants are unknown. Studies using genetic rescue and knockin mouse systems have established that GATA2 mutations differentially affect molecular processes in distinct target genes and within a single target cell. Considering that target genes for a transcription factor can differ in sensitivity to altered levels of the factor, and transcriptional mechanisms are often cell type specific, the context-dependent consequences of GATA2 mutations in experimental systems portend the complex phenotypes and interindividual variation of GATA2 deficiency syndrome. This review documents GATA2 human genetics and the state of efforts to traverse from physiological insights to pathogenic mechanisms.

Insights into the pathogenicity of missense variants in the forkhead domain of FOX proteins underlying Mendelian disorders

Forkhead box (FOX) proteins are members of a conserved family of transcription factors. Pathogenic variants in FOX genes have been shown to be responsible for several human genetic diseases. Here, we have studied the molecular and structural features of germline pathogenic variants in seven FOX proteins involved in Mendelian disorders and compared them with those of variants present in the general population (gnomAD). Our study shows that the DNA-binding domain of FOX proteins is particularly sensitive to damaging variation, although some family members show greater mutational tolerance than others. Next, we set to demonstrate that this tolerance depends on the inheritance mode of FOX-linked disorders. Accordingly, genes whose variants underlie recessive conditions are supposed to have a greater tolerance to variation. This is what we found. As expected, variants responsible for disorders with a dominant inheritance pattern show a higher degree of pathogenicity compared to those segregating in the general population. Moreover, we show that pathogenic and likely pathogenic variants tend to affect mutually exclusive sites with respect to those reported in gnomAD. The former also tend to affect sites with lower solvent exposure and a higher degree of conservation. Our results show the value of using publicly available databases and bioinformatics to gain insights into the molecular and structural bases of disease-causing genetic variation.

An empirical bayesian approach for testing gene expression fold change and its application in detecting global dosage effects

We are motivated by biological studies intended to understand global gene expression fold change. Biologists have generally adopted a fixed cutoff to determine the significance of fold changes in gene expression studies (e.g. by using an observed fold change equal to two as a fixed threshold). Scientists can also use a t-test or a modified differential expression test to assess the significance of fold changes. However, these methods either fail to take advantage of the high dimensionality of gene expression data or fail to test fold change directly. Our research develops a new empirical Bayesian approach to substantially improve the power and accuracy of fold-change detection. Specifically, we more accurately estimate gene-wise error variation in the log of fold change. We then adopt a t-test with adjusted degrees of freedom for significance assessment. We apply our method to a dosage study in Arabidopsis and a Down syndrome study in humans to illustrate the utility of our approach. We also present a simulation study based on real datasets to demonstrate the accuracy of our method relative to error variance estimation and power in fold-change detection. Our developed R package with a detailed user manual is publicly available on GitHub at https://github.com/cuiyingbeicheng/Foldseq.

Identification of haploinsufficient genes from epigenomic data using deep forest

Haploinsufficiency, wherein a single allele is not enough to maintain normal functions, can lead to many diseases including cancers and neurodevelopmental disorders. Recently, computational methods for identifying haploinsufficiency have been developed. However, most of those computational methods suffer from study bias, experimental noise and instability, resulting in unsatisfactory identification of haploinsufficient genes. To address those challenges, we propose a deep forest model, called HaForest, to identify haploinsufficient genes. The multiscale scanning is proposed to extract local contextual representations from input features under Linear Discriminant Analysis. After that, the cascade forest structure is applied to obtain the concatenated features directly by integrating decision-tree-based forests. Meanwhile, to exploit the complex dependency structure among haploinsufficient genes, the LightGBM library is embedded into HaForest to reveal the highly expressive features. To validate the effectiveness of our method, we compared it to several computational methods and four deep learning algorithms on five epigenomic data sets. The results reveal that HaForest achieves superior performance over the other algorithms, demonstrating its unique and complementary performance in identifying haploinsufficient genes. The standalone tool is available at https://github.com/yangyn533/HaForest.

GmKIX8-1 regulates organ size in soybean and is the causative gene for the major seed weight QTL qSw17-1

Seed weight is one of the most important agronomic traits in soybean for yield improvement and food production. Several quantitative trait loci (QTLs) associated with the trait have been identified in soybean. However, the genes underlying the QTLs and their functions remain largely unknown. Using forward genetic methods and CRISPR/Cas9 gene editing, we identified and characterized the role of GmKIX8-1 in the control of organ size in soybean. GmKIX8-1 belongs to a family of KIX domain-containing proteins that negatively regulate cell proliferation in plants. Consistent with this predicted function, we found that loss-of-function GmKIX8-1 mutants showed a significant increase in the size of aerial plant organs, such as seeds and leaves. Likewise, the increase in organ size is due to increased cell proliferation, rather than cell expansion, and increased expression of CYCLIN D3;1-10. Lastly, molecular analysis of soybean germplasms harboring the qSw17-1 QTL for the big-seeded phenotype indicated that reduced expression of GmKIX8-1 is the genetic basis of the qSw17-1 phenotype.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

Functional effects of protein variants

Genetic and other variations frequently affect protein functions. Scientific articles can contain confusing descriptions about which function or property is affected, and in many cases the statements are pure speculation without any experimental evidence. To clarify functional effects of protein variations of genetic or non-genetic origin, a systematic conceptualisation and framework are introduced. This framework describes protein functional effects on abundance, activity, specificity and affinity, along with countermeasures, which allow cells, tissues and organisms to tolerate, avoid, repair, attenuate or resist (TARAR) the effects. Effects on abundance discussed include gene dosage, restricted expression, mis-localisation and degradation. Enzymopathies, effects on kinetics, allostery and regulation of protein activity are subtopics for the effects of variants on activity. Variation outcomes on specificity and affinity comprise promiscuity, specificity, affinity and moonlighting. TARAR mechanisms redress variations with active and passive processes including chaperones, redundancy, robustness, canalisation and metabolic and signalling rewiring. A framework for pragmatic protein function analysis and presentation is introduced. All of the mechanisms and effects are described along with representative examples, most often in relation to diseases. In addition, protein function is discussed from evolutionary point of view. Application of the presented framework facilitates unambiguous, detailed and specific description of functional effects and their systematic study.

Promoter CpG Density Predicts Downstream Gene Loss-of-Function Intolerance

The aggregation and joint analysis of large numbers of exome sequences has recently made it possible to derive estimates of intolerance to loss-of-function (LoF) variation for human genes. Here, we demonstrate strong and widespread coupling between genic LoF intolerance and promoter CpG density across the human genome. Genes downstream of the most CpG-rich promoters (top 10% CpG density) have a 67.2% probability of being highly LoF intolerant, using the LOEUF metric from gnomAD. This is in contrast to 7.4% of genes downstream of the most CpG-poor (bottom 10% CpG density) promoters. Combining promoter CpG density with exonic and promoter conservation explains 33.4% of the variation in LOEUF, and the contribution of CpG density exceeds the individual contributions of exonic and promoter conservation. We leverage this to train a simple and easily interpretable predictive model that outperforms other existing predictors and allows us to classify 1,760 genes-which are currently unascertained in gnomAD-as highly LoF intolerant or not. These predictions have the potential to aid in the interpretation of novel variants in the clinical setting. Moreover, our results reveal that high CpG density is not merely a generic feature of human promoters but is preferentially encountered at the promoters of the most selectively constrained genes, calling into question the prevailing view that CpG islands are not subject to selection.

Computational Approaches for Unraveling the Effects of Variation in the Human Genome and Microbiome

The past two decades of analytical efforts have highlighted how much more remains to be learned about the human genome and, particularly, its complex involvement in promoting disease development and progression. While numerous computational tools exist for the assessment of the functional and pathogenic effects of genome variants, their precision is far from satisfactory, particularly for clinical use. Accumulating evidence also suggests that the human microbiome's interaction with the human genome plays a critical role in determining health and disease states. While numerous microbial taxonomic groups and molecular functions of the human microbiome have been associated with disease, the reproducibility of these findings is lacking. The human microbiome–genome interaction in healthy individuals is even less well understood. This review summarizes the available computational methods built to analyze the effect of variation in the human genome and microbiome. We address the applicability and precision of these methods across their possible uses. We also briefly discuss the exciting, necessary, and now possible integration of the two types of data to improve the understanding of pathogenicity mechanisms.

Trait-associated noncoding variant regions affect TBX3 regulation and cardiac conduction

Genome-wide association studies have implicated common genomic variants in the gene desert upstream of TBX3 in cardiac conduction velocity. Whether these noncoding variants affect expression of TBX3 or neighboring genes and how they affect cardiac conduction is not understood. Here, we use high-throughput STARR-seq to test the entire 1.3 Mb human and mouse TBX3 locus, including two cardiac conduction-associated variant regions, for regulatory function. We identified multiple accessible and functional regulatory DNA elements that harbor variants affecting their activity. Both variant regions drove gene expression in the cardiac conduction tissue in transgenic reporter mice. Genomic deletion from the mouse genome of one of the regions caused increased cardiac expression of only Tbx3, PR interval shortening and increased QRS duration. Combined, our findings address the mechanistic link between trait-associated variants in the gene desert, TBX3 regulation and cardiac conduction.

Loss of the transcription factor MAFB limits β-cell derivation from human PSCs

Next generation sequencing studies have highlighted discrepancies in β-cells which exist between mice and men. Numerous reports have identified MAF BZIP Transcription Factor B (MAFB) to be present in human β-cells postnatally, while its expression is restricted to embryonic and neo-natal β-cells in mice. Using CRISPR/Cas9-mediated gene editing, coupled with endocrine cell differentiation strategies, we dissect the contribution of MAFB to β-cell development and function specifically in humans. Here we report that MAFB knockout hPSCs have normal pancreatic differentiation capacity up to the progenitor stage, but favor somatostatin- and pancreatic polypeptide–positive cells at the expense of insulin- and glucagon-producing cells during endocrine cell development. Our results describe a requirement for MAFB late in the human pancreatic developmental program and identify it as a distinguishing transcription factor within islet cell subtype specification. We propose that hPSCs represent a powerful tool to model human pancreatic endocrine development and associated disease pathophysiology.

The MAF bZIP transcription factor B (MAFB) is present in postnatal human beta cells but its role is unclear. Here, the authors show that MAFB regulates endocrine pancreatic cell fate specification.

Functional Characterization of Rare Variants in the SHOX2 Gene Identified in Sinus Node Dysfunction and Atrial Fibrillation

Sinus node dysfunction (SND) and atrial fibrillation (AF) often coexist; however, the molecular mechanisms linking both conditions remain elusive. Mutations in the homeobox-containing SHOX2 gene have been recently associated with early-onset and familial AF. Shox2 is a key regulator of sinus node development, and its deficiency leads to bradycardia, as demonstrated in animal models. To provide an extended SHOX2 gene analysis in patients with distinct arrhythmias, we investigated SHOX2 as a susceptibility gene for SND and AF by screening 98 SND patients and 450 individuals with AF. The functional relevance of the novel mutations was investigated in vivo and in vitro, together with the previously reported p.H283Q variant. A heterozygous missense mutation (p.P33R) was identified in the SND cohort and four heterozygous variants (p.G77D, p.L129=, p.L130F, p.A293=) in the AF cohort. Overexpression of the pathogenic predicted mutations in zebrafish revealed pericardial edema for p.G77D and the positive control p.H283Q, whereas the p.P33R and p.A293= variants showed no effect. In addition, a dominant-negative effect with reduced heart rates was detected for p.G77D and p.H283Q. In vitro reporter assays demonstrated for both missense variants p.P33R and p.G77D significantly impaired transactivation activity, similar to the described p.H283Q variant. Also, a reduced Bmp4 target gene expression was revealed in zebrafish hearts upon overexpression of the p.P33R mutant. This study associates additional rare variants in the SHOX2 gene implicated in the susceptibility to distinct arrhythmias and allows frequency estimations in the AF cohort (3/990). We also demonstrate for the first time a genetic link between SND and AF involving SHOX2. Moreover, our data highlight the importance of functional investigations of rare variants.

PUMILIO1 Links Epilepsy to Spinocerebellar Ataxia

A Mild PUM1 Mutation Is Associated With Adult-Onset Ataxia, Whereas Haploinsufficiency Causes Developmental Delay and Seizures

Gennarino VA, Palmer EE, McDonell LM, et al. Cell. 2018;172(5):924-936.e11. doi:10.1016/j.cell.2018.02.006.

Certain mutations can cause proteins to accumulate in neurons, leading to neurodegeneration. We recently showed, however, that upregulation of a wild-type protein, Ataxin1, caused by haploinsufficiency of its repressor, the RNA-binding protein Pumilio1 (PUM1), also causes neurodegeneration in mice. We therefore searched for human patients with PUM1 mutations. We identified 11 individuals with either PUM1 deletions or de novo missense variants who suffer a developmental syndrome (PUM1-associated developmental disability, ataxia, and seizure). We also identified a milder missense mutation in a family with adult-onset ataxia with incomplete penetrance (PUM1-related cerebellar ataxia). Studies in patient-derived cells revealed that the missense mutations reduced PUM1 protein levels by ∼25% in the adult-onset cases and by ∼50% in the infantile-onset cases; levels of known PUM1 targets increased accordingly. Changes in protein levels thus track with phenotypic severity, and identifying posttranscriptional modulators of protein expression should identify new candidate disease genes.

miR-17∼92 in lymphocyte development and lymphomagenesis

microRNAs (miRNAs) down-modulate the levels of proteins by sequence-specific binding to their respective target mRNAs, causing translational repression or mRNA degradation. The miR-17∼92 cluster encodes for six miRNAs whose target recognition specificities are determined by their distinct sequence. In mice, the four miRNA families generated from the miR-17∼92 cluster coordinate to allow for proper lymphocyte development and effective adaptive immune responses following infection or immunization. Lymphocyte development and homeostasis rely on tight regulation of PI3K signaling to avoid autoimmunity or immunodeficiency, and the miR-17∼92 miRNAs appear as key mediators to appropriately tune PI3K activity. On the other hand, in lymphoid tumors overexpression of the miR-17∼92 miRNAs is a common oncogenic event. In this review, we touch on what we have learned so far about the miR-17∼92 miRNAs, particularly with respect to their role in lymphocyte development, homeostasis and pathology.

C9orf72 Protein Plasmatic Concentrations Are Similar between C9ORF72 Expansion Carriers and Noncarriers in Frontotemporal Dementia

BACKGROUND/AIMS

The aim of the study was to assess the theory of haploinsufficiency in C9ORF72 expansion carriers, the most frequent causative gene of frontotemporal dementia.

METHODS

Plasmatic concentrations of C9orf72 protein were measured in 33 patients suspected of familial frontotemporal dementia using an enzyme-linked immunosorbent assay.

RESULTS

No difference was observed between C9ORF72 expansion carriers (21.2% of patients) and noncarriers (78.8% of patients). C9orf72 protein determination is not a suitable biomarker for screening C9ORF72 expansion carriers.

CONCLUSION

Our results provide new evidence against the hypothesis of haploinsufficiency leading to frontotemporal dementia in C9ORF72 expansion carriers.

Important cardiac transcription factor genes are accompanied by bidirectional long non-coding RNAs

Background

Heart development is a relatively fragile process in which many transcription factor genes show dose-sensitive characteristics such as haploinsufficiency and lower penetrance. Despite efforts to unravel the genetic mechanism for overcoming the fragility under normal conditions, our understanding still remains in its infancy. Recent studies on the regulatory mechanisms governing gene expression in mammals have revealed that long non-coding RNAs (lncRNAs) are important modulators at the transcriptional and translational levels. Based on the hypothesis that lncRNAs also play important roles in mouse heart development, we attempted to comprehensively identify lncRNAs by comparing the embryonic and adult mouse heart and brain.

Results

We have identified spliced lncRNAs that are expressed during development and found that lncRNAs that are expressed in the heart but not in the brain are located close to genes that are important for heart development. Furthermore, we found that many important cardiac transcription factor genes are located in close proximity to lncRNAs. Importantly, many of the lncRNAs are divergently transcribed from the promoter of these genes. Since the lncRNA divergently transcribed from Tbx5 is highly evolutionarily conserved, we focused on and analyzed the transcript. We found that this lncRNA exhibits a different expression pattern than that of Tbx5, and knockdown of this lncRNA leads to embryonic lethality.

Conclusion

These results suggest that spliced lncRNAs, particularly bidirectional lncRNAs, are essential regulators of mouse heart development, potentially through the regulation of neighboring transcription factor genes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5233-5) contains supplementary material, which is available to authorized users.

De novo variants in congenital diaphragmatic hernia identify MYRF as a new syndrome and reveal genetic overlaps with other developmental disorders

Congenital diaphragmatic hernia (CDH) is a severe birth defect that is often accompanied by other congenital anomalies. Previous exome sequencing studies for CDH have supported a role of de novo damaging variants but did not identify any recurrently mutated genes. To investigate further the genetics of CDH, we analyzed de novo coding variants in 362 proband-parent trios including 271 new trios reported in this study. We identified four unrelated individuals with damaging de novo variants in MYRF (P = 5.3x10(-8)), including one likely gene-disrupting (LGD) and three deleterious missense (D-mis) variants. Eight additional individuals with de novo LGD or missense variants were identified from our other genetic studies or from the literature. Common phenotypes of MYRF de novo variant carriers include CDH, congenital heart disease and genitourinary abnormalities, suggesting that it represents a novel syndrome. MYRF is a membrane associated transcriptional factor highly expressed in developing diaphragm and is depleted of LGD variants in the general population. All de novo missense variants aggregated in two functional protein domains. Analyzing the transcriptome of patient-derived diaphragm fibroblast cells suggest that disease associated variants abolish the transcription factor activity. Furthermore, we showed that the remaining genes with damaging variants in CDH significantly overlap with genes implicated in other developmental disorders. Gene expression patterns and patient phenotypes support pleiotropic effects of damaging variants in these genes on CDH and other developmental disorders. Finally, functional enrichment analysis implicates the disruption of regulation of gene expression, kinase activities, intra-cellular signaling, and cytoskeleton organization as pathogenic mechanisms in CDH.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Coupling of DNA Replication and Negative Feedback Controls Gene Expression for Cell-Fate Decisions

Cellular decision-making arises from the expression of genes along a regulatory cascade, which leads to a choice between distinct phenotypic states. DNA dosage variations, often introduced by replication, can significantly affect gene expression to ultimately bias decision outcomes. The bacteriophage lambda system has long served as a paradigm for cell-fate determination, yet the effect of DNA replication remains largely unknown. Here, through single-cell studies and mathematical modeling we show that DNA replication drastically boosts cI expression to allow lysogenic commitment by providing more templates. Conversely, expression of CII, the upstream regulator of cI, is surprisingly robust to DNA replication due to the negative autoregulation of the Cro repressor. Our study exemplifies how living organisms can not only utilize DNA replication for gene expression control but also implement mechanisms such as negative feedback to allow the expression of certain genes to be robust to dosage changes resulting from DNA replication.

Dosage-Dependent Expression Variation Suppressed on the Drosophila Male X Chromosome

DNA copy number variation is associated with many high phenotypic heterogeneity disorders. We systematically examined the impact of Drosophila melanogaster deletions on gene expression profiles to ask whether increased expression variability owing to reduced gene dose might underlie this phenotypic heterogeneity. Indeed, we found that one-dose genes have higher gene expression variability relative to two-dose genes. We then asked whether this increase in variability could be explained by intrinsic noise within cells due to stochastic biochemical events, or whether expression variability is due to extrinsic noise arising from more complex interactions. Our modeling showed that intrinsic gene expression noise averages at the organism level and thus cannot explain increased variation in one-dose gene expression. Interestingly, expression variability was related to the magnitude of expression compensation, suggesting that regulation, induced by gene dose reduction, is noisy. In a remarkable exception to this rule, the single X chromosome of males showed reduced expression variability, even compared with two-dose genes. Analysis of sex-transformed flies indicates that X expression variability is independent of the male differentiation program. Instead, we uncovered a correlation between occupancy of the chromatin-modifying protein encoded by males absent on the first (mof) and expression variability, linking noise suppression to the specialized X chromosome dosage compensation system. MOF occupancy on autosomes in both sexes also lowered transcriptional noise. Our results demonstrate that gene dose reduction can lead to heterogeneous responses, which are often noisy. This has implications for understanding gene network regulatory interactions and phenotypic heterogeneity. Additionally, chromatin modification appears to play a role in dampening transcriptional noise.

Mechanisms of Mendelian dominance

Genetic dominance has long been considered as a qualitative reflection of interallelic interactions. Dominance arises from many multiple sources whose unifying theme is the existence of non-linear relationships between the genotypic and phenotypic values. One of the clearest examples are dominant negative mutations (DNMs) in which a defective subunit poisons a macromolecular complex. Dominance can also be due to the presence of a heterozygous null allele, as is the case of haploinsufficiency. Dominance can also be influenced by epistatic (interloci) interactions. For instance, a pre-existing genetic variant can make possible the expression of a pathogenic variant in a seemingly "dominant" fashion. Such interactions, which can make an individual more or less sensitive to a particular pathogenic variant, will also be discussed here.

Decoding disease-causing mechanisms of missense mutations from supramolecular structures

The inheritance modes of pathogenic missense mutations are known to be highly associated with protein structures; recessive mutations are mainly observed in the buried region of protein structures, whereas dominant mutations are significantly enriched in the interfaces of molecular interactions. However, the differences in phenotypic impacts among various dominant mutations observed in individuals are not fully understood. In the present study, the functional effects of pathogenic missense mutations on three-dimensional macromolecular complex structures were explored in terms of dominant mutation types, namely, haploinsufficiency, dominant-negative, or toxic gain-of-function. The major types of dominant mutation were significantly associated with the different types of molecular interactions, such as protein-DNA, homo-oligomerization, or intramolecular domain-domain interactions, affected by mutations. The dominant-negative mutations were biased toward molecular interfaces for cognate protein or DNA. The haploinsufficiency mutations were enriched on the DNA interfaces. The gain-of-function mutations were localized to domain-domain interfaces. Our results demonstrate a novel use of macromolecular complex structures for predicting the disease-causing mechanisms through inheritance modes.

BACH2 immunodeficiency illustrates an association between super-enhancers and haploinsufficiency

The transcriptional programs that guide lymphocyte differentiation depend on the precise expression and timing of transcription factors (TFs). The TF BACH2 is essential for T and B lymphocytes and is associated with an archetypal super-enhancer (SE). Single-nucleotide variants in the BACH2 locus are associated with several autoimmune diseases, but BACH2 mutations that cause Mendelian monogenic primary immunodeficiency have not previously been identified. Here we describe a syndrome of BACH2-related immunodeficiency and autoimmunity (BRIDA) that results from BACH2 haploinsufficiency. Affected subjects had lymphocyte-maturation defects that caused immunoglobulin deficiency and intestinal inflammation. The mutations disrupted protein stability by interfering with homodimerization or by causing aggregation. We observed analogous lymphocyte defects in Bach2-heterozygous mice. More generally, we observed that genes that cause monogenic haploinsufficient diseases were substantially enriched for TFs and SE architecture. These findings reveal a previously unrecognized feature of SE architecture in Mendelian diseases of immunity: heterozygous mutations in SE-regulated genes identified by whole-exome/genome sequencing may have greater significance than previously recognized.

Estimating the selective effects of heterozygous protein-truncating variants from human exome data

The dispensability of individual genes for viability has interested generations of geneticists. For some genes it is essential to maintain two functional chromosomal copies, while others may tolerate the loss of one or both copies. Exome sequence data from 60,706 individuals provide sufficient observations of rare protein truncating variants (PTVs) to make genome-wide estimates of selection against heterozygous loss of gene function. The cumulative frequency of rare deleterious PTVs is primarily determined by the balance between incoming mutations and purifying selection rather than genetic drift. This enables the estimation of the genome-wide distribution of selection coefficients for heterozygous PTVs and corresponding Bayesian estimates for individual genes. The strength of selection can discriminate the severity, age of onset, and mode of inheritance in Mendelian exome sequencing cases. We find that genes under the strongest selection are enriched in embryonic lethal mouse knockouts, putatively cell-essential genes, Mendelian disease genes, and regulators of transcription. Screening by essentiality, we find a large set of genes under strong selection that likely have critical function but have not yet been extensively annotated in published literature.

Differential Sensitivity of Target Genes to Translational Repression by miR-17~92

MicroRNAs (miRNAs) are thought to exert their functions by modulating the expression of hundreds of target genes and each to a small degree, but it remains unclear how small changes in hundreds of target genes are translated into the specific function of a miRNA. Here, we conducted an integrated analysis of transcriptome and translatome of primary B cells from mutant mice expressing miR-17~92 at three different levels to address this issue. We found that target genes exhibit differential sensitivity to miRNA suppression and that only a small fraction of target genes are actually suppressed by a given concentration of miRNA under physiological conditions. Transgenic expression and deletion of the same miRNA gene regulate largely distinct sets of target genes. miR-17~92 controls target gene expression mainly through translational repression and 5’UTR plays an important role in regulating target gene sensitivity to miRNA suppression. These findings provide molecular insights into a model in which miRNAs exert their specific functions through a small number of key target genes.

MicroRNAs (miRNAs) are small RNAs encoded by our genome. Each miRNA binds hundreds of target mRNAs and performs specific functions. It is thought that miRNAs exert their function by reducing the expression of all these target genes and each to a small degree. However, these target genes often have very diverse functions. It has been unclear how small changes in hundreds of target genes with diverse functions are translated into the specific function of a miRNA. Here we take advantage of recent technical advances to globally examine the mRNA and protein levels of 868 target genes regulated by miR-17~92, the first oncogenic miRNA, in mutant mice with transgenic overexpression or deletion of this miRNA gene. We show that miR-17~92 regulates target gene expression mainly at the protein level, with little effect on mRNA. Surprisingly, only a small fraction of target genes respond to miR-17~92 expression changes. Further studies show that the sensitivity of target genes to miR-17~92 is determined by a non-coding region of target mRNA. Our findings demonstrate that not every target gene is equal, and suggest that the function of a miRNA is mediated by a small number of key target genes.

Genome Editing in hPSCs Reveals GATA6 Haploinsufficiency and a Genetic Interaction with GATA4 in Human Pancreatic Development

Human disease phenotypes associated with haploinsufficient gene requirements are often not recapitulated well in animal models. Here, we have investigated the association between human GATA6 haploinsufficiency and a wide range of clinical phenotypes that include neonatal and adult-onset diabetes using CRISPR (clustered regularly interspaced short palindromic repeat)/Cas9-mediated genome editing coupled with human pluripotent stem cell (hPSC) directed differentiation. We found that loss of one GATA6 allele specifically affects the differentiation of human pancreatic progenitors from the early PDX1+ stage to the more mature PDX1+NKX6.1+ stage, leading to impaired formation of glucose-responsive β-like cells. In addition to this GATA6 haploinsufficiency, we also identified dosage-sensitive requirements for GATA6 and GATA4 in the formation of both definitive endoderm and pancreatic progenitor cells. Our work expands the application of hPSCs from studying the impact of individual gene loci to investigation of multigenic human traits, and it establishes an approach for identifying genetic modifiers of human disease.

MLPA analysis in a cohort of patients with autism

Background

Autism is a global neurodevelopmental disorder which generally manifests during the first 2 years and continues throughout life, with a range of symptomatic variations. Epidemiological studies show an important role of genetic factors in autism and several susceptible regions and genes have been identified. The aim of our study was to validate a cost-effective set of commercial Multiplex Ligation dependent Probe Amplification (MLPA) and methylation specific multiplex ligation dependent probe amplification (MS-MLPA) test in autistic children refered by the neurodevelopmental center and autism unit of a Paediatric Hospital.

Results

In this study 150 unrelated children with autism spectrum disorders were analysed for copy number variation in specific regions of chromosomes 15, 16 and 22, using MLPA. All the patients had been previously studied by conventional karyotype and fluorescence in situ hybridization (FISH) analysis for 15(q11.2q13) and, with these techniques, four alterations were identified. The MLPA technique confirmed these four and identified further six alterations by the combined application of the two different panels.

Conclusions

Our data show that MLPA is a cost effective straightforward and rapid method for detection of imbalances in a clinically characterized population with autism. It contributes to strengthen the relationship between genotype and phenotype of children with autism, showing the clinical difference between deletions and duplications.

Gain-of-function mutations in GATA6 lead to atrial fibrillation

BACKGROUND

The genetic basis of atrial fibrillation (AF) and congenital heart disease remains incompletely understood.

OBJECTIVE

We sought to determine the causative mutation in a family with AF, atrial septal defects, and ventricular septal defects.

METHODS

We evaluated a pedigree with 16 family members, 1 with an atrial septal defect, 1 with a ventricular septal defect, and 3 with AF; we performed whole exome sequencing in 3 affected family members. Given that early-onset AF was prominent in the family, we then screened individuals with early-onset AF, defined as an age of onset <66 years, for mutations in GATA6. Variants were functionally characterized using reporter assays in a mammalian cell line.

RESULTS

Exome sequencing in 3 affected individuals identified a conserved mutation, R585L, in the transcription factor gene GATA6. In the Massachusetts General Hospital Atrial Fibrillation (MGH AF) Study, the mean age of AF onset was 47.1 ± 10.9 years; 79% of the participants were men; and there was no evidence of structural heart disease. We identified 3 GATA6 variants (P91S, A177T, and A543G). Using wild-type and mutant GATA6 constructs driving atrial natriuretic peptide promoter reporter, we found that 3 of the 4 variants had a marked upregulation of luciferase activity (R585L: 4.1-fold, P < .0001; P91S: 2.5-fold, P = .0002; A177T; 1.7-fold, P = .03). In addition, when co-overexpressed with GATA4 and MEF2C, GATA6 variants exhibited upregulation of the alpha myosin heavy chain and atrial natriuretic peptide reporter activity.

CONCLUSION

Overall, we found gain-of-function mutations in GATA6 in both a family with early-onset AF and atrioventricular septal defects as well as in a family with sporadic, early-onset AF.