The molecular basis of genetic dominance

Yes, polygenic sex determination is a thing!

The process of sexual development in animals is modulated by a variety of mechanisms. Some species respond to environmental cues, while, in others, sex determination is thought to be controlled by a single 'master regulator' gene. However, many animals respond to a combination of environmental cues (e.g., temperature) and genetic factors (e.g., sex chromosomes). Even among species in which genetic factors predominate, there is a continuum between monofactorial and polygenic systems. The perception that polygenic systems are rare may result from experiments that lack the statistical power to detect multiple loci. Intellectual biases against the existence of polygenic sex determination (PSD) may further arise from misconceptions about the regulation of developmental processes and a misreading of theoretical results on the stability of polygenic systems of sex determination.

Late-Onset Slowly Progressing Cone/Macular Dystrophy in Patients With the Biallelic Hypomorphic Variant p.Arg1933Ter in RP1

Purpose

Homozygous hypomorphic variants of the RP1 gene, including c.5797C>T, p.Arg1933Ter, have traditionally been considered non-pathogenic. This study aimed to elucidate the clinical manifestations of late-onset, slowly progressive cone/macular dystrophy in patients homozygous for p.Arg1933Ter in the RP1 gene.

Methods

Five patients with biallelic p.Arg1933Ter in RP1 were retrospectively recruited, and their clinical profiles were analyzed. Copy number variation analysis and Alu insertion assessment of genes associated with inherited retinal diseases were conducted. The results of comprehensive ophthalmological examinations, multimodal imaging, and full-field electroretinogram tests were analyzed.

Results

No specific sequencing errors or structural variations associated with the clinical phenotypes were identified. Alu element insertion in RP1 was not detected. The mean ± SD age at the first visit was 62.2 ± 9.8 years, with symptoms typically starting between 45 and 50 years of age. Two patients exhibited a mild form of cone/macular dystrophy, characterized by a relatively preserved fundus appearance and blurring of the ellipsoid zone on optical coherence tomography. Three patients had late-onset cone/macular dystrophy with significant atrophy.

Conclusions

To our knowledge, this study is the first to report that a homozygous hypomorphic variant of RP1, previously considered non-pathogenic, leads to cone/macular dystrophy.

Translational Relevance

The study introduces novel possibilities suggesting that the homozygous hypomorphic variant of RP1 may be linked to variant pathogenicity.

Ancestral role of Pax6 in chordate brain regionalization

The Pax6 gene is essential for eye and brain development across various animal species. Here, we investigate the function of Pax6 in the development of the anterior central nervous system (CNS) of the invertebrate chordate amphioxus using CRISPR/Cas9-induced genome editing. Specifically, we examined Pax6 mutants featuring a 6 bp deletion encompassing two invariant amino acids in the conserved paired domain, hypothesized to impair Pax6 DNA-binding capacity and gene regulatory functions. Although this mutation did not result in gross morphological changes in amphioxus larvae, it demonstrated a reduced ability to activate Pax6-responsive reporter gene, suggesting a hypomorphic effect. Expression analysis in mutant larvae revealed changes in gene expression within the anterior CNS, supporting the conserved role of Pax6 gene in brain regionalization across chordates. Additionally, our findings lend support to the hypothesis of a zona limitans intrathalamica (ZLI)-like region in amphioxus, suggesting evolutionary continuity in brain patterning mechanisms. ZLI region, found in both hemichordates and vertebrates, functions as a key signaling center and serves as a restrictive boundary between major thalamic regions.

Identification of potential candidate genes for the Huoyan trait in developing Wulong goose embryos by transcriptomic analysis

1. The Wulong goose is a Chinese breed and a source of high-quality meat and eggs. A characteristic of the Wulong goose is that a proportion of the birds do not have eyelids, known as the Huoyon trait.2. Wulong geese exhibiting the Huoyan trait at embryonic stages of 9 days (E9), 12 days (E12) and 14 days (E14) were selected alongside those with normal eyelids for comprehensive transcriptome sequencing. Differentially expressed gene (DEG) and functional enrichment analyses were performed and finally, eight DEG were chosen to verify the accuracy of qPCR sequencing.3. Overall, 466, 962 and 550 DEG were obtained from the three control groups, D9 vs. N9, D12 vs. N12 and D14 vs. N14, respectively, by differential analysis (p < 0.05). CDKN1C, CRH, CROCC and TYSND1 were significantly expressed in the three groups. Enrichment analysis revealed the enrichment of CROCC and TYSND1 in pathways of cell cycle process, endocytosis, microtubule-based process, microtubule organising centre organisation, protein processing and protein maturation. CDKN1C and CRH were enriched in the cell cycle and cAMP signalling pathway.4. Some collagen family genes were detected among the DEGs, including COL3A1, COL4A5, COL4A2 and COL4A1. FREM1 and FREM2 genes were detected in both Huoyan and normal eyelids. There was a significant difference (p < 0.01) in FREM1 expression between ED9 and ED14 in female embryos, but this difference was not observed in male embryos.

Hypomorphic variants in inherited retinal and ocular diseases: A review of the literature with clinical cases

Hypomorphic variants decrease, but do not eliminate, gene function via a reduction in the amount of mRNA or protein product produced by a gene or by production of a gene product with reduced function. Many hypomorphic variants have been implicated in inherited retinal diseases (IRDs) and other genetic ocular conditions; however, there is heterogeneity in the use of the term "hypomorphic" in the scientific literature. We searched for all hypomorphic variants reported to cause IRDs and ocular disorders. We also discuss the presence of hypomorphic variants in the patient population of our ocular genetics department over the past decade. We propose that standardized criteria should be adopted for use of the term "hypomorphic" to describe gene variants to improve genetic counseling and patient care outcomes.

Penetrance interactions of colour pattern loci in the African Monarch and their implications for the evolution of dominance

Scoring the penetrance of heterozygotes in complex phenotypes, like colour pattern, is difficult and complicates the analysis of systems in which dominance is incomplete or evolving. The African Monarch (Danaus chrysippus) represents an example where colour pattern heterozygotes, formed in the contact zone between the different subspecies, show such intermediate dominance. Colour pattern in this aposematic butterfly is controlled by three loci A, B and C. The B and C loci are closely linked in a B/C supergene and significant interaction of B and C phenotypes is therefore expected via linkage alone. The A locus, however, is not linked to B/C and is found on a different chromosome. To study interactions between these loci we generated colour pattern heterozygotes by crossing males and females bearing different A and B/C genotypes, collected from different parts of Africa. We derived a novel scoring system for the expressivity of the heterozygotes and, as predicted, we found significant interactions between the genotypes of the closely linked B and C loci. Surprisingly, however, we also found highly significant interactions between C and the unlinked A locus, modifications that generally increased the resemblance of heterozygotes to homozygous ancestors. In contrast, we found no difference in the penetrance of any of the corresponding heterozygotes from crosses conducted either in allopatry or sympatry, in reciprocal crosses of males and females, or in the presence or absence of endosymbiont mediated male-killing or its associated neoW mediated sex-linkage of colour pattern. Together, this data supports the idea that the different colour morphs of the African Monarch meet transiently in the East African contact zone and that genetic modifiers act to mask inappropriate expression of colour patterns in the incorrect environments.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Dominance is common in mammals and is associated with trans-acting gene expression and alternative splicing

BACKGROUND

Dominance and other non-additive genetic effects arise from the interaction between alleles, and historically these phenomena play a major role in quantitative genetics. However, most genome-wide association studies (GWAS) assume alleles act additively.

RESULTS

We systematically investigate both dominance-here representing any non-additive within-locus interaction-and additivity across 574 physiological and gene expression traits in three mammalian stocks: F2 intercross pigs, rat heterogeneous stock, and mice heterogeneous stock. Dominance accounts for about one quarter of heritable variance across all physiological traits in all species. Hematological and immunological traits exhibit the highest dominance variance, possibly reflecting balancing selection in response to pathogens. Although most quantitative trait loci (QTLs) are detectable as additive QTLs, we identify 154, 64, and 62 novel dominance QTLs in pigs, rats, and mice respectively that are undetectable as additive QTLs. Similarly, even though most cis-acting expression QTLs are additive, gene expression exhibits a large fraction of dominance variance, and trans-acting eQTLs are enriched for dominance. Genes causal for dominance physiological QTLs are less likely to be physically linked to their QTLs but instead act via trans-acting dominance eQTLs. In addition, thousands of eQTLs are associated with alternatively spliced isoforms with complex additive and dominant architectures in heterogeneous stock rats, suggesting a possible mechanism for dominance.

CONCLUSIONS

Although heritability is predominantly additive, many mammalian genetic effects are dominant and likely arise through distinct mechanisms. It is therefore advantageous to consider both additive and dominance effects in GWAS to improve power and uncover causality.

A single amino acid substitution in MdLAZY1A dominantly impairs shoot gravitropism in Malus

Plant architecture is 1 of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is 1 of the 4 paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that colocalizes with Region III, 1 of the 5 conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Overexpressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Dominance vs epistasis: the biophysical origins and plasticity of genetic interactions within and between alleles

An important challenge in genetics, evolution and biotechnology is to understand and predict how mutations combine to alter phenotypes, including molecular activities, fitness and disease. In diploids, mutations in a gene can combine on the same chromosome or on different chromosomes as a "heteroallelic combination". However, a direct comparison of the extent, sign, and stability of the genetic interactions between variants within and between alleles is lacking. Here we use thermodynamic models of protein folding and ligand-binding to show that interactions between mutations within and between alleles are expected in even very simple biophysical systems. Protein folding alone generates within-allele interactions and a single molecular interaction is sufficient to cause between-allele interactions and dominance. These interactions change differently, quantitatively and qualitatively as a system becomes more complex. Altering the concentration of a ligand can, for example, switch alleles from dominant to recessive. Our results show that intra-molecular epistasis and dominance should be widely expected in even the simplest biological systems but also reinforce the view that they are plastic system properties and so a formidable challenge to predict. Accurate prediction of both intra-molecular epistasis and dominance will require either detailed mechanistic understanding and experimental parameterization or brute-force measurement and learning.

Stretch-activated ion channel TMEM63B associates with developmental and epileptic encephalopathies and progressive neurodegeneration

By converting physical forces into electrical signals or triggering intracellular cascades, stretch-activated ion channels allow the cell to respond to osmotic and mechanical stress. Knowledge of the pathophysiological mechanisms underlying associations of stretch-activated ion channels with human disease is limited. Here, we describe 17 unrelated individuals with severe early-onset developmental and epileptic encephalopathy (DEE), intellectual disability, and severe motor and cortical visual impairment associated with progressive neurodegenerative brain changes carrying ten distinct heterozygous variants of TMEM63B, encoding for a highly conserved stretch-activated ion channel. The variants occurred de novo in 16/17 individuals for whom parental DNA was available and either missense, including the recurrent p.Val44Met in 7/17 individuals, or in-frame, all affecting conserved residues located in transmembrane regions of the protein. In 12 individuals, hematological abnormalities co-occurred, such as macrocytosis and hemolysis, requiring blood transfusions in some. We modeled six variants (p.Val44Met, p.Arg433His, p.Thr481Asn, p.Gly580Ser, p.Arg660Thr, and p.Phe697Leu), each affecting a distinct transmembrane domain of the channel, in transfected Neuro2a cells and demonstrated inward leak cation currents across the mutated channel even in isotonic conditions, while the response to hypo-osmotic challenge was impaired, as were the Ca2+ transients generated under hypo-osmotic stimulation. Ectopic expression of the p.Val44Met and p.Gly580Cys variants in Drosophila resulted in early death. TMEM63B-associated DEE represents a recognizable clinicopathological entity in which altered cation conductivity results in a severe neurological phenotype with progressive brain damage and early-onset epilepsy associated with hematological abnormalities in most individuals.

Spinocerebellar ataxia 38: structure-function analysis shows ELOVL5 G230V is proteotoxic, conformationally altered and a mutational hotspot

Fatty acid elongase ELOVL5 is part of a protein family of multipass transmembrane proteins that reside in the endoplasmic reticulum where they regulate long-chain fatty acid elongation. A missense variant (c.689G>T p.Gly230Val) in ELOVL5 causes Spinocerebellar Ataxia subtype 38 (SCA38), a neurodegenerative disorder characterized by autosomal dominant inheritance, cerebellar Purkinje cell demise and adult-onset ataxia. Having previously showed aberrant accumulation of p.G230V in the Golgi complex, here we further investigated the pathogenic mechanisms triggered by p.G230V, integrating functional studies with bioinformatic analyses of protein sequence and structure. Biochemical analysis showed that p.G230V enzymatic activity was normal. In contrast, SCA38-derived fibroblasts showed reduced expression of ELOVL5, Golgi complex enlargement and increased proteasomal degradation with respect to controls. By heterologous overexpression, p.G230V was significantly more active than wild-type ELOVL5 in triggering the unfolded protein response and in decreasing viability in mouse cortical neurons. By homology modelling, we generated native and p.G230V protein structures whose superposition revealed a shift in Loop 6 in p.G230V that altered a highly conserved intramolecular disulphide bond. The conformation of this bond, connecting Loop 2 and Loop 6, appears to be elongase-specific. Alteration of this intramolecular interaction was also observed when comparing wild-type ELOVL4 and the p.W246G variant which causes SCA34. We demonstrate by sequence and structure analyses that ELOVL5 p.G230V and ELOVL4 p.W246G are position-equivalent missense variants. We conclude that SCA38 is a conformational disease and propose combined loss of function by mislocalization and gain of toxic function by ER/Golgi stress as early events in SCA38 pathogenesis.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Dominance mechanisms in supergene alleles controlling butterfly wing pattern variation: insights from gene expression in Heliconius numata

Loci under balancing selection, where multiple alleles are maintained, offer a relevant opportunity to investigate the role of natural selection in shaping genetic dominance: the high frequency of heterozygotes at these loci has been shown to enable the evolution of dominance among alleles. In the butterfly Heliconius numata, mimetic wing color variations are controlled by an inversion polymorphism of a circa 2 Mb genomic region (supergene P), with strong dominance between sympatric alleles. To test how differences in dominance observed on wing patterns correlate with variations in expression levels throughout the supergene region, we sequenced the complete transcriptome of heterozygotes at the prepupal stage and compared it to corresponding homozygotes. By defining dominance based on non-overlapping ranges of transcript expression between genotypes, we found contrasting patterns of dominance between the supergene and the rest of the genome; the patterns of transcript expression in the heterozygotes were more similar to the expression observed in the dominant homozygotes in the supergene region. Dominance also differed among the three subinversions of the supergene, suggesting possible epistatic interactions among their gene contents underlying dominance evolution. We found the expression pattern of the melanization gene cortex located in the P-region to predict wing pattern phenotype in the heterozygote. We also identify new candidate genes that are potentially involved in mimetic color pattern variations highlighting the relevance of transcriptomic analyses in heterozygotes to pinpoint candidate genes in non-recombining regions.

Therapeutic landscape for inherited ocular diseases: Current and emerging therapies

Inherited ocular diseases comprise a heterogeneous group of rare and complex diseases, including inherited retinal diseases (IRDs) and inherited optic neuropathies. Recent success in adeno-associated virus-based gene therapy, voretigene neparvovec (Luxturna^®) for RPE65-related IRDs, has heralded rapid evolution in gene therapy platform technologies and strategies, from gene augmentation to RNA editing, as well as gene agnostic approaches such as optogenetics. This review discusses the fundamentals underlying the mode of inheritance, natural history studies and clinical trial outcomes, as well as current and emerging therapies covering gene therapy strategies, cell-based therapies and bionic vision.

Mono- and biallelic variant effects on disease at biobank scale

Identifying causal factors for Mendelian and common diseases is an ongoing challenge in medical genetics1. Population bottleneck events, such as those that occurred in the history of the Finnish population, enrich some homozygous variants to higher frequencies, which facilitates the identification of variants that cause diseases with recessive inheritance2,3. Here we examine the homozygous and heterozygous effects of 44,370 coding variants on 2,444 disease phenotypes using data from the nationwide electronic health records of 176,899 Finnish individuals. We find associations for homozygous genotypes across a broad spectrum of phenotypes, including known associations with retinal dystrophy and novel associations with adult-onset cataract and female infertility. Of the recessive disease associations that we identify, 13 out of 20 would have been missed by the additive model that is typically used in genome-wide association studies. We use these results to find many known Mendelian variants whose inheritance cannot be adequately described by a conventional definition of dominant or recessive. In particular, we find variants that are known to cause diseases with recessive inheritance with significant heterozygous phenotypic effects. Similarly, we find presumed benign variants with disease effects. Our results show how biobanks, particularly in founder populations, can broaden our understanding of complex dosage effects of Mendelian variants on disease.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

Mutagenesis alters sperm swimming velocity in Astyanax cave fish

We investigated the hypothesis that intra ejaculate sperm competition screens against the transmission of deleterious alleles, including new mutants, from male parent to offspring. Recent investigations have established that sperm haploid genotypes can have major effects on sperm traits such as cellular robustness, longevity, and fertilization success. However, there is no evidence that new mutations can meaningfully affect sperm phenotypes. We tested this directly by comparing sperm from mutagenized and non-mutagenized control males in Astyanax fish. We used N-ethyl-N-nitrosourea (ENU) to induce single base substitutions in spermatogonial stem cells. We looked at swimming velocity, an important factor contributing to fertilization success, and flagellar length. Variability in swimming velocity was significantly higher in sperm from mutagenized males than in control sperm, reflecting their increased allelic diversity. In contrast, flagellar length, which is fixed during diploid stages of spermatogenesis, was unaffected by ENU treatment. We briefly discuss the implications of intra-ejaculate screening for maintenance of anisogamy and for outcomes of assisted reproductive technology.

Evolution and dysfunction of human cognitive and social traits: A transcriptional regulation perspective

Evolutionary changes in brain and craniofacial development have endowed humans with unique cognitive and social skills, but also predisposed us to debilitating disorders in which these traits are disrupted. What are the developmental genetic underpinnings that connect the adaptive evolution of our cognition and sociality with the persistence of mental disorders with severe negative fitness effects? We argue that loss of function of genes involved in transcriptional regulation represents a crucial link between the evolution and dysfunction of human cognitive and social traits. The argument is based on the haploinsufficiency of many transcriptional regulator genes, which makes them particularly sensitive to loss-of-function mutations. We discuss how human brain and craniofacial traits evolved through partial loss of function (i.e. reduced expression) of these genes, a perspective compatible with the idea of human self-domestication. Moreover, we explain why selection against loss-of-function variants supports the view that mutation-selection-drift, rather than balancing selection, underlies the persistence of psychiatric disorders. Finally, we discuss testable predictions.

Human HspB1, HspB3, HspB5 and HspB8: Shaping these disease factors during vertebrate evolution

Small heat shock proteins (sHSPs) emerged early in evolution and occur in all domains of life and nearly in all species, including humans. Mutations in four sHSPs (HspB1, HspB3, HspB5, HspB8) are associated with neuromuscular disorders. The aim of this study is to investigate the evolutionary forces shaping these sHSPs during vertebrate evolution. We performed comparative evolutionary analyses on a set of orthologous sHSP sequences, based on the ratio of non-synonymous: synonymous substitution rates for each codon. We found that these sHSPs had been historically exposed to different degrees of purifying selection, decreasing in this order: HspB8 > HspB1, HspB5 > HspB3. Within each sHSP, regions with different degrees of purifying selection can be discerned, resulting in characteristic selective pressure profiles. The conserved α-crystallin domains were exposed to the most stringent purifying selection compared to the flanking regions, supporting a 'dimorphic pattern' of evolution. Thus, during vertebrate evolution the different sequence partitions were exposed to different and measurable degrees of selective pressures. Among the disease-associated mutations, most are missense mutations primarily in HspB1 and to a lesser extent in the other sHSPs. Our data provide an explanation for this disparate incidence. Contrary to the expectation, most missense mutations cause dominant disease phenotypes. Theoretical considerations support a connection between the historic exposure of these sHSP genes to a high degree of purifying selection and the unusual prevalence of genetic dominance of the associated disease phenotypes. Our study puts the genetics of inheritable sHSP-borne diseases into the context of vertebrate evolution.

Missplicing suppressor alleles of Arabidopsis PRE-MRNA PROCESSING FACTOR 8 increase splicing fidelity by reducing the use of novel splice sites

Efficient splicing requires a balance between high-fidelity splice-site (SS) selection and speed. In Saccharomyces cerevisiae, Pre-mRNA processing factor 8 (Prp8) helps to balance precise SS selection and rapid, efficient intron excision and exon joining. argonaute1-52 (ago1-52) and incurvata13 (icu13) are hypomorphic alleles of the Arabidopsis thaliana genes ARGONAUTE1 (AGO1) and AUXIN RESISTANT6 (AXR6) that harbor point mutations creating a novel 3'SS and 5'SS, respectively. The spliceosome recognizes these novel SSs, as well as the intact genuine SSs, producing a mixture of wild-type and aberrant mature mRNAs. Here, we characterized five novel mutant alleles of PRP8 (one of the two Arabidopsis co-orthologs of yeast Prp8), naming these alleles morphology of ago1-52 suppressed5 (mas5). In the mas5-1 background, the spliceosome preferentially recognizes the intact genuine 3'SS of ago1-52 and 5'SS of icu13. Since point mutations that damage genuine SSs make the spliceosome prone to recognizing cryptic SSs, we also tested alleles of four genes carrying damaged genuine SSs, finding that mas5-1 did not suppress their missplicing. The mas5-1 and mas5-3 mutations represent a novel class of missplicing suppressors that increase splicing fidelity by hampering the use of novel SSs, but do not alter general pre-mRNA splicing.

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.

Novel Compound Heterozygous Variants in the CDC6 Gene in a Russian Patient with Meier-Gorlin Syndrome

Background

Meier-Gorlin syndrome (MGS) is a rare genetic syndrome inherited in an autosomal dominant or autosomal recessive manner. The disorder is characterized by bilateral microtia, absence or hypoplasia of the patella, and an intrauterine growth retardation as well as a number of other characteristic features. The cause of the disease is mutations in genes encoding proteins involved in the regulation of the cell cycle (ORC1, ORC4, ORC6, CDT1, CDC6, GMNN, CDC45L, MCM3, MCM5, MCM7, GINS2, and DONSON). Meier-Gorlin syndrome 5 due to mutations in the CDC6 gene is difficult to diagnose, and few clinical data have been described to date. Only one patient (male) with a missense mutation in a homozygous state has been previously reported. This report describes a new clinical case of Meier-Gorlin syndrome 5. This is also the first report of a Russian patient with Meier-Gorlin syndrome.

Case Presentation

The patient, a female, had extremely low physical development, neonatal progeroid appearance, lipodystrophy, thin skin, partial alopecia, cyanosis of the face, triangular face, microgenia, arachnodactyly, delayed bone age, hepatomegaly, hypoplasia of the labia majora, and hypertrophy of the clitoris in addition to known clinical signs. Differential diagnosis was performed with chromosomal abnormalities and Hutchinson-Gilford progeria. According to the results of sequencing of the clinical exome, the patient had two previously undescribed variants in the CDC6 gene, c.230A>G (p.(Lys77Arg)) and c.232C>T (p.(Gln78Ter)), NM_001254.3, in a compound heterozygous state.

Conclusion

This case allows us to learn more about the clinical features and nature of MGS 5 and improve the speed of diagnostics and quality of genetic counseling for such families.

Identification of discriminative gene-level and protein-level features associated with pathogenic gain-of-function and loss-of-function variants

Identifying whether a given genetic mutation results in a gene product with increased (gain-of-function; GOF) or diminished (loss-of-function; LOF) activity is an important step toward understanding disease mechanisms because they may result in markedly different clinical phenotypes. Here, we generated an extensive database of documented germline GOF and LOF pathogenic variants by employing natural language processing (NLP) on the available abstracts in the Human Gene Mutation Database. We then investigated various gene- and protein-level features of GOF and LOF variants and applied machine learning and statistical analyses to identify discriminative features. We found that GOF variants were enriched in essential genes, for autosomal-dominant inheritance, and in protein binding and interaction domains, whereas LOF variants were enriched in singleton genes, for protein-truncating variants, and in protein core regions. We developed a user-friendly web-based interface that enables the extraction of selected subsets from the GOF/LOF database by a broad set of annotated features and downloading of up-to-date versions. These results improve our understanding of how variants affect gene/protein function and may ultimately guide future treatment options.

Multiple mutations in the predicted transmembrane domains of the cellulose synthase 3 (CesA3) of Phytophthora capsici can confer semi-dominant resistance to carboxylic acid amide fungicides

The current study established a clearer understanding of the molecular basis for resistance to carboxylic acid amide (CAA) fungicides. Although four cellulose synthase (CesA) genes were investigated, only F1073L, G1105A, V1109L in CesA3 were found to link to CAA-resistance in Phytophthora capsici. Back-transformation experiments confirmed the role of the three mutations in CAA-resistance. Inheritance studies also confirmed the link and indicated the resistance was semi-dominant with the heterozygous F₁ and F₂ progeny exhibiting intermediate resistance levels compared to the homozygous parents, which was validated by the pyrosequencing results. The semi-dominant nature of CAA-resistance implies that it could be easy for resistance to spread once resistance emerged, being facilitated by both sexual and asexual reproduction. Bioinformatic analysis indicated all mutations occurred in either the first or second of the predicted transmembrane domains at C-terminus of CesA3. Resistant isolates bearing different combinations of mutations were found to exhibit different resistance levels to different CAAs, indicating that each mutation could make different contributions to resistance phenotype depending on structural differences in different CAAs. The current results highlight the complex combinations of mutations and resistance phenotype, and further reinforces the research necessity to completely characterize CAA-resistance to develop appropriate strategies to manage resistance development.

The tumor suppression theory of aging

Despite continued increases in human life expectancy, the factors determining the rate of human biological aging remain unknown. Without understanding the molecular mechanisms underlying aging, efforts to prevent aging are unlikely to succeed. The tumor suppression theory of aging introduced here proposes somatic mutation as the proximal cause of aging, but postulates that oncogenic transformation and clonal expansion, not functional impairment, are the relevant consequences of somatic mutation. Obesity and caloric restriction accelerate and decelerate aging due to their effect on cell proliferation, during which most mutations arise. Most phenotypes of aging are merely tumor-suppressive mechanisms that evolved to limit malignant growth, the dominant age-related cause of death in early and middle life. Cancer limits life span for most long-lived mammals, a phenomenon known as Peto's paradox. Its conservation across species demonstrates that mutation is a fundamental but hard limit on mammalian longevity. Cell senescence and apoptosis and differentiation induced by oncogenes, telomere shortening or DNA damage evolved as a second line of defense to limit the tumorigenic potential of clonally expanding cells, but accumulating senescent cells, senescence-associated secretory phenotypes and stem cell exhaustion eventually cause tissue dysfunction and the majority, if not most, phenotypes of aging.

A fresh look to the phenotype in mono-allelic likely pathogenic variants of the leptin and the leptin receptor gene

Leptin (LEP) and leptin receptor (LEPR) play a major role in energy homeostasis, metabolism, and reproductive function. While effects of biallelic likely pathogenic variants (-/-) on the phenotype are well characterized, effects of mono-allelic likely pathogenic variants (wt/-) in the LEP and LEPR gene on the phenotype compared to wild-type homozygosity (wt/wt) have not been systematically investigated. We identified in our systematic review 44 animal studies (15 on Lep, 29 on Lepr) and 39 studies in humans reporting on 130 mono-allelic likely pathogenic variant carriers with 20 distinct LEP variants and 108 heterozygous mono-allelic likely pathogenic variant carriers with 35 distinct LEPR variants. We found indications for a higher weight status in carriers of mono-allelic likely pathogenic variant in the leptin and in the leptin receptor gene compared to wt/wt, in both animal and human studies. In addition, animal studies showed higher body fat percentage in Lep and Lepr wt/- vs wt/wt. Animal studies provided indications for lower leptin levels in Lep wt/- vs. wt/wt and indications for higher leptin levels in Lepr wt/- vs wt/wt. Data on leptin levels in human studies was limited. Evidence for an impaired metabolism in mono-allelic likely pathogenic variants of the leptin and in leptin receptor gene was not conclusive (animal and human studies). Mono-allelic likely pathogenic variants in the leptin and in leptin receptor gene have phenotypic effects disposing to increased body weight and fat accumulation.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

N471D WASH complex subunit strumpellin knock-in mice display mild motor and cardiac abnormalities and BPTF and KLHL11 dysregulation in brain tissue

AIMS

We investigated N471D WASH complex subunit strumpellin (Washc5) knock-in and Washc5 knock-out mice as models for hereditary spastic paraplegia type 8 (SPG8).

METHODS

We generated heterozygous and homozygous N471D Washc5 knock-in mice and subjected them to a comprehensive clinical, morphological and laboratory parameter screen, and gait analyses. Brain tissue was used for proteomic analysis. Furthermore, we generated heterozygous Washc5 knock-out mice. WASH complex subunit strumpellin expression was determined by qPCR and immunoblotting.

RESULTS

Homozygous N471D Washc5 knock-in mice showed mild dilated cardiomyopathy, decreased acoustic startle reactivity, thinner eye lenses, increased alkaline phosphatase and potassium levels and increased white blood cell counts. Gait analyses revealed multiple aberrations indicative of locomotor instability. Similarly, the clinical chemistry, haematology and gait parameters of heterozygous mice also deviated from the values expected for healthy animals, albeit to a lesser extent. Proteomic analysis of brain tissue depicted consistent upregulation of BPTF and downregulation of KLHL11 in heterozygous and homozygous knock-in mice. WASHC5-related protein interaction partners and complexes showed no change in abundancies. Heterozygous Washc5 knock-out mice showing normal WASHC5 levels could not be bred to homozygosity.

CONCLUSIONS

While biallelic ablation of Washc5 was prenatally lethal, expression of N471D mutated WASHC5 led to several mild clinical and laboratory parameter abnormalities, but not to a typical SPG8 phenotype. The consistent upregulation of BPTF and downregulation of KLHL11 suggest mechanistic links between the expression of N471D mutated WASHC5 and the roles of both proteins in neurodegeneration and protein quality control, respectively.

Further evidence for functional recovery of AQP2 mutations associated with nephrogenic diabetes insipidus

Aquaporin-2 (AQP2) is a homotetrameric water channel responsible for the final water reuptake in the kidney. Disease-causing AQP2 mutations induce nephrogenic diabetes insipidus (NDI), a condition that challenges the bodily water balance by producing large urinary volumes. In this study, we characterize three new AQP2 mutations identified in our lab from NDI patients (A120D, A130V, T179N) along the previously reported A47V variant. Using Xenopus oocytes, we compared the key functional and biochemical features of these mutations against classical recessive (R187C) and dominant (R254Q) forms, and once again found clear functional recovery features (increased protein stability and function) for all mutations under study. This behaviour, attributed to heteromerization to wt-AQP2, challenge the classical model to NDI which often depicts recessive mutations as ill-structured proteins unable to oligomerize. Consequently, we propose a revised model to the cell pathophysiology of AQP2-related NDI which accounts for the functional recovery of recessive AQP2 mutations.

The mitochondrial ADP/ATP carrier exists and functions as a monomer

For more than 40 years, the oligomeric state of members of the mitochondrial carrier family (SLC25) has been the subject of debate. Initially, the consensus was that they were dimeric, based on the application of a large number of different techniques. However, the structures of the mitochondrial ADP/ATP carrier, a member of the family, clearly demonstrated that its structural fold is monomeric, lacking a conserved dimerisation interface. A re-evaluation of previously published data, with the advantage of hindsight, concluded that technical errors were at the basis of the earlier dimer claims. Here, we revisit this topic, as new claims for the existence of dimers of the bovine ADP/ATP carrier have emerged using native mass spectrometry of mitochondrial membrane vesicles. However, the measured mass does not agree with previously published values, and a large number of post-translational modifications are proposed to account for the difference. Contrarily, these modifications are not observed in electron density maps of the bovine carrier. If they were present, they would interfere with the structure and function of the carrier, including inhibitor and substrate binding. Furthermore, the reported mass does not account for three tightly bound cardiolipin molecules, which are consistently observed in other studies and are important stabilising factors for the transport mechanism. The monomeric carrier has all of the required properties for a functional transporter and undergoes large conformational changes that are incompatible with a stable dimerisation interface. Thus, our view that the native mitochondrial ADP/ATP carrier exists and functions as a monomer remains unaltered.

Novel REST Truncation Mutations Causing Hereditary Gingival Fibromatosis

Hereditary gingival fibromatosis (HGF) is a rare genetic disorder featured by nonsyndromic pathological overgrowth of gingiva. The excessive gingival tissues can cause dental, masticatory, and phonetic problems, which impose severe functional and esthetic burdens on affected individuals. Due to its high recurrent rate, patients with HGF have to undergo repeated surgical procedures of gingival resection, from childhood to adulthood, which significantly compromises their quality of life. Unraveling the genetic etiology and molecular pathogenesis of HGF not only gains insight into gingival physiology and homeostasis but also opens avenues for developing potential therapeutic strategies for this disorder. Recently, mutations in REST (OMIM *600571), encoding a transcription repressor, were reported to cause HGF (GINGF5; OMIM #617626) in 3 Turkish families. However, the functions of REST in gingival homeostasis and pathogenesis of REST-associated HGF remain largely unknown. In this study, we characterized 2 HGF families and identified 2 novel REST mutations, c.2449C>T (p.Arg817*) and c.2771_2793dup (p.Glu932Lysfs*3). All 5 mutations reported to date are nonsenses or frameshifts in the last exon of REST and would presumably truncate the protein. In vitro reporter gene assays demonstrated a partial or complete loss of repressor activity for these truncated RESTs. When coexpressed with the full-length protein, the truncated RESTs impaired the repressive ability of wild-type REST, suggesting a dominant negative effect. Immunofluorescent studies showed nuclear localization of overexpressed wild-type and truncated RESTs in vitro, indicating preservation of the nuclear localization signal in shortened proteins. Immunohistochemistry demonstrated a comparable pattern of ubiquitous REST expression in both epithelium and lamina propria of normal and HGF gingival tissues despite a reduced reactivity in HGF gingiva. Results of this study confirm the pathogenicity of REST truncation mutations occurring in the last exon causing HGF and suggest the pathosis is caused by an antimorphic (dominant negative) disease mechanism.

MVP predicts the pathogenicity of missense variants by deep learning

Accurate pathogenicity prediction of missense variants is critically important in genetic studies and clinical diagnosis. Previously published prediction methods have facilitated the interpretation of missense variants but have limited performance. Here, we describe MVP (Missense Variant Pathogenicity prediction), a new prediction method that uses deep residual network to leverage large training data sets and many correlated predictors. We train the model separately in genes that are intolerant of loss of function variants and the ones that are tolerant in order to take account of potentially different genetic effect size and mode of action. We compile cancer mutation hotspots and de novo variants from developmental disorders for benchmarking. Overall, MVP achieves better performance in prioritizing pathogenic missense variants than previous methods, especially in genes tolerant of loss of function variants. Finally, using MVP, we estimate that de novo coding variants contribute to 7.8% of isolated congenital heart disease, nearly doubling previous estimates.

Homozygosity for the pathogenic RET hotspot variant p.Cys634Trp: A consanguineous family with MEN2A

Multiple endocrine neoplasia type 2 (MEN2) is a dominantly inherited condition with defined correlations between the genetic variant and clinical presentations. The location of pathogenic variants in the RET gene is a significant determinant of disease presentation and is associated with variable gene activation. Heterozygous pathogenic variants in codon 634 result in earlier onset of medullary thyroid carcinoma and higher incidence of phaeochromocytoma. Here we describe a consanguineous family with MEN2A that includes two children homozygous for the established pathogenic variant p. Cys634Trp. Both parents and a sibling were confirmed to being heterozygotes. Previous reports of biallelic or multiple RET variants have been limited to weakly activating variants. We present the first report of individuals homozygous for the highly activating RET p. Cys634Trp pathogenic variant and discuss disease severity and onset in this rare occurrence.

Mutant VAPB: Culprit or Innocent Bystander of Amyotrophic Lateral Sclerosis?

Nearly twenty years ago a mutation in the VAPB gene, resulting in a proline to serine substitution (p.P56S), was identified as the cause of a rare, slowly progressing, familial form of the motor neuron degenerative disease Amyotrophic Lateral Sclerosis (ALS). Since then, progress in unravelling the mechanistic basis of this mutation has proceeded in parallel with research on the VAP proteins and on their role in establishing membrane contact sites between the ER and other organelles. Analysis of the literature on cellular and animal models reviewed here supports the conclusion that P56S-VAPB, which is aggregation-prone, non-functional and unstable, is expressed at levels that are insufficient to support toxic gain-of-function or dominant negative effects within motor neurons. Instead, insufficient levels of the product of the single wild-type allele appear to be required for pathological effects, and may be the main driver of the disease. In light of the multiple interactions of the VAP proteins, we address the consequences of specific VAPB depletion and highlight various affected processes that could contribute to motor neuron degeneration. In the future, distinction of specific roles of each of the two VAP paralogues should help to further elucidate the basis of p.P56S familial ALS, as well as of other more common forms of the disease.

Predicting functional effects of missense variants in voltage-gated sodium and calcium channels

Malfunctions of voltage-gated sodium and calcium channels (encoded by SCNxA and CACNA1x family genes, respectively) have been associated with severe neurologic, psychiatric, cardiac, and other diseases. Altered channel activity is frequently grouped into gain or loss of ion channel function (GOF or LOF, respectively) that often corresponds not only to clinical disease manifestations but also to differences in drug response. Experimental studies of channel function are therefore important, but laborious and usually focus only on a few variants at a time. On the basis of known gene-disease mechanisms of 19 different diseases, we inferred LOF (n = 518) and GOF (n = 309) likely pathogenic variants from the disease phenotypes of variant carriers. By training a machine learning model on sequence- and structure-based features, we predicted LOF or GOF effects [area under the receiver operating characteristics curve (ROC) = 0.85] of likely pathogenic missense variants. Our LOF versus GOF prediction corresponded to molecular LOF versus GOF effects for 87 functionally tested variants in SCN1/2/8A and CACNA1I (ROC = 0.73) and was validated in exome-wide data from 21,703 cases and 128,957 controls. We showed respective regional clustering of inferred LOF and GOF nucleotide variants across the alignment of the entire gene family, suggesting shared pathomechanisms in the SCNxA/CACNA1x family genes.

Phylogenetic Modeling of Regulatory Element Turnover Based on Epigenomic Data

Evolutionary changes in gene expression are often driven by gains and losses of cis-regulatory elements (CREs). The dynamics of CRE evolution can be examined using multispecies epigenomic data, but so far such analyses have generally been descriptive and model-free. Here, we introduce a probabilistic modeling framework for the evolution of CREs that operates directly on raw chromatin immunoprecipitation and sequencing (ChIP-seq) data and fully considers the phylogenetic relationships among species. Our framework includes a phylogenetic hidden Markov model, called epiPhyloHMM, for identifying the locations of multiply aligned CREs, and a combined phylogenetic and generalized linear model, called phyloGLM, for accounting for the influence of a rich set of genomic features in describing their evolutionary dynamics. We apply these methods to previously published ChIP-seq data for the H3K4me3 and H3K27ac histone modifications in liver tissue from nine mammals. We find that enhancers are gained and lost during mammalian evolution at about twice the rate of promoters, and that turnover rates are negatively correlated with DNA sequence conservation, expression level, and tissue breadth, and positively correlated with distance from the transcription start site, consistent with previous findings. In addition, we find that the predicted dosage sensitivity of target genes positively correlates with DNA sequence constraint in CREs but not with turnover rates, perhaps owing to differences in the effect sizes of the relevant mutations. Altogether, our probabilistic modeling framework enables a variety of powerful new analyses.

Loss-of-function variants in myocardin cause congenital megabladder in humans and mice

Myocardin (MYOCD) is the founding member of a class of transcriptional coactivators that bind the serum-response factor to activate gene expression programs critical in smooth muscle (SM) and cardiac muscle development. Insights into the molecular functions of MYOCD have been obtained from cell culture studies, and to date, knowledge about in vivo roles of MYOCD comes exclusively from experimental animals. Here, we defined an often lethal congenital human disease associated with inheritance of pathogenic MYOCD variants. This disease manifested as a massively dilated urinary bladder, or megabladder, with disrupted SM in its wall. We provided evidence that monoallelic loss-of-function variants in MYOCD caused congenital megabladder in males only, whereas biallelic variants were associated with disease in both sexes, with a phenotype additionally involving the cardiovascular system. These results were supported by cosegregation of MYOCD variants with the phenotype in 4 unrelated families by in vitro transactivation studies in which pathogenic variants resulted in abrogated SM gene expression and by the finding of megabladder in 2 distinct mouse models with reduced Myocd activity. In conclusion, we have demonstrated that variants in MYOCD result in human disease, and the collective findings highlight a vital role for MYOCD in mammalian organogenesis.

Sequencing of lipoprotein lipase gene in the Mediterranean river buffalo identified novel variants affecting gene expression

Lipoprotein lipase (LPL) is a key enzyme for lipid metabolism, playing a fundamental role in the composition of fat in adipose tissue and milk. The LPL gene has been seldom investigated in dairy ruminants and barely studied in river buffalo (Bubalus bubalis). The aim of this work was to explore the genetic diversity of LPL and its promoter and to identify functional mutations, using a combined approach based on sequencing, dual-color electrophoretic mobility shift assay, and quantitative PCR. Thirteen consensus sequences for transcription factors were found in the promoter. Eleven SNP were detected, and the attention was focused on the SNP with potential functional effects: g.-446A>G, because the presence of G created a consensus motif for the transcription factor Sp1, and g.107A>G, which was the only exonic SNP. We developed PCR-RFLP methods for genotyping the 2 SNP and calculated the allele frequencies. A strong linkage disequilibrium (D' = 1; r² = 0.903) was found between the 2 SNP. The dual-color electrophoretic mobility shift assay demonstrated that only genotype g.-446GG allowed the binding of the Sp1 transcription factor, resulting in overexpression of the gene (~2.5 fold), as confirmed by the quantitative PCR results. Haploinsufficiency is proposed as a regulation mechanism. This study adds further knowledge on the structure of the LPL gene and its expression in river buffalo, with potential effects on milk qualitative and quantitative production.

Identification of two novel mutations in POU4F3 gene associated with autosomal dominant hearing loss in Chinese families

Autosomal dominant non‐syndromic hearing loss is genetically heterogeneous with 47 genes identified to date, including POU4F3. In this study, by using a next‐generation sequencing panel targeting 127 deafness genes, we identified a pathogenic frameshift mutation c.704_705del and a missense mutation c.593G>A in two three‐generation Chinese families with late‐onset progressive ADNSHL, respectively. The novel mutations of POU4F3 co‐segregated with the deafness phenotype in these two families. c.704_705del caused a frameshift p.T235fs and c.593G>A caused an amino acid substitution of p.R198H. Both mutations led to an abnormal and incomplete protein structure. POU4F3 with either of the two mutations was transiently transfected into HEI‐OC1 and HEK 293 cell lines and immunofluorescence assay was performed to investigate the subcellular localization of mutated protein. The results indicated that both c.704_705del (p.T235fs) and c.593G>A (p.R198H) could impair the nuclear localization function of POU4F3. The p.R198H POU4F3 protein was detected as a weak band of the correct molecular weight, indicating that the stability of p.R198H POU4F3 differed from that of the wild‐type protein. While, the p.T235fs POU4F3 protein was expressed with a smaller molecular weight, implying this mutation result in a frameshift and premature termination of the POU4F3 protein. In summary, we report two novel mutations of POU4F3 associated with progressive ADNSHL and explored their effects on POU4F3 nuclear localization. These findings expanded the mutation spectrum of POU4F3 and provided new knowledge for the pathogenesis of POU4F3 in hearing loss.

Genomic analysis of a riboflavin-overproducing Ashbya gossypii mutant isolated by disparity mutagenesis

BACKGROUND

Ashbya gossypii naturally overproduces riboflavin and has been utilized for industrial riboflavin production. To improve riboflavin production, various approaches have been developed. In this study, to investigate the change in metabolism of a riboflavin-overproducing mutant, namely, the W122032 strain (MT strain) that was isolated by disparity mutagenesis, genomic analysis was carried out.

RESULTS

In the genomic analysis, 33 homozygous and 1377 heterozygous mutations in the coding sequences of the genome of MT strain were detected. Among these heterozygous mutations, the proportion of mutated reads in each gene was different, ranging from 21 to 75%. These results suggest that the MT strain may contain multiple nuclei containing different mutations. We tried to isolate haploid spores from the MT strain to prove its ploidy, but this strain did not sporulate under the conditions tested. Heterozygous mutations detected in genes which are important for sporulation likely contribute to the sporulation deficiency of the MT strain. Homozygous and heterozygous mutations were found in genes encoding enzymes involved in amino acid metabolism, the TCA cycle, purine and pyrimidine nucleotide metabolism and the DNA mismatch repair system. One homozygous mutation in AgILV2 gene encoding acetohydroxyacid synthase, which is also a flavoprotein in mitochondria, was found. Gene ontology (GO) enrichment analysis showed heterozygous mutations in all 22 DNA helicase genes and genes involved in oxidation-reduction process.

CONCLUSION

This study suggests that oxidative stress and the aging of cells were involved in the riboflavin over-production in A. gossypii riboflavin over-producing mutant and provides new insights into riboflavin production in A. gossypii and the usefulness of disparity mutagenesis for the creation of new types of mutants for metabolic engineering.

Enhancer architecture sensitizes cell specific responses to Notch gene dose via a bind and discard mechanism

Notch pathway haploinsufficiency can cause severe developmental syndromes with highly variable penetrance. Currently, we have a limited mechanistic understanding of phenotype variability due to gene dosage. Here, we unexpectedly found that inserting an enhancer containing pioneer transcription factor sites coupled to Notch dimer sites can induce a subset of Notch haploinsufficiency phenotypes in Drosophila with wild type Notch gene dose. Using Drosophila genetics, we show that this enhancer induces Notch phenotypes in a Cdk8-dependent, transcription-independent manner. We further combined mathematical modeling with quantitative trait and expression analysis to build a model that describes how changes in Notch signal production versus degradation differentially impact cellular outcomes that require long versus short signal duration. Altogether, these findings support a 'bind and discard' mechanism in which enhancers with specific binding sites promote rapid Cdk8-dependent Notch turnover, and thereby reduce Notch-dependent transcription at other loci and sensitize tissues to gene dose based upon signal duration.

A framework for high-resolution phenotyping of candidate male infertility mutants: from human to mouse

Male infertility is a heterogeneous condition of largely unknown etiology that affects at least 7% of men worldwide. Classical genetic approaches and emerging next-generation sequencing studies support genetic variants as a frequent cause of male infertility. Meanwhile, the barriers to transmission of this disease mean that most individual genetic cases will be rare, but because of the large percentage of the genome required for spermatogenesis, the number of distinct causal mutations is potentially large. Identifying bona fide causes of male infertility thus requires advanced filtering techniques to select for high-probability candidates, including the ability to test causality in animal models. The mouse remains the gold standard for defining the genotype-phenotype connection in male fertility. Here, we present a best practice guide consisting of (a) major points to consider when interpreting next-generation sequencing data performed on infertile men, and, (b) a systematic strategy to categorize infertility types and how they relate to human male infertility. Phenotyping infertility in mice can involve investigating the function of multiple cell types across the testis and epididymis, as well as sperm function. These findings will feed into the diagnosis and treatment of male infertility as well as male health broadly.

Identification of six novel variants in Waardenburg syndrome type II by next-generation sequencing

Background

Waardenburg syndrome (WS) is a dominantly inherited, genetically heterogeneous auditory‐pigmentary syndrome characterized by nonprogressive sensorineural hearing loss and iris discoloration. This study aimed to investigate the underlying molecular pathology in Chinese WS families.

Methods

A total of 13 patients with Waardenburg syndrome type II (WS2) from six unrelated Chinese families were enrolled. We investigated the mutation profile of genes related to congenital deafness in these families through a targeted sequencing technology and validated the candidate variants by Sanger sequencing.

Results

We identified six novel variants in microphthalmia‐associated transcription factor (MITF) and SRY‐box 10 (SOX10), which were predicted to be disease causing by in silico analysis. Our results showed that mutations in SOX10 and MITF are two major causes of deafness associated with WS, and de novo mutations were frequently found in probands with SOX10 mutations but not in those with MITF mutations.

Conclusion

Results showed that targeted next‐generation sequencing (NGS) enabled us to detect disease‐causing mutations with high accuracy, stability, speed and throughput. Our study extends the pathogenic mutation spectrum of MITF and SOX10.

A Drosophila model of neuronal ceroid lipofuscinosis CLN4 reveals a hypermorphic gain of function mechanism

The autosomal dominant neuronal ceroid lipofuscinoses (NCL) CLN4 is caused by mutations in the synaptic vesicle (SV) protein CSPα. We developed animal models of CLN4 by expressing CLN4 mutant human CSPα (hCSPα) in Drosophila neurons. Similar to patients, CLN4 mutations induced excessive oligomerization of hCSPα and premature lethality in a dose-dependent manner. Instead of being localized to SVs, most CLN4 mutant hCSPα accumulated abnormally, and co-localized with ubiquitinated proteins and the prelysosomal markers HRS and LAMP1. Ultrastructural examination revealed frequent abnormal membrane structures in axons and neuronal somata. The lethality, oligomerization and prelysosomal accumulation induced by CLN4 mutations was attenuated by reducing endogenous wild type (WT) dCSP levels and enhanced by increasing WT levels. Furthermore, reducing the gene dosage of Hsc70 also attenuated CLN4 phenotypes. Taken together, we suggest that CLN4 alleles resemble dominant hypermorphic gain of function mutations that drive excessive oligomerization and impair membrane trafficking.

Importance of the Conserved Carboxyl-Terminal CNOT1 Binding Domain to Tristetraprolin Activity In Vivo

Tristetraprolin (TTP) is an anti-inflammatory protein that modulates the stability of certain cytokine/chemokine mRNAs. After initial high-affinity binding to AU-rich elements in 3' untranslated regions of target mRNAs, mediated through its tandem zinc finger (TZF) domain, TTP promotes the deadenylation and ultimate decay of target transcripts. These transcripts and their encoded proteins accumulate abnormally in TTP knockout (KO) mice, leading to a severe inflammatory syndrome. To assess the importance of the highly conserved C-terminal CNOT1 binding domain (CNBD) of TTP to the TTP deficiency phenotype in mice, we created a mouse model in which TTP lacked its CNBD. CNBD deletion mice exhibited a less severe phenotype than the complete TTP KO mice. In macrophages, the stabilization of target transcripts seen in KO mice was partially normalized in the CNBD deletion mice. In cell-free experiments, recombinant TTP lacking its CNBD could still activate target mRNA deadenylation by purified recombinant Schizosaccharomyces pombe CCR4/NOT complexes, although to a lesser extent than full-length TTP. Thus, TTP lacking its CNBD can still act to promote target mRNA instability in vitro and in vivo These data have implications for TTP family members throughout the eukarya, since species from all four kingdoms contain proteins with linked TZF and CNOT1 binding domains.

A Model of Hereditary Sensory and Autonomic Neuropathy Type 1 Reveals a Role of Glycosphingolipids in Neuronal Polarity

Hereditary sensory and autonomic neuropathy Type 1 (HSAN1) is a rare autosomal dominantly inherited neuropathy, clinically characterized by a loss of distal peripheral sensory and motoneuronal function. Mutations in subunits of serine palmitoyltransferase (SPT) have been linked to the majority of HSAN1 cases. SPTs catalyze the condensation of l-serine with palmitoyl-CoA, the first committed and rate-limiting step in de novo sphingolipid biosynthesis. Despite extensive investigation, the molecular pathogenesis of HSAN1 remains controversial. Here, we established a Caenorhabditis elegans (C. elegans) model of HSAN1 by generating a sptl-1(c363g) mutation, encoding SPTL-1(C121W) and equivalent to human SPTLC1^C133W, at the C. elegans genomic locus through CRISPR. The sptl-1(c363g) homozygous mutants exhibited the same larval lethality and epithelial polarity defect as observed in sptl-1(RNAi) animals, suggesting a loss-of-function effect of the SPTL-1(C121W) mutation. sptl-1(c363g)/+ heterozygous mutants displayed sensory dysfunction with concomitant neuronal morphology and axon-dendrite polarity defects, demonstrating that the C. elegans model recapitulates characteristics of the human disease. sptl-1(c363g)-derived neuronal defects were copied in animals with defective sphingolipid biosynthetic enzymes downstream of SPTL-1, including ceramide glucosyltransferases, suggesting that SPTLC1^C133W contributes to the HSAN1 pathogenesis by limiting the production of complex sphingolipids, including glucosylceramide. Overexpression of SPTL-1(C121W) led to similar epithelial and neuronal defects and to reduced levels of complex sphingolipids, specifically glucosylceramide, consistent with a dominant-negative effect of SPTL-1(C121W) that is mediated by loss of this downstream product. Genetic interactions between SPTL-1(C121W) and components of directional trafficking in neurons suggest that the neuronal polarity phenotype could be caused by glycosphingolipid-dependent defects in polarized vesicular trafficking.SIGNIFICANCE STATEMENT The symptoms of inherited metabolic diseases are often attributed to the accumulation of toxic intermediates or byproducts, no matter whether the disease-causing enzyme participates in a biosynthetic or a degradation pathway. By showing that the phenotypes observed in a C. elegans model of HSAN1 disease could be caused by loss of a downstream product (glucosylceramide) rather than the accumulation of a toxic byproduct, our work provides new insights into the origins of the symptoms of inherited metabolic diseases while expanding the repertoire of sphingolipid functions, specifically, of glucosylceramides. These findings not only have their most immediate relevance for neuroprotective treatments for HSAN1, they may also have implications for a much broader range of neurologic conditions.

De Novo Missense Substitutions in the Gene Encoding CDK8, a Regulator of the Mediator Complex, Cause a Syndromic Developmental Disorder

The Mediator is an evolutionarily conserved, multi-subunit complex that regulates multiple steps of transcription. Mediator activity is regulated by the reversible association of a four-subunit module comprising CDK8 or CDK19 kinases, together with cyclin C, MED12 or MED12L, and MED13 or MED13L. Mutations in MED12, MED13, and MED13L were previously identified in syndromic developmental disorders with overlapping phenotypes. Here, we report CDK8 mutations (located at 13q12.13) that cause a phenotypically related disorder. Using whole-exome or whole-genome sequencing, and by international collaboration, we identified eight different heterozygous missense CDK8 substitutions, including 10 shown to have arisen de novo, in 12 unrelated subjects; a recurrent mutation, c.185C>T (p.Ser62Leu), was present in five individuals. All predicted substitutions localize to the ATP-binding pocket of the kinase domain. Affected individuals have overlapping phenotypes characterized by hypotonia, mild to moderate intellectual disability, behavioral disorders, and variable facial dysmorphism. Congenital heart disease occurred in six subjects; additional features present in multiple individuals included agenesis of the corpus callosum, ano-rectal malformations, seizures, and hearing or visual impairments. To evaluate the functional impact of the mutations, we measured phosphorylation at STAT1-Ser727, a known CDK8 substrate, in a CDK8 and CDK19 CRISPR double-knockout cell line transfected with wild-type (WT) or mutant CDK8 constructs. These experiments demonstrated a reduction in STAT1 phosphorylation by all mutants, in most cases to a similar extent as in a kinase-dead control. We conclude that missense mutations in CDK8 cause a developmental disorder that has phenotypic similarity to syndromes associated with mutations in other subunits of the Mediator kinase module, indicating probable overlap in pathogenic mechanisms.

Chemical Biology Gateways to Mapping Location, Association, and Pathway Responsivity

Here we discuss, how by applying chemical concepts to biological problems, methods have been developed to map spatiotemporal regulation of proteins and small-molecule modulation of proteome signaling responses. We outline why chemical-biology platforms are ideal for such purposes. We further discuss strengths and weaknesses of chemical-biology protocols, contrasting them against classical genetic and biochemical approaches. We make these evaluations based on three parameters: occupancy; functional information; and spatial restriction. We demonstrate how the specific choice of chemical reagent and experimental set-up unite to resolve biological problems. Potential improvements/extensions as well as specific controls that in our opinion are often overlooked or employed incorrectly are also considered. Finally, we discuss some of the latest emerging methods to illuminate how chemical-biology innovations provide a gateway toward information hitherto inaccessible by conventional genetic/biochemical means. Finally, we also caution against solely relying on chemical-biology strategies and urge the field to undertake orthogonal validations to ensure robustness of results.