Search

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

BackgroundThousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Moreover, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed.MethodsTo derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including n = 229 African American and n = 381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, n = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (n = 27,955) and Hispanic/Latino (n = 28,324) ancestry participants.ResultsOur results revealed 24 suggestive signals (p < 1 × 10-4) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (n = 802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN.ConclusionsThese analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel.

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders1. They are heritable2,3 and etiologically related4,5 behaviors that have been resistant to gene discovery efforts6-11. In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures.

Body mass index (BMI), a proxy measure for obesity, is determined by both environmental (including ethnicity, age, and sex) and genetic factors, with > 400 BMI-associated loci identified to date. However, the impact, interplay, and underlying biological mechanisms among BMI, environment, genetics, and ancestry are not completely understood. To further examine these relationships, we utilized 427,509 calendar year-averaged BMI measurements from 100,418 adults from the single large multiethnic Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We observed substantial independent ancestry and nationality differences, including ancestry principal component interactions and nonlinear effects. To increase the list of BMI-associated variants before assessing other differences, we conducted a genome-wide association study (GWAS) in GERA, with replication in the Genetic Investigation of Anthropomorphic Traits (GIANT) consortium combined with the UK Biobank (UKB), followed by GWAS in GERA combined with GIANT, with replication in the UKB. We discovered 30 novel independent BMI loci (P < 5.0 × 10-8) that replicated. We then assessed the proportion of BMI variance explained by sex in the UKB using previously identified loci compared to previously and newly identified loci and found slight increases: from 3.0 to 3.3% for males and from 2.7 to 3.0% for females. Further, the variance explained by previously and newly identified variants decreased with increasing age in the GERA and UKB cohorts, echoed in the variance explained by the entire genome, which also showed gene-age interaction effects. Finally, we conducted a tissue expression QTL enrichment analysis, which revealed that GWAS BMI-associated variants were enriched in the cerebellum, consistent with prior work in humans and mice.

Primary open-angle glaucoma (POAG) is a leading cause of irreversible vision loss, yet much of the genetic risk remains unaccounted for, especially in African-Americans who have a higher risk for developing POAG. We conduct a multiethnic genome-wide association study (GWAS) of POAG in the GERA cohort, with replication in the UK Biobank (UKB), and vice versa, GWAS in UKB with replication in GERA. We identify 24 loci (P

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma, a leading cause of blindness. IOP heritability has been estimated to up to 67%, and to date only 11 IOP loci have been reported, accounting for 1.5% of IOP variability. Here, we conduct a genome-wide association study of IOP in 69,756 untreated individuals of European, Latino, Asian, and African ancestry. Multiple longitudinal IOP measurements were collected through electronic health records and, in total, 356,987 measurements were included. We identify 47 genome-wide significant IOP-associated loci (P

BackgroundFamilial cerebral cavernous malformation type 1 (CCM1) is an autosomal dominant disease caused by mutations in the Krev Interaction Trapped 1 (KRIT1/CCM1) gene, and characterized by multiple brain lesions that often result in intracerebral hemorrhage (ICH), seizures, and neurological deficits. Carriers of the same genetic mutation can present with variable symptoms and severity of disease, suggesting the influence of modifier factors. Evidence is emerging that inflammation and immune response play a role in the pathogenesis of CCM. The purpose of this study was to investigate whether common variants in inflammatory and immune response genes influence the severity of familial CCM1 disease, as manifested by ICH and greater brain lesion count.MethodsHispanic CCM1 patients (n=188) harboring the founder Q455X 'common Hispanic mutation' (CHM) in the KRIT1 gene were analyzed at baseline. Participants were enrolled between June 2010 and March 2014 either through the Brain Vascular Malformation Consortium (BVMC) study or through the Angioma Alliance organization. Clinical assessment and cerebral susceptibility-weighted magnetic resonance imaging were performed to determine ICH as well as total and large (≥5 mm in diameter) lesion counts. Samples were genotyped on the Affymetrix Axiom Genome-Wide LAT1 Human Array. We analyzed 830 variants in 56 inflammatory and immune response genes for association with ICH and residuals of log-transformed total or large lesion count adjusted for age at enrollment and gender. Variants were analyzed individually or grouped by sub-pathways or whole pathways.ResultsAt baseline, 30.3% of CCM1-CHM subjects had ICH, with a mean ± standard deviation (SD) of 60.1±115.0 (range 0-713) for total lesions and 4.9±8.7 (range 0-104) for large lesions. The heritability estimates explained by all autosomal variants were 0.20 (SE=0.31), 0.81 (SE=0.17), and 0.48 (SE=0.19), for ICH, total lesion count, and large lesion count, respectively. TGFBR2 rs9823731 was significantly associated with ICH as well as with the total and large lesion counts (p≤0.017). Further, IL-4 rs9327638, CD14 rs778588, IL-6R rs114660934 and MSR1 rs62489577 were associated with two markers of disease severity. Finally, the whole pathway was associated with total lesion count (p=0.005) with TLR-4 rs10759930, CD14 rs778588, IL-6R rs114660934 and IGH rs57767447 mainly bearing this association. Eicosanoid signaling, extracellular pattern recognition, and immune response sub-pathways were also associated with the total lesion count.ConclusionsThese results suggest that polymorphisms in inflammatory and immune response pathways contribute to variability in CCM1 disease severity and might be used as predictors of disease severity. In particular, TGFBR2 rs9823731 was associated with all three markers of CCM1 disease severity tested, suggesting that TGFBR2 might be a key participant in the mechanism underlying CCM1 disease severity and phenotype variability. However, further longitudinal studies in larger sample sizes are needed to confirm these findings.

Abstract Background Considering that a portion of the heterogeneity amongst previous replication studies may be due to a variable proportion of obese subjects in case-control designs, we assessed the association of genetic variants with type 2 diabetes (T2D) in large groups of obese and non-obese subjects. Methods We genotyped RETN, KCNJ11, HNF4A, HNF1A, GCK, SLC30A8, ENPP1, ADIPOQ, PPARG, and TCF7L2 polymorphisms in 1,283 normoglycemic (NG) and 1,581 T2D obese individuals as well as in 3,189 NG and 1,244 T2D non-obese subjects of European descent, allowing us to examine T2D risk over a wide range of BMI. Results Amongst non-obese individuals, we observed significant T2D associations with HNF1A I27L [odds ratio (OR) = 1.14, P = 0.04], GCK -30G>A (OR = 1.23, P = 0.01), SLC30A8 R325W (OR = 0.87, P = 0.04), and TCF7L2 rs7903146 (OR = 1.89, P = 4.5 × 10-23), and non-significant associations with PPARG Pro12Ala (OR = 0.85, P = 0.14), ADIPOQ -11,377C>G (OR = 1.00, P = 0.97) and ENPP1 K121Q (OR = 0.99, P = 0.94). In obese subjects, associations with T2D were detected with PPARG Pro12Ala (OR = 0.73, P = 0.004), ADIPOQ -11,377C>G (OR = 1.26, P = 0.02), ENPP1 K121Q (OR = 1.30, P = 0.003) and TCF7L2 rs7903146 (OR = 1.30, P = 1.1 × 10-4), and non-significant associations with HNF1A I27L (OR = 0.96, P = 0.53), GCK -30G>A (OR = 1.15, P = 0.12) and SLC30A8 R325W (OR = 0.95, P = 0.44). However, a genotypic heterogeneity was only found for TCF7L2 rs7903146 (P = 3.2 × 10-5) and ENPP1 K121Q (P = 0.02). No association with T2D was found for KCNJ11, RETN, and HNF4A polymorphisms in non-obese or in obese individuals. Conclusion Genetic variants modulating insulin action may have an increased effect on T2D susceptibility in the presence of obesity, whereas genetic variants acting on insulin secretion may have a greater impact on T2D susceptibility in non-obese individuals.

Abstract Background In adults, the TCF7L2 rs7903146 T allele, commonly associated with type 2 diabetes (T2D), has been also associated with a lower body mass index (BMI) in T2D individuals and with a smaller waist circumference in subjects with impaired glucose tolerance. Methods The present association study aimed at analyzing the contribution of the rs7903146 SNP to smallness for gestational age (SGA) and metabolic profiles in subjects with SGA or appropriate for gestational age birth weight (AGA). Two groups of French Caucasian subjects were selected on birth data: SGA (birth weight < 10th percentile; n = 764), and AGA (25th < birth weight < 75th percentile; n = 627). Family-based association tests were also performed in 3,012 subjects from 628 SGA and AGA pedigrees. Results The rs7903146 genotypic distributions between AGA (30.7%) and SGA (29.0%) were not statistically different (allelic OR = 0.92 [0.78–1.09], p = 0.34). Family association-based studies did not show a distortion of T allele transmission in SGA subjects (p = 0.52). No significant effect of the T allele was detected on any of the metabolic parameters in the SGA group. However, in the AGA group, trends towards a lower insulin secretion (p = 0.03) and a higher fasting glycaemia (p = 0.002) were detected in carriers of the T allele. Conclusion The TCF7L2 rs7903146 variant neither increases the risk for SGA nor modulates birth weight and young adulthood glucose homeostasis in French Caucasian subjects born with SGA.

Background Common variants in the CHRNA5-A3-B4 gene cluster have been shown to be associated with nicotine dependence and alcohol use disorders ( AUDs) and related traits, including the level of response ( LR) to alcohol. Recently, rare variants ( MAF < 0.05) in CHRNB4 have been reported to be associated with a decreased risk of developing nicotine dependence. However, the role of rare variants in the CHRNA5-A3-B4 gene cluster to the LR to alcohol has not yet been established. Methods To determine whether rare variants in the CHRNA5-A3-B4 gene cluster contribute to the LR to alcohol, the coding regions of these 3 genes were sequenced in 538 subjects from the San Diego Sibling Pair study. Results The analyses identified 16 rare missense variants, 9 of which were predicted to be damaging using in silico analysis tools. Carriers of these variants were compared to noncarriers using a family-based design for each gene and for the gene cluster as a whole. In these analyses, a CHRNA5 carrier status was significantly associated with the phenotype related to the feeling of intoxication experienced during the alcohol challenge ( p = 0.039). Conclusions These results indicate that rare genetic variation in the CHRNA5-A3-B4 gene cluster contributes modestly to the LR to alcohol in the San Diego Sibling Pair study and may protect against AUDs. However, replication studies are needed to confirm our findings. [ABSTRACT FROM AUTHOR]

Common PCSK1 variants (notably rs6232 and rs6235) have been shown to be associated with obesity in European, Asian and Mexican populations. To determine whether common PCSK1 variants contribute to obesity in American population, we conducted association analyses in 8,359 subjects using two multi-ethnic American studies: the Coronary Artery Risk Development in Young Adults (CARDIA) study and the Multi-Ethnic Study of Atherosclerosis (MESA). By evaluating the contribution of rs6232 and rs6235 in each ethnic group, we found that in European-American subjects from CARDIA, only rs6232 was associated with BMI (P = 0.006) and obesity (P = 0.018) but also increased the obesity incidence during the 20 years of follow-up (HR = 1.53 [1.07–2.19], P = 0.019). Alternatively, in African-American subjects from CARDIA, rs6235 was associated with BMI (P = 0.028) and obesity (P = 0.018). Further, by combining the two case-control ethnic groups from the CARDIA study in a meta-analysis, association between rs6235 and obesity risk remained significant (OR = 1.23 [1.05–1.45], P = 9.5×10−3). However, neither rs6232 nor rs6235 was associated with BMI or obesity in the MESA study. Interestingly, rs6232 was associated with BMI (P = 4.2×10−3) and obesity (P = 3.4×10−3) in the younger European-American group combining samples from the both studies [less than median age (53 years)], but not among the older age group (P = 0.756 and P = 0.935 for BMI and obesity, respectively). By combining all the case-control ethnic groups from CARDIA and MESA in a meta-analysis, we found no significant association for the both variants and obesity risk. Finally, by exploring the full PCSK1 locus, we observed that no variant remained significant after correction for multiple testing. These results indicate that common PCSK1 variants (notably rs6232 and rs6235) contribute modestly to obesity in multi-ethnic American population. Further, these results suggest that the association of rs6232 with obesity may be age-dependent in European-Americans. However, multiple replication studies in multi-ethnic American population are needed to confirm our findings. [ABSTRACT FROM AUTHOR]