Search

Human cell line models, including the neuronal precursor line LUHMES, are important for investigating developmental transcriptional dynamics within imprinted regions, particularly the 15q11-q13 Angelman (AS) and Prader-Willi (PWS) syndrome locus. AS results from loss of maternal UBE3A in neurons, where the paternal allele is silenced by a convergent antisense transcript UBE3A-ATS, a lncRNA that terminates at PWAR1 in non-neurons. qRT-PCR analysis confirmed the exclusive and progressive increase in UBE3A-ATS in differentiating LUHMES neurons, validating their use for studying UBE3A silencing. Genome-wide transcriptome analyses revealed changes to 11 834 genes during neuronal differentiation, including the upregulation of most genes within the 15q11-q13 locus. To identify dynamic changes in chromatin loops linked to transcriptional activity, we performed a HiChIP validated by 4C, which identified two neuron-specific CTCF loops between MAGEL2-SNRPN and PWAR1-UBE3A. To determine if allele-specific differentially methylated regions (DMR) may be associated with CTCF loop anchors, whole genome long-read nanopore sequencing was performed. We identified a paternally hypomethylated DMR near the SNRPN upstream loop anchor exclusive to neurons and a paternally hypermethylated DMR near the PWAR1 CTCF anchor exclusive to undifferentiated cells, consistent with increases in neuronal transcription. Additionally, DMRs near CTCF loop anchors were observed in both cell types, indicative of allele-specific differences in chromatin loops regulating imprinted transcription. These results provide an integrated view of the 15q11-q13 epigenetic landscape during LUHMES neuronal differentiation, underscoring the complex interplay of transcription, chromatin looping, and DNA methylation. They also provide insights for future therapeutic approaches for AS and PWS.

Duchenne Muscular Dystrophy (DMD) is a progressive and fatal neuromuscular disease. Cycles of myofibre degeneration and regeneration are hallmarks of the disease where immune cells infiltrate to repair damaged skeletal muscle. Benfotiamine is a lipid soluble precursor to thiamine, shown clinically to reduce inflammation in diabetic related complications. We assessed whether benfotiamine administration could reduce inflammation related dystrophic pathology. Benfotiamine (10 mg/kg/day) was fed to male mdx mice (n = 7) for 15 weeks from 4 weeks of age. Treated mice had an increased growth weight (5-7 weeks) and myofibre size at treatment completion. Markers of dystrophic pathology (area of damaged necrotic tissue, central nuclei) were reduced in benfotiamine mdx quadriceps. Grip strength was increased and improved exercise capacity was found in mdx treated with benfotiamine for 12 weeks, before being placed into individual cages and allowed access to an exercise wheel for 3 weeks. Global gene expression profiling (RNAseq) in the gastrocnemius revealed benfotiamine regulated signalling pathways relevant to dystrophic pathology (Inflammatory Response, Myogenesis) and fibrotic gene markers (Col1a1, Col1a2, Col4a5, Col5a2, Col6a2, Col6a2, Col6a3, Lum) towards wildtype levels. In addition, we observed a reduction in gene expression of inflammatory gene markers in the quadriceps (Emr1, Cd163, Cd4, Cd8, Ifng). Overall, these data suggest that benfotiamine reduces dystrophic pathology by acting on inflammatory and fibrotic gene markers and signalling pathways. Given benfotiamine's excellent safety profile and current clinical use, it could be used in combination with glucocorticoids to treat DMD patients.

Craniosynostosis, defined by premature fusion of one or multiple cranial sutures, is a common congenital defect affecting more than 1/2000 infants and results in restricted brain expansion. Single gene mutations account for 15%-20% of cases, largely as part of a syndrome, but the majority are nonsyndromic with complex underlying genetics. We hypothesized that the two noncoding genomic regions identified by a GWAS for craniosynostosis contain distal regulatory elements for the risk genes BMPER and BMP2. To identify such regulatory elements, we surveyed conserved noncoding sequences from both risk loci for enhancer activity in transgenic Danio rerio. We identified enhancers from both regions that direct expression to skeletal tissues, consistent with the endogenous expression of bmper and bmp2. For each locus, we also found a skeletal enhancer that also contains a sequence variant associated with craniosynostosis risk. We examined the activity of each enhancer during craniofacial development and found that the BMPER-associated enhancer is active in the restricted region of cartilage closely associated with frontal bone initiation. The same enhancer is active in mouse skeletal tissues, demonstrating evolutionarily conserved activity. Using enhanced yeast one-hybrid assays, we identified transcription factors that bind each enhancer and observed differential binding between alleles, implicating multiple signaling pathways. Our findings help unveil the genetic mechanism of the two craniosynostosis risk loci. More broadly, our combined in vivo approach is applicable to many complex genetic diseases to build a link between association studies and specific genetic mechanisms.

Astrocytes and brain endothelial cells are components of the neurovascular unit that comprises the blood-brain barrier (BBB) and their dysfunction contributes to pathogenesis in Huntington's disease (HD). Defining the contribution of these cells to disease can inform cell-type-specific effects and uncover new disease-modifying therapeutic targets. These cells express integrin (ITG) adhesion receptors that anchor the cells to the extracellular matrix (ECM) to maintain the integrity of the BBB. We used HD patient-derived induced pluripotent stem cell (iPSC) modeling to study the ECM-ITG interface in astrocytes and brain microvascular endothelial cells and found ECM-ITG dysregulation in human iPSC-derived cells that may contribute to the dysfunction of the BBB in HD. This disruption has functional consequences since reducing ITG expression in glia in an HD Drosophila model suppressed disease-associated CNS dysfunction. Since ITGs can be targeted therapeutically and manipulating ITG signaling prevents neurodegeneration in other diseases, defining the role of ITGs in HD may provide a novel strategy of intervention to slow CNS pathophysiology to treat HD.

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.

NADK2 encodes the mitochondrial form of nicotinamide adenine dinucleotide (NAD) kinase, which phosphorylates NAD. Rare recessive mutations in human NADK2 are associated with a syndromic neurological mitochondrial disease that includes metabolic changes, such as hyperlysinemia and 2,4 dienoyl CoA reductase (DECR) deficiency. However, the full pathophysiology resulting from NADK2 deficiency is not known. Here, we describe two chemically induced mouse mutations in Nadk2-S326L and S330P-which cause severe neuromuscular disease and shorten lifespan. The S330P allele was characterized in detail and shown to have marked denervation of neuromuscular junctions by 5 weeks of age and muscle atrophy by 11 weeks of age. Cerebellar Purkinje cells also showed progressive degeneration in this model. Transcriptome profiling on brain and muscle was performed at early and late disease stages. In addition, metabolomic profiling was performed on the brain, muscle, liver and spinal cord at the same ages and on plasma at 5 weeks. Combined transcriptomic and metabolomic analyses identified hyperlysinemia, DECR deficiency and generalized metabolic dysfunction in Nadk2 mutant mice, indicating relevance to the human disease. We compared findings from the Nadk model to equivalent RNA sequencing and metabolomic datasets from a mouse model of infantile neuroaxonal dystrophy, caused by recessive mutations in Pla2g6. This enabled us to identify disrupted biological processes that are common between these mouse models of neurological disease, as well as those processes that are gene-specific. These findings improve our understanding of the pathophysiology of neuromuscular diseases and describe mouse models that will be useful for future preclinical studies.

RationaleGenetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear.ObjectivesWe sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.MethodsUsing the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants.Measurements and main resultsIn European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value

Coiled-coil-helix-coiled-coil-helix domain containing 2 (CHCHD2) is a mitochondrial protein that plays important roles in cristae structure, oxidative phosphorylation and apoptosis. Multiple mutations in CHCHD2 have been associated with Lewy body disorders (LBDs), such as Parkinson's disease (PD) and dementia with Lewy bodies, with the CHCHD2-T61I mutation being the most widely studied. However, at present, only CHCHD2 knockout or CHCHD2/CHCHD10 double knockout mouse models have been investigated. They do not recapitulate the pathology seen in patients with CHCHD2 mutations. We generated the first transgenic mouse model expressing the human PD-linked CHCHD2-T61I mutation driven by the mPrP promoter. We show that CHCHD2-T61I Tg mice exhibit perinuclear mitochondrial aggregates, neuroinflammation, and have impaired long-term synaptic plasticity associated with synaptic dysfunction. Dopaminergic neurodegeneration, a hallmark of PD, is also observed along with α-synuclein pathology. Significant motor dysfunction is seen with no changes in learning and memory at 1 year of age. A minor proportion of the CHCHD2-T61I Tg mice (~10%) show a severe motor phenotype consistent with human Pisa Syndrome, an atypical PD phenotype. Unbiased proteomics analysis reveals surprising increases in many insoluble proteins predominantly originating from mitochondria and perturbing multiple canonical biological pathways as assessed by ingenuity pathway analysis, including neurodegenerative disease-associated proteins such as tau, cofilin, SOD1 and DJ-1. Overall, CHCHD2-T61I Tg mice exhibit pathological and motor changes associated with LBDs, indicating that this model successfully captures phenotypes seen in human LBD patients with CHCHD2 mutations and demonstrates changes in neurodegenerative disease-associated proteins, which delineates relevant pathological pathways for further investigation.

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

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

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

Many neurodevelopmental disorders (NDDs) are the result of mutations on the X chromosome. One severe NDD resulting from mutations on the X chromosome is CDKL5 deficiency disorder (CDD). CDD is an epigenetic, X-linked NDD characterized by intellectual disability (ID), pervasive seizures and severe sleep disruption, including recurring hospitalizations. CDD occurs at a 4:1 ratio, with a female bias. CDD is driven by the loss of cyclin-dependent kinase-like 5 (CDKL5), a serine/threonine kinase that is essential for typical brain development, synapse formation and signal transmission. Previous studies focused on male subjects from animal models, likely to avoid the complexity of X mosaicism. For the first time, we report translationally relevant behavioral phenotypes in young adult (8-20 weeks) females and males with robust signal size, including impairments in learning and memory, substantial hyperactivity and increased susceptibility to seizures/reduced seizure thresholds, in both sexes, and in two models of CDD preclinical mice, one with a general loss-of-function mutation and one that is a patient-derived mutation.

Genetic variants in SCN2A, encoding the NaV1.2 voltage-gated sodium channel, are associated with a range of neurodevelopmental disorders with overlapping phenotypes. Some variants fit into a framework wherein gain-of-function missense variants that increase neuronal excitability lead to developmental and epileptic encephalopathy, while loss-of-function variants that reduce neuronal excitability lead to intellectual disability and/or autism spectrum disorder (ASD) with or without co-morbid seizures. One unique case less easily classified using this framework is the de novo missense variant SCN2A-p.K1422E, associated with infant-onset developmental delay, infantile spasms and features of ASD. Prior structure-function studies demonstrated that K1422E substitution alters ion selectivity of NaV1.2, conferring Ca2+ permeability, lowering overall conductance and conferring resistance to tetrodotoxin (TTX). Based on heterologous expression of K1422E, we developed a compartmental neuron model incorporating variant channels that predicted reductions in peak action potential (AP) speed. We generated Scn2aK1422E mice and characterized effects on neurons and neurological/neurobehavioral phenotypes. Cultured cortical neurons from heterozygous Scn2aK1422E/+ mice exhibited lower current density with a TTX-resistant component and reversal potential consistent with mixed ion permeation. Recordings from Scn2aK1442E/+ cortical slices demonstrated impaired AP initiation and larger Ca2+ transients at the axon initial segment during the rising phase of the AP, suggesting complex effects on channel function. Scn2aK1422E/+ mice exhibited rare spontaneous seizures, interictal electroencephalogram abnormalities, altered induced seizure thresholds, reduced anxiety-like behavior and alterations in olfactory-guided social behavior. Overall, Scn2aK1422E/+ mice present with phenotypes similar yet distinct from other Scn2a models, consistent with complex effects of K1422E on NaV1.2 channel function.

Rheumatoid arthritis (RA) and osteoporosis (OP) are two comorbid complex inflammatory conditions with evidence of shared genetic background and causal relationships. We aimed to clarify the genetic architecture underlying RA and various OP phenotypes while additionally considering an inflammatory component, C-reactive protein (CRP). Genome-wide association study summary statistics were acquired from the GEnetic Factors for OSteoporosis Consortium, Cohorts for Heart and Aging Research Consortium and UK Biobank. Mendelian randomization (MR) was used to detect the presence of causal relationships. Colocalization analysis was performed to determine shared genetic variants between CRP and OP phenotypes. Analysis of pleiotropy between traits owing to shared causal single nucleotide polymorphisms (SNPs) was performed using PL eiotropic A nalysis under CO mposite null hypothesis (PLACO). MR analysis was suggestive of horizontal pleiotropy between RA and OP traits. RA was a significant causal risk factor for CRP (β = 0.027, 95% confidence interval = 0.016-0.038). There was no evidence of CRP→OP causal relationship, but horizontal pleiotropy was apparent. Colocalization established shared genomic regions between CRP and OP, including GCKR and SERPINA1 genes. Pleiotropy arising from shared causal SNPs revealed through the colocalization analysis was all confirmed by PLACO. These genes were found to be involved in the same molecular function 'protein binding' (GO:0005515) associated with RA, OP and CRP. We identified three major components explaining the epidemiological relationship among RA, OP and inflammation: (1) Pleiotropy explains a portion of the shared genetic relationship between RA and OP, albeit polygenically; (2) RA contributes to CRP elevation and (3) CRP, which is influenced by RA, demonstrated pleiotropy with OP.

Mutations in the mitochondrial protein CHCHD2 cause autosomal dominant Parkinson's disease characterized by the preferential loss of substantia nigra dopamine (DA) neurons. Therefore, understanding the function of CHCHD2 in neurons may provide vital insights into how mitochondrial dysfunction contributes to neurodegeneration in PD. To investigate the normal requirement and function of CHCHD2 in neurons, we first examined CHCHD2 levels and showed that DA neurons have higher CHCHD2 levels than other neuron types, both in vivo and in co-culture. We then generated mice with either a targeted deletion of CHCHD2 in DA neurons or a deletion in the brain or total body. All three models were viable, and loss of CHCHD2 in the brain did not cause degeneration of DA neurons. Mice lacking CHCHD2 in DA neurons did display sex-specific changes to locomotor activity, but we did not observe differences in assays of muscle strength, exercise endurance or motor coordination. Furthermore, mitochondria derived from mice lacking CHCHD2 did not display abnormalities in OXPHOS function. Lastly, resilience to CHCHD2 deletion could not be explained by functional complementation by its paralog CHCHD10, as deletion of both CHCHD10 and CHCHD2 did not cause degeneration of DA neurons in the midbrain. These findings support the hypothesis that pathogenic CHCHD2 mutations cause PD through a toxic gain-of-function, rather than loss-of-function mechanism.

The molecular mechanisms leading to high-altitude pulmonary hypertension (HAPH) remains poorly understood. We previously analyzed the whole genome sequence of Kyrgyz highland population and identified eight genomic intervals having a potential role in HAPH. Tropomodulin 3 gene (TMOD3), which encodes a protein that binds and caps the pointed ends of actin filaments and inhibits cell migration, was one of the top candidates. Here we systematically sought additional evidence to validate the functional role of TMOD3. In-silico analysis reveals that some of the SNPs in HAPH associated genomic intervals were positioned in a regulatory region that could result in alternative splicing of TMOD3. In order to functionally validate the role of TMOD3 in HAPH, we exposed Tmod3-/+ mice to 4 weeks of constant hypoxia, i.e. 10% O2 and analyzed both functional (hemodynamic measurements) and structural (angiography) parameters related to HAPH. The hemodynamic measurements, such as right ventricular systolic pressure, a surrogate measure for pulmonary arterial systolic pressure, and right ventricular contractility (RV- ± dP/dt), increases with hypoxia did not separate between Tmod3-/+ and control mice. Remarkably, there was a significant increase in the number of lung vascular branches and total length of pulmonary vascular branches (P

BackgroundProprotein convertase subtilisin/kexin type 9 (PCSK9) is a key player in lipid metabolism, as it degrades low-density lipoprotein (LDL) receptors from hepatic cell membranes. So far, only variants of the PCSK9 gene locus were found to be associated with PCSK9 levels. Here we aimed to identify novel genetic loci that regulate PCSK9 levels and how they relate to other lipid traits. Additionally, we investigated to what extend the causal effect of PCSK9 on coronary artery disease (CAD) is mediated by low-density lipoprotein-cholesterol (LDL-C).Methods and resultsWe performed a genome-wide association study meta-analysis of PCSK9 levels in up to 12 721 samples of European ancestry. The estimated heritability was 10.3%, which increased to 12.6% using only samples from patients without statin treatment. We successfully replicated the known PCSK9 hit consisting of three independent signals. Interestingly, in a study of 300 African Americans, we confirmed the locus with a different PCSK9 variant. Beyond PCSK9, our meta-analysis detected three novel loci with genome-wide significance. Co-localization analysis with cis-eQTLs and lipid traits revealed biologically plausible candidate genes at two of them: APOB and TM6SF2. In a bivariate Mendelian Randomization analysis, we detected a strong effect of PCSK9 on LDL-C, but not vice versa. LDL-C mediated 63% of the total causal effect of PCSK9 on CAD.ConclusionOur study identified novel genetic loci with plausible candidate genes affecting PCSK9 levels. Ethnic heterogeneity was observed at the PCSK9 locus itself. Although the causal effect of PCSK9 on CAD is mainly mediated by LDL-C, an independent direct effect also occurs.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

We conducted cohort- and race-specific epigenome-wide association analyses of mitochondrial deoxyribonucleic acid (mtDNA) copy number (mtDNA CN) measured in whole blood from participants of African and European origins in five cohorts (n = 6182, mean age = 57-67 years, 65% women). In the meta-analysis of all the participants, we discovered 21 mtDNA CN-associated DNA methylation sites (CpG) (P

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P

Nonalcoholic fatty liver disease (NAFLD) is a leading cause of chronic liver disease and is highly correlated with metabolic disease. NAFLD results from environmental exposures acting on a susceptible polygenic background. This study performed the largest multiethnic investigation of exonic variation associated with NAFLD and correlated metabolic traits and diseases. An exome array meta-analysis was carried out among eight multiethnic population-based cohorts (n = 16 492) with computed tomography (CT) measured hepatic steatosis. A fixed effects meta-analysis identified five exome-wide significant loci (P

Angelman syndrome (AS) is a rare neurodevelopmental disorder characterized by impaired communication skills, ataxia, motor and balance deficits, intellectual disabilities, and seizures. The genetic cause of AS is the neuronal loss of UBE3A expression in the brain. A novel approach, described here, is a stem cell gene therapy which uses lentivector-transduced hematopoietic stem and progenitor cells to deliver functional UBE3A to affected cells. We have demonstrated both the prevention and reversal of AS phenotypes upon transplantation and engraftment of human CD34+ cells transduced with a Ube3a lentivector in a novel immunodeficient Ube3amat-/pat+ IL2rg-/y mouse model of AS. A significant improvement in motor and cognitive behavioral assays as well as normalized delta power measured by electroencephalogram was observed in neonates and adults transplanted with the gene modified cells. Human hematopoietic profiles observed in the lymphoid organs by detection of human immune cells were normal. Expression of UBE3A was detected in the brains of the adult treatment group following immunohistochemical staining illustrating engraftment of the gene-modified cells expressing UBE3A in the brain. As demonstrated with our data, this stem cell gene therapy approach offers a promising treatment strategy for AS, not requiring a critical treatment window.

Although several high-resolution recombination maps exist for European-descent populations, the recombination landscape of African populations remains relatively understudied. Given that there is high genetic divergence among groups in Africa, it is possible that recombination hotspots also diverge significantly. Both limitations and opportunities exist for developing recombination maps for these populations. In this review, we discuss various recombination inference methods, and the strengths and weaknesses of these methods in analyzing recombination in African-descent populations. Furthermore, we provide a decision tree and recommendations for which inference method to use in various research contexts. Establishing an appropriate methodology for recombination rate inference in a particular study will improve the accuracy of various downstream analyses including but not limited to local ancestry inference, haplotype phasing, fine-mapping of GWAS loci and genome assemblies.

Oligodendrocytes exist in a heterogenous state and are implicated in multiple neuropsychiatric diseases including dementia. Cortical oligodendrocytes are a glial population uniquely positioned to play a key role in neurodegeneration by synchronizing circuit connectivity but molecular pathways specific to this role are lacking. We utilized oligodendrocyte-specific translating ribosome affinity purification and RNA-seq (TRAP-seq) to transcriptionally profile adult mature oligodendrocytes from different regions of the central nervous system. Weighted gene co-expression network analysis reveals distinct region-specific gene networks. Two of these mature myelinating oligodendrocyte gene networks uniquely define cortical oligodendrocytes and differentially regulate cortical myelination (M8) and synaptic signaling (M4). These two cortical oligodendrocyte gene networks are enriched for genes associated with dementia including MAPT and include multiple gene targets of the regulatory microRNA, miR-142-3p. Using a combination of TRAP-qPCR, miR-142-3p overexpression in vitro, and miR-142-null mice, we show that miR-142-3p negatively regulates cortical myelination. In rTg4510 tau-overexpressing mice, cortical myelination is compromised, and tau-mediated neurodegeneration is associated with gene co-expression networks that recapitulate both the M8 and M4 cortical oligodendrocyte gene networks identified from normal cortex. We further demonstrate overlapping gene networks in mature oligodendrocytes present in normal cortex, rTg4510 and miR-142-null mice, and existing datasets from human tauopathies to provide evidence for a critical role of miR-142-3p-regulated cortical myelination and oligodendrocyte-mediated synaptic signaling in neurodegeneration.

Knock-in homozygote VCPR155H/R155H mutant mice are a lethal model of valosin-containing protein (VCP)-associated inclusion body myopathy associated with Paget disease of bone, frontotemporal dementia and amyotrophic lateral sclerosis. Ceramide (d18:1/16:0) levels are elevated in skeletal muscle of the mutant mice, compared to wild-type controls. Moreover, exposure to a lipid-enriched diet reverses lethality, improves myopathy and normalizes ceramide levels in these mutant mice, suggesting that dysfunctions in lipid-derived signaling are critical to disease pathogenesis. Here, we investigated the potential role of ceramide in VCP disease using pharmacological agents that manipulate the ceramide levels in myoblast cultures from VCP mutant mice and VCP patients. Myoblasts from wild-type, VCPR155H/+ and VCPR155H/R155H mice, as well as patient-induced pluripotent stem cells (iPSCs), were treated with an inhibitor of ceramide degradation to increase ceramide via acid ceramidase (ARN082) for proof of principle. Three chemically distinct inhibitors of ceramide biosynthesis via serine palmitoyl-CoA transferase (L-cycloserine, myriocin or ARN14494) were used as a therapeutic strategy to reduce ceramide in myoblasts. Acid ceramidase inhibitor, ARN082, elevated cellular ceramide levels and concomitantly enhanced pathology. Conversely, inhibitors of ceramide biosynthesis L-cycloserine, myriocin and ARN14494 reduced ceramide production. The results point to ceramide-mediated signaling as a key contributor to pathogenesis in VCP disease and suggest that manipulating this pathway by blocking ceramide biosynthesis might exert beneficial effects in patients with this condition. The ceramide pathway appears to be critical in VCP pathogenesis, and small-molecule inhibitors of ceramide biosynthesis might provide therapeutic benefits in VCP and related neurodegenerative diseases.

Previously we showed that dimethyl fumarate (DMF) dose-dependently increased mitochondrial gene expression and function in cells and might be considered as a therapeutic for inherited mitochondrial disease, including Friedreich's ataxia (FA). Here we tested DMF's ability to dose-dependently increase mitochondrial function, mitochondrial gene expression (frataxin and cytochrome oxidase protein) and mitochondrial copy number in C57BL6 wild-type mice and the FXNKD mouse model of FA. We first dosed DMF at 0-320 mg/kg in C57BL6 mice and observed significant toxicity above 160 mg/kg orally, defining the maximum tolerated dose. Oral dosing of C57BL6 mice in the range 0-160 mg/kg identified a maximum increase in aconitase activity and mitochondrial gene expression in brain and quadriceps at 110 mg/kg DMF, thus defining the maximum effective dose (MED). The MED of DMF in mice overlaps the currently approved human-equivalent doses of DMF prescribed for multiple sclerosis (480 mg/day) and psoriasis (720 mg/day). In the FXNKD mouse model of FA, which has a doxycycline-induced deficit of frataxin protein, we observed significant decreases of multiple mitochondrial parameters, including deficits in brain mitochondrial Complex 2, Complex 4 and aconitase activity, supporting the idea that frataxin deficiency reduces mitochondrial gene expression, mitochondrial functions and biogenesis. About 110 mg/kg of oral DMF rescued these enzyme activities in brain and rescued frataxin and cytochrome oxidase expression in brain, cerebellum and quadriceps muscle of the FXNKD mouse model. Taken together, these results support the idea of using fumarate-based molecules to treat FA or other mitochondrial diseases.

Neonatal dried blood spots (NDBS) are a widely banked sample source that enables retrospective investigation into early life molecular events. Here, we performed low-pass whole genome bisulfite sequencing (WGBS) of 86 NDBS DNA to examine early life Down syndrome (DS) DNA methylation profiles. DS represents an example of genetics shaping epigenetics, as multiple array-based studies have demonstrated that trisomy 21 is characterized by genome-wide alterations to DNA methylation. By assaying over 24 million CpG sites, thousands of genome-wide significant (q

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

The RNA exosome is an essential ribonuclease complex required for processing and/or degradation of both coding and non-coding RNAs. We identified five patients with biallelic variants in EXOSC5, which encodes a structural subunit of the RNA exosome. The clinical features of these patients include failure to thrive, short stature, feeding difficulties, developmental delays that affect motor skills, hypotonia and esotropia. Brain MRI revealed cerebellar hypoplasia and ventriculomegaly. While we ascertained five patients, three patients with distinct variants of EXOSC5 were studied in detail. The first patient had a deletion involving exons 5-6 of EXOSC5 and a missense variant, p.Thr114Ile, that were inherited in trans, the second patient was homozygous for p.Leu206His and the third patient had paternal isodisomy for chromosome 19 and was homozygous for p.Met148Thr. The additional two patients ascertained are siblings who had an early frameshift mutation in EXOSC5 and the p.Thr114Ile missense variant that were inherited in trans. We employed three complementary approaches to explore the requirement for EXOSC5 in brain development and assess consequences of pathogenic EXOSC5 variants. Loss of function for exosc5 in zebrafish results in shortened and curved tails/bodies, reduced eye/head size and edema. We modeled pathogenic EXOSC5 variants in both budding yeast and mammalian cells. Some of these variants cause defects in RNA exosome function as well as altered interactions with other RNA exosome subunits. These findings expand the number of genes encoding RNA exosome subunits linked to human disease while also suggesting that disease mechanism varies depending on the specific pathogenic variant.

Hundreds of loci have been associated with blood pressure (BP) traits from many genome-wide association studies. We identified an enrichment of these loci in aorta and tibial artery expression quantitative trait loci in our previous work in ~100 000 Genetic Epidemiology Research on Aging study participants. In the present study, we sought to fine-map known loci and identify novel genes by determining putative regulatory regions for these and other tissues relevant to BP. We constructed maps of putative cis-regulatory elements (CREs) using publicly available open chromatin data for the heart, aorta and tibial arteries, and multiple kidney cell types. Variants within these regions may be evaluated quantitatively for their tissue- or cell-type-specific regulatory impact using deltaSVM functional scores, as described in our previous work. We aggregate variants within these putative CREs within 50 Kb of the start or end of 'expressed' genes in these tissues or cell types using public expression data and use deltaSVM scores as weights in the group-wise sequence kernel association test to identify candidates. We test for association with both BP traits and expression within these tissues or cell types of interest and identify the candidates MTHFR, C10orf32, CSK, NOV, ULK4, SDCCAG8, SCAMP5, RPP25, HDGFRP3, VPS37B and PPCDC. Additionally, we examined two known QT interval genes, SCN5A and NOS1AP, in the Atherosclerosis Risk in Communities Study, as a positive control, and observed the expected heart-specific effect. Thus, our method identifies variants and genes for further functional testing using tissue- or cell-type-specific putative regulatory information.

Altered cellular metabolism is believed to be an important contributor to pathogenesis of the neurodegenerative disorder Huntington's disease (HD). Research has primarily focused on mitochondrial toxicity, which can cause death of the vulnerable striatal neurons, but other aspects of metabolism have also been implicated. Most previous studies have been carried out using postmortem human brain or non-human cells. Here, we studied bioenergetics in an induced pluripotent stem cell-based model of the disease. We found decreased adenosine triphosphate (ATP) levels in HD cells compared to controls across differentiation stages and protocols. Proteomics data and multiomics network analysis revealed normal or increased levels of mitochondrial messages and proteins, but lowered expression of glycolytic enzymes. Metabolic experiments showed decreased spare glycolytic capacity in HD neurons, while maximal and spare respiratory capacities driven by oxidative phosphorylation were largely unchanged. ATP levels in HD neurons could be rescued with addition of pyruvate or late glycolytic metabolites, but not earlier glycolytic metabolites, suggesting a role for glycolytic deficits as part of the metabolic disturbance in HD neurons. Pyruvate or other related metabolic supplements could have therapeutic benefit in HD.

Altered cellular metabolism is believed to be an important contributor to pathogenesis of the neurodegenerative disorder Huntington's disease (HD). Research has primarily focused on mitochondrial toxicity, which can cause death of the vulnerable striatal neurons, but other aspects of metabolism have also been implicated. Most previous studies have been carried out using postmortem human brain or non-human cells. Here, we studied bioenergetics in an induced pluripotent stem cell-based model of the disease. We found decreased adenosine triphosphate (ATP) levels in HD cells compared to controls across differentiation stages and protocols. Proteomics data and multiomics network analysis revealed normal or increased levels of mitochondrial messages and proteins, but lowered expression of glycolytic enzymes. Metabolic experiments showed decreased spare glycolytic capacity in HD neurons, while maximal and spare respiratory capacities driven by oxidative phosphorylation were largely unchanged. ATP levels in HD neurons could be rescued with addition of pyruvate or late glycolytic metabolites, but not earlier glycolytic metabolites, suggesting a role for glycolytic deficits as part of the metabolic disturbance in HD neurons. Pyruvate or other related metabolic supplements could have therapeutic benefit in HD.

Pre-B cell leukemia factor 1 (PBX1) is an essential developmental transcription factor, mutations in which have recently been associated with CAKUTHED syndrome, characterized by multiple congenital defects including congenital heart disease (CHD). During analysis of a whole-exome-sequenced cohort of heterogeneous CHD patients, we identified a de novo missense variant, PBX1:c.551G>C p.R184P, in a patient with tetralogy of Fallot with absent pulmonary valve and extra-cardiac phenotypes. Functional analysis of this variant by creating a CRISPR-Cas9 gene-edited mouse model revealed multiple congenital anomalies. Congenital heart defects (persistent truncus arteriosus and ventricular septal defect), hypoplastic lungs, hypoplastic/ectopic kidneys, aplastic adrenal glands and spleen, as well as atretic trachea and palate defects were observed in the homozygous mutant embryos at multiple stages of development. We also observed developmental anomalies in a proportion of heterozygous embryos, suggestive of a dominant mode of inheritance. Analysis of gene expression and protein levels revealed that although Pbx1 transcripts are higher in homozygotes, amounts of PBX1 protein are significantly decreased. Here, we have presented the first functional model of a missense PBX1 variant and provided strong evidence that p.R184P is disease-causal. Our findings also expand the phenotypic spectrum associated with pathogenic PBX1 variants in both humans and mice.

Huntington's disease (HD) is caused by an expansion of a poly glutamine (polyQ) stretch in the huntingtin protein (HTT) that is necessary to cause pathology and formation of HTT aggregates. Here we ask whether expanded polyQ is sufficient to cause pathology and aggregate formation. By addressing the sufficiency question, one can identify cellular processes and structural parameters that influence HD pathology and HTT subcellular behavior (i.e. aggregation state and subcellular location). Using Drosophila, we compare the effects of expressing mutant full-length human HTT (fl-mHTT) to the effects of mutant human HTTexon1 and to two commonly used synthetic fragments, HTT171 and shortstop (HTT118). Expanded polyQ alone is not sufficient to cause inclusion formation since full-length HTT and HTTex1 with expanded polyQ are both toxic although full-length HTT remains diffuse while HTTex1 forms inclusions. Further, inclusions are not sufficient to cause pathology since HTT171-120Q forms inclusions but is benign and co-expression of HTT171-120Q with non-aggregating pathogenic fl-mHTT recruits fl-mHTT to aggregates and rescues its pathogenicity. Additionally, the influence of sequences outside the expanded polyQ domain is revealed by finding that small modifications to the HTT118 or HTT171 fragments can dramatically alter their subcellular behavior and pathogenicity. Finally, mutant HTT subcellular behavior is strongly modified by different cell and tissue environments (e.g. fl-mHTT appears as diffuse nuclear in one tissue and diffuse cytoplasmic in another but toxic in both). These observations underscore the importance of cellular and structural context for the interpretation and comparison of experiments using different fragments and tissues to report the effects of expanded polyQ.

Frontotemporal dementia (FTD) is an early onset dementia characterized by progressive atrophy of the frontal and/or temporal lobes. FTD is highly heritable with mutations in progranulin accounting for 5-26% of cases in different populations. Progranulin is involved in endocytosis, secretion and lysosomal processes, but its functions under physiological and pathological conditions remains to be defined. Many FTD-causing non-sense progranulin mutations contain a premature termination codon (PTC), thus progranulin haploinsufficiency has been proposed as a major disease mechanism. Currently, there is no effective FTD treatment or therapy. Aminoglycosides are a class of antibiotics that possess a less-known function to induce eukaryotic ribosomal readthrough of PTCs to produce a full-length protein. The aminoglycoside-induced readthrough strategy has been utilized to treat multiple human diseases caused by PTCs. In this study, we tested the only clinically approved readthrough small molecule PTC124 and 11 aminoglycosides in a cell culture system on four PTCs responsible for FTD or a related neurodegenerative disease amyotrophic lateral sclerosis. We found that the aminoglycosides G418 and gentamicin rescued the expression of the progranulin R493X mutation. G418 was more effective than gentamicin (~50% rescue versus

Mitochondrial disorders are the result of nuclear and mitochondrial DNA mutations that affect multiple organs, with the central and peripheral nervous system often affected. Currently, there is no cure for mitochondrial disorders. Currently, gene therapy offers a novel approach for treating monogenetic disorders, including nuclear genes associated with mitochondrial disorders. We utilized a mouse model carrying a knockout of the mitochondrial fusion-fission-related gene solute carrier family 25 member 46 (Slc25a46) and treated them with neurotrophic AAV-PHP.B vector carrying the mouse Slc25a46 coding sequence. Thereafter, we used immunofluorescence staining and western blot to test the transduction efficiency of this vector. Toluidine blue staining and electronic microscopy were utilized to assess the morphology of optic and sciatic nerves following treatment, and the morphology and respiratory chain activity of mitochondria within these tissues were determined as well. The adeno-associated virus (AAV) vector effectively transduced in the cerebrum, cerebellum, heart, liver and sciatic nerves. AAV-Slc25a46 treatment was able to rescue the premature death in the mutant mice (Slc25a46-/-). The treatment-improved electronic conductivity of the peripheral nerves increased mobility and restored mitochondrial complex activities. Most notably, mitochondrial morphology inside the tissues of both the central and peripheral nervous systems was normalized, and the neurodegeneration, chronic neuroinflammation and loss of Purkinje cell dendrites observed within the mutant mice were alleviated. Overall, our study shows that AAV-PHP.B's neurotrophic properties are plausible for treating conditions where the central nervous system is affected, such as many mitochondrial diseases, and that AAV-Slc25a46 could be a novel approach for treating SLC25A46-related mitochondrial disorders.

Transcriptional and epigenetic alterations occur early in Huntington's disease (HD), and treatment with epigenetic modulators is beneficial in several HD animal models. The drug JQ1, which inhibits histone acetyl-lysine reader bromodomains, has shown promise for multiple cancers and neurodegenerative disease. We tested whether JQ1 could improve behavioral phenotypes in the R6/2 mouse model of HD and modulate HD-associated changes in transcription and epigenomics. R6/2 and non-transgenic (NT) mice were treated with JQ1 daily from 5 to 11 weeks of age and behavioral phenotypes evaluated over this period. Following the trial, cortex and striatum were isolated and subjected to mRNA-seq and ChIP-seq for the histone marks H3K4me3 and H3K27ac. Initially, JQ1 enhanced motor performance in NT mice. In R6/2 mice, however, JQ1 had no effect on rotarod or grip strength but exacerbated weight loss and worsened performance on the pole test. JQ1-induced gene expression changes in NT mice were distinct from those in R6/2 and primarily involved protein translation and bioenergetics pathways. Dysregulation of HD-related pathways in striatum was exacerbated by JQ1 in R6/2 mice, but not in NTs, and JQ1 caused a corresponding increase in the formation of a mutant huntingtin protein-dependent high molecular weight species associated with pathogenesis. This study suggests that drugs predicted to be beneficial based on their mode of action and effects in wild-type or in other neurodegenerative disease models may have an altered impact in the HD context. These observations have important implications in the development of epigenetic modulators as therapies for HD.

Genetic diversity provides a rich repository for understanding the role of proteostasis in the management of the protein fold in human biology. Failure in proteostasis can trigger multiple disease states, affecting both human health and lifespan. Niemann-Pick C1 (NPC1) disease is a rare genetic disorder triggered by mutations in NPC1, a multi-spanning transmembrane protein that is trafficked through the exocytic pathway to late endosomes (LE) and lysosomes (Ly) (LE/Ly) to globally manage cholesterol homeostasis. Defects triggered by >300 NPC1 variants found in the human population inhibit export of NPC1 protein from the endoplasmic reticulum (ER) and/or function in downstream LE/Ly, leading to cholesterol accumulation and onset of neurodegeneration in childhood. We now show that the allosteric inhibitor JG98, that targets the cytosolic Hsp70 chaperone/co-chaperone complex, can significantly improve the trafficking and post-ER protein level of diverse NPC1 variants. Using a new approach to model genetic diversity in human disease, referred to as variation spatial profiling, we show quantitatively how JG98 alters the Hsp70 chaperone/co-chaperone system to adjust the spatial covariance (SCV) tolerance and set-points on an amino acid residue-by-residue basis in NPC1 to differentially regulate variant trafficking, stability, and cholesterol homeostasis, results consistent with the role of BCL2-associated athanogene family co-chaperones in managing the folding status of NPC1 variants. We propose that targeting the cytosolic Hsp70 system by allosteric regulation of its chaperone/co-chaperone based client relationships can be used to adjust the SCV tolerance of proteostasis buffering capacity to provide an approach to mitigate systemic and neurological disease in the NPC1 population.

Ubiquitin E3 ligase 3A (UBE3A) encodes an E3 ubiquitin ligase whose loss from the maternal allele causes the neurodevelopmental disorder Angelman syndrome (AS). Previous studies of UBE3A function have not examined full Ube3a deletion in mouse, the complexity of imprinted gene networks in brain nor the molecular basis of systems-level cognitive dysfunctions in AS. We therefore utilized a systems biology approach to elucidate how UBE3A loss impacts the early postnatal brain in a novel CRISPR/Cas9-engineered rat Angelman model of a complete Ube3a deletion. Strand-specific transcriptome analysis of offspring from maternally or paternally inherited Ube3a deletions revealed the expected parental expression patterns of Ube3a sense and antisense transcripts by postnatal day 2 (P2) in hypothalamus and day 9 (P9) in cortex, compared to wild-type littermates. The dependency of genome-wide effects on parent-of-origin, Ube3a genotype and time (P2 and P9) was investigated through transcriptome (RNA sequencing of cortex and hypothalamus) and methylome (whole-genome bisulfite sequencing of hypothalamus). Weighted gene co-expression and co-methylation network analyses identified co-regulated networks in maternally inherited Ube3a deletion offspring enriched in postnatal developmental processes including Wnt signaling, synaptic regulation, neuronal and glial functions, epigenetic regulation, ubiquitin, circadian entrainment and splicing. Furthermore, we showed that loss of the paternal Ube3a antisense transcript resulted in both unique and overlapping dysregulated gene pathways with maternal loss, predominantly at the level of differential methylation. Together, these results provide a holistic examination of the molecular impacts of UBE3A loss in brain, supporting the existence of interactive epigenetic networks between maternal and paternal transcripts at the Ube3a locus.

Humans are social animals that experience intense suffering when they perceive a lack of social connection. Modern societies are experiencing an epidemic of loneliness. Although the experience of loneliness is universally human, some people report experiencing greater loneliness than others. Loneliness is more strongly associated with mortality than obesity, emphasizing the need to understand the nature of the relationship between loneliness and health. Although it is intuitive that circumstantial factors such as marital status and age influence loneliness, there is also compelling evidence of a genetic predisposition toward loneliness. To better understand the genetic architecture of loneliness and its relationship with associated outcomes, we extended the genome-wide association study meta-analysis of loneliness to 511 280 subjects, and detect 19 significant genetic variants from 16 loci, including four novel loci, as well as 58 significantly associated genes. We investigated the genetic overlap with a wide range of physical and mental health traits by computing genetic correlations and by building loneliness polygenic scores in an independent sample of 18 498 individuals with EHR data to conduct a PheWAS with. A genetic predisposition toward loneliness was associated with cardiovascular, psychiatric, and metabolic disorders and triglycerides and high-density lipoproteins. Mendelian randomization analyses showed evidence of a causal, increasing, the effect of both BMI and body fat on loneliness. Our results provide a framework for future studies of the genetic basis of loneliness and its relationship to mental and physical health.

Single germline or somatic activating mutations of mammalian target of rapamycin (mTOR) pathway genes are emerging as a major cause of type II focal cortical dysplasia (FCD), hemimegalencephaly (HME) and tuberous sclerosis complex (TSC). A double-hit mechanism, based on a primary germline mutation in one allele and a secondary somatic hit affecting the other allele of the same gene in a small number of cells, has been documented in some patients with TSC or FCD. In a patient with HME, severe intellectual disability, intractable seizures and hypochromic skin patches, we identified the ribosomal protein S6 (RPS6) p.R232H variant, present as somatic mosaicism at ~15.1% in dysplastic brain tissue and ~11% in blood, and the MTOR p.S2215F variant, detected as ~8.8% mosaicism in brain tissue, but not in blood. Overexpressing the two variants independently in animal models, we demonstrated that MTOR p.S2215F caused neuronal migration delay and cytomegaly, while RPS6 p.R232H prompted increased cell proliferation. Double mutants exhibited a more severe phenotype, with increased proliferation and migration defects at embryonic stage and, at postnatal stage, cytomegalic cells exhibiting eccentric nuclei and binucleation, which are typical features of balloon cells. These findings suggest a synergistic effect of the two variants. This study indicates that, in addition to single activating mutations and double-hit inactivating mutations in mTOR pathway genes, severe forms of cortical dysplasia can also result from activating mutations affecting different genes in this pathway. RPS6 is a potential novel disease-related gene.

Duchenne muscular dystrophy (DMD) is caused by mutations in the gene encoding dystrophin. Prior work has shown that DMD progression can vary, depending on the genetic makeup of the patient. Several modifier alleles have been identified including LTBP4 and SPP1. We previously showed that Spp1 exacerbates the DMD phenotype in the mdx mouse model by promoting fibrosis and by skewing macrophage polarization. Here, we studied the mechanisms involved in Spp1's promotion of fibrosis by using both isolated fibroblasts and genetically modified mice. We found that Spp1 upregulates collagen expression in mdx fibroblasts by enhancing TGFβ signaling. Spp1's effects on TGFβ signaling are through induction of MMP9 expression. MMP9 is a protease that can release active TGFβ ligand from its latent complex. In support for activation of this pathway in our model, we showed that treatment of mdx fibroblasts with MMP9 inhibitor led to accumulation of the TGFβ latent complex, decreased levels of active TGFβ and reduced collagen expression. Correspondingly, we found reduced active TGFβ in Spp1-/-mdxB10 and Mmp9-/-mdxB10 muscles in vivo. Taken together with previous observations of reduced fibrosis in both models, these data suggest that Spp1 acts upstream of TGFβ to promote fibrosis in mdx muscles. We found that in the context of constitutively upregulated TGFβ signaling (such as in the mdxD2 model), ablation of Spp1 has very little effect on fibrosis. Finally, we performed proof-of-concept studies showing that postnatal pharmacological inhibition of Spp1 reduces fibrosis and improves muscle function in mdx mice.

Agrin is a large extracellular matrix protein whose isoforms differ in their tissue distribution and function. Motoneuron-derived y+z+ agrin regulates the formation of the neuromuscular junction (NMJ), while y-z- agrin is widely expressed and has diverse functions. Previously we identified a missense mutation (V1727F) in the second laminin globular (LG2) domain of agrin that causes severe congenital myasthenic syndrome. Here, we define pathogenic effects of the agrin V1727F mutation that account for the profound dysfunction of the NMJ. First, by expressing agrin variants in heterologous cells, we show that the V1727F mutation reduces the secretion of y+z+ agrin compared to wild type, whereas it has no effect on the secretion of y-z- agrin. Second, we find that the V1727F mutation significantly impairs binding of y+z+ agrin to both heparin and the low-density lipoprotein receptor-related protein 4 (LRP4) coreceptor. Third, molecular modeling of the LG2 domain suggests that the V1727F mutation primarily disrupts the y splice insert, and consistent with this we find that it partially occludes the contribution of the y splice insert to agrin binding to heparin and LRP4. Together, these findings identify several pathogenic effects of the V1727F mutation that reduce its expression and ability to bind heparan sulfate proteoglycan and LRP4 coreceptors involved in the muscle-specific kinase signaling pathway. These defects primarily impair the function of neural y+z+ agrin and combine to cause a severe CMS phenotype, whereas y-z- agrin function in other tissues appears preserved.

DNA methylation acts at the interface of genetic and environmental factors relevant for autism spectrum disorder (ASD). Placenta, normally discarded at birth, is a potentially rich source of DNA methylation patterns predictive of ASD in the child. Here, we performed whole methylome analyses of placentas from a prospective study MARBLES (Markers of Autism Risk in Babies-Learning Early Signs) of high-risk pregnancies. A total of 400 differentially methylated regions (DMRs) discriminated placentas stored from children later diagnosed with ASD compared to typically developing controls. These ASD DMRs were significantly enriched at promoters, mapped to 596 genes functionally enriched in neuronal development, and overlapped genetic ASD risk. ASD DMRs at CYP2E1 and IRS2 reached genome-wide significance, replicated by pyrosequencing and correlated with expression differences in brain. Methylation at CYP2E1 associated with both ASD diagnosis and genotype within the DMR. In contrast, methylation at IRS2 was unaffected by within DMR genotype but modified by preconceptional maternal prenatal vitamin use. This study therefore identified two potentially useful early epigenetic markers for ASD in placenta.

Global gene-expression analysis has shown remarkable difference between males and females in response to exposure to many diseases. Nevertheless, gene expression studies in asthmatics have so far focused on sex-combined analysis, ignoring inherent variabilities between the sexes, which potentially drive disparities in asthma prevalence. The objectives of this study were to identify (1) sex-specific differentially expressed genes (DEGs), (2) genes that show sex-interaction effects and (3) sex-specific pathways and networks enriched in asthma risk. We analyzed 711 males and 689 females and more than 2.8 million transcripts covering 20 000 genes leveraged from five different tissues and cell types (i.e. epithelial, blood, induced sputum, T cells and lymphoblastoids). Using tissue-specific meta-analysis, we identified 439 male- and 297 female-specific DEGs in all cell types, with 32 genes in common. By linking DEGs to the genome-wide association study (GWAS) catalog and the lung and blood eQTL annotation data from GTEx, we identified four male-specific genes (FBXL7, ITPR3 and RAD51B from epithelial tissue and ALOX15 from blood) and one female-specific gene (HLA-DQA1 from epithelial tissue) that are disregulated during asthma. The hypoxia-inducible factor 1 signaling pathway was enriched only in males, and IL-17 and chemokine signaling pathways were enriched in females. The cytokine-cytokine signaling pathway was enriched in both sexes. The presence of sex-specific genes and pathways demonstrates that sex-combined analysis does not identify genes preferentially expressed in each sex in response to diseases. Linking DEG and molecular eQTLs to GWAS catalog represents an important avenue for identifying biologically and clinically relevant genes.

Autism spectrum disorder (ASD) encompasses a collection of complex neuropsychiatric disorders characterized by deficits in social functioning, communication and repetitive behaviour. Building on recent studies supporting a role for developmentally moderated regulatory genomic variation in the molecular aetiology of ASD, we quantified genome-wide patterns of DNA methylation in 223 post-mortem tissues samples isolated from three brain regions [prefrontal cortex, temporal cortex and cerebellum (CB)] dissected from 43 ASD patients and 38 non-psychiatric control donors. We identified widespread differences in DNA methylation associated with idiopathic ASD (iASD), with consistent signals in both cortical regions that were distinct to those observed in the CB. Individuals carrying a duplication on chromosome 15q (dup15q), representing a genetically defined subtype of ASD, were characterized by striking differences in DNA methylationacross a discrete domain spanning an imprinted gene cluster within the duplicated region. In addition to the dramatic cis-effects on DNA methylation observed in dup15q carriers, we identified convergent methylomic signatures associated with both iASD and dup15q, reflecting the findings from previous studies of gene expression and H3K27ac. Cortical co-methylation network analysis identified a number of co-methylated modules significantly associated with ASD that are enriched for genomic regions annotated to genes involved in the immune system, synaptic signalling and neuronal regulation. Our study represents the first systematic analysis of DNA methylation associated with ASD across multiple brain regions, providing novel evidence for convergent molecular signatures associated with both idiopathic and syndromic autism.

Mutations in the microtubule-associated protein tau (MAPT) underlie multiple neurodegenerative disorders, yet the pathophysiological mechanisms are unclear. A novel variant in MAPT resulting in an alanine to threonine substitution at position 152 (A152T tau) has recently been described as a significant risk factor for both frontotemporal lobar degeneration and Alzheimer's disease. Here we use complementary computational, biochemical, molecular, genetic and imaging approaches in Caenorhabditis elegans and mouse models to interrogate the effects of the A152T variant on tau function. In silico analysis suggests that a threonine at position 152 of tau confers a new phosphorylation site. This finding is borne out by mass spectrometric survey of A152T tau phosphorylation in C. elegans and mouse. Optical pulse-chase experiments of Dendra2-tau demonstrate that A152T tau and phosphomimetic A152E tau exhibit increased diffusion kinetics and the ability to traverse across the axon initial segment more efficiently than wild-type (WT) tau. A C. elegans model of tauopathy reveals that A152T and A152E tau confer patterns of developmental toxicity distinct from WT tau, likely due to differential effects on retrograde axonal transport. These data support a role for phosphorylation of the variant threonine in A152T tau toxicity and suggest a mechanism involving impaired retrograde axonal transport contributing to human neurodegenerative disease.

Obstructive sleep apnea (OSA) is a common disorder associated with increased risk of cardiovascular disease and mortality. Its prevalence and severity vary across ancestral background. Although OSA traits are heritable, few genetic associations have been identified. To identify genetic regions associated with OSA and improve statistical power, we applied admixture mapping on three primary OSA traits [the apnea hypopnea index (AHI), overnight average oxyhemoglobin saturation (SaO2) and percentage time SaO2

Birdshot Uveitis (Birdshot) is a rare eye condition that affects HLA-A29-positive individuals and could be considered a prototypic member of the recently proposed 'MHC-I (major histocompatibility complex class I)-opathy' family. Genetic studies have pinpointed the endoplasmic reticulum aminopeptidase (ERAP1) and (ERAP2) genes as shared associations across MHC-I-opathies, which suggests ERAP dysfunction may be a root cause for MHC-I-opathies. We mapped the ERAP1 and ERAP2 haplotypes in 84 Dutch cases and 890 controls. We identified association at variant rs10044354, which mediated a marked increase in ERAP2 expression. We also identified and cloned an independently associated ERAP1 haplotype (tagged by rs2287987) present in more than half of the cases; this ERAP1 haplotype is also the primary risk and protective haplotype for other MHC-I-opathies. We show that the risk ERAP1 haplotype conferred significantly altered expression of ERAP1 isoforms in transcriptomic data (n = 360), resulting in lowered protein expression and distinct enzymatic activity. Both the association for rs10044354 (meta-analysis: odds ratio (OR) [95% CI]=2.07[1.58-2.71], P = 1.24 × 10(-7)) and rs2287987 (OR[95% CI]: =2.01[1.51-2.67], P = 1.41 × 10(-6)) replicated and showed consistent direction of effect in an independent Spanish cohort of 46 cases and 2103 controls. In both cohorts, the combined rs2287987-rs10044354 haplotype associated with Birdshot more strongly than either variant alone [meta-analysis: P=3.9 × 10(-9)]. Finally, we observed that ERAP2 protein expression is dependent on the ERAP1 background across three European populations (n = 3353). In conclusion, a functionally distinct combination of ERAP1 and ERAP2 are a hallmark of Birdshot and provide rationale for strategies designed to correct ERAP function for treatment of Birdshot and MHC-I-opathies more broadly.

The genetic basis of human neural tube defects (NTDs), such as anencephaly and spina bifida (SB), is complex and heterogeneous. Grainyhead-like genes represent candidates for involvement in NTDs based on the presence of SB and exencephaly in mice carrying loss-of-function alleles of Grhl2 or Grhl3. We found that reinstatement of Grhl3 expression, by bacterial artificial chromosome (BAC)-mediated transgenesis, prevents SB in Grhl3-null embryos, as in the Grhl3 hypomorphic curly tail strain. Notably, however, further increase in expression of Grhl3 causes highly penetrant SB. Grhl3 overexpression recapitulates the spinal NTD phenotype of loss-of-function embryos, although the underlying mechanism differs. However, it does not phenocopy other defects of Grhl3-null embryos such as abnormal axial curvature, cranial NTDs (exencephaly) or skin barrier defects, the latter being rescued by the Grhl3-transgene. Grhl2 and Grhl3 can form homodimers and heterodimers, suggesting a possible model in which defects arising from overexpression of Grhl3 result from sequestration of Grhl2 in heterodimers, mimicking Grhl2 loss of function. This hypothesis predicts that increased abundance of Grhl2 would have an ameliorating effect in Grhl3 overexpressing embryo. Instead, we observed a striking additive genetic interaction between Grhl2 and Grhl3 gain-of-function alleles. Severe SB arose in embryos in which both genes were expressed at moderately elevated levels that individually do not cause NTDs. Furthermore, moderate Grhl3 overexpression also interacted with the Vangl2Lp allele to cause SB, demonstrating genetic interaction with the planar cell polarity signalling pathway that is implicated in mouse and human NTDs.