Neurodevelopmental disorder with microcephaly, hypotonia and variable brain anomalies (NMIHBA) is an autosomal recessive neurodevelopmental and neurodegenerative disorder characterized by global developmental delay and severe intellectual disability. Microcephaly, progressive cortical atrophy, cerebellar hypoplasia and delayed myelination are neurological hallmarks in affected individuals. NMIHBA is caused by biallelic variants in PRUNE1 encoding prune exopolyphosphatase 1. We provide in-depth clinical description of two affected siblings harboring compound heterozygous variant alleles, c.383G > A (p.Arg128Gln), c.520G > T (p.Gly174*) in PRUNE1. To gain insights into disease biology, we biochemically characterized missense variants within the conserved N-terminal aspartic acid-histidine-histidine (DHH) motif and provide evidence that they result in the destabilization of protein structure and/or loss of exopolyphosphatase activity. Genetic ablation of Prune1 results in midgestational lethality in mice, associated with perturbations to embryonic growth and vascular development. Our findings suggest that NMIHBA results from hypomorphic variant alleles in humans and underscore the potential key role of PRUNE1 exopolyphoshatase activity in neurodevelopment.

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma. Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness. To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness. The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry. We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10-8) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

TNNI3K expression worsens disease progression in several mouse heart pathology models. TNNI3K expression also reduces the number of diploid cardiomyocytes, which may be detrimental to adult heart regeneration. However, the gene is evolutionarily conserved, suggesting a beneficial function that has remained obscure. Here, we show that C57BL/6J-inbred Tnni3k mutant mice develop concentric remodeling, characterized by ventricular wall thickening and substantial reduction of cardiomyocyte aspect ratio. This pathology occurs in mice carrying a Tnni3k null allele, a K489R point mutation rendering the protein kinase-dead, or an allele corresponding to human I686T, the most common human non-synonymous TNNI3K variant, which is hypomorphic for kinase activity. Mutant mice develop these conditions in the absence of fibrosis or hypertension, implying a primary cardiomyocyte etiology. In culture, mutant cardiomyocytes were impaired in contractility and calcium dynamics and in protein kinase A signaling in response to isoproterenol, indicating diminished contractile reserve. These results demonstrate a beneficial function of TNNI3K in the adult heart that might explain its evolutionary conservation and imply that human TNNI3K variants, in particular the widespread I686T allele, may convey elevated risk for altered heart geometry and hypertrophy.

Myotonic dystrophy type 1 (DM1) is caused by expansion of a CTG repeat in the DMPK gene, where expansion size and somatic mosaicism correlates with disease severity and age of onset. While it is known that the mismatch repair protein MSH2 contributes to the unstable nature of the repeat, its role on other disease-related features, such as CpG methylation upstream of the repeat, is unknown. In this study, we investigated the effect of an MSH2 knock-down (MSH2KD) on both CTG repeat dynamics and CpG methylation pattern in human embryonic stem cells (hESC) carrying the DM1 mutation. Repeat size in MSH2 wild-type (MSH2WT) and MSH2KD DM1 hESC was determined by PacBio sequencing and CpG methylation by bisulfite massive parallel sequencing. We found stabilization of the CTG repeat concurrent with a gradual loss of methylation upstream of the repeat in MSH2KD cells, while the repeat continued to expand and upstream methylation remained unchanged in MSH2WT control lines. Repeat instability was re-established and biased towards expansions upon MSH2 transgenic re-expression in MSH2KD lines while upstream methylation was not consistently re-established. We hypothesize that the hypermethylation at the mutant DM1 locus is promoted by the MMR machinery and sustained by a constant DNA repair response, establishing a potential mechanistic link between CTG repeat instability and upstream CpG methylation. Our work represents a first step towards understanding how epigenetic alterations and repair pathways connect and contribute to the DM1 pathology.

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q < 0.05) differentially methylated regions (DMRs) that distinguished DS from typical development and idiopathic developmental delay were identified. Machine learning feature selection refined these DMRs to 22 loci. The DS DMRs mapped to genes involved in neurodevelopment, metabolism, and transcriptional regulation. Based on comparisons with previous DS methylation studies and reference epigenomes, the hypermethylated DS DMRs were significantly (q < 0.05) enriched across tissues while the hypomethylated DS DMRs were significantly (q < 0.05) enriched for blood-specific chromatin states. A ~28 kb block of hypermethylation was observed on chromosome 21 in the RUNX1 locus, which encodes a hematopoietic transcription factor whose binding motif was the most significantly enriched (q < 0.05) overall and specifically within the hypomethylated DMRs. Finally, we also identified DMRs that distinguished DS NDBS based on the presence or absence of congenital heart disease (CHD). Together, these results not only demonstrate the utility of low-pass WGBS on NDBS samples for epigenome-wide association studies, but also provide new insights into the early life mechanisms of epigenomic dysregulation resulting from trisomy 21.

The glycogenin knockout mouse is a model of Glycogen Storage Disease type XV. These animals show high perinatal mortality (90%) due to respiratory failure. The lungs of glycogenin-deficient embryos and P0 mice have a lower glycogen content than that of wild-type counterparts. Embryonic lungs were found to have decreased levels of mature surfactant proteins SP-B and SP-C, together with incomplete processing of precursors. Furthermore, non-surviving pups showed collapsed sacculi, which may be linked to a significantly reduced amount of surfactant proteins. A similar pattern was observed in glycogen synthase1-deficient mice, which are devoid of glycogen in the lungs and are also affected by high perinatal mortality due to atelectasis. These results indicate that glycogen availability is a key factor for the burst of surfactant production required to ensure correct lung expansion at the establishment of air breathing. Our findings confirm that glycogen deficiency in lungs can cause respiratory distress syndrome and suggest that mutations in glycogenin and glycogen synthase 1 genes may underlie cases of idiopathic neonatal death.

Autism spectrum disorder (ASD) is characterized by a triad of behavioural impairments including social behaviour. Neuroligin, a trans-synaptic adhesion molecule, has emerged as a penetrant genetic determinant of behavioural traits that signature the neuroatypical behaviours of autism. However, the function of neuroligin in social circuitry and the impact of genetic variation to this gene is not fully understood. Indeed, in animal studies designed to model autism, there remains controversy regarding the role of neuroligin dysfunction in the expression of disrupted social behaviours. The model organism, Caenorhabditis elegans, offers an informative experimental platform to investigate the impact of genetic variants on social behaviour. In a number of paradigms, it has been shown that inter-organismal communication by chemical cues regulates C. elegans social behaviour. We utilize this social behaviour to investigate the effect of autism-associated genetic variants within the social domain of the research domain criteria. We have identified neuroligin as an important regulator of social behaviour and segregate the importance of this gene to the recognition and/or processing of social cues. We also use CRISPR/Cas9 to edit an R-C mutation that mimics a highly penetrant human mutation associated with autism. C. elegans carrying this mutation phenocopy the behavioural dysfunction of a C. elegans neuroligin null mutant, thus confirming its significance in the regulation of animal social biology. This highlights that quantitative behaviour and precision genetic intervention can be used to manipulate discrete social circuits of the worm to provide further insight into complex social behaviour.

Gliomas have substantial mortality to incidence rate ratio and a dismal clinical course. Newer molecular insights, therefore, are imperative to refine glioma diagnosis, prognosis and therapy. Meningioma 1 (MN1) gene is a transcriptional co-regulator implicated in other malignancies, albeit its significance in glioma pathology remains to be explored. IGFBP5 is regulated transcriptionally by MN1 and IGF1 and is associated with higher glioma grade and shorter survival time, prompting us to ascertain their correlation in these tumors. We quantified the expression of MN1, IGFBP5 and IGF1 in 40 glioma samples and examined their interrelatedness. MN1 mRNA-protein inter-correlation and the gene's copy number were evaluated in these tumors. Publicly available TCGA datasets were used to examine the association of MN1 expression levels with patient survival and for validating our findings. We observed MN1 overexpression correlated with low-grade (LGGs) and not high-grade gliomas and is not determined by the copy number alteration of the gene. Notably, gliomas with upregulated MN1 have better overall survival (OS) and progression-free survival (PFS). IGFBP5 expression associated inversely with MN1 expression levels in gliomas but correlated positively with IGF1 expression in only LGGs. This suggests a potential grade-specific interplay between repressive and activating roles of MN1 and IGF1, respectively, in the regulation of IGFBP5. Thus, MN1 overexpression, a promising predictor of OS and PFS in gliomas, may serve as a prognostic biomarker in clinical practice to categorize patients with survival advantage.

Integrating findings from genome-wide association studies with molecular datasets can develop insight into the underlying functional mechanisms responsible for trait-associated genetic variants. We have applied the principles of Mendelian randomization (MR) to investigate whether brain-derived gene expression (n = 1194) may be responsible for mediating the effect of genetic variants on eight cognitive and psychological outcomes (attention deficit hyperactivity disorder (ADHD), Alzheimer's disease, bipolar disorder, depression, intelligence, insomnia, neuroticism and schizophrenia). Transcriptome-wide analyses identified 83 genes associated with at least one outcome (PBonferroni < 6.72 × 10-6), with multiple-trait colocalization also implicating changes to brain-derived DNA methylation at nine of these loci. Comparing effects between outcomes identified evidence of enrichment which may reflect putative causal relationships, such as an inverse relationship between genetic liability towards schizophrenia risk and cognitive ability in later life. Repeating these analyses in whole blood (n = 31 684), we replicated 58.2% of brain-derived effects (based on P < 0.05). Finally, we undertook phenome-wide evaluations at associated loci to investigate pleiotropic effects with 700 complex traits. This highlighted pleiotropic loci such as FURIN (initially implicated in schizophrenia risk (P = 1.05 × 10-7)) which had evidence of an effect on 28 other outcomes, as well as genes which may have a more specific role in disease pathogenesis (e.g. SLC12A5 which only provided evidence of an effect on depression (P = 7.13 × 10-10)). Our results support the utility of whole blood as a valuable proxy for informing initial target identification but also suggest that gene discovery in a tissue-specific manner may be more informative. Finally, non-pleiotropic loci highlighted by our study may be of use for therapeutic translational endeavours.

Erratum in Hum Mol Genet. 2021 Dec 17;31(1):156.

OPA1 mutations are the major cause of dominant optic atrophy (DOA) and the syndromic form DOA plus, pathologies for which there is no established cure. We used a 'drug repurposing' approach to identify FDA-approved molecules able to rescue the mitochondrial dysfunctions induced by OPA1 mutations. We screened two different chemical libraries by using two yeast strains carrying the mgm1I322M and the chim3P646L mutations, identifying 26 drugs able to rescue their oxidative growth phenotype. Six of them, able to reduce the mitochondrial DNA instability in yeast, have been then tested in Opa1 deleted mouse embryonic fibroblasts expressing the human OPA1 isoform 1 bearing the R445H and D603H mutations. Some of these molecules were able to ameliorate the energetic functions and/or the mitochondrial network morphology, depending on the type of OPA1 mutation. The final validation has been performed in patients' fibroblasts, allowing to select the most effective molecules. Our current results are instrumental to rapidly translating the findings of this drug repurposing approach into clinical trial for DOA and other neurodegenerations caused by OPA1 mutations.

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is the most common pediatric musculoskeletal disorder, affecting ~3% of the population worldwide. However, its genetic bases and tissues of origin remain largely unknown. Several genome-wide association studies (GWAS) have implicated nucleotide variants in non-coding sequences that control genes with important roles in cartilage, muscle, bone, connective tissue and intervertebral disks (IVDs) as drivers of AIS susceptibility. Here, we set out to define the expression of AIS-associated genes and active regulatory elements by performing RNA-seq and chromatin immunoprecipitation-sequencing against H3 lysine 27 acetylation in these tissues in mouse and human. Our study highlights genetic pathways involving AIS-associated loci that regulate chondrogenesis, IVD development and connective tissue maintenance and homeostasis. In addition, we identify thousands of putative AIS-associated regulatory elements which may orchestrate tissue-specific expression in musculoskeletal tissues of the spine. Quantification of enhancer activity of several candidate regulatory elements from our study identifies three functional enhancers carrying AIS-associated GWAS SNPs at the ADGRG6 and BNC2 loci. Our findings provide a novel genome-wide catalog of AIS-relevant genes and regulatory elements and aid in the identification of novel targets for AIS causality and treatment.

Adrenocortical Carcinoma (ACC) is a rare endocrine tumor with poor overall prognosis and 1.5-fold overrepresentation in females. In children, ACC is associated with inherited cancer syndromes with 50-80% of childhood-ACC associated with TP53 germline variants. ACC in adolescents and young adults (AYA) is rarely due to germline TP53, IGF2, PRKAR1A and MEN1 variants. We analyzed exome sequencing data from 21 children (<15y), 32 AYA (15-39y), and 60 adults (>39y) with ACC, and retained all pathogenic, likely pathogenic, and highly prioritized variants of uncertain significance. We engineered a stable lentiviral-mutant ACC cell line, harboring an EGFR variant (p.Asp1080Asn) from a 21-year-old female without germline-TP53-variant and with aggressive ACC. We found that 4.8% of the children (P = 0.004) and 6.2% of AYA (P < 0.0001), all-female participants, harbored germline EGFR variants, compared to only 0.3% of the control group. Expanding our analysis to the RTK-RAS-MAPK pathway, we found that the RTK genes have the highest number of highly prioritized germline variants in these individuals amongst all three arms of this pathway. We showed EGFR mutant cells migrate faster and are characterized by a stem-like phenotype compared to wild type cells. While EGFR inhibitors did not affect the stemness of mutant cells, Sunitinib, a multireceptor tyrosine kinase inhibitor, significantly reduced their stem-like behavior. Our data suggest that EGFR could be a novel underlying germline predisposition factor for ACC, especially in the Childhood-AYA (C-AYA) population. Further clinical validation can improve precision oncology management of this disease, which is known to have limited therapeutic options.

Spastic paraplegia 35 (SPG35) (OMIM: 612319) or fatty acid hydroxylase-associated neurodegeneration (FAHN) is caused by deficiency of fatty acid 2-hydroxylase (FA2H). This enzyme synthesizes sphingolipids containing 2-hydroxylated fatty acids, which are particularly abundant in myelin. Fa2h-deficient (Fa2h-/-) mice develop symptoms reminiscent of the human disease and therefore serve as animal model of SPG35. In order to understand further the pathogenesis of SPG35, we compared the proteome of purified CNS myelin isolated from wild type and Fa2h-/- mice at different time points of disease progression using tandem mass tag labeling. Data analysis with a focus on myelin membrane proteins revealed a significant increase of the oligodendrocytic myelin paranodal and inner loop protein (Opalin) in Fa2h-/- mice, whereas the concentration of other major myelin proteins was not significantly changed. Western blot analysis revealed an almost 6-fold increase of Opalin in myelin of Fa2h-/- mice aged 21-23 months. A concurrent unaltered Opalin gene expression suggested a decreased turnover of the Opalin protein in Fa2h-/- mice. Supporting this hypothesis, Opalin protein half-life was reduced significantly when expressed in CHO cells synthesizing 2-hydroxylated sulfatide, compared to cells synthesizing only non-hydroxylated sulfatide. Degradation of Opalin was inhibited by inhibitors of lysosomal degradation but unaffected by proteasome inhibitors. Taken together, these results reveal a new function of 2-hydroxylated sphingolipids namely affecting the turnover of a myelin membrane protein. This may play a role in the pathogenesis of SPG35.

Mutations in the GDAP1 gene cause Charcot-Marie-Tooth (CMT) neuropathy. GDAP1 is an atypical glutathione S-transferase (GST) of the outer mitochondrial membrane and the mitochondrial membrane contacts with the endoplasmic reticulum (MAMs). Here, we investigate the role of this GST in the autophagic flux and the membrane contact sites (MCSs) between mitochondria and lysosomes in the cellular pathophysiology of GDAP1 deficiency. We demonstrate that GDAP1 participates in basal autophagy and that its depletion affects LC3 and PI3P biology in autophagosome biogenesis and membrane trafficking from MAMs. GDAP1 also contributes to the maturation of lysosome by interacting with PYKfyve kinase, a pH-dependent master lysosomal regulator. GDAP1 deficiency causes giant lysosomes with hydrolytic activity, a delay in the autophagic lysosome reformation, and TFEB activation. Notably, we found that GDAP1 interacts with LAMP-1, which supports that GDAP1-LAMP-1 is a new tethering pair of mitochondria and lysosome membrane contacts. We observed mitochondria-lysosome MCSs in soma and axons of cultured mouse embryonic motor neurons and human neuroblastoma cells. GDAP1 deficiency reduces the MCSs between these organelles, causes mitochondrial network abnormalities, and decreases levels of cellular glutathione (GSH). The supply of GSH-MEE suffices to rescue the lysosome membranes and the defects of the mitochondrial network, but not the interorganelle MCSs nor early autophagic events. Overall, we show that GDAP1 enables the proper function of mitochondrial MCSs in both degradative and nondegradative pathways, which could explain primary insults in GDAP1-related CMT pathophysiology, and highlights new redox-sensitive targets in axonopathies where mitochondria and lysosomes are involved.

Mutations in retinitis pigmentosa GTPase regulator (RPGR) cause severe retinal ciliopathy, X-linked retinitis pigmentosa. Although two major alternatively spliced isoforms, RPGRex1-19 and RPGRORF15, are expressed, the relative importance of these isoforms in disease pathogenesis is unclear. Here, we analyzed fibroblast samples from eight patients and found that all of them form longer cilia than normal controls, albeit to different degrees. Although all mutant RPGRORF15 messenger RNAs (mRNAs) are unstable, their steady-state levels were similar or higher than those in the control cells, suggesting there may be increased transcription. Three of the fibroblasts that had higher levels of mutant RPGRORF15 mRNA also exhibited significantly higher levels of RPGRex1-19 mRNA. Four samples with unaltered RPGRex1-19 levels carried mutations in RPGRORF15 that resulted in this isoform being relatively less stable. Thus, in all cases, the RPGRex1-19/RPGRORF15 isoform ratio was increased, and this was highly correlative to the cilia extension defect. Moreover, overexpression of RPGRex1-19 (mimicking the increase in RPGRex1-19 to RPGRORF15 isoform ratio) or RPGRORF15 (mimicking reduction of the ratio) resulted in significantly longer or shorter cilia, respectively. Notably, the cilia length defect appears to be attributable to both the loss of the wild-type RPGRORF15 protein and to the higher levels of the RPGRex1-19 isoform, indicating that the observed defect is due to the altered isoform ratios. These results suggest that maintaining the optimal RPGRex1-9 to RPGRORF15 ratio is critical for cilia growth and that designing strategies that focus on the best ways to restore the RPGRex1-19/RPGRORF15 ratio may lead to better therapeutic outcomes.

Monoallelic gene expression at the Igf2/H19 locus is controlled by paternal allele-specific DNA methylation of the imprinting control region (H19 ICR) that is established during spermatogenesis. We demonstrated that the H19 ICR fragment in transgenic mice acquires allele-specific methylation only after fertilization, which is essential for maintaining its allelic methylation during early embryogenesis. We identified a DNA element required for establishing postfertilization methylation within a 118 bp (m118) region. A previously generated knock-in mouse whose endogenous H19 ICR was substituted with the human H19 ICR (hIC1; 4.8 kb) sequence revealed that the hIC1 sequence was partially methylated in sperm, although this methylation was lost by the blastocyst stage, which we assume is due to a lack of an m118-equivalent sequence in the hIC1 transgene. To identify a cis sequence involved in postfertilization methylation within the hIC1 region, we generated three transgenic mouse lines (TgM): one carrying an 8.8 kb hIC1 sequence joined to m118 (hIC1+m118), one with the 8.8 kb hIC1 and one with the 5.8 kb hIC1 sequence joined to m118 (hIC1-3'+m118). We found that the hIC1-3' region was resistant to de novo DNA methylation throughout development. In contrast, the 5' portion of the hIC1 (hIC1-5') in both hIC1+m118 and hIC1 TgM were preferentially methylated on the paternal allele only during preimplantation. As DNA methylation levels were higher in hIC1+m118, the m118 sequence could also induce imprinted methylation of the human sequence. Most importantly, the hIC1-5' sequence appears to possess an activity equivalent to that of m118.

Talipes equinovarus (clubfoot, TEV) is a congenital rotational foot deformity occurring in 1 per 1000 births with increased prevalence in males compared with females. The genetic etiology of isolated clubfoot (iTEV) remains unclear. Using a genome-wide association study, we identified a locus within FSTL5, encoding follistatin-like 5, significantly associated with iTEV. FSTL5 is an uncharacterized gene whose potential role in embryonic and postnatal development was previously unstudied. Utilizing multiple model systems, we found that Fstl5 was expressed during later stages of embryonic hindlimb development, and, in mice, expression was restricted to the condensing cartilage anlage destined to form the limb skeleton. In the postnatal growth plate, Fstl5 was specifically expressed in prehypertrophic chondrocytes. As Fstl5 knockout rats displayed no gross malformations, we engineered a conditional transgenic mouse line (Fstl5LSL) to overexpress Fstl5 in skeletal osteochondroprogenitors. We observed that hindlimbs were slightly shorter and that bone mineral density was reduced in adult male, but not female, Prrx1-cre;Fstl5LSL mice compared with control. No overt clubfoot-like deformity was observed in Prrx1-cre;Fstl5LSL mice, suggesting FSTL5 may function in other cell types to contribute to iTEV pathogenesis. Interrogating published mouse embryonic single-cell expression data showed that Fstl5 was expressed in cell lineage subclusters whose transcriptomes were associated with neural system development. Moreover, our results suggest that lineage-specific expression of the Fstl genes correlates with their divergent roles as modulators of transforming growth factor beta and bone morphogenetic protein signaling. Results from this study associate FSTL5 with iTEV and suggest a potential sexually dimorphic role for Fstl5 in vivo.