Cardiometabolic diseases, such as type 2 diabetes and cardiovascular disease, have a high public health burden. Understanding the genetically-determined regulation of proteins that are dysregulated in disease can help to dissect the complex biology underpinning them. Here, we perform a protein quantitative trait locus (pQTL) analysis of 255 serum proteins relevant to cardiometabolic processes in 2893 individuals. Meta-analysing whole-genome sequencing (WGS) data from two Greek cohorts, MANOLIS (n = 1356; 22.5x WGS) and Pomak (n = 1537; 18.4x WGS), we detect 302 independently-associated pQTL variants for 171 proteins, including 12 rare variants (minor allele frequency [MAF] < 1%). We additionally find 15 pQTL variants that are rare in non-Finnish European populations, but have drifted up in frequency in the discovery cohorts here. We identify proteins causally associated with cardiometabolic traits, including MEP1B for high-density lipoprotein levels; and describe a knock-out Mep1b mouse model. Our findings furnish insights into the genetic architecture of the serum proteome, identify new protein-disease relationships, and demonstrate the importance of isolated populations in pQTL analysis.

Loss of heterozygosity (LOH) is a genetic alteration that results from the loss of one allele at a heterozygous locus. In particular, copy neutral LOH (CN-LOH) events are generated, for example, by mitotic homologous recombination after monoallelic defection or gene conversion, resulting in novel homozygous locus having two copies of the normal counterpart allele. This phenomenon can serve as a source of genome diversity and is associated with various diseases. To clarify the nature of the CN-LOH such as the frequency, genomic distribution, and inheritance pattern, we made use of whole-genome sequencing data of the three-generation CEPH/Utah family cohort, with the pedigree consisting of grandparents, parents, and offspring. We identified an average of 40.7 CN-LOH events per individual taking advantage of 285 healthy individuals from 33 families in the cohort. On average 65% of them were classified as gonosomal-mosaicism-associated CN-LOH, which exists in both germline and somatic cells. We also confirmed that the incidence of the CN-LOH has little to do with the parents' age and sex. Furthermore, through the analysis of the genomic region including the CN-LOH, we found that the chance of the occurrence of the CN-LOH tends to increase at the GC-rich locus and/or on the chromosome having a relatively close inter-homolog distance. We expect that these results provide significant insights into the association between genetic alteration and spatial position of chromosomes as well as the intrinsic genetic property of the CN-LOH.

ADP-ribosylation factor 1 (ARF1) is a small GTPase that regulates membrane traffic at the Golgi apparatus and endosomes through recruitment of several coat proteins and lipid-modifying enzymes. Here, we report a pediatric patient with an ARF1-related disorder due to a monoallelic de novo missense variant (c.296 G > A; p.R99H) in the ARF1 gene, associated with developmental delay, hypotonia, intellectual disability, and motor stereotypies. Neuroimaging revealed a hypoplastic corpus callosum and subcortical white-matter abnormalities. Notably, this patient did not exhibit periventricular heterotopias previously observed in other patients with ARF1 variants (including p.R99H). Functional analysis of the R99H-ARF1 variant protein revealed that it was expressed at normal levels and properly localized to the Golgi apparatus; however, expression of this variant caused swelling of the Golgi apparatus, increased the recruitment of coat proteins such as COPI, AP-1, and GGA3, and altered the morphology of recycling endosomes. In addition, we observed that expression of R99H-ARF1 prevented dispersal of the Golgi apparatus by the ARF1-inhibitor brefeldin A. Finally, protein-interaction analyses showed that R99H-ARF1 bound more tightly to the ARF1-effector GGA3 relative to wild-type ARF1. These properties were similar to those of the well-characterized constitutively-active Q71L-ARF1 mutant, indicating that the pathogenetic mechanism of the R99H-ARF1 variant involves constitutive activation with resultant Golgi and endosomal alterations. The absence of periventricular nodular heterotopias in this R99H-ARF1 subject also indicates this finding may not be a consistent phenotypic expression of all ARF1-related disorders.

The principal component of the protein homeostasis network is the ubiquitin-proteasome system. Ubiquitination is mediated by an enzymatic cascade involving, i.e. E3 ubiquitin ligases, many of which belong to the cullin-RING ligases family. Genetic defects in the ubiquitin-proteasome system components, including cullin-RING ligases, are known causes of neurodevelopmental disorders. Using exome sequencing to diagnose a pediatric patient with developmental delay, pyramidal signs, and limb ataxia, we identified a de novo missense variant c.376G > C; p.(Asp126His) in the FEM1C gene encoding a cullin-RING ligase substrate receptor. This variant alters a conserved amino acid located within a highly constrained coding region and is predicted as pathogenic by most in silico tools. In addition, a de novo FEM1C mutation of the same residue p.(Asp126Val) was associated with an undiagnosed developmental disorder, and the relevant variant (FEM1CAsp126Ala) was found to be functionally compromised in vitro. Our computational analysis showed that FEM1CAsp126His hampers protein substrate binding. To further assess its pathogenicity, we used the nematode Caenorhabditis elegans. We found that the FEM-1Asp133His animals (expressing variant homologous to the FEM1C p.(Asp126Val)) had normal muscle architecture yet impaired mobility. Mutant worms were sensitive to the acetylcholinesterase inhibitor aldicarb but not levamisole (acetylcholine receptor agonist), showing that their disabled locomotion is caused by synaptic abnormalities and not muscle dysfunction. In conclusion, we provide the first evidence from an animal model suggesting that a mutation in the evolutionarily conserved FEM1C Asp126 position causes a neurodevelopmental disorder in humans.

Mitochondrial dynamics is essential for maintaining the physiological function of the mitochondrial network, and the disorders of it leads to a variety of diseases. Our previous study identified mitochondrial dynamics controlled anti-tumor immune responses and anxiety symptoms. However, how mitochondrial dynamics affects auditory function in the inner ear remains unclear. Here, we show that deficiency of FAM73a or FAM73b, two mitochondrial outer membrane proteins that mediate mitochondrial fusion, lead to outer hair cells damage and progressive hearing loss in FVB/N mice. Abnormal mitochondrial fusion causes elevated oxidative stress and apoptosis of hair cells in early stage. Thereafter, the activation of macrophages and CD4+ T cell is found in the mutant mice with the increased expression of the inflammatory cytokines IL-12 and IFN-γ compared to control mice. Strikingly, dramatically decreased number of macrophages by LCCA treatment alleviates the hearing loss of mutant mice. Collectively, our finding highlights that FAM73a or FAM73b deficiency affects hair cells survival by disturbing mitochondrial function, and the subsequent immune response in the cochleae worsens the damage of hair cells.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease (CKD) and diabetes. Our two-stage whole-exome sequencing study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort (CRIC) and Atherosclerosis Risk in Communities (ARIC) studies (stage-1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine (TOPMed) participants (stage-2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex, and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test (SKAT-O) implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds (95% confidence interval: 33.6, 1105) of DKD compared with non-carriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% confidence interval: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.

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.

Birt-Hogg-Dubé syndrome (BHD) is an autosomal dominant disorder characterized by fibrofolliculomas, pulmonary cysts, pneumothoraces and renal cell carcinomas. Here, we reveal a novel hereditary disorder in a family with skin and mucosal lesions, extensive lipomatosis and renal cell carcinomas. The proband was initially diagnosed with BHD based on the presence of fibrofolliculomas, but no pathogenic germline variant was detected in FLCN, the gene associated with BHD. By whole exome sequencing we identified a heterozygous missense variant (p.(Cys677Tyr)) in a zinc-finger encoding domain of the PRDM10 gene which co-segregated with the phenotype in the family. We show that PRDM10Cys677Tyr loses affinity for a regulatory binding motif in the FLCN promoter, abrogating cellular FLCN mRNA and protein levels. Overexpressing inducible PRDM10Cys677Tyr in renal epithelial cells altered the transcription of multiple genes, showing overlap but also differences with the effects of knocking out FLCN. We propose that PRDM10 controls an extensive gene program and acts as a critical regulator of FLCN gene transcription in human cells. The germline variant PRDM10Cys677Tyr curtails cellular folliculin expression and underlies a distinguishable syndrome characterized by extensive lipomatosis, fibrofolliculomas and renal cell carcinomas.

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.

High-throughput circular RNA (circRNA) sequencing identified circRNA_001678 (circ_001678) as an upregulated circRNA in non-small cell lung cancer (NSCLC) tissues. Hence, the current study sought to investigate the function and the underlying mechanism of circRNA_001678 in immune escape of NSCLC. Briefly, commercially purchased NSCLC cell lines were adopted for in vitro experiment to evaluate the effects of circ_001678 over-expression or knockdown on cell biological functions, including proliferation, migration and invasive abilities. In addition, the effects of circ_001678 on the in vivo tumorigenicity ability were evaluated for verification. Accordingly, we uncovered that circ_001678 over-expression augmented NSCLC progression in vitro and enhanced tumorigenicity ability in vivo. The interaction between circ_001678 and miR-326 predicted online was verified by means of luciferase and RNA pull-down assays. Furthermore, circ_001678 could sponge miR-326 to up-regulate ZEB1. On the other hand, the tumor-promoting effects of circ_001678 could be inhibited by anti-PD-L1/PD-1 treatment. Mechanistically, circ_001678 led to the activation of the PD-1/PD-L1 pathway to promote CD8+ T cell apoptosis, thereby inducing NSCLC cell immune escape via regulation of the miR-326/ZEB1 axis. To conclude, our findings revealed that circ_001678 sponges miR-326 to up-regulate ZEB1 expression and induce the PD-1/PD-L1 pathway-dependent immune escape, thereby promoting the malignant progression of NSCLC.

All mutations in exon 9 of the Cullin3 gene associated with pseudohypoaldosteronism type II (PHA II) contribute to exon skipping to different degrees, but the specific molecular mechanism of this aberrant splicing is still unclear. The aims of this study were to investigate the regulatory mechanism underlying two synonymous splicing events, c.1221A > G (p. Glu407Glu) and c.1236G > A (p. Leu412Leu), and to discover a therapeutic strategy for correcting this aberrant splicing by targeting potential regulatory sites. Through a series of RNA pulldown, silver staining, western blotting, small interfering RNA knockdown, in vitro overexpression and single or double site-directed mutagenesis experiments, we first explored the pathogenesis of exon 9 skipping caused by mutations in the CUL3 gene and verified that the main splicing regulators associated with the synonymous c.1221A > G and c.1236G > A mutations were heterogeneous nuclear ribonucleoproteins. In addition, we verified that introducing another synonymous mutation, c.1224A > G (A18G), significantly rescued the abnormal splicing caused by c.1221A > G and c.1236G > A, highlighting the therapeutic potential for the treatment of PHA II.

To nominate novel disease genes for obesity and type 2 diabetes (T2D), we recently generated two mouse backcross populations of the T2D-susceptible New Zealand Obese (NZO/HI) mouse strain and two genetically different, lean and T2D-resistant strains, 129P2/OlaHsd and C3HeB/FeJ. Comparative linkage analysis of our two female backcross populations identified seven novel body fat-associated quantitative trait loci (QTL). Only the locus Nbw14 (NZO body weight on chromosome 14) showed linkage to obesity-related traits in both backcross populations, indicating that the causal gene variant is likely specific for the NZO strain as NZO allele carriers in both crosses displayed elevated body weight and fat mass. To identify candidate genes for Nbw14, we used a combined approach of gene expression and haplotype analysis to filter for NZO-specific gene variants in gonadal white adipose tissue, defined as the main QTL-target tissue. Only two genes, Arl11 and Sgcg, fulfilled our candidate criteria. In addition, expression QTL analysis revealed cis-signals for both genes within the Nbw14 locus. Moreover, retroviral overexpression of Sgcg in 3T3-L1 adipocytes resulted in increased insulin-stimulated glucose uptake. In humans, mRNA levels of SGCG correlated with body mass index and body fat mass exclusively in diabetic subjects, suggesting that SGCG may present a novel marker for metabolically unhealthy obesity. In conclusion, our comparative-cross analysis could substantially improve the mapping resolution of the obesity locus Nbw14. Future studies will throw light on the mechanism by which Sgcg may protect from the development of obesity.

Congenital hypothyroidism due to thyroid dysgenesis (TD), presented as thyroid aplasia, hypoplasia or ectopia, is one of the most prevalent rare diseases with an isolated organ malformation. The pathogenesis of TD is largely unknown, although a genetic predisposition has been suggested. We performed a genome-wide association study (GWAS) with 142 Japanese TD cases and 8380 controls and found a significant locus at 2q33.3 (top single nucleotide polymorphism, rs9789446: P = 4.4 × 10-12), which was replicated in a German patient cohort (P = 0.0056). A subgroup analysis showed that rs9789446 confers a risk for thyroid aplasia (per allele odds ratio = 3.17) and ectopia (3.12) but not for hypoplasia. Comprehensive epigenomic characterization of the 72-kb disease-associated region revealed that it was enriched for active enhancer signatures in human thyroid. Analysis of chromosome conformation capture data showed long-range chromatin interactions of this region with promoters of two genes, FZD5 and CCNYL1, mediating Wnt signaling. Moreover, rs9789446 was found to be a thyroid-specific quantitative trait locus, adding further evidence for a cis-regulatory function of this region in thyroid tissue. Specifically, because the risk rs9789446 allele is associated with increased thyroidal expression of FDZ5 and CCNYL1 and given the recent demonstration of perturbed early thyroid development following overactivation of Wnt signaling in zebrafish embryos, an enhanced Wnt signaling in risk allele carriers provides a biologically plausible TD mechanism. In conclusion, our work found the first risk locus for TD, exemplifying that in rare diseases with relatively low biological complexity, GWAS may provide mechanistic insights even with a small sample size.

The ClinGen malignant hyperthermia susceptibility (MHS) variant curation expert panel specified the American College of Medical Genetics and Genomics/Association of Molecular Pathologists (ACMG/AMP) criteria for RYR1-related MHS and a pilot analysis of 84 variants was published. We have now classified an additional 251 variants for RYR1-related MHS according to current ClinGen standards and updated the criteria where necessary. Criterion PS4 was modified such that individuals with multiple RYR1 variants classified as pathogenic (P), likely pathogenic (LP), or variant of uncertain significance (VUS) were not considered as providing evidence for pathogenicity. Criteria PS1 and PM5 were revised to consider LP variants at the same amino-acid residue as providing evidence for pathogenicity at reduced strength. Finally, PM1 was revised such that if PS1 or PM5 are used PM1, if applicable, should be downgraded to supporting. Of the 251 RYR1 variants, 42 were classified as P/LP, 16 as B/LB, and 193 as VUS. The primary driver of 175 VUS classifications was insufficient evidence supporting pathogenicity, rather than evidence against pathogenicity. Functional data supporting PS3/BS3 was identified for only 13 variants. Based on the posterior probabilities of pathogenicity and variant frequencies in gnomAD, we estimated the prevalence of individuals with RYR1-related MHS pathogenic variants to be between 1/300 and 1/1075, considerably higher than current estimates. We have updated ACMG/AMP criteria for RYR1-related MHS and classified 251 variants. We suggest that prioritization of functional studies is needed to resolve the large number of VUS classifications and allow for appropriate risk assessment. RYR1-related MHS pathogenic variants are likely to be more common than currently appreciated.

Pathogenic variants in ANKRD11 or microdeletions at 16q24.3 are the cause of KBG syndrome (KBGS), a neurodevelopmental syndrome characterized by intellectual disability, dental and skeletal anomalies, and characteristic facies. The ANKRD11 gene encodes the ankyrin repeat-containing protein 11A transcriptional regulator, expressed in the brain and implicated in neural development. Syndromic conditions caused by pathogenic variants in epigenetic regulatory genes show unique patterns of DNA methylation (DNAm) in peripheral blood termed DNAm signatures. Given ANKRD11's role in chromatin modification, we tested whether pathogenic ANKRD11 variants underlying KBGS are associated with a DNAm signature. We profiled whole-blood DNAm in 21 individuals with ANKRD11 variants, 2 individuals with microdeletions at 16q24.3, and 28 typically developing individuals, using Illumina's Infinium EPIC array. We identified 95 differentially methylated CpG sites that distinguished individuals with KBGS and pathogenic variants in ANKRD11 (n = 14) from typically developing controls (n = 28) and was also validated in seven additional individuals with pathogenic ANKRD11 variants. We generated a machine learning model from the KBGS DNAm signature and classified the DNAm profiles of 4 individuals with variants of uncertain significance (VUS) in ANKRD11. We identified an intermediate classification score for an inherited missense variant transmitted from a clinically unaffected mother to her affected child. In conclusion, we show the DNAm profile of two individuals with 16q24.3 microdeletions, was indistinguishable from the DNAm profile of individuals with pathogenic variants in ANKRD11 and we demonstrate the diagnostic utility of the new KBGS signature by classifying the DNAm profiles of individuals with VUS in ANKRD11.

Nineteen genetic susceptibility loci for esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) have been identified through genome-wide association studies (GWAS). Clinical translation of such discoveries, however, has been hindered by the slow pace of discovery of functional/causal variants and gene targets at these loci. We previously developed a systematic informatics pipeline to prioritize candidate functional variants using functional potential scores, applied the pipeline to select high-scoring BE/EAC risk loci and validated a functional variant at chr19p13.11 (rs10423674). Here, we selected two additional prioritized loci for experimental interrogation: chr3p13/rs1522552 and chr8p23.1/rs55896564. Candidate enhancer regions encompassing these variants were evaluated using luciferase reporter assays in two EAC cell lines. One of the two regions tested exhibited allele-specific enhancer activity - 8p23.1/rs55896564. CRISPR-mediated deletion of the putative enhancer in EAC cell lines correlated with reduced expression of three candidate gene targets: B lymphocyte kinase (BLK), nei like DNA glycosylase 2 (NEIL2) and cathepsin B (CTSB). Expression quantitative trait locus (eQTL) mapping in normal esophagus and stomach revealed strong associations between the BE/EAC risk allele at rs55896564 (G) and lower expression of CTSB, a protease gene implicated in epithelial wound repair. These results further support the utility of functional potential scores for GWAS variant prioritization, and provide the first experimental evidence of a functional variant and risk enhancer at the 8p23.1 GWAS locus. Identification of CTSB, BLK and NEIL2 as candidate gene targets suggests that altered expression of these genes may underlie the genetic risk association at 8p23.1 with BE/EAC.

Despite early interest, the evidence linking fatty acids to cardiovascular diseases (CVDs) remains controversial. We used Mendelian randomization to explore the involvement of polyunsaturated (PUFA) and monounsaturated (MUFA) fatty acids biosynthesis in the etiology of several CVD endpoints in up to 1 153 768 European (maximum 123 668 cases) and 212 453 East Asian (maximum 29 319 cases) ancestry individuals. As instruments, we selected single nucleotide polymorphisms mapping to genes with well-known roles in PUFA (i.e. FADS1/2 and ELOVL2) and MUFA (i.e. SCD) biosynthesis. Our findings suggest that higher PUFA biosynthesis rate (proxied by rs174576 near FADS1/2) is related to higher odds of multiple CVDs, particularly ischemic stroke, peripheral artery disease and venous thromboembolism, whereas higher MUFA biosynthesis rate (proxied by rs603424 near SCD) is related to lower odds of coronary artery disease among Europeans. Results were unclear for East Asians as most effect estimates were imprecise. By triangulating multiple approaches (i.e. uni-/multi-variable Mendelian randomization, a phenome-wide scan, genetic colocalization and within-sibling analyses), our results are compatible with higher low-density lipoprotein (LDL) cholesterol (and possibly glucose) being a downstream effect of higher PUFA biosynthesis rate. Our findings indicate that PUFA and MUFA biosynthesis are involved in the etiology of CVDs and suggest LDL cholesterol as a potential mediating trait between PUFA biosynthesis and CVDs risk.

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate 'on average', creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an 'arching' profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Cyclin-dependent kinase-like 5 (Cdkl5) deficiency disorder (CDD) is a severe neurodevelopmental condition caused by mutations in the X-linked Cdkl5 gene. CDD is characterized by early-onset seizures in the first month of life, intellectual disability, motor and social impairment. No effective treatment is currently available and medical management is only symptomatic and supportive. Recently, mouse models of Cdkl5 disorder have demonstrated that mice lacking Cdkl5 exhibit autism-like phenotypes, hyperactivity and dysregulations of the arousal system, suggesting the possibility to use these features as translational biomarkers. In this study, we tested Cdkl5 male and female mutant mice in an appetitive operant conditioning chamber to assess cognitive and motor abilities, and performed pupillometry to assess the integrity of the arousal system. Then, we evaluated the performance of artificial intelligence models to classify the genotype of the animals from the behavioral and physiological phenotype. The behavioral results show that CDD mice display impulsivity, together with low levels of cognitive flexibility and perseverative behaviors. We assessed arousal levels by simultaneously recording pupil size and locomotor activity. Pupillometry reveals in CDD mice a smaller pupil size and an impaired response to unexpected stimuli associated with hyperlocomotion, demonstrating a global defect in arousal modulation. Finally, machine learning reveals that both behavioral and pupillometry parameters can be considered good predictors of CDD. Since early diagnosis is essential to evaluate treatment outcomes and pupillary measures can be performed easily, we proposed the monitoring of pupil size as a promising biomarker for CDD.

The mitochondrial DNA mutation m.9032 T > C was previously identified in patients presenting with NARP (Neuropathy Ataxia Retinitis Pigmentosa). Their clinical features had a maternal transmission and patient's cells showed a reduced oxidative phosphorylation capacity, elevated reactive oxygen species (ROS) production and hyperpolarization of the mitochondrial inner membrane, providing evidence that m.9032 T > C is truly pathogenic. This mutation leads to replacement of a highly conserved leucine residue with proline at position 169 of ATP synthase subunit a (L169P). This protein and a ring of identical c-subunits (c-ring) move protons through the mitochondrial inner membrane coupled to ATP synthesis. We herein investigated the consequences of m.9032 T > C on ATP synthase in a strain of Saccharomyces cerevisiae with an equivalent mutation (L186P). The mutant enzyme assembled correctly but was mostly inactive as evidenced by a > 95% drop in the rate of mitochondrial ATP synthesis and absence of significant ATP-driven proton pumping across the mitochondrial membrane. Intragenic suppressors selected from L186P yeast restoring ATP synthase function to varying degrees (30-70%) were identified at the original mutation site (L186S) or in another position of the subunit a (H114Q, I118T). In light of atomic structures of yeast ATP synthase recently described, we conclude from these results that m.9032 T > C disrupts proton conduction between the external side of the membrane and the c-ring, and that H114Q and I118T enable protons to access the c-ring through a modified pathway.

OBJECTIVE: The objective of this study was to discuss the mechanism of artesunate (ART) in improving cartilage damage in osteoarthritis (OA) by regulating the expression levels of MTA1, LXA4, and the downstream JAK2/STAT3 signaling pathway. METHODS: The OA model in vitro was constructed by stimulating chondrocytes for 24 h with 10 ng/mL IL-1β, and cell proliferation and apoptosis, expression levels of Aggrecan, MTA1, LXA4, MMP3, MMP13, and Collagen II, and inflammatory cytokines in the culture supernatants were examined. Histopathological changes, inflammatory response, and chondrocyte apoptosis of the cartilage tissues of OA mice were performed. RESULTS: In vitro cell experiments, ART enhanced cell proliferation capacity, accompanied by decreased apoptosis rate, decreased expression of MMP-3 and MMP-13, elevated expression of Collagen II and Aggrecan, as well as reduced levels of IL-6 and TNF-α in the cell supernatant. ART also ameliorated IL-1β-induced chondrocyte damage by upregulating MTA1. The LXA4 promoter region had two potential binding sites for MTA1. There was a positive correlation between MTA1 and LXA4. MTA1 enhanced the expression of LXA4 through transcription and blocked the activation of the JAK2/STAT3 signaling pathway. In vivo animal model experiments further showed that ART treatment alleviated cartilage tissue damage in OA model mice by up-regulating MTA1. CONCLUSION: Our study demonstrates that ART improves the cartilage damage of OA by up-regulating MTA1 expression and promoting the transcriptional activation of LXA4, and further blocking the JAK2/STAT3 signaling pathway.

XIST is a 17-19 kb long non-coding RNA critical for X-chromosome inactivation. Tandem repeats within the RNA serve as functional domains involved in the cis-limited recruitment of heterochromatic changes and silencing. To explore the sufficiency of these domains while generating a functional miniXIST for targeted silencing approaches, we tested inducible constructs integrated into 8p in a male cell line. Previous results suggested silencing could be accomplished with a transgene comprised of the repeat A which is highly conserved and critical for silencing; the repeat F which overlaps regulatory elements and the repeat E which contributes to XIST localization by binding proteins such as CIZ1. As PRC1 is recruited through HNRNPK binding of repeats B-C-D we included a second 'miniXIST' comprising AFE with the mouse PID, a 660-nucleotide region known to recruit PRC1. Silencing of an adjacent gene was possible with and without PID, however, silencing more distally required the addition of PID. The recruitment of heterochromatic marks, evaluated by IF combined with RNA FISH, revealed that the AFE domains were sufficient only for CIZ1 recruitment. However, miniXIST transgene recruited all marks, albeit not to full XIST levels. The ability of the PID domain to facilitate silencing and heterochromatic marks recruitment was unexpected, and inhibition of PRC1 suggested that many of these are PRC1-independent. These results suggest that the addition of this small region allowed the partial recruitment of all the changes induced by a full XIST, demonstrating the feasibility of finding a minimal functional XIST.

Fukuyama congenital muscular dystrophy (FCMD) is an autosomal recessive disorder caused by fukutin (FKTN) gene mutations. FCMD is the second most common form of childhood muscular dystrophy in Japan, and the most patients possess a homozygous retrotransposal SINE-VNTR-Alu insertion in the 3'-untranslated region of FKTN. A deep-intronic variant (DIV) was previously identified as the second most prevalent loss-of-function mutation in Japanese patients with FCMD. The DIV creates a new splicing donor site in intron 5 that causes aberrant splicing and the formation of a 64-base pair pseudoexon in the mature mRNA, resulting in a truncated nonfunctional protein. Patients with FCMD carrying the DIV present a more severe symptoms, and currently, there is no radical therapy available for this disorder. In the present study, we describe in vitro evaluation of antisense oligonucleotide mediated skipping of pseudoexon inclusion and restoration of functional FKTN protein. A total of 16 19-26-mer antisense oligonucleotide sequences were designed with a 2'-O-methyl backbone and were screened in patient-derived fibroblasts, lymphoblast cells, and minigene splice assays. One antisense oligonucleotide targeting the exonic splice enhancer region significantly induced pseudoexon skipping and increased the expression of normal mRNA. It also rescued FKTN protein production in lymphoblast cells and restored functional O-mannosyl glycosylation of alpha-dystroglycan in patient-derived myotubes. Based on our results, ASO-based splicing correction should be investigated further as a potential treatment for patients with FCMD carrying the DIV.

Mitochondrial involvement in neurodegenerative diseases is widespread and multifactorial. Targeting mitochondrial pathology is therefore of interest. The recent development of bioactive molecules that modulate mitochondrial dynamics (fusion, fission and motility) offers a new therapeutic approach for neurodegenerative diseases with either indirect or direct mitochondrial involvement. Here, we asked: 1. Can enhanced mitochondrial fusion and motility improve secondary mitochondrial pathology in SOD1 mutant amyotrophic lateral sclerosis (ALS)? And: 2. What is the impact of enhancing mitochondria fitness on in vivo manifestations of SOD1 mutant ALS? We observed that small molecule mitofusin activators corrected mitochondrial fragmentation, depolarization and dysmotility in genetically diverse ALS patient reprogrammed motor neurons and fibroblasts, and in motor neurons, sensory neurons and fibroblasts from SOD1 G93A mice. Continuous, but not intermittent, pharmacologic mitofusin activation delayed phenotype progression and lethality in SOD1 G93A mice, reducing neuron loss and improving neuromuscular connectivity. Mechanistically, mitofusin activation increased mitochondrial motility, fitness and residency within neuromuscular synapses, reduced mitochondrial ROS production, and diminished apoptosis in SOD1 mutant neurons. These benefits were accompanied by improved mitochondrial respiratory coupling, despite characteristic SOD1 mutant ALS-associated downregulation of mitochondrial respiratory complexes. Targeting mitochondrial dysdynamism is a promising approach to alleviate pathology caused by secondary mitochondrial dysfunction in some neurodegenerative diseases.

Charcot-Marie-Tooth disease is an inherited peripheral neuropathy that is clinically and genetically heterogenous. Mutations in IGHMBP2, a ubiquitously expressed DNA/RNA helicase, have been shown to cause the infantile motor neuron disease spinal muscular atrophy with respiratory distress type 1 (SMARD1), and, more recently, juvenile-onset Charcot-Marie-Tooth disease Type 2S (CMT2S). Using CRISPR-cas9 mutagenesis we developed the first mouse models of CMT2S (p.Glu365del (E365del) and p.Tyr918Cys (Y918C)). E365del is the first CMT2S mouse model to be discovered and Y918C is the first human CMT2S allele knock-in model. Phenotypic characterization of the homozygous models found progressive peripheral motor and sensory axonal degeneration. Neuromuscular and locomotor assays indicate that both E365del and Y918C mice have motor deficits, while neurobehavioral characterization of sensory function found that E365del mutants have mechanical allodynia. Analysis of femoral motor and sensory nerves identified axonal degeneration, which does not impact nerve conduction velocities in E365del mice, but does in the Y918C model. Based on these results, the E365del mutant mouse, as well as the human allele knock-in, Y918C, represent mouse models with the hallmark phenotypes of CMT2S, which will be critical for understanding the pathogenic mechanisms of IGHMBP2. These mice will complement existing Ighmbp2 alleles modeling SMARD1 to help understand the complex phenotypic and genotypic heterogeneity that is observed in patients with IGHMBP2 variants.

G4C2 repeat expansion in C9orf72 causes the most common familial frontotemporal dementia and amyotrophic lateral sclerosis (C9FTD/ALS). The pathogenesis includes haploinsufficiency of C9orf72, which forms a protein complex with Smcr8, as well as G4C2 repeat-induced gain-of-function including toxic dipeptide repeats (DPRs). The key in vivo disease-driving mechanisms and how loss- and gain-of-function interplay remain poorly understood. Here we identified a lysosome-ribosome biogenesis circuit dysregulation as an early and key disease mechanism using a physiologically relevant mouse model with combined loss- and gain-of-function across the ageing process. C9orf72 deficiency exacerbates FTD/ALS-like pathologies and behaviors in C9ORF72 bacterial artificial chromosome (C9-BAC) mice with G4C2 repeats under endogenous regulatory elements from patients. Single nucleus RNA sequencing (snRNA-seq) and bulk RNA-seq revealed that C9orf72 depletion disrupts lysosomes in neurons and leads to transcriptional dysregulation of ribosomal protein (RP) genes, which are likely due to the proteotoxic stress response and resemble ribosomopathy defects. Importantly, ectopic expression of C9orf72 or its partner Smcr8 in C9FTD/ALS mutant mice promotes lysosomal functions and restores ribosome biogenesis gene transcription, resulting in the mitigation of DPR accumulation, neurodegeneration, as well as FTD/ALS-like motor and cognitive behaviors. Therefore, we conclude that loss- and gain-of-function crosstalk in C9FTD/ALS converges on neuronal dysregulation of a lysosome-ribosome biogenesis circuit leading to the proteotoxicity, neurodegeneration, and behavioral defects.

Autosomal dominant polycystic kidney disease (ADPKD) is driven by mutations in the PKD1 and PKD2 genes and is characterized by renal cyst formation, inflammation and fibrosis. Forkhead box protein M1 (FoxM1), a transcription factor of the Forkhead box (Fox) protein super family, has been reported to promote tumor formation, inflammation and fibrosis in many organs. However, the role and mechanism of FoxM1 in regulation of ADPKD progression is still poorly understood. Here, we show that FoxM1 is an important regulator of cyst growth in ADPKD. FoxM1 is upregulated in cyst-lining epithelial cells in Pkd1 mutant mouse kidneys and human ADPKD kidneys. FoxM1 promotes cystic renal epithelial cell proliferation by increasing the expression of Akt and Stat3 and the activation of ERK and Rb. FoxM1 also regulates cystic renal epithelial cell apoptosis through NF-κB signaling pathways. In addition, FoxM1 regulates the recruitment and retention of macrophages in Pkd1 mutant mouse kidneys, a process that is associated with FoxM1-mediated upregulation of monocyte chemotactic protein 1 (MCP-1). Targeting FoxM1 with its specific inhibitor, FDI-6, delays cyst growth in rapidly progressing and slowly progressing Pkd1 mutant mouse kidneys. This study suggests that FoxM1 is a central and upstream regulator of ADPKD pathogenesis and provides a rationale for targeting FoxM1 as a therapeutic strategy for ADPKD treatment.

Nemaline myopathy 8 (NEM8) is typically a severe autosomal recessive disorder associated with variants in the kelch-like family member 40 gene (KLHL40). Common features include fetal akinesia, fractures, contractures, dysphagia, respiratory failure, and neonatal death. Here, we describe a 26-year-old man with relatively mild NEM8. He presented with hypotonia and bilateral femur fractures at birth, later developing bilateral Achilles' contractures, scoliosis, and elbow and knee contractures. He had walking difficulties throughout childhood and became wheelchair bound from age 13 after prolonged immobilisation. Muscle MRI at age 13 indicated prominent fat replacement in his pelvic girdle, posterior compartments of thighs, and vastus intermedius. Muscle biopsy revealed nemaline bodies and intranuclear rods. RNA sequencing and western blotting of patient skeletal muscle indicated significant reduction in KLHL40 mRNA and protein respectively. Using gene panel screening, exome sequencing and RNA sequencing, we identified compound heterozygous variants in KLHL40; a truncating 10.9 kb deletion in trans with a likely pathogenic variant (c.*152G > T) in the 3' untranslated region (UTR). Computational tools SpliceAI and Introme predicted the c.*152G > T variant created a cryptic donor splice site. RNA-seq and in vitro analyses indicated that the c.*152G > T variant induces multiple de novo splicing events that likely provoke nonsense mediated decay of KLHL40 mRNA explaining the loss of mRNA expression and protein abundance in the patient. Analysis of 3' UTR variants in ClinVar suggests variants that introduce aberrant 3' UTR splicing may be underrecognised in Mendelian disease. We encourage consideration of this mechanism during variant curation.

More than 250 million people in the world are chronically infected with hepatitis B virus (HBV), which causes serious complications. Host genetic susceptibility is essential for CHB, and our previous genome-wide association study identified a single-nucleotide polymorphism (SNP), rs1883832, in the 5' untranslated region of CD40 predisposing to chronic HBV infection, but the underlying mechanism remains undefined. This study aimed to investigate whether rs1883832 was the real functional SNP (fSNP) of CD40 and how it modulated HBV clearance in hepatocytes. We determined the fSNP of CD40 and its regulatory protein(s) using luciferase reporter assays, electrophoretic mobility shift assay (EMSA), flanking restriction enhanced pulldown (FREP) and chromatin immunoprecipitation (ChIP). The potential anti-HBV activity of CD40 and its downstream molecule BST2 was assessed in HBV-transfected and HBV-infected hepatoma cells and HBV-infected primary human hepatocytes (PHHs). Moreover, the mechanism of CD40 was investigated by mRNA sequencing, quantitative real-time polymerase chain reaction (qRT-PCR), immunofluorescence and western blot. We revealed rs1883832 as the true fSNP of CD40 and identified ANXA2 as a negative regulatory protein that preferentially bound to the risk allele T of rs1883832 and hence reduced CD40 expression. Furthermore, CD40 suppressed HBV replication and transcription in hepatocytes via activating the JAK-STAT pathway. BST2 was identified to be the key IFN-stimulated gene regulated by CD40 after activating JAK-STAT pathway. Inhibition of JAK/STAT/BST2 axis attenuated CD40-induced antiviral effect. In conclusion, a functional variant of CD40 modulates HBV clearance via regulation of the ANXA2/CD40/BST2 axis, which may shed new light on HBV personalized therapy.

BACE1 is largely expressed by neurons and is the sole β-secretase for initiating the production of neuronal β-amyloid peptides (Aβ). To fully understand the physiological functions of neuronal BACE1, we used mouse genetic approach coupled with unbiased single nucleus RNA sequencing (snRNAseq) to investigate how targeted deletion of Bace1 in neurons, driven by Thy-1-Cre recombinase, would affect functions in the nervous system. Our transcriptome results revealed that BACE1 is essential for maturation of neural precursor cells and oligodendrocytes in mice. RNA velocity analysis confirmed deficit in the trajectory of neuroblasts in reaching the immature granule neuron state in young Bace1fl/fl; Thy1-cre mice. Further analysis of differential gene expression indicated changes in genes important for SNARE signaling, tight junction signaling, synaptogenesis and insulin secretion pathways. Morphological studies revealed a hypomyelination in Bace1fl/fl;Thy1-cre sciatic nerves, but no detectable myelination changes in the corpus callosum, even though clear reduction in myelination proteins in the brain. Functional study showed reduction in long term potential, defects in synaptogenesis, and learning behavioral. Together, our results show that neuronal BACE1 is critical for optimal development of central and peripheral nervous system, and inhibition of neuronal BACE1 will result in deficits in synaptic functions and cognitive behaviors.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (COVID-19) positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (P < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (P = 2.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, P = 4.1 × 10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

Congenital hearing impairment is a genetically highly heterogenous disorder in which prompt recognition and intervention are crucial to optimize outcomes. In this study, we used exome sequencing to investigate a large consanguineous Pakistani family with eight affected individuals showing bilateral severe-to-profound hearing impairment. This identified a homozygous splice region variant in STX4 (c.232 + 6 T > C), which causes exon skipping and a frameshift, that segregated with hearing impairment (two-point LOD score = 5.9). STX4, a member of the syntaxin family, is a component of the SNARE machinery involved in several vesicle transport and recycling pathways. In silico analysis showed that murine orthologue Stx4a is highly and widespread expressed in the developing and adult inner ear. Immunofluorescent imaging revealed localization of STX4A in the cell body, cell membrane and stereocilia of inner and outer hair cells. Furthermore, a morpholino based knockdown of stx4 in zebrafish showed an abnormal startle response, morphological and developmental defects, and a disrupted mechanotransduction function in neuromast hair cells measured via FM1-43 uptake. Our findings indicate that STX4 dysfunction leads to hearing impairment in humans and zebrafish and supports the evolutionary conserved role of STX4 in inner ear development and hair cell functioning.

Mutations in the fibrillin-1 (FBN1) gene are responsible for the autosomal dominant form of geleophysic dysplasia (GD), which is characterized by short stature and extremities, thick skin and cardiovascular disease. All known FBN1 mutations in patients with GD are localized within the region encoding the transforming growth factor-β binding protein-like 5 (TB5) domain of this protein. Herein, we generated a knock-in mouse model, Fbn1Y1698C by introducing the p.Tyr1696Cys mutation from a patient with GD into the TB5 domain of murine Fbn1 to elucidate the specific role of this domain in endochondral ossification. We found that both Fbn1Y1698C/+ and Fbn1Y1698C/Y1698C mice exhibited a reduced stature reminiscent of the human GD phenotype. The Fbn1 point mutation introduced in these mice affected the growth plate formation owing to abnormal chondrocyte differentiation such that mutant chondrocytes failed to establish a dense microfibrillar network composed of FBN1. This original Fbn1 mutant mouse model offers new insight into the pathogenic events underlying GD. Our findings suggest that the etiology of GD involves the dysregulation of the extracellular matrix composed of an abnormal FBN1 microfibril network impacting the differentiation of the chondrocytes.

Fear and anxiety have proven to be essential during the evolutionary process. However, the mechanisms involved in recognizing and categorizing threat probability (i.e. low to high) to elicit the appropriate defensive behavior are yet to be determined. In this study, we investigated the cerebellar contribution in evoking appropriate defensive escape behavior using a purely cerebellar, neurodegenerative mouse model for spinocerebellar ataxia type 6 which is caused by an expanded CAG repeat in exon 47 of the P/Q type calcium channel α1A subunit. These mice overexpress the carboxy terminus (CT) of the P/Q type calcium channel containing an expanded 27 CAG repeat specifically in cerebellar Purkinje cells (CT-longQ27PC). We found that our CT-longQ27PC mice exhibit anxiolytic behavior in the open field, elevated plus maze and light/dark place preference tests, which could be recovered with more threatening conditions such as brighter lighting, meowing sounds and an ultrasound repellent. Their innate fear to find safety in the Barnes maze and visual cliff tests was also diminished with subsequent trials, which could be partially recovered with an ultrasound repellent in the Barnes maze. However, under higher threat conditions such as in the light/dark place preference with ultrasound repellent and in the looming tests, CT-longQ27PC mice responded with higher defensive escape behaviors as controls. Moreover, CT-longQ27PC mice displayed increased levels of CT-labeled aggregates compared with controls. Together these data suggest that cerebellar degeneration by overexpression of CT-longQ27PC is sufficient to impair defensive escape responses in those mice.

The D620N mutation in vacuolar protein sorting protein 35 (VPS35) gene has been identified to be linked to late onset familial Parkinson disease (PD). However, the pathophysiological roles of VPS35-D620N in PD remain unclear. Here, we generated the transgenic Caenorhabditis elegans overexpressing either human wild type or PD-linked mutant VPS35-D620N in neurons. C. elegans expressing VPS35-D620N, compared with non-transgenic controls, showed movement disorders and dopaminergic neuron loss. VPS35-D620N worms displayed more swimming induced paralysis but showed no defects in BSR assays, thus indicating the disruption of dopamine (DA) recycling back inside neurons. Moreover, VPS35 formed a protein interaction complex with DA transporter (DAT), RAB5, RAB11 and FAM21. In contrast, the VPS35-D620N mutant destabilized these interactions, thus disrupting DAT transport from early endosomes to recycling endosomes, and decreasing DAT at the cell surface. These effects together increased DA in synaptic clefts, and led to dopaminergic neuron degeneration and motor dysfunction. Treatment with reserpine significantly decreased the swimming induced paralysis in VPS35-D620N worms, as compared with vehicle treated VPS35-D620N worms. Our studies not only provide novel insights into the mechanisms of VPS35-D620N-induced dopaminergic neuron degeneration and motor dysfunction via disruption of DAT function and the DA signaling pathway but also indicate a potential strategy to treat VPS35-D620N-related PD and other disorders.

In vitro fertilization (IVF) is associated with DNA methylation abnormalities and a higher incidence of adverse pregnancy outcomes. However, which exposure(s), among the many IVF interventions, contributes to these outcomes remains unknown. Frozen embryo transfer (ET) is increasingly utilized as an alternative to fresh ET, but reports suggest a higher incidence of pre-eclampsia and large for gestational age infants. This study examines DNA methylation in human placentas using the 850K Infinium MethylationEPIC BeadChip array obtained after 65 programmed frozen ET cycles, 82 fresh ET cycles and 45 unassisted conceptions. Nine patients provided placentas following frozen and fresh ET from consecutive pregnancies for a paired subgroup analysis. In parallel, eight mouse placentas from fresh and frozen ET were analyzed using the Infinium Mouse Methylation BeadChip array. Human and mouse placentas were significantly hypermethylated after frozen ET compared with fresh. Paired analysis showed similar trends. Sex-specific analysis revealed that these changes were driven by male placentas in humans and mice. Frozen and fresh ET placentas were significantly different from controls, with frozen samples hypermethylated compared with controls driven by males and fresh samples being hypomethylated compared with controls, driven by females. Sexually dimorphic epigenetic changes could indicate differential susceptibility to IVF-associated perturbations, which highlights the importance of sex-specific evaluation of adverse outcomes. Similarities between changes in mice and humans underscore the suitability of the mouse model in evaluating how IVF impacts the epigenetic landscape, which is valuable given limited access to human tissue and the ability to isolate specific interventions in mice.

Genome-wide association studies (GWAS) of Crohn's disease (CD) in European and leprosy in Chinese population have shown that CD and leprosy share genetic risk loci. As these shared loci were identified through cross-comparisons across different ethnic populations, we hypothesized that meta-analysis of GWAS on CD and leprosy in East Asian populations would increase power to identify additional shared loci. We performed a cross-disease meta-analysis of GWAS data from CD (1621 cases and 4419 controls) and leprosy (2901 cases 3801 controls) followed by replication in additional datasets comprising 738 CD cases and 488 controls and 842 leprosy cases and 925 controls. We identified one novel locus at 7p22.3, rs77992257 in intron 2 of ADAP1, shared between CD and leprosy with genome-wide significance (P = 3.80 × 10-11) and confirmed 10 previously established loci in both diseases: IL23R, IL18RAP, IL12B, RIPK2, TNFSF15, ZNF365-EGR2, CCDC88B, LACC1, IL27, NOD2. Phenotype variance explained by the polygenic risk scores derived from Chinese leprosy data explained up to 5.28% of variance of Korean CD, supporting similar genetic structures between the two diseases. Although CD and leprosy shared a substantial number of genetic susceptibility loci in East Asians, the majority of shared susceptibility loci showed allelic effects in the opposite direction. Investigation of the genetic correlation using cross-trait linkage disequilibrium score regression also showed a negative genetic correlation between CD and leprosy (rg [SE] = -0.40[0.13], P = 2.6 × 10-3). These observations implicate the possibility that CD might be caused by hyper-sensitive reactions toward pathogenic stimuli.

Mental disorders present a global health concern and have limited treatment options. In today's medical practice, medications such as antidepressants are prescribed not only for depression but also for conditions such as anxiety and attention deficit hyperactivity disorder (ADHD). Therefore, identifying gene targets for specific disorders is important and offers improved precision. In this study, we performed a genetic analysis of six common mental disorders-ADHD, anxiety, depression, delays in mental development, intellectual disabilities (IDs) and speech/language disorder-in the ethnic minority of African Americans (AAs) using whole genome sequencing (WGS). WGS data were generated from blood-derived DNA from 4178 AA individuals, including 1384 patients with the diagnosis of at least one mental disorder. Mutation burden analysis was applied based on rare and deleterious mutations in the AA population between cases and controls, and further analyzed in the context of patients with single mental disorder diagnosis. Certain genes uncovered demonstrated significant P-values in mutation burden analysis. In addition, exclusive recurrences in specific type of disorder were scanned through gene-drug interaction databases to assess for availability of potential medications. We uncovered 15 genes harboring deleterious mutations, including 3-Hydroxy-3-Methylglutaryl-CoA Reductase (HMGCR) and Uronyl 2-Sulfotransferase (UST) for ADHD; Farnesyltransferase, CAAX Box, Beta (FNTB) for anxiety; Xin Actin Binding Repeat Containing 2 (XIRP2), Natriuretic Peptide C (NPPC), Serine/Threonine Kinase 33 (STK33), Pannexin 1 (PANX1) and Neurotensin (NTS) for depression; RUNX Family Transcription Factor 3 (RUNX3), Tachykinin Receptor 1 (TACR1) and NADH:Ubiquinone Oxidoreductase Core Subunit S7 (NDUFS7) for delays in mental development; Hepsin (HPN) for ID and Collagen Type VI Alpha 3 Chain (COL6A3), Damage Specific DNA Binding Protein 1 (DDB1) and NADH:Ubiquinone Oxidoreductase Subunit A11 (NDUFA11) for speech/language disorder. Taken together, we have established critical insights into the development of new precision medicine approaches for mental disorders in AAs.

The aim of this study is to investigate if extracellular vesicles (EVs) from bone marrow mesenchymal stem cells (BMSCs) deliver microRNA (miR)-331-3p to regulate LIM zinc finger domain containing 2 (LIMS2) methylation in cervical cancer cells. Cervical cancer cells were incubated with EVs from BMSCs with altered expression of miR-331-3p, DNA methyltransferase 3 alpha (DNMT3A) or/and LIMS2 and then subjected to 5-ethynyl-2'-deoxyuridine, Transwell, flow cytometry and western blotting analyses. Dual-luciferase reporter assay was conducted to verify the binding between miR-331-3p and DNMT3A. A xenograft model was established to evaluate the effect of BMSC-derived EV-miR-331-3p on cervical tumor growth. miR-331-3p was lowly and DNMT3A was highly expressed in cervical cancer. BMSC-derived EVs delivered miR-331-3p to control the behaviors of cervical cancer cells. miR-331-3p inhibited the expression of DNMT3A by binding DNMT3A mRNA. DNMT3A promoted LIMS2 methylation and reduced the expression of LIMS2. Overexpression of DNMT3A or silencing of LIMS2 in BMSCs counteracted the tumor suppressive effects of miR-331-3p. BMSC-derived EV-miR-331-3p also inhibited the growth of cervical tumors in vivo. BMSC-derived EVs alleviate cervical cancer partially by delivering miR-331-3p to reduce DNMT3A-dependent LIMS2 methylation in tumor cells.

Peroxiredoxin 3 (PRDX3) encodes a mitochondrial antioxidant protein, which is essential for the control of reactive oxygen species homeostasis. So far, PRDX3 mutations are involved in mild-to-moderate progressive juvenile onset cerebellar ataxia. We aimed to unravel the molecular bases underlying the disease in an infant suffering from cerebellar ataxia that started at 19 months old and presented severe cerebellar atrophy and peripheral neuropathy early in the course of disease. By whole exome sequencing, we identified a novel homozygous mutation, PRDX3 p.D163E, which impaired the mitochondrial ROS defense system. In mouse primary cortical neurons, the exogenous expression of PRDX3 p.D163E was reduced and triggered alterations in neurite morphology and in mitochondria. Mitochondrial computational parameters showed that p.D163E led to serious mitochondrial alterations. In transfected HeLa cells expressing the mutation, mitochondria accumulation was detected by correlative light electron microscopy. Mitochondrial morphology showed severe changes, including extremely damaged outer and inner membranes with a notable cristae disorganization. Moreover, spherical structures compatible with lipid droplets were identified, which can be associated with a generalized response to stress and can be involved in the removal of unfolded proteins. In the patient's fibroblasts, PRDX3 expression was nearly absent. The biochemical analysis suggested that the mutation p.D163E would result in an unstable structure tending to form aggregates that trigger unfolded protein responses via mitochondria and endoplasmic reticulum. Altogether, our findings broaden the clinical spectrum of the recently described PRDX3-associated neurodegeneration and provide new insight into the pathological mechanisms underlying this new form of cerebellar ataxia.

RATIONALE: Genetic 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. OBJECTIVES: We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program. METHODS: Using 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 RESULTS: In 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 < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability. CONCLUSIONS: Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

CTR9 is one of five genes that form the PAF1 complex, which binds to RNA polymerase II and plays critical roles in transcriptional elongation and transcription-coupled histone modifications including histones H3K4me3 and H3K36me3. In this study, de novo CTR9 non-synonymous variants (p.(Glu15Asp) and p.(Pro25Arg)) were detected in two unrelated patients with macrocephaly, motor delay, and intellectual disability. A pull-down assay showed that the mutant CTR9 proteins had stronger affinities to the PAF1 protein than the wild-type protein. Functional analyses using zebrafish showed that the knockout of the ctr9 gene caused motor defects and enlargement of the telencephalon, which is homologous to the mammalian cerebrum. The rescue experiment, in which the human CTR9 mutants were introduced into ctr9-knockout zebrafish, failed to maintain the swimming posture of the ctr9-knockout fish, suggesting that the human CTR9 mutant proteins do not function normally in vivo. In addition, the overexpression of human CTR9 mutant mRNA caused telencephalon enlargement in zebrafish larvae, suggesting that the human CTR9 mutant proteins interfered with normal endogenous CTR9 function. We concluded that the two missense variants in CTR9 (p.(Glu15Asp) and p.(Pro25Arg)) cause a new syndrome involving macrocephaly, motor delay and intellectual disability through the loss of the normal function of CTR9 and the inhibition of the normal intrinsic CTR9 function of the contralateral allele.

The basic motif-leucine zipper (bZIP) transcription factor neural retina leucine zipper (NRL) determines rod photoreceptor cell fate during retinal development, and its loss leads to cone-only retina in mice. NRL works synergistically with homeodomain protein Cone-Rod Homeobox and other regulatory factors to control the transcription of most genes associated with rod morphogenesis and functional maturation, which span over a period of several weeks in the mammalian retina. We predicted that NRL gradually establishes rod cell identity and function by temporal and dynamic regulation of stage-specific transcriptional targets. Therefore, we mapped the genomic occupancy of NRL at four stages of mouse photoreceptor differentiation by CUT&RUN analysis. Dynamics of NRL binding revealed concordance with the corresponding changes in transcriptome of the developing rods. Notably, we identified c-Jun proto-oncogene as one of the targets of NRL, which could bind to specific cis-elements in the c-Jun promoter and modulate its activity in HEK293 cells. Coimmunoprecipitation studies showed the association of NRL with c-Jun, also a bZIP protein, in transfected cells as well as in developing mouse retina. Additionally, shRNA-mediated knockdown of c-Jun in the mouse retina in vivo resulted in altered expression of almost 1000 genes, with reduced expression of phototransduction genes and many direct targets of NRL in rod photoreceptors. We propose that c-Jun-NRL heterodimers prime the NRL-directed transcriptional program in neonatal rod photoreceptors before high NRL expression suppresses c-Jun at later stages. Our study highlights a broader cooperation among cell-type restricted and widely expressed bZIP proteins, such as c-Jun, in specific spatiotemporal contexts during cellular differentiation.