The shaker rat carries a naturally occurring mutation leading to progressive ataxia characterized by Purkinje cell (PC) loss. We previously reported on fine-mapping the shaker locus to the long arm of the rat X chromosome. In this work, we sought to identify the mutated gene underlying the shaker phenotype and confirm its identity by functional complementation. We fine-mapped the candidate region and analyzed cerebellar transcriptomes, identifying a XM_217630.9 (Slc9a6):c.[191_195delinsA] variant in the Slc9a6 gene that segregated with disease. We generated an adeno-associated virus (AAV) targeting Slc9a6 expression to PCs using the mouse L7-6 (L7) promoter. We administered the AAV prior to the onset of PC degeneration through intracerebroventricular injection and found that it reduced the shaker motor, molecular, and cellular phenotypes. Therefore, Slc9a6 is mutated in shaker and AAV-based gene therapy may be a viable therapeutic strategy for Christianson syndrome, also caused by Slc9a6 mutation.

Autism spectrum disorder (ASD) affects 1 in 44 children. Chromatin regulatory proteins are overrepresented among genes that contain high risk variants in ASD. Disruption of the chromatin environment leads to widespread dysregulation of gene expression, which is traditionally thought of as a mechanism of disease pathogenesis associated with ASD. Alternatively, alterations in chromatin dynamics could also lead to dysregulation of alternative splicing, which is understudied as a mechanism of ASD pathogenesis. The anticonvulsant valproic acid (VPA) is a well-known environmental risk factor for ASD that acts as a class I histone deacetylase (HDAC) inhibitor. However, the precise molecular mechanisms underlying defects in human neuronal development associated with exposure to VPA are understudied. To dissect how VPA exposure and subsequent chromatin hyperacetylation influence molecular signatures involved in ASD pathogenesis, we conducted RNA sequencing (RNA-seq) in human cortical neurons that were treated with VPA. We observed that differentially expressed genes (DEGs) were enriched for mRNA splicing, mRNA processing, histone modification, and metabolism related gene sets. Furthermore, we observed widespread increase in the number and the type of alternative splicing events. Analysis of differential transcript usage (DTU) showed that exposure to VPA induces extensive alterations in transcript isoform usage across neurodevelopmentally important genes. Finally, we find that DEGs and genes that display DTU overlap with known ASD-risk genes. Together, these findings suggest that, in addition to differential gene expression, changes in alternative splicing correlated with alterations in the chromatin environment could act as an additional mechanism of disease in ASD.

Dominant mutations in ubiquitously expressed mitofusin 2 gene (MFN2) cause Charcot-Marie-Tooth type 2A (CMT2A; OMIM 609260), an inherited sensory-motor neuropathy that affects peripheral nerve axons. Mitofusin 2 protein has been found to take part in mitochondrial fusion, mitochondria-endoplasmic reticulum tethering, mitochondrial trafficking along axons, mitochondrial quality control and various types of cancer, in which MFN2 has been indicated as a tumor suppressor gene. Discordant data on the mitochondrial altered phenotypes in patient-derived fibroblasts harboring MFN2 mutations and in animal models have been reported. We addressed some of these issues by focusing on mitochondria behavior during autophagy and mitophagy in fibroblasts derived from a CMT2AMFN2 patient with an MFN2650G > T/C217F mutation in the GTPase domain. This study investigated mitochondrial dynamics, respiratory capacity and autophagy/mitophagy, to tackle the multifaceted MFN2 contribution to CMT2A pathogenesis. We found that MFN2 mutated fibroblasts showed impairment of mitochondrial morphology, bioenergetics capacity, and impairment of the early stages of autophagy, but not mitophagy. Unexpectedly, transcriptomic analysis of mutated fibroblasts highlighted marked differentially expressed pathways related to cell population proliferation and extracellular matrix organization. We consistently found the activation of mTORC2/AKT signaling and accelerated proliferation in the CMT2AMFN2 fibroblasts. In conclusion, our evidence indicates that MFN2 mutation can positively drive cell proliferation in CMT2AMFN2 fibroblasts.

Mutations in LMNA, the gene encoding A-type lamins, cause laminopathies-diseases of striated muscle and other tissues. The aetiology of laminopathies has been attributed to perturbation of chromatin organization or structural weakening of the nuclear envelope (NE) such that the nucleus becomes more prone to mechanical damage. The latter model requires a conduit for force transmission to the nucleus. NE-associated Linker of Nucleoskeleton and Cytoskeleton (LINC) complexes are one such pathway. Using clustered regularly interspaced short palindromic repeats to disrupt the Nesprin-1 KASH (Klarsicht, ANC-1, Syne Homology) domain, we identified this LINC complex protein as the predominant NE anchor for microtubule cytoskeleton components, including nucleation activities and motor complexes, in mouse cardiomyocytes. Loss of Nesprin-1 LINC complexes resulted in loss of microtubule cytoskeleton proteins at the nucleus and changes in nuclear morphology and positioning in striated muscle cells, but with no overt physiological defects. Disrupting the KASH domain of Nesprin-1 suppresses Lmna-linked cardiac pathology, likely by reducing microtubule cytoskeleton activities at the nucleus. Nesprin-1 LINC complexes thus represent a potential therapeutic target for striated muscle laminopathies.

As the auditory and balance receptor cells in the inner ear, hair cells are responsible for converting mechanical stimuli into electrical signals, a process referred to as mechano-electrical transduction (MET). Hair cell development and function are tightly regulated, and hair cell deficits are the main reasons for hearing loss and balance disorders. TMCC2 is an endoplasmic reticulum (ER)-residing transmembrane protein whose physiological function largely remains unknown. In the present work, we show that Tmcc2 is specifically expressed in the auditory hair cells of mouse inner ear. Tmcc2 knockout mice were then established to investigate its physiological role in hearing. Auditory brainstem responses (ABR) measurements show that Tmcc2 knockout mice suffer from congenital hearing loss. Further investigations reveal progressive auditory hair cell loss in Tmcc2 knockout mice. The general morphology and function of ER is unaffected in Tmcc2 knockout hair cells. However, increased ER stress was observed in Tmcc2 knockout mice and knockdown cells, suggesting that loss of TMCC2 leads to auditory hair cell death through elevated ER stress.

Type Ib glycogen storage disease (GSD-Ib) is caused by a deficiency in the glucose-6-phosphate (G6P) transporter (G6PT) that translocates G6P from the cytoplasm into the endoplasmic reticulum lumen, where the intraluminal G6P is hydrolyzed to glucose by glucose-6-phosphatase-α (G6Pase-α). Clinically, GSD-Ib patients manifest a metabolic phenotype of impaired blood glucose homeostasis and a long-term risk of hepatocellular adenoma/carcinoma (HCA/HCC). Studies have shown that autophagy deficiency contributes to hepatocarcinogenesis. In this study, we show that G6PT deficiency leads to impaired hepatic autophagy evident from attenuated expression of many components of the autophagy network, decreased autophagosome formation and reduced autophagy flux. The G6PT-deficient liver displayed impaired sirtuin 1 (SIRT1) and AMP-activated protein kinase (AMPK) signaling, along with reduced expression of SIRT1, forkhead boxO3a (FoxO3a), liver kinase B-1 (LKB1) and the active p-AMPK. Importantly, we show that overexpression of either SIRT1 or LKB1 in G6PT-deficient liver restored autophagy and SIRT1/FoxO3a and LKB1/AMPK signaling. The hepatosteatosis in G6PT-deficient liver decreased SIRT1 expression. LKB1 overexpression reduced hepatic triglyceride levels, providing a potential link between LKB1/AMPK signaling upregulation and the increase in SIRT1 expression. In conclusion, downregulation of SIRT1/FoxO3a and LKB1/AMPK signaling underlies impaired hepatic autophagy which may contribute to HCA/HCC development in GSD-Ib. Understanding this mechanism may guide future therapies.

Shorter gestational age (GA) is a risk factor of developmental delay. GA is usually estimated clinically from last menstrual period and ultrasound. DNA methylation (DNAm) estimates GA using sets of cytosine-guanine-sites coupled with a clock algorithm. Therefore, DNAm estimated GA may better reflect biological maturation. A DNAm GA greater than clinical GA, known as gestational age acceleration (GAA), may indicate epigenetic maturity and holds potential as an early biomarker for developmental delay risk. We used data from the Upstate KIDS Study to examine associations of DNAm GA and developmental delay within the first 3 years based on the Ages & Stages Questionnaire® (n = 1010). We estimated DNAm GA using two clocks specific to the Illumina Methylation EPIC 850 K the Haftorn clock and one developed from the Effects of Aspirin in Gestation and Reproduction study, in which women were followed to detect pregnancy at the earliest time possible. Among singletons, each week increase in DNAm GA was protective for overall delay (OR:0.74; 95%CI:0.61-0.90) and delay in all domains except for problem solving skills. Among twins, we observed similar point estimates but lower precision. Results were similar for clinical GA. GAA was largely not associated with developmental delays. In summary, either DNAm GA or clinical GA at birth, but not epigenetic maturity (i.e. GAA), was associated with decreased odds of developmental delay in early childhood. Our study does not support using DNAm GA or GAA as separate risk factors for future risk of developmental delay within the first 3 years of age.

DNA methylation plays a critical function in establishing and maintaining cell identity in brain. Disruption of DNA methylation-related processes leads to diverse neurological disorders. However, the role of DNA methylation characteristics in neuronal diversity remains underexplored. Here, we report detailed context-specific DNA methylation maps for GABAergic, glutamatergic (Glu) and Purkinje neurons, together with matched transcriptome profiles. Genome-wide mCH levels are distinguishable, while the mCG levels are similar among the three cell types. Substantial CG-differentially methylated regions (DMRs) are also seen, with Glu neurons experiencing substantial hypomethylation events. The relationship between mCG levels and gene expression displays cell type-specific patterns, while genic CH methylation exhibits a negative effect on transcriptional abundance. We found that cell type-specific CG-DMRs are informative in terms of represented neuronal function. Furthermore, we observed that the identified Glu-specific hypo-DMRs have a high level of consistency with the chromatin accessibility of excitatory neurons and the regions enriched for histone modifications (H3K27ac and H3K4me1) of active enhancers, suggesting their regulatory potential. Hypomethylation regions specific to each cell type are predicted to bind neuron type-specific transcription factors. Finally, we show that the DNA methylation changes in a mouse model of Rett syndrome, a neurodevelopmental disorder caused by the de novo mutations in MECP2, are cell type- and brain region-specific. Our results suggest that cell type-specific DNA methylation signatures are associated with the functional characteristics of the neuronal subtypes. The presented results emphasize the importance of DNA methylation-mediated epigenetic regulation in neuronal diversity and disease.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Somatic activating Kirsten rat sarcoma viral oncogene homologue (KRAS) mutations have been reported in patients with arteriovenous malformations. By producing LSL-Kras (G12D); Cdh5 (PAC)-CreERT2 [iEC-Kras (G12D*)] mice, we hoped to activate KRAS within vascular endothelial cells (ECs) to generate an arteriovenous malformation mouse model. Neonatal mice were treated daily with tamoxifen from postnatal (PN) days 1-3. Mortality and phenotypes varied amongst iEC-Kras (G12D*) pups, with only 31.5% surviving at PN14. Phenotypes (focal lesions, vessel dilations) developed in a consistent manner, although with unpredictable severity within multiple soft tissues (such as the brain, liver, heart and brain). Overall, iEC-Kras (G12D*) pups developed significantly larger vascular lumen areas compared with control littermates, beginning at PN8. We subsequently tested whether the MEK inhibitor trametinib could effectively alleviate lesion progression. At PN16, iEC-Kras (G12D*) pup survival improved to 76.9%, and average vessel sizes were closer to controls than in untreated and vehicle-treated mutants. In addition, trametinib treatment helped normalize iEC-Kras (G12D*) vessel morphology in PN14 brains. Thus, trametinib could act as an effective therapy for KRAS-induced vascular malformations in patients.

Responding effectively to external stress is crucial for neurons. Defective stress granule dynamics has been hypothesized as one of the pathways that renders motor neurons in amyotrophic lateral sclerosis (ALS) more prone to early death. Specifically, it is thought that stress granules seed the cytoplasmic TDP-43 inclusions that are observed in the neurons of most ALS patients, as well as ~50% of all frontotemporal dementia (FTD) patients. In this study, we tested this hypothesis in an intact mammalian nervous system. We established an in vivo heat stress paradigm in mice that effectively triggers the eIF2α pathway and the formation of stress granules in the CNS. In non-transgenic mice, we report an age-dependent decline in the formation of heat-induced stress granules, with 18-month-old animals showing a significant impairment. Furthermore, although neuronal stress granules were robustly observed in non-transgenic mice and SOD1G93A mice, they were largely absent in age-matched TDP-43M337V animals. The observed defect in stress granule formation in TDP-43M337V mice correlated with deficits in expression of key protein components typically required for phase separation. Lastly, while TDP-43 was not localized to stress granules, we observed complete nuclear depletion of TDP-43 in a subset of neurons, with the highest proportion being in the TDP-43M337V mice. Overall, our results indicate that mutant TDP-43 expression is associated with defective stress granule assembly and increased TDP-43 nuclear depletion in the mammalian nervous system, which could be relevant to ALS/FTD pathogenesis.

Heterozygous germline missense variants in the HRAS gene underlie Costello syndrome (CS). The molecular basis for cutaneous manifestations in CS is largely unknown. We used an immortalized human cell line, HaCaT keratinocytes, stably expressing wild-type or CS-associated (p.Gly12Ser) HRAS and defined RIN1 as quantitatively most prominent, high-affinity effector of active HRAS in these cells. As an exchange factor for RAB5 GTPases, RIN1 is involved in endosomal sorting of cell-adhesion integrins. RIN1-dependent RAB5A activation was strongly increased by HRASGly12Ser, and HRAS-RIN1-ABL1/2 signaling was induced in HRASWT- and HRASGly12Ser-expressing cells. Along with that, HRASGly12Ser expression decreased total integrin levels and enriched β1 integrin in RAB5- and EEA1-positive early endosomes. The intracellular level of active β1 integrin was increased in HRASGly12Ser HaCaT keratinocytes due to impaired recycling, whereas RIN1 disruption raised β1 integrin cell surface distribution. HRASGly12Ser induced co-localization of β1 integrin with SNX17 and RAB7 in early/sorting and late endosomes, respectively. Thus, by retaining β1 integrin in intracellular endosomal compartments, HRAS-RIN1 signaling affects the subcellular availability of β1 integrin. This may interfere with integrin-dependent processes as we detected for HRASGly12Ser cells spreading on fibronectin. We conclude that dysregulation of receptor trafficking and integrin-dependent processes such as cell adhesion are relevant in the pathobiology of CS.

The neuronal-specific SNORD115 has gathered interest because its deficiency may contribute to the pathophysiology of Prader-Willi syndrome (PWS), possibly by altering post-transcriptional regulation of the gene encoding the serotonin (HTR2C) receptor. Yet, Snord115-KO mice do not resume the main symptoms of PWS, and only subtle-altered A-to-I RNA editing of Htr2c mRNAs was uncovered. Because HTR2C signaling fine-tunes the activity of monoaminergic neurons, we addressed the hypothesis that lack of Snord115 alters monoaminergic systems. We first showed that Snord115 was expressed in both monoaminergic and non-monoaminergic cells of the ventral tegmental area (VTA) and the dorsal raphe nucleus (DRN) harboring cell bodies of dopaminergic and serotonergic neurons, respectively. Measuring the tissue level of monoamines and metabolites, we found very few differences except that the content of homovanillic acid-a metabolite of dopamine-was decreased in the orbitofrontal and prefrontal cortex of Snord115-KO mice. The latter effects were, however, associated with a few changes in monoamine tissue content connectivity across the 12 sampled brain regions. Using in vivo single-cell extracellular recordings, we reported that the firing rate of VTA dopaminergic neurons and DRN serotonergic neurons was significantly increased in Snord115-KO mice. These neural circuit dysfunctions were not, however, associated with apparent defects in binge eating, conditioned place preference to cocaine, cocaine-induced hyperlocomotion or compulsive behavior. Altogether, our multiscale study shows that the absence of Snord115 impacts central monoaminergic circuits to an extent that does not elicit gross behavioral abnormalities.

Participant overlap can induce overfitting bias into Mendelian randomization (MR) and polygenic risk score (PRS) studies. Here, we evaluated a block jackknife resampling framework for genome-wide association studies (GWAS) and PRS construction to mitigate overfitting bias in MR analyses and implemented this study design in a causal inference setting using data from the UK Biobank. We simulated PRS and MR under three scenarios: (1) using weighted SNP estimates from an external GWAS, (2) using weighted SNP estimates from an overlapping GWAS sample and (3) using a block jackknife resampling framework. Based on a P-value threshold to derive genetic instruments for MR studies (P < 5 × 10-8) and a 10% variance in the exposure explained by all SNPs, block-jackknifing PRS did not suffer from overfitting bias (mean R2 = 0.034) compared with the externally weighted PRS (mean R2 = 0.040). In contrast, genetic instruments derived from overlapping samples explained a higher variance (mean R2 = 0.048) compared with the externally derived score. Overfitting became considerably more severe when using a more liberal P-value threshold to construct PRS (e.g. P < 0.05, overlapping sample PRS mean R2 = 0.103, externally weighted PRS mean R2 = 0.086), whereas estimates using jackknife score remained robust to overfitting (mean R2 = 0.084). Using block jackknife resampling MR in an applied analysis, we examined the effects of body mass index on circulating biomarkers which provided comparable estimates to an externally weighted instrument, whereas the overfitted scores typically provided narrower confidence intervals. Furthermore, we extended this framework into sex-stratified, multivariate and bidirectional settings to investigate the effect of childhood body size on adult testosterone levels.

Leber's hereditary optic neuropathy (LHON) is a maternally inherited eye disease that results from degeneration of retinal ganglion cells (RGC). Mitochondrial ND4 11778G > A mutation, which affects structural components of complex I, is the most prevalent LHON-associated mitochondrial DNA (mtDNA) mutation worldwide. The m.11778G > A mutation is the primary contributor underlying the development of LHON and X-linked PRICKLE3 allele (c.157C > T, p.Arg53Trp) linked to biogenesis of ATPase interacts with m.11778G > A mutation to cause LHON. However, the lack of appropriate cell and animal models of LHON has been significant obstacles for deep elucidation of disease pathophysiology, specifically the tissue-specific effects. Using RGC-like cells differentiated from induced pluripotent stem cells (iPSCs) from members of one Chinese family (asymptomatic subjects carrying only m.11778G > A mutation or PRICKLE3 p.Arg53Trp mutation, symptomatic individuals bearing both m.11778G > A and PRICKLE3 p.Arg53Trp mutations and control lacking these mutations), we demonstrated the deleterious effects of mitochondrial dysfunctions on the morphology and functions of RGCs. Notably, iPSCs bearing only m.11778G > A or p.Arg53Trp mutation exhibited mild defects in differentiation to RGC-like cells. The RGC-like cells carrying only m.11778G > A or p.Arg53Trp mutation displayed mild defects in RGC morphology, including the area of soma and numbers of neurites, electrophysiological properties, ATP contents and apoptosis. Strikingly, those RGC-like cells derived from symptomatic individuals harboring both m.11778G > A and p.Arg53Trp mutations displayed greater defects in the development, morphology and functions than those in cells bearing single mutation. These findings provide new insights into pathophysiology of LHON arising from RGC deficiencies caused by synergy between m.11778G > A and PRICKLE3 p.Arg53Trp mutation.

EFEMP1 R345W is a dominant mutation causing Doyne honeycomb retinal dystrophy/malattia leventinese (DHRD/ML), a rare blinding disease with clinical pathology similar to age-related macular degeneration (AMD). Aged Efemp1  R345W/R345W knock-in mice (Efemp1ki/ki) develop microscopic deposits on the basal side of retinal pigment epithelial cells (RPE), an early feature in DHRD/ML and AMD. Here, we assessed the role of alternative complement pathway component factor B (FB) in the formation of these deposits. RNA-seq analysis of the posterior eyecups revealed increased unfolded protein response, decreased mitochondrial function in the neural retina (by 3 months of age) and increased inflammatory pathways in both neural retina and posterior eyecups (at 17 months of age) of Efemp1ki/ki mice compared with wild-type littermate controls. Proteomics analysis of eye lysates confirmed similar dysregulated pathways as detected by RNA-seq. Complement activation was increased in aged Efemp1ki/ki eyes with an approximately 2-fold elevation of complement breakdown products iC3b and Ba (P < 0.05). Deletion of the Cfb gene in female Efemp1ki/ki mice partially normalized the above dysregulated biological pathway changes and oral dosing of a small molecule FB inhibitor from 10 to 12 months of age reduced sub-RPE deposits by 65% (P = 0.029). In contrast, male Efemp1ki/ki mice had fewer sub-RPE deposits than age-matched females, no elevation of ocular complement activation and no effect of FB inhibition on sub-RPE deposits. The effects of FB deletion or inhibition on Efemp1ki/ki mice supports systemic inhibition of the alternative complement pathway as a potential treatment of dry AMD and DHRD/ML.

Stress granules are membrane-less ribonucleoprotein organelles that assemble upon exposure to stress conditions, but rapidly disassemble upon removal of stress. However, chronic stress can lead to persistent stress granules, a feature of distinct age-related neurodegenerative disorders. Among them, Huntington's disease (HD) which is caused by mutant expansion of the polyglutamine (polyQ) repeats of huntingtin protein (HTT), leading to its aggregation. To identify modulators of mutant HTT aggregation, we define its interactome in striatal neurons from patient-derived induced pluripotent stem cells (HD-iPSCs). We find that HTT interacts with G3BP1, a characteristic component of stress granules. Knockdown of G3BP1 increases mutant HTT protein levels and abolishes the ability of iPSCs as well as their differentiated neural counterparts to suppress mutant HTT aggregation. Moreover, loss of G3BP1 hastens polyQ-expanded aggregation and toxicity in the neurons of HD C. elegans models. Likewise, the assembly of G3BP1 into stress granules upon distinct stress conditions also reduces its interaction with HTT in human cells, promoting mutant HTT aggregation. Notably, enhancing the levels of G3BP1 is sufficient to induce proteasomal degradation of mutant HTT and prevent its aggregation, whereas the formation of stress granules blocks these ameliorative effects. In contrast, a mutant G3BP1 variant that cannot accumulate into granules retains its capacity to prevent mutant HTT aggregation even when the cells assemble stress granules. Thus, our findings indicate a direct role of G3BP1 and stress granule assembly in mutant HTT aggregation that may have implications for HD.

Leber's hereditary optic neuropathy (LHON) is a maternally transmitted eye disease due to the degeneration of retinal ganglion cells (RGC). Mitochondrial 11778G > A mutation is the most common LHON-associated mitochondrial DNA (mtDNA) mutation. Our recent studies demonstrated some LHON families manifested by synergic interaction between m.11778G > A mutation and YARS2 allele (c.572G > T, p.Gly191Val) encoding mitochondrial tyrosyl-tRNA synthetase. However, the RGC-specific effects of LHON-associated mtDNA mutations remains elusive and there is no highly effective therapy for LHON. Here, we generated patients-derived induced pluripotent stem cells (iPSCs) from fibroblasts derived from a Chinese LHON family (both m.11778G > A and c.572G > T mutations, only m.11778G > A mutation, and control subject). The c.572G > T mutation in iPSC lines from a syndromic individual was corrected by CRISPR/Cas9. Those iPSCs were differentiated into neural progenitor cells (NPCs) and subsequently induced RGC-like cells using a stepwise differentiation procedure. Those RGC-like cells derived from symptomatic individual harboring both m.11778G > A and c.572G > T mutations exhibited greater defects in neuronal differentiation, morphology including reduced area of soma, numbers of neurites, and shortened length of axons, electrophysiological properties than those in cells bearing only m.11778G > A mutation. Furthermore, these RGC-like cells revealed more drastic reductions in oxygen consumption rates, levels of mitochondrial ATP and increasing productions of reactive oxygen species than those in other cell models. These mitochondrial dysfunctions promoted the apoptotic process for RGC degenerations. Correction of YARS2 c.572G > T mutation rescued deficiencies of patient-derived RGC-like cells. These findings provide new insights into pathophysiology of LHON arising from RGC-specific mitochondrial dysfunctions and step toward therapeutic intervention for this disease.

The abnormal expansion of GGGGCC hexanucleotide repeats within the C9orf72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The accumulation of GGGGCC repeat-containing RNAs as RNA foci, and the deposition of dipeptide repeat proteins (DPR) produced from these repeat RNAs by unconventional translation are major pathological hallmarks of C9orf72-linked ALS/FTD (C9-ALS/FTD), and are both thought to play a crucial role in the pathogenesis of these diseases. Because GGGGCC repeat RNA is likely to be the most upstream therapeutic target in the pathogenic cascade of C9-ALS/FTD, lowering the cellular level of GGGGCC repeat RNA is expected to mitigate repeat RNA toxicity, and will therefore be a disease-modifying therapeutic strategy for the treatment of C9-ALS/FTD. In this study, we demonstrated using a Drosophila model of C9-ALS/FTD that elevated expression of a subset of human RNA-binding proteins that bind to GGGGCC repeat RNA, including hnRNPA3, IGF2BP1, hnRNPA2B1, hnRNPR, and SF3B3, reduces the level of GGGGCC repeat RNA, resulting in the suppression of neurodegeneration. We further showed that hnRNPA3-mediated reduction of GGGGCC repeat RNA suppresses disease pathology, such as RNA foci and DPR accumulation. These results demonstrate that hnRNPA3 and other RNA-binding proteins negatively regulate the level of GGGGCC repeat RNA, and mitigate repeat RNA toxicity in vivo, indicating the therapeutic potential of the repeat RNA-lowering approach mediated by endogenous RNA-binding proteins for the treatment of C9-ALS/FTD.

Drosophila is an excellent model organism for studying human neurodegenerative diseases (NDs). However, there is still almost no experimental system which could directly observe the degeneration of neurons and automatically quantify axonal degeneration. In this study, we created MeDUsA (a 'method for the quantification of degeneration using fly axons'), a standalone executable computer program based on Python that combines a pre-trained deep-learning masking tool with an axon terminal counting tool. This software automatically quantifies the number of retinal R7 axons in Drosophila from a confocal z-stack image series. Using this software, we were able to directly demonstrate that axons were degenerated by the representative causative genes of NDs for the first time in Drosophila. The fly retinal axon is an excellent experimental system that is capable of mimicking the pathology of axonal degeneration in human NDs. MeDUsA rapidly and accurately quantifies axons in Drosophila photoreceptor neurons. It enables large-scale research into axonal degeneration, including screening to identify genes or drugs that mediate axonal toxicity caused by ND proteins and diagnose the pathological significance of novel variants of human genes in axons.

Congenital myasthenic syndrome (CMS) is a heterogeneous condition associated with 34 different genes, including SLC5A7, which encodes the high affinity choline transporter 1 (CHT1). CHT1 is expressed in presynaptic neurons of the neuromuscular junction where it uses the inward sodium gradient to re-uptake choline. Bi-allelic CHT1 mutations often lead to neonatal lethality, and less commonly to non-lethal motor weakness and developmental delays. Here, we report detailed biochemical characterization of two novel mutations in CHT1, p.I294T and p.D349N, that we identified in an 11 year-old patient with a history of neonatal respiratory distress, and subsequent hypotonia and global developmental delay. Heterologous expression of each CHT1 mutant in human embryonic kidney cells showed two different mechanisms of reduced protein function. The p.I294T CHT1 mutant transporter function was detectable, but its abundance and half-life were significantly reduced. In contrast, the p.D349N CHT1 mutant was abundantly expressed at the cell membrane, but transporter function was absent. The residual function of the p.I294T CHT1 mutant may explain the non-lethal form of CMS in this patient, and the divergent mechanisms of reduced CHT1 function that we identified may guide future functional studies of the CHT1 myasthenic syndrome. Based on these in vitro studies that provided a diagnosis, treatment with cholinesterase inhibitor together with physical and occupational therapy significantly improved the patient's strength and quality of life.

Biallelic germline mutations in BRCA2 occur in the FA-D1 subtype of the rare pediatric disorder, Fanconi anemia (FA), characterized clinically by congenital abnormalities and a high propensity to develop malignancies early in life. Clinical and genetic data from 96 FA-D1 patients with biallelic BRCA2 mutations were collected and used to develop a new cancer risk prediction score system based on the specific mutations in BRCA2. This score takes into account the location of frameshift/stop and missense mutations relative to exon 11 of BRCA2, which encodes the major sites for interaction with the RAD51 recombinase, and uses the MaxEnt and HBond splicing scores to analyze potential splice site pertubations. Among 75 FA-D1 patients with ascertained BRCA2 mutations, 66 patients developed 102 malignancies, ranging from one to three independent tumors per individual. The median age at the manifestation of peripheral embryonal tumors was 1.0, of hematologic malignancies 1.8, and of CNS tumors 2.7 years, respectively. Patients who received treatment lived longer than those without. Using our novel scoring system, we could distinguish three distinct cancer risk groups among FA-D1 patients: in the first, patients developed their initial malignancy at a median age of 1.3 years (n = 36, 95% CI 0.9-1.8), in the second group at 2.3 years (n = 17, 95% CI 1.4-4.4) and in the 3rd group at 23.0 years (n = 22, 95% CI 4.3-n/a). Therefore, this scoring system allows, for the first time, to predict the cancer manifestation of FA-D1 patients simply based on the type and position of the mutations in BRCA2.

Tandem CAG repeat sizes of 36 or more in the huntingtin gene (HTT) cause Huntington disease. Apart from neuropsychiatric complications, the disease is also accompanied by metabolic dysregulation and weight loss, which contribute to a progressive functional decline. Recent studies also reported an association between repeats below the pathogenic threshold (<36) for Huntington's disease and body mass index (BMI), suggesting that HTT repeat sizes in the non-pathogenic range are associated with metabolic dysregulation. In this study we hypothesized that HTT repeat sizes < 36 are associated with metabolite levels, possibly mediated through reduced BMI. We pooled data from three European cohorts (n = 10 228) with genotyped HTT CAG repeat size and metabolomic measurements. All 145 metabolites were measured on the same targeted platform in all studies. Multilevel mixed-effects analysis using the CAG repeat size in HTT identified 67 repeat size-metabolite associations. Overall, the metabolomic profile associated with larger CAG repeat sizes in HTT were unfavorable-similar to those of higher risk of coronary artery disease and type 2 diabetes-and included elevated levels of amino acids, fatty acids, LDL, VLDL and IDL related metabolites whilst with decreased levels of very large HDL related metabolites. Furthermore, the associations of 50 metabolites, in particular specific very large HDL related metabolites, were mediated by lower BMI. However, no mediation effect was found for 17 metabolites related to LDL and IDL. In conclusion, our findings indicate that large non-pathogenic CAG repeat sizes in HTT are associated with an unfavorable metabolomic profile despite their association with a lower BMI.

Pathogenic variations in the sodium voltage-gated channel alpha subunit 1 (SCN1A) gene are responsible for multiple epilepsy phenotypes, including Dravet syndrome (DS), febrile seizures (FS), and genetic epilepsy with febrile seizures plus (GEFS+). Phenotypic heterogeneity is a hallmark of SCN1A-related epilepsies, the causes of which are yet to be clarified. Genetic variation in the non-coding regulatory regions of SCN1A could be one potential causal factor. However, a comprehensive understanding of the SCN1A regulatory landscape is currently lacking. Here, we summarised the current state of knowledge of SCN1A regulation, providing details of its promoter and enhancer regions. We then integrated currently available data on SCN1A promoters by extracting information related to the SCN1A locus from genome-wide repositories, and clearly defined the promoter and enhancer regions of SCN1A. Further, we explored the cellular specificity of differential SCN1A promoter usage. We also reviewed and integrated the available human brain-derived enhancer databases and mouse-derived data to provide a comprehensive computationally-developed summary of SCN1A brain-active enhancers. By querying genome-wide data repositories, extracting SCN1A-specific data and integrating the different types of independent evidence, we created a comprehensive catalogue that better defines the regulatory landscape of SCN1A, which could be used to explore the role of SCN1A regulatory regions in disease.

Berardinelli-Seip Congenital Lipodystrophy type 2 (CGL2) is a very rare human genetic disorder with potential significance to the understanding of the pathobiology of aging. CGL2 patients display characteristic progeroid features and suffer from type 2 diabetes, insulin resistance, and fatty liver. In this study, we profiled genome-wide DNA methylation levels in CGL2 patients with BSCL2 mutations to study epigenetic age acceleration and DNA methylation alterations. This analysis revealed significant age acceleration in blood DNA of CGL2 patients using both first and second generation epigenetic clocks. We also observed a shortened lifespan of C. elegans following knockdown of the BSCL2 homolog seip-1 on a daf-16/Foxo mutant background. DNA methylation analysis revealed significant differentially methylated sites enriched for lyase activity, kinase regulator activity, protein kinase regulator activity, and kinase activator activity. We could also observe significant hypomethylation in the promoter of the Dual Specificity Phosphatase 22 (DUSP22) gene when comparing CGL2 patients vs controls. We conclude that in line with the observed progeroid features, CGL2 patients exhibit significant epigenetic age acceleration and DNA methylation alterations that might affect pathways/genes of potential relevance to the disease.

Oxidative stress is a common feature of inflammation-driven cancers, and it promotes genomic instability and aggressive tumour phenotypes. It is known that oxidative stress transiently modulates gene expression through the oxidation of transcription factors and associated regulatory proteins. Neutrophils are our most abundant white blood cells and accumulate at sites of infection and inflammation. Activated neutrophils produce hypochlorous acid and chloramines, which can disrupt DNA methylation by oxidizing methionine. The goal of the current study was to determine whether chloramine exposure results in sequence-specific modifications in DNA methylation that enable long-term alterations in transcriptional output. Proliferating Jurkat T-lymphoma cells were exposed to sublethal doses of glycine chloramine and differential methylation patterns were compared using Illumina EPIC 850 K bead chip arrays. There was a substantial genome-wide decrease in methylation 4 h after exposure that correlated with altered RNA expression for 24 and 48 h, indicating sustained impacts on exposed cells. A large proportion of the most significant differentially methylated CpG sites were situated towards chromosomal ends, suggesting that these regions are most susceptible to inhibition of maintenance DNA methylation. This may contribute to epigenetic instability of chromosomal ends in rapidly dividing cells, with potential implications for the regulation of telomere length and cellular longevity.

DEPDC5 (DEP Domain-Containing Protein 5) encodes an inhibitory component of the mammalian target of rapamycin (mTOR) pathway and is commonly implicated in sporadic and familial focal epilepsies, both non-lesional and in association with focal cortical dysplasia. Germline pathogenic variants are typically heterozygous and inactivating. We describe a novel phenotype caused by germline biallelic missense variants in DEPDC5. Cases were identified clinically. Available records, including magnetic resonance imaging and electroencephalography, were reviewed. Genetic testing was performed by whole exome and whole-genome sequencing and cascade screening. In addition, immunohistochemistry was performed on skin biopsy. The phenotype was identified in nine children, eight of which are described in detail herein. Six of the children were of Irish Traveller, two of Tunisian and one of Lebanese origin. The Irish Traveller children shared the same DEPDC5 germline homozygous missense variant (p.Thr337Arg), whereas the Lebanese and Tunisian children shared a different germline homozygous variant (p.Arg806Cys). Consistent phenotypic features included extensive bilateral polymicrogyria, congenital macrocephaly and early-onset refractory epilepsy, in keeping with other mTOR-opathies. Eye and cardiac involvement and severe neutropenia were also observed in one or more patients. Five of the children died in infancy or childhood; the other four are currently aged between 5 months and 6 years. Skin biopsy immunohistochemistry was supportive of hyperactivation of the mTOR pathway. The clinical, histopathological and genetic evidence supports a causal role for the homozygous DEPDC5 variants, expanding our understanding of the biology of this gene.

Crohn's disease (CD) and ulcerative colitis (UC), two major subtypes of inflammatory bowel disease, show substantial differences in their clinical course and treatment response. To identify the genetic factors underlying the distinct characteristics of these two diseases, we performed a genome-wide association study (GWAS) between CD (n = 2359) and UC (n = 2175) in a Korean population, followed by replication in an independent sample of 772 CD and 619 UC cases. Two novel loci were identified with divergent effects on CD and UC: rs9842650 in CD200 and rs885026 in NCOR2. In addition, the seven established susceptibility loci [major histocompatibility complex (MHC), TNFSF15, OTUD3, USP12, IL23R, FCHSD2 and RIPK2] reached genome-wide significance. Of the nine loci, six (MHC, TNFSF15, OTUD3, USP12, IL23R and CD200) were replicated in the case-case GWAS of European populations. The proportion of variance explained in CD-UC status by polygenic risk score analysis was up to 22.6%. The area under the receiver-operating characteristic curve value was 0.74, suggesting acceptable discrimination between CD and UC. This CD-UC GWAS provides new insights into genetic differences between the two diseases with similar symptoms and might be useful in improving their diagnosis and treatment.

BACKGROUND: The relations of alcohol consumption and gene expression remain to be elucidated. MATERIALS AND METHODS: We examined cross-sectional associations between alcohol consumption and whole blood derived gene expression levels and between alcohol-associated genes and obesity, hypertension, and diabetes in 5531 Framingham Heart Study (FHS) participants. RESULTS: We identified 25 alcohol-associated genes. We further showed cross-sectional associations of 16 alcohol-associated genes with obesity, nine genes with hypertension, and eight genes with diabetes at P < 0.002. For example, we observed decreased expression of PROK2 (β = -0.0018; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) and PAX5 (β = -0.0014; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) per 1 g/day increase in alcohol consumption. Consistent with our previous observation on the inverse association of alcohol consumption with obesity and positive association of alcohol consumption with hypertension, we found that PROK2 was positively associated with obesity (OR = 1.42; 95%CI: 1.17, 1.72; P = 4.5e - 4) and PAX5 was negatively associated with hypertension (OR = 0.73; 95%CI: 0.59, 0.89; P = 1.6e - 3). We also observed that alcohol consumption was positively associated with expression of ABCA13 (β = 0.0012; 95%CI: 0.0007, 0.0017; P = 1.3e - 6) and ABCA13 was positively associated with diabetes (OR = 2.57; 95%CI: 1.73, 3.84; P = 3.5e - 06); this finding, however, was inconsistent with our observation of an inverse association between alcohol consumption and diabetes. CONCLUSIONS: We showed strong cross-sectional associations between alcohol consumption and expression levels of 25 genes in FHS participants. Nonetheless, complex relationships exist between alcohol-associated genes and CVD risk factors.

Genome-wide association studies (GWAS) have markedly advanced our understanding of the genetics of Parkinson's disease (PD), but they currently do not account for the full heritability of PD. In many cases it is difficult to unambiguously identify a specific gene within each locus because GWAS does not provide functional information on the identified candidate loci. Here we present an integrative approach that combines transcriptome-wide association study (TWAS) with high-throughput neuronal dysfunction analyses in Drosophila to discover and validate candidate PD genes. We identified 160 candidate genes whose misexpression is associated with PD risk via TWAS. Candidates were validated using orthogonal in silico methods and found to be functionally related to PD-associated pathways (i.e. endolysosome). We then mimicked these TWAS-predicted transcriptomic alterations in a Drosophila PD model and discovered that 50 candidates can modulate α-Synuclein(α-Syn)-induced neurodegeneration, allowing us to nominate new genes in previously known PD loci. We also uncovered additional novel PD candidate genes within GWAS suggestive loci (e.g. TTC19, ADORA2B, LZTS3, NRBP1, HN1L), which are also supported by clinical and functional evidence. These findings deepen our understanding of PD, and support applying our integrative approach to other complex trait disorders.

Mutations and aberrant gene expression during cellular differentiation lead to neurodevelopmental disorders, such as Prader-Willi syndrome (PWS), which results from the deletion of an imprinted locus on paternally inherited chromosome 15. We analyzed chromatin-associated RNA in human induced pluripotent cells (iPSCs) upon depletion of hybrid small nucleolar long non-coding RNAs (sno-lncRNAs) and 5' snoRNA capped and polyadenylated long non-coding RNAs (SPA-lncRNAs) transcribed from the locus deleted in PWS. We found that rapid ablation of these lncRNAs affects transcription of specific gene classes. Downregulated genes contribute to neurodevelopment and neuronal maintenance, while upregulated genes are predominantly involved in the negative regulation of cellular metabolism and apoptotic processes. Our data reveal the importance of SPA-lncRNAs and sno-lncRNAs in controlling gene expression in iPSCs and provide a platform for synthetic experimental approaches in PWS studies. We conclude that ncRNAs transcribed from the PWS locus are critical regulators of a transcriptional signature, which is important for neuronal differentiation and development.

Myotonic dystrophy type 1 is a complex disease caused by a genetically unstable CTG repeat expansion in the 3'-untranslated region of the DMPK gene. Age-dependent, tissue-specific somatic instability has confounded genotype-phenotype associations, but growing evidence suggests that it also contributes directly toward disease progression. Using a well-characterized clinical cohort of DM1 patients from Costa Rica, we quantified somatic instability in blood, buccal cells, skin and skeletal muscle. Whilst skeletal muscle showed the largest expansions, modal allele lengths in skin were also very large and frequently exceeded 2000 CTG repeats. Similarly, the degree of somatic expansion in blood, muscle and skin were associated with each other. Notably, we found that the degree of somatic expansion in skin was highly predictive of that in skeletal muscle. More importantly, we established that individuals whose repeat expanded more rapidly than expected in one tissue (after correction for progenitor allele length and age) also expanded more rapidly than expected in other tissues. We also provide evidence suggesting that individuals in whom the repeat expanded more rapidly than expected in skeletal muscle have an earlier age at onset than expected (after correction for the progenitor allele length). Pyrosequencing analyses of the genomic DNA flanking the CTG repeat revealed that the degree of methylation in muscle was well predicted by the muscle modal allele length and age, but that neither methylation of the flanking DNA nor levels of DMPK sense and anti-sense transcripts could obviously explain individual- or tissue-specific patterns of somatic instability.

Human spermatogenesis requires an orchestrated expression of numerous genes in various germ cell subtypes. Therefore, the genetic landscape of male infertility is highly complex. Known genetic factors alone account for at least 15% of male infertility. However, ~40% of infertile men remain undiagnosed and are classified as idiopathic infertile men. We performed exome sequencing in 47 idiopathic infertile men (discovery cohort), followed by replication study (40 variants in 33 genes) in 844 infertile men and 709 controls using Sequenom MassARRAY® based genotyping. We report 17 variants in twelve genes that comprise both previously reported (DNAH8, DNAH17, FISP2 and SPEF2) and novel candidate genes (BRDT, CETN1, CATSPERD, GMCL1, SPATA6, TSSK4, TSKS and ZNF318) for male infertility. The latter have a strong biological nexus to human spermatogenesis and their respective mouse knockouts are concordant with human phenotypes. One candidate gene CETN1, identified in this study, was sequenced in another independent cohort of 840 infertile and 689 fertile men. Further, CETN1 variants were functionally characterized using biophysical and cell biology approaches. We demonstrate that CETN1 variant- p.Met72Thr leads to multipolar cells, fragmented nuclei during mitosis leading to cell death and show significantly perturbed ciliary disassembly dynamics. Whereas CETN1-5' UTR variant; rs367716858 leads to loss of a methylation site and increased reporter gene expression in vitro. We report a total of eight novel candidate genes identified by exome sequencing, which may have diagnostic relevance and can contribute to improved diagnostic workup and clinical management of male infertility.

Erratum in Hum Mol Genet. 2022 Dec 26;:

The purpose of this paper is to identify likely pathogenic non-coding variants in inherited retinal dystrophy (IRD) genes, using genome sequencing (GS). Patients with IRD were recruited to the study and underwent comprehensive ophthalmological evaluation and GS. The results of GS were investigated through virtual gene panel analysis, and plausible pathogenic variants and clinical phenotype evaluated by the multidisciplinary team (MDT) discussion. For unsolved patients in whom a specific gene was suspected to harbor a missed pathogenic variant, targeted re-analysis of non-coding regions was performed on GS data. Candidate variants were functionally tested by messenger RNA analysis, minigene or luciferase reporter assays. Previously unreported, likely pathogenic, non-coding variants in 7 genes (PRPF31, NDP, IFT140, CRB1, USH2A, BBS10 and GUCY2D), were identified in 11 patients. These were shown to lead to mis-splicing (PRPF31, IFT140, CRB1 and USH2A) or altered transcription levels (BBS10 and GUCY2D). MDT-led, phenotype-driven, non-coding variant re-analysis of GS is effective in identifying the missing causative alleles.

Targeting AMP-activated protein kinase (AMPK) is emerging as a promising strategy for treating myotonic dystrophy type 1 (DM1), the most prevalent form of adult-onset muscular dystrophy. We previously demonstrated that 5-aminomidazole-4-carboxamide-1-β-D-ribofuranoside (AICAR) and exercise, two potent AMPK activators, improve disease features in DM1 mouse skeletal muscles. Here, we employed a combinatorial approach with these AMPK activators and examined their joint impact on disease severity in male and female DM1 mice. Our data reveal that swimming exercise additively enhances the effect of AICAR in mitigating the nuclear accumulation of toxic CUGexp RNA foci. In addition, our findings show a trend towards an enhanced reversal of MBNL1 sequestration and correction in pathogenic alternative splicing events. Our results further demonstrate that the combinatorial impact of exercise and AICAR promotes muscle fiber hypertrophy in DM1 skeletal muscle. Importantly, these improvements occur in a sex-specific manner with greater benefits observed in female DM1 mice. Our findings demonstrate that combining AMPK-activating interventions may prove optimal for rescuing the DM1 muscle phenotype and uncover important sex differences in the response to AMPK-based therapeutic strategies in DM1 mice.

Oligoasthenoteratozoospermia (OAT) can result in male infertility due to reduced sperm motility and abnormal spermatozoan morphology. The Tektins are a family of highly conserved filamentous proteins expressed in the axoneme and associated structures in many different metazoan species. Earlier studies on mice identified Tektin3 (Tekt3) as a testis-enriched gene, and knockout of Tekt3 resulted in asthenozoospermia in the mice. Here, whole exome sequencing of 100 males with asthenozoospermia from unrelated families was performed, followed by Sanger sequencing, leading to the identification of TEKT3 as a candidate gene in two of these patients and their associated family members. In total, three mutations in the TEKT3 gene were identified in both these patients, including one homozygous deletion-insertion mutation (c.543_547delinsTTGAT: p.Glu182*) and one compound heterozygous mutation (c.[548G > A]; [752A > C], p.[(Arg183Gln)]; [(Gln251Pro)]). Both of these mutations resulted in the complete loss of TEKT3 expression. The patients were both found to produce sperm that, although they showed no apparent defects in the flagellar structure, had reduced progressive motility. In contrast to mice, most sperm from these two patients exhibited acrosomal hypoplasia, although this did not prevent the use of the sperm for in vitro fertilization through an ICSI approach. TEKT3 was found to bind to other TEKT proteins, suggesting that these proteins form a complex within human spermatozoa. Overall, these results suggest that a loss of TEKT3 function can contribute to OAT incidence in humans. TEKT3 deficiencies can reduce sperm motility and contribute to severe acrosomal hypoplasia in spermatozoa, compromising their normal function.

The testis-specific adenosine deaminase domain-containing (ADAD) protein family, including ADAD1 and ADAD2, has been confirmed to be essential in mouse male fertility. However, the roles of ADAD1 and ADAD2 in human reproductive biology are unclear. Herein, whole exome sequencing (WES) was conducted for 337 infertile patients to detect pathogenic variants in ADAD1 and ADAD2. Importantly, a novel deleterious biallelic variant of NM_001159285.2:c.1408G > T (p.V470F) and NM_001159285.2:c.1418A > G (p.E473G) in ADAD1 and a pathogenic homozygous missense variant of NM_001145400.2:c.1381C > T (p.R461W) in ADAD2 were identified in this infertile cohort with frequencies of 0.29% (1/337) and 0.59% (2/337), respectively. Electron microscopy revealed an abnormal morphology and severely disorganized ultrastructure of sperm from the patients. Immunofluorescence and western blotting showed a sharp decrease in ADAD1 and ADAD2 expression in sperm from the patients. Mechanistically, bioinformatics analysis suggested that ADAD2 interacts with DNAH17. Furthermore, we demonstrated that the expression of DNAH17 was markedly downregulated in the sperm of patients harboring ADAD2 variants. In addition, the expression of several autophagy regulators was significantly disrupted in the sperm of patients harboring ADAD2 variants. In conclusion, we identified novel ADAD1 and ADAD2 variants in three infertile patients from a large infertile cohort, first providing evidence that ADAD1 and ADAD2 variants might be a candidate genetic cause of human male infertility. Moreover, an important new dimension to our understanding of the genotype-phenotype correlations between the ADAD gene family and male infertility in humans has been uncovered, providing valuable information for the genetic diagnosis of male infertility.

The large DMD gene encodes a group of dystrophin proteins in brain and retina, produced from multiple promoters and alternative splicing events. Dystrophins are core components of different scaffolding complexes in distinct cell types. Their absence may thus alter several cellular pathways, which might explain the heterogeneous genotype-phenotype relationships underlying central comorbidities in Duchenne muscular dystrophy (DMD). However, the cell-specific expression of dystrophins and associated proteins (DAPs) is still largely unknown. The present study provides a first RNA-Seq-based reference showing tissue- and cell-specific differential expression of dystrophins, splice variants and DAPs in mouse brain and retina. We report that a cell type may express several dystrophin complexes, perhaps due to expression in separate cell subdomains and/or subpopulations, some of which with differential expression at different maturation stages. We also identified new splicing events in addition to the common exon-skipping events. These include a new exon within intron 51 (E51b) in frame with the flanking exons in retina, as well as inclusions of intronic sequences with stop codons leading to the presence of transcripts with elongated exons 40 and/or 41 (E40e, E41e) in both retina and brain. PCR validations revealed that the new exons may affect several dystrophins. Moreover, immunoblot experiments using a combination of specific antibodies and dystrophin-deficient mice unveiled that the transcripts with stop codons are translated into truncated proteins lacking their C-terminus, which we called N-Dp427 and N-Dp260. This study thus uncovers a range of new findings underlying the complex neurobiology of DMD.

DNA damage response (DDR) is an important signaling-transduction network that promotes the repair of DNA lesions which can induce and/or support diseases. However, the mechanisms involved in its regulation are not fully understood. Recent studies suggest that the peroxiredoxin 5 (Prdx5) enzyme, which detoxifies reactive oxygen species, is associated to genomic instability and signal transduction. Its role in the regulation of DDR, however, is not well characterized. In this study, we demonstrate a role of Prdx5 in the regulation of the DDR signaling pathway. Knockdown of Prdx5 resulted in DNA damage manifested by the induction of phosphorylated histone H2AX (γ-H2AX) and p53-binding protein 1 (53BP1). We show that Prdx5 regulates DDR through (1) polo-like kinase 1 (Plk1) mediated phosphorylation of ataxia telangiectasia mutated (ATM) kinase to further trigger downstream mediators Chek1 and Chek2; (2) the increase of the acetylation of p53 at lysine 382, stabilizing p53 in the nucleus and enhancing transcription and (3) the induction of autophagy, which regulates the recycling of molecules involved in DDR. We identified Sirt2 as a novel deacetylase of p53 at lysine 382, and Sirt2 regulated the acetylation status of p53 at lysine 382 in a Prdx5-dependent manner. Furthermore, we found that exogenous expression of Prdx5 decreased DNA damage and the activation of ATM in Pkd1 mutant renal epithelial cells, suggesting that Prdx5 may play a protective role from DNA damage in cystic renal epithelial cells. This study identified a novel mechanism of Prdx5 in the regulation of DDR through the ATM/p53/Sirt2 signaling cascade.

BACKGROUND: Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations. METHODS: We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants. RESULTS: Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22. CONCLUSION: Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

Isolated hypogonadotropic hypogonadism (IHH) is a rare disease with hypogonadism and infertility caused by the defects in embryonic migration of hypothalamic gonadotropin-releasing hormone (GnRH) neurons, hypothalamic GnRH secretion or GnRH signal transduction. PROKR2 gene, encoding a G-protein coupled receptor PROKR2, is one of the most frequently mutated genes identified in IHH patients. However, the functional consequences of several PROKR2 mutants remain elusive. In this study, we systematically analyzed the Gαq, Gαs and ERK1/2 signaling of 23 IHH-associated PROKR2 mutations which are yet to be functionally characterized. We demonstrate that blockage of Gαq, instead of MAPK/ERK pathway, inhibited PROK2-induced migration of PROKR2-expressing cells, implying that PROKR2-related IHH results primarily due to Gαq signaling pathway disruption. Combined with previous reports, we categorized a total of 63 IHH-associated PROKR2 mutations into four distinct groups according Gαq pathway functionality: (1) Neutral (N, > 80% activity); (2) Low pathogenicity (L, 50-80% activity); (3) Medium pathogenicity (M, 20-50% activity) and (4) High pathogenicity (H, < 20% activity). We further compared the cell-based functional results with in silico mutational prediction programs. Our results indicated that while SIFT predictions were accurate for transmembrane region mutations, mutations localized in the intracellular and extracellular domains were accurately predicted by the CADD prediction tool. Our results thus provide a functional database that can be used to guide diagnosis and appropriate genetic counseling in IHH patients with PROKR2 mutations.

Nemaline myopathy is a rare neuromuscular disorder associated with congenital or childhood-onset of skeletal muscle weakness and hypotonia that results in limited motor function. Nemaline myopathy is a genetic disorder and mutations in 12 genes are known to contribute to autosomal dominant or recessive forms of the disease. Recessive mutations in nebulin (NEB) are the most common cause of nemaline myopathy affecting about 50% of patients. Due to the large size of the NEB gene and lack of mutational hot spots, developing therapies that can benefit a wide group of patients is challenging. Although there are several promising therapies under investigation, there is no cure for nemaline myopathy. Therefore, targeting disease modifiers that can stabilize or improve skeletal muscle function may represent alternative therapeutic strategies. Our studies have identified Nrap upregulation in nebulin deficiency that contributes to structural and functional deficits in nemaline myopathy. We show that genetic ablation of nrap in nebulin deficiency restored sarcomeric disorganization, reduced protein aggregates and improved skeletal muscle function in zebrafish. Our findings suggest that Nrap is a disease modifier that affects skeletal muscle structure and function in nemaline myopathy, thus therapeutic targeting of Nrap in nebulin related nemaline myopathy and related diseases may be beneficial for patients.

Neuroinflammation and immune dysregulation play a key role in Alzheimer's disease (ad) and are also associated with severe Covid-19 and neurological symptoms. Also, genome-wide association studies found many risk SNPs for ad and Covid-19. However, our understanding of underlying gene regulatory mechanisms from risk SNPs to ad, Covid-19 and phenotypes is still limited. To this end, we performed an integrative multi-omics analysis to predict gene regulatory networks for major brain regions from population data in ad. Our networks linked transcription factors (TFs) to TF binding sites (TFBSs) on regulatory elements to target genes. Comparative network analyses revealed cross-region-conserved and region-specific regulatory networks, in which many immunological genes are present. Furthermore, we identified a list of ad-Covid genes using our networks involving known ad and Covid-19 genes. Our machine learning analysis prioritized 36 ad-Covid candidate genes for predicting Covid severity. Our independent validation analyses found that these genes outperform known genes for classifying Covid-19 severity and ad. Finally, we mapped GWAS SNPs of ad and severe Covid that interrupt TFBSs on our regulatory networks, revealing potential mechanistic insights of those disease risk variants. Our analyses and results are open-source available, providing an ad-Covid functional genomic resource at the brain-region level.

Membrane fusion is mediated by soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins. During neurotransmitter exocytosis, SNARE proteins on a synaptic vesicle and the target membrane form a complex, resulting in neurotransmitter release. N-ethylmaleimide-sensitive factor (NSF), a homohexameric ATPase, disassembles the complex, allowing individual SNARE proteins to be recycled. Recently, the association between pathogenic NSF variants and developmental and epileptic encephalopathy (DEE) was reported; however, the molecular pathomechanism of NSF-related DEE remains unclear. Here, three patients with de novo heterozygous NSF variants were presented, of which two were associated with DEE and one with a very mild phenotype. One of the DEE patients also had hypocalcemia from parathyroid hormone deficiency and neuromuscular junction impairment. Using PC12 cells, a neurosecretion model, we show that NSF with DEE-associated variants impaired the recycling of vesicular membrane proteins and vesicle enlargement in response to exocytotic stimulation. In addition, DEE-associated variants caused neurodegenerative change and defective autophagy through overactivation of the mTOR pathway. Treatment with rapamycin, an mTOR inhibitor, or overexpression of wild-type NSF ameliorated these phenotypes. Furthermore, neurons differentiated from patient-derived induced pluripotent stem cells showed neurite degeneration, which was also alleviated by rapamycin treatment or gene correction using genome editing. Protein structure analysis of NSF revealed that DEE-associated variants might disrupt the transmission of the conformational change of NSF monomers and consequently halt the rotation of ATP hydrolysis, indicating a dominant negative mechanism. In conclusion, this study elucidates the pathomechanism underlying NSF-related DEE and identifies a potential therapeutic approach.

Age related macular degeneration (AMD) is the most prevalent cause of blindness in the developed world. Vision loss in the advanced stages of the disease is caused by atrophy of retinal photoreceptors, overlying retinal pigment epithelium (RPE) and choroidal endothelial cells (CECs). The molecular events that underline the development of these cell types from in utero to adult as well as the progression to intermediate and advanced stages AMD are not yet fully understood. We performed single cell RNA-Sequencing (Seq) of human fetal and adult RPE-choroidal tissues, profiling in detail all the cell types and elucidating cell type specific proliferation, differentiation, and immunomodulation events that occur up to midgestation. Our data demonstrate that progression from the fetal to adult state is characterised by an increase in expression of genes involved in the oxidative stress response and detoxification from heavy metals, suggesting a better defence against oxidative stress in the adult RPE-choroid tissue. Single cell comparative transcriptional analysis between a patient with intermediate AMD and an unaffected subject revealed a reduction in the number of RPE and melanocytes in the macular region of the AMD patient. Together these findings may suggest a macular loss of RPE cells and melanocytes in the AMD patients, but given the complex processing of tissues required for single cell RNA-Seq which is prone to technical artefacts, these findings need to be validated by additional techniques in a larger number of AMD patients and controls.

DNA drug molecules are not only widely used in gene therapy, but also play an important role in controlling the electrical properties of molecular electronics. Covalent binding, groove binding and intercalation are all important forms of drug-DNA interaction. But its applications are limited due to a lack of understanding of the electron transport mechanisms after different drug-DNA interaction modes. Here, we used a combination of density functional theory calculations and nonequilibrium Green's function formulation with decoherence to study the effect of drug molecules on the charge transport property of DNA under three different binding modes. Conductance of DNA is found to decrease from 2.35E-5 G0 to 1.95E-6 G0 upon doxorubicin intercalation due to modifications of the density of states in the near-highest occupied molecular orbital region, δG = 1105.13%. Additionally, the conductance of DNA after cis-[Pt(NH3)2(py)Cl]+ covalent binding increases from 1.02E-6 G0 to 5.25E-5 G0, δG = 5047.06%. However, in the case of pentamidine groove binding, because there is no direct change in DNA molecular structure during drug binding, the conductance changes before and after drug binding is much smaller than in the two above cases, δG = 90.43%. Our theoretical calculations suggest that the conductance of DNA can be regulated by different drug molecules or switching the interaction modes between small molecules and DNA. This regulation opens new possibilities for their potential applications in controllable modulation of the electron transport property of DNA.

Nutritional influences have been discussed as potential modulators of Parkinson's disease (PD) pathology through various epidemiological and physiological studies. In animal models, a high-fat diet (HFD) with greater intake of lipid-derived calories leads to accelerated disease onset and progression. The underlying molecular mechanisms of HFD-induced aggravated pathology, however, remain largely unclear. In this study, we aimed to further illuminate the effects of a fat-enriched diet in PD by examining the brainstem and hippocampal transcriptome of alpha-synuclein transgenic mice exposed to a life-long HFD. Investigating individual transcript isoforms, differential gene expression and co-expression clusters, we observed that transcriptional differences between wild-type (WT) and transgenic animals intensified in both regions under HFD. Both brainstem and hippocampus displayed strikingly similar transcriptomic perturbation patterns. Interestingly, expression differences resulted mainly from responses in WT animals to HFD, while these genes remained largely unchanged or were even slightly oppositely regulated by diet in transgenic animals. Genes and co-expressed gene groups exhibiting this dysregulation were linked to metabolic and mitochondrial pathways. Our findings propose the failure of metabolic adaptions as the potential explanation for accelerated disease unfolding under exposure to HFD. From the identified clusters of co-expressed genes, several candidates lend themselves to further functional investigations.

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Atopic dermatitis is a chronically recurrent dermatologic disease affected by complex pathophysiology with limited therapeutic options. To identify promising biomarkers for atopic dermatitis, we conducted a Mendelian randomization (MR) study to systematically screen blood metabolome for potential causal mediators of atopic dermatitis and further predict target-mediated side effects. We selected 128 unique blood metabolites from 3 European-descent metabolome genome-wide association studies (GWASs) with a total of 147 827 participants. Atopic dermatitis dataset originated from a large-scale GWAS including 10 788 cases and 30 047 controls of European ancestry. MR analyses were performed to estimate the associations of blood metabolites with atopic dermatitis. We then applied a phenome-wide MR analysis to ascertain potential on-target side effects of metabolite intervention. Three metabolites were identified as potential causal mediators for atopic dermatitis, including docosahexaenoic acid (odds ratio [OR], 0.87; 95% CI, 0.81-0.94; P = 3.45 × 10-4), arachidonate (OR, 0.30; 95% CI, 0.17-0.53; P = 4.09 × 10-5), and 1-arachidonoylglycerophosphoethanolamine (1-arachidonoyl-GPE)(OR, 0.25; 95% CI, 0.12-0.53; P = 2.58 × 10-4). In the phenome-wide MR analysis, docosahexaenoic acid and arachidonate were also identified to have beneficial or detrimental effects on multiple diseases beyond atopic dermatitis, respectively. No adverse side effects were found for 1-arachidonoyl-GPE. In this systematic MR study, docosahexaenoic acid, arachidonate, and 1-arachidonoyl-GPE were identified as potential causal and beneficial mediators in the development of atopic dermatitis. Side-effect profiles were characterized to help inform drug target prioritization, and 1-arachidonoyl-GPE was a promising target for prevention and treatment of atopic dermatitis with no predicted adverse side effects.

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Congenital hypopituitarism is a genetically heterogeneous condition that is part of a spectrum disorder that can include holoprosencephaly. Heterozygous mutations in SIX3 cause variable holoprosencephaly in humans and mice. We identified two children with neonatal hypopituitarism and thin pituitary stalk who were doubly heterozygous for rare, likely deleterious variants in the transcription factors SIX3 and POU1F1. We used genetically engineered mice to understand the disease pathophysiology. Pou1f1 loss-of-function heterozygotes are unaffected; Six3 heterozygotes have pituitary gland dysmorphology and incompletely ossified palate; and the Six3+/-; Pou1f1+/dw double heterozygote mice have a pronounced phenotype, including pituitary growth through the palate. The interaction of Pou1f1 and Six3 in mice supports the possibility of digenic pituitary disease in children. Disruption of Six3 expression in the oral ectoderm completely ablated anterior pituitary development, and deletion of Six3 in the neural ectoderm blocked the development of the pituitary stalk and both anterior and posterior pituitary lobes. Six3 is required in both oral and neural ectodermal tissues for the activation of signaling pathways and transcription factors necessary for pituitary cell fate. These studies clarify the mechanism of SIX3 action in pituitary development and provide support for a digenic basis for hypopituitarism.

YWHAZ encodes an adapter protein 14-3-3ζ, which is involved in many signaling pathways that control cellular proliferation, migration and differentiation. It has not been definitely correlated to any phenotype in OMIM. To investigate the role of YWHAZ gene in intellectual disability and global developmental delay, we conducted whole-exon sequencing in all of the available members from a large three-generation family and we discovered that a novel variant of the YWHAZ gene was associated with intellectual disability and global developmental delay. This variant is a missense mutation of YWHAZ, p.Lys49Asn/c.147A > T, which was found in all affected members but not found in other unaffected members. We also conducted computational modeling and knockdown/knockin with Drosophila to confirm the role of the YWHAZ variant in intellectual disability. Computational modeling showed that the binding energy was increased in the mutated protein combining with the ligand indicating that the c147A > T variation was a loss-of-function variant. Cognitive defects and mushroom body morphological abnormalities were observed in YWHAZ c.147A > T knockin flies. The YWHAZ knockdown flies also manifested serious cognitive defects with hyperactivity behaviors, which is consistent with the clinical features. Our clinical and experimental results consistently suggested that YWHAZ was a novel intellectual disability pathogenic gene.

Kinesins are motor proteins involved in microtubule (MT)-mediated intracellular transport. They contribute to key cellular processes, including intracellular trafficking, organelle dynamics and cell division. Pathogenic variants in kinesin-encoding genes underlie several human diseases characterized by an extremely variable clinical phenotype, ranging from isolated neurodevelopmental/neurodegenerative disorders to syndromic phenotypes belonging to a family of conditions collectively termed as 'ciliopathies.' Among kinesins, kinesin-1 is the most abundant MT motor for transport of cargoes towards the plus end of MTs. Three kinesin-1 heavy chain isoforms exist in mammals. Different from KIF5A and KIF5C, which are specifically expressed in neurons and established to cause neurological diseases when mutated, KIF5B is an ubiquitous protein. Three de novo missense KIF5B variants were recently described in four subjects with a syndromic skeletal disorder characterized by kyphomelic dysplasia, hypotonia and DD/ID. Here, we report three dominantly acting KIF5B variants (p.Asn255del, p.Leu498Pro and p.Leu537Pro) resulting in a clinically wide phenotypic spectrum, ranging from dilated cardiomyopathy with adult-onset ophthalmoplegia and progressive skeletal myopathy to a neurodevelopmental condition characterized by severe hypotonia with or without seizures. In vitro and in vivo analyses provide evidence that the identified disease-associated KIF5B variants disrupt lysosomal, autophagosome and mitochondrial organization, and impact cilium biogenesis. All variants, and one of the previously reported missense changes, were shown to affect multiple developmental processes in zebrafish. These findings document pleiotropic consequences of aberrant KIF5B function on development and cell homeostasis, and expand the phenotypic spectrum resulting from altered kinesin-mediated processes.

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Usher syndrome (USH) is the most common form of hereditary deaf-blindness in humans. USH is a complex genetic disorder, assigned to three clinical subtypes differing in onset, course and severity, with USH1 being the most severe. Rodent USH1 models do not reflect the ocular phenotype observed in human patients to date; hence, little is known about the pathophysiology of USH1 in the human eye. One of the USH1 genes, USH1C, exhibits extensive alternative splicing and encodes numerous harmonin protein isoforms that function as scaffolds for organizing the USH interactome. RNA-seq analysis of human retinae uncovered harmonin_a1 as the most abundant transcript of USH1C. Bulk RNA-seq analysis and immunoblotting showed abundant expression of harmonin in Müller glia cells (MGCs) and retinal neurons. Furthermore, harmonin was localized in the terminal endfeet and apical microvilli of MGCs, presynaptic region (pedicle) of cones and outer segments (OS) of rods as well as at adhesive junctions between MGCs and photoreceptor cells (PRCs) in the outer limiting membrane (OLM). Our data provide evidence for the interaction of harmonin with OLM molecules in PRCs and MGCs and rhodopsin in PRCs. Subcellular expression and colocalization of harmonin correlate with the clinical phenotype observed in USH1C patients. We also demonstrate that primary cilia defects in USH1C patient-derived fibroblasts could be reverted by the delivery of harmonin_a1 transcript isoform. Our studies thus provide novel insights into PRC cell biology, USH1C pathophysiology and development of gene therapy treatment(s).

INTRODUCTION: Since 2005, disease-related human genetic diversity has been intensively characterized using genome-wide association studies (GWAS). Understanding how and by whom this work was performed may yield valuable insights into the generalizability of GWAS discoveries to global populations and how high-impact genetics research can be equitably sustained in the future. MATERIALS AND METHODS: We mined the NHGRI-EBI GWAS Catalog (2005-2022) for the most burdensome non-communicable causes of death worldwide. We then compared (i) the geographic, ethnic and socioeconomic characteristics of study populations; (ii) the geographic and socioeconomic characteristics of the regions within which researchers were located and (iii) the extent to which male and female investigators undertook and led the research. RESULTS: The research institutions leading the work are often US-based (37%), while the origin of samples is more diverse, with the Nordic countries having contributed as much data to GWAS as the United States (~17% of data). The majority of first (60%), senior (75%) and all (66%) authors are male; although proportions vary by disease and leadership level, male co-authors are the ubiquitous majority. The vast majority (91%) of complex trait GWAS has been performed in European ancestry populations, with cohorts and scientists predominantly located in medium-to-high socioeconomically ranked countries; apart from East Asians (~5%), other ethnicities rarely feature in published GWAS. See: https://hugofitipaldi.shinyapps.io/gwas_results/ to browse all results. CONCLUSION: Most GWAS cohorts are of European ancestry residing outside the United States, with a smaller yet meaningful proportion of East Asian ancestry. Papers describing GWAS research are predominantly authored by male scientists based in medium-to-high income countries.

Folic acid (synthetic folate, FA) is consumed in excess in North America and may interact with common pathogenic variants in methylenetetrahydrofolate reductase (MTHFR); the most prevalent inborn error of folate metabolism with wide ranging obesity-related comorbidities. While preclinical murine models have been valuable to inform on diet-gene interactions, a recent Folate Expert panel has encouraged validation of new animal models. In this study, we characterized a novel zebrafish model of mthfr deficiency and evaluated the effects of genetic loss of mthfr function and FA supplementation during embryonic development on energy homeostasis and metabolism. Mthfr-deficient zebrafish were generated using CRISPR mutagenesis and supplemented with no FA (control, 0FA) or 100 μM FA (100FA) throughout embryonic development (0-5 days post-fertilization). We show that genetic loss of mthfr function in zebrafish recapitulates key biochemical hallmarks reported in MTHFR deficiency in humans and leads to greater lipid accumulation and aberrant cholesterol metabolism as reported in the Mthfr murine model. In mthfr-deficient zebrafish, energy homeostasis was also impaired as indicated by altered food intake, reduced metabolic rate, and lower expression of central energy-regulatory genes. Microglia abundance, involved in healthy neuronal development, was also reduced. FA supplementation to control zebrafish mimicked many of the adverse effects of mthfr deficiency, some of which were also exacerbated in mthfr-deficient zebrafish. Together, these findings support the translatability of the mthfr-deficient zebrafish as a preclinical model in folate research.

Mucopolysaccharidosis type IIIB (MPS IIIB) is an autosomal recessive lysosomal storage disease caused by mutations in the gene that encodes the protein N-acetyl-glucosaminidase (NAGLU). Defective NAGLU activity results in aberrant retention of heparan sulfate within lysosomes leading to progressive central nervous system (CNS) degeneration. Intravenous treatment options are limited by the need to overcome the blood-brain barrier and gain successful entry into the CNS. Additionally, we have demonstrated that AAV8 provides a broader transduction area in the MPS IIIB mouse brain compared with AAV5, 9 or rh10. A triple-capsid mutant (tcm) modification of AAV8 further enhanced GFP reporter expression and distribution. Using the MPS IIIB mouse model, we performed a study using either intracranial six site or intracisterna magna injection of AAVtcm8-codon-optimized (co)-NAGLU using untreated MPS IIIB mice as controls to assess disease correction. Disease correction was evaluated based on enzyme activity, heparan sulfate storage levels, CNS lysosomal signal intensity, coordination, activity level, hearing and survival. Both histologic and enzymatic assessments show that each injection method results in supranormal levels of NAGLU expression in the brain. In this study, we have shown correction of lifespan and auditory deficits, increased CNS NAGLU activity and reduced lysosomal storage levels of heparan sulfate following AAVtcm8-coNAGLU administration and partial correction of NAGLU activity in several peripheral organs in the murine model of MPS IIIB.

Accumulating toxic protein assemblies, including Aβ and tau, and dysfunctional mitochondria are associated with synaptic and neuronal loss in Alzheimer's disease (ad). Such accumulations are thought to be due to clearance defects in the autophagy-lysosome pathway. Mitochondrial dysfunction is evident in ad brains and animal models at multiple levels, such as mitochondrial genomic mutations, disrupted bioenergetics, deregulated mitochondrial dynamics, and impaired clearance of damaged mitochondria (mitophagy). Slingshot Homolog-1 (SSH1) is a phosphatase activated by oxidative stress, high intracellular levels of Ca+2, and Aβ42 oligomers (Aβ42O), known for its function to dephosphorylate/activate cofilin through the N-terminal region. SSH1-mediated cofilin dephosphorylation results in Ab42O-induced severing of F-actin and translocation of cofilin to mitochondria, which promotes mitochondria-mediated apoptosis, synaptic loss, and synaptic deficits. On the other hand, SSH1-mediated dephosphorylation/deactivation of the autophagy-cargo receptor p62 (SQSTM1), through its C-terminal region, inhibits p62 autophagy flux. However, the interplay between these two different activities of SSH1 in Aβ42O-induced mitochondrial toxicity remains unclear. In this study, we assessed the role of endogenous SSH1 and different regions of SSH1 in regulating mitochondrial health, mitochondrial respiration, clearance of damaged mitochondria, and synaptic integrity in vitro and in vivo. Our results indicate that SSH1 suppresses mitochondrial health and respiration through the cofilin-binding N-terminal region, whereas SSH1 impairs mitophagy through a newly identified ~ 100 residue p62-binding domain in the C-terminal region. These results indicate that both N-terminal and C-terminal regions negatively impact mitochondria by distinct and independent modalities to amplify mitochondrial abnormalities, making SSH1 an excellent target to mitigate ad pathogenesis.

De novo deleterious and heritable biallelic mutations in the DNA binding domain (DBD) of the transcription factor deformed epidermal autoregulatory factor 1 (DEAF1) result in a phenotypic spectrum of disorders termed DEAF1-associated neurodevelopmental disorders (DAND). RNA-sequencing using hippocampal RNA from mice with conditional deletion of Deaf1 in the central nervous system indicate that loss of Deaf1 activity results in the altered expression of genes involved in neuronal function, dendritic spine maintenance, development, and activity, with reduced dendritic spines in hippocampal regions. Since DEAF1 is not a dosage-sensitive gene, we assessed the dominant negative activity of previously identified de novo variants and a heritable recessive DEAF1 variant on selected DEAF1-regulated genes in 2 different cell models. While no altered gene expression was observed in cells over-expressing the recessive heritable variant, the gene expression profiles of cells over-expressing de novo variants resulted in similar gene expression changes as observed in CRISPR-Cas9-mediated DEAF1-deleted cells. Altered expression of DEAF1-regulated genes was rescued by exogenous expression of WT-DEAF1 but not by de novo variants in cells lacking endogenous DEAF1. De novo heterozygous variants within the DBD of DEAF1 were identified in 10 individuals with a phenotypic spectrum including autism spectrum disorder, developmental delays, sleep disturbance, high pain tolerance, and mild dysmorphic features. Functional assays demonstrate these variants alter DEAF1 transcriptional activity. Taken together, this study expands the clinical phenotypic spectrum of individuals with DAND, furthers our understanding of potential roles of DEAF1 on neuronal function, and demonstrates dominant negative activity of identified de novo variants.

Epilepsy is a chronic neurological disorder featuring recurrent, unprovoked seizures, which affect more than 65 million people worldwide. Here, we discover that the PKHD1L1, which is encoded by polycystic kidney and hepatic disease1-like 1 (Pkhd1l1), wildly distributes in neurons in the central nervous system (CNS) of mice. Disruption of PKHD1L1 in the dentate gyrus region of the hippocampus leads to increased susceptibility to pentylenetetrazol-induced seizures in mice. The disturbance of PKHD1L1 leads to the overactivation of the mitogen-activated protein kinase (MAPK)/extracellular regulated kinase (ERK)-Calpain pathway, which is accompanied by remarkable degradation of cytoplasmic potassium chloride co-transporter 2 (KCC2) level together with the impaired expression and function of membrane KCC2. However, the reduction of membrane KCC2 is associated with the damaged inhibitory ability of the vital GABA receptors, which ultimately leads to the significantly increased susceptibility to epileptic seizures. Our data, thus, indicate for the first time that Pkhd1l1, a newly discovered polycystic kidney disease (PKD) association gene, is required in neurons to maintain neuronal excitability by regulation of KCC2 expression in CNS. A new mechanism of the clinical association between genetic PKD and seizures has been built, which could be a potential therapeutic target for treating PKD-related seizures.

Non-small cell lung cancer (NSCLC) is one of the most malignant epithelial tumors. Studies have suggested that DNA hypermethylation of promoters and abnormal histone modifications could induce tumor suppressor genes (TSGs) downregulation in NSCLC. However, the exact mechanism of TSGs downregulation remains unclear. In this study, we found that there is no difference in the regions of most TSGs promoters in NSCLC. Moreover, we found that there is no DNA methylation difference in the region of VILL promoter in NSCLC compared with adjacent tissue samples by pyrosequencing. We further demonstrated that VILL was markedly reactivated in A549 and H1703 cells infected with miR-26A1 lentivirus while this activation was inhibited by JQ1, an enhancer inhibitor. In addition, we identified that miR-26A1 could function as a tumor suppressor to inhibit proliferation and metastasis of NSCLC cells. Chromatin immunoprecipitation assays revealed that overexpression of miR-26A1 could significantly induce the enrichment of H3K27ac at the enhancer regions in A549 cells. To sum up, our findings revealed that enhancer-mediated TSGs regulation occured in NSCLC, suggesting that miR-26A1 could serve as a key regulator and may be a potential therapeutic target for NSCLC.

Pathogenic variants in ATL1 are a known cause of autosomal-dominantly inherited hereditary spastic paraplegia (HSP-ATL1, SPG3A) with a predominantly 'pure' HSP phenotype. Although a relatively large number of patients have been reported, no genotype-phenotype correlations have been established for specific ATL1 variants. Confronted with five children carrying de novo ATL1 variants showing early, complex and severe symptoms, we systematically investigated the molecular and phenotypic spectrum of HSP-ATL1. Through a cross-sectional analysis of 537 published and novel cases, we delineate a distinct phenotype observed in patients with de novo variants. Guided by this systematic phenotyping approach and structural modelling of disease-associated variants in atlastin-1, we demonstrate that this distinct phenotypic signature is also prevalent in a subgroup of patients with inherited ATL1 variants and is largely explained by variant localization within a three-dimensional mutational cluster. Establishing genotype-phenotype correlations, we find that symptoms that extend well beyond the typical pure HSP phenotype (i.e. neurodevelopmental abnormalities, upper limb spasticity, bulbar symptoms, peripheral neuropathy and brain imaging abnormalities) are prevalent in patients with variants located within this mutational cluster.

One of the most significant risk variants for Parkinson's disease (PD), rs356182, is located at the PD-associated locus near the alpha-synuclein (α-syn) encoding gene, SNCA. SNCA-proximal variants, including rs356182, are thought to function in PD risk through enhancers via allele-specific regulatory effects on SNCA expression. However, this interpretation discounts the complex activity of genetic enhancers and possible non-conical functions of α-syn. Here we investigated a novel risk mechanism for rs356182. We use CRISPR-Cas9 in LUHMES cells, a model for dopaminergic midbrain neurons, to generate precise hemizygous lesions at rs356182. The PD-protective (A/-), PD-risk (G/-) and wild-type (A/G) clones were neuronally differentiated and then compared transcriptionally and morphologically. Among the affected genes was SNCA, whose expression was promoted by the PD-protective allele (A) and repressed in its absence. In addition to SNCA, hundreds of genes were differentially expressed and associated with neurogenesis and axonogenesis-an effect not typically ascribed to α-syn. We also found that the transcription factor FOXO3 specifically binds to the rs356182 A-allele in differentiated LUHMES cells. Finally, we compared the results from the rs356182-edited cells to our previously published knockouts of SNCA and found only minimal overlap between the sets of significant differentially expressed genes. Together, the data implicate a risk mechanism for rs356182 in which the risk-allele (G) is associated with abnormal neuron development, independent of SNCA expression. We speculate that these pathological effects manifest as a diminished population of dopaminergic neurons during development leading to the predisposition for PD later in life.

Genetic variation in genes regulating metabolism may be advantageous in some settings but not others. The non-failing adult heart relies heavily on fatty acids as a fuel substrate and source of ATP. In contrast, the failing heart favors glucose as a fuel source. A bootstrap analysis for genes with deviant allele frequencies in cardiomyopathy cases versus controls identified the MTCH2 gene as having unusual variation. MTCH2 encodes an outer mitochondrial membrane protein, and prior genome-wide studies associated MTCH2 variants with body mass index, consistent with its role in metabolism. We identified the referent allele of rs1064608 (p.Pro290) as being overrepresented in cardiomyopathy cases compared to controls, and linkage disequilibrium analysis associated this variant with the MTCH2 cis eQTL rs10838738 and lower MTCH2 expression. To evaluate MTCH2, we knocked down Mtch in Drosophila heart tubes which produced a dilated and poorly functioning heart tube, reduced adiposity and shortened life span. Cardiac Mtch mutants generated more lactate at baseline, and they displayed impaired oxygen consumption in the presence of glucose but not palmitate. Treatment of cardiac Mtch mutants with dichloroacetate, a pyruvate dehydrogenase kinase inhibitor, reduced lactate and rescued lifespan. Deletion of MTCH2 in human cells similarly impaired oxygen consumption in the presence of glucose but not fatty acids. These data support a model in which MTCH2 reduction may be favorable when fatty acids are the major fuel source, favoring lean body mass. However, in settings like heart failure, where the heart shifts toward using more glucose, reduction of MTCH2 is maladaptive.

Author information: (1)Institute for Molecular Medicine Finland, FIMM, HiLIFE, University of Helsinki, Helsinki, Finland. (2)Department of Pulmonology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland. (3)Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA. (4)Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, USA.

Bietti crystalline corneoretinal dystrophy (BCD) is an autosomal recessive retinal degenerative disease characterized by yellow-white crystal deposits in the posterior pole, degeneration of the retinal pigment epithelium (RPE), and sclerosis of the choroid. Mutations in the cytochrome P450 4V2 gene (CYP4V2) cause BCD, which is associated with lipid metabolic disruption. The use of gene-replacement therapy in BCD has been hampered by the lack of disease models. To advance CYP4V2 gene-replacement therapy, we generated BCD patient-specific induced pluripotent stem cell (iPSC)-RPE cells and Cyp4v3 knockout (KO) mice as disease models and AAV2/8-CAG-CYP4V2 as treatment vectors. We demonstrated that after adeno-associated virus (AAV)-mediated CYP4V2 gene-replacement therapy BCD-iPSC-RPE cells presented restored cell survival and reduced lipid droplets accumulation; restoration of vision in Cyp4v3 KO mice was revealed by elevated electroretinogram amplitude and ameliorated RPE degeneration. These results suggest that AAV-mediated gene-replacement therapy in BCD patients is a promising strategy.

Filamin A (FLNA) is a cytoplasmic actin binding protein, recently shown to be expressed as a long and short isoform. Mutations in FLNA are associated with a wide spectrum of disorders, including an X-linked form of chronic intestinal pseudo-obstruction (CIPO). However, the role of FLNA in intestinal development and function is largely unknown. In this study, we show that FLNA is expressed in the muscle layer of the small intestine from early human fetal stages. Expression of FLNA variants associated with CIPO, blocked expression of the long flna isoform and led to an overall reduction of RNA and protein levels. As a consequence, contractility of human intestinal smooth muscle cells was affected. Lastly, our transgenic zebrafish line showed that the flna long isoform is required for intestinal elongation and peristalsis. Histological analysis revealed structural and architectural changes in the intestinal smooth muscle of homozygous fish, likely triggered by the abnormal expression of intestinal smooth muscle markers. No defect in the localization or numbers of enteric neurons was observed. Taken together, our study demonstrates that the long FLNA isoform contributes to intestinal development and function. Since loss of the long FLNA isoform does not seem to affect the enteric nervous system, it likely results in a myopathic form of CIPO, bringing new insights to disease pathogenesis.

Sox9 plays an essential role in mammalian testis formation. It has been reported that gene expression in the testes is regulated by enhancers. Among them, mXYSRa/Enh13-which is located at far upstream of the transcription start site-plays a critical role, wherein its deletion causes complete male-to-female sex reversal in mice. It has been proposed that the binding sites (BSs) of SOX9 and SRY, the latter of which is the sex determining gene on the Y chromosome, are associated with mXYSRa/Enh13. They function as an enhancer, whereby the sequences are evolutionarily conserved and in vivo binding of SOX9 and SRY to mXYSRa/Enh13 has been demonstrated previously. However, their precise in vivo functions have not been examined to date. To this end, this study generated mice with substitutions on the SOX9 and SRY BSs to reveal their in vivo functions. Homozygous mutants of SOX9 and SRY BS were indistinguishable from XY males, whereas double mutants had small testes, suggesting that these functions are redundant and that there is another functional sequence on mXYSRa/Enh13, since mXYSRa/Enh13 deletion mice are XY females. In addition, the majority of hemizygous mice with substitutions in SOX9 BS and SRY BS were female and male, respectively, suggesting that SOX9 BS contributes more to SRY BS for mXYSRa/Enh13 to function. The additive effect of SOX9 and SRY via these BSs was verified using an in vitro assay. In conclusion, SOX9 BS and SRY BS function redundantly in vivo, and at least one more functional sequence should exist in mXYSRa/Enh13.

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1-P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10-5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10-8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

Eukaryotic RNA polymerase I (Pol I) products play fundamental roles in ribosomal assembly, protein synthesis, metabolism and cell growth. Abnormal expression of both Pol I transcription-related factors and Pol I products causes a range of diseases, including ribosomopathies and cancers. However, the factors and mechanisms governing Pol I-dependent transcription remain to be elucidated. Here, we report that transcription factor IIB-related factor 1 (BRF1), a subunit of transcription factor IIIB required for RNA polymerase III (Pol III)-mediated transcription, is a nucleolar protein and modulates Pol I-mediated transcription. We showed that BRF1 can be localized to the nucleolus in several human cell types. BRF1 expression correlates positively with Pol I product levels and tumour cell growth in vitro and in vivo. Pol III transcription inhibition assays confirmed that BRF1 modulates Pol I-directed transcription in an independent manner rather than through a Pol III product-to-45S pre-rRNA feedback mode. Mechanistically, BRF1 binds to the Pol I transcription machinery components and can be recruited to the rDNA promoter along with them. Additionally, alteration of BRF1 expression affects the recruitment of Pol I transcription machinery components to the rDNA promoter and the expression of TBP and TAF1A. These findings indicate that BRF1 modulates Pol I-directed transcription by controlling the expression of selective factor 1 subunits. In summary, we identified a novel role of BRF1 in Pol I-directed transcription, suggesting that BRF1 can independently regulate both Pol I- and Pol III-mediated transcription and act as a key coordinator of Pol I and Pol III.

Glutaminase deficiency has recently been associated with ataxia and developmental delay due to repeat expansions in the 5'UTR of the glutaminase (GLS) gene. Patients with the described GLS repeat expansion may indeed remain undiagnosed due to the rarity of this variant, the challenge of its detection and the recency of its discovery. In this study, we combined advanced bioinformatics screening of ~3000 genomes and ~1500 exomes with optical genome mapping and long-read sequencing for confirmation studies. We identified two GLS families, previously intensely and unsuccessfully analyzed. One family carries an unusual and complex structural change involving a homozygous repeat expansion nested within a quadruplication event in the 5'UTR of GLS. Glutaminase deficiency and its metabolic consequences were validated by in-depth biochemical analysis. The identified GLS patients showed progressive early-onset ataxia, cognitive deficits, pyramidal tract damage and optic atrophy, thus demonstrating susceptibility of several specific neuron populations to glutaminase deficiency. This large-scale screening study demonstrates the ability of bioinformatics analysis-validated by latest state-of-the-art technologies (optical genome mapping and long-read sequencing)-to effectively flag complex repeat expansions using short-read datasets and thus facilitate diagnosis of ultra-rare disorders.

Women of African ancestry have the highest mortality from triple-negative breast cancer (TNBC) of all racial groups. To understand the genomic basis of breast cancer in the populations, we previously conducted genome-wide association studies and identified single nucleotide polymorphisms (SNPs) associated with breast cancer in Black women. In this study, we investigated the functional significance of the top associated SNP rs13074711. We found the SNP served as an enhancer variant and regulated TNFSF10 (TRAIL) expression in TNBC cells, with a significant association between the SNP genotype and TNFSF10 expression in breast tumors. Mechanistically, rs13074711 modulated the binding activity of c-MYB at the motif and thereby controlled TNFSF10 expression. Interestingly, TNFSF10 expression in many cancers was consistently lower in African Americans compared with European Americans. Furthermore, TNFSF10 expression in TNBC was significantly correlated with the expression of antiviral immune genes and was regulated by type I interferons (IFNs). Accordingly, loss of TNFSF10 resulted in a profound decrease in apoptosis of TNBC cells in response to type I IFNs and poly(I:C), a synthetic analogue of double stranded virus. Lastly, in a syngeneic mouse model of breast cancer, TNFSF10-deficiency in breast tumors decreased tumor-infiltrated CD4+ and CD8+ T cell quantities. Collectively, our results suggested that TNFSF10 plays an important role in the regulation of antiviral immune responses in TNBC, and the expression is in part regulated by a genetic variant associated with breast cancer in Black women. Our results underscore the important contributions of genetic variants to immune defense mechanisms.

Huntington's disease (HD) is a neurodegenerative disorder caused by an inherited unstable HTT CAG repeat that expands further, thereby eliciting a disease process that may be initiated by polyglutamine-expanded huntingtin or a short polyglutamine-product. Phosphorylation of selected candidate residues is reported to mediate polyglutamine-fragment degradation and toxicity. Here to support the discovery of phosphosites involved in the life-cycle of (full-length) huntingtin, we employed mass spectrometry-based phosphoproteomics to systematically identify sites in purified huntingtin and in the endogenous protein by proteomic and phosphoproteomic analyses of members of an HD neuronal progenitor cell panel. Our results bring total huntingtin phosphosites to 95, with more located in the N-HEAT domain relative to numbers in the Bridge and C-HEAT domains. Moreover, phosphorylation of C-HEAT Ser2550 by cAMP-dependent protein kinase (PKA), the top hit in kinase activity screens, was found to hasten huntingtin degradation, such that levels of the catalytic subunit (PRKACA) were inversely related to huntingtin levels. Taken together, these findings highlight categories of phosphosites that merit further study and provide a phosphosite kinase pair (pSer2550-PKA) with which to investigate the biological processes that regulate huntingtin degradation and thereby influence the steady state levels of huntingtin in HD cells.