Osteoarthritis (OA) is a polygenic disease of older people resulting in the breakdown of cartilage within articular joints. Although a leading cause of disability, there are no disease-modifying therapies. Evidence is emerging to support the origins of OA in skeletogenesis. Whilst methylation QTLs (mQTLs) co-localizing with OA GWAS signals have been identified in aged human cartilage and used to identify effector genes and variants, such analyses have never been conducted during human development. Here, for the first time, we have investigated the developmental origins of OA genetic risk at seven well-characterized OA risk loci, comprising 39 OA-mQTL CpGs, in human foetal limb (FL) and cartilage (FC) tissues using a range of molecular genetic techniques. We compared our results to aged cartilage samples (AC) and identified significant OA-mQTLs at 14 CpGs and 29 CpGs in FL and FC tissues, respectively. Differential methylation was observed at 26 sites between foetal and aged cartilage, with the majority becoming actively hypermethylated in old age. Notably, 6/9 OA effector genes showed allelic expression imbalances during foetal development. Finally, we conducted ATAC-sequencing in cartilage from the developing and aged hip and knee to identify accessible chromatin regions, and found enrichment for transcription factor binding motifs including SOX9 and FOS/JUN. For the first time, we have demonstrated the activity of OA-mQTLs and expression imbalance of OA effector genes during skeletogenesis. We show striking differences in the spatiotemporal function of these loci, contributing to our understanding of OA aetiology, with implications for the timing and strategy of pharmacological interventions.

Germline gain-of-function missense variants in the pore-forming Cav1.3 α1-subunit (CACNA1D gene) confer high risk for a severe neurodevelopmental disorder with or without endocrine symptoms. Here we report a four weeks old new-born with the novel de novo missense variant F747S with a so far not described prominent jittering phenotype in addition to symptoms previously reported for CACNA1D mutations including developmental delay, elevated aldosterone level, and transient hypoglycemia. We confirmed the pathogenicity of this variant in whole-cell patch-clamp experiments with wild-type and F747S mutant channels heterologously expressed together with α2δ1 and cytosolic β3 or membrane-bound β2a subunits. Mutation F747S caused the quantitatively largest shift in the voltage dependence of activation (-28 mV) reported so far for CACNA1D germline mutations. It also shifted inactivation to more negative voltages, slowed the time course of current inactivation and slowed current deactivation upon repolarization with both co-expressed β-subunits. In silico modelling and molecular docking simulations revealed that this gain-of-function phenotype can be explained by formation of a novel inter-domain hydrogen bond between mutant residues S747 (IIS6) with N1145 (IIIS6) stabilizing selectively the activated open channel state. F747S displayed 2-6-fold increased sensitivity for the L-type Ca2+ channel blocker isradipine compared to wild-type. Our data confirm the pathogenicity of the F747S variant with very strong gain-of-function gating changes, which may contribute to the novel jittering phenotype. Increased sensitivity for isradipine suggests this drug for potential symptomatic off-label treatment for carriers of this mutation.

Olig2 is a basic helix-loop-helix transcription factor that plays a critical role in the central nervous system. It directs the specification of motor neurons and oligodendrocyte precursor cells (OPCs) from neural progenitors and the subsequent maturation of OPCs into myelin-forming oligodendrocytes (OLs). It is also required for the development of astrocytes. Despite a decade-long search, enhancers that regulate the expression of Olig2 remain elusive. We have recently developed an innovative method that maps promoter-distal enhancers to genes in a principled manner. Here, we applied it to Olig2 in the context of OL lineage cells, uncovering an OL enhancer for it (termed Olig2-E1). Silencing Olig2-E1 by CRISPRi epigenome editing significantly downregulated Olig2 expression. Luciferase assay and ATAC-seq and ChIP-seq data show that Olig2-E1 is an OL-specific enhancer that is conserved across human, mouse, and rat. Hi-C data reveal that Olig2-E1 physically interacts with OLIG2 and suggest that this interaction is specific to OL lineage cells. In sum, Olig2-E1 is an evolutionarily conserved OL-specific enhancer that drives the expression of Olig2.

Maintaining protein lipoylation is vital for cell metabolism. The H-protein encoded by GCSH has a dual role in protein lipoylation required for bioenergetic enzymes including pyruvate dehydrogenase and 2-ketoglutarate dehydrogenase, and in the one-carbon metabolism through its involvement in glycine cleavage enzyme system, intersecting two vital roles for cell survival. Here we report six patients with biallelic pathogenic variants in GCSH and a broad clinical spectrum ranging from neonatal fatal glycine encephalopathy to an attenuated phenotype of developmental delay, behavioral problems, limited epilepsy, and variable movement problems. The mutational spectrum includes one insertion c.293-2_293-1insT, one deletion c.122_(228 + 1_229-1) del, one duplication of exons 4 and 5, one nonsense variant p.Gln76*and four missense p.His57Arg, p.Pro115Leu, and p.Thr148Pro and the previously described p.Met1?. Via functional studies in patient's fibroblasts, molecular modelling, expression analysis in GCSH knock-down COS7 cells and yeast, and in vitro protein studies, we demonstrate for the first time that most variants identified in our cohort produced a hypomorphic effect on both mitochondrial activities, protein lipoylation and glycine metabolism, causing combined deficiency whereas some missense variants affect primarily one function only. The clinical features of the patients reflect the impact of the GCSH changes on any of the two functions analyzed. Our analysis illustrates the complex interplay of functional and clinical impact when pathogenic variants affect a multifunctional protein involved in two metabolic pathways and emphasizes the value of the functional assays to select the treatment and investigate new personalized options.

TFIIH is a complex essential for transcription of protein-coding genes by RNA polymerase II, DNA repair of UV-lesions and transcription of rRNA by RNA polymerase I. Mutations in TFIIH cause the cancer prone DNA-repair disorder xeroderma pigmentosum (XP) and the developmental and premature aging disorders trichothiodystrophy (TTD) and Cockayne syndrome (CS). 50% of TTD cases are caused by TFIIH mutations. Using TFIIH mutant patient cells from TTD and XP subjects we can show that the stress-sensitivity of the proteome is reduced in TTD, but not in XP. Using three different methods to investigate the accuracy of protein synthesis by the ribosome, we demonstrate that translational fidelity of the ribosomes of TTD, but not XP cells, is decreased. The process of ribosomal synthesis and maturation is affected in TTD cells and can lead to instable ribosomes. Isolated ribosomes from TTD patients show an elevated error rate when challenged with oxidized mRNA, explaining the oxidative hypersensitivity of TTD cells. Treatment of TTD cells with N-acetyl cysteine normalized the increased translational error-rate and restored translational fidelity. Here we describe a pathomechanism that might be relevant for our understanding of impaired development and aging-associated neurodegeneration.

Inflammatory bowel disease (IBD) is an immune-mediated chronic intestinal disorder with major phenotypes; ulcerative colitis (UC), and Crohn's disease (CD). Multiple studies have identified over 240 IBD susceptibility loci. However, most studies have centered on European (EUR) and East Asian (EAS) populations. The prevalence of IBD in non-EUR, including African Americans (AAs), has risen in recent years. Here we present the first attempt to identify loci in AAs using a trans-ancestry Bayesian approach (MANTRA) accounting for heterogeneity between diverse ancestries while allowing for the similarity between closely related populations. We meta-analyzed GWAS and Immunochip data from a 2015 EUR meta-analysis of 38 155 IBD cases and 48 485 controls and EAS Immunochip study of 2824 IBD cases and 3719 controls, and our recent AA IBD GWAS of 2345 cases and 5002 controls. Across the major IBD phenotypes, we found significant evidence for 92% of the 2015 loci but also 3 IBD loci only established in latter studies. We detected 20 novel loci, all containing immunity-related genes or genes with other evidence for IBD or immune-mediated disease relevance: PLEKHG5;TNFSFR25 (encoding death receptor 3, receptor for TNFSF15 gene product TL1A), XKR6, ELMO1, BC021024;PI4KB;PSMD4, and APLP1 for IBD; AUTS2, XKR6, OSER1, TET2;AK094561, BCAP29, and APLP1 for CD; and GABBR1;MOG, DQ570892, SPDEF;ILRUN, SMARCE1;CCR7;KRT222;KRT24;KRT25, ANKS1A;TCP11, IL7, LRRC18;WDFY4, XKR6 and TNFSF4 for UC. Our study highlights the value of combining low-powered genomic studies from understudied populations of diverse ancestral backgrounds together with a high-powered study to enable novel locus discovery, including potentially important therapeutic IBD gene targets.

Obesity has a highly complex genetic architecture, making it difficult to understand the genetic mechanisms, despite the large number of discovered loci via genome-wide association studies (GWAS). Omics techniques have provided a better resolution to view this problem. As a proxy of cell-level biology, extracellular vesicles (EVs) are useful for studying cellular regulation of complex phenotypes such as obesity. Here, in a well-established Scottish cohort, we utilized a novel technology to detect surface proteins across millions of single EVs in each individual's plasma sample. Integrating the results with established obesity GWAS, we inferred 78 types of EVs carrying one or two of 12 surface proteins to be associated with adiposity-related traits such as waist circumference. We then verified that particular EVs' abundance is negatively correlated with body adiposity, while no association with lean body mass. We also revealed that genetic variants associated with protein-specific EVs capture 2-4-fold heritability enrichment for blood cholesterol levels. Our findings provide evidence that EVs with specific surface proteins have phenotypic and genetic links to obesity and blood lipids, respectively, guiding future EV biomarker research.

Parkinson's disease (PD) is a neurological disorder with complex interindividual etiology that is becoming increasingly prevalent worldwide. Elevated alpha-synuclein levels can increase risk of PD and may influence epigenetic regulation of PD pathways. Here, we report genome-wide DNA methylation and hydroxymethylation alterations associated with overexpression of two PD-linked alpha-synuclein variants (wild-type and A30P) in LUHMES cells differentiated to dopaminergic neurons. Alpha-synuclein altered DNA methylation at thousands of CpGs and DNA hydroxymethylation at hundreds of CpGs in both genotypes, primarily in locomotor behavior and glutamate signaling pathway genes. In some cases, epigenetic changes were associated with transcription. SMITE network analysis incorporating H3K4me1 ChIP-seq to score DNA methylation and hydroxymethylation changes across promoters, enhancers, and gene bodies confirmed epigenetic and transcriptional deregulation of glutamate signaling modules in both genotypes. Our results identify distinct and shared impacts of alpha-synuclein variants on the epigenome, and associate alpha-synuclein with the epigenetic etiology of PD.

Rabenosyn (RBSN) is a conserved endosomal protein necessary for regulating internalized cargo. Here, we present clinical, genetic, cellular and biochemical evidence that two distinct RBSN missense variants are responsible for a novel Mendelian disorder consisting of progressive muscle weakness, facial dysmorphisms, ophthalmoplegia and intellectual disability. Using exome sequencing, we identified recessively acting germline alleles p.Arg180Gly and p.Gly183Arg, which are both situated in the FYVE domain of RBSN. We find that these variants abrogate binding to its cognate substrate phosphatidylinositol 3-phosphate (PI3P) and thus prevent its translocation to early endosomes. Although the endosomal recycling pathway was unaltered, mutant p.Gly183Arg patient fibroblasts show accumulation of cargo tagged for lysosomal degradation. Our results suggest that these variants are separation-of-function alleles, which cause a delay in endosomal maturation without affecting cargo recycling. We conclude that distinct germline mutations in RBSN cause non-overlapping phenotypes with specific and discrete endolysosomal cellular defects.

Kabuki syndrome is frequently caused by loss-of-function mutations in one allele of histone 3 lysine 4 (H3K4) methyltransferase KMT2D and is associated with problems in neurological, immunological and skeletal system development. We generated heterozygous KMT2D knockout and Kabuki patient-derived cell models to investigate the role of reduced dosage of KMT2D in stem cells. We discovered chromosomal locus-specific alterations in gene expression, specifically a 110 Kb region containing Synaptotagmin 3 (SYT3), C-Type Lectin Domain Containing 11A (CLEC11A), Chromosome 19 Open Reading Frame 81 (C19ORF81) and SH3 And Multiple Ankyrin Repeat Domains 1 (SHANK1), suggesting locus-specific targeting of KMT2D. Using whole genome histone methylation mapping, we confirmed locus-specific changes in H3K4 methylation patterning coincident with regional decreases in gene expression in Kabuki cell models. Significantly reduced H3K4 peaks aligned with regions of stem cell maps of H3K27 and H3K4 methylation suggesting KMT2D haploinsufficiency impact bivalent enhancers in stem cells. Preparing the genome for subsequent differentiation cues may be of significant importance for Kabuki-related genes. This work provides a new insight into the mechanism of action of an important gene in bone and brain development and may increase our understanding of a specific function of a human disease-relevant H3K4 methyltransferase family member.

Wilson's disease (WD) is a copper metabolic disorder caused by a defective ATP7B function. Conventional therapies cause severe side effects and significant variation in efficacy, according to cohort studies. Thus, exploring new therapeutic approaches to prevent progression to liver failure is urgent. To study the physiology and pathology of WD, immortalized cell lines and rodent WD models have been used conventionally; however, a large gap remains among different species as well as in genetic backgrounds among individuals. We generated induced pluripotent stem cells (iPSCs) from four WD patients carrying compound heterozygous mutations in the ATP7B gene. ATP7B loss- and gain-of-functions were further manifested with ATP7B-deficient iPSCs and heterozygously corrected R778L WD patient-derived iPSCs using CRISPR-Cas9-based gene editing. Although the expression of ATP7B protein varied among WD-specific hepatocytes differentiated from these iPSCs, the expression and secretion of ceruloplasmin (Cp), a downstream copper carrier in plasma, were consistently decreased in WD patient-derived and ATP7B-deficient hepatocytes. A transcriptome analysis detected abnormalities in the retinoid signaling pathway and lipid metabolism in WD-specific hepatocytes. Drug screening using WD patient-derived hepatocytes identified retinoids as promising candidates for rescuing Cp secretion. All-trans retinoic acid also alleviates reactive oxygen species production induced by lipid accumulation in WD-specific hepatocytes treated with oleic acid. These patient-derived iPSC-based hepatic models function as effective platforms for the development of potential therapeutics for hepatic steatosis in WD and other fatty liver diseases.

Gout is of particularly high prevalence in the Māori and Pacific (Polynesian) populations of Aotearoa New Zealand (NZ). Here, we investigated the contribution of common population-specific copy number variation (CNV) to gout in the Aotearoa NZ Polynesian population. Microarray-generated genome-wide genotype data from Aotearoa NZ Polynesian individuals with (n = 1196) and without (n = 1249) gout were analyzed. Comparator population groups were 552 individuals of European ancestry and 1962 of Han Chinese ancestry. Levels of circulating major histocompatibility complex (MHC) class I polypeptide-related sequence A (MICA) were measured by enzyme-linked immunosorbent assay. Fifty-four CNV regions (CNVRs) appearing in at least 10 individuals were detected, of which seven common (>2%) CNVRs were specific to or amplified in Polynesian people. A burden test of these seven revealed associations of insertion/deletion with gout (odds ratio (OR) 95% confidence interval [CI] = 1.80 [1.01; 3.22], P = 0.046). Individually testing of the seven CNVRs for association with gout revealed nominal association of CNVR1 with gout in Western Polynesian (Chr6: 31.36-31.45 Mb, OR = 1.72 [1.03; 2.92], P = 0.04), CNVR6 in the meta-analyzed Polynesian sample sets (Chr1: 196.75-196.92 Mb, OR = 1.86 [1.16; 3.00], P = 0.01) and CNVR9 in Western Polynesian (Chr1: 189.35-189.54 Mb, OR = 2.75 [1.15; 7.13], P = 0.03). Analysis of European gout genetic association data demonstrated a signal of association at the CNVR1 locus that was an expression quantitative trait locus for MICA. The most common CNVR (CNVR1) includes deletion of the MICA gene, encoding an immunomodulatory protein. Expression of MICA was reduced in the serum of individuals with the deletion. In summary, we provide evidence for the association of CNVR1 containing MICA with gout in Polynesian people, implicating class I MHC-mediated antigen presentation in gout.

Genome-wide association studies have identified a growing number of single nucleotide polymorphisms (SNPs) associated with childhood acute lymphoblastic leukemia (ALL), yet the functional roles of most SNPs are unclear. Multiple lines of evidence suggest that epigenetic mechanisms may mediate the impact of heritable genetic variation on phenotypes. Here, we investigated whether DNA methylation mediates the effect of genetic risk loci for childhood ALL. We performed an epigenome-wide association study (EWAS) including 808 childhood ALL cases and 919 controls from California-based studies using neonatal blood DNA. For differentially methylated CpG positions (DMPs), we next conducted association analysis with 23 known ALL risk SNPs followed by causal mediation analyses addressing the significant SNP-DMP pairs. DNA methylation at CpG cg01139861, in the promoter region of IKZF1, mediated the effects of the intronic IKZF1 risk SNP rs78396808, with the average causal mediation effect (ACME) explaining ~30% of the total effect (ACME P = 0.0031). In analyses stratified by self-reported race/ethnicity, the mediation effect was only significant in Latinos, explaining ~41% of the total effect of rs78396808 on ALL risk (ACME P = 0.0037). Conditional analyses confirmed the presence of at least three independent genetic risk loci for childhood ALL at IKZF1, with rs78396808 unique to non-European populations. We also demonstrated that the most significant DMP in the EWAS, CpG cg13344587 at gene ARID5B (P = 8.61 × 10-10), was entirely confounded by the ARID5B ALL risk SNP rs7090445. Our findings provide new insights into the functional pathways of ALL risk SNPs and the DNA methylation differences associated with risk of childhood ALL.

Isolated thoracic aortic aneurysms (TAAs) are asymptomatic before dissection or rupture and heterogeneous in clinical phenotype. It is urgent need but difficult to identify individuals at high risk to enable enhanced screening or preventive therapies. Because TAAs have a genetic component, one possible approach is to stratify individuals based on inherited DNA variations. Here, we constructed an integrated exome risk score (ERS) based on both common and rare variants found in whole-exome sequencing through a machine-learning framework in discovery population consisting of 551 cases and 1071 controls. We evaluated the performance of the ERS in an independent population including 151 cases and 779 controls with a raw odds ratio (OR) per 1 standard deviation (SD) = 1.95 and area under the receiver operating characteristic curve (AUC) = 0.680. When adjusted by gender and the first four principal components, OR per SD = 1.68 and AUC reached 0.783. Individuals in the top 20% of ERS distribution had an OR of 3.20 compared with others. Finally, we found that individuals with top 20% ERS developed TAA at a younger age (P = 0.002) and with a larger diameter (P = 0.016) compared with lower ERS, and were more likely to suffer from aortic root aneurysms (P = 0.009). Our analysis provides a global view of the genetic components of isolated TAA. The exome score developed and evaluated here is the first polygenic risk score for TAA and is a promising predictor of disease risk and severity, which will facilitate the implementation of the risk-reduction strategies.

Factor I (FI) is a serine protease inhibitor of the complement system. Heterozygous rare genetic variants in complement factor I (CFI) are associated with advanced age-related macular degeneration (AMD). The clinical impact of these variants is unknown since a majority have not been functionally characterized and are classified as 'variants of uncertain significance' (VUS). This study assessed the functional significance of VUS in CFI. Our previous cross-sectional study using a serum-based assay demonstrated that CFI variants in advanced AMD can be categorized into three types. Type 1 variants cause a quantitative deficiency of FI. Type 2 variants demonstrate a qualitative deficiency. However, Type 3 variants consist of VUS that are less dysfunctional than Types 1 and 2 but are not as biologically active as wild type (WT). In this study, we employed site-directed mutagenesis followed by expression of the recombinant variant and a comprehensive set of functional assays to characterize nine Type 3 variants that were identified in 37 individuals. Our studies establish that the expression of the recombinant protein compared with WT is reduced for R202I, Q217H, S221Y and G263V. Further, G362A and N536K, albeit expressed normally, have significantly less cofactor activity. These results led to re-categorization of CFI variants R202I, Q217H, S221Y and G263V as Type 1 variants and to reclassification of N536K and G362A as Type 2. The variants K441R, Q462H and I492L showed no functional defect and remained as Type 3. This study highlights the utility of an in-depth biochemical analysis in defining the pathologic and clinical implications of complement variants underlying AMD.

Erratum in Hum Mol Genet. 2022 Dec 22;:

Humans present remarkable diversity in their mitochondrial DNA (mtDNA) in terms of variants across individuals as well as across tissues and even cells within one person. We have investigated the timing of the first appearance of this variant-driven mosaicism. For this, we deep-sequenced the mtDNA of 254 oocytes from 85 donors, 158 single blastomeres of 25 day-3 embryos, 17 inner cell mass and trophectoderm samples of 7 day-5 blastocysts, 142 bulk DNA and 68 single cells of different adult tissues. We found that day-3 embryos present blastomeres that carry variants only detected in that cell, showing that mtDNA mosaicism arises very early in human development. We classified the mtDNA variants based on their recurrence or uniqueness across different samples. Recurring variants had higher heteroplasmic loads and more frequently resulted in synonymous changes or were located in non-coding regions than variants unique to one oocyte or single embryonic cell. These differences were maintained through development, suggesting that the mtDNA mosaicism arising in the embryo is maintained into adulthood. We observed a decline in potentially pathogenic variants between day 3 and day 5 of development, suggesting early selection. We propose a model in which closely clustered mitochondria carrying specific mtDNA variants in the ooplasm are asymmetrically distributed throughout the cell divisions of the preimplantation embryo, resulting in the earliest form of mtDNA mosaicism in human development.

A growing body of evidence suggests that nuclear alpha-synuclein (αSyn) plays a role in the pathogenesis of Parkinson's disease (PD). However, this question has been difficult to address as controlling the localization of αSyn in experimental systems often requires protein overexpression, which affects its aggregation propensity. To overcome this, we engineered SncaNLS mice, which localize endogenous αSyn to the nucleus. We characterized these mice on a behavioral, histological and biochemical level to determine whether the increase of nuclear αSyn is sufficient to elicit PD-like phenotypes. SncaNLS mice exhibit age-dependent motor deficits and altered gastrointestinal function. We found that these phenotypes were not linked to αSyn aggregation or phosphorylation. Through histological analyses, we observed motor cortex atrophy in the absence of midbrain dopaminergic neurodegeneration. We sampled cortical proteomes of SncaNLS mice and controls to determine the molecular underpinnings of these pathologies. Interestingly, we found several dysregulated proteins involved in dopaminergic signaling, including Darpp32, Pde10a and Gng7, which we further confirmed was decreased in cortical samples of the SncaNLS mice compared with controls. These results suggest that chronic endogenous nuclear αSyn can elicit toxic phenotypes in mice, independent of its aggregation. This model raises key questions related to the mechanism of αSyn toxicity in PD and provides a new model to study an underappreciated aspect of PD pathogenesis.

Pathogenesis of the inherited neurodegenerative disorder Huntington's disease (HD) is progressive with a long presymptomatic phase in which subtle changes occur up to 15 years before the onset of symptoms. Thus, there is a need for early, functional biomarker to better understand disease progression and to evaluate treatment efficacy far from onset. Recent studies have shown that white matter may be affected early in mutant HTT gene carriers. A previous study performed on 12 months old Ki140CAG mice showed reduced glutamate level measured by Chemical Exchange Saturation Transfer of glutamate (gluCEST), especially in the corpus callosum. In this study, we scanned longitudinally Ki140CAG mice with structural MRI, diffusion tensor imaging, gluCEST and magnetization transfer imaging, in order to assess white matter integrity over the life of this mouse model characterized by slow progression of symptoms. Our results show early defects of diffusion properties in the anterior part of the corpus callosum at 5 months of age, preceding gluCEST defects in the same region at 8 and 12 months that spread to adjacent regions. At 12 months, frontal and piriform cortices showed reduced gluCEST, as well as the pallidum. MT imaging showed reduced signal in the septum at 12 months. Cortical and striatal atrophy then appear at 18 months. Vulnerability of the striatum and motor cortex, combined with alterations of anterior corpus callosum, seems to point out the potential role of white matter in the brain dysfunction that characterizes HD and the pertinence of gluCEST and DTI as biomarkers in HD.

Auditory synaptopathy/neuropathy (AS/AN) is a distinct type of sensorineural hearing loss in which the cochlear sensitivity to sound (i.e. active cochlear amplification by outer hair cells) is preserved whereas sound encoding by inner hair cells and/or auditory nerve fibers is disrupted due to genetic or environmental factors. Autosomal-dominant auditory neuropathy type 2 (AUNA2) was linked either to chromosomal bands 12q24 or 13q34 in a large German family in 2017. By whole genome sequencing, we now detected a 5500 bp deletion in ATP11A on chromosome 13q34 segregating with the phenotype in this family. ATP11A encodes a P-type ATPase that translocates phospholipids from the exoplasmic to the cytoplasmic leaflet of the plasma membrane. The deletion affects both isoforms of ATP11A and activates a cryptic splice site leading to the formation of an alternative last exon. ATP11A carrying the altered C-terminus loses its flippase activity for phosphatidylserine. Atp11a is expressed in fibers and synaptic contacts of the auditory nerve and in the cochlear nucleus in mice and conditional Atp11a knockout mice show a progressive reduction of the spiral ganglion neuron compound action potential, recapitulating the human phenotype of auditory neuropathy. By combining whole genome sequencing, immunohistochemistry, in vitro functional assays and generation of a mouse model, we could thus identify a partial deletion of ATP11A as the genetic cause of AUNA2.

BACKGROUND: Cysteamine is currently the only therapy for nephropathic cystinosis. It significantly improves life expectancy and delays progression to end-stage kidney disease, however, it cannot prevent it. Unfortunately, compliance to therapy is often weak, particularly during adolescence. Therefore, finding better treatments is a priority in the field of cystinosis. Previously, we found that genistein, an isoflavone particularly enriched in soy, can revert part of the cystinotic cellular phenotype that is not sensitive to cysteamine in vitro. METHODS: To test the effects of genistein in vivo, we fed 2-month-old wild type and Ctns-/- female mice with either a control diet, a genistein-containing diet or a cysteamine-containing diet for 14 months. RESULTS: Genistein (160 mg/kg/day) did not affect the growth of the mice or hepatic functionality. Compared with untreated mice at 16 months, Ctns-/- mice fed with genistein had lower cystine concentrations in their kidneys, reduced formation of cystine crystals, a smaller number of LAMP1-positive structures, and an overall better-preserved parenchymal architecture. Cysteamine (400 mg/kg/day) was efficient in reverting the lysosomal phenotype and in preventing the development of renal lesions. CONCLUSIONS: These preclinical data indicate that genistein ameliorates kidney injury resulting from cystinosis with no side effects. Genistein therapy represents a potential treatment to improve the outcome for patients with cystinosis.

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is a congenital condition characterized by aplasia or hypoplasia of the uterus and vagina in women with a 46,XX karyotype. This condition can occur as type I when isolated or as type II when associated with extragenital anomalies including kidney and skeletal abnormalities. The genetic basis of MRKH syndrome remains unexplained and several candidate genes have been proposed to play a role in its etiology, including HNF1B, LHX1, and WNT4. Here, we conducted a microarray analysis of 13 women affected by MRKH syndrome, resulting in the identification of chromosomal changes, including the deletion at 17q12, which contains both HNF1B and LHX1. We focused on HNF1B for further investigation due to its known association with, but unknown etiological role in, MRKH syndrome. We ablated Hnf1b specifically in the epithelium of the Müllerian ducts in mice and found that this caused hypoplastic development of the uterus, as well as kidney anomalies, closely mirroring the MRKH type II phenotype. Using single-cell RNA sequencing of uterine tissue in the Hnf1b-ablated embryos, we analyzed the molecules and pathways downstream of Hnf1b, revealing a dysregulation of processes associated with cell proliferation, migration, and differentiation. Thus, we establish that loss of Hnf1b function leads to an MRKH phenotype, and generate the first mouse model of MRKH syndrome type II. Our results support the investigation of HNF1B in clinical genetic settings of MRKH syndrome, and shed new light on the molecular mechanisms underlying this poorly understood condition in women's reproductive health.

Emery-Dreifuss muscular dystrophy (EDMD) is a genetically and clinically variable disorder. Previous attempts to use gene expression changes to find its pathomechanism were unavailing, so we engaged a functional pathway analysis. RNA-Seq was performed on cells from 10 patients diagnosed with an EDMD spectrum disease with different mutations in 7 genes. Upon comparing to controls, the pathway analysis revealed that multiple genes involved in fibrosis, metabolism, myogenic signaling, and splicing were affected in all patients. Splice variant analysis revealed alterations of muscle-specific variants for several important muscle genes. Deeper analysis of metabolic pathways revealed a reduction in glycolytic and oxidative metabolism and reduced numbers of mitochondria across a larger set of 14 EDMD spectrum patients and 7 controls. Intriguingly, the gene expression signatures segregated the patients into three subgroups whose distinctions could potentially relate to differences in clinical presentation. Finally, differential expression analysis of miRNAs changing in the patients similarly highlighted fibrosis, metabolism, and myogenic signaling pathways. This pathway approach revealed a transcriptome profile that can both be used as a template for establishing a biomarker panel for EDMD and direct further investigation into its pathomechanism. Furthermore, the segregation of specific gene changes into distinct groups that appear to correlate with clinical presentation may template development of prognostic biomarkers, though this will first require their testing in a wider set of patients with more clinical information.

Multiple sclerosis (MS) is a complex dysimmune disorder of the central nervous system. Genome-wide association studies (GWAS) have identified 233 genetic variations associated with MS at the genome-wide significant level. Epigenetic studies have pinpointed differentially methylated CpG sites in MS patients. However, the interplay between genetic risk factors and epigenetic regulation remains elusive. Here, we employed a network model to integrate GWAS summary statistics of 14 802 MS cases and 26 703 controls with DNA methylation profiles from 140 MS cases and 139 controls and the human interactome. We identified differentially methylated genes by aggregating additive effects of differentially methylated CpG sites within promoter regions. We reconstructed a gene regulatory network (GRN) using literature-curated transcription factor knowledge. Colocalization of the MS GWAS and methylation quantitative trait loci (mQTL) was performed to assess the GRN. The resultant MS-associated GRN highlighted several single nucleotide polymorphisms (SNPs) with GWAS-mQTL colocalization: rs6032663, rs6065926, and rs2024568 of CD40 locus, rs9913597 of STAT3 locus, and rs887864 and rs741175 of CIITA locus. Moreover, synergistic mQTL and eQTL signals were identified in CD40, suggesting gene expression alteration was likely induced by epigenetic changes. Web-based Cell-type Specific Enrichment Analysis of Genes (WebCSEA) indicated that the GRN was enriched in T follicular helper cells (p-value = 0.0016). Drug target enrichment analysis of annotations from the Therapeutic Target Database (TTD) revealed the GRN was also enriched with drug target genes (p-value = 3.89 × 10-4), revealing repurposable candidates for MS treatment. These candidates included vorinostat (HDAC1 inhibitor) and sivelestat (ELANE inhibitor), which warrant further investigation.

Autosomal dominant adult-onset neuronal ceroid lipofuscinosis (ANCL) is a rare neurodegenerative disorder characterised by progressive dementia and premature death. Four ANCL-causing mutations have been identified, all mapping to the DNAJC5 gene that encodes cysteine string protein α (CSPα). Here, using Caenorhabditis elegans, we describe an animal model of ANCL in which disease-causing mutations are introduced into their endogenous chromosomal locus, thereby mirroring the human genetic disorder. This was achieved through CRISPR/Cas9 mediated gene editing of dnj-14, the C. elegans ortholog of DNAJC5. The resultant homozygous ANCL mutant worms exhibited reduced lifespans and severely impaired chemotaxis, similar to isogenic dnj-14 null mutants. Importantly, these phenotypes were also seen in balanced heterozygotes carrying one wild-type and one ANCL mutant dnj-14 allele, mimicking the heterozygosity of ANCL patients. We observed a more severe chemotaxis phenotype in heterozygous ANCL mutant worms compared to haploinsufficient worms lacking one copy of CSP, consistent with a dominant-negative mechanism of action. Additionally, we provide evidence of CSP haploinsufficiency in longevity, as heterozygous null mutants exhibited significantly shorter lifespan than wild-type controls. The chemotaxis phenotype of dnj-14 null mutants was fully rescued by transgenic human CSPα, confirming the translational relevance of the worm model. Finally, a focused compound screen revealed that the anti-epileptic drug ethosuximide could restore chemotaxis in dnj-14 ANCL mutants to wild-type levels. This suggests that ethosuximide may have therapeutic potential for ANCL and demonstrates the utility of this C. elegans model for future larger-scale drug screening.

BACKGROUND: Variants in KCNH2, encoding the hERG channel which is responsible for the rapid component of the cardiac delayed rectifier K+ current (IKr), are causal to Long QT Syndrome type 2 (LQTS2). We identified eight index patients with a new variant of unknown significance (VUS), KCNH2:c.2717C > T:p.(Ser906Leu). We aimed to elucidate the biophysiological effect of this variant, to enable reclassification and consequent clinical decision-making. METHODS: A genotype-phenotype overview of the patients and relatives was created. The biophysiological effects were assessed independently by manual and automated calibrated patch-clamp. HEK293a cells expressing: (I) wild type (WT) KCNH2, (II) KCNH2-p.S906L alone (homozygous, Hm) or (III) KCNH2-p.S906L in combination with WT (1:1) (heterozygous, Hz), were used for manual patching. Automated patch-clamp measured the variants function against known benign and pathogenic variants, using Flp-In T-rex HEK293 KCNH2-variant cell lines. RESULTS: Incomplete penetrance of LQTS2 in KCNH2:p.(Ser906Leu) carriers was observed. In addition, some patients were heterozygous for other VUSs in CACNA1C, PKP2, RYR2, or AKAP9. The phenotype of carriers of KCNH2:p.(Ser906Leu) ranged from asymptomatic to life-threatening arrhythmic events. Manual patch-clamp showed a reduced current density by 69.8%, and 60.4% in KCNH2-p.S906L-Hm and KCNH2-p.S906L-Hz, respectively. The time constant of activation was significantly increased with 80.1% in KCNH2-p.S906L-Hm compared to KCNH2-WT. Assessment of KCNH2-p.S906L-Hz, by calibrated automatic patch-clamp assay showed a reduction in current density by 35.6%. CONCLUSION: The reduced current density in the KCNH2-p.S906L-Hz indicates a moderate loss-of-function. Combined with the reduced penetrance and variable phenotype, we conclude that KCNH2:p.(Ser906Leu) is a low penetrant likely pathogenic variant for LQTS2.

Enhancers are pivotal for regulating gene transcription that occurs at promoters. Identification of the interacting enhancer-promoter pairs and understanding the mechanisms behind how they interact and how enhancers modulate transcription can provide fundamental insight into gene regulatory networks. Recently, advances in high-throughput methods in three major areas-chromosome conformation capture assay, such as Hi-C to study basic chromatin architecture, ectopic reporter experiments such as self-transcribing active regulatory region sequencing (STARR-seq) to quantify promoter and enhancer activity, and endogenous perturbations such as clustered regularly interspaced short palindromic repeat interference (CRISPRi) to identify enhancer-promoter compatibility-have further our knowledge about transcription. In this review, we will discuss the major method developments and key findings from these assays.

Fourteen years after the first genome-wide association study (GWAS) of lung cancer was published, approximately 45 genomic loci have now been significantly associated with lung cancer risk. While functional characterization was performed for several of these loci, a comprehensive summary of the current molecular understanding of lung cancer risk has been lacking. Further, many novel computational and experimental tools now became available to accelerate the functional assessment of disease-associated variants, moving beyond locus-by-locus approaches. In this review, we first highlight the heterogeneity of lung cancer GWAS findings across histological subtypes, ancestries and smoking status, which poses unique challenges to follow-up studies. We then summarize the published lung cancer post-GWAS studies for each risk-associated locus to assess the current understanding of biological mechanisms beyond the initial statistical association. We further summarize strategies for GWAS functional follow-up studies considering cutting-edge functional genomics tools and providing a catalog of available resources relevant to lung cancer. Overall, we aim to highlight the importance of integrating computational and experimental approaches to draw biological insights from the lung cancer GWAS results beyond association.

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

Genome-wide association studies (GWAS) have found the majority of disease-associated variants to be non-coding. Major efforts into the charting of the non-coding regulatory landscapes have allowed for the development of tools and methods which aim to aid in the identification of causal variants and their mechanism of action. In this review, we give an overview of current tools and methods for the analysis of non-coding GWAS variants in disease. We provide a workflow that allows for the accumulation of in silico evidence to generate novel hypotheses on mechanisms underlying disease and prioritize targets for follow-up study using non-coding GWAS variants. Lastly, we discuss the need for comprehensive benchmarks and novel tools for the analysis of non-coding variants.

Non-coding genetic variants outside of protein-coding genome regions play an important role in genetic and epigenetic regulation. It has become increasingly important to understand their roles, as non-coding variants often make up the majority of top findings of genome-wide association studies (GWAS). In addition, the growing popularity of disease-specific whole-genome sequencing (WGS) efforts expands the library of and offers unique opportunities for investigating both common and rare non-coding variants, which are typically not detected in more limited GWAS approaches. However, the sheer size and breadth of WGS data introduce additional challenges to predicting functional impacts in terms of data analysis and interpretation. This review focuses on the recent approaches developed for efficient, at-scale annotation and prioritization of non-coding variants uncovered in WGS analyses. In particular, we review the latest scalable annotation tools, databases and functional genomic resources for interpreting the variant findings from WGS based on both experimental data and in silico predictive annotations. We also review machine learning-based predictive models for variant scoring and prioritization. We conclude with a discussion of future research directions which will enhance the data and tools necessary for the effective functional analyses of variants identified by WGS to improve our understanding of disease etiology.

Obesity and type 2 diabetes are major contributors to the growing prevalence of non-alcoholic fatty liver disease (NAFLD), a chronic liver condition characterized by the accumulation of fat in individuals without a significant amount of alcohol intake. The NAFLD spectrum ranges from simple steatosis (early stages, known as NAFL) to non-alcoholic steatohepatitis, which can progress to fibrosis and cirrhosis or hepatocellular carcinoma. Obesity, type 2 diabetes and NAFLD are strongly associated with insulin resistance. In the liver, insulin resistance increases hepatic glucose output, lipogenesis and very-low-density lipoprotein secretion, leading to a combination of hyperglycemia and hypertriglyceridemia. Aberrant gene expression is a hallmark of insulin resistance. Non-coding RNAs (ncRNAs) have emerged as prominent regulators of gene expression that operate at the transcriptional, post-transcriptional and post-translational levels. In the last couple of decades, a wealth of studies have provided evidence that most processes of liver metabolism are orchestrated by ncRNAs. This review focuses on the role of microRNAs, long non-coding RNAs and circular RNAs as coordinators of hepatic function, as well as the current understanding on how their dysregulation contributes to abnormal metabolism and pathophysiology in animal models of insulin resistance and NAFLD. Moreover, ncRNAs are emerging as useful biomarkers that may be able to discriminate between the different stages of NAFLD. The potential of ncRNAs as therapeutic drugs for NAFLD treatment and as biomarkers is discussed.

Linkage disequilibrium and the incomplete regulatory annotation of the noncoding genome complicates the identification of functional noncoding genetic variants and their causal association with disease. Current computational methods for variant prioritization have limited predictive value, necessitating the application of highly parallelized experimental assays to efficiently identify functional noncoding variation. Here, we summarize two distinct approaches, massively parallel reporter assays and CRISPR-based pooled screens and describe their flexible implementation to characterize human noncoding genetic variation at unprecedented scale. Each approach provides unique advantages and limitations, highlighting the importance of multimodal methodological integration. These multiplexed assays of variant effects are undoubtedly poised to play a key role in the experimental characterization of noncoding genetic risk, informing our understanding of the underlying mechanisms of disease-associated loci and the development of more robust predictive classification algorithms.

Advances in genomics have led to the identification of many risk loci with hundreds of genes and thousands of DNA variants associated with neuropsychiatric disorders. A significant barrier to understanding the genetic underpinnings of complex diseases is the lack of functional characterization of risk genes and variants in biological systems relevant to human health and connecting disease-associated variants to pathological phenotypes. Characterizing gene and DNA variant functions requires genetic perturbations followed by molecular and cellular assays of neurobiological phenotypes. However, generating null or mutant alleles is low throughput, making it impossible to characterize disease-associated variants in large quantities efficiently. CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) screens can be leveraged to dissect the biological consequences of the tested genes and variants in their native context. Nevertheless, testing non-coding variants associated with complex diseases remains non-trivial. In this review, we first discuss the current challenges of interpreting the function of the non-coding genome and approaches to prioritizing disease-associated variants in the context of the 3D epigenome. Second, we provide a brief overview of high-throughput CRISPRi and CRISPRa screening strategies applicable for characterizing non-coding sequences in appropriate biological systems. Lastly, we discuss the promising prospects of using CRISPR-based technologies to dissect DNA sequences associated with neuropsychiatric diseases.

Advances and reduction of costs in various sequencing technologies allow for a closer look at variations present in the non-coding regions of the human genome. Correlating non-coding variants with large-scale multi-omic data holds the promise not only of a better understanding of likely causal connections between non-coding DNA and expression of traits but also identifying potential disease-modifying medicines. Genome-phenome association studies have created large datasets of DNA variants that are associated with multiple traits or diseases, such as Alzheimer's disease; yet, the functional consequences of variants, in particular of non-coding variants, remain largely unknown. Recent advances in functional genomics and computational approaches have led to the identification of potential roles of DNA variants, such as various quantitative trait locus (xQTL) techniques. Multi-omic assays and analytic approaches toward xQTL have identified links between genetic loci and human transcriptomic, epigenomic, proteomic and metabolomic data. In this review, we first discuss the recent development of xQTL from multi-omic findings. We then highlight multimodal analysis of xQTL and genetic data for identification of risk genes and drug targets using Alzheimer's disease as an example. We finally discuss challenges and future research directions (e.g. artificial intelligence) for annotation of non-coding variants in complex diseases.

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

The exponential accumulation of DNA sequencing data has opened new avenues for discovering the causative roles of single-nucleotide polymorphisms (SNPs) in neurological diseases. The opportunities emerging from this are staggering, yet only as good as our abilities to glean insights from this surplus of information. Whereas computational biology continues to improve with respect to predictions and molecular modeling, the differences between in silico and in vivo analysis remain substantial. Invertebrate in vivo model systems represent technically advanced, experimentally mature, high-throughput, efficient and cost-effective resources for investigating a disease. With a decades-long track record of enabling investigators to discern function from DNA, fly (Drosophila) and worm (Caenorhabditis elegans) models have never been better poised to serve as living engines of discovery. Both of these animals have already proven useful in the classification of genetic variants as either pathogenic or benign across a range of neurodevelopmental and neurodegenerative disorders-including autism spectrum disorders, ciliopathies, amyotrophic lateral sclerosis, Alzheimer's and Parkinson's disease. Pathogenic SNPs typically display distinctive phenotypes in functional assays when compared with null alleles and frequently lead to protein products with gain-of-function or partial loss-of-function properties that contribute to neurological disease pathogenesis. The utility of invertebrates is logically limited by overt differences in anatomical and physiological characteristics, and also the evolutionary distance in genome structure. Nevertheless, functional annotation of disease-SNPs using invertebrate models can expedite the process of assigning cellular and organismal consequences to mutations, ascertain insights into mechanisms of action, and accelerate therapeutic target discovery and drug development for neurological conditions.

Neurodegenerative diseases such as Alzheimer's disease (AD) are characterized by the progressive loss of neurons in the brain and the spinal cord. The pathophysiology of AD is multifactorial with heterogeneous molecular manifestations. The lack of efficacious therapies for AD reinforces the importance of exploring in depth multifaceted disease mechanisms. Recent progresses on AD have generated a large amount of RNA-sequencing data at both bulk and single cell levels and revealed thousands of genes with expression changes in AD. However, the upstream regulators of such gene expression changes are largely unknown. Non-coding RNAs (ncRNAs) represent the majority of the human transcriptome, and regulatory ncRNAs have been found to play an important role in regulating gene expression. A single miRNA usually targets a number of mRNAs and thus such ncRNAs are particular important for understanding disease mechanisms and developing novel therapeutics. This review aims to summarize the recent findings on the roles of ncRNAs in AD from ncRNA-omics studies with a focus on ncRNA signatures, interactions between ncRNAs and mRNAs, and ncRNA-regulated pathways in AD. We also review the potential of specific ncRNAs to serve as biomarkers and therapeutic targets for AD. In the end, we point out future directions for studying ncRNAs in AD.

How ancestry-associated genetic variance affects disparities in the risk for polygenic diseases and influences the identification of disease-associated genes warrant a deeper understanding. We hypothesized that the discovery of genes associated with polygenic diseases may be limited by overreliance on single-nucleotide polymorphism (SNP)-based genomic investigation, since most significant variants identified in genome-wide SNP association studies map to introns and intergenic regions of the genome. To overcome such potential limitation, we developed a gene-constrained and function-based analytical method centered on high-risk variants (hrV) that encode frameshifts, stopgains, or splice site disruption. We analyzed the total number of hrV per gene in populations of different ancestry, representing a total of 185 934 subjects. Using this analysis, we developed a quantitative index of hrV (hrVI) across 20 428 genes within each population. We then applied hrVI analysis to the discovery of genes associated with type 2 diabetes mellitus (T2DM), a polygenic disease with ancestry-related disparity. HrVI profiling and gene-to-gene comparisons of ancestry-specific hrV between the case (20 781 subjects) and control (24 440 subjects) populations in the T2DM national repository identified 57 genes associated with T2DM, 40 of which were discoverable only by ancestry-specific analysis. These results illustrate how function-based and ancestry-specific analysis of genetic variations can accelerate the identification of genes associated with polygenic diseases. Besides T2DM, such analysis may facilitate our understanding of the genetic basis for other polygenic diseases that are also greatly influenced by environmental and behavioral factors, such as obesity, hypertension, and Alzheimer's disease.

Spinal muscular atrophy (SMA) is a fatal neuromuscular disease caused by homozygous deletions or mutations of the SMN1 gene. SMN2 is a paralogous gene of SMN1 and a modifying gene of SMA. A better understanding of how SMN2 exon 7 splicing is regulated helps discover new therapeutic targets for SMA therapy. Based on an antisense walk method to map exonic and intronic splicing silencers (ESSs and ISSs) in SMN2 exon 7 and the proximal regions of its flanking introns, we identified one ISS (ISS6-KH) at upstream of the branch point site in intron 6. By using mutagenesis-coupled RT-PCR with SMN1/2 minigenes, immunochromatography, over-expression, and siRNA-knockdown, we found this ISS consists of a bi-partite hnRNP A1 binding cis-element and a poly-U sequence located between the proximal hnRNP A1 binding site (UAGCUA) and the branch site. Both HuR and hnRNP C1 proteins promote exon 7 skipping through the poly-U stretch. Mutations or deletions of these motifs lead to efficient SMN2 exon 7 inclusion comparable to SMN1 gene. Furthermore, we identified an optimal antisense oligonucleotide (ASO) that binds the intron 6 ISS and causes striking exon 7 inclusion in the SMN2 gene in patient fibroblasts and SMA mouse model. Our findings demonstrate that this novel ISS plays an important role in SMN2 exon 7 skipping and highlight a new therapeutic target for SMA therapy.

SUPT16H encodes the large subunit of the FACT complex, which functions as a nucleosome organizer during transcription. We identified two individuals from unrelated families carrying de novo missense variants in SUPT16H. The probands exhibit global developmental delay, intellectual disability, epilepsy, facial dysmorphism, and brain structural abnormalities. We used Drosophila to characterize these two variants, p.T171I and p.G808R. Loss of the fly ortholog, dre4, causes lethality at an early developmental stage. RNAi-mediated knockdown of dre4 in either glia or neurons causes severely reduced eclosion and longevity. Tissue-specific knockdown of dre4 in the eye or wing leads to the loss of these tissues whereas overexpression of SUPT16H has no dominant effect. Moreover, expression of the reference SUPT16H significantly rescues the loss-of-function phenotypes in the nervous system as well as wing and eye. In contrast, expression of SUPT16H p.T171I or p.G808R rescue the phenotypes poorly, indicating that the variants are partial loss-of-function alleles. While previous studies argued that the developmental arrest caused by loss of dre4 is due to impaired ecdysone production in the prothoracic gland, our data show that dre4 is required for proper cell growth and survival in multiple tissues in a cell-autonomous manner. Altogether, our data indicate that the de novo loss-of-function variants in SUPT16H are indeed associated with developmental and neurological defects observed in the probands.

Haploinsufficiency of TGF-Beta Activated Kinase 1 (MAP3K7) Binding Protein 2 (TAB2) has been associated with congenital heart disease and more recently multiorgan structural abnormalities. Missense variant represents a major proportion of non-synonymous TAB2 variants reported in gnomAD (295/576) and Clinvar (16/73), most of which are Variants of uncertain significance (VUSs). However, interpretation of TAB2 missense variants remains challenging due to lack of functional assay. To address this issue, we established a cell-based luciferase assay that enables high-throughput screening of TAB2 variants to assess the functional consequence for predicting variant pathogenicity. Using this platform, we screened 47 TAB2 variants including 5 pathogenic controls and 1 benign control, and the results showed that the transcriptional activity of AP-1 but not NF-κB predicts the TAB2 variant pathogenicity. This assay provides accurate functional readout for both loss-of-function (LOF) and gain-of-function (GOF) variants, which are associated with distinct phenotypes. 22 out of 32 tested VUSs were reclassified. Genotype-phenotype association showed that most patients with partial LOF variants do not exhibit congenital heart disease but high frequency of developmental delay, hypotonia, and dysmorphic features, which suggests genetic testing for TAB2 is needed for a broader spectrum of patients with more diverse phenotypes. Molecular modeling with Npl4 zinc finger (NZF) domain variants revealed that the stability of the NZF domain in TAB2 protein is crucial for AP-1 activation. In conclusion, we developed a highly effective functional assay for TAB2 variant prediction and interpretation.

NRSF/REST plays a key role in neuronal homeostasis as a transcriptional repressor of neuronal genes. NRSF/REST relates to cognitive preservation and longevity of humans, but its specific functions in age-dependent and Alzheimer's disease (ad)-related memory deficits remain unclear. Here we show that conditional NRSF/REST knockout either in the dorsal telencephalon or specially in neurons induced an age-dependently diminished retrieval performance in spatial or fear conditioning memory tasks and altered hippocampal synaptic transmission and activity-dependent synaptic plasticity. The NRSF/REST deficient mice were also characterized by an increase of activated glial cells, complement C3 protein, and the transcription factor C/EBPβ in the cortex and hippocampus. Reduction of NRSF/REST by conditional depletion upregulated the activation of astrocytes in APP/PS1 mice, and increased the C3-positive glial cells, but did not alter the Aβ loads and memory retrieval performances of 6- and 12-month-old APP/PS1 mice. Simultaneously, overexpression of NRSF/REST improved cognitive abilities of aged WT, but not in ad mice. These findings demonstrated that NRSF/REST is essential for the preservation of memory performance and activity-dependent synaptic plasticity during aging and takes potential roles in the onset of age-related memory impairments. However, while altering the glial activation, NRSF/REST deficiency does not interfere with the Aβ deposits, and the electrophysiological and cognitive ad-like pathologies.

Maf-family basic motif leucine zipper protein NRL specifies rod photoreceptor cell fate during retinal development and, in concert with homeodomain protein CRX and other regulatory factors, controls the expression of most rod-expressed genes including the visual pigment gene Rhodospin (Rho). Transcriptional regulatory activity of NRL is modulated by post-translational modifications (PTMs), especially phosphorylation, and mutations at specific phosphosites can lead to retinal degeneration. During our studies to elucidate NRL-mediated transcriptional regulation, we identified Protein Kinase CK2 in NRL-enriched complexes bound to Rho promoter-enhancer regions and in NRL-enriched high molecular mass fractions from the bovine retina. The presence of CK2 in NRL complexes was confirmed by co-immunoprecipitation from developing and adult mouse retinal extracts. In vitro kinase assay and bioinformatic analysis indicated phosphorylation of NRL at Ser117 residue by CK2. Co-transfection of Csnk2a1 cDNA encoding murine CK2 with human NRL and CRX reduced the bovine Rho promoter-driven luciferase expression in HEK293 cells and mutagenesis of NRL-Ser117 residue to Ala restored the reporter gene activity. In concordance, over-expression of CK2 in the mouse retina in vivo by electroporation resulted in reduction of Rho promoter-driven DsRed reporter expression as well as the transcript level of many phototransduction genes. Thus, our studies demonstrate that CK2 can phosphorylate Ser117 of NRL. Modulation of NRL activity by CK2 suggests intricate inter-dependence of transcriptional and signaling pathways in maintaining rod homeostasis.

Chromatin regulators constitute a fundamental means of transcription regulation, which have been implicated in neurodevelopment and neurodevelopment disorders (NDDs). Supt16, one of candidate genes for NDDs, encodes the large subunit of facilitates chromatin transcription (FACT). However, the underlying mechanisms remain poorly understood. Here, Supt16+/- mice was generated, modeling the neurodevelopment disorder. Abnormal cognitive and social behavior were observed in the Supt16 +/- mice. Simultaneously, the number of neurocytes in the cerebral cortex and hippocampus is decreased, which might be resulted from the impairment of mouse neural stem cells (mNSCs) in the SVZ. Supt16 haploinsufficiency affects the proliferation and apoptosis of mNSCs. As the RNA-seq and ChIP-seq assays showed, Supt16 haploinsufficiency disrupts the stemness of mNSCs by inhibiting MAPK signal pathway. Thus, this study demonstrates a critical role of Supt16 gene in the proliferation and apoptosis of mNSCs and provids a novel insight in the pathogenesis of NDDs.

Myotonic dystrophy type 1 (DM1) is a multi-systemic disorder caused by expansion of CTG microsatellite repeats within DMPK. The most severe form, congenital myotonic dystrophy (CDM), has symptom onset at birth due to large intergenerational repeat expansions. Despite a common mutation, CDM individuals present with a distinct clinical phenotype and absence of common DM1 symptoms. Given the clinical divergence, it is unknown if the hallmark of DM1 pathology, dysregulation of alternative splicing (AS) due to sequestration of MBNL proteins within toxic CUG repeat RNAs, contributes to disease throughout pediatric development. To evaluate global transcriptomic dysregulation RNA-seq was performed on 36 CDM skeletal muscle biopsies ages 2 weeks to 16 years, including two longitudinal samples. 50 DM1 and adult/pediatric controls were also sequenced as comparative groups. Despite a large CTG expansion and shared age of onset, CDM individuals presented with a heterogenous, MBNL-dependent mis-splicing signature. Estimation of intracellular MBNL concentrations from splicing responses of select events correlated with total spliceopathy and revealed a distinct, triphasic pattern of AS dysregulation across pediatric development. CDM infants (< 2 years) possess severe mis-splicing which significantly improves in early childhood (2-8 years) independent of sex or CTG repeat load. Adolescent individuals (8-16 years) stratified into two populations with a full range of global splicing dysregulation. DMPK expression changes correlated with alterations in splicing severity during development. This study reveals the complex dynamics of the CDM muscle transcriptome and provides insights into new therapeutic strategies, timing of therapeutic intervention, and biomarker development.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. Its complex pathogenesis and phenotypic heterogeneity hinder therapeutic development and early diagnosis. Altered RNA metabolism is a recurrent pathophysiologic theme, including distinct microRNA (miRNA) profiles in ALS tissues. We profiled miRNAs in accessible biosamples, including skin fibroblasts and whole blood and compared them in age- and sex-matched healthy controls versus ALS participants with and without repeat expansions to chromosome 9 open reading frame 72 (C9orf72; C9-ALS and nonC9-ALS), the most frequent ALS mutation. We identified unique and shared profiles of differential miRNA (DmiRNA) levels in each C9-ALS and nonC9-ALS tissues versus controls. Fibroblast DmiRNAs were validated by qPCR and their target mRNAs by 5-bromouridine and 5-bromouridine-chase sequencing. We also performed pathway analysis to infer biological meaning, revealing anticipated, tissue-specific pathways and pathways previously linked to ALS, as well as novel pathways that could inform future research directions. Overall, we report a comprehensive study of a miRNA profile dataset from C9-ALS and nonC9-ALS participants across two accessible biosamples, providing evidence of dysregulated miRNAs in ALS and possible targets of interest. Distinct miRNA patterns in accessible tissues may also be leveraged to distinguish ALS participants from healthy controls for earlier diagnosis. Future directions may look at potential correlations of miRNA profiles with clinical parameters.

Parkinson's disease (PD) is the second most common neurodegenerative disorder, and aging and genetic and environmental exposure can contribute to its pathogenesis. DNA methylation has been suggested to play a pivotal role in neurodevelopment and neurodegenerative diseases. 5-hydroxymethylcytosine (5hmC) is generated through 5-methylcytosine (5mC) oxidization by ten-eleven translocation proteins and is particularly enriched in the brain. Although 5hmC has been linked to multiple neurological disorders, little is known about 5hmC alterations in the substantia nigra of patients with PD. To determine the specific alterations in DNA methylation and hydroxymethylation in PD brain samples, we examined the genome-wide profiles of 5mC and 5hmC in the substantia nigra of patients with PD and Alzheimer's disease (ad). We identified 4119 differentially hydroxymethylated regions (DhMRs) and no differentially methylated regions (DMRs) in the postmortem brains of patients with PD compared with those of controls. These DhMRs were PD-specific when compared with the results of AD. Gene ontology analysis revealed that several signaling pathways, such as neurogenesis and neuronal differentiation, were significantly enriched in PD DhMRs. KEGG enrichment analysis revealed substantial alterations in multiple signaling pathways, including phospholipase D (PLD), cAMP and Rap1. In addition, using a PD Drosophila model, we found that one of the 5hmC-modulated genes, PLD1, modulated α-synuclein toxicity. Our analysis suggested that 5hmC may act as an independent epigenetic marker and contribute to the pathogenesis of PD.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

Mutations in genes encoding subunits of the BAF (BRG1/BRM-associated factor) complex cause various neurodevelopmental diseases. However, the underlying pathophysiology remains largely unknown. Here, we analyzed the function of Brahma-related gene 1 (Brg1), a core ATPase of BAF complexes, in the developing cerebral cortex. Loss of Brg1 causes several morphological defects resembling human malformations of cortical developments (MCDs), including microcephaly, cortical dysplasia, cobblestone lissencephaly and periventricular heterotopia. We demonstrated that neural progenitor cell renewal, neuronal differentiation, neuronal migration, apoptotic cell death, pial basement membrane and apical junctional complexes, which are associated with MCD formation, were impaired after Brg1 deletion. Furthermore, transcriptome profiling indicated that a large number of genes were deregulated. The deregulated genes were closely related to MCD formation, and most of these genes were bound by Brg1. Cumulatively, our study indicates an essential role of Brg1 in cortical development and provides a new possible pathogenesis underlying Brg1-based BAF complex-related neurodevelopmental disorders.

Autosomal dominant optic atrophy (DOA) is the most common inherited optic neuropathy, characterized by the preferential loss of retinal ganglion cells (RGCs), resulting in optic nerve degeneration and progressive bilateral central vision loss. More than 60% of genetically confirmed patients with DOA carry variants in the nuclear OPA1 gene, which encodes for a ubiquitously expressed, mitochondrial GTPase protein. OPA1 has diverse functions within the mitochondrial network, facilitating inner membrane fusion and cristae modelling, regulating mitochondrial DNA maintenance and coordinating mitochondrial bioenergetics. There are currently no licensed disease-modifying therapies for DOA and the disease mechanisms driving RGC degeneration are poorly understood. Here, we describe the generation of isogenic, heterozygous OPA1 null induced pluripotent stem cell (iPSC) (OPA1+/-) through clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing of a control cell line, in conjunction with the generation of DOA patient-derived iPSC carrying OPA1 variants, namely, the c.2708_2711delTTAG variant (DOA iPSC), and previously reported missense variant iPSC line (c.1334G>A, DOA plus [DOA]+ iPSC) and CRISPR/Cas9 corrected controls. A two-dimensional (2D) differentiation protocol was used to study the effect of OPA1 variants on iPSC-RGC differentiation and mitochondrial function. OPA1+/-, DOA and DOA+ iPSC showed no differentiation deficit compared to control iPSC lines, exhibiting comparable expression of all relevant markers at each stage of differentiation. OPA1+/- and OPA1 variant iPSC-RGCs exhibited impaired mitochondrial homeostasis, with reduced bioenergetic output and compromised mitochondrial DNA maintenance. These data highlight mitochondrial deficits associated with OPA1 dysfunction in human iPSC-RGCs, and establish a platform to study disease mechanisms that contribute to RGC loss in DOA, as well as potential therapeutic interventions.

Author information: (1)Department of Neurology, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA. (2)Department of Psychiatry and Behavioral Sciences, Stanford University, Stanford, CA 94305, USA. (3)Department of Psychology, John Carroll University, University Heights, OH 44118, USA. (4)Department of Pediatrics, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA. (5)Rosamund Stone Zander Translational Neuroscience Center, Boston Children's Hospital, Harvard Medical School, Boston, MA 02115, USA. (6)Genomic Medicine Institute, Cleveland Clinic, Cleveland, OH 44195, USA. (7)Department of Pediatrics, Hematology, Oncology, Blood and Marrow Transplantation, Cleveland Clinic, Cleveland, OH 44195, USA.

Although autism is typically characterized by differences in language, social interaction and restrictive, repetitive behaviors, it is becoming more well known in the field that alterations in energy metabolism and mitochondrial function are comorbid disorders in autism. The synaptic cell adhesion molecule, neurexin-1 (NRXN1), has previously been implicated in autism, and here we show that in Drosophila melanogaster, the homologue of NRXN1, called Nrx-1, regulates energy metabolism and nutrient homeostasis. First, we show that Nrx-1-null flies exhibit decreased resistance to nutrient deprivation and heat stress compared to controls. Additionally, Nrx-1 mutants exhibit a significantly altered metabolic profile characterized by decreased lipid and carbohydrate stores. Nrx-1-null Drosophila also exhibit diminished levels of nicotinamide adenine dinucleotide (NAD+), an important coenzyme in major energy metabolism pathways. Moreover, loss of Nrx-1 resulted in striking abnormalities in mitochondrial morphology in the flight muscle of Nrx-1-null Drosophila and impaired flight ability in these flies. Further, following a mechanical shock Nrx-1-null flies exhibited seizure-like activity, a phenotype previously linked to defects in mitochondrial metabolism and a common symptom of patients with NRXN1 deletions. The current studies indicate a novel role for NRXN1 in the regulation of energy metabolism and uncover a clinically relevant seizure phenotype in Drosophila lacking Nrx-1.

Heterozygous variants in bromodomain and plant homeodomain containing transcription factor (BPTF) cause the neurodevelopmental disorder with dysmorphic facies and distal limb anomalies (NEDDFL) syndrome (MIM#617755) characterized by intellectual disability, speech delay and postnatal microcephaly. BPTF functions within nucleosome and remodeling factor (NURF), a complex comprising sucrose non-fermenting like (SNF2L), an Imitation SWItching (ISWI) chromatin remodeling protein encoded by the SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily a, member 1 (SMARCA1) gene. Surprisingly, ablation of Smarca1 resulted in mice with enlarged brains, a direct contrast to the phenotype of NEDDFL patients. To model the NEDDFL syndrome, we generated forebrain-specific Bptf knockout (Bptf conditional Knockout (cKO)) mice. The Bptf cKO mice were born in normal Mendelian ratios, survived to adulthood but were smaller in size with severe cortical hypoplasia. Prolonged progenitor cell cycle length and a high incidence of cell death reduced the neuronal output. Cortical lamination was also disrupted with reduced proportions of deep layer neurons, and neuronal maturation defects that impaired the acquisition of distinct cell fates (e.g. COUP-TF-interacting protein 2 (Ctip2)+ neurons). RNAseq and pathway analysis identified altered expression of fate-determining transcription factors and the biological pathways involved in neural development, apoptotic signaling and amino acid biosynthesis. Dysregulated genes were enriched for MYC Proto-Oncogene, BHLH Transcription Factor (Myc)-binding sites, a known BPTF transcriptional cofactor. We propose the Bptf cKO mouse as a valuable model for further study of the NEDDFL syndrome.

Frataxin deficiency in Friedreich's ataxia results from transcriptional downregulation of the FXN gene caused by expansion of the intronic trinucleotide guanine-adenine-adenine (GAA) repeats. We used multiple transcriptomic approaches to determine the molecular mechanism of transcription inhibition caused by long GAAs. We uncovered that transcription of FXN in patient cells is prematurely terminated upstream of the expanded repeats leading to the formation of a novel, truncated and stable RNA. This FXN early terminated transcript (FXN-ett) undergoes alternative, non-productive splicing and does not contribute to the synthesis of functional frataxin. The level the FXN-ett RNA directly correlates with the length of the longer of the two expanded GAA tracts. Targeting GAAs with antisense oligonucleotides or excision of the repeats eliminates the transcription impediment, diminishes expression of the aberrant FXN-ett, while increasing levels of FXN mRNA and frataxin. Non-productive transcription may represent a common phenomenon and attractive therapeutic target in diseases caused by repeat-mediated transcription aberrations.

Progressive dilation of the infrarenal aortic diameter is a consequence of the ageing process and is considered the main determinant of abdominal aortic aneurysm (AAA). We aimed to investigate the genetic and clinical determinants of abdominal aortic diameter (AAD). We conducted a meta-analysis of genome-wide association studies in 10 cohorts (n = 13 542) imputed to the 1000 Genome Project reference panel including 12 815 subjects in the discovery phase and 727 subjects [Partners Biobank cohort 1 (PBIO)] as replication. Maximum anterior-posterior diameter of the infrarenal aorta was used as AAD. We also included exome array data (n = 14 480) from seven epidemiologic studies. Single-variant and gene-based associations were done using SeqMeta package. A Mendelian randomization analysis was applied to investigate the causal effect of a number of clinical risk factors on AAD. In genome-wide association study (GWAS) on AAD, rs74448815 in the intronic region of LDLRAD4 reached genome-wide significance (beta = -0.02, SE = 0.004, P-value = 2.10 × 10-8). The association replicated in the PBIO1 cohort (P-value = 8.19 × 10-4). In exome-array single-variant analysis (P-value threshold = 9 × 10-7), the lowest P-value was found for rs239259 located in SLC22A20 (beta = 0.007, P-value = 1.2 × 10-5). In the gene-based analysis (P-value threshold = 1.85 × 10-6), PCSK5 showed an association with AAD (P-value = 8.03 × 10-7). Furthermore, in Mendelian randomization analyses, we found evidence for genetic association of pulse pressure (beta = -0.003, P-value = 0.02), triglycerides (beta = -0.16, P-value = 0.008) and height (beta = 0.03, P-value < 0.0001), known risk factors for AAA, consistent with a causal association with AAD. Our findings point to new biology as well as highlighting gene regions in mechanisms that have previously been implicated in the genetics of other vascular diseases.

Metabolic alterations shared between the nervous system and skin fibroblasts have emerged in amyotrophic lateral sclerosis (ALS). Recently, we found that a subgroup of sporadic ALS (sALS) fibroblasts (sALS1) is characterized by metabolic profiles distinct from other sALS cases (sALS2) and controls, suggesting that metabolic therapies could be effective in sALS. The metabolic modulators nicotinamide riboside and pterostilbene (EH301) are under clinical development for the treatment of ALS. Here, we studied the transcriptome and metabolome of sALS cells to understand the molecular bases of sALS metabotypes and the impact of EH301. Metabolomics and transcriptomics were investigated at baseline and after EH301 treatment. Moreover, weighted gene coexpression network analysis (WGCNA) was used to investigate the association of the metabolic and clinical features. We found that the sALS1 transcriptome is distinct from sALS2 and that EH301 modifies gene expression differently in sALS1, sALS2 and the controls. Furthermore, EH301 had strong protective effects against metabolic stress, an effect linked to the antiinflammatory and antioxidant pathways. WGCNA revealed that the ALS functional rating scale and metabotypes are associated with gene modules enriched for the cell cycle, immunity, autophagy and metabolic genes, which are modified by EH301. The meta-analysis of publicly available transcriptomic data from induced motor neurons by Answer ALS confirmed the functional associations of genes correlated with disease traits. A subset of genes differentially expressed in sALS fibroblasts was used in a machine learning model to predict disease progression. In conclusion, multiomic analyses highlighted the differential metabolic and transcriptomic profiles in patient-derived fibroblast sALS, which translate into differential responses to the investigational drug EH301.

Recent research on familial dysautonomia (FD) has focused on the development of therapeutics that facilitate the production of the correctly spliced, exon 20-containing, transcript in cells and individuals bearing the splice-altering, FD-causing mutation in the elongator acetyltransferase complex subunit I (ELP1) gene. We report here the ability of carnosol, a diterpene present in plant species of the Lamiaceae family, including rosemary, to enhance the cellular presence of the correctly spliced ELP1 transcript in FD patient-derived fibroblasts by upregulating transcription of the ELP1 gene and correcting the aberrant splicing of the ELP1 transcript. Carnosol treatment also elevates the level of the RNA binding motif protein 24 (RBM24) and RNA binding motif protein 38 (RBM38) proteins, two multifunctional RNA-binding proteins. Transfection-mediated expression of either of these RNA binding motif (RBMs) facilitates the inclusion of exon 20 sequence into the transcript generated from a minigene-bearing ELP1 genomic sequence containing the FD-causing mutation. Suppression of the carnosol-mediated induction of either of these RBMs, using targeting siRNAs, limited the carnosol-mediated inclusion of the ELP1 exon 20 sequence. Carnosol treatment of FD patient peripheral blood mononuclear cells facilitates the inclusion of exon 20 into the ELP1 transcript. The increased levels of the ELP1 and RBM38 transcripts and the alternative splicing of the sirtuin 2 (SIRT2) transcript, a sentinel for exon 20 inclusion in the FD-derived ELP1 transcript, are observed in RNA isolated from whole blood of healthy adults following the ingestion of carnosol-containing rosemary extract. These findings and the excellent safety profile of rosemary together justify an expedited clinical study of the impact of carnosol on the FD patient population.

The correct expression of folded, functional rhodopsin (Rho) is critical for visual perception. However, this seven-transmembrane helical G protein-coupled receptor is prone to mutations with pathological consequences of retinal degeneration in retinitis pigmentosa (RP) due to Rho misfolding. Pharmacological chaperones that stabilize the inherited Rho variants by assisting their folding and membrane targeting could slow the progression of RP. In this study, we employed virtual screening of synthetic compounds with a natural product scaffold in conjunction with in vitro and in vivo evaluations to discover a novel chromenone-containing small molecule with favorable pharmacological properties that stabilize rod opsin. This compound reversibly binds to unliganded bovine rod opsin with an EC50 value comparable to the 9-cis-retinal chromophore analog and partially rescued membrane trafficking of multiple RP-related rod opsin variants in vitro. Importantly, this novel ligand of rod opsin was effective in vivo in murine models, protecting photoreceptors from deterioration caused by either bright light or genetic insult. Together, our current study suggests potential broad therapeutic implications of the new chromenone-containing non-retinoid small molecule against retinal diseases associated with photoreceptor degeneration.