Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function

BACKGROUND

Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function.

OBJECTIVES

We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions.

METHODS

We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings.

RESULTS

Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers.

CONCLUSION

Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Key variants via the Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases = 2184, N controls = 2383) and targeted analyses in subpopulations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Association of large core middle cerebral artery stroke and hemorrhagic transformation with hospitalization outcomes

Historically, investigators have not differentiated between patients with and without hemorrhagic transformation (HT) in large core ischemic stroke at risk for life-threatening mass effect (LTME) from cerebral edema. Our objective was to determine whether LTME occurs faster in those with HT compared to those without. We conducted a two-center retrospective study of patients with ≥ 1/2 MCA territory infarct between 2006 and 2021. We tested the association of time-to-LTME and HT subtype (parenchymal, petechial) using Cox regression, controlling for age, mean arterial pressure, glucose, tissue plasminogen activator, mechanical thrombectomy, National Institute of Health Stroke Scale, antiplatelets, anticoagulation, temperature, and stroke side. Secondary and exploratory outcomes included mass effect-related death, all-cause death, disposition, and decompressive hemicraniectomy. Of 840 patients, 358 (42.6%) had no HT, 403 (48.0%) patients had petechial HT, and 79 (9.4%) patients had parenchymal HT. LTME occurred in 317 (37.7%) and 100 (11.9%) had mass effect-related deaths. Parenchymal (HR 8.24, 95% CI 5.46-12.42, p < 0.01) and petechial HT (HR 2.47, 95% CI 1.92-3.17, p < 0.01) were significantly associated with time-to-LTME and mass effect-related death. Understanding different risk factors and sequelae of mass effect with and without HT is critical for informed clinical decisions.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals

MOTIVATION

Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus.

RESULTS

We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modelling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 (RSPO3) GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals.

AVAILABILITY

SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes

OBJECTIVE

To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D).

RESEARCH DESIGN AND METHODS

We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD.

RESULTS

Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs11444867 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Later Midline Shift Is Associated with Better Outcomes after Large Middle Cerebral Artery Stroke

Background/Objective

Space occupying cerebral edema is the most feared early complication after large ischemic stroke, occurring in up to 30% of patients with middle cerebral artery (MCA) occlusion, and is reported to peak 2-4 days after injury. Little is known about the factors and outcomes associated with peak edema timing, especially when it occurs after 96 hours. We aimed to characterize differences between patients who experienced maximum midline shift (MLS) or decompressive hemicraniectomy (DHC) in the acute (<48 hours), average (48-96 hours), and subacute (>96 hours) groups and determine whether patients with subacute peak edema timing have improved discharge dispositions.

Methods

We performed a two-center, retrospective study of patients with ≥1/2 MCA territory infarct and MLS. We constructed a multivariable model to test the association of subacute peak edema and favorable discharge disposition, adjusting for age, admission Alberta Stroke Program Early CT Score (ASPECTS), National Institute of Health Stroke Scale (NIHSS), acute thrombolytic intervention, cerebral atrophy, maximum MLS, parenchymal hemorrhagic transformation, DHC, and osmotic therapy receipt.

Results

Of 321 eligible patients with MLS, 32%, 36%, and 32% experienced acute, average, and subacute peak edema. Subacute peak edema was significantly associated with higher odds of favorable discharge than non-subacute swelling, adjusting for confounders (aOR, 1.85; 95% CI, 1.05-3.31).

Conclusions

Subacute peak edema after large MCA stroke is associated with better discharge disposition compared to earlier peak edema courses. Understanding how the timing of cerebral edema affects risk of unfavorable discharge has important implications for treatment decisions and prognostication.

Frequency of Variants in Mendelian Alzheimer's Disease Genes within the Alzheimer's Disease Sequencing Project (ADSP)

Background

Prior studies using the ADSP data examined variants within presenilin-2 ( PSEN2 ), presenilin-1 ( PSEN1 ), and amyloid precursor protein ( APP ) genes. However, previously-reported clinically-relevant variants and other predicted damaging missense (DM) variants have not been characterized in a newer release of the Alzheimer's Disease Sequencing Project (ADSP).

Objective

To characterize previously-reported clinically-relevant variants and DM variants in PSEN2, PSEN1, APP within the participants from the ADSP.

Methods

We identified rare variants (MAF <1%) previously-reported in PSEN2 , PSEN1, and APP in the available ADSP sample of 14,641 individuals with whole genome sequencing and 16,849 individuals with whole exome sequencing available for research-use (N total = 31,490). We additionally curated variants in these three genes from ClinVar, OMIM, and Alzforum and report carriers of variants in clinical databases as well as predicted DM variants in these genes.

Results

We detected 31 previously-reported clinically-relevant variants with alternate alleles observed within the ADSP: 4 variants in PSEN2 , 25 in PSEN1 , and 2 in APP . The overall variant carrier rate for the 31 clinically-relevant variants in the ADSP was 0.3%. We observed that 79.5% of the variant carriers were cases compared to 3.9% were controls. In those with AD, the mean age of onset of AD among carriers of these clinically-relevant variants was 19.6 ± 1.4 years earlier compared with non-carriers (p-value=7.8×10 -57 ).

Conclusion

A small proportion of individuals in the ADSP are carriers of a previously-reported clinically-relevant variant allele for AD and these participants have significantly earlier age of AD onset compared to non-carriers.

Genetically Predicted Body Mass Index and Mortality in COPD

RATIONALE

BMI is associated with COPD mortality, but the underlying mechanisms are unclear. The effect of genetic variants aggregated into a polygenic score may elucidate causal mechanisms and predict risk.

OBJECTIVES

To examine the associations of genetically predicted BMI with all-cause and cause-specific mortality in COPD.

METHODS

We developed a polygenic score for BMI (PGSBMI) and tested for associations of the PGSBMI with all-cause, respiratory, and cardiovascular mortality in participants with COPD from the COPDGene, ECLIPSE, and Framingham Heart studies. We calculated the difference between measured BMI and PGS-predicted BMI (BMIdiff) and categorized participants into groups of discordantly low (BMIdiff < 20th percentile), concordant (BMIdiff between 20th - 80th percentile), and discordantly high (BMIdiff > 80th percentile) BMI. We applied Cox models, examined potential non-linear associations of the PGSBMI and BMIdiff with mortality, and summarized results with meta-analysis.

MEASUREMENTS AND MAIN RESULTS

We observed significant non-linear associations of measured BMI and BMIdiff, but not PGSBMI, with all-cause mortality. In meta-analyses, a one standard deviation increase in the PGSBMI was associated with an increased hazard for cardiovascular mortality (HR=1.29, 95% CI=1.12-1.49), but not with respiratory or all-cause mortality. Compared to participants with concordant measured and genetically predicted BMI, those with discordantly low BMI had higher mortality risk for all-cause (HR=1.57, CI=1.41-1.74) and respiratory death (HR=2.01, CI=1.61-2.51).

CONCLUSIONS

In people with COPD, higher genetically predicted BMI is associated with higher cardiovascular mortality but not respiratory mortality. Individuals with discordantly low BMI have higher all-cause and respiratory mortality compared to those with concordant BMI.

Genome-wide association study of preserved ratio impaired spirometry (PRISm)

BACKGROUND

Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities.

METHODS

We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed.

RESULTS

22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits.

CONCLUSION

This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Follow-up ASPECTS improves prediction of potentially lethal malignant edema in patients with large middle cerebral artery stroke

BACKGROUND

Recent studies have shown that follow-up head CT is a strong predictor of functional outcomes in patients with middle cerebral artery stroke and mechanical thrombectomy. We sought to determine whether total and/or regional follow-up Alberta Stroke Program Early CT Score (ASPECTSfu) are associated with important clinical outcomes during hospitalization and improve the performance of clinical prediction models of potentially lethal malignant edema (PLME).

METHODS

We conducted a retrospective study of patients at three medical centers in a major North American metropolitan area with baseline and follow-up head CTs after large middle cerebral artery stroke between 2006 and 2022. We used multivariable logistic regression to test the association of total and regional ASPECTSfu with PLME (cerebral edema related death or surgery), adjusting for total baseline ASPECTS, age, sex, admission glucose, tissue plasminogen activator, and mechanical thrombectomy. We compared existing clinical risk models with and without total or regional ASPECTSfu using area under the curve.

RESULTS

In our 560 patient cohort, lower total ASPECTSfu was significantly associated with higher odds of PLME when adjusting for confounders (OR 1.69, 95% CI 1.49 to 2.0), and improved model discrimination compared with existing models and models using baseline ASPECTS. Deep territory involvement (OR 2.46, 95% CI 1.53 to 4.01) and anterior territory involvement (OR 3.23, 95% CI 1.88 to 5.71) were significantly associated with PLME.

CONCLUSIONS

Lower ASPECTSfu and certain locations on regional ASPECTSfu, including deep and anterior areas, were significantly associated with PLME. Including ASPECTSfu information improved discrimination of established edema prediction models and could be used immediately to help facilitate clinical management decisions and prognostication.

Bik promotes proteasomal degradation to control low-grade inflammation

Although chronic low-grade inflammation does not cause immediate clinical symptoms, over the longer term, it can enhance other insults or age-dependent damage to organ systems and thereby contribute to age-related disorders, such as respiratory disorders, heart disease, metabolic disorders, autoimmunity, and cancer. However, the molecular mechanisms governing low-level inflammation are largely unknown. We discovered that Bcl-2-interacting killer (Bik) deficiency causes low-level inflammation even at baseline and the development of spontaneous emphysema in female but not male mice. Similarly, a single nucleotide polymorphism that reduced Bik levels was associated with increased inflammation and enhanced decline in lung function in humans. Transgenic expression of Bik in the airways of Bik-deficient mice inhibited allergen- or LPS-induced lung inflammation and reversed emphysema in female mice. Bik deficiency increased nuclear but not cytosolic p65 levels because Bik, by modifying the BH4 domain of Bcl-2, interacted with regulatory particle non-ATPase 1 (RPN1) and RPN2 and enhanced proteasomal degradation of nuclear proteins. Bik deficiency increased inflammation primarily in females because Bcl-2 and Bik levels were reduced in lung tissues and airway cells of female compared with male mice. Therefore, controlling low-grade inflammation by modifying the unappreciated role of Bik and Bcl-2 in facilitating proteasomal degradation of nuclear proteins may be crucial in treating chronic age-related diseases.

Dynamic trajectories of life-threatening mass effect in patients with large middle cerebral artery stroke

Background Life-threatening, space-occupying mass effect due to cerebral edema and/or hemorrhagic transformation is an early complication of patients with middle cerebral artery ( MCA ) stroke. Little is known about longitudinal trajectories of laboratory and vital signs leading up to radiographic and clinical deterioration related to this mass effect. Methods We curated a granular retrospective dataset of 635 patients with large middle cerebral artery ( MCA ) stroke totaling 108,547 data points for repeated measurements of 10 covariates, and 40 time-independent covariates. We assessed longitudinal trajectories of the 10 longitudinal variables during the 72 hours preceding three outcomes representative of life-threatening mass effect: midline shift ( MLS ) \(\ge\)5mm, pineal gland shift ( PGS ) \(>\)4mm, and decompressive hemicraniectomy ( DHC ). We used a "backward looking" trajectory approach. Patients were aligned according to the time of outcome occurrence and the trajectory of each variable was assessed prior to that outcome by accounting for both cases and non-cases. Results Of 635 patients, 49% were female, and mean age was 69 years. Thirty five percent of patients had MLS \(\ge\)5mm, 24.1% had PGS \(>\)4mm, and DHC occurred in 10.7%. For the three outcomes of interest, backward-looking trajectories showed mild increases in white blood cell count (10 up to 11 K/UL within 72 hours), temperature (up to half a degree within 24 hours), and sodium (1-3 mEq/L within 24 hours) leading up to outcomes. We also observed a decrease in heart rate (75 - 65 beats per minute) 24 hours prior to DHC. Conclusions Univariable longitudinal profiling showed that temperature, white blood cell count, and sodium increase prior to radiographic and clinical indicators of space-occupying mass effect. These findings will inform development of multivariable dynamic risk models to aid prediction of life-threatening space-occupying mass effect.

Clonal Hematopoiesis of Indeterminate Potential (CHIP) and Incident Type 2 Diabetes Risk

OBJECTIVE

Clonal hematopoiesis of indeterminate potential (CHIP) is an aging-related accumulation of somatic mutations in hematopoietic stem cells, leading to clonal expansion. CHIP presence has been implicated in atherosclerotic coronary heart disease (CHD) and all-cause mortality, but its association with incident type 2 diabetes (T2D) is unknown. We hypothesized that CHIP is associated with elevated risk of T2D.

RESEARCH DESIGN AND METHODS

CHIP was derived from whole-genome sequencing of blood DNA in the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine (TOPMed) prospective cohorts. We performed analysis for 17,637 participants from six cohorts, without prior T2D, cardiovascular disease, or cancer. We evaluated baseline CHIP versus no CHIP prevalence with incident T2D, including associations with DNMT3A, TET2, ASXL1, JAK2, and TP53 variants. We estimated multivariable-adjusted hazard ratios (HRs) and 95% CIs with adjustment for age, sex, BMI, smoking, alcohol, education, self-reported race/ethnicity, and combined cohorts' estimates via fixed-effects meta-analysis.

RESULTS

Mean (SD) age was 63.4 (11.5) years, 76% were female, and CHIP prevalence was 6.0% (n = 1,055) at baseline. T2D was diagnosed in n = 2,467 over mean follow-up of 9.8 years. Participants with CHIP had 23% (CI 1.04, 1.45) higher risk of T2D than those with no CHIP. Specifically, higher risk was for TET2 (HR 1.48; CI 1.05, 2.08) and ASXL1 (HR 1.76; CI 1.03, 2.99) mutations; DNMT3A was nonsignificant (HR 1.15; CI 0.93, 1.43). Statistical power was limited for JAK2 and TP53 analyses.

CONCLUSIONS

CHIP was associated with higher incidence of T2D. CHIP mutations located on genes implicated in CHD and mortality were also related to T2D, suggesting shared aging-related pathology.

Multi-tissue epigenetic analysis identifies distinct associations underlying insulin resistance and Alzheimer's disease at CPT1A locus

BACKGROUND

Insulin resistance (IR) is a major risk factor for Alzheimer's disease (AD) dementia. The mechanisms by which IR predisposes to AD are not well-understood. Epigenetic studies may help identify molecular signatures of IR associated with AD, thus improving our understanding of the biological and regulatory mechanisms linking IR and AD.

METHODS

We conducted an epigenome-wide association study of IR, quantified using the homeostatic model assessment of IR (HOMA-IR) and adjusted for body mass index, in 3,167 participants from the Framingham Heart Study (FHS) without type 2 diabetes at the time of blood draw used for methylation measurement. We identified DNA methylation markers associated with IR at the genome-wide level accounting for multiple testing (P < 1.1 × 10-7) and evaluated their association with neurological traits in participants from the FHS (N = 3040) and the Religious Orders Study/Memory and Aging Project (ROSMAP, N = 707). DNA methylation profiles were measured in blood (FHS) or dorsolateral prefrontal cortex (ROSMAP) using the Illumina HumanMethylation450 BeadChip. Linear regressions (ROSMAP) or mixed-effects models accounting for familial relatedness (FHS) adjusted for age, sex, cohort, self-reported race, batch, and cell type proportions were used to assess associations between DNA methylation and neurological traits accounting for multiple testing.

RESULTS

We confirmed the strong association of blood DNA methylation with IR at three loci (cg17901584-DHCR24, cg17058475-CPT1A, cg00574958-CPT1A, and cg06500161-ABCG1). In FHS, higher levels of blood DNA methylation at cg00574958 and cg17058475 were both associated with lower IR (P = 2.4 × 10-11 and P = 9.0 × 10-8), larger total brain volumes (P = 0.03 and P = 9.7 × 10-4), and smaller log lateral ventricular volumes (P = 0.07 and P = 0.03). In ROSMAP, higher levels of brain DNA methylation at the same two CPT1A markers were associated with greater risk of cognitive impairment (P = 0.005 and P = 0.02) and higher AD-related indices (CERAD score: P = 5 × 10-4 and 0.001; Braak stage: P = 0.004 and P = 0.01).

CONCLUSIONS

Our results suggest potentially distinct epigenetic regulatory mechanisms between peripheral blood and dorsolateral prefrontal cortex tissues underlying IR and AD at CPT1A locus.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk

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<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P=0.11) or in the reverse direction (β=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (β=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Investigating Associations of Omega-3 Fatty Acids, Lung Function Decline, and Airway Obstruction

Rationale: Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have antiinflammatory properties and may benefit lung health. Objectives: To investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in a diverse sample of adults from general-population cohorts. Methods: Complementary study designs: 1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV1 and FVC measures in the NHLBI Pooled Cohorts Study and 2) two-sample Mendelian randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters. Measurements and Main Results: The longitudinal study found that higher omega-3 fatty acid levels were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for the most metabolically downstream omega-3 fatty acid, docosahexaenoic acid (DHA). An increase in DHA of 1% of total fatty acids was associated with attenuations of 1.4 ml/yr for FEV1 (95% confidence interval [CI], 1.1-1.8) and 2.0 ml/yr for FVC (95% CI, 1.6-2.4) and a 7% lower incidence of spirometry-defined airway obstruction (95% CI, 0.89-0.97). DHA associations persisted across sexes and smoking histories and in Black, White, and Hispanic participants, with associations of the largest magnitude in former smokers and Hispanic participants. The MR study showed similar trends toward positive associations of genetically predicted downstream omega-3 fatty acids with FEV1 and FVC. Conclusions: The longitudinal and MR studies provide evidence supporting beneficial effects of higher levels of downstream omega-3 fatty acids, especially DHA, on lung health.

Bronchial epithelial gene expression and interstitial lung abnormalities

INTRODUCTION

Interstitial lung abnormalities (ILA) often represent early fibrotic changes that can portend a progressive fibrotic phenotype. In particular, the fibrotic subtype of ILA is associated with increased mortality and rapid decline in lung function. Understanding the differential gene expression that occurs in the lungs of participants with fibrotic ILA may provide insight into development of a useful biomarker for early detection and therapeutic targets for progressive pulmonary fibrosis.

METHODS

Measures of ILA and gene expression data were available in 213 participants in the Detection of Early Lung Cancer Among Military Personnel (DECAMP1 and DECAMP2) cohorts. ILA was defined using Fleischner Society guidelines and determined by sequential reading of computed tomography (CT) scans. Primary analysis focused on comparing gene expression in ILA with usual interstitial pneumonia (UIP) pattern with those with no ILA.

RESULTS

ILA was present in 51 (24%) participants, of which 16 (7%) were subtyped as ILA with a UIP pattern. One gene, pro platelet basic protein (PPBP) and seventeen pathways (e.g. TNF-α signalling) were significantly differentially expressed between those with a probable or definite UIP pattern of ILA compared to those without ILA. 16 of these 17 pathways, but no individual gene, met significance when comparing those with ILA to those without ILA.

CONCLUSION

Our study demonstrates that abnormal inflammatory processes are apparent in the bronchial airway gene expression profiles of smokers with and without lung cancer with ILA. Future studies with larger and more diverse populations will be needed to confirm these findings.

A Type 1 Diabetes Polygenic Score Is Not Associated With Prevalent Type 2 Diabetes in Large Population Studies

Context

Both type 1 diabetes (T1D) and type 2 diabetes (T2D) have significant genetic contributions to risk and understanding their overlap can offer clinical insight.

Objective

We examined whether a T1D polygenic score (PS) was associated with a diagnosis of T2D in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium.

Methods

We constructed a T1D PS using 79 known single nucleotide polymorphisms associated with T1D risk. We analyzed 13 792 T2D cases and 14 169 controls from CHARGE cohorts to determine the association between the T1D PS and T2D prevalence. We validated findings in an independent sample of 2256 T2D cases and 27 052 controls from the Mass General Brigham Biobank (MGB Biobank). As secondary analyses in 5228 T2D cases from CHARGE, we used multivariable regression models to assess the association of the T1D PS with clinical outcomes associated with T1D.

Results

The T1D PS was not associated with T2D both in CHARGE (P = .15) and in the MGB Biobank (P = .87). The partitioned human leukocyte antigens only PS was associated with T2D in CHARGE (OR 1.02 per 1 SD increase in PS, 95% CI 1.01-1.03, P = .006) but not in the MGB Biobank. The T1D PS was weakly associated with insulin use (OR 1.007, 95% CI 1.001-1.012, P = .03) in CHARGE T2D cases but not with other outcomes.

Conclusion

In large biobank samples, a common variant PS for T1D was not consistently associated with prevalent T2D. However, possible heterogeneity in T2D cannot be ruled out and future studies are needed do subphenotyping.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Association of Common and Rare Variants with Alzheimer's Disease in over 13,000 Diverse Individuals with Whole-Genome Sequencing from the Alzheimer's Disease Sequencing Project

Alzheimer's Disease (AD) is a common disorder of the elderly that is both highly heritable and genetically heterogeneous. Here, we investigated the association between AD and both common variants and aggregates of rare coding and noncoding variants in 13,371 individuals of diverse ancestry with whole genome sequence (WGS) data. Pooled-population analyses identified genetic variants in or near APOE, BIN1, and LINC00320 significantly associated with AD (p < 5×10-8). Population-specific analyses identified a haplotype on chromosome 14 including PSEN1 associated with AD in Hispanics, further supported by aggregate testing of rare coding and noncoding variants in this region. Finally, we observed suggestive associations (p < 5×10-5) of aggregates of rare coding rare variants in ABCA7 among non-Hispanic Whites (p=5.4×10-6), and rare noncoding variants in the promoter of TOMM40 distinct of APOE in pooled-population analyses (p=7.2×10-8). Complementary pooled-population and population-specific analyses offered unique insights into the genetic architecture of AD.

GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification.

Key variants via Alzheimer's Disease Sequencing Project whole genome sequence data

INTRODUCTION

Genome-wide association studies (GWAS) have identified loci associated with Alzheimer's disease (AD) but did not identify specific causal genes or variants within those loci. Analysis of whole genome sequence (WGS) data, which interrogates the entire genome and captures rare variations, may identify causal variants within GWAS loci.

METHODS

We performed single common variant association analysis and rare variant aggregate analyses in the pooled population (N cases=2,184, N controls=2,383) and targeted analyses in sub-populations using WGS data from the Alzheimer's Disease Sequencing Project (ADSP). The analyses were restricted to variants within 100 kb of 83 previously identified GWAS lead variants.

RESULTS

Seventeen variants were significantly associated with AD within five genomic regions implicating the genes OARD1/NFYA/TREML1, JAZF1, FERMT2, and SLC24A4. KAT8 was implicated by both single variant and rare variant aggregate analyses.

DISCUSSION

This study demonstrates the utility of leveraging WGS to gain insights into AD loci identified via GWAS.

Systemic Markers of Lung Function and Forced Expiratory Volume in 1 Second Decline across Diverse Cohorts

Rationale: Chronic obstructive pulmonary disease (COPD) is a complex disease characterized by airway obstruction and accelerated lung function decline. Our understanding of systemic protein biomarkers associated with COPD remains incomplete. Objectives: To determine what proteins and pathways are associated with impaired pulmonary function in a diverse population. Methods: We studied 6,722 participants across six cohort studies with both aptamer-based proteomic and spirometry data (4,566 predominantly White participants in a discovery analysis and 2,156 African American cohort participants in a validation). In linear regression models, we examined protein associations with baseline forced expiratory volume in 1 second (FEV1) and FEV1/forced vital capacity (FVC). In linear mixed effects models, we investigated the associations of baseline protein levels with rate of FEV1 decline (ml/yr) in 2,777 participants with up to 7 years of follow-up spirometry. Results: We identified 254 proteins associated with FEV1 in our discovery analyses, with 80 proteins validated in the Jackson Heart Study. Novel validated protein associations include kallistatin serine protease inhibitor, growth differentiation factor 2, and tumor necrosis factor-like weak inducer of apoptosis (discovery β = 0.0561, Q = 4.05 × 10-10; β  = 0.0421, Q = 1.12 × 10-3; and β = 0.0358, Q = 1.67 × 10-3, respectively). In longitudinal analyses within cohorts with follow-up spirometry, we identified 15 proteins associated with FEV1 decline (Q < 0.05), including elafin leukocyte elastase inhibitor and mucin-associated TFF2 (trefoil factor 2; β = -4.3 ml/yr, Q = 0.049; β = -6.1 ml/yr, Q = 0.032, respectively). Pathways and processes highlighted by our study include aberrant extracellular matrix remodeling, enhanced innate immune response, dysregulation of angiogenesis, and coagulation. Conclusions: In this study, we identify and validate novel biomarkers and pathways associated with lung function traits in a racially diverse population. In addition, we identify novel protein markers associated with FEV1 decline. Several protein findings are supported by previously reported genetic signals, highlighting the plausibility of certain biologic pathways. These novel proteins might represent markers for risk stratification, as well as novel molecular targets for treatment of COPD.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus

BACKGROUND

Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD.

METHODS

From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D.

RESULTS

A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler

A primary obstacle in translating genetic associations with disease into therapeutic strategies is elucidating the cellular programs affected by genetic risk variants and effector genes. Here, we introduce LipocyteProfiler, a cardiometabolic-disease-oriented high-content image-based profiling tool that enables evaluation of thousands of morphological and cellular profiles that can be systematically linked to genes and genetic variants relevant to cardiometabolic disease. We show that LipocyteProfiler allows surveillance of diverse cellular programs by generating rich context- and process-specific cellular profiles across hepatocyte and adipocyte cell-state transitions. We use LipocyteProfiler to identify known and novel cellular mechanisms altered by polygenic risk of metabolic disease, including insulin resistance, fat distribution, and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse deep phenotypic profiling in lipocytes and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

The power of TOPMed imputation for the discovery of Latino-enriched rare variants associated with type 2 diabetes

AIMS/HYPOTHESIS

The Latino population has been systematically underrepresented in large-scale genetic analyses, and previous studies have relied on the imputation of ungenotyped variants based on the 1000 Genomes (1000G) imputation panel, which results in suboptimal capture of low-frequency or Latino-enriched variants. The National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) released the largest multi-ancestry genotype reference panel representing a unique opportunity to analyse rare genetic variations in the Latino population. We hypothesise that a more comprehensive analysis of low/rare variation using the TOPMed panel would improve our knowledge of the genetics of type 2 diabetes in the Latino population.

METHODS

We evaluated the TOPMed imputation performance using genotyping array and whole-exome sequence data in six Latino cohorts. To evaluate the ability of TOPMed imputation to increase the number of identified loci, we performed a Latino type 2 diabetes genome-wide association study (GWAS) meta-analysis in 8150 individuals with type 2 diabetes and 10,735 control individuals and replicated the results in six additional cohorts including whole-genome sequence data from the All of Us cohort.

RESULTS

Compared with imputation with 1000G, the TOPMed panel improved the identification of rare and low-frequency variants. We identified 26 genome-wide significant signals including a novel variant (minor allele frequency 1.7%; OR 1.37, p=3.4 × 10-9). A Latino-tailored polygenic score constructed from our data and GWAS data from East Asian and European populations improved the prediction accuracy in a Latino target dataset, explaining up to 7.6% of the type 2 diabetes risk variance.

CONCLUSIONS/INTERPRETATION

Our results demonstrate the utility of TOPMed imputation for identifying low-frequency variants in understudied populations, leading to the discovery of novel disease associations and the improvement of polygenic scores.

DATA AVAILABILITY

Full summary statistics are available through the Common Metabolic Diseases Knowledge Portal ( https://t2d.hugeamp.org/downloads.html ) and through the GWAS catalog ( https://www.ebi.ac.uk/gwas/ , accession ID: GCST90255648). Polygenic score (PS) weights for each ancestry are available via the PGS catalog ( https://www.pgscatalog.org , publication ID: PGP000445, scores IDs: PGS003443, PGS003444 and PGS003445).

Genome-wide association study and functional characterization identifies candidate genes for insulin-stimulated glucose uptake

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10-8) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits.

Identifying important gene signatures of BMI using network structure-aided nonparametric quantile regression

We focus on identifying genomics risk factors of higher body mass index (BMI) incorporating a priori information, such as biological pathways. However, the commonly used methods to incorporate prior information provide a model for the mean function of the outcome and rely on unmet assumptions. To address these concerns, we propose a method for nonparametric additive quantile regression with network regularization to incorporate the information encoded by known networks. To account for nonlinear associations, we approximate the unknown additive functional effect of each predictor with the expansion of a B-spline basis. We implement the group Lasso penalty to obtain a sparse model. We define the network-constrained penalty by the totalℓ 2 $$ {\ell}_2 $$ norm of the difference between the effect functions of any two linked genes in the known network. We further propose an efficient computation procedure to solve the optimization problem that arises in our model. Simulation studies show that our proposed method performs well in identifying more truly associated genes and less falsely associated genes than alternative approaches. We apply the proposed method to analyze the microarray gene-expression dataset in the Framingham Heart Study and identify several 75 percentile BMI associated genes. In conclusion, our proposed approach efficiently identifies the outcome-associated variables in a nonparametric additive quantile regression framework by leveraging known network information.

Rationale and Design of the Diet, CKD, and Apolipoprotein L1 Study in Low-Income and Middle-Income Countries

Introduction

Diet, chronic kidney disease (CKD), and Apolipoprotein L1 (APOL1) (DCA) Study is examining the role of dietary factors in CKD progression and APOL1 nephropathy. We describe enrollment and retention efforts and highlight facilitators and barriers to enrollment and operational challenges, as well as accommodations made in the study protocol.

Methods

The DCA study is enrolling participants in 7 centers in West Africa. Participants who consented were invited to complete dietary recalls and 24-hour urine collections in year 1. We conducted focus groups and semistructured interviews among study personnel to identify facilitators and barriers to enrollment as well as retention and operational challenges in the execution of the study protocol. We analyzed emerging themes using content analyses.

Results

A total of 712 participants were enrolled in 18 months with 1256 24-hour urine and 1260 dietary recalls. Barriers to enrollment were the following: (i) a lack of understanding of research, (ii) the burden of research visits, and (iii) incorporating cultural and traditional nuances when designing research protocols. Factors facilitating enrollment were the following: (i) designing convenient research visits, (ii) building rapport and increased communication between the research team and participants, and (iii) cultural sensitivity - adapting research protocols for the populations involved. Offering home visits, providing free dietary counseling, reducing the volume of study blood collection, and reducing the frequency of visits were some changes made in the study protocol that increased participant satisfaction.

Conclusion

Adopting a participant-centered approach with accommodations in the protocol for cultural adaptability and incorporating participant feedback is vital for carrying out research in low-income and middle-income regions.

Multi-ancestry genome-wide association analyses improve resolution of genes and pathways influencing lung function and chronic obstructive pulmonary disease risk

Lung-function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry genome-wide association meta-analysis of lung function to date, comprising 580,869 participants, we identified 1,020 independent association signals implicating 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies for selected associated variants as well as trait and pathway-specific genetic risk scores to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Genetic Effect on Body Mass Index and Cardiovascular Disease Across Generations

BACKGROUND

Whether genetics contribute to the rising prevalence of obesity or its cardiovascular consequences in today's obesogenic environment remains unclear. We sought to determine whether the effects of a higher aggregate genetic burden of obesity risk on body mass index (BMI) or cardiovascular disease (CVD) differed by birth year.

METHODS

We split the FHS (Framingham Heart Study) into 4 equally sized birth cohorts (birth year before 1932, 1932 to 1946, 1947 to 1959, and after 1960). We modeled a genetic predisposition to obesity using an additive genetic risk score (GRS) of 941 BMI-associated variants and tested for GRS-birth year interaction on log-BMI (outcome) when participants were around 50 years old (N=7693). We repeated the analysis using a GRS of 109 BMI-associated variants that increased CVD risk factors (type 2 diabetes, blood pressure, total cholesterol, and high-density lipoprotein) in addition to BMI. We then evaluated whether the effects of the BMI GRSs on CVD risk differed by birth cohort when participants were around 60 years old (N=5493).

RESULTS

Compared with participants born before 1932 (mean age, 50.8 yrs [2.4]), those born after 1960 (mean age, 43.3 years [4.5]) had higher BMI (median, 25.4 [23.3-28.0] kg/m2 versus 26.9 [interquartile range, 23.7-30.6] kg/m2). The effect of the 941-variant BMI GRS on BMI and CVD risk was stronger in people who were born in later years (GRS-birth year interaction: P=0.0007 and P=0.04 respectively).

CONCLUSIONS

The significant GRS-birth year interactions indicate that common genetic variants have larger effects on middle-age BMI and CVD risk in people born more recently. These findings suggest that the increasingly obesogenic environment may amplify the impact of genetics on the risk of obesity and possibly its cardiovascular consequences.

Investigating associations of omega-3 fatty acids, lung function decline, and airway obstruction

Rationale

Inflammation contributes to lung function decline and the development of chronic obstructive pulmonary disease. Omega-3 fatty acids have anti-inflammatory properties and may benefit lung health.

Objectives

Investigate associations of omega-3 fatty acids with lung function decline and incident airway obstruction in adults of diverse races/ethnicities from general population cohorts.

Methods

Complementary study designs: (1) longitudinal study of plasma phospholipid omega-3 fatty acids and repeated FEV 1 and FVC measures in the National Heart, Lung, and Blood Institute Pooled Cohorts Study, and (2) two-sample Mendelian Randomization (MR) study of genetically predicted omega-3 fatty acids and lung function parameters.

Measurements and Main Results

The longitudinal study found that higher omega-3 fatty acid concentrations were associated with attenuated lung function decline in 15,063 participants, with the largest effect sizes for docosahexaenoic acid (DHA). One standard deviation higher DHA was associated with an attenuation of 1.8 mL/year for FEV 1 (95% confidence interval [CI] 1.3-2.2) and 2.4 mL/year for FVC (95% CI 1.9-3.0). One standard deviation higher DHA was also associated with a 9% lower incidence of spirometry-defined airway obstruction (95% CI 0.86-0.97). DHA associations persisted across sexes, smoking histories, and Black, white and Hispanic participants, with the largest magnitude associations in former smokers and Hispanics. The MR study showed positive associations of genetically predicted omega-3 fatty acids with FEV 1 and FVC, with statistically significant findings across multiple MR methods.

Conclusions

The longitudinal and MR studies provide evidence supporting beneficial effects of higher circulating omega-3 fatty acids, especially DHA, on lung health.

Assisted clustering of gene expression data using regulatory data from partially overlapping sets of individuals

BACKGROUND

As omics measurements profiled on different molecular layers are interconnected, integrative approaches that incorporate the regulatory effect from multi-level omics data are needed. When the multi-level omics data are from the same individuals, gene expression (GE) clusters can be identified using information from regulators like genetic variants and DNA methylation. When the multi-level omics data are from different individuals, the choice of integration approaches is limited.

METHODS

We developed an approach to improve GE clustering from microarray data by integrating regulatory data from different but partially overlapping sets of individuals. We achieve this through (1) decomposing gene expression into the regulated component and the other component that is not regulated by measured factors, (2) optimizing the clustering goodness-of-fit objective function. We do not require the availability of different omics measurements on all individuals. A certain amount of individual overlap between GE data and the regulatory data is adequate for modeling the regulation, thus improving GE clustering.

RESULTS

A simulation study shows that the performance of the proposed approach depends on the strength of the GE-regulator relationship, degree of missingness, data dimensionality, sample size, and the number of clusters. Across the various simulation settings, the proposed method shows competitive performance in terms of accuracy compared to the alternative K-means clustering method, especially when the clustering structure is due mostly to the regulated component, rather than the unregulated component. We further validate the approach with an application to 8,902 Framingham Heart Study participants with data on up to 17,873 genes and regulation information of DNA methylation and genotype from different but partially overlapping sets of participants. We identify clustering structures of genes associated with pulmonary function while incorporating the predicted regulation effect from the measured regulators. We further investigate the over-representation of these GE clusters in pathways of other diseases that may be related to lung function and respiratory health.

CONCLUSION

We propose a novel approach for clustering GE with the assistance of regulatory data that allowed for different but partially overlapping sets of individuals to be included in different omics data.

Quantitative pupillometry and radiographic markers of intracranial midline shift: A pilot study

Background

Asymmetric pupil reactivity or size can be early clinical indicators of midbrain compression due to supratentorial ischemic stroke or primary intraparenchymal hemorrhage (IPH). Radiographic midline shift is associated with worse functional outcomes and life-saving interventions. Better understanding of quantitative pupil characteristics would be a non-invasive, safe, and cost-effective way to improve identification of life-threatening mass effect and resource utilization of emergent radiographic imaging. We aimed to better characterize the association between midline shift at various anatomic levels and quantitative pupil characteristics.

Methods

We conducted a multicenter retrospective study of brain CT images within 75 min of a quantitative pupil observation from patients admitted to Neuro-ICUs between 2016 and 2020 with large (>1/3 of the middle cerebral artery territory) acute supratentorial ischemic stroke or primary IPH > 30 mm³. For each image, we measured midline shift at the septum pellucidum (MLS-SP), pineal gland shift (PGS), the ratio of the ipsilateral to contralateral midbrain width (IMW/CMW), and other exploratory markers of radiographic shift/compression. Pupil reactivity was measured using an automated infrared pupillometer (NeurOptics^®, Inc.), specifically the proprietary algorithm for Neurological Pupil Index^® (NPi). We used rank-normalization and linear mixed-effects models, stratified by diagnosis and hemorrhagic conversion, to test associations of radiographic markers of shift and asymmetric pupil reactivity (Diff NPi), adjusting for age, lesion volume, Glasgow Coma Scale, and osmotic medications.

Results

Of 53 patients with 74 CT images, 26 (49.1%) were female, and median age was 67 years. MLS-SP and PGS were greater in patients with IPH, compared to patients with ischemic stroke (6.2 v. 4.0 mm, 5.6 v. 3.4 mm, respectively). We found no significant associations between pupil reactivity and the radiographic markers of shift when adjusting for confounders. However, we found potentially relevant relationships between MLS-SP and Diff NPi in our IPH cohort (β = 0.11, SE 0.04, P = 0.01), and PGS and Diff NPi in the ischemic stroke cohort (β = 0.16, SE 0.09, P = 0.07).

Conclusion

We found the relationship between midline shift and asymmetric pupil reactivity may differ between IPH and ischemic stroke. Our study may serve as necessary preliminary data to guide further prospective investigation into how clinical manifestations of radiographic midline shift differ by diagnosis and proximity to the midbrain.

Exploiting family history in aggregation unit-based genetic association tests

The development of sequencing technology calls for new powerful methods to detect disease associations and lower the cost of sequencing studies. Family history (FH) contains information on disease status of relatives, adding valuable information about the probands' health problems and risk of diseases. Incorporating data from FH is a cost-effective way to improve statistical evidence in genetic studies, and moreover, overcomes limitations in study designs with insufficient cases or missing genotype information for association analysis. We proposed family history aggregation unit-based test (FHAT) and optimal FHAT (FHAT-O) to exploit available FH for rare variant association analysis. Moreover, we extended liability threshold model of case-control status and FH (LT-FH) method in aggregated unit-based methods and compared that with FHAT and FHAT-O. The computational efficiency and flexibility of the FHAT and FHAT-O were demonstrated through both simulations and applications. We showed that FHAT, FHAT-O, and LT-FH methods offer reasonable control of the type I error unless case/control ratio is unbalanced, in which case they result in smaller inflation than that observed with conventional methods excluding FH. We also demonstrated that FHAT and FHAT-O are more powerful than LT-FH and conventional methods in many scenarios. By applying FHAT and FHAT-O to the analysis of all cause dementia and hypertension using the exome sequencing data from the UK Biobank, we showed that our methods can improve significance for known regions. Furthermore, we replicated the previous associations in all cause dementia and hypertension and detected novel regions through the exome-wide analysis.

Assessing the contribution of rare genetic variants to phenotypes of chronic obstructive pulmonary disease using whole-genome sequence data

RATIONALE

Genetic 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.

OBJECTIVES

We sought to assess the heritability of COPD and lung function using whole-genome sequence data from the Trans-Omics for Precision Medicine program.

METHODS

Using 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 RESULTS

In 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 < 0.05). In African-Americans, common variant heritability was similar to European ancestry participants, but lower sample size precluded calculation of rare variant heritability.

CONCLUSIONS

Our study provides updated and unbiased estimates of heritability for COPD and lung function, and suggests an important contribution of rare variants. Larger studies of more diverse ancestry will improve accuracy of these estimates.

Leveraging family history in genetic association analyses of binary traits

BACKGROUND

Considering relatives' health history in logistic regression for case-control genome-wide association studies (CC-GWAS) may provide new information that increases accuracy and power to detect disease associated genetic variants. We conducted simulations and analyzed type 2 diabetes (T2D) data from the Framingham Heart Study (FHS) to compare two methods, liability threshold model conditional on both case-control status and family history (LT-FH) and Fam-meta, which incorporate family history into CC-GWAS.

RESULTS

In our simulation scenario of trait with modest T2D heritability (h² = 0.28), variant minor allele frequency ranging from 1% to 50%, and 1% of phenotype variance explained by the genetic variants, Fam-meta had the highest overall power, while both methods incorporating family history were more powerful than CC-GWAS. All three methods had controlled type I error rates, while LT-FH was the most conservative with a lower-than-expected error rate. In addition, we observed a substantial increase in power of the two familial history methods compared to CC-GWAS when the prevalence of the phenotype increased with age. Furthermore, we showed that, when only the phenotypes of more distant relatives were available, Fam-meta still remained more powerful than CC-GWAS, confirming that leveraging disease history of both close and distant relatives can increase power of association analyses. Using FHS data, we confirmed the well-known association of TCF7L2 region with T2D at the genome-wide threshold of P-value < 5 × 10^-8, and both familial history methods increased the significance of the region compared to CC-GWAS. We identified two loci at 5q35 (ADAMTS2) and 5q23 (PRR16), not previously reported for T2D using CC-GWAS and Fam-meta; both genes play a role in cardiovascular diseases. Additionally, CC-GWAS detected one more significant locus at 13q31 (GPC6) reported associated with T2D-related traits.

CONCLUSIONS

Overall, LT-FH and Fam-meta had higher power than CC-GWAS in simulations, especially using phenotypes that were more prevalent in older age groups, and both methods detected known genetic variants with lower P-values in real data application, highlighting the benefits of including family history in genetic association studies.

Apolipoprotein L1 Genotypes and the Association of Urinary Potassium Excretion with CKD Progression

BACKGROUND AND OBJECTIVES

Progressive CKD in Black individuals is strongly associated with polymorphisms in the APOL1 gene, but it is unknown whether dietary risk factors for CKD progression vary in high- versus low-risk APOL1 genotypes. We investigated if APOL1 genotypes modify associations of dietary potassium and sodium with CKD progression and death.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

We analyzed 1399 self-identified Black participants enrolled in the Chronic Renal Insufficiency Cohort from April 2003 to September 2008. Exposures were calibrated 24-hour urine potassium and sodium excretion. The primary outcome was CKD progression defined as the time to 50% decline in eGFR or kidney failure. The secondary outcome was CKD progression or death. We tested for an interaction between urinary potassium and sodium excretion and APOL1 genotypes.

RESULTS

Median 24-hour urinary sodium and potassium excretions in Black participants were 150 mmol (interquartile range, 118-188) and 43 mmol (interquartile range, 35-54), respectively. Individuals with high- and low-risk APOL1 genotypes numbered 276 (20%) and 1104 (79%), respectively. After a median follow-up of 5.23 years, CKD progression events equaled 605, and after 7.29 years, CKD progression and death events equaled 868. There was significant interaction between APOL1 genotypes and urinary potassium excretion with CKD progression and CKD progression or death (P=0.003 and P=0.03, respectively). In those with high-risk APOL1 genotypes, higher urinary potassium excretion was associated with a lower risk of CKD progression (quartiles 2-4 versus 1: hazard ratio, 0.83; 95% confidence interval, 0.50 to 1.39; hazard ratio, 0.54; 95% confidence interval, 0.31 to 0.93; and hazard ratio, 0.50; 95% confidence interval, 0.27 to 0.93, respectively). In the low-risk APOL1 genotypes, higher urinary potassium excretion was associated with a higher risk of CKD progression (quartiles 2-4 versus 1: hazard ratio, 1.01; 95% confidence interval, 0.75 to 1.36; hazard ratio, 1.23; 95% confidence interval, 0.91 to 1.66; and hazard ratio, 1.53; 95% confidence interval, 1.12 to 2.09, respectively). We found no interaction between APOL1 genotypes and urinary sodium excretion with CKD outcomes.

CONCLUSIONS

Higher urinary potassium excretion was associated with lower versus higher risk of CKD progression in APOL1 high-risk and low-risk genotypes, respectively.

A polygenic risk score and age of diagnosis of COPD

BACKGROUND

Genetic susceptibility may be associated with earlier onset of chronic obstructive pulmonary disease (COPD). We hypothesised that a polygenic risk score (PRS) for COPD would be associated with earlier age of diagnosis of COPD.

METHODS

In 6647 non-Hispanic White (NHW) and 2464 African American (AA) participants from COPDGene, and 6812 participants from the Framingham Heart Study (FHS), we tested the relationship of the PRS and age of COPD diagnosis. Age at diagnosis was determined by: 1) self-reported age at COPD diagnosis or 2) age at visits when moderate-to-severe airflow limitation (Global Initiative for Chronic Obstructive Lung Disease (GOLD) grade 2-4) was observed on spirometry. We used Cox regression to examine the overall and time-dependent effects of the PRS on incident COPD. In the COPDGene study, we also examined the PRS's predictive value for COPD at age <50 years (COPD50) using logistic regression and area under the curve (AUC) analyses, with and without the addition of other risk factors present at early life (e.g. childhood asthma).

RESULTS

In Cox models, the PRS demonstrated age-dependent associations with incident COPD, with larger effects at younger ages in both cohorts. The PRS was associated with COPD50 (OR 1.55 (95% CI 1.41-1.71) for NHW, OR 1.23 (95% CI 1.05-1.43) for AA and OR 2.47 (95% CI 2.12-2.88) for FHS participants). In COPDGene, adding the PRS to known early-life risk factors improved prediction of COPD50 in NHW (AUC 0.69 versus 0.74; p<0.0001) and AA (AUC 0.61 versus 0.64; p=0.04) participants.

CONCLUSIONS

A COPD PRS is associated with earlier age of diagnosis of COPD and retains predictive value when added to known early-life risk factors.

Epidemiology and impact of frailty in patients with atrial fibrillation in Europe

BACKGROUND

Frailty is a medical syndrome characterised by reduced physiological reserve and increased vulnerability to stressors. Data regarding the relationship between frailty and atrial fibrillation (AF) are still inconsistent.

OBJECTIVES

We aim to perform a comprehensive evaluation of frailty in a large European cohort of AF patients.

METHODS

A 40-item frailty index (FI) was built according to the accumulation of deficits model in the AF patients enrolled in the ESC-EHRA EORP-AF General Long-Term Registry. Association of baseline characteristics, clinical management, quality of life, healthcare resources use and risk of outcomes with frailty was examined.

RESULTS

Among 10,177 patients [mean age (standard deviation) 69.0 (11.4) years, 4,103 (40.3%) females], 6,066 (59.6%) were pre-frail and 2,172 (21.3%) were frail, whereas only 1,939 (19.1%) were considered robust. Baseline thromboembolic and bleeding risks were independently associated with increasing FI. Frail patients with AF were less likely to be treated with oral anticoagulants (OACs) (odds ratio 0.70, 95% confidence interval 0.55-0.89), especially with non-vitamin K antagonist OACs and managed with a rhythm control strategy, compared with robust patients. Increasing frailty was associated with a higher risk for all outcomes examined, with a non-linear exponential relationship. The use of OAC was associated with a lower risk of outcomes, except in patients with very/extremely high frailty.

CONCLUSIONS

In this large cohort of AF patients, there was a high burden of frailty, influencing clinical management and risk of adverse outcomes. The clinical benefit of OAC is maintained in patients with high frailty, but not in very high/extremely frail ones.

Sugar-Sweetened Beverage Consumption and Plasma Lipoprotein Cholesterol, Apolipoprotein, and Lipoprotein Particle Size Concentrations in US Adults

BACKGROUND

Prospective cohort studies have found a relation between sugar-sweetened beverage (SSB) consumption (sodas and fruit drinks) and dyslipidemia. There is limited evidence linking SSB consumption to emerging features of dyslipidemia, which can be characterized by variation in lipoprotein particle size, remnant-like particle (RLP), and apolipoprotein concentrations.

OBJECTIVES

To examine the association between SSB consumption and plasma lipoprotein cholesterol, apolipoprotein, and lipoprotein particle size concentrations among US adults.

METHODS

We examined participants from the Framingham Offspring Study (FOS; 1987-1995, n = 3047) and the Women's Health Study (1992, n = 26,218). Concentrations of plasma LDL cholesterol, apolipoprotein B (apoB), HDL cholesterol, apolipoprotein A1 (apoA1), triglyceride (TG), and non-HDL cholesterol, as well as total cholesterol:HDL cholesterol ratio and apoB:apoA1 ratio, were quantified in both cohorts; concentrations of apolipoprotein E, apolipoprotein C3, RLP-TG, and RLP cholesterol (RLP-C) were measured in the FOS only. Lipoprotein particle sizes were calculated from nuclear magnetic resonance signals for lipoprotein particle subclass concentrations (TG-rich lipoprotein particles [TRL-Ps]: very large, large, medium, small, and very small; LDL particles [LDL-Ps]: large, medium, and small; HDL particles [HDL-Ps]: large, medium, and small). SSB consumption was estimated from food frequency questionnaire data. We examined the associations between SSB consumption and all lipoprotein and apoprotein measures in linear regression models, adjusting for confounding factors such as lifestyle, diet, and traditional lipoprotein risk factors.

RESULTS

SSB consumption was positively associated with LDL cholesterol, apoB, TG, RLP-TG, RLP-C, and non-HDL cholesterol concentrations and total cholesterol:HDL cholesterol and apoB:apoA1 ratios; and negatively associated with HDL cholesterol and apoA1 concentrations (P-trend range: <0.0001 to 0.008). After adjustment for traditional lipoprotein risk factors, SSB consumers had smaller LDL-P and HDL-P sizes; lower concentrations of large LDL-Ps and medium HDL-Ps; and higher concentrations of small LDL-Ps, small HDL-Ps, and large TRL-Ps (P-trend range: <0.0001 to 0.001).

CONCLUSIONS

Higher SSB consumption was associated with multiple emerging features of dyslipidemia that have been linked to higher cardiometabolic risk in US adults.

Interstitial lung abnormalities are associated with decreased mean telomere length

BACKGROUND

Interstitial lung abnormalities (ILA) share many features with idiopathic pulmonary fibrosis; however, it is not known if ILA are associated with decreased mean telomere length (MTL).

METHODS

Telomere length was measured with quantitative PCR in the Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) and Age Gene/Environment Susceptibility Reykjavik (AGES-Reykjavik) cohorts and Southern blot analysis was used in the Framingham Heart Study (FHS). Logistic and linear regression were used to assess the association between ILA and MTL; Cox proportional hazards models were used to assess the association between MTL and mortality.

RESULTS

In all three cohorts, ILA were associated with decreased MTL. In the COPDGene and AGES-Reykjavik cohorts, after adjustment there was greater than twofold increase in the odds of ILA when comparing the shortest quartile of telomere length to the longest quartile (OR 2.2, 95% CI 1.5-3.4, p=0.0001, and OR 2.6, 95% CI 1.4-4.9, p=0.003, respectively). In the FHS, those with ILA had shorter telomeres than those without ILA (-767 bp, 95% CI 76-1584 bp, p=0.03). Although decreased MTL was associated with chronic obstructive pulmonary disease (OR 1.3, 95% CI 1.1-1.6, p=0.01) in COPDGene, the effect estimate was less than that noted with ILA. There was no consistent association between MTL and risk of death when comparing the shortest quartile of telomere length in COPDGene and AGES-Reykjavik (HR 0.82, 95% CI 0.4-1.7, p=0.6, and HR 1.2, 95% CI 0.6-2.2, p=0.5, respectively).

CONCLUSION

ILA are associated with decreased MTL.

Pulmonary Function and Blood DNA Methylation: A Multiancestry Epigenome-Wide Association Meta-analysis

Rationale: Methylation integrates factors present at birth and modifiable across the lifespan that can influence pulmonary function. Studies are limited in scope and replication. Objectives: To conduct large-scale epigenome-wide meta-analyses of blood DNA methylation and pulmonary function. Methods: Twelve cohorts analyzed associations of methylation at cytosine-phosphate-guanine probes (CpGs), using Illumina 450K or EPIC/850K arrays, with FEV₁, FVC, and FEV₁/FVC. We performed multiancestry epigenome-wide meta-analyses (total of 17,503 individuals; 14,761 European, 2,549 African, and 193 Hispanic/Latino ancestries) and interpreted results using integrative epigenomics. Measurements and Main Results: We identified 1,267 CpGs (1,042 genes) differentially methylated (false discovery rate, <0.025) in relation to FEV₁, FVC, or FEV₁/FVC, including 1,240 novel and 73 also related to chronic obstructive pulmonary disease (1,787 cases). We found 294 CpGs unique to European or African ancestry and 395 CpGs unique to never or ever smokers. The majority of significant CpGs correlated with nearby gene expression in blood. Findings were enriched in key regulatory elements for gene function, including accessible chromatin elements, in both blood and lung. Sixty-nine implicated genes are targets of investigational or approved drugs. One example novel gene highlighted by integrative epigenomic and druggable target analysis is TNFRSF4. Mendelian randomization and colocalization analyses suggest that epigenome-wide association study signals capture causal regulatory genomic loci. Conclusions: We identified numerous novel loci differentially methylated in relation to pulmonary function; few were detected in large genome-wide association studies. Integrative analyses highlight functional relevance and potential therapeutic targets. This comprehensive discovery of potentially modifiable, novel lung function loci expands knowledge gained from genetic studies, providing insights into lung pathogenesis.

Whole genome sequence association analysis of fasting glucose and fasting insulin levels in diverse cohorts from the NHLBI TOPMed program

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.

The relationship between interstitial lung abnormalities, mortality, and multimorbidity: a cohort study

BACKGROUND

Interstitial lung abnormalities (ILAs) are associated with increased mortality. It is unclear whether multimorbidity accounts for the mortality association or how strongly ILA is associated with mortality relative to other common age-associated diseases. We determined the association of ILA with all-cause mortality adjusted for multimorbidity, compared mortality associated with ILA and prevalent cardiovascular disease (CVD), diabetes mellitus, chronic kidney disease, chronic obstructive pulmonary disease and cancer and also determined the association between ILA and these diseases.

METHODS

We measured ILA (none, indeterminant, definite) using blinded reads of CT images, prevalent chronic diseases and potential confounders in two observational cohorts, the Framingham Heart Study (FHS) (n=2449) and Age, Gene/Environment Susceptibility - Reykjavik Study (AGES-Reykjavik) (n=5180). We determined associations with mortality using Cox proportional hazards models and between ILA and diseases with multinomial logistic regression.

RESULTS

Over a median (IQR) follow-up of 8.8 (1.4) years in FHS and 12.0 (7.7) years in AGES-Reykjavik, in adjusted models, ILAs were significantly associated with increased mortality (HR, 95% CI 1.95, 1.23 to 3.08, p=0.0042, in FHS; HR 1.60, 1.41 to 1.82, p<0.0001, in AGES-Reykjavik) adjusted for multimorbidity. In both cohorts, the association of ILA with mortality was of similar magnitude to the association of most other diseases. In adjusted models, ILAs were associated only with prevalent kidney disease (OR, 95% CI 1.90, 1.01 to 3.57, p=0.0452) in FHS and with prevalent CVD (OR 1.42, 1.12 to 1.81, p=0.0040) in AGES-Reykjavik.

CONCLUSIONS

ILAs were associated with mortality adjusted for multimorbidity and were similarly associated with increased mortality compared with several common chronic diseases. ILAs were not consistently associated with the prevalence of these diseases themselves.