A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats.

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Increased soluble urokinase plasminogen activator levels modulate monocyte function to promote atherosclerosis.

People with kidney disease are disproportionately affected by atherosclerosis for unclear reasons. Soluble urokinase plasminogen activator receptor (suPAR) is an immune-derived mediator of kidney disease, levels of which are strongly associated with cardiovascular outcomes. We assessed suPAR's pathogenic involvement in atherosclerosis using epidemiologic, genetic, and experimental approaches. We found serum suPAR levels to be predictive of coronary artery calcification and cardiovascular events in 5,406 participants without known coronary disease. In a genome-wide association meta-analysis including over 25,000 individuals, we identified a missense variant in the plasminogen activator, urokinase receptor (PLAUR) gene (rs4760), confirmed experimentally to lead to higher suPAR levels. Mendelian randomization analysis in the UK Biobank using rs4760 indicated a causal association between genetically predicted suPAR levels and atherosclerotic phenotypes. In an experimental model of atherosclerosis, proprotein convertase subtilisin/kexin-9 (Pcsk9) transfection in mice overexpressing suPAR (suPARTg) led to substantially increased atherosclerotic plaques with necrotic cores and macrophage infiltration compared with those in WT mice, despite similar cholesterol levels. Prior to induction of atherosclerosis, aortas of suPARTg mice excreted higher levels of CCL2 and had higher monocyte counts compared with WT aortas. Aortic and circulating suPARTg monocytes exhibited a proinflammatory profile and enhanced chemotaxis. These findings characterize suPAR as a pathogenic factor for atherosclerosis acting at least partially through modulation of monocyte function.

The phenotypic spectrum associated with OTX2 mutations in humans.

Objective: The transcription factor OTX2 is implicated in ocular, craniofacial, and pituitary development. Design: We aimed to establish the contribution of OTX2 mutations in congenital hypopituitarism patients with/without eye abnormalities, study functional consequences, and establish OTX2 expression in the human brain, with a view to investigate the mechanism of action. Methods: We screened patients from the UK (n = 103), international centres (n = 24), and Brazil (n = 282); 145 were within the septo-optic dysplasia spectrum, and 264 had no eye phenotype. Transactivation ability of OTX2 variants was analysed in murine hypothalamic GT1-7 neurons. In situ hybridization was performed on human embryonic brain sections. Genetically engineered mice were generated with a series of C-terminal OTX2 variants. Results: Two chromosomal deletions and six haploinsufficient mutations were identified in individuals with eye abnormalities; an affected relative of one patient harboured the same mutation without an ocular phenotype. OTX2 truncations led to significant transactivation reduction. A missense variant was identified in another patient without eye abnormalities; however, studies revealed it was most likely not causative. In the mouse, truncations proximal to aa219 caused anophthalmia, while distal truncations and the missense variant were tolerated. During human embryogenesis, OTX2 was expressed in the posterior pituitary, retina, ear, thalamus, choroid plexus, and partially in the hypothalamus, but not in the anterior pituitary. Conclusions: OTX2 mutations are rarely associated with hypopituitarism in isolation without eye abnormalities, and may be variably penetrant, even within the same pedigree. Our data suggest that the endocrine phenotypes in patients with OTX2 mutations are of hypothalamic origin.

FGL1 as a modulator of plasma D-dimer levels: Exome-wide marker analysis of plasma tPA, PAI-1, and D-dimer.

BACKGROUND: Use of targeted exome-arrays with common, rare variants and functionally enriched variation has led to discovery of new genes contributing to population variation in risk factors. Plasminogen activator-inhibitor 1 (PAI-1), tissue plasminogen activator (tPA), and the plasma product D-dimer are important components of the fibrinolytic system. There have been few large-scale genome-wide or exome-wide studies of PAI-1, tPA, and D-dimer. OBJECTIVES: We sought to discover new genetic loci contributing to variation in these traits using an exome-array approach. METHODS: Cohort-level analyses and fixed effects meta-analyses of PAI-1 (n = 15 603), tPA (n = 6876,) and D-dimer (n = 19 306) from 12 cohorts of European ancestry with diverse study design were conducted, including single-variant analyses and gene-based burden testing. RESULTS: Five variants located in NME7, FGL1, and the fibrinogen locus, all associated with D-dimer levels, achieved genome-wide significance (P < 5 × 10-8 ). Replication was sought for these 5 variants, as well as 45 well-imputed variants with P < 1 × 10-4 in the discovery using an independent cohort. Replication was observed for three out of the five significant associations, including a novel and uncommon (0.013 allele frequency) coding variant p.Trp256Leu in FGL1 (fibrinogen-like-1) with increased plasma D-dimer levels. Additionally, a candidate-gene approach revealed a suggestive association for a coding variant (rs143202684-C) in SERPINB2, and suggestive associations with consistent effect in the replication analysis include an intronic variant (rs11057830-A) in SCARB1 associated with increased D-dimer levels. CONCLUSION: This work provides new evidence for a role of FGL1 in hemostasis.

The Endoplasmic Reticulum Cargo Receptor SURF4 Facilitates Efficient Erythropoietin Secretion.

Erythropoietin (EPO) stimulates erythroid differentiation and maturation. Though the transcriptional regulation of EPO has been well studied, the molecular determinants of EPO secretion remain unknown. Here, we generated a HEK293T reporter cell line that provides a quantifiable and selectable readout of intracellular EPO levels and performed a genome-scale CRISPR screen that identified SURF4 as an important mediator of EPO secretion. Targeting SURF4 with multiple independent single guide RNAs (sgRNAs) resulted in intracellular accumulation and extracellular depletion of EPO. Both of these phenotypes were rescued by expression of SURF4 cDNA. Additionally, we found that disruption of SURF4 resulted in accumulation of EPO in the endoplasmic reticulum (ER) compartment and that SURF4 and EPO physically interact. Furthermore, SURF4 disruption in Hep3B cells also caused a defect in the secretion of endogenous EPO under conditions mimicking hypoxia, ruling out an artifact of heterologous overexpression. This work demonstrates that SURF4 functions as an ER cargo receptor that mediates the efficient secretion of EPO. Our findings also suggest that modulating SURF4 may be an effective treatment for disorders of erythropoiesis that are driven by aberrant EPO levels. Finally, we show that SURF4 overexpression results in increased secretion of EPO, suggesting a new strategy for more efficient production of recombinant EPO.

Variants in myelin regulatory factor (MYRF) cause autosomal dominant and syndromic nanophthalmos in humans and retinal degeneration in mice.

Nanophthalmos is a rare, potentially devastating eye condition characterized by small eyes with relatively normal anatomy, a high hyperopic refractive error, and frequent association with angle closure glaucoma and vision loss. The condition constitutes the extreme of hyperopia or farsightedness, a common refractive error that is associated with strabismus and amblyopia in children. NNO1 was the first mapped nanophthalmos locus. We used combined pooled exome sequencing and strong linkage data in the large family used to map this locus to identify a canonical splice site alteration upstream of the last exon of the gene encoding myelin regulatory factor (MYRF c.3376-1G>A), a membrane bound transcription factor that undergoes autoproteolytic cleavage for nuclear localization. This variant produced a stable RNA transcript, leading to a frameshift mutation p.Gly1126Valfs*31 in the C-terminus of the protein. In addition, we identified an early truncating MYRF frameshift mutation, c.769dupC (p.S264QfsX74), in a patient with extreme axial hyperopia and syndromic features. Myrf conditional knockout mice (CKO) developed depigmentation of the retinal pigment epithelium (RPE) and retinal degeneration supporting a role of this gene in retinal and RPE development. Furthermore, we demonstrated the reduced expression of Tmem98, another known nanophthalmos gene, in Myrf CKO mice, and the physical interaction of MYRF with TMEM98. Our study establishes MYRF as a nanophthalmos gene and uncovers a new pathway for eye growth and development.

Extended regions of suspected mis-assembly in the rat reference genome.

We performed whole-genome sequencing for eight inbred rat strains commonly used in genetic mapping studies. They are the founders of the NIH heterogeneous stock (HS) outbred colony. We provide their sequences and variant calls to the rat genomics community. When analyzing the variant calls we identified regions with unusually high levels of heterozygosity. These regions are consistent across the eight inbred strains, including Brown Norway, which is the basis of the rat reference genome. These regions show higher read depths than other regions in the genome and contain higher rates of apparent tri-allelic variant sites. The evidence suggests that these regions may correspond to duplicated segments that were incorrectly overlaid as a single segment in the reference genome. We provide masks for these regions of suspected mis-assembly as a resource for the community to flag potentially false interpretations of mapping or functional results.

Analyses of LMNA-negative juvenile progeroid cases confirms biallelic POLR3A mutations in Wiedemann-Rautenstrauch-like syndrome and expands the phenotypic spectrum of PYCR1 mutations.

Juvenile segmental progeroid syndromes are rare, heterogeneous disorders characterized by signs of premature aging affecting more than one tissue or organ starting in childhood. Hutchinson-Gilford progeria syndrome (HGPS), caused by a recurrent de novo synonymous LMNA mutation resulting in aberrant splicing and generation of a mutant product called progerin, is a prototypical example of such disorders. Here, we performed a joint collaborative study using massively parallel sequencing and targeted Sanger sequencing, aimed at delineating the underlying genetic cause of 14 previously undiagnosed, clinically heterogeneous, non-LMNA-associated juvenile progeroid patients. The molecular diagnosis was achieved in 11 of 14 cases (~ 79%). Furthermore, we firmly establish biallelic mutations in POLR3A as the genetic cause of a recognizable, neonatal, Wiedemann-Rautenstrauch-like progeroid syndrome. Thus, we suggest that POLR3A mutations are causal for a portion of under-diagnosed early-onset segmental progeroid syndromes. We additionally expand the clinical spectrum associated with PYCR1 mutations by showing that they can somewhat resemble HGPS in the first year of life. Moreover, our results lead to clinical reclassification in one single case. Our data emphasize the complex genetic and clinical heterogeneity underlying progeroid disorders.

Mutations in VPS13D lead to a new recessive ataxia with spasticity and mitochondrial defects.

OBJECTIVE: To identify novel causes of recessive ataxias, including spinocerebellar ataxia with saccadic intrusions, spastic ataxias, and spastic paraplegia. METHODS: In an international collaboration, we independently performed exome sequencing in 7 families with recessive ataxia and/or spastic paraplegia. To evaluate the role of VPS13D mutations, we evaluated a Drosophila knockout model and investigated mitochondrial function in patient-derived fibroblast cultures. RESULTS: Exome sequencing identified compound heterozygous mutations in VPS13D on chromosome 1p36 in all 7 families. This included a large family with 5 affected siblings with spinocerebellar ataxia with saccadic intrusions (SCASI), or spinocerebellar ataxia, recessive, type 4 (SCAR4). Linkage to chromosome 1p36 was found in this family with a logarithm of odds score of 3.1. The phenotypic spectrum in our 12 patients was broad. Although most presented with ataxia, additional or predominant spasticity was present in 5 patients. Disease onset ranged from infancy to 39 years, and symptoms were slowly progressive and included loss of independent ambulation in 5. All but 2 patients carried a loss-of-function (nonsense or splice site) mutation on one and a missense mutation on the other allele. Knockdown or removal of Vps13D in Drosophila neurons led to changes in mitochondrial morphology and impairment in mitochondrial distribution along axons. Patient fibroblasts showed altered morphology and functionality including reduced energy production. INTERPRETATION: Our study demonstrates that compound heterozygous mutations in VPS13D cause movement disorders along the ataxia-spasticity spectrum, making VPS13D the fourth VPS13 paralog involved in neurological disorders. Ann Neurol 2018.

Genetic variants in ADAMTS13 as well as smoking are major determinants of plasma ADAMTS13 levels.

The metalloprotease ADAMTS13 cleaves von Willebrand factor (VWF) in circulating blood, limiting the size of VWF multimers and regulating VWF activity. Abnormal regulation of VWF contributes to bleeding and to thrombotic disorders. ADAMTS13 levels in plasma are highly variable among healthy individuals, although the heritability and the genetic determinants of this variation are unclear. We performed genome-wide association studies of plasma ADAMTS13 concentrations in 3244 individuals from 2 independent cohorts of healthy individuals. The heritability of ADAMTS13 levels was between 59.1% (all individuals) and 83.5% (siblings only), whereas tobacco smoking was associated with a decrease in plasma ADAMTS13 levels. Meta-analysis identified common variants near the ADAMTS13 locus on chromosome 9q34.2 that were significantly associated with ADAMTS13 levels and collectively explained 20.0% of the variance. The top single nucleotide polymorphism (SNP), rs28673647, resides in an intron of ADAMTS13 (ß, 6.7%; P = 1.3E-52). Conditional analysis revealed 3 additional independent signals represented by rs3739893 (ß, -22.3%; P = 1.2E-30) and rs3124762 (ß, 3.5%; P = 8.9E-9) close to ADAMTS13 and rs4075970 (ß, 2.4%; P = 6.8E-9) on 21q22.3. Linkage analysis also identified the region around ADAMTS13 (9q34.2) as the top signal (LOD 3.5), consistent with our SNP association analyses. Two nonsynonymous ADAMTS13 variants in the top 2 independent linkage disequilibrium blocks (Q448E and A732V) were identified and characterized in vitro. This study uncovered specific common genetic polymorphisms that are key genetic determinants of the variation in plasma ADAMTS13 levels in healthy individuals.

Genome-wide analysis of central corneal thickness in primary open-angle glaucoma cases in the NEIGHBOR and GLAUGEN consortia.

PURPOSE: To investigate the effects of central corneal thickness (CCT)-associated variants on primary open-angle glaucoma (POAG) risk using single nucleotide polymorphisms (SNP) data from the Glaucoma Genes and Environment (GLAUGEN) and National Eye Institute (NEI) Glaucoma Human Genetics Collaboration (NEIGHBOR) consortia. METHODS: A replication analysis of previously reported CCT SNPs was performed in a CCT dataset (n = 1117) and these SNPs were then tested for association with POAG using a larger POAG dataset (n = 6470). Then a CCT genome-wide association study (GWAS) was performed. Top SNPs from this analysis were selected and tested for association with POAG. cDNA libraries from fetal and adult brain and ocular tissue samples were generated and used for candidate gene expression analysis. RESULTS: Association with one of 20 previously published CCT SNPs was replicated: rs12447690, near the ZNF469 gene (P = 0.001; ß = -5.08 µm/allele). None of these SNPs were significantly associated with POAG. In the CCT GWAS, no SNPs reached genome-wide significance. After testing 50 candidate SNPs for association with POAG, one SNP was identified, rs7481514 within the neurotrimin (NTM) gene, that was significantly associated with POAG in a low-tension subset (P = 0.00099; Odds Ratio [OR] = 1.28). Additionally, SNPs in the CNTNAP4 gene showed suggestive association with POAG (top SNP = rs1428758; P = 0.018; OR = 0.84). NTM and CNTNAP4 were shown to be expressed in ocular tissues. CONCLUSIONS: The results suggest previously reported CCT loci are not significantly associated with POAG susceptibility. By performing a quantitative analysis of CCT and a subsequent analysis of POAG, SNPs in two cell adhesion molecules, NTM and CNTNAP4, were identified and may increase POAG susceptibility in a subset of cases.

Target concentration dependence of DNA melting temperature on oligonucleotide microarrays.

The design of microarrays is currently based on studies focusing on DNA hybridization reaction in bulk solution. However, the presence of a surface to which the probe strand is attached can make the solution-based approximations invalid, resulting in sub-optimum hybridization conditions. To determine the effect of surfaces on DNA duplex formation, the authors studied the dependence of DNA melting temperature (T(m)) on target concentration. An automated system was developed to capture the melting profiles of a 25-mer perfect-match probe-target pair initially hybridized at 23°C. Target concentrations ranged from 0.0165 to 15 nM with different probe amounts (0.03-0.82 pmol on a surface area of 10(18) Å(2)), a constant probe density (5 × 10(12) molecules/cm(2)) and spacer length (15 dT). The authors found that T(m) for duplexes anchored to a surface is lower than in-solution, and this difference increases with increasing target concentration. In a representative set, a target concentration increase from 0.5 to 15 nM with 0.82 pmol of probe on the surface resulted in a T(m) decrease of 6°C when compared with a 4°C increase in solution. At very low target concentrations, a multi-melting process was observed in low temperature domains of the curves. This was attributed to the presence of truncated or mismatch probes.