SMAD6 variants in craniosynostosis : genotype and phenotype evaluation

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism near BMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure

Identifying low density lipoprotein cholesterol associated variants in the Annexin A2 (ANXA2) gene.

BACKGROUND AND AIMS: Annexin-A2 (AnxA2) is an endogenous inhibitor of proprotein convertase subtilisin/kexin type-9 (PCSK9). The repeat-one (R1) domain of AnxA2 binds to PCSK9, blocking its ability to promote degradation of low-density lipoprotein cholesterol-receptors (LDL-R) and thereby regulate low-density lipoprotein cholesterol (LDL-C) levels. Here we identify variants in ANXA2 influencing LDL-C levels and we determine the molecular mechanisms of their effects. RESULTS: The ANXA2 single nucleotide polymorphism (SNP) genotype-phenotype association was examined using the Second-Northwick-Park Heart Study (NPHSII) (n∼2700) and the UCL-LSHTM-Edinburgh-Bristol (UCLEB) consortium (n∼14,600). The ANXA2-R1 domain coding-SNP rs17845226 (V98L) associated with LDL-C, homozygotes for the minor allele having ≈18.8% higher levels of LDL-C (p = 0.004), and higher risk of coronary heart disease (CHD) (p = 0.04). The SNP is in modest linkage disequilibrium (r(2) > 0.5) with two intergenic SNPs, rs17191344 and rs11633032. Both SNPs showed allele-specific protein binding, and the minor alleles caused significant reduction in reporter gene expression (≈18%, p < 0.001). In the expression quantitative trait loci (eQTL) study, minor allele homozygotes have significantly lower levels of ANXA2-mRNA expression (p = 1.36 × 10(-05)). CONCLUSIONS: Both rs11633032 and rs17191344 SNPs are functional variants, where the minor alleles create repressor-binding protein sites for transcription factors that contribute to reduced ANXA2 gene expression. Lower AnxA2 levels could increase plasma levels of PCSK9 and thus increase LDL-C levels and risk of CHD. This supports, for the first time in humans, previous observations in mouse models that changes in the levels of AnxA2 directly influence plasma LDL-C levels, and thus implicate this protein as a potential therapeutic target for LDL-C lowering

The inner junction protein CFAP20 functions in motile and non-motile cilia and is critical for vision

Motile and non-motile cilia are associated with mutually-exclusive genetic disorders. Motile cilia propel sperm or extracellular fluids, and their dysfunction causes primary ciliary dyskinesia. Non-motile cilia serve as sensory/signalling antennae on most cell types, and their disruption causes single-organ ciliopathies such as retinopathies or multi-system syndromes. CFAP20 is a ciliopathy candidate known to modulate motile cilia in unicellular eukaryotes. We demonstrate that in zebrafish, cfap20 is required for motile cilia function, and in C. elegans, CFAP-20 maintains the structural integrity of non-motile cilia inner junctions, influencing sensory-dependent signalling and development. Human patients and zebrafish with CFAP20 mutations both exhibit retinal dystrophy. Hence, CFAP20 functions within a structural/functional hub centered on the inner junction that is shared between motile and non-motile cilia, and is distinct from other ciliopathy-associated domains or macromolecular complexes. Our findings suggest an uncharacterised pathomechanism for retinal dystrophy, and potentially for motile and non-motile ciliopathies in general.</p

Assessing the digenic model in rare disorders using population sequencing data

An important fraction of patients with rare disorders remains with no clear genetic diagnostic, even after whole-exome or whole-genome sequencing, posing a difficulty in giving adequate treatment and genetic counseling. The analysis of genomic data in rare disorders mostly considers the presence of single gene variants in coding regions that follow a concrete monogenic mode of inheritance. A digenic inheritance, with variants in two functionally-related genes in the same individual, is a plausible alternative that might explain the genetic basis of the disease in some cases. In this case, digenic disease combinations should be absent or underrepresented in healthy individuals. We develop a framework to evaluate the significance of digenic combinations and test its statistical power in different scenarios. We suggest that this approach will be relevant with the advent of new sequencing efforts including hundreds of thousands of samples

Identification of 4 novel human ocular coloboma genes ANK3, BMPR1B, PDGFRA, and CDH4 through evolutionary conserved vertebrate gene analysis

Purpose: Ocular coloboma arises from genetic or environmental perturbations that inhibit optic fissure (OF) fusion during early eye development. Despite high genetic heterogeneity, 70% to 85% of patients remain molecularly undiagnosed. In this study, we have identified new potential causative genes using cross-species comparative meta-analysis. Methods: Evolutionarily conserved differentially expressed genes were identified through in silico analysis, with in situ hybridization, gene knockdown, and rescue performed to confirm spatiotemporal gene expression and phenotype. Interrogation of the 100,000 Genomes Project for putative pathogenic variants was performed. Results: Nine conserved differentially expressed genes between zebrafish and mouse were identified. Expression of zebrafish ank3a, bmpr1ba/b, cdh4, and pdgfaa was localized to the OF, periocular mesenchyme cells, or ciliary marginal zone, regions traversed by the OF. Knockdown of ank3, bmpr1b, and pdgfaa revealed a coloboma and/or microphthalmia phenotype. Novel pathogenic variants in ANK3, BMPR1B, PDGFRA, and CDH4 were identified in 8 unrelated coloboma families. We showed BMPR1B rescued the knockdown phenotype but variant messenger RNAs failed, providing evidence of pathogenicity. Conclusion: We show the utility of cross-species meta-analysis to identify several novel coloboma disease-causing genes. There is a potential to increase the diagnostic yield for new and unsolved patients while adding to our understanding of the genetic basis of OF morphogenesis

Author Correction: Nuclear-mitochondrial DNA segments resemble paternally inherited mitochondrial DNA in humans

An amendment to this paper has been published and can be accessed via a link at the top of the paper

Biallelic variants in KARS1 are associated with neurodevelopmental disorders and hearing loss recapitulated by the knockout zebrafish

Purpose: Pathogenic variants in Lysyl-tRNA synthetase 1 (KARS1) have increasingly been recognized as a cause of early-onset complex neurological phenotypes. To advance the timely diagnosis of KARS1-related disorders, we sought to delineate its phenotype and generate a disease model to understand its function in vivo. Methods: Through international collaboration, we identified 22 affected individuals from 16 unrelated families harboring biallelic likely pathogenic or pathogenic in KARS1 variants. Sequencing approaches ranged from disease-specific panels to genome sequencing. We generated loss-of-function alleles in zebrafish. Results: We identify ten new and four known biallelic missense variants in KARS1 presenting with a moderate-to-severe developmental delay, progressive neurological and neurosensory abnormalities, and variable white matter involvement. We describe novel KARS1-associated signs such as autism, hyperactive behavior, pontine hypoplasia, and cerebellar atrophy with prevalent vermian involvement. Loss of kars1 leads to upregulation of p53, tissue-specific apoptosis, and downregulation of neurodevelopmental related genes, recapitulating key tissue-specific disease phenotypes of patients. Inhibition of p53 rescued several defects of kars1−/− knockouts. Conclusion: Our work delineates the clinical spectrum associated with KARS1 defects and provides a novel animal model for KARS1-related human diseases revealing p53 signaling components as potential therapeutic targets

Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14

TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara−/−;tspearb−/− double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants

DYNC2H1 hypomorphic or retina-predominant variants cause nonsyndromic retinal degeneration

Purpose: Determining the role of DYNC2H1 variants in nonsyndromic inherited retinal disease (IRD). Methods: Genome and exome sequencing were performed for five unrelated cases of IRD with no identified variant. In vitro assays were developed to validate the variants identified (fibroblast assay, induced pluripotent stem cell [iPSC] derived retinal organoids, and a dynein motility assay). Results: Four novel DYNC2H1 variants (V1, g.103327020_103327021dup; V2, g.103055779A>T; V3, g.103112272C>G; V4, g.103070104A>C) and one previously reported variant (V5, g.103339363T>G) were identified. In proband 1 (V1/V2), V1 was predicted to introduce a premature termination codon (PTC), whereas V2 disrupted the exon 41 splice donor site causing incomplete skipping of exon 41. V1 and V2 impaired dynein-2 motility in vitro and perturbed IFT88 distribution within cilia. V3, homozygous in probands 2–4, is predicted to cause a PTC in a retina-predominant transcript. Analysis of retinal organoids showed that this new transcript expression increased with organoid differentiation. V4, a novel missense variant, was in trans with V5, previously associated with Jeune asphyxiating thoracic dystrophy (JATD). Conclusion: The DYNC2H1 variants discussed herein were either hypomorphic or affecting a retina-predominant transcript and caused nonsyndromic IRD. Dynein variants, specifically DYNC2H1 variants are reported as a cause of non syndromic IRD

Human and mouse essentiality screens as a resource for disease gene discovery.

The identification of causal variants in sequencing studies remains a considerable challenge that can be partially addressed by new gene-specific knowledge. Here, we integrate measures of how essential a gene is to supporting life, as inferred from viability and phenotyping screens performed on knockout mice by the International Mouse Phenotyping Consortium and essentiality screens carried out on human cell lines. We propose a cross-species gene classification across the Full Spectrum of Intolerance to Loss-of-function (FUSIL) and demonstrate that genes in five mutually exclusive FUSIL categories have differing biological properties. Most notably, Mendelian disease genes, particularly those associated with developmental disorders, are highly overrepresented among genes non-essential for cell survival but required for organism development. After screening developmental disorder cases from three independent disease sequencing consortia, we identify potentially pathogenic variants in genes not previously associated with rare diseases. We therefore propose FUSIL as an efficient approach for disease gene discovery