Search

BackgroundHeterogeneity in the long-term metabolic response to dietary macronutrient composition can be partially explained by genetic factors. However, few studies have demonstrated reproducible gene-diet interactions (GDIs), likely due in part to measurement error in dietary intake estimation as well as insufficient capture of rare genetic variation. Discovery analyses in ancestry-diverse cohorts that include rare genetic variants from whole-genome sequencing (WGS) could help identify genetic variants modifying the effects of dietary macronutrient composition on glycemic phenotypes.ObjectiveWe aimed to identify macronutrient GDIs across the genetic frequency spectrum associated with continuous glycemic traits in genetically and culturally diverse cohorts.MethodsWe analyzed N=33,187 diabetes-free participants from 10 cohorts in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program with WGS, self-reported diet, and glycemic traits (fasting glucose [FG], insulin [FI], and hemoglobin A1c [HbA1c]). We fit multivariable-adjusted linear mixed models for the main effect of diet, modeled as an isocaloric substitution of carbohydrate for fat, and for its interactions with genetic variants genome-wide. Tests were performed for both common variants and gene-based rare variant sets in each cohort followed by a combined cohort meta-analysis.ResultsIn main effect models, participants consuming more calories from carbohydrate at the expense of fat had modestly lower glycemic trait values (β per 250 kcal substitution for FG: −0.030 mmol/L, p=2.7×10−6; lnFI: −0.008 log(pmol/L), p=0.17; HbA1c: −0.013 %, p=0.025). In GDI analyses, a common African ancestry-enriched variant (rs79762542; 78 kb upstream of the FRAS1 gene) reached study-wide significance (p = 1.14×10−8) indicating a higher HbA1c with greater proportion of calories from carbohydrate vs. fat among minor allele carriers only. This interaction was replicated in the UK Biobank cohort. Simulations revealed that there is (1) a substantial impact of measurement error on statistical power for GDI discovery at these sample sizes, especially for rare genetic variants, and (2) over 150,000 samples may be necessary to identify similar macronutrient GDIs under realistic assumptions about effect size and measurement error.ConclusionsOur analysis identified a potential genetic interaction modifying the dietary macronutrient-HbA1c association while highlighting the importance of precise exposure measurement and significantly increased sample size to identify additional similar effects.

Plasma lipids are heritable modifiable causal factors for coronary artery disease, the leading cause of death globally. Despite the well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing, partly due to limited sample sizes, ancestral diversity, and interpretation of potential clinical significance. Increasingly larger whole genome sequence datasets with plasma lipids coupled with methodologic advances enable us to more fully catalog the allelic spectrum for lipids. Here, among 66,329 ancestrally diverse (56% non-European ancestry) participants, we associate 428M variants from deep-coverage whole genome sequences with plasma lipids. Approximately 400M of these variants were not studied in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with plasma lipids through analysis of common and rare coding variants. We additionally discover several significantly associated rare non-coding variants largely at Mendelian lipid genes. Notably, we detect rareLDLRintronic variants associated with markedly increased LDL-C, similar to rareLDLRexonic variants. In conclusion, we conducted a systematic whole genome scan for plasma lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite notable anthropometric variability with ancestry proportions, and a high burden of growth stunting and overweight/obesity in Hispanic/Latino populations. This address this knowledge gap, we analyzed densely-imputed genetic data in a sample of Hispanic/Latino adults, to identify and fine-map common genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (Stage 1, n=59,769) and validated our findings in 9 additional studies (HISLA Stage 2, n=9,336). We conducted a trans-ethnic GWAS with summary statistics from HISLA Stage 1 and existing consortia of European and African ancestries. In our HISLA Stage 1+2 analyses, we discovered one novel BMI locus, as well two novel BMI signals and another novel height signal, each within established anthropometric loci. In our trans-ethnic meta- analysis, we identified three additional novel BMI loci, one novel height locus, and one novel WHRadjBMI locus. We also identified three secondary signals for BMI, 28 for height, and two for WHRadjBMI. We replicated >60 established anthropometric loci in Hispanic/Latino populations at genome-wide significance—representing up to 30% of previously-reported index SNP anthropometric associations. Trans-ethnic meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our novel findings demonstrate that future studies may also benefit from leveraging differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

Since 2005, genome-wide association (GWA) datasets have been largely biased toward sampling European ancestry individuals, and recent studies have shown that GWA results estimated from self-identified European individuals are not transferable to non-European individuals due to various confounding challenges. Here, we demonstrate that enrichment analyses which aggregate SNP-level association statistics at multiple genomic scales—from genes to genomic regions and pathways—have been underutilized in the GWA era and can generate biologically interpretable hypotheses regarding the genetic basis of complex trait architecture. We illustrate examples of the robust associations generated by enrichment analyses while studying 25 continuous traits assayed in 566,786 individuals from seven diverse self-identified human ancestries in the UK Biobank and the Biobank Japan, as well as 44,348 admixed individuals from the PAGE consortium including cohorts of African-American, Hispanic and Latin American, Native Hawaiian, and American Indian/Alaska Native individuals. We identify 1,000 gene-level associations that are genome-wide significant in at least two ancestry cohorts across these 25 traits, as well as highly conserved pathway associations with triglyceride levels in European, East Asian, and Native Hawaiian cohorts.

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in Trans-Omics for Precision Medicine (TOPMed) program (N=63,302) and UK Biobank (N=48,658). Bidirectional Mendelian randomization studies were consistent with LTL lengthening increasing propensity to develop CHIP, but CHIP then in turn hastening LTL shortening. We also demonstrated evidence of modest mediation between CHIP and CAD by LTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, and especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aim to: 1) identify novel loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC), and hip circumference (HIP), all adjusted for body mass index (adjBMI), using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); 2) determine if differences in genetic associations differ by background group within HCHS/SOL; 3) determine whether previously reported association regions generalize to HCHS/SOL. Our analyses included 7,472 women and 5,200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban, and Dominican) background residing in the US, with genome-wide array data imputed to the 1000 genomes Phase I multiethnic reference panel. We analyzed associations stratified by sex in addition to sexes combined using linear mixed-model regression. We identified 16 variants for WHRadjBMI, 22 for WCadjBMI, and 28 for HIPadjBMI that reached suggestive significance (P−6). Many of the loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in nine cohort studies (N=34,161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (ancestry-specific P−8 after meta-analysis with HCHS/SOL), including rs13301996 in the sexes-combined analysis, and rs79478137 for women-only for WHRadjBMI; rs28692724 in women-only for HIPadjBMI; and rs3168072 in the sexes combined analysis for WCadjBMI. Also, a total of eight previously reported WHRadjBMI association regions, 12 for HIPadjBMI, and 10 for WCadjBMI generalized to HCHS/SOL. Our study findings highlight the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity-susceptibility that may be ancestry-specific.

Hematological measures are important intermediate clinical phenotypes for many acute and chronic diseases. Hematological measures are highly heritable, and although genome-wide association studies (GWAS) have identified thousands of loci containing trait-associated variants, the causal genes underlying these associations are often uncertain. To better understand the underlying genetic regulatory mechanisms, we performed a transcriptome-wide association study (TWAS) using PrediXcan to systematically investigate the association between genetically-predicted gene expression and hematological measures in 54,542 individuals of European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We found 239 significant gene-trait associations with hematological measures. Among this set of 239 associations, we replicated 71 at p < 0.05 with same direction of effect for the blood cell trait in a meta-analysis of TWAS results consisting of up to 35,900 European ancestry individuals from the Women’s Health Initiative (WHI), the Atherosclerosis Risk in Communities Study (ARIC), and BioMe Biobank. We further attempted to refine this list of candidate genes by performing conditional analyses, adjusting for individual variants previously associated with these hematological measures, and performed further fine-mapping of TWAS loci. To assist with the interpretation of TWAS findings, we designed an R Shiny application to interactively visualize TWAS results, one genomic locus at a time, by integrating our TWAS results with additional genetic data sources (GWAS, TWAS from other gene expression reference panels, conditional analyses, known GWAS variants, etc.). Our results and R Shiny application highlight frequently overlooked challenges with TWAS and illustrate the complexity of TWAS fine-mapping efforts.Author SummaryTranscriptome-wide association studies (TWAS) have shown great promise in furthering our understanding of the genetic regulatory mechanisms underlying complex trait variation. However, interpreting TWAS results can be incredibly complex, especially in large-scale analyses where hundreds of signals appear throughout the genome, with multiple genes often identified in a single chromosomal region. Our research demonstrates this complexity through real data examples from our analysis of hematological traits, and we provide a useful web application to visualize TWAS results in a broadly approachable format. Together, our results and web application illustrate the importance of interpreting TWAS studies in context and highlight the need to carefully examine results in a region-wide context to draw reasonable conclusions and formulate mechanistic hypotheses.

We assembled an ancestrally diverse collection of genome-wide association studies of type 2 diabetes (T2D) in 180,834 cases and 1,159,055 controls (48.9% non-European descent). We identified 277 loci at genome-wide significance (p-8), including 237 attaining a more stringent trans-ancestry threshold (p-9), which were delineated to 338 distinct association signals. Trans-ancestry meta-regression offered substantial enhancements to fine-mapping, with 58.6% of associations more precisely localised due to population diversity, and 54.4% of signals resolved to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying foundations for functional investigations. Trans-ancestry genetic risk scores enhanced transferability across diverse populations, providing a step towards more effective clinical translation to improve global health.

Despite the dramatic underrepresentation of non-European populations in human genetics studies, researchers continue to exclude participants of non-European ancestry, even when these data are available. This practice perpetuates existing research disparities and can lead to important and large effect size associations being missed. Here, we conducted genome-wide association studies (GWAS) of 31 serum and urine biomarker quantitative traits in African (n=9354), East Asian (n=2559) and South Asian (n=9823) UK Biobank participants ancestry. We adjusted for all known GWAS catalog variants for each trait, as well as novel signals identified in European ancestry UK Biobank participants alone. We identify 12 novel signals in African ancestry and 3 novel signals in South Asian participants (p−10). Many of these signals are highly plausible and rare in Europeans (1% or lower minor allele frequency), includingcispQTLs for the genes encoding serum biomarkers like gamma-glutamyl transferase and apolipoprotein A,PIEZ01andG6PDvariants with impacts on HbA1c through likely erythocytic mechanisms, and a coding variant inGPLD1, a gene which cleaves GPI-anchors, associated with normally GPI-anchored protein alkaline phosphatase in serum. This work illustrates the importance of using the genetic data we already have in diverse populations, with many novel discoveries possible in even modest sample sizes.

Glycaemic traits are used to diagnose and monitor type 2 diabetes, and cardiometabolic health. To date, most genetic studies of glycaemic traits have focused on individuals of European ancestry. Here, we aggregated genome-wide association studies in up to 281,416 individuals without diabetes (30% non-European ancestry) with fasting glucose, 2h-glucose post-challenge, glycated haemoglobin, and fasting insulin data. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P-8), 80% with no significant evidence of between-ancestry heterogeneity. Analyses restricted to European ancestry individuals with equivalent sample size would have led to 24 fewer new loci. Compared to single-ancestry, equivalent sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase understanding of diabetes pathophysiology by use of trans-ancestry studies for improved power and resolution.

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase power. Harmonization of the data usually requires substantial effort due to heterogeneity in phenotype definitions, study design, data collection procedures, and data set organization. Here we describe a centralized system for phenotype harmonization that includes input from phenotype domain and study experts, quality control, documentation, reproducible results, and data sharing mechanisms. This system was developed for the National Heart, Lung and Blood Institute’s Trans-Omics for Precision Medicine (TOPMed) program, which is generating genomic and other omics data for >80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants from up to 17 TOPMed studies per phenotype. We discuss the challenges faced in this undertaking and how they were addressed. The harmonized phenotype data and associated documentation have been submitted to National Institutes of Health data repositories for controlled-access by the scientific community. We also provide materials to facilitate future harmonization efforts by the community, which include (1) the code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify or extend these harmonizations to additional studies; and (2) results of labeling thousands of phenotype variables with controlled vocabulary terms.

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.1 Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations was associated with both hematologic cancer incidence2–4 and coronary heart disease prevalence.5 Somatic mutations with leukemogenic potential may confer selective cellular advantages leading to clonal expansion, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).6 Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer. Aggregates of rare germline loss-of-function variants in CHEK2, a DNA damage repair gene, predisposed to CHIP acquisition. Overall, we observe that germline genetic variation altering hematopoietic stem cell function and the fidelity of DNA-damage repair increase the likelihood of somatic mutations leading to CHIP.

We conducted a methylome-wide association study to examine associations between DNA methylation in whole blood and central adiposity and body fat distribution, measured as waist circumference, waist- to-hip ratio, and waist-to-height ratio adjusted for body mass index, in 2684 African American adults in the Atherosclerosis Risk in Communities study. We validated significantly associated Cytosine-phosphate-Guanine methylation sites (CpGs) among adults using the Women’s Health Initiative and Framingham Heart Study participants (combined N=5743) and generalized associations in adolescents from The Raine Study (N=820). We identified 11 CpGs that were robustly associated with one or more central adiposity trait in adults and 2 in adolescents, including CpG site associations near TXNIP, ADCY7, SREBF1, and RAP1GAP2 that had not previously been associated with obesity-related traits.

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., Author summary Admixed African and Hispanic/Latino populations remain understudied in genetic studies of complex diseases. These populations have more complex linkage disequilibrium (LD) structure that can impair mapping of variants. Genotype imputation represents an approach to improve genome coverage, especially for rare or ancestry-specific variation; however, these understudied populations also have smaller relevant imputation reference panels. In this study, we leveraged >100,000 phased sequences generated from the multi-ethnic NHLBI TOPMed project for imputation in ~21,600 individuals of African ancestry (AAs) and ~21,700 Hispanics/Latinos. We demonstrated substantially higher imputation quality for low frequency and rare variants in comparison to the 1000 Genomes Project and Haplotype Reference Consortium reference panels. Analysis of quantitative hematological traits led to the discovery of associations with two rare variants in the HBB gene; one of these variants was replicated in an independent sample, and the other is known to cause anemia in the homozygous state. By comparison, the same HBB variants would not have been genome-wide significant using current reference panels due to lower imputation quality. Our findings demonstrate the power of TOPMed whole genome sequencing data for imputation and subsequent association analysis in admixed African and Hispanic/Latino populations.

Heritability, the proportion of phenotypic variance explained by genetic factors, can be estimated from pedigree data 1, but such estimates are uninformative with respect to the underlying genetic architecture. Analyses of data from genome-wide association studies (GWAS) on unrelated individuals have shown that for human traits and disease, approximately one-third to two-thirds of heritability is captured by common SNPs 2–5. It is not known whether the remaining heritability is due to the imperfect tagging of causal variants by common SNPs, in particular if the causal variants are rare, or other reasons such as overestimation of heritability from pedigree data. Here we show that pedigree heritability for height and body mass index (BMI) appears to be largely recovered from whole-genome sequence (WGS) data on 25,465 unrelated individuals of European ancestry. We assigned 33.7 million genetic variants to groups based upon their minor allele frequencies (MAF) and linkage disequilibrium (LD) with variants nearby, and estimated and partitioned genetic variance accordingly. The estimated heritability was 0.68 (SE 0.10) for height and 0.30 (SE 0.10) for BMI, with a range of ~0.60 – 0.71 for height and ~0.25 – 0.35 for BMI, depending on quality control and analysis strategies. Low-MAF variants in low LD with neighbouring variants were enriched for heritability, to a greater extent for protein-altering variants, consistent with negative selection thereon. Cumulatively variants with 0.0001 < MAF < 0.1 explained 0.47 (SE 0.07) and 0.30 (SE 0.10) of heritability for height and BMI, respectively. Our results imply that rare variants, in particular those in regions of low LD, is a major source of the still missing heritability of complex traits and disease.

Summary paragraphThe Trans-Omics for Precision Medicine (TOPMed) program seeks to elucidate the genetic architecture and disease biology of heart, lung, blood, and sleep disorders, with the ultimate goal of improving diagnosis, treatment, and prevention. The initial phases of the program focus on whole genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here, we describe TOPMed goals and design as well as resources and early insights from the sequence data. The resources include a variant browser, a genotype imputation panel, and sharing of genomic and phenotypic data via dbGaP. In 53,581 TOPMed samples, >400 million single-nucleotide and insertion/deletion variants were detected by alignment with the reference genome. Additional novel variants are detectable through assembly of unmapped reads and customized analysis in highly variable loci. Among the >400 million variants detected, 97% have frequency

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To provide new insights in the biology of sleep-associated adverse lipid profile, we conducted multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identified 49 novel lipid loci, and 10 additional novel lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identified new gene-sleep interactions for known lipid loci such asLPLandPCSK9. The novel gene-sleep interactions had a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explained 4.25% of the variance in triglyceride concentration. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N=60,552), African (N=8,429), Asian (N=9,959), and Hispanic/Latino (N=11,775) ethnicities. We identified over 50 novel loci at genome-wide significance in ancestry-specific and/or multiethnic meta-analyses. Recent fine mapping methods incorporating functional annotation, gene expression, and/or differences in linkage disequilibrium between ethnicities identified potential causal variants and genes at known and newly identified loci. Sixteen of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12.

Summary/AbstractGenome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development, and clinical guidelines. However, the dominance of European-ancestry populations in GWAS creates a biased view of the role of human variation in disease, and hinders the equitable translation of genetic associations into clinical and public health applications. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioral phenotypes in 49,839 non-European individuals. Using strategies designed for analysis of multi-ethnic and admixed populations, we confirm 574 GWAS catalog variants across these traits, and find 38 secondary signals in known loci and 27 novel loci. Our data shows strong evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts, and insights into clinical implications. We strongly advocate for continued, large genome-wide efforts in diverse populations to reduce health disparities.

BackgroundSmoking and alcohol use behaviors in humans have been associated with common genetic variants within multiple genomic loci. Investigation of rare variation within these loci holds promise for identifying causal variants impacting biological mechanisms in the etiology of disordered behavior. Microarrays have been designed to genotype rare nonsynonymous and putative loss of function variants. Such variants are expected to have greater deleterious consequences on gene function than other variants, and significantly contribute to disease risk.MethodsIn the present study, we analyzed ∼250,000 rare variants from 17 independent studies. Each variant was tested for association with five addiction-related phenotypes: cigarettes per day, pack years, smoking initiation, age of smoking initiation, and alcoholic drinks per week. We conducted single variant tests of all variants, and gene-based burden tests of nonsynonymous or putative loss of function variants with minor allele frequency less than 1%.ResultsMeta-analytic sample sizes ranged from 70,847 to 164,142 individuals, depending on the phenotype. Known loci tagged by common variants replicated, but there was no robust evidence for individually associated rare variants, either in gene based or single variant tests. Using a modified method-of-moment approach, we found that all low frequency coding variants, in aggregate, contributed 1.7% to 3.6% of the phenotypic variation for the five traits (pConclusionsThe findings indicate that rare coding variants contribute to phenotypic variation, but that much larger samples and/or denser genotyping of rare variants will be required to successfully identify associations with these phenotypes, whether individual variants or gene‐ based associations.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV1), forced vital capacity (FVC) and the ratio of FEV1 to FVC (FEV1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P-7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs (SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

The role that vitamin D plays in pulmonary function remains uncertain. Epidemiological studies reported mixed findings for the association of serum 25-hydroxyvitamin D [25(OH)D] and pulmonary function. We conducted the largest cross-sectional meta-analysis of the 25(OH)D– pulmonary function association to date, based on nine European ancestry (EA) cohorts (n=22,838) and five African ancestry (AA) cohorts (n=4,290) in the CHARGE Consortium. Data were analyzed using linear models by cohort and ancestry. Effect modification by smoking status (current/former/never) was tested. Results were combined using fixed-effects meta-analysis. Mean (SD) serum 25(OH)D was 68 (29) nmol/L for EAs and 49 (21) nmol/L for AAs. For each 1 nmol/L higher 25(OH)D, forced expiratory volume in the first second (FEV1) was higher by 1.1 mL in EAs (95% CI: 0.9,1.3; P=2.5×10-21) and 1.8 mL (95% CI: 1.1,2.5; P=1.6×10-7) in Aas (Prace difference=0.06), and forced vital capacity (FVC) was higher by 1.3 mL in EAs (95% CI: 1.0,1.6; P=1.1×10-20) and 1.5 mL (95% CI: 0.8,2.3; P=1.2×10-4) in AAs (Prace difference=0.56). Among EAs, the 25(OH)D–FVC association was stronger in smokers: per 1nmol/L higher 25(OH)D, FVC was higher by 1.7 mL (95% CI: 1.1,2.3) for current smokers and 1.7 mL (95% CI: 1.2,2.1) for former smokers, compared to 0.8 mL (95% CI: 0.4,1.2) for never smokers. In summary, the 25(OH)D associations with FEV1 and FVC were positive in both ancestries. In EAs, a stronger association was observed for smokers compared to never smokers, which supports the importance of vitamin D in vulnerable populations.Cohort FundingThis work was supported by National Institutes of Health (NIH) grant number R21 HL125574 funded by the National Heart, Lung, and Blood Institute (NHLBI) and the NIH Office of Dietary Supplements (ODS) (co-Principal Investigators [co-PIs]: DBH and PAC). The corresponding author (PAC) had full access to the data for the meta-analysis, and had final responsibility for the decision to submit for publication. No funding source had any role in the analysis of the data, the writing of the manuscript, or the decision to submit it. This work was also supported in part by R01HL077612 (PI: RGB) and by the Intramural Research Program of the National Institutes of Health (NIH), National Institute of Environmental Health Sciences (ZO1 ES043012, PI: SJL). SJL is supported by the Intramural Research Program of NIH, National Institute of Environmental Health Sciences. Infrastructure for the CHARGE Consortium is supported in part by the NHLBI grant R01HL105756.The Age, Gene/Environment Susceptibility (AGES)–Reykjavik Study has been funded by NIH contracts N01-AG-1-2100 and 271201200022C, the National Institute on Aging (NIA) Intramural Research Program, Hjartavernd (the Icelandic Heart Association), and the Althingi (the Icelandic Parliament). The study is approved by the Icelandic National Bioethics Committee, VSN: 00-063. The researchers are indebted to the participants for their willingness to participate in the study.The Atherosclerosis Risk in Communities Study is carried out as a collaborative study supported by NHLBI contracts HHSN268201100005C, HHSN268201100006C, HHSN268201100007C, HHSN268201100008C, HHSN268201100009C, HHSN268201100010C, HHSN268201100011C, and HHSN268201100012C. 25(OH)D measurements were conducted with the support of R01 HL103706 from the NHLBI and R01 HL103706-S1 from the NIH ODS. The authors thank the staff and participants of the ARIC study for their important contributions.This Cardiovascular Health Study (CHS) research was supported by NHLBI contracts HHSN268201200036C, HHSN268200800007C, N01HC55222, N01HC85079, N01HC85080, N01HC85081, N01HC85082, N01HC85083, N01HC85086; and NHLBI grants U01HL080295, R01HL085251, R01HL087652, R01HL105756, R01HL103612, R01HL120393, and R01HL130114 with additional contribution from the National Institute of Neurological Disorders and Stroke (NINDS). Additional support was provided through R01AG023629 from NIA. A full list of principal CHS investigators and institutions can be found at CHS-NHLBI.org. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Vitamin D measurements were made possible by NHLBI (R01HL084443-01A2).This work in Framingham Heart Study was supported by NHLBI’s Framingham Heart Study contract (N01-HC-25195 and HHSN268201500001I). Vitamin D measurements in the Framingham study were made possible by NIA (R01 AG14759 to SLB.).The Health Aging and Body Composition cohort study was supported by NIA contracts N01AG62101, N01AG2103, and N01AG62106, NIA grant R01-AG028050, NINR grant R01-NR012459, and in part by the Intramural Research Program of the NIA, NIH. This research was further supported by RC1AG035835, and the serum vitamin D assays were supported by R01AG029364.The Multi-Ethnic Study of Atherosclerosis (MESA) study is conducted and supported by NHLBI in collaboration with MESA investigators. Support for MESA is provided by contracts HHSN268201500003I, N01-HC-95159, N01-HC-95160, N01-HC-95161, N01-HC-95162, N01-HC-95163, N01-HC-95164, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, and N01-HC-95169 from NHLBI, UL1-TR-000040, UL1-TR-001079, and UL1-TR-001881 from NCRR, and DK063491 from the NIDDK. The MESA Lung study was supported by grants R01 HL077612, RC1 HL100543 and R01 HL093081 from NHLBI. Support for the Mineral Metabolite dataset was provided by grant HL096875.The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, the Netherlands; the Organization for the Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly (RIDE); the Dutch Ministry of Education, Culture, and Science; the Dutch Ministry for Health, Welfare, and Sports; the European Commission (DG XII), and the Municipality of Rotterdam. LL was a postdoctoral fellow of the Research Foundation—Flanders (FWO) in Brussels, Belgium. Part of this work was supported by a FWO-grant G035014N. DSM Nutritional Products AG, Kaiseraugst, Switzerland, sponsored the Vitamin D serum analyses. The authors are grateful to the study participants, the staff from the Rotterdam Study, and the participating general practitioners and pharmacists.The Coronary Artery Risk Development in Young Adults Study (CARDIA) is supported by contracts HHSN268201300025C, HHSN268201300026C, HHSN268201300027C, HHSN268201300028C, HHSN268201300029C, and HHSN268200900041C from the National Heart, Lung, and Blood Institute (NHLBI), the Intramural Research Program of the National Institute on Aging (NIA), and an intra-agency agreement between NIA and NHLBI (AG0005).Author DisclosureDr. Psaty serves on the DSMB of a clinical trial funded by the manufacturer (Zoll LifeCor) and on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson.All other authors have no conflicts of interest. There is no commercial support or financial interest from the tobacco industry for the research presented.The study sponsors were not involved in study design, data collection, data analysis, data interpretation, report writing, or decisions to submit the paper for publication. PAC and DBH had final responsibility for the decision to submit for publication.Online Supporting MaterialSupplemental table, figures, and methods are available.Abbreviation Footnote25(OH)D25-Hydroxyvitamin DAAAfrican AncestryAGESAge, Gene, Environment, Susceptibility Study—Reykjavik, IcelandARICAtherosclerosis Risk in Communities StudyCARDIACoronary Artery Risk Development in Young Adults StudyCHARGECohorts for Heart and Aging Research in Genomic Epidemiology ConsortiumCHSCardiovascular Health StudyCLIAChemiluminescence ImmunoassayCOPDChronic Obstructive Pulmonary DiseaseEAEuropean AncestryFEV1Forced Expiratory Volume in the First SecondFHS (Offspring)Framingham Heart Study—Offspring CohortFHS (Gen3)Framingham Heart Study—Generation 3 CohortFVCForced Vital CapacityHABCHealth, Aging, and Body Composition StudyLC-MS/MSLiquid Chromatography in Tandem with Mass SpectrometryMESAMulti-Ethnic Study of AtherosclerosisNHANESNational Health and Nutrition Examination SurveyPFTPulmonary Function TestRIARadioimmunoassayRSRotterdam (Netherlands) Study