Search

Rationale & objectiveTubulointerstitial damage is a feature of early chronic kidney disease (CKD), but current clinical tests capture it poorly. Urine biomarkers of tubulointerstitial health may identify risk of CKD.Study designProspective cohort (Atherosclerosis Risk in Communities [ARIC]) and case-cohort (Multi-Ethnic Study of Atherosclerosis [MESA] and Reasons for Geographic and Racial Differences in Stroke [REGARDS]).Setting & participantsAdults with estimated glomerular filtration rate (eGFR) ≥60 mL/min/1.73 m2 and without diabetes in the ARIC, REGARDS, and MESA studies.ExposuresBaseline urine monocyte chemoattractant protein-1 (MCP-1), alpha-1-microglobulin (α1m), kidney injury molecule-1, epidermal growth factor, and chitinase-3-like protein 1.OutcomeIncident CKD or end-stage kidney disease.Analytical approachMultivariable Cox proportional hazards regression for each cohort; meta-analysis of results from all 3 cohorts.Results872 ARIC participants (444 cases of incident CKD), 636 MESA participants (158 cases), and 924 REGARDS participants (488 cases) were sampled. Across cohorts, mean age ranged from 60 ± 10 to 63 ± 8 years, and baseline eGFR ranged from 88 ± 13 to 91 ± 14 mL/min/1.73 m2. In ARIC, higher concentrations of urine MCP-1, α1m, and kidney injury molecule-1 were associated with incident CKD. In MESA, higher concentration of urine MCP-1 and lower concentration of epidermal growth factor were each associated with incident CKD. In REGARDS, none of the biomarkers were associated with incident CKD. In meta-analysis of all 3 cohorts, each 2-fold increase α1m concentration was associated with incident CKD (HR, 1.19; 95% CI, 1.08-1.31).LimitationsObservational design susceptible to confounding; competing risks during long follow-up period; meta-analysis limited to 3 cohorts.ConclusionsIn 3 combined cohorts of adults without prevalent CKD or diabetes, higher urine α1m concentration was independently associated with incident CKD. 4 biomarkers were associated with incident CKD in at least 1 of the cohorts when analyzed individually. Kidney tubule health markers might inform CKD risk independent of eGFR and albuminuria.

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P

Markers of glomerular disease, estimated glomerular filtration rate (eGFR) and albuminuria, are associated with cardiac structural abnormalities and incident cardiovascular disease (CVD). We aimed to determine whether biomarkers of kidney tubule injury, function, and systemic inflammation are associated with cardiac structural abnormalities. Among 393 Multi-Ethnic Study of Atherosclerosis participants without diabetes, CVD, or chronic kidney disease, we assessed the association of 12 biomarkers of kidney tubule injury, function, and systemic inflammation with the left ventricular mass/volume ratio (LVmvr) and left ventricular ejection fraction (LVEF) on cardiac magnetic resonance imaging using linear regression. The average age was 60 ± 10 years; 48% were men; mean eGFR was 96±16 ml/min/1.73 m2; mean LVmvr was 0.93±0.18 g/ml, and mean LVEF was 62±6%. Each twofold greater concentration of plasma soluble urokinase plasminogen activator receptor was associated with a 0.04 g/ml (95% confidence interval [CI] 0.01 to 0.08 g/ml) higher LVmvr and 2.1% (95% CI 0.6 to 3.5%) lower LVEF, independent of risk factors for CVD, eGFR, and albuminuria. Each twofold greater plasma monocyte chemoattractant protein 1 was associated with higher LVmvr with a similar coefficient to that of plasma soluble urokinase plasminogen activator receptor. Each twofold greater concentration of plasma chitinase-3-like protein 1 and urine alpha-1-microglobulin was associated with a 1.1% (95% CI 0.4 to 1.7%) and 1.2% (95% CI 0.2 to 2.2%) lower LVEF, respectively. In conclusion, abnormal kidney tubule health may lead to cardiac dysfunction above and beyond eGFR and albuminuria.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

BackgroundAlzheimer's disease and related dementias (ADRD) involve biological processes that begin years to decades before onset of clinical symptoms. The plasma proteome can offer insight into brain aging and risk of incident dementia among cognitively healthy adults.ObjectiveTo identify biomarkers and biological pathways associated with neuroimaging measures and incident dementia in two large community-based cohorts by applying a correlation-based network analysis to the plasma proteome.MethodsWeighted co-expression network analysis of 1,305 plasma proteins identified four modules of co-expressed proteins, which were related to MRI brain volumes and risk of incident dementia over a median 20-year follow-up in Framingham Heart Study (FHS) Offspring cohort participants (n = 1,861). Analyses were replicated in the Cardiovascular Health Study (CHS) (n = 2,117, mean 6-year follow-up).ResultsTwo proteomic modules, one related to protein clearance and synaptic maintenance (M2) and a second to inflammation (M4), were associated with total brain volume in FHS (M2: p = 0.014; M4: p = 4.2×10-5). These modules were not significantly associated with hippocampal volume, white matter hyperintensities, or incident all-cause or AD dementia. Associations with TCBV did not replicate in CHS, an older cohort with a greater burden of comorbidities.ConclusionsProteome networks implicate an early role for biological pathways involving inflammation and synaptic function in preclinical brain atrophy, with implications for clinical dementia.

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p

A bidirectional communication exists between the brain and the gut, in which the gut microbiota influences cognitive function and vice-versa. Gut dysbiosis has been linked to several diseases, including Alzheimer's disease and related dementias (ADRD). However, the relationship between gut dysbiosis and markers of cerebral small vessel disease (cSVD), a major contributor to ADRD, is unknown. In this cross-sectional study, we examined the connection between the gut microbiome, cognitive, and neuroimaging markers of cSVD in the Framingham Heart Study (FHS). Markers of cSVD included white matter hyperintensities (WMH), peak width of skeletonized mean diffusivity (PSMD), and executive function (EF), estimated as the difference between the trail-making tests B and A. We included 972 FHS participants with MRI scans, neurocognitive measures, and stool samples and quantified the gut microbiota composition using 16S rRNA sequencing. We used multivariable association and differential abundance analyses adjusting for age, sex, BMI, and education level to estimate the association between gut microbiota and WMH, PSMD, and EF measures. Our results suggest an increased abundance of Pseudobutyrivibrio and Ruminococcus genera was associated with lower WMH and PSMD (p values

We develop a closed-form Haseman-Elston estimator for genetic and environmental correlation coefficients between complex phenotypes, which we term HEc, that is as precise as GCTA yet ∼20× faster. We estimate genetic and environmental correlations between over 7,000 phenotype pairs in subgroups from the Trans-Omics in Precision Medicine (TOPMed) program. We demonstrate substantial differences in both heritabilities and genetic correlations for multiple phenotypes and phenotype pairs between individuals of self-reported Black, Hispanic/Latino, and White backgrounds. We similarly observe differences in many of the genetic and environmental correlations between genders. To estimate the contribution of genetics to the observed phenotypic correlation, we introduce "fractional genetic correlation" as the fraction of phenotypic correlation explained by genetics. Finally, we quantify the enrichment of correlations between phenotypic domains, each of which is comprised of multiple phenotypes. Altogether, we demonstrate that the observed correlations between complex human phenotypes depend on the genetic background of the individuals, their gender, and their environment.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

BackgroundWhile large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries.ResultsAssociations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10- 7).ConclusionsLow frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

Rationale: Obstructive sleep apnea (OSA) is a common disorder associated with increased risk for cardiovascular disease, diabetes, and premature mortality. There is strong clinical and epidemiologic evidence supporting the importance of genetic factors influencing OSA but limited data implicating specific genes. Objectives: To search for rare variants contributing to OSA severity. Methods: Leveraging high-depth genomic sequencing data from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and imputed genotype data from multiple population-based studies, we performed linkage analysis in the CFS (Cleveland Family Study), followed by multistage gene-based association analyses in independent cohorts for apnea-hypopnea index (AHI) in a total of 7,708 individuals of European ancestry. Measurements and Main Results: Linkage analysis in the CFS identified a suggestive linkage peak on chromosome 7q31 (LOD = 2.31). Gene-based analysis identified 21 noncoding rare variants in CAV1 (Caveolin-1) associated with lower AHI after accounting for multiple comparisons (P = 7.4 × 10-8). These noncoding variants together significantly contributed to the linkage evidence (P

RationaleGenetic variation has a substantial contribution to chronic obstructive pulmonary disease (COPD) and lung function measurements. Heritability estimates using genome-wide genotyping data can be biased if analyses do not appropriately account for the nonuniform distribution of genetic effects across the allele frequency and linkage disequilibrium (LD) spectrum. In addition, the contribution of rare variants has been unclear.ObjectivesWe sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.MethodsUsing the genome-based restricted maximum likelihood method, we partitioned the genome into bins based on minor allele frequency and LD scores and estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio in 11 051 European ancestry and 5853 African-American participants.Measurements and main resultsIn European ancestry participants, the estimated heritability of COPD, FEV1% predicted and FEV1/FVC ratio were 35.5%, 55.6% and 32.5%, of which 18.8%, 19.7%, 17.8% were from common variants, and 16.6%, 35.8%, and 14.6% were from rare variants. These estimates had wide confidence intervals, with common variants and some sets of rare variants showing a statistically significant contribution (P-value

Common genetic variants explain less variation in complex phenotypes than inferred from family-based studies, and there is a debate on the source of this 'missing heritability'. We investigated the contribution of rare genetic variants to tobacco use with whole-genome sequences from up to 26,257 unrelated individuals of European ancestries and 11,743 individuals of African ancestries. Across four smoking traits, single-nucleotide-polymorphism-based heritability ([Formula: see text]) was estimated from 0.13 to 0.28 (s.e., 0.10-0.13) in European ancestries, with 35-74% of it attributable to rare variants with minor allele frequencies between 0.01% and 1%. These heritability estimates are 1.5-4 times higher than past estimates based on common variants alone and accounted for 60% to 100% of our pedigree-based estimates of narrow-sense heritability ([Formula: see text], 0.18-0.34). In the African ancestry samples, [Formula: see text] was estimated from 0.03 to 0.33 (s.e., 0.09-0.14) across the four smoking traits. These results suggest that rare variants are important contributors to the heritability of smoking.

BackgroundEpigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD.MethodsBlood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE-/-) mouse model of atherosclerosis.ResultsA total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic.ConclusionsDifferential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

The Collaborative Cohort of Cohorts for COVID-19 Research (C4R) is a national prospective study of adults comprising 14 established US prospective cohort studies. Starting as early as 1971, investigators in the C4R cohort studies have collected data on clinical and subclinical diseases and their risk factors, including behavior, cognition, biomarkers, and social determinants of health. C4R links this pre-coronavirus disease 2019 (COVID-19) phenotyping to information on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and acute and postacute COVID-related illness. C4R is largely population-based, has an age range of 18-108 years, and reflects the racial, ethnic, socioeconomic, and geographic diversity of the United States. C4R ascertains SARS-CoV-2 infection and COVID-19 illness using standardized questionnaires, ascertainment of COVID-related hospitalizations and deaths, and a SARS-CoV-2 serosurvey conducted via dried blood spots. Master protocols leverage existing robust retention rates for telephone and in-person examinations and high-quality event surveillance. Extensive prepandemic data minimize referral, survival, and recall bias. Data are harmonized with research-quality phenotyping unmatched by clinical and survey-based studies; these data will be pooled and shared widely to expedite collaboration and scientific findings. This resource will allow evaluation of risk and resilience factors for COVID-19 severity and outcomes, including postacute sequelae, and assessment of the social and behavioral impact of the pandemic on long-term health trajectories.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p

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 the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

ObjectiveTrans fatty acids (TFAs) have harmful biologic effects that could increase the risk of type 2 diabetes (T2D), but evidence remains uncertain. We aimed to investigate the prospective associations of TFA biomarkers and T2D by conducting an individual participant-level pooled analysis.Research design and methodsWe included data from an international consortium of 12 prospective cohorts and nested case-control studies from six nations. TFA biomarkers were measured in blood collected between 1990 and 2008 from 25,126 participants aged ≥18 years without prevalent diabetes. Each cohort conducted de novo harmonized analyses using a prespecified protocol, and findings were pooled using inverse-variance weighted meta-analysis. Heterogeneity was explored by prespecified between-study and within-study characteristics.ResultsDuring a mean follow-up of 13.5 years, 2,843 cases of incident T2D were identified. In multivariable-adjusted pooled analyses, no significant associations with T2D were identified for trans/trans-18:2, relative risk (RR) 1.09 (95% CI 0.94-1.25); cis/trans-18:2, 0.89 (0.73-1.07); and trans/cis-18:2, 0.87 (0.73-1.03). Trans-16:1n-9, total trans-18:1, and total trans-18:2 were inversely associated with T2D (RR 0.81 [95% CI 0.67-0.99], 0.86 [0.75-0.99], and 0.84 [0.74-0.96], respectively). Findings were not significantly different according to prespecified sources of potential heterogeneity (each P ≥ 0.1).ConclusionsCirculating individual trans-18:2 TFA biomarkers were not associated with risk of T2D, while trans-16:1n-9, total trans-18:1, and total trans-18:2 were inversely associated. Findings may reflect the influence of mixed TFA sources (industrial vs. natural ruminant), a general decline in TFA exposure due to policy changes during this period, or the relatively limited range of TFA levels.

Metabolomics genome wide association study (GWAS) help outline the genetic contribution to human metabolism. However, studies to date have focused on relatively healthy, population-based samples of White individuals. Here, we conducted a GWAS of 537 blood metabolites measured in the Chronic Renal Insufficiency Cohort (CRIC) Study, with separate analyses in 822 White and 687 Black study participants. Trans-ethnic meta-analysis was then applied to improve fine-mapping of potential causal variants. Mean estimated glomerular filtration rate was 44.4 and 41.5 mL/min/1.73m2 in the White and Black participants, respectively. There were 45 significant metabolite associations at 19 loci, including novel associations at PYROXD2, PHYHD1, FADS1-3, ACOT2, MYRF, FAAH, and LIPC. The strength of associations was unchanged in models additionally adjusted for estimated glomerular filtration rate and proteinuria, consistent with a direct biochemical effect of gene products on associated metabolites. At several loci, trans-ethnic meta-analysis, which leverages differences in linkage disequilibrium across populations, reduced the number and/or genomic interval spanned by potentially causal single nucleotide polymorphisms compared to fine-mapping in the White participant cohort alone. Across all validated associations, we found strong concordance in effect sizes of the potentially causal single nucleotide polymorphisms between White and Black study participants. Thus, our study identifies novel genetic determinants of blood metabolites in chronic kidney disease, demonstrates the value of diverse cohorts to improve causal inference in metabolomics GWAS, and underscores the shared genetic basis of metabolism across race.

Background and purposeStroke is the leading cause of death and long-term disability worldwide. Previous genome-wide association studies identified 51 loci associated with stroke (mostly ischemic) and its subtypes among predominantly European populations. Using whole-genome sequencing in ancestrally diverse populations from the Trans-Omics for Precision Medicine (TOPMed) Program, we aimed to identify novel variants, especially low-frequency or ancestry-specific variants, associated with all stroke, ischemic stroke and its subtypes (large artery, cardioembolic, and small vessel), and hemorrhagic stroke and its subtypes (intracerebral and subarachnoid).MethodsWhole-genome sequencing data were available for 6833 stroke cases and 27 116 controls, including 22 315 European, 7877 Black, 2616 Hispanic/Latino, 850 Asian, 54 Native American, and 237 other ancestry participants. In TOPMed, we performed single variant association analysis examining 40 million common variants and aggregated association analysis focusing on rare variants. We also combined TOPMed European populations with over 28 000 additional European participants from the UK BioBank genome-wide array data through meta-analysis.ResultsIn the single variant association analysis in TOPMed, we identified one novel locus 13q33 for large artery at whole-genome-wide significance (P

BackgroundGuidelines emphasize screening people with elevated BP for the presence of end-organ damage.MethodsWe characterized the prevalence, correlates, and prognosis of hypertension-mediated organ damage (HMOD) in the community-based Framingham Study. 7898 participants (mean age 51.6 years, 54% women) underwent assessment for the following HMOD: electrocardiographic and echocardiographic left ventricular hypertrophy, abnormal brain imaging findings consistent with vascular injury, increased carotid intima-media thickness, elevated carotid-femoral pulse wave velocity, reduced kidney function, microalbuminuria, and low ankle-brachial index. We characterized HMOD prevalence according to blood pressure (BP) categories defined by four international BP guidelines. Participants were followed up for incidence of cardiovascular disease.ResultsThe prevalence of HMOD varied positively with systolic BP and pulse pressure but negatively with diastolic BP; it increased with age, was similar in both sexes, and varied across BP guidelines based on their thresholds defining hypertension. Among participants with hypertension, elevated carotid-femoral pulse wave velocity was the most prevalent HMOD (40%-60%), whereas low ankle-brachial index was the least prevalent (

BackgroundObesity and cardiometabolic dysfunction have been associated with cancer risk and severity. Underlying mechanisms remain unclear.ObjectivesThe aim of this study was to examine associations of obesity and related cardiometabolic traits with incident cancer.MethodsFHS (Framingham Heart Study) and PREVEND (Prevention of Renal and Vascular End-Stage Disease) study participants without prevalent cancer were studied, examining associations of obesity, body mass index (BMI), waist circumference, visceral adipose tissue (VAT) and subcutaneous adipose tissue depots, and C-reactive protein (CRP) with future cancer in Cox models.ResultsAmong 20,667 participants (mean age 50 years, 53% women), 2,619 cancer events were observed over a median follow-up duration of 15 years. Obesity was associated with increased risk for future gastrointestinal (HR: 1.30; 95% CI: 1.05-1.60), gynecologic (HR: 1.62; 95% CI: 1.08-2.45), and breast (HR: 1.32; 95% CI: 1.05-1.66) cancer and lower risk for lung cancer (HR: 0.62; 95% CI: 0.44-0.87). Similarly, waist circumference was associated with increased risk for overall, gastrointestinal, and gynecologic but not lung cancer. VAT but not subcutaneous adipose tissue was associated with risk for overall cancer (HR: 1.22; 95% CI: 1.05-1.43), lung cancer (HR: 1.92; 95% CI: 1.01-3.66), and melanoma (HR: 1.56; 95% CI: 1.02-2.38) independent of BMI. Last, higher CRP levels were associated with higher risk for overall, colorectal, and lung cancer (P < 0.05 for all).ConclusionsObesity and abdominal adiposity are associated with future risk for specific cancers (eg, gastrointestinal, gynecologic). Although obesity was associated with lower risk for lung cancer, greater VAT and CRP were associated with higher lung cancer risk after adjusting for BMI.

Analyses of data from genome-wide association studies on unrelated individuals have shown that, for human traits and diseases, approximately one-third to two-thirds of heritability is captured by common SNPs. However, it is not known whether the remaining heritability is due to the imperfect tagging of causal variants by common SNPs, in particular whether the causal variants are rare, or whether it is overestimated due to bias in inference from pedigree data. Here we estimated heritability for height and body mass index (BMI) from whole-genome sequence data on 25,465 unrelated individuals of European ancestry. The estimated heritability was 0.68 (standard error 0.10) for height and 0.30 (standard error 0.10) for body mass index. Low minor allele frequency variants in low linkage disequilibrium (LD) with neighboring variants were enriched for heritability, to a greater extent for protein-altering variants, consistent with negative selection. Our results imply that rare variants, in particular those in regions of low linkage disequilibrium, are a major source of the still missing heritability of complex traits and disease.

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.

Rationale & objectivePlasma kidney injury molecule 1 (KIM-1) is a sensitive marker of proximal tubule injury, but its association with risks of adverse clinical outcomes across a spectrum of kidney diseases is unknown.Study designProspective, observational cohort study.Setting & participants524 individuals enrolled into the Boston Kidney Biopsy Cohort (BKBC) Study undergoing clinically indicated native kidney biopsy with biopsy specimens adjudicated for semiquantitative scores of histopathology by 2 kidney pathologists and 3,800 individuals with common forms of chronic kidney disease (CKD) enrolled into the Chronic Renal Insufficiency Cohort (CRIC) Study.ExposureHistopathologic lesions and clinicopathologic diagnosis in cross-sectional analyses, baseline plasma KIM-1 levels in prospective analyses.OutcomesBaseline plasma KIM-1 levels in cross-sectional analyses, kidney failure (defined as initiation of kidney replacement therapy) and death in prospective analyses.Analytical approachMultivariable-adjusted linear regression models tested associations of plasma KIM-1 levels with histopathologic lesions and clinicopathologic diagnoses. Cox proportional hazards models tested associations of plasma KIM-1 levels with future kidney failure and death.ResultsIn the BKBC Study, higher plasma KIM-1 levels were associated with more severe acute tubular injury, tubulointerstitial inflammation, and more severe mesangial expansion after multivariable adjustment. Participants with diabetic nephropathy, glomerulopathies, and tubulointerstitial disease had significantly higher plasma KIM-1 levels after multivariable adjustment. In the BKBC Study, CKD in 124 participants progressed to kidney failure and 85 participants died during a median follow-up time of 5 years. In the CRIC Study, CKD in 1,153 participants progressed to kidney failure and 1,356 participants died during a median follow-up time of 11.5 years. In both cohorts, each doubling of plasma KIM-1 level was associated with an increased risk of kidney failure after multivariable adjustment (hazard ratios of 1.19 [95% CI, 1.03-1.38] and 1.10 [95% CI, 1.06-1.15] for BKBC and CRIC, respectively). There was no statistically significant association of plasma KIM-1 levels with death in either cohort.LimitationsGeneralizability and unmeasured confounding.ConclusionsPlasma KIM-1 is associated with underlying tubulointerstitial and mesangial lesions and progression to kidney failure in 2 cohort studies of individuals with kidney diseases.

Blood lipids are heritable modifiable causal factors for coronary artery disease. Despite well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing (WGS), partly due to limited sample sizes, ancestral diversity, and interpretation of clinical significance. Among 66,329 ancestrally diverse (56% non-European) participants, we associate 428M variants from deep-coverage WGS with lipid levels; ~400M variants were not assessed in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with blood lipids through analysis of common and rare coding variants. We discover several associated rare non-coding variants, largely at Mendelian lipid genes. Notably, we observe rare LDLR intronic variants associated with markedly increased LDL-C, similar to rare LDLR exonic variants. In conclusion, we conducted a systematic whole genome scan for blood lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.

Deficiency of the immune checkpoint lymphocyte activation gene-3 (LAG3) protein is significantly associated with both elevated HDL-cholesterol (HDL-C) and myocardial infarction risk. We determined the association of genetic variants within ±500 kb of LAG3 with plasma LAG3 and defined LAG3-associated plasma proteins with HDL-C and clinical outcomes. Whole genome sequencing and plasma proteomics were obtained from the Multi-Ethnic Study of Atherosclerosis (MESA) and the Framingham Heart Study (FHS) cohorts as part of the Trans-Omics for Precision Medicine program. In situ Hi-C chromatin capture was performed in EBV-transformed cell lines isolated from four MESA participants. Genetic association analyses were performed in MESA using multivariate regression models, with validation in FHS. A LAG3-associated protein network was tested for association with HDL-C, coronary heart disease, and all-cause mortality. We identify an association between the LAG3 rs3782735 variant and plasma LAG3 protein. Proteomics analysis reveals 183 proteins significantly associated with LAG3 with four proteins associated with HDL-C. Four proteins discovered for association with all-cause mortality in FHS shows nominal associations in MESA. Chromatin capture analysis reveals significant cis interactions between LAG3 and C1S, LRIG3, TNFRSF1A, and trans interactions between LAG3 and B2M. A LAG3-associated protein network has significant associations with HDL-C and mortality.

BackgroundDysregulation of compensatory mechanisms to regulate blood pressure (BP) upon postural change is a phenotype of BP variability and an emerging risk factor for cardiovascular outcomes.Materials and methodsWe assessed postural change in BP (starting 2 min after standing from a supine position), carotid-femoral pulse wave velocity (cfPWV), and markers of hypertension-mediated organ damage (HMOD) in the heart, kidney, and brain in Framingham Third Generation, Omni-2, and New Offspring Spouse Cohort participants. We related vascular measures (postural change in BP measures and cfPWV) with HMOD in 3,495 participants (mean age 47 years, 53% women) using multivariable logistic and linear regression models.ResultsIn multivariable-adjusted models, we did not observe significant associations of vascular measures with presence of left ventricular hypertrophy, albuminuria, covert brain infarcts, or white matter hyperintensities (Bonferroni-adjusted P-values > 0.05/20 > 0.0025). In multivariable models, greater cfPWV (est. β = 0.11 ± 0.03; P < 0.001), but not postural change in BP measures (Bonferroni-adjusted P-values > 0.05/20 > 0.0025), was associated with higher white matter free water using brain magnetic resonance imaging. In multivariable models, greater postural change in pulse pressure was associated with higher urinary albumin-creatinine ratio (est. β = 0.07 ± 0.02; P < 0.001). No other postural change in BP measure was associated with urinary albumin-creatinine ratio (Bonferroni-adjusted P-values > 0.05/20 > 0.0025). In sex-specific analyses, higher cfPWV was associated with higher urinary albumin-creatinine ratio in men (est. β: 0.11 ± 0.04; P = 0.002) but not in women (est. β: 0.03 ± 0.03; P = 0.44). We also observed marginal to strong effect modification by above vs. at/below median postural change in BP for the association of cfPWV with urinary albumin-creatinine ratio (Bonferroni-adjusted interaction P < 0.001-0.01). Vascular measures were not related to left ventricular mass index or fractional anisotropy (Bonferroni-adjusted P-values > 0.05/20 > 0.0025).ConclusionBaroreflex dysfunction is associated with greater subclinical kidney damage. Additionally, relations of higher aortic stiffness with greater kidney damage may be modified by associated baroreflex dysregulation.

BackgroundPrevious studies suggest poor pulmonary function is associated with increased burden of cerebral white matter hyperintensities and brain atrophy among elderly individuals, but the results are inconsistent.ObjectiveTo study the cross-sectional associations of pulmonary function with structural brain variables.MethodsData from six large community-based samples (N = 11,091) were analyzed. Spirometric measurements were standardized with respect to age, sex, height, and ethnicity using reference equations of the Global Lung Function Initiative. Associations of forced expiratory volume in 1 second (FEV1), forced vital capacity (FVC), and their ratio FEV1/FVC with brain volume, gray matter volume, hippocampal volume, and volume of white matter hyperintensities were investigated using multivariable linear regressions for each study separately and then combined using random-effect meta-analyses.ResultsFEV1 and FVC were positively associated with brain volume, gray matter volume, and hippocampal volume, and negatively associated with white matter hyperintensities volume after multiple testing correction, with little heterogeneity present between the studies. For instance, an increase of FVC by one unit was associated with 3.5 ml higher brain volume (95% CI: [2.2, 4.9]). In contrast, results for FEV1/FVC were more heterogeneous across studies, with significant positive associations with brain volume, gray matter volume, and hippocampal volume, but not white matter hyperintensities volume. Associations of brain variables with both FEV1 and FVC were consistently stronger than with FEV1/FVC, specifically with brain volume and white matter hyperintensities volume.ConclusionIn cross-sectional analyses, worse pulmonary function is associated with smaller brain volumes and higher white matter hyperintensities burden.

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value

BackgroundInsulin-like growth factor 1 (IGF-1) signaling has been implicated in Alzheimer's disease pathogenesis, and further evidence suggests inflammation can be a moderator of this association. However, most research to date has been conducted on older adults.ObjectiveTo investigate the association of serum IGF-1 and IGF binding protein 3 (IGFBP-3) concentrations with MRI markers of Alzheimer's disease in predominantly middle-aged adults, and further assess moderation by chronic inflammation.MethodsWe included participants from the Framingham Heart Study (n = 1,852, mean age 46±8, 46% men) and the Study of Health in Pomerania (n = 674, mean age 50±13, 42% men) with available serum IGF-1, IFGBP-3, as well as brain MRI. IGF-1 and IFGBP-3 were related to MRI outcomes (i.e., total brain, cortical gray matter, white matter, white matter hyperintensities (WMH), and hippocampal volumes) using multivariable regression models adjusting for potential confounders. Subgroup analyses by C-reactive protein (CRP) concentrations were also performed. Cohort-specific summary statistics were meta-analyzed using random-effects models and corrected for multiple comparisons.ResultsMeta-analysis results revealed that higher IGF-1 concentrations were associated with lower WMH (estimate [β] [95% CI], -0.05 [-0.09, -0.02], p = 0.006) and larger hippocampal volumes (0.07 [0.02, 0.12], p = 0.01), independent of vascular risk factors. These associations occurred predominantly in individuals with CRP concentrations

Circulating total-tau levels can be used as an endophenotype to identify genetic risk factors for tauopathies and related neurological disorders. Here, we confirmed and better characterized the association of the 17q21 MAPT locus with circulating total-tau in 14,721 European participants and identified three novel loci in 953 African American participants (4q31, 5p13, and 6q25) at P

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome arrays, resulting in over 100 associated variants. We extended this work with high-coverage whole genome sequencing analyses from fifteen cohorts in NHLBI's Trans-Omics for Precision Medicine (TOPMed) program. Over 23,000 non-diabetic individuals from five race-ethnicities/populations (African, Asian, European, Hispanic and Samoan) were included. Eight variants were significantly associated with FG or FI across previously identified regions MTNR1B, G6PC2, GCK, GCKR and FOXA2. We additionally characterize suggestive associations with FG or FI near previously identified SLC30A8, TCF7L2, and ADCY5 regions as well as APOB, PTPRT, and ROBO1. Functional annotation resources including the Diabetes Epigenome Atlas were compiled for each signal (chromatin states, annotation principal components, and others) to elucidate variant-to-function hypotheses. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions creating a foundation for future sequencing-based investigations of glycemic traits.

Although immunoassays are the most widely used protein measurement method, aptamer-based methods such as the SomaScan platform can quantify up to 7000 proteins per biosample, creating new opportunities for unbiased discovery. However, there is limited research comparing the consistency of biomarker-disease associations between immunoassay and aptamer-based platforms. In a substudy of the TRIBE-AKI cohort, preoperative and postoperative plasma samples from 294 patients with previous immunoassay measurements were analyzed using the SomaScan platform. Inter-platform Spearman correlations (rs) and biomarker-AKI associations were compared across 30 preoperative and 34 postoperative immunoassay-aptamer pairs. Possible factors contributing to inter-platform differences were examined including target protein characteristics, immunoassay, and SomaScan coefficients of variation, other assay characteristics, and sample storage time. The median rs was 0.54 (interquartile range [IQR] 0.34-0.83) in postoperative samples and 0.41 (IQR 0.21-0.69) in preoperative samples. We observed a trend of greater rs in biomarkers with greater concentrations; the Spearman correlation between the concentration of protein and the inter-platform correlation was 0.64 in preoperative pairs and 0.53 in postoperative pairs. Of proteins measured by immunoassays, we observed significant biomarker-AKI associations for 13 proteins preop and 24 postop; of all corresponding aptamers, 8 proteins preop and 12 postop. All proteins significantly associated with AKI as measured by SomaScan were also significantly associated with AKI as measured by immunoassay. All biomarker-AKI odds ratios were significantly different (P < 0.05) between platforms in 14% of aptamer-immunoassay pairs, none of which had high (rs > 0.50) inter-platform correlations. Although similar biomarker-disease associations were observed overall, biomarkers with high physiological concentrations tended to have the highest-confidence inter-platform operability in correlations and biomarker-disease associations. Aptamer assays provide excellent precision and an unprecedented coverage and promise for disease associations but interpretation of results should keep in mind a broad range of correlations with immunoassays.

BackgroundTotal serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE.ObjectiveWe aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901).MethodsWe performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data.ResultsWe identified 6 loci at genome-wide significance (P < 5 × 10-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (Pdiscovery = 2 × 10-4; Preplication = 5 × 10-4; Pdiscovery+replication = 4 × 10-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (PHispanic/Latino discovery+replication = 8 × 10-8).ConclusionWe performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.

BackgroundSleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing.MethodsThe study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (

Genotype-phenotype association studies often combine phenotype data from multiple studies to increase statistical 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 more than 80 studies with extensive phenotype data. To date, 63 phenotypes have been harmonized across thousands of participants (recruited in 1948-2012) from up to 17 studies per phenotype. Here we discuss challenges 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 software code used to generate the 63 harmonized phenotypes, enabling others to reproduce, modify, or extend these harmonizations to additional studies, and 2) the results of labeling thousands of phenotype variables with controlled vocabulary terms.

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.