Rare and low-frequency coding variants alter human adult height

Genome-wide compound heterozygote analysis highlights alleles associated with adult height in Europeans

Adult height is the most widely genetically studied common trait in humans; however, the trait variance explainable by currently known height-associated single nucleotide polymorphisms (SNPs) identified from the previous genome-wide association studies (GWAS) is yet far from complete given the high heritability of this complex trait. To exam if compound heterozygotes (CH) may explain extra height variance, we conducted a genome-wide analysis to screen for CH in association with adult height in 10,631 Dutch Europeans enriched with extremely tall people, using our recently developed method implemented in the software package CollapsABEL. The analysis identified six regions (3q23, 5q35.1, 6p21.31, 6p21.33, 7q21.2, and 9p24.3), where multiple pairs of SNPs as CH showed genome-wide significant association with height (P < 1.67 × 10⁻¹⁰). Of those, 9p24.3 represents a novel region influencing adult height, whereas the others have been highlighted in the previous GWAS on height based on analysis of individual SNPs. A replication analysis in 4080 Australians of European ancestry confirmed the significant CH-like association at 9p24.3 (P < 0.05). Together, the collapsed genotypes at these six loci explained 2.51% of the height variance (after adjusting for sex and age), compared with 3.23% explained by the 14 top-associated SNPs at 14 loci identified by traditional GWAS in the same data set (P < 5 × 10⁻⁸). Overall, our study empirically demonstrates that CH plays an important role in adult height and may explain a proportion of its “missing heritability”. Moreover, our findings raise promising expectations for other highly polygenic complex traits to explain missing heritability identifiable through CH-like associations.

Genetics of hypertension: an assessment of progress in the spontaneously hypertensive rat

The application of gene mapping methods to uncover the genetic basis of hypertension in the inbred spontaneously hypertensive rat (SHR) began over 25 yr ago. This animal provides a useful model of genetic high blood pressure, and some of its features are described. In particular, it appears to be a polygenic model of disease, and polygenes participate in human hypertension genetic risk. The SHR hypertension alleles were fixed rapidly by selective breeding in just a few generations and so are presumably common genetic variants present in the outbred Wistar strain from which SHR was created. This review provides a background to the origins and genesis of this rat line. It considers its usefulness as a model organism for a common cardiovascular disease. The progress and obstacles facing mapping are considered in depth, as are the emergence and application of other genome-wide genetic discovery approaches that have been applied to investigate this model. Candidate genes, their identification, and the evidence to support their potential role in blood pressure elevation are considered. The review assesses the progress that has arisen from this work has been limited. Consideration is given to some of the factors that have impeded progress, and prospects for advancing understanding of the genetic basis of hypertension in this model are discussed.

The Rotterdam Study: 2018 update on objectives, design and main results

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, otolaryngological, locomotor, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. Since 2016, the cohort is being expanded by persons aged 40 years and over. The findings of the Rotterdam Study have been presented in over 1500 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.

Association analyses based on false discovery rate implicate new loci for coronary artery disease

Genome-wide association studies (GWAS) in coronary artery disease (CAD) had identified 66 loci at 'genome-wide significance' (P < 5 × 10-8) at the time of this analysis, but a much larger number of putative loci at a false discovery rate (FDR) of 5% (refs. 1,2,3,4). Here we leverage an interim release of UK Biobank (UKBB) data to evaluate the validity of the FDR approach. We tested a CAD phenotype inclusive of angina (SOFT; ncases = 10,801) as well as a stricter definition without angina (HARD; ncases = 6,482) and selected cases with the former phenotype to conduct a meta-analysis using the two most recent CAD GWAS. This approach identified 13 new loci at genome-wide significance, 12 of which were on our previous list of loci meeting the 5% FDR threshold, thus providing strong support that the remaining loci identified by FDR represent genuine signals. The 304 independent variants associated at 5% FDR in this study explain 21.2% of CAD heritability and identify 243 loci that implicate pathways in blood vessel morphogenesis as well as lipid metabolism, nitric oxide signaling and inflammation.

Cholangiocytes in the pathogenesis of primary sclerosing cholangitis and development of cholangiocarcinoma

Primary sclerosing cholangitis (PSC) is an idiopathic cholangiopathy strongly associated with inflammatory bowel disease (IBD) and characterized by cholestasis, chronic immune infiltration and progressive fibrosis of the intrahepatic and extrahepatic bile ducts. PSC confers a high risk of cholangiocarcinoma (CCA) with PSC-CCA representing the leading cause of PSC-associated mortality. PSC-CCA is derived from cholangiocytes and associated progenitor cells - a heterogeneous group of dynamic epithelial cells lining the biliary tree that modulate the composition and volume of bile production by the liver. Infection, inflammation and cholestasis can trigger cholangiocyte activation leading to an increased expression of adhesion and antigen-presenting molecules as well as the release of various inflammatory and fibrogenic mediators. As a result, activated cholangiocytes engage in a myriad of cellular processes, including hepatocellular proliferation, apoptosis, angiogenesis and fibrosis. Cholangiocytes can also regulate the recruitment of immune cells, mesenchymal cells, and endothelial cells that participate in tissue repair and destruction in settings of persistent inflammation. In PSC, the role of cholangiocytes and the mechanisms governing their transformation to PSC-CCA are unclear however localization of disease suggests that cholangiocytes are a key target and potential regulator of hepatobiliary immunity, fibrogenesis and tumorigenesis. Herein, we summarize mechanisms of cholangiocyte activation in PSC and highlight new insights into disease pathways that may contribute to the development of PSC-CCA. This article is part of a Special Issue entitled: Cholangiocytes in Health and Disease edited by Jesus Banales, Marco Marzioni, Nicholas LaRusso and Peter Jansen.

Imaging IGF-I uptake in growth plate cartilage using in vivo multiphoton microscopy

Bones elongate through endochondral ossification in cartilaginous growth plates located at ends of primary long bones. Linear growth ensues from a cascade of biochemical signals initiated by actions of systemic and local regulators on growth plate chondrocytes. Although cellular processes are well defined, there is a fundamental gap in understanding how growth regulators are physically transported from surrounding blood vessels into and through dense, avascular cartilage matrix. Intravital imaging using in vivo multiphoton microscopy is one promising strategy to overcome this barrier by quantitatively tracking molecular delivery to cartilage from the vasculature in real time. We previously used in vivo multiphoton imaging to show that hindlimb heating increases vascular access of large molecules to growth plates using 10-, 40-, and 70-kDa dextran tracers. To comparatively evaluate transport of similarly sized physiological regulators, we developed and validated methods for measuring uptake of biologically active IGF-I into proximal tibial growth plates of live 5-wk-old mice. We demonstrate that fluorescently labeled IGF-I (8.2 kDa) is readily taken up in the growth plate and localizes to chondrocytes. Bioactivity tests performed on cultured metatarsal bones confirmed that the labeled protein is functional, assessed by phosphorylation of its signaling kinase, Akt. This methodology, which can be broadly applied to many different proteins and tissues, is relevant for understanding factors that affect delivery of biologically relevant molecules to the skeleton in real time. Results may lead to the development of drug-targeting strategies to treat a wide range of bone and cartilage pathologies.NEW & NOTEWORTHY This paper describes and validates a novel method for imaging transport of biologically active, fluorescently labeled IGF-I into skeletal growth plates of live mice using multiphoton microscopy. Cellular patterns of fluorescence in the growth plate were completely distinct from our prior publications using biologically inert probes, demonstrating for the first time in vivo localization of IGF-I in chondrocytes and perichondrium. These results form important groundwork for future studies aimed at targeting therapeutics into growth plates.

Rare coding variants pinpoint genes that control human hematological traits

The identification of rare coding or splice site variants remains the most straightforward strategy to link genes with human phenotypes. Here, we analyzed the association between 137,086 rare (minor allele frequency (MAF) <1%) coding or splice site variants and 15 hematological traits in up to 308,572 participants. We found 56 such rare coding or splice site variants at P<5x10-8, including 31 that are associated with a blood-cell phenotype for the first time. All but one of these 31 new independent variants map to loci previously implicated in hematopoiesis by genome-wide association studies (GWAS). This includes a rare splice acceptor variant (rs146597587, MAF = 0.5%) in interleukin 33 (IL33) associated with reduced eosinophil count (P = 2.4x10-23), and lower risk of asthma (P = 2.6x10-7, odds ratio [95% confidence interval] = 0.56 [0.45-0.70]) and allergic rhinitis (P = 4.2x10-4, odds ratio = 0.55 [0.39-0.76]). The single new locus identified in our study is defined by a rare p.Arg172Gly missense variant (rs145535174, MAF = 0.05%) in plasminogen (PLG) associated with increased platelet count (P = 6.8x10-9), and decreased D-dimer concentration (P = 0.018) and platelet reactivity (P<0.03). Finally, our results indicate that searching for rare coding or splice site variants in very large sample sizes can help prioritize causal genes at many GWAS loci associated with complex human diseases and traits.

Rare genetic variants in CX3CR1 and their contribution to the increased risk of schizophrenia and autism spectrum disorders

CX3CR1, a G protein-coupled receptor solely expressed by microglia in the brain, has been repeatedly reported to be associated with neurodevelopmental disorders including schizophrenia (SCZ) and autism spectrum disorders (ASD) in transcriptomic and animal studies but not in genetic studies. To address the impacts of variants in CX3CR1 on neurodevelopmental disorders, we conducted coding exon-targeted resequencing of CX3CR1 in 370 Japanese SCZ and 192 ASD patients using next-generation sequencing technology, followed by a genetic association study in a sample comprising 7054 unrelated individuals (2653 SCZ, 574 ASD and 3827 controls). We then performed in silico three-dimensional (3D) structural modeling and in vivo disruption of Akt phosphorylation to determine the impact of the detected variant on CX3CR1-dependent signal transduction. We detected a statistically significant association between the variant Ala55Thr in CX3CR1 with SCZ and ASD phenotypes (odds ratio=8.3, P=0.020). A 3D structural model indicated that Ala55Thr could destabilize the conformation of the CX3CR1 helix 8 and affect its interaction with a heterotrimeric G protein. In vitro functional analysis showed that the CX3CR1-Ala55Thr mutation inhibited cell signaling induced by fractalkine, the ligand for CX3CR1. The combined data suggested that the variant Ala55Thr in CX3CR1 might result in the disruption of CX3CR1 signaling. Our results strengthen the association between microglia-specific genes and neurodevelopmental disorders.

Contribution of rare and low-frequency whole-genome sequence variants to complex traits variation in dairy cattle

Background

Whole-genome sequencing and imputation methodologies have enabled the study of the effects of genomic variants with low to very low minor allele frequency (MAF) on variation in complex traits. Our objective was to estimate the proportion of variance explained by imputed sequence variants classified according to their MAF compared with the variance explained by the pedigree-based additive genetic relationship matrix for 17 traits in Nordic Holstein dairy cattle.

Results

Imputed sequence variants were grouped into seven classes according to their MAF (0.001–0.01, 0.01–0.05, 0.05–0.1, 0.1–0.2, 0.2–0.3, 0.3–0.4 and 0.4–0.5). The total contribution of all imputed sequence variants to variance in deregressed estimated breeding values or proofs (DRP) for different traits ranged from 0.41 [standard error (SE) = 0.026] for temperament to 0.87 (SE = 0.011) for milk yield. The contribution of rare variants (MAF < 0.01) to the total DRP variance explained by all imputed sequence variants was relatively small (a maximum of 12.5% for the health index). Rare and low-frequency variants (MAF < 0.05) contributed a larger proportion of the explained DRP variances (>13%) for health-related traits than for production traits (<11%). However, a substantial proportion of these variance estimates across different MAF classes had large SE, especially when the variance explained by a MAF class was small. The proportion of DRP variance that was explained by all imputed whole-genome sequence variants improved slightly compared with variance explained by the 50 k Illumina markers, which are routinely used in bovine genomic prediction. However, the proportion of DRP variance explained by imputed sequence variants was lower than that explained by pedigree relationships, ranging from 1.5% for milk yield to 37.9% for the health index.

Conclusions

Imputed sequence variants explained more of the variance in DRP than the 50 k markers for most traits, but explained less variance than that captured by pedigree-based relationships. Although in humans partitioning variants into groups based on MAF and linkage disequilibrium was used to estimate heritability without bias, many of our bovine estimates had a high SE. For a reliable estimate of the explained DRP variance for different MAF classes, larger sample sizes are needed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0336-z) contains supplementary material, which is available to authorized users.

Modeling human disease in rodents by CRISPR/Cas9 genome editing

Modeling human disease has proven to be a challenge for the scientific community. For years, generating an animal model was complicated and restricted to very few species. With the rise of CRISPR/Cas9, it is now possible to generate more or less any animal model. In this review, we will show how this technology is and will change our way to obtain relevant disease animal models and how it should impact human health.

The role of human host genetics in tuberculosis resistance

INTRODUCTION

Tuberculosis (TB) remains a public health problem: the latest estimate of new incident cases per year is a staggering 10.4 million. Despite this overwhelming number, the majority of the immunocompetent population can control infection with Mycobacterium tuberculosis. The human genome underlies the immune response and contributes to the outcome of TB infection. Areas covered: Investigations of TB resistance in the general population have closely mirrored those of other infectious diseases and initially involved epidemiological observations. Linkage and association studies, including studies of VDR, SLC11A1 and HLA-DRB1 followed. Genome-wide association studies of common variants, not necessarily sufficient for disease, became possible after technological advancements. Other approaches involved the identification of those individuals with rare disease-causing mutations that strongly predispose to TB, epistasis and the role of ethnicity in disease. Despite these efforts, infection outcome, on an individual basis, cannot yet be predicted. Expert commentary: The early identification of future disease progressors is necessary to stem the TB epidemic. Human genetics may contribute to this endeavour and could in future suggest pathways to target for disease prevention. This will however require concerted efforts to establish large, well-phenotyped cohorts from different ethnicities, improved genomic resources and a better understanding of the human genome architecture.

10 Years of GWAS Discovery: Biology, Function, and Translation

Application of the experimental design of genome-wide association studies (GWASs) is now 10 years old (young), and here we review the remarkable range of discoveries it has facilitated in population and complex-trait genetics, the biology of diseases, and translation toward new therapeutics. We predict the likely discoveries in the next 10 years, when GWASs will be based on millions of samples with array data imputed to a large fully sequenced reference panel and on hundreds of thousands of samples with whole-genome sequencing data.

Personalized medicine: Genetic risk prediction of drug response

Pharmacogenomics (PGx), a substantial component of "personalized medicine", seeks to understand each individual's genetic composition to optimize drug therapy -- maximizing beneficial drug response, while minimizing adverse drug reactions (ADRs). Drug responses are highly variable because innumerable factors contribute to ultimate phenotypic outcomes. Recent genome-wide PGx studies have provided some insight into genetic basis of variability in drug response. These can be grouped into three categories. [a] Monogenic (Mendelian) traits include early examples mostly of inherited disorders, and some severe (idiosyncratic) ADRs typically influenced by single rare coding variants. [b] Predominantly oligogenic traits represent variation largely influenced by a small number of major pharmacokinetic or pharmacodynamic genes. [c] Complex PGx traits resemble most multifactorial quantitative traits -- influenced by numerous small-effect variants, together with epigenetic effects and environmental factors. Prediction of monogenic drug responses is relatively simple, involving detection of underlying mutations; due to rarity of these events and incomplete penetrance, however, prospective tests based on genotype will have high false-positive rates, plus pharmacoeconomics will require justification. Prediction of predominantly oligogenic traits is slowly improving. Although a substantial fraction of variation can be explained by limited numbers of large-effect genetic variants, uncertainty in successful predictions and overall cost-benefit ratios will make such tests elusive for everyday clinical use. Prediction of complex PGx traits is almost impossible in the foreseeable future. Genome-wide association studies of large cohorts will continue to discover relevant genetic variants; however, these small-effect variants, combined, explain only a small fraction of phenotypic variance -- thus having limited predictive power and clinical utility.

Genome-wide association analysis identifies a GLUL haplotype for familial hepatitis B virus-related hepatocellular carcinoma

BACKGROUND

A family history of liver cancer increases the risk of developing hepatocellular carcinoma (HCC) by 2-fold to 10-fold among patients with chronic hepatitis B virus (HBV). Previous genome-wide association studies have identified many possible susceptible loci associated with sporadic HBV-related HCC. However, despite family history being a well-known risk factor for HBV-related HCC, to the authors' knowledge its genetic mechanisms and associating loci remain largely unknown or unexplored, most likely due to the relative rarity of familial HCC and the difficulty of sample collection.

METHODS

The authors conducted a genome-wide association study with 139 male cases with familial HBV-related HCC and 139 non-HCC male controls with chronic HBV. The results were corroborated further with an independent cohort of 101 patients with familial HBV-related HCC and comparison with both the 1000 Genomes Project and the Taiwan Biobank.

RESULTS

A total of 51 risk single-nucleotide polymorphisms (P≤1E-04) were identified in the association analyses, which included 2 clusters of associated single-nucleotide polymorphisms and haplotypes at 1q25.3 (glutamate-ammonia ligase [GLUL]/transmembrane epididymal protein 1 [TEDDM1]/long intergenic non-protein-coding RNA 272 [LINC00272]/regulator of G-protein signaling-like 1 [RGSL1]) and 17q11.2 (solute carrier family 13 member 2 [SLC13A2]/forkhead box N1 [FOXN1]). Both the GLUL and SLC13A2/FOXN1 haplotypes have large effect sizes and were found to be different from those found from genome-wide association studies of sporadic HCCs.

CONCLUSIONS

To the authors' knowledge, the current study is the first genome-wide association study to identify genetic factors for familial HBV-related HCC. The results identified 2 large effect susceptible haplotypes located at GLUL and SLC13A2/FOXN1. The current study findings also suggest different genetic susceptibility between familial and sporadic HBV-related HCC. Cancer 2017;123:3966-76. © 2017 American Cancer Society.

An Expanded View of Complex Traits: From Polygenic to Omnigenic

A central goal of genetics is to understand the links between genetic variation and disease. Intuitively, one might expect disease-causing variants to cluster into key pathways that drive disease etiology. But for complex traits, association signals tend to be spread across most of the genome-including near many genes without an obvious connection to disease. We propose that gene regulatory networks are sufficiently interconnected such that all genes expressed in disease-relevant cells are liable to affect the functions of core disease-related genes and that most heritability can be explained by effects on genes outside core pathways. We refer to this hypothesis as an "omnigenic" model.

An Expanded View of Complex Traits: From Polygenic to Omnigenic

A central goal of genetics is to understand the links between genetic variation and disease. Intuitively, one might expect disease-causing variants to cluster into key pathways that drive disease etiology. But for complex traits, association signals tend to be spread across most of the genome-including near many genes without an obvious connection to disease. We propose that gene regulatory networks are sufficiently interconnected such that all genes expressed in disease-relevant cells are liable to affect the functions of core disease-related genes and that most heritability can be explained by effects on genes outside core pathways. We refer to this hypothesis as an "omnigenic" model.

The impact of rare and low-frequency genetic variants in common disease

Despite thousands of genetic loci identified to date, a large proportion of genetic variation predisposing to complex disease and traits remains unaccounted for. Advances in sequencing technology enable focused explorations on the contribution of low-frequency and rare variants to human traits. Here we review experimental approaches and current knowledge on the contribution of these genetic variants in complex disease and discuss challenges and opportunities for personalised medicine.