topr: an R package for viewing and annotating genetic association results.

BACKGROUND: The successful identification of genetic loci for complex traits in genome-wide association studies (GWAS) has resulted in thousands of GWAS summary statistics becoming publicly available for hundreds of complex traits from multiple cohorts and studies. Visualisation is an important aid for interpreting, comparing, validating, and obtaining an overview of large amounts of data. However, the current software is limited in its ability and flexibility to annotate and simultaneously display multiple GWAS results which is useful when interpreting and comparing association results. Therefore, I created the topr R package to facilitate visualisation, annotation, and comparisons of single or multiple GWAS results. It contains functions tailored for viewing and analysing GWAS results. RESULTS: topr provides a fast and elegant visual display of association results, along with the annotation of association peaks with their nearest gene. Association results from multiple analyses can be viewed simultaneously over the entire genome or in a more detailed regional view along with gene information. Users can perform the essential steps of visually exploring and annotating association results and generating elegant publication-ready plots. CONCLUSIONS: topr is developed as a package for the R statistical computing environment, released under the GNU General Public License, and is freely available on the Comprehensive R Archive Network ( http://cran.r-project.org/package=topr ). The source code is available at GitHub ( https://github.com/totajuliusd/topr ). topr provides several advantages and advances over the current alternatives, particularly in its gene annotation functionality and customisable display of single- or multiple-association results. With topr, I provide a flexible tool with multiple features to aid in the analysis and evaluation of GWAS association results.

Genetic effects on the timing of parturition and links to fetal birth weight.

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Genes encoding agrin (AGRN) and neurotrypsin (PRSS12) are associated with muscle mass, strength and plasma C-terminal agrin fragment concentration.

Although physiological data suggest that neuromuscular junction (NMJ) dysfunction is a principal mechanism underpinning sarcopenia, genetic studies have implicated few genes involved in NMJ function. Accordingly, we explored whether genes encoding agrin (AGRN) and neurotrypsin (PRSS12) were associated with sarcopenia phenotypes: muscle mass, strength and plasma C-terminal agrin fragment (CAF). PhenoScanner was used to determine if AGRN and/or PRSS12 variants had previously been implicated with sarcopenia phenotypes. For replication, we combined genotype from whole genome sequencing with phenotypic data from 6715 GenoFit participants aged 18-83 years. Dual energy X-ray absorptiometry assessed whole body lean mass (WBLM) and appendicular lean mass (ALM), hand dynamometry determined grip strength and ELISA measured plasma CAF in a subgroup (n = 260). Follow-up analyses included eQTL analyses, carrier analyses, single-variant and gene-burden tests. rs2710873 (AGRN) and rs71608359 (PRSS12) associate with muscle mass and strength phenotypes, respectively, in the UKBB (p = 8.9 × 10-6 and p = 8.4 × 10-6) and GenoFit cohort (p = 0.019 and p = 0.014). rs2710873 and rs71608359 are eQTLs for AGRN and PRSS12, respectively, in ≥ three tissues. Compared to non-carriers, carriers of rs2710873 had 4.0% higher WBLM and ALM (both p < 0.001), and 9.5% lower CAF concentrations (p < 0.001), while carriers of rs71608359 had 2.3% lower grip strength (p = 0.034). AGRN and PRSS12 are associated with muscle strength and mass in single-variant analyses, while PRSS12 has further associations with muscle strength in gene-burden tests. Our findings provide novel evidence of the relevance of AGRN and PRSS12 to sarcopenia phenotypes and support existing physiological data illustrating the importance of the NMJ in maintaining muscle health during ageing.

Distinction between the effects of parental and fetal genomes on fetal growth.

Birth weight is a common measure of fetal growth that is associated with a range of health outcomes. It is directly affected by the fetal genome and indirectly by the maternal genome. We performed genome-wide association studies on birth weight in the genomes of the child and parents and further analyzed birth length and ponderal index, yielding a total of 243 fetal growth variants. We clustered those variants based on the effects of transmitted and nontransmitted alleles on birth weight. Out of 141 clustered variants, 22 were consistent with parent-of-origin-specific effects. We further used haplotype-specific polygenic risk scores to directly test the relationship between adult traits and birth weight. Our results indicate that the maternal genome contributes to increased birth weight through blood-glucose-raising alleles while blood-pressure-raising alleles reduce birth weight largely through the fetal genome.

The genetic architecture of age-related hearing impairment revealed by genome-wide association analysis.

Age-related hearing impairment (ARHI) is the most common sensory disorder in older adults. We conducted a genome-wide association meta-analysis of 121,934 ARHI cases and 591,699 controls from Iceland and the UK. We identified 21 novel sequence variants, of which 13 are rare, under either additive or recessive models. Of special interest are a missense variant in LOXHD1 (MAF = 1.96%) and a tandem duplication in FBF1 covering 4 exons (MAF = 0.22%) associating with ARHI (OR = 3.7 for homozygotes, P = 1.7 × 10-22 and OR = 4.2 for heterozygotes, P = 5.7 × 10-27, respectively). We constructed an ARHI genetic risk score (GRS) using common variants and showed that a common variant GRS can identify individuals at risk comparable to carriers of rare high penetrance variants. Furthermore, we found that ARHI and tinnitus share genetic causes. This study sheds a new light on the genetic architecture of ARHI, through several rare variants in both Mendelian deafness genes and genes not previously linked to hearing.

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease.

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable1. Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported2-7. A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10-24). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10-4), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10-4) and coeliac disease (OR = 1.62, P = 1.20 × 10-4). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia8 with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10-3). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Eighty-eight variants highlight the role of T cell regulation and airway remodeling in asthma pathogenesis.

Asthma is one of the most common chronic diseases affecting both children and adults. We report a genome-wide association meta-analysis of 69,189 cases and 702,199 controls from Iceland and UK biobank. We find 88 asthma risk variants at 56 loci, 19 previously unreported, and evaluate their effect on other asthma and allergic phenotypes. Of special interest are two low frequency variants associated with protection against asthma; a missense variant in TNFRSF8 and 3' UTR variant in TGFBR1. Functional studies show that the TNFRSF8 variant reduces TNFRSF8 expression both on cell surface and in soluble form, acting as loss of function. eQTL analysis suggests that the TGFBR1 variant acts through gain of function and together with an intronic variant in a downstream gene, SMAD3, points to defective TGFßR1 signaling as one of the biological perturbations increasing asthma risk. Our results increase the number of asthma variants and implicate genes with known role in T cell regulation, inflammation and airway remodeling in asthma pathogenesis.

A homozygous loss-of-function mutation leading to CYBC1 deficiency causes chronic granulomatous disease.

Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase's main subunit (gp91phox) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chip-genotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10-8; OR = 67.6), as well as reduced height (P = 3.3 × 10-4; -8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction.

Toppar: an interactive browser for viewing association study results.

Summary: Data integration and visualization help geneticists make sense of large amounts of data. To help facilitate interpretation of genetic association data we developed Toppar, a customizable visualization tool that stores results from association studies and enables browsing over multiple results, by combining features from existing tools and linking to appropriate external databases. Availability and implementation: Detailed information on Toppar's features and functionality are on our website http://mccarthy.well.ox.ac.uk/toppar/docs along with instructions on how to download, install and run Toppar. Our online version of Toppar is accessible from the website and can be test-driven using Firefox, Safari or Chrome on sub-sets of publicly available genome-wide association study anthropometric waist and body mass index data (Locke et al., 2015; Shungin et al., 2015) from the Genetic Investigation of ANthropometric Traits consortium. Contact: totajuliusd@gmail.com.

The Genetic Landscape of Renal Complications in Type 1 Diabetes.

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

POPE--a tool to aid high-throughput phylogenetic analysis.

UNLABELLED: POPE (Phylogeny, Ortholog and Paralog Extractor) provides an integrated platform for automatic ortholog identification. Intermediate steps can be visualized, modified and analyzed in order to assess and improve the underlying quality of orthology and paralogy assignments. AVAILABILITY: POPE is available for download from the website: http://www.well.ox.ac.uk/~tota/pope.

Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida.

Deuterostomes comprise vertebrates, the related invertebrate chordates (tunicates and cephalochordates) and three other invertebrate taxa: hemichordates, echinoderms and Xenoturbella. The relationships between invertebrate and vertebrate deuterostomes are clearly important for understanding our own distant origins. Recent phylogenetic studies of chordate classes and a sea urchin have indicated that urochordates might be the closest invertebrate sister group of vertebrates, rather than cephalochordates, as traditionally believed. More remarkable is the suggestion that cephalochordates are closer to echinoderms than to vertebrates and urochordates, meaning that chordates are paraphyletic. To study the relationships among all deuterostome groups, we have assembled an alignment of more than 35,000 homologous amino acids, including new data from a hemichordate, starfish and Xenoturbella. We have also sequenced the mitochondrial genome of Xenoturbella. We support the clades Olfactores (urochordates and vertebrates) and Ambulacraria (hemichordates and echinoderms). Analyses using our new data, however, do not support a cephalochordate and echinoderm grouping and we conclude that chordates are monophyletic. Finally, nuclear and mitochondrial data place Xenoturbella as the sister group of the two ambulacrarian phyla. As such, Xenoturbella is shown to be an independent phylum, Xenoturbellida, bringing the number of living deuterostome phyla to four.