Mono- and biallelic variants in ALDH18A1 cause a spectrum of human disorders associated with cutaneous and neurological findings that overlap with both cutis laxa and spastic paraplegia. ALDH18A1 encodes the bifunctional enzyme pyrroline-5-carboxylate synthetase (P5CS) that plays a role in the de novo biosynthesis of proline and ornithine. Here we characterize a previously unreported homozygous ALDH18A1 variant (p.Thr331Pro) in four affected probands from two unrelated families, and demonstrate broad-based alterations in amino acid and antioxidant metabolism. These four patients exhibit variable developmental delay, neurological deficits, and loose skin. Functional characterization of the p.Thr331Pro variant demonstrated a lack of any impact on the steady-state level of the P5CS monomer or mitochondrial localization of the enzyme, but reduced incorporation of the monomer into P5CS oligomers. Using an unlabeled NMR-based metabolomics approach in patient fibroblasts and ALDH18A1-null human embryonic kidney cells expressing the variant P5CS, we identified reduced abundance of glutamate and several metabolites derived from glutamate, including proline and glutathione. Biosynthesis of the polyamine putrescine, derived from ornithine, was also decreased in patient fibroblasts, highlighting the functional consequence on another metabolic pathway involved in antioxidant responses in the cell. RNA sequencing of patient fibroblasts revealed transcript abundance changes in several metabolic and ECM-related genes, adding further insight into pathogenic processes associated with impaired P5CS function. Together these findings shed new light on amino acid and antioxidant pathways associated with ALDH18A1-related disorders, and underscore the value of metabolomic and transcriptomic profiling to discover new pathways that impact disease pathogenesis.

The Transcriptome-Wide Association Study (TWAS) is a widely used approach integrating gene expression and GWAS data to study the role of cis-regulated gene expression (GEx) in complex traits. However, the genetic architecture of GEx varies across populations, and recent findings point to possible ancestral heterogeneity in the effects of GEx on complex traits, which may be amplified in TWAS by modeling GEx as a function of cis-eQTLs. Here, we present a novel extension to TWAS to account for heterogeneity in the effects of cis-regulated GEx which are correlated with ancestry. Our proposed Multi- Ancestry TwaS (MATS) framework jointly analyzes samples from multiple populations and distinguishes between shared, ancestry-specific, and/or subject-specific expression-trait associations. As such, MATS amplifies power to detect shared GEx associations over ancestry-stratified TWAS through increased sample sizes, and facilitates detection of genes with subgroup-specific associations which may be masked by standard TWAS. Our simulations highlight the improved Type-I error conservation and power of MATS compared to competing approaches. Our real data applications to ad case-control genotypes from the Alzheimer's Disease Sequencing Project (ADSP) and continuous phenotypes from the UK Biobank (UKBB) identify a number of unique gene-trait associations which were not discovered through standard and/or ancestry-stratified TWAS. Ultimately, these findings promote MATS as a powerful method for detecting and estimating significant gene expression effects on complex traits within multi-ancestry cohorts, and corroborates the mounting evidence for inter-population heterogeneity in gene-trait associations.

Gestational diabetes mellitus (GDM) is associated with increased risk of pregnancy complications and adverse perinatal outcomes. GDM often reoccurs and is associated with increased risk of subsequent diagnosis of type 2 diabetes (T2D). To improve our understanding of the aetiological factors and molecular processes driving the occurrence of GDM, including the extent to which these overlap with T2D pathophysiology, the GENetics of Diabetes In Pregnancy Consortium assembled genome-wide association studies of diverse ancestry in a total of 5485 women with GDM and 347 856 without GDM. Through multi-ancestry meta-analysis, we identified five loci with genome-wide significant association (P < 5 × 10-8) with GDM, mapping to/near MTNR1B (P = 4.3 × 10-54), TCF7L2 (P = 4.0 × 10-16), CDKAL1 (P = 1.6 × 10-14), CDKN2A-CDKN2B (P = 4.1 × 10-9) and HKDC1 (P = 2.9 × 10-8). Multiple lines of evidence pointed to the shared pathophysiology of GDM and T2D: (i) four of the five GDM loci (not HKDC1) have been previously reported at genome-wide significance for T2D; (ii) significant enrichment for associations with GDM at previously reported T2D loci; (iii) strong genetic correlation between GDM and T2D and (iv) enrichment of GDM associations mapping to genomic annotations in diabetes-relevant tissues and transcription factor binding sites. Mendelian randomization analyses demonstrated significant causal association (5% false discovery rate) of higher body mass index on increased GDM risk. Our results provide support for the hypothesis that GDM and T2D are part of the same underlying pathology but that, as exemplified by the HKDC1 locus, there are genetic determinants of GDM that are specific to glucose regulation in pregnancy.

In the era of personalized medicine with more and more patient-specific targeted therapies being used, we need reliable, dynamic, faster and sensitive biomarkers both to track the causes of disease and to develop and evolve therapies during the course of treatment. Metabolomics recently has shown substantial evidence to support its emerging role in disease diagnosis and prognosis. Aside from biomarkers and development of therapies, it is also an important goal to understand the involvement of mitochondrial DNA (mtDNA) in metabolic regulation, aging and disease development. Somatic mutations of the mitochondrial genome are also heavily implicated in age-related disease and aging. The general hypothesis is that an alteration in the concentration of metabolite profiles (possibly conveyed by lifestyle and environmental factors) influences the increase of mutation rate in the mtDNA and thereby contributes to a range of pathophysiological alterations observed in complex diseases. We performed an inverted mitochondrial genome-wide association analysis between mitochondrial nucleotide variants (mtSNVs) and concentration of metabolites. We used 151 metabolites and the whole sequenced mitochondrial genome from 2718 individuals to identify the genetic variants associated with metabolite profiles. Because of the high coverage, next-generation sequencing-based analysis of the mitochondrial genome allows for an accurate detection of mitochondrial heteroplasmy and for the identification of variants associated with the metabolome. The strongest association was found for mt715G > A located in the MT-12SrRNA with the metabolite ratio of C2/C10:1 (P-value = 6.82*10-09, β = 0.909). The second most significant mtSNV was found for mt3714A > G located in the MT-ND1 with the metabolite ratio of phosphatidylcholine (PC) ae C42:5/PC ae C44:5 (P-value = 1.02*10-08, β = 3.631). A large number of significant metabolite ratios were observed involving PC aa C36:6 and the variant mt10689G > A, located in the MT-ND4L gene. These results show an important interconnection between mitochondria and metabolite concentrations. Considering that some of the significant metabolites found in this study have been previously related to complex diseases, such as neurological disorders and metabolic conditions, these associations found here might play a crucial role for further investigations of such complex diseases. Understanding the mechanisms that control human health and disease, in particular, the role of genetic predispositions and their interaction with environmental factors is a prerequisite for the development of safe and efficient therapies for complex disorders.

Axonal degeneration is observed in early stages of several neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). This degeneration generally precedes apoptosis and therefore may be a promising therapeutic target. An increasing number of genes have been identified to actively regulate axonal degeneration and regeneration; however, only a few potential therapeutic targets have been identified in the context of neurodegenerative diseases. Here we investigate DLK-1, a major axonal regeneration pathway and its contribution to axonal degeneration phenotypes in several Caenorhabditis elegans ALS models. From this pathway, we identified the poly (ADP-ribose) (PAR) polymerases (PARP) PARP-1 and PARP-2 as the most consistent modifiers of axonal degeneration in our models of ALS. Genetic and pharmacological inhibition of PARP-1 and PARP-2 reduces axonal degeneration and improves related motor phenotypes.

Recent studies identifying expression quantitative trait loci (eQTLs) in immune cells have uncovered important links between disease risk alleles and gene expression trends in monocytes, T cells and other cell types. However, these studies are generally done with young, healthy subjects, limiting the utility of their findings for age-related conditions such as Alzheimer's disease (AD). We have performed RNA sequencing on four T-cell subsets in genome-wide genotyped and well-characterized AD subjects and age- and sex-matched controls from the Religious Orders Study/Memory and Aging Project. We correlated gene expression data with AD neuropathological traits and with single-nucleotide polymorphisms to detect eQTLs. We identified several significant genes involved in T-cell senescence and cytotoxicity, consistent with T-cell RNA sequencing studies in aged/AD cohorts. We identified unexpected eQTLs previously associated with neuropsychiatric disease traits. Finally, we discovered that pathways related to axon guidance and synaptic function were enriched among trans-eQTLs in coding regions of the genome. Our data strengthen the potential link between T-cell senescence and age-related neurodegenerative disease. In addition, our eQTL data suggest that T-cell phenotypes may influence neuropsychiatric disorders and can be influenced by genes involved in neurodevelopmental processes.

Transcriptome-wide association studies increase the yield of loci associated with disease phenotypes by focusing on expression quantitative trait loci (eQTL). The major source of eQTL data for is the Gene and Tissue Expression (GTEx) project, which is comprised entirely of adults, mainly those >50 years of age at death. Since gene expression levels differ by developmental stage, it is not clear whether eQTLs derived from adult data sources are best suited for use in young-onset diseases such as pediatric cancers. To fill in this knowledge gap, we performed a large-scale eQTL mapping analysis in the GenCord study with newborn samples and compared it with GTEx. Under matched conditions, we found around 80% of the eQTLs in one study can be replicated in the other. However, among all eQTLs identified in GenCord (GTEx), 584 (1045) showed statistically significant differences in effect sizes in GTEx (GenCord). We further investigated how using fetal eQTL data can facilitate the genetic association study of acute lymphoblastic leukemia. GenCord and GTEx identified the same genetic loci with statistical significance; however, the overall association pattern was only weakly correlated. Our paper demonstrates age-differential eQTLs and shows their potential influence on childhood leukemia research.

Intellectual disability (ID) is a neurodevelopmental disorder frequently caused by monogenic defects. In this study, we collected 14 SEMA6B heterozygous variants in 16 unrelated patients referred for ID to different centers. Whereas, until now, SEMA6B variants have mainly been reported in patients with progressive myoclonic epilepsy, our study indicates that the clinical spectrum is wider and also includes non-syndromic ID without epilepsy or myoclonus. To assess the pathogenicity of these variants, selected mutated forms of Sema6b were overexpressed in Human Embryonic Kidney 293T (HEK293T) cells and in primary neuronal cultures. shRNAs targeting Sema6b were also used in neuronal cultures to measure the impact of the decreased Sema6b expression on morphogenesis and synaptogenesis. The overexpression of some variants leads to a subcellular mislocalization of SEMA6B protein in HEK293T cells and to a reduced spine density owing to loss of mature spines in neuronal cultures. Sema6b knockdown also impairs spine density and spine maturation. In addition, we conducted in vivo rescue experiments in chicken embryos with the selected mutated forms of Sema6b expressed in commissural neurons after knockdown of endogenous SEMA6B. We observed that expression of these variants in commissural neurons fails to rescue the normal axon pathway. In conclusion, identification of SEMA6B variants in patients presenting with an overlapping phenotype with ID and functional studies highlight the important role of SEMA6B in neuronal development, notably in spine formation and maturation and in axon guidance. This study adds SEMA6B to the list of ID-related genes.

Mutations in the γ-aminobutyric acid type A (GABAA) receptor γ2 subunit gene, GABRG2, have been associated with a variety of epilepsy syndromes. A de novo mutation (c.T1027C, p.F343L) in GABRG2 was identified in a patient with early onset epileptic encephalopathy. Zebrafish overexpressing mutant human GABRG2 (F343L) subunits displayed spontaneous seizure activity and convulsive behaviors. In this study, we demonstrated that Tg (hGABRG2F343L) zebrafish displayed hyperactivity during light phase with normal circadian rhythm, as well as increased drug-induced locomotor activity. Real-time quantitative PCR, whole mount in situ hybridization and western blotting showed that Tg(hGABRG2F343L) zebrafish had altered expression of GABAA receptor subunits. Furthermore, investigation of synaptic protein expression and synapse ultrastructure uncovered a robust synaptic phenotype that is causally linked to GABRG2(F343L) mutation. Strikingly, Tg(hGABRG2F343L) zebrafish not only had postsynaptic defects, but also displayed an unanticipated deficit at the presynaptic level. Overall, our Tg(hGABRG2F343L) overexpression zebrafish model has expanded the GABAergic paradigm in epileptic encephalopathy from channelopathy to synaptopathy.

High myopia [refractive error ≤ -6 diopters (D)] is a heterogeneous condition, and without clear accompanying features, it can be difficult to pinpoint a genetic cause. This observational study aimed to evaluate the utility of whole exome sequencing (WES) using an eye disorder gene panel in European patients with high myopia. Patients with high myopia were recruited by ophthalmologists and clinical geneticists. Clinical features were categorized into isolated high myopia, high myopia with other ocular involvement or with systemic involvement. WES was performed and an eye disorder gene panel of ~500 genes was evaluated. Hundred and thirteen patients with high myopia [mean (SD) refractive error - 11.8D (5.2)] were included. Of these, 53% were children younger than 12 years of age (53%), 13.3% were aged 12-18 years and 34% were adults (aged > 18 years). Twenty-three out of 113 patients (20%) received a genetic diagnosis of which 11 patients displayed additional ocular or systemic involvement. Pathogenic variants were identified in retinal dystrophy genes (e.g. GUCY2D and CACNA1F), connective tissue disease genes (e.g. COL18A1 and COL2A1), non-syndromic high myopia genes (ARR3), ocular development genes (e.g. PAX6) and other genes (ASPH and CNNM4). In 20% of our high myopic study population, WES using an eye gene panel enabled us to diagnose the genetic cause for this disorder. Eye genes known to cause retinal dystrophy, developmental or syndromic disorders can cause high myopia without apparent clinical features of other pathology.

Mutations in Talpid3, a basal body protein essential for the assembly of primary cilia, have been reported to be causative for Joubert Syndrome (JS). Herein, we report prominent developmental defects in the hippocampus of a conditional knockout mouse lacking the conserved exons 11 and 12 of Talpid3. At early postnatal stages, the Talpid3 mutants exhibit a reduction in proliferation in the dentate gyrus and a disrupted glial scaffold. The occurrence of mis-localized progenitors in the granule cell layer suggests a role for the disrupted glial scaffold in cell migration resulting in defective subpial neurogenic zone-to-hilar transition. Neurospheres derived from the hippocampus of Talpid3fl/flUbcCre mouse, in which Talpid3 was conditionally deleted, lacked primary cilia and were smaller in size. In addition, neurosphere cells showed a disrupted actin cytoskeleton and defective migration. Our findings suggest a link between the hippocampal defects and the learning/memory deficits seen in JS patients.

A disproportionate tall stature is the most evident manifestation in Marfan syndrome (MFS), a multisystem condition caused by mutations in the extracellular protein and TGFβ modulator, fibrillin-1. Unlike cardiovascular manifestations, there has been little effort devoted to unravel the molecular mechanism responsible for long bone overgrowth in MFS. By combining the Cre-LoxP recombination system with metatarsal bone cultures, here we identify the outer layer of the perichondrium as the tissue responsible for long bone overgrowth in MFS mice. Analyses of differentially expressed genes in the fibrillin-1-deficient perichondrium predicted that loss of TGFβ signaling may influence chondrogenesis in the neighboring epiphyseal growth plate (GP). Immunohistochemistry revealed that fibrillin-1 deficiency in the outer perichondrium is associated with decreased accumulation of latent TGFβ-binding proteins (LTBPs)-3 and -4, and reduced levels of phosphorylated (activated) Smad2. Consistent with these findings, mutant metatarsal bones grown in vitro were longer and released less TGFβ than the wild-type counterparts. Moreover, addition of recombinant TGFβ1 normalized linear growth of mutant metatarsal bones. We conclude that longitudinal bone overgrowth in MFS is accounted for by diminished sequestration of LTBP-3 and LTBP-4 into the fibrillin-1-deficient matrix of the outer perichondrium, which results in less TGFβ signaling locally and improper GP differentiation distally.

Most of the early Hungarian tribes originated from the Volga-Kama and South-Ural regions, where they were composed of a mixed population based on historical, philological and archaeological data. We present here the uniparental genetic makeup of the mediaeval era of these regions that served as a melting pot for ethnic groups with different linguistic and historical backgrounds. Representing diverse cultural contexts, the new genetic data originate from ancient proto-Ob-Ugric people from Western Siberia (6th-13th century), the pre-Conquest period and subsisting Hungarians from the Volga-Ural region (6th-14th century) and their neighbours. By examining the eastern archaeology traits of Hungarian prehistory, we also study their genetic composition and origin in an interdisciplinary framework. We analyzed 110 deep-sequenced mitogenomes and 42 Y-chromosome haplotypes from 18 archaeological sites in Russia. The results support the studied groups' genetic relationships regardless of geographical distances, suggesting large-scale mobility. We detected long-lasting genetic connections between the sites representing the Kushnarenkovo and Chiyalik cultures and the Carpathian Basin Hungarians and confirmed the Uralic transmission of several East Eurasian uniparental lineages in their gene pool. Based on phylogenetics, we demonstrate and model the connections and splits of the studied Volga-Ural and conqueror groups. Early Hungarians and their alliances conquered the Carpathian Basin around 890 AD. Re-analysis of the Hungarian conquerors' maternal gene pool reveals numerous surviving maternal relationships in both sexes; therefore, we conclude that men and women came to the Carpathian Basin together, and although they were subsequently genetically fused into the local population, certain eastern lineages survived for centuries.

In human autosomal recessive woolly hair/hypotrichosis (ARWH/HT), many mutations have been identified in a gene encoding LPA6, a G protein-coupled receptor (GPCR) for lysophosphatidic acid (LPA). However, information regarding the effects of such mutations on receptor function is limited. In this study, we examined functional impacts of selected amino acid changes in LPA6 identified in ARWH/HT patients. In our exogenous expression experiments, all mutants except S3T failed to respond to LPA, indicating that they are loss-of-function mutants. Among the nine mutants, five (D63V, G146R, N246D, L277P, and C278Y) displayed impaired expression at the cell surface due to endoplasmic reticulum (ER) retention, indicating that these mutants are trafficking-defective, as reported in other disease-associated GPCRs. Notably, alkyl-OMPT, a potent synthetic agonist for LPA6 restored the defective cell surface expression of two of the ER-retained mutants, D63V and N246D, possibly by its chaperoning function that allows them to escape intracellular retention as well as proteasomal degradation. Furthermore, the alkyl-OMPT-rescued N246D mutant was shown be functional. Our findings encourage future application of pharmacoperone therapy for ARWH/HT patients with specific LPA6 mutations.

BACKGROUND: The endoplasmic reticulum (ER)-membrane protein complex (EMC) is a multi-protein transmembrane complex composed of 10 subunits that functions as a membrane-protein chaperone. Variants in EMC1 lead to neurodevelopmental delay and cerebellar degeneration. Multiple families with biallelic variants have been published, yet to date, only a single report of a monoallelic variant has been described, and functional evidence is sparse. METHODS: Exome sequencing was used to investigate the genetic cause underlying severe developmental delay in three unrelated children. EMC1 variants were modeled in Drosophila, using loss-of-function (LoF) and overexpression studies. Glial-specific and neuronal-specific assays were used to determine whether the dysfunction was specific to one cell type. RESULTS: Exome sequencing identified de novo variants in EMC1 in three individuals affected by global developmental delay, hypotonia, seizures, visual impairment and cerebellar atrophy. All variants were located at Pro582 or Pro584. Drosophila studies indicated that imbalance of EMC1-either overexpression or knockdown-results in pupal lethality and suggest that the tested homologous variants are LoF alleles. In addition, glia-specific gene dosage, overexpression or knockdown, of EMC1 led to lethality, whereas neuron-specific alterations were tolerated. DISCUSSION: We establish de novo monoallelic EMC1 variants as causative of a neurological disease trait by providing functional evidence in a Drosophila model. The identified variants failed to rescue the lethality of a null allele. Variations in dosage of the wild-type EMC1, specifically in glia, lead to pupal lethality, which we hypothesize results from the altered stoichiometry of the multi-subunit protein complex EMC.

Genome-wide association studies (GWAS) have identified more than 75 genetic variants associated with Alzheimer's disease (ad). However, how these variants function and impact protein expression in brain regions remain elusive. Large-scale proteomic datasets of ad postmortem brain tissues have become available recently. In this study, we used these datasets to investigate brain region-specific molecular pathways underlying ad pathogenesis and explore their potential drug targets. We applied our new network-based tool, Edge-Weighted Dense Module Search of GWAS (EW_dmGWAS), to integrate ad GWAS statistics of 472 868 individuals with proteomic profiles from two brain regions from two large-scale ad cohorts [parahippocampal gyrus (PHG), sample size n = 190; dorsolateral prefrontal cortex (DLPFC), n = 192]. The resulting network modules were evaluated using a scale-free network index, followed by a cross-region consistency evaluation. Our EW_dmGWAS analyses prioritized 52 top module genes (TMGs) specific in PHG and 58 TMGs in DLPFC, of which four genes (CLU, PICALM, PRRC2A and NDUFS3) overlapped. Those four genes were significantly associated with ad (GWAS gene-level false discovery rate < 0.05). To explore the impact of these genetic components on TMGs, we further examined their differentially co-expressed genes at the proteomic level and compared them with investigational drug targets. We pinpointed three potential drug target genes, APP, SNCA and VCAM1, specifically in PHG. Gene set enrichment analyses of TMGs in PHG and DLPFC revealed region-specific biological processes, tissue-cell type signatures and enriched drug signatures, suggesting potential region-specific drug repurposing targets for ad.

Leber's hereditary optic neuropathy (LHON) is a maternally inherited eye disease due to mitochondrial DNA (mtDNA) mutations. LHON-linked ND6 14484T > C (p.M64V) mutation affected structural components of complex I but its pathophysiology is poorly understood. The structural analysis of complex I revealed that the M64 forms a nonpolar interaction Y59 in the ND6, Y59 in the ND6 interacts with E34 of ND4L, and L60 of ND6 interacts with the Y114 of ND1. These suggested that the m.14484T > C mutation may perturb the structure and function of complex I. Mutant cybrids constructed by transferring mitochondria from lymphoblastoid cell lines of one Chinese LHON family into mtDNA-less (ρo) cells revealed decreases in the levels of ND6, ND1 and ND4L. The m.14484T > C mutation may affect mitochondrial mRNA homeostasis, supported by reduced levels of SLIRP and SUPV3L1 involved in mRNA degradation and increasing expression of ND6, ND1 and ND4L genes. These alterations yielded decreased activity of complex I, respiratory deficiency, diminished mitochondrial ATP production and reduced membrane potential, and increased production of reactive oxygen species in the mutant cybrids. Furthermore, the m.14484T > C mutation promoted apoptosis, evidenced by elevating Annexin V-positive cells, release of cytochrome c into cytosol, levels in apoptotic proteins BAX, caspases 3, 7, 9 and decreasing levels in anti-apoptotic protein Bcl-xL in the mutant cybrids. Moreover, the cybrids bearing the m.14484T > C mutation exhibited the reduced levels of autophagy protein LC3, increased levels of substrate P62 and impaired PINK1/Parkin-dependent mitophagy. Our findings highlighted the critical role of m.14484T > C mutation in the pathogenesis of LHON.

Polygenic scores (PGS) can identify individuals at risk of adverse health events and guide genetics-based personalized medicine. However, it is not clear how well PGS translate between different populations, limiting their application to well-studied ethnicities. Proteins are intermediate traits linking genetic predisposition and environmental factors to disease, with numerous blood circulating protein levels representing functional readouts of disease-related processes. We hypothesized that studying the genetic architecture of a comprehensive set of blood-circulating proteins between a European and an Arab population could shed fresh light on the translatability of PGS to understudied populations. We therefore conducted a genome-wide association study with whole-genome sequencing data using 1301 proteins measured on the SOMAscan aptamer-based affinity proteomics platform in 2935 samples of Qatar Biobank and evaluated the replication of protein quantitative traits (pQTLs) from European studies in an Arab population. Then, we investigated the colocalization of shared pQTL signals between the two populations. Finally, we compared the performance of protein PGS derived from a Caucasian population in a European and an Arab cohort. We found that the majority of shared pQTL signals (81.8%) colocalized between both populations. About one-third of the genetic protein heritability was explained by protein PGS derived from a European cohort, with protein PGS performing ~ 20% better in Europeans when compared to Arabs. Our results are relevant for the translation of PGS to non-Caucasian populations, as well as for future efforts to extend genetic research to understudied populations.

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes responsible for charging amino acids onto cognate tRNAs during protein synthesis. In histidyl-tRNA synthetase (HARS), autosomal dominant mutations V133F, V155G, Y330C, S356N in the HARS catalytic domain cause Charcot Marie Tooth Disease Type 2 W (CMT2W), while RNA-binding domain mutation Y454S causes recessive Usher Syndrome Type IIIB (USH3B). In a yeast model, all human HARS variants complemented a genomic deletion of the yeast ortholog hts1 at high expression levels. CMT2W associated mutations, but not Y454S, resulted in reduced growth. We show mistranslation of histidine to glutamine and threonine in V155G and S356N but not Y330C mutants in yeast. Mistranslating V155G and S356N mutants lead to accumulation of insoluble proteins, which was rescued by histidine. Mutants V133F and Y330C showed the most significant growth defect and decreased HARS abundance in cells. Here, histidine supplementation led to insoluble protein aggregation and further reduced viability, indicating histidine toxicity associated with these mutants. V133F proteins displayed reduced thermal stability in vitro, which was rescued by tRNA. Our data will inform future treatment options for HARS patients, where histidine supplementation may either have a toxic or compensating effect depending on the nature of the causative HARS variant.

Rod and cone photoreceptors in the retina mediate dim light and day light vision, respectively. Despite their distinctive functions, rod and cone visual pigments utilize the same vitamin A-derived chromophore. To sustain vision, vitamin A precursors must be acquired in the gut, metabolized, and distributed to the eyes. Deficiencies in this pathway in inherited ocular disease states deplete cone photoreceptors from chromophore and eventually lead to cell death, whereas the more abundant rod photoreceptors are less affected. However, pathways that support cone function and survival under such conditions are largely unknown. Using biochemical, histological, and physiological approaches, we herein show that intervention with β-carotene in STRA6-deficient mice improved chromophore supply to cone photoreceptors. Relieving the inherent negative feedback regulation of β-carotene-oxygenase-1 activity in the intestine by genetic means further bolstered cone photoreceptor functioning in the STRA6-deficient eyes. A vitamin A-rich diet, however, did not improve cone photoreceptor function in STRA6-deficiency. We provide evidence that the beneficial effect of β-carotene on cones results from favorable serum kinetics of retinyl esters in lipoproteins. The respective alterations in lipoprotein metabolism maintained a steady supply of retinoids to the STRA6-deficient eyes that ameliorated the competition for chromophore between rod and cone photoreceptors. Together, our study elucidates a cone photoreceptor-survival pathway and unravels an unexpected metabolic connection between the gut and the retina.

Few GWAS studies analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of twelve morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity, and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (pvalue = 8.33x10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

Sequence variants or mutations in the GBA gene are numerically the most important risk factor for Parkinson disease (PD). The GBA gene encodes for the lysosomal hydrolase enzyme, glucocerebrosidase (GCase). GBA mutations often reduce GCase activity and lead to impairment of the autophagy-lysosomal pathway, which is important in the turnover of alpha-synuclein, accumulation of which is a key pathological hallmark of PD. Although the E326K variant is one of the most common GBA variants associated with PD, there is limited understanding of its biochemical effects. We have characterised homozygous and heterozygous E326K variants in human fibroblasts. We found that E326K variants did not cause significant loss of GCase protein or activity, endoplasmic reticulum (ER) retention or ER stress, in contrast to the L444P GBA mutation. This was confirmed in human dopaminergic SH-SY5Y neuroblastoma cell lines over-expressing GCase with either E326K or L444P protein. Despite no loss of GCase activity, a significant increase of insoluble alpha-synuclein aggregates in E326K and L444P mutants was observed. Notably, SH-SY5Y over-expressing E326K demonstrated a significant increase in lipid droplet number under basal conditions, which was exacerbated following treatment with the fatty acid oleic acid. Similarly, a significant increase in lipid droplet formation following lipid loading was observed in heterozygous and homozygous E326K fibroblasts. In conclusion, the work presented here demonstrates that the E326K mutation behaves differently to common loss of function GBA mutations, however lipid dyshomeostasis and alpha-synuclein pathology is still evident.

Hereditary hearing loss is a highly genetically heterogeneous disorder. More than 150 genes have been identified to link to human nonsyndromic hearing impairment. To identify genetic mutations and underlying molecular mechanisms in affected individuals and families with congenital hearing loss, we recruited a cohort of 389 affected individuals in 354 families for whole-exome sequencing analysis. In the present study, we report a novel homozygous missense variant (c.233A > G, p.Lys78Arg) in the OXR1 gene, which was identified in a 4-year-old girl with sensorineural hearing loss (SNHL). OXR1 encodes Oxidation Resistance 1 and is evolutionarily conserved from zebrafish to human. We found that the ortholog oxr1b gene is expressed in the statoacoustic ganglion (SAG, a sensory ganglion of ear) and posterior lateral line ganglion (pLL) in zebrafish. Knockdown of oxr1b in zebrafish resulted in a significant developmental defect of SAG and pLL. This phenotype can be rescued by co-injection of wild-type human OXR1 mRNAs, but not mutant OXR1 (c.233A > G) mRNAs. OXR1-associated pathway analysis revealed that mutations of TBC1D24, a TLDc-domain-containing homolog gene of OXR1, have previously been identified in patients with hearing loss. Interestingly, mutations or knockout of OXR1 interacting molecules such as ATP6V1B1 and ESR1 are also associated with hearing loss in patients or animal models, hinting an important role of OXR1 and associated partners in cochlear development and hearing function.

Author information: (1)Department of Medical Genetics, Hospital Universitario de Navarra, IdiSNA, Pamplona 31008, Spain. (2)Department of Health Sciences, Universidad Pública de Navarra, IdiSNA, Pamplona 31008, Spain. (3)Fundación Miguel Servet-Navarrabiomed, IdiSNA, Pamplona 31008, Spain. (4)Department of Statistics, Informatics and Mathematics, Universidad Pública de Navarra, IdiSNA, Pamplona 31006, Spain. (5)School of Molecular Biosciences, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow G12 8QQ, UK.

Lamin B Receptor (LBR) is an inner nuclear membrane protein that assembles the nuclear envelope post mitosis. Here we show that LBR depletion induces mitotic defects accompanied by recurrent chromosomal losses. In addition, LBR knockdown results in nuclear aberrations such as nuclear blebs and micronuclei, with chromosomes showing a higher frequency of losses, being enriched within the micronucleus. Furthermore, doxycycline-induced conditional depletion of LBR significantly increased tumor volumes that form within the subcutaneous xenografts of mice. Of note, the tumor-derived primary cells recapitulated chromosomal losses and gains, revealing a novel role for LBR as a tumor suppressor. Co-immunoprecipitation (Co-IP) of LBR uncovered an association of LBR with telomere-associated factors. Interestingly, qPCR array-based gene expression profiling showed a significant upregulation of TRF1 upon LBR depletion. Remarkably, TRF1 knockdown in the background of LBR depletion maintains chromosomal stability, unraveling a novel mechanism involving LBR and TRF in the maintenance of chromosomal stability in colorectal cancer cells.

BACKGROUND: To summarize the clinical and genetic characteristics of patients with pyridox(am)ine-5'-phosphate oxidase (PNPO) deficiency. METHODS: Clinical and genetic data of the patients were collected and analyzed. RESULTS: 18 patients from 17 families with variants in PNPO were collected, and 15 cases survived to date. The age of onset ranged from 1 day to 5 months (median age 6.5 days) and seven of them presented with seizures less than 24 hours. 7/18 (39%) of patients showed seizure free with pyridoxine (PN) or pyridoxal-5'-phosphate (PLP) treatment. Two patients showed surprised therapeutic responses to antiseizure medications (ASMs) therapy: one could be controlled for up to 1 year and 5 months, and the other showed seizure free for more than eight years. The neurodevelopment was normal in one patient, mild delay in four, in whom responded well to PN. Severe delay could be seen in the remaining 10 surviving patients. Genetic analysis revealed 14 variants of PNPO, seven of which were novel. Five pairs of unrelated patients were observed to carry the same variants respectively, and had similar developmental status and onset age of seizures in some degree in each pair, while also had differences. CONCLUSIONS: The clinical characteristics, including age of onset, treatment response, and prognosis, were variable and difficult to classify into different types clearly. Patients with PNPO deficiency who used PN as their main treatment and being able to control seizures seemed to be associated with better outcomes. Patients with the same genotype tended to show the correlation of phenotype-genotype.

Asians are underrepresented across many omics databases, thereby limiting the potential of precision medicine in nearly 60% of the global population. As such, there is a pressing need for multi-omics derived quantitative trait loci (QTLs) to fill the knowledge gap of complex traits in populations of Asian ancestry. Here, we provide the first blood-based multi-omics analysis of Asian pregnant women, constituting high-resolution genotyping (N = 1079), DNA methylation (N = 915) and transcriptome profiling (N = 238). Integrative omics analysis identified 219 154 CpGs associated with cis-DNA methylation QTLs (meQTLs) and 3703 RNAs associated with cis-RNA expression QTLs (eQTLs). Ethnicity was the largest contributor of inter-individual variation across all omics datasets, with 2561 genes identified as hotspots of this variation; 395 of these hotspot genes also contained both ethnicity-specific eQTLs and meQTLs. Gene set enrichment analysis of these ethnicity QTL hotspots showed pathways involved in lipid metabolism, adaptive immune system and carbohydrate metabolism. Pathway validation by profiling the lipidome (~480 lipids) of antenatal plasma (N = 752) and placenta (N = 1042) in the same cohort showed significant lipid differences among Chinese, Malay and Indian women, validating ethnicity-QTL gene effects across different tissue types. To develop deeper insights into the complex traits and benefit future precision medicine research in Asian pregnant women, we developed iMOMdb, an open-access database.

BACKGROUND: TASP1 encodes an endopeptidase activating histone methyltransferases of the KMT2 family. Homozygous loss-of-function variants in TASP1 have recently been associated with Suleiman-El-Hattab syndrome. We report six individuals with Suleiman-El-Hattab syndrome and provide functional characterization of this novel histone modification disorder in a multi-omics approach. METHODS: Chromosomal microarray/exome sequencing in all individuals. Western blotting from fibroblasts in two individuals. RNA sequencing and proteomics from fibroblasts in one individual. Methylome analysis from blood in two individuals. Knock-out of tasp1 orthologue in zebrafish and phenotyping. RESULTS: All individuals had biallelic TASP1 loss-of-function variants and a phenotype including developmental delay, multiple congenital anomalies (including cardiovascular and posterior fossa malformations), a distinct facial appearance and happy demeanor. Western blot revealed absence of TASP1. RNA sequencing/proteomics showed HOX gene downregulation (HOXA4, HOXA7, HOXA1 and HOXB2) and dysregulation of transcription factor TFIIA. A distinct methylation profile intermediate between control and Kabuki syndrome (KMT2D) profiles could be produced. Zebrafish tasp1 knock-out revealed smaller head size and abnormal cranial cartilage formation in tasp1 crispants. CONCLUSION: This work further delineates Suleiman-El-Hattab syndrome, a recognizable neurodevelopmental syndrome. Possible downstream mechanisms of TASP1 deficiency include perturbed HOX gene expression and dysregulated TFIIA complex. Methylation pattern suggests that Suleiman-El-Hattab syndrome can be categorized into the group of histone modification disorders including Wiedemann-Steiner and Kabuki syndrome.

Large-scale genomic studies of schizophrenia implicate genes involved in the epigenetic regulation of transcription by histone methylation and genes encoding components of the synapse. However, the interactions between these pathways in conferring risk to psychiatric illness are unknown. Loss-of-function (LoF) mutations in the gene encoding histone methyltransferase, SETD1A, confer substantial risk to schizophrenia. Among several roles, SETD1A is thought to be involved in the development and function of neuronal circuits. Here, we employed a multi-omics approach to study the effects of heterozygous Setd1a LoF on gene expression and synaptic composition in mouse cortex across five developmental timepoints from embryonic day 14 to postnatal day 70. Using RNA sequencing, we observed that Setd1a LoF resulted in the consistent downregulation of genes enriched for mitochondrial pathways. This effect extended to the synaptosome, in which we found age-specific disruption to both mitochondrial and synaptic proteins. Using large-scale patient genomics data, we observed no enrichment for genetic association with schizophrenia within differentially expressed transcripts or proteins, suggesting they derive from a distinct mechanism of risk from that implicated by genomic studies. This study highlights biological pathways through which SETD1A LOF may confer risk to schizophrenia. Further work is required to determine whether the effects observed in this model reflect human pathology.

Spinal muscular atrophy (SMA) is a childhood motor neuron disease caused by anomalies in the SMN1 gene. Although therapeutics have been approved for the treatment of SMA, there is a therapeutic time window, after which efficacy is reduced. Hallmarks of motor unit pathology in SMA include loss of motor-neurons and neuromuscular junction (NMJs). Following an increase in Smn levels, it is unclear how much damage can be repaired and the degree to which normal connections are re-established. Here, we perform a detailed analysis of motor unit pathology before and after restoration of Smn levels. Using a Smn-inducible mouse model of SMA, we show that genetic restoration of Smn results in a dramatic reduction in NMJ pathology, with restoration of innervation patterns, preservation of axon and endplate number and normalized expression of P53-associated transcripts. Notably, presynaptic swelling and elevated Pmaip levels remained. We analysed the effect of either early or delayed treated of an antisense oligonucleotide (ASO) targeting SMN2 on a range of differentially vulnerable muscles. Following ASO administration, the majority of endplates appeared fully occupied. However, there was an underlying loss of axons and endplates, which was more prevalent following a delay in treatment. There was an increase in average motor unit size following both early and delayed treatment. Together this work demonstrates the remarkably regenerative capacity of the motor neuron following Smn restoration, but highlights that recovery is incomplete. This work suggests that there is an opportunity to enhance neuromuscular junction recovery following administration of Smn-enhancing therapeutics.

The coronavirus disease 2019 (COVID-19) pandemic has changed the paradigms for disease surveillance and rapid deployment of scientific-based evidence for understanding disease biology, susceptibility and treatment. We have organized a large-scale genome-wide association study (GWAS) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infected individuals in Sao Paulo, Brazil, one of the most affected areas of the pandemic in the country, itself one of the most affected in the world. Here, we present the results of the initial analysis in the first 5233 participants of the BRACOVID study. We have conducted a GWAS for COVID-19 hospitalization enrolling 3533 cases (hospitalized COVID-19 participants) and 1700 controls (non-hospitalized COVID-19 participants). Models were adjusted by age, sex and the 4 first principal components. A meta-analysis was also conducted merging BRACOVID hospitalization data with the Human Genetic Initiative (HGI) Consortia results. BRACOVID results validated most loci previously identified in the HGI meta-analysis. In addition, no significant heterogeneity according to ancestral group within the Brazilian population was observed for the two most important COVID-19 severity associated loci: 3p21.31 and Chr21 near IFNAR2. Using only data provided by BRACOVID, a new genome-wide significant locus was identified on Chr1 near the genes DSTYK and RBBP5. The associated haplotype has also been previously associated with a number of blood cell related traits and might play a role in modulating the immune response in COVID-19 cases.

Most epigenetic epidemiology to date has utilized microarrays to identify positions in the genome where variation in DNA methylation is associated with environmental exposures or disease. However, these profile less than 3% of DNA methylation sites in the human genome, potentially missing affected loci and preventing the discovery of disrupted biological pathways. Third generation sequencing technologies, including Nanopore sequencing, have the potential to revolutionize the generation of epigenetic data, not only by providing genuine genome-wide coverage but profiling epigenetic modifications direct from native DNA. Here we assess the viability of using Nanopore sequencing for epidemiology by performing a comparison with DNA methylation quantified using the most comprehensive microarray available, the Illumina EPIC array. We implemented a CRISPR-Cas9 targeted sequencing approach in concert with Nanopore sequencing to profile DNA methylation in three genomic regions to attempt to rediscover genomic positions that existing technologies have shown are differentially methylated in tobacco smokers. Using Nanopore sequencing reads, DNA methylation was quantified at 1779 CpGs across three regions, providing a finer resolution of DNA methylation patterns compared to the EPIC array. The correlation of estimated levels of DNA methylation between platforms was high. Furthermore, we identified 12 CpGs where hypomethylation was significantly associated with smoking status, including 10 within the AHRR gene. In summary, Nanopore sequencing is a valid option for identifying genomic loci where large differences in DNAm are associated with a phenotype and has the potential to advance our understanding of the role differential methylation plays in the etiology of complex disease.

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.

Myotonic dystrophy (DM) is caused by expansions of C(C)TG repeats in the non-coding regions of the DMPK and CNBP genes, and DM patients often suffer from sudden cardiac death due to lethal conduction block or arrhythmia. Specific molecular changes that underlie DM cardiac pathology have been linked to repeat-associated depletion of Muscleblind-like (MBNL) 1 and 2 proteins and upregulation of CUGBP, Elav-like family member 1 (CELF1). Hypothesis solely targeting MBNL1 or CELF1 pathways that could address all the consequences of repeat expansion in heart remained inconclusive, particularly when the direct cause of mortality and results of transcriptome analyses remained undetermined in Mbnl compound knockout (KO) mice with cardiac phenotypes. Here, we develop Myh6-Cre double KO (DKO) (Mbnl1-/-; Mbnl2cond/cond; Myh6-Cre+/-) mice to eliminate Mbnl1/2 in cardiomyocytes and observe spontaneous lethal cardiac events under no anesthesia. RNA sequencing recapitulates DM heart spliceopathy and shows gene expression changes that were previously undescribed in DM heart studies. Notably, immunoblotting reveals a nearly 6-fold increase of Calsequestrin 1 and 50% reduction of epidermal growth factor proteins. Our findings demonstrate that complete ablation of MBNL1/2 in cardiomyocytes is essential for generating sudden death due to lethal cardiac rhythms and reveal potential mechanisms for DM heart pathogenesis.

Type 2 diabetes is a complex, systemic disease affected by both genetic and environmental factors. Previous research has identified genetic variants associated with type 2 diabetes risk; however, gene regulatory changes underlying progression to metabolic dysfunction are still largely unknown. We investigated RNA expression changes that occur during diabetes progression using a two-stage approach. In our discovery stage, we compared changes in gene expression using two longitudinally collected blood samples from subjects whose fasting blood glucose transitioned to a level consistent with type 2 diabetes diagnosis between the time points against those who did not with a novel analytical network approach. Our network methodology identified 17 networks, one of which was significantly associated with transition status. This 822-gene network harbors many genes novel to the type 2 diabetes literature but is also significantly enriched for genes previously associated with type 2 diabetes. In the validation stage, we queried associations of genetically determined expression with diabetes-related traits in a large biobank with linked electronic health records. We observed a significant enrichment of genes in our identified network whose genetically determined expression is associated with type 2 diabetes and other metabolic traits and validated 31 genes that are not near previously reported type 2 diabetes loci. Finally, we provide additional functional support, which suggests that the genes in this network are regulated by enhancers that operate in human pancreatic islet cells. We present an innovative and systematic approach that identified and validated key gene expression changes associated with type 2 diabetes transition status and demonstrated their translational relevance in a large clinical resource.

Polygenic risk scores (PRSs) are useful for predicting breast cancer risk, but the prediction accuracy of existing PRSs in women of African ancestry (AA) remains relatively low. We aim to develop optimal PRSs for the prediction of overall and estrogen receptor (ER) subtype-specific breast cancer risk in AA women. The AA dataset comprised 9235 cases and 10 184 controls from four genome-wide association study (GWAS) consortia and a GWAS study in Ghana. We randomly divided samples into training and validation sets. We built PRSs using individual-level AA data by a forward stepwise logistic regression and then developed joint PRSs that combined (1) the PRSs built in the AA training dataset and (2) a 313-variant PRS previously developed in women of European ancestry. PRSs were evaluated in the AA validation set. For overall breast cancer, the odds ratio per standard deviation of the joint PRS in the validation set was 1.34 [95% confidence interval (CI): 1.27-1.42] with the area under receiver operating characteristic curve (AUC) of 0.581. Compared with women with average risk (40th-60th PRS percentile), women in the top decile of the PRS had a 1.98-fold increased risk (95% CI: 1.63-2.39). For PRSs of ER-positive and ER-negative breast cancer, the AUCs were 0.608 and 0.576, respectively. Compared with existing methods, the proposed joint PRSs can improve prediction of breast cancer risk in AA women.

Sensorineural hearing loss often results from damaged or deficient inner ear hair cells. Mitochondrial 12S rRNA 1555A>G mutation has been associated with hearing loss in many families. The m.1555A>G mutation is a primary factor underlying the development of hearing loss and TRMU allele (c.28G>T, p.Ala10Sser) encoding tRNA thiouridylase interact with m.1555A>G mutation to cause hearing loss. However, the tissue specificity of mitochondrial dysfunction remains elusive and there is no highly effective therapy for mitochondrial deafness. We report here the generation of induced pluripotent stem cells (iPSCs) from lymphoblastoid cell lines derived from members of an Arab-Israeli family (asymptomatic individual carrying only m.1555A>G mutation, symptomatic individual bearing both m.1555A>G and c.28G>T mutations, and control subject). The c.28G>T mutation in iPSC lines from a hearing-impaired subject was corrected by CRISPR/Cas9. These iPSCs were differentiated into otic epithelial progenitor (OEP) cells and subsequent inner ear hair cell (HC)-like cells. The iPSCs bearing m.1555A>G mutation exhibited mildly deficient differentiation into OEP and resultant HC-like cells displayed mild defects in morphology and electrophysiological properties. Strikingly, those HC-like cells harboring m.1555A>G and TRMU c.28G>T mutations displayed greater defects in the development, morphology and functions than those in cells bearing only m.1555A>G mutation. Transcriptome analysis of patients-derived HC-like cells revealed altered expressions of genes vital for mechanotransduction of hair cells. Genetic correction of TRMU c.28G>T mutation yielded morphologic and functional recovery of patient derived HC-like cells. These findings provide new insights into pathophysiology of maternally inherited hearing loss and a step toward therapeutic interventions for this disease.

RTL1/PEG11, which has been associated with anxiety disorders, is a retrotransposon-derived imprinted gene in the placenta. However, imprinting patterns and functions of RTL1 in the brain have not been well-investigated. We found Rtl1 was paternally, but not maternally, expressed in brain stem, thalamus, and hypothalamus of mice, and imprinting status of RTL1 was maintained in human brain. Paternal Rtl1 knockout (Rtl1m+/p-) mice had higher neonatal death rates due to impaired suckling, and low body weights beginning on embryonic day 16.5. High paternal expression of Rtl1 was detected in the locus coeruleus (LC) and Rtl1m+/p- mice showed an increased delay in time of onset for action potentials and inward currents with decreased neuronal excitability of LC neurons. Importantly, Rtl1m+/p- mice exhibited behaviors associated with anxiety, depression, fear-related learning and memory, social dominance, and low locomotor activity. Taken together, our findings demonstrate RTL1 is imprinted in brain, mediates emotional and social behaviors, and regulates excitability in LC neurons.

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Precision medicine carries great potential for management of all tumor types. The aim of this retrospective study was to investigate if the two most common genetically distinct uterine fibroid subclasses, driven by aberrations in MED12 and HMGA2 genes, respectively, influence response to treatment with the progesterone receptor modulator ulipristal acetate. Changes in diameter and mutation status were derived for 101 uterine fibroids surgically removed after ulipristal acetate treatment. A significant difference in treatment response between the two major subclasses was detected. MED12 mutant fibroids had 4.4 times higher odds of shrinking in response to ulipristal acetate treatment as compared to HMGA2 driven fibroids (95% confidence interval 1.37-13.9; P = 0.013), and in a multivariate analysis molecular subclassification was an independent predictive factor. Compatible with this finding, gene expression and DNA methylation analyses revealed subclass specific differences in progesterone receptor signaling. The work provides a proof-of-principle that uterine fibroid treatment response is influenced by molecular subclass and that the genetic subclasses should be taken into account when evaluating current and future uterine fibroid therapies.

Hereditary hearing loss has a genetic and phenotypic heterogeneity. However, it is still difficult to explain this heterogeneity perfectly with known deafness genes. Here, we report a novel causative gene EPHA10 as well as its non-coding variant in 5' untranslated region identified in a family with post-lingual autosomal dominant non-syndromic hearing loss from southern China. One affected member of this family had an ideal hearing restoration after cochlear implantation. We speculated that there were probable deafness-causing abnormalities in the cochlea according to clinical imaging and auditory evaluations. A heterozygous variant c.-81_-73delinsAGC was found co-segregating with hearing loss. Epha10 was expressed in mouse cochlea at both transcription and translation levels. The variant caused up-regulation of EPHA10 which may result from promoter activity enhancement after sequence change. Over-expression of Eph (the homologue of human EPHA10) exerted effects on the structure and function of chordotonal organ in fly model. In summary, our study linked pseudo-kinase EPHA10 to hearing loss in humans for the first time.