Cardiovascular Burden of the V142I Transthyretin Variant.

Importance: Individual cohort studies concur that the amyloidogenic V142I variant of the transthyretin (TTR) gene, present in 3% to 4% of US Black individuals, increases heart failure (HF) and mortality risk. Precisely defining carrier risk across relevant clinical outcomes and estimating population burden of disease are important given established and emerging targeted treatments. Objectives: To better define the natural history of disease in carriers across mid to late life, assess variant modifiers, and estimate cardiovascular burden to the US population. Design, Setting, and Participants: A total of 23 338 self-reported Black participants initially free from HF were included in 4 large observational studies across the US (mean [SD], 15.5 [8.2] years of follow-up). Data analysis was performed between May 2023 and February 2024. Exposure: V142I carrier status (n = 754, 3.2%). Main Outcomes and Measures: Hospitalizations for HF (including subtypes of reduced and preserved ejection fraction) and all-cause mortality. Outcomes were analyzed by generating 10-year hazard ratios for each age between 50 and 90 years. Using actuarial methods, mean survival by carrier status was estimated and applied to the 2022 US population using US Census data. Results: Among the 23 338 participants, the mean (SD) age at baseline was 62 (9) years and 76.7% were women. Ten-year carrier risk increased for HF hospitalization by age 63 years, predominantly driven by HF with reduced ejection fraction, and 10-year all-cause mortality risk increased by age 72 years. Only age (but not sex or other select variables) modified risk with the variant, with estimated reductions in longevity ranging from 1.9 years (95% CI, 0.6-3.1) at age 50 to 2.8 years (95% CI, 2.0-3.6) at age 81. Based on these data, 435 851 estimated US Black carriers between ages 50 and 95 years are projected to cumulatively lose 957 505 years of life (95% CI, 534 475-1 380 535) due to the variant. Conclusions and Relevance: Among self-reported Black individuals, male and female V142I carriers faced similar and substantial risk for HF hospitalization, predominantly with reduced ejection fraction, and death, with steep age-dependent penetrance. Delineating the individual contributions of, and complex interplay among, the V142I variant, ancestry, the social construct of race, and biological or social determinants of health to cardiovascular disease merits further investigation.

AtheroSpectrum Reveals Novel Macrophage Foam Cell Gene Signatures Associated With Atherosclerotic Cardiovascular Disease Risk.

BACKGROUND: Whereas several interventions can effectively lower lipid levels in people at risk for atherosclerotic cardiovascular disease (ASCVD), cardiovascular event risks remain, suggesting an unmet medical need to identify factors contributing to cardiovascular event risk. Monocytes and macrophages play central roles in atherosclerosis, but studies have yet to provide a detailed view of macrophage populations involved in increased ASCVD risk. METHODS: A novel macrophage foaming analytics tool, AtheroSpectrum, was developed using 2 quantitative indices depicting lipid metabolism and the inflammatory status of macrophages. A machine learning algorithm was developed to analyze gene expression patterns in the peripheral monocyte transcriptome of MESA participants (Multi-Ethnic Study of Atherosclerosis; set 1; n=911). A list of 30 genes was generated and integrated with traditional risk factors to create an ASCVD risk prediction model (30-gene cardiovascular disease risk score [CR-30]), which was subsequently validated in the remaining MESA participants (set 2; n=228); performance of CR-30 was also tested in 2 independent human atherosclerotic tissue transcriptome data sets (GTEx [Genotype-Tissue Expression] and GSE43292). RESULTS: Using single-cell transcriptomic profiles (GSE97310, GSE116240, GSE97941, and FR-FCM-Z23S), AtheroSpectrum detected 2 distinct programs in plaque macrophages-homeostatic foaming and inflammatory pathogenic foaming-the latter of which was positively associated with severity of atherosclerosis in multiple studies. A pool of 2209 pathogenic foaming genes was extracted and screened to select a subset of 30 genes correlated with cardiovascular event in MESA set 1. A cardiovascular disease risk score model (CR-30) was then developed by incorporating this gene set with traditional variables sensitive to cardiovascular event in MESA set 1 after cross-validation generalizability analysis. The performance of CR-30 was then tested in MESA set 2 (P=2.60×10-4; area under the receiver operating characteristic curve, 0.742) and 2 independent data sets (GTEx: P=7.32×10-17; area under the receiver operating characteristic curve, 0.664; GSE43292: P=7.04×10-2; area under the receiver operating characteristic curve, 0.633). Model sensitivity tests confirmed the contribution of the 30-gene panel to the prediction model (likelihood ratio test; df=31, P=0.03). CONCLUSIONS: Our novel computational program (AtheroSpectrum) identified a specific gene expression profile associated with inflammatory macrophage foam cells. A subset of 30 genes expressed in circulating monocytes jointly contributed to prediction of symptomatic atherosclerotic vascular disease. Incorporating a pathogenic foaming gene set with known risk factors can significantly strengthen the power to predict ASCVD risk. Our programs may facilitate both mechanistic investigations and development of therapeutic and prognostic strategies for ASCVD risk.

Whole Genome Sequence Analysis of the Plasma Proteome in Black Adults Provides Novel Insights Into Cardiovascular Disease.

BACKGROUND: Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants. METHODS: Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics). RESULTS: We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10-11. These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, ß=0.61±0.05, P=3.27×10-30) and MMP-3 (ß=-0.60±0.05, P=1.67×10-32), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, ß=0.34±0.04, P=1.34×10-17) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure. CONCLUSIONS: Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

MUC5B, telomere length and longitudinal quantitative interstitial lung changes: the MESA Lung Study.

BACKGROUND: The MUC5B promoter variant (rs35705950) and telomere length are linked to pulmonary fibrosis and CT-based qualitative assessments of interstitial abnormalities, but their associations with longitudinal quantitative changes of the lung interstitium among community-dwelling adults are unknown. METHODS: We used data from participants in the Multi-Ethnic Study of Atherosclerosis with high-attenuation areas (HAAs, Examinations 1-6 (2000-2018)) and MUC5B genotype (n=4552) and telomere length (n=4488) assessments. HAA was defined as the per cent of imaged lung with attenuation of -600 to -250 Hounsfield units. We used linear mixed-effects models to examine associations of MUC5B risk allele (T) and telomere length with longitudinal changes in HAAs. Joint models were used to examine associations of longitudinal changes in HAAs with death and interstitial lung disease (ILD). RESULTS: The MUC5B risk allele (T) was associated with an absolute change in HAAs of 2.60% (95% CI 0.36% to 4.86%) per 10 years overall. This association was stronger among those with a telomere length below an age-adjusted percentile of 5% (p value for interaction=0.008). A 1% increase in HAAs per year was associated with 7% increase in mortality risk (rate ratio (RR)=1.07, 95% CI 1.02 to 1.12) for overall death and 34% increase in ILD (RR=1.34, 95% CI 1.20 to 1.50). Longer baseline telomere length was cross-sectionally associated with less HAAs from baseline scans, but not with longitudinal changes in HAAs. CONCLUSIONS: Longitudinal increases in HAAs were associated with the MUC5B risk allele and a higher risk of death and ILD.

Pulmonary emphysema subtypes defined by unsupervised machine learning on CT scans.

BACKGROUND: Treatment and preventative advances for chronic obstructive pulmonary disease (COPD) have been slow due, in part, to limited subphenotypes. We tested if unsupervised machine learning on CT images would discover CT emphysema subtypes with distinct characteristics, prognoses and genetic associations. METHODS: New CT emphysema subtypes were identified by unsupervised machine learning on only the texture and location of emphysematous regions on CT scans from 2853 participants in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS), a COPD case-control study, followed by data reduction. Subtypes were compared with symptoms and physiology among 2949 participants in the population-based Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study and with prognosis among 6658 MESA participants. Associations with genome-wide single-nucleotide-polymorphisms were examined. RESULTS: The algorithm discovered six reproducible (interlearner intraclass correlation coefficient, 0.91-1.00) CT emphysema subtypes. The most common subtype in SPIROMICS, the combined bronchitis-apical subtype, was associated with chronic bronchitis, accelerated lung function decline, hospitalisations, deaths, incident airflow limitation and a gene variant near DRD1, which is implicated in mucin hypersecretion (p=1.1 ×10-8). The second, the diffuse subtype was associated with lower weight, respiratory hospitalisations and deaths, and incident airflow limitation. The third was associated with age only. The fourth and fifth visually resembled combined pulmonary fibrosis emphysema and had distinct symptoms, physiology, prognosis and genetic associations. The sixth visually resembled vanishing lung syndrome. CONCLUSION: Large-scale unsupervised machine learning on CT scans defined six reproducible, familiar CT emphysema subtypes that suggest paths to specific diagnosis and personalised therapies in COPD and pre-COPD.

Incidence of Interstitial Lung Abnormalities: The MESA Lung Study.

The incidence of newly developed interstitial lung abnormalities (ILA) and fibrotic ILA have not been previously reported.Trained thoracic radiologists evaluated 13 944 cardiac CT scans for the presence of ILA in 6197 Multi-Ethnic Study of Atherosclerosis longitudinal cohort study participants >45 years of age from 2000 to 2012. 5% of the scans were re-read by the same or a different observer in a blinded fashion. After exclusion of participants with ILA at baseline, incidence rates and incidence rate ratios for ILA and fibrotic ILA were calculated.The intra-reader agreement of ILA was 92.0% (Gwet AC1=0.912, ICC=0.982) and the inter-reader agreement of ILA was 83.5% (Gwet AC1=0.814; ICC=0.969). Incidence of ILA and fibrotic ILA was estimated to be 13.1 cases/1000 person-years and 3.5/1000 person-years, respectively. In multivariable analyses, age (HR 1.06 (1.05, 1.08), p <0.001; HR 1.08 (1.06, 1.11), p <0.001), high attenuation area (HAA) at baseline (HR 1.05 (1.03, 1.07), p <0.001; HR 1.06 (1.02, 1.10), p=0.002), and the MUC5B promoter SNP (HR 1.73 (1.17, 2.56) p=0.01; HR 4.96 (2.68, 9.15), p <0.001) were associated with incident ILA and fibrotic ILA, respectively. Ever smoking (HR 2.31 (1.34, 3.96), p= 0.002) and an IPF polygenic risk score (HR 2.09 (1.61-2.71), p<0.001) were associated only with incident fibrotic ILA.Incident ILA and fibrotic ILA were estimated by review of cardiac imaging studies. These findings may lead to wider application of a screening tool for atherosclerosis to identify preclinical lung disease.

Genetic Regulation of Atherosclerosis-Relevant Phenotypes in Human Vascular Smooth Muscle Cells.

RATIONALE: Coronary artery disease (CAD) is a major cause of morbidity and mortality worldwide. Recent genome-wide association studies revealed 163 loci associated with CAD. However, the precise molecular mechanisms by which the majority of these loci increase CAD risk are not known. Vascular smooth muscle cells (VSMCs) are critical in the development of CAD. They can play either beneficial or detrimental roles in lesion pathogenesis, depending on the nature of their phenotypic changes. OBJECTIVE: To identify genetic variants associated with atherosclerosis-relevant phenotypes in VSMCs. METHODS AND RESULTS: We quantified 12 atherosclerosis-relevant phenotypes related to calcification, proliferation, and migration in VSMCs isolated from 151 multiethnic heart transplant donors. After genotyping and imputation, we performed association mapping using 6.3 million genetic variants. We demonstrated significant variations in calcification, proliferation, and migration. These phenotypes were not correlated with each other. We performed genome-wide association studies for 12 atherosclerosis-relevant phenotypes and identified 4 genome-wide significant loci associated with at least one VSMC phenotype. We overlapped the previously identified CAD loci with our data set and found nominally significant associations at 79 loci. One of them was the chromosome 1q41 locus, which harbors MIA3. The G allele of the lead risk single nucleotide polymorphism (SNP) rs67180937 was associated with lower VSMC MIA3 expression and lower proliferation. Lentivirus-mediated silencing of MIA3 (melanoma inhibitory activity protein 3) in VSMCs resulted in lower proliferation, consistent with human genetics findings. Furthermore, we observed a significant reduction of MIA3 protein in VSMCs in thin fibrous caps of late-stage atherosclerotic plaques compared to early fibroatheroma with thick and protective fibrous caps in mice and humans. CONCLUSIONS: Our data demonstrate that genetic variants have significant influences on VSMC function relevant to the development of atherosclerosis. Furthermore, high MIA3 expression may promote atheroprotective VSMC phenotypic transitions, including increased proliferation, which is essential in the formation or maintenance of a protective fibrous cap.

Characterisation of gas exchange in COPD with dissolved-phase hyperpolarised xenon-129 MRI.

To investigate whether hyperpolarised xenon-129 MRI (HXeMRI) enables regional and physiological resolution of diffusing capacity limitations in chronic obstructive pulmonary disease (COPD), we evaluated 34 COPD subjects and 11 healthy volunteers. We report significant correlations between airflow abnormality quantified by HXeMRI and per cent predicted forced expiratory volume in 1 s; HXeMRI gas transfer capacity to red blood cells and carbon monoxide diffusion capacity (%DLCO); and HXeMRI gas transfer capacity to interstitium and per cent emphysema quantified by multidetector chest CT. We further demonstrate the capability of HXeMRI to distinguish varying pathology underlying COPD in subjects with low %DLCO and minimal emphysema.

Genome-Wide Analysis of Left Ventricular Image-Derived Phenotypes Identifies Fourteen Loci Associated With Cardiac Morphogenesis and Heart Failure Development.

BACKGROUND: The genetic basis of left ventricular (LV) image-derived phenotypes, which play a vital role in the diagnosis, management, and risk stratification of cardiovascular diseases, is unclear at present. METHODS: The LV parameters were measured from the cardiovascular magnetic resonance studies of the UK Biobank. Genotyping was done using Affymetrix arrays, augmented by imputation. We performed genome-wide association studies of 6 LV traits-LV end-diastolic volume, LV end-systolic volume, LV stroke volume, LV ejection fraction, LV mass, and LV mass to end-diastolic volume ratio. The replication analysis was performed in the MESA study (Multi-Ethnic Study of Atherosclerosis). We identified the candidate genes at genome-wide significant loci based on the evidence from extensive bioinformatic analyses. Polygenic risk scores were constructed from the summary statistics of LV genome-wide association studies to predict the heart failure events. RESULTS: The study comprised 16 923 European UK Biobank participants (mean age 62.5 years; 45.8% men) without prevalent myocardial infarction or heart failure. We discovered 14 genome-wide significant loci (3 loci each for LV end-diastolic volume, LV end-systolic volume, and LV mass to end-diastolic volume ratio; 4 loci for LV ejection fraction, and 1 locus for LV mass) at a stringent P<1×10-8. Three loci were replicated at Bonferroni significance and 7 loci at nominal significance (P<0.05 with concordant direction of effect) in the MESA study (n=4383). Follow-up bioinformatic analyses identified 28 candidate genes that were enriched in the cardiac developmental pathways and regulation of the LV contractile mechanism. Eight genes (TTN, BAG3, GRK5, HSPB7, MTSS1, ALPK3, NMB, and MMP11) supported by at least 2 independent lines of in silico evidence were implicated in the cardiac morphogenesis and heart failure development. The polygenic risk scores of LV phenotypes were predictive of heart failure in a holdout UK Biobank sample of 3106 cases and 224 134 controls (odds ratio 1.41, 95% CI 1.26 - 1.58, for the top quintile versus the bottom quintile of the LV end-systolic volume risk score). CONCLUSIONS: We report 14 genetic loci and indicate several candidate genes that not only enhance our understanding of the genetic architecture of prognostically important LV phenotypes but also shed light on potential novel therapeutic targets for LV remodeling.

A Genetic Risk Score Associated with Chronic Obstructive Pulmonary Disease Susceptibility and Lung Structure on Computed Tomography.

Rationale: Chronic obstructive pulmonary disease (COPD) has been associated with numerous genetic variants, yet the extent to which its genetic risk is mediated by variation in lung structure remains unknown.Objectives: To characterize associations between a genetic risk score (GRS) associated with COPD susceptibility and lung structure on computed tomography (CT).Methods: We analyzed data from MESA Lung (Multi-Ethnic Study of Atherosclerosis Lung Study), a U.S. general population-based cohort, and SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study). A weighted GRS was calculated from 83 SNPs that were previously associated with lung function. Lung density, spatially matched airway dimensions, and airway counts were assessed on full-lung CT. Generalized linear models were adjusted for age, age squared, sex, height, principal components of genetic ancestry, smoking status, pack-years, CT model, milliamperes, and total lung volume.Measurements and Main Results: MESA Lung and SPIROMICS contributed 2,517 and 2,339 participants, respectively. Higher GRS was associated with lower lung function and increased COPD risk, as well as lower lung density, smaller airway lumens, and fewer small airways, without effect modification by smoking. Adjustment for CT lung structure, particularly small airway measures, attenuated associations between the GRS and FEV1/FVC by 100% and 60% in MESA and SPIROMICS, respectively. Lung structure (P < 0.0001), but not the GRS (P > 0.10), improved discrimination of moderate-to-severe COPD cases relative to clinical factors alone.Conclusions: A GRS associated with COPD susceptibility was associated with CT lung structure. Lung structure may be an important mediator of heritability and determinant of personalized COPD risk.

Omega-3 Fatty Acids and Genome-Wide Interaction Analyses Reveal DPP10-Pulmonary Function Association.

RATIONALE: Omega-3 polyunsaturated fatty acids (n-3 PUFAs) have anti-inflammatory properties that could benefit adults with comprised pulmonary health. OBJECTIVE: To investigate n-3 PUFA associations with spirometric measures of pulmonary function tests (PFTs) and determine underlying genetic susceptibility. METHODS: Associations of n-3 PUFA biomarkers (α-linolenic acid, eicosapentaenoic acid, docosapentaenoic acid [DPA], and docosahexaenoic acid [DHA]) were evaluated with PFTs (FEV1, FVC, and FEV1/FVC) in meta-analyses across seven cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (N = 16,134 of European or African ancestry). PFT-associated n-3 PUFAs were carried forward to genome-wide interaction analyses in the four largest cohorts (N = 11,962) and replicated in one cohort (N = 1,687). Cohort-specific results were combined using joint 2 degree-of-freedom (2df) meta-analyses of SNP associations and their interactions with n-3 PUFAs. RESULTS: DPA and DHA were positively associated with FEV1 and FVC (P < 0.025), with evidence for effect modification by smoking and by sex. Genome-wide analyses identified a novel association of rs11693320-an intronic DPP10 SNP-with FVC when incorporating an interaction with DHA, and the finding was replicated (P2df = 9.4 × 10-9 across discovery and replication cohorts). The rs11693320-A allele (frequency, ∼80%) was associated with lower FVC (PSNP = 2.1 × 10-9; ßSNP = -161.0 ml), and the association was attenuated by higher DHA levels (PSNP×DHA interaction = 2.1 × 10-7; ßSNP×DHA interaction = 36.2 ml). CONCLUSIONS: We corroborated beneficial effects of n-3 PUFAs on pulmonary function. By modeling genome-wide n-3 PUFA interactions, we identified a novel DPP10 SNP association with FVC that was not detectable in much larger studies ignoring this interaction.

D-Dimer in African Americans: Whole Genome Sequence Analysis and Relationship to Cardiovascular Disease Risk in the Jackson Heart Study.

OBJECTIVE: Plasma levels of the fibrinogen degradation product D-dimer are higher among African Americans (AAs) compared with those of European ancestry and higher among women compared with men. Among AAs, little is known of the genetic architecture of D-dimer or the relationship of D-dimer to incident cardiovascular disease. APPROACH AND RESULTS: We measured baseline D-dimer in 4163 AAs aged 21 to 93 years from the prospective JHS (Jackson Heart Study) cohort and assessed association with incident cardiovascular disease events. In participants with whole genome sequencing data (n=2980), we evaluated common and rare genetic variants for association with D-dimer. Each standard deviation higher baseline D-dimer was associated with a 20% to 30% increased hazard for incident coronary heart disease, stroke, and all-cause mortality. Genetic variation near F3 was associated with higher D-dimer (rs2022030, ß=0.284, P=3.24×10-11). The rs2022030 effect size was nearly 3× larger among women (ß=0.373, P=9.06×10-13) than among men (ß=0.135, P=0.06; P interaction =0.009). The sex by rs2022030 interaction was replicated in an independent sample of 10 808 multiethnic men and women (P interaction =0.001). Finally, the African ancestral sickle cell variant (HBB rs334) was significantly associated with higher D-dimer in JHS (ß=0.507, P=1.41×10-14), and this association was successfully replicated in 1933 AAs (P=2.3×10-5). CONCLUSIONS: These results highlight D-dimer as an important predictor of cardiovascular disease risk in AAs and suggest that sex-specific and African ancestral genetic effects of the F3 and HBB loci contribute to the higher levels of D-dimer among women and AAs.

Meta-analysis across Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium provides evidence for an association of serum vitamin D with pulmonary function.

The role that vitamin D plays in pulmonary function remains uncertain. Epidemiological studies reported mixed findings for serum 25-hydroxyvitamin D (25(OH)D)-pulmonary function association. We conducted the largest cross-sectional meta-analysis of the 25(OH)D-pulmonary function association to date, based on nine European ancestry (EA) cohorts (n 22 838) and five African ancestry (AA) cohorts (n 4290) in the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium. Data were analysed using linear models by cohort and ancestry. Effect modification by smoking status (current/former/never) was tested. Results were combined using fixed-effects meta-analysis. Mean serum 25(OH)D was 68 (sd 29) nmol/l for EA and 49 (sd 21) nmol/l for AA. For each 1 nmol/l higher 25(OH)D, forced expiratory volume in the 1st second (FEV1) was higher by 1·1 ml in EA (95 % CI 0·9, 1·3; P<0·0001) and 1·8 ml (95 % CI 1·1, 2·5; P<0·0001) in AA (P race difference=0·06), and forced vital capacity (FVC) was higher by 1·3 ml in EA (95 % CI 1·0, 1·6; P<0·0001) and 1·5 ml (95 % CI 0·8, 2·3; P=0·0001) in AA (P race difference=0·56). Among EA, the 25(OH)D-FVC association was stronger in smokers: per 1 nmol/l higher 25(OH)D, FVC was higher by 1·7 ml (95 % CI 1·1, 2·3) for current smokers and 1·7 ml (95 % CI 1·2, 2·1) for former smokers, compared with 0·8 ml (95 % CI 0·4, 1·2) for never smokers. In summary, the 25(OH)D associations with FEV1 and FVC were positive in both ancestries. In EA, a stronger association was observed for smokers compared with never smokers, which supports the importance of vitamin D in vulnerable populations.

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.

Proteomic Networks and Related Genetic Variants Associated with Smoking and Chronic Obstructive Pulmonary Disease.

Background: Studies have identified individual blood biomarkers associated with chronic obstructive pulmonary disease (COPD) and related phenotypes. However, complex diseases such as COPD typically involve changes in multiple molecules with interconnections that may not be captured when considering single molecular features. Methods: Leveraging proteomic data from 3,173 COPDGene Non-Hispanic White (NHW) and African American (AA) participants, we applied sparse multiple canonical correlation network analysis (SmCCNet) to 4,776 proteins assayed on the SomaScan v4.0 platform to derive sparse networks of proteins associated with current vs. former smoking status, airflow obstruction, and emphysema quantitated from high-resolution computed tomography scans. We then used NetSHy, a dimension reduction technique leveraging network topology, to produce summary scores of each proteomic network, referred to as NetSHy scores. We next performed genome-wide association study (GWAS) to identify variants associated with the NetSHy scores, or network quantitative trait loci (nQTLs). Finally, we evaluated the replicability of the networks in an independent cohort, SPIROMICS. Results: We identified networks of 13 to 104 proteins for each phenotype and exposure in NHW and AA, and the derived NetSHy scores significantly associated with the variable of interests. Networks included known (sRAGE, ALPP, MIP1) and novel molecules (CA10, CPB1, HIS3, PXDN) and interactions involved in COPD pathogenesis. We observed 7 nQTL loci associated with NetSHy scores, 4 of which remained after conditional analysis. Networks for smoking status and emphysema, but not airflow obstruction, demonstrated a high degree of replicability across race groups and cohorts. Conclusions: In this work, we apply state-of-the-art molecular network generation and summarization approaches to proteomic data from COPDGene participants to uncover protein networks associated with COPD phenotypes. We further identify genetic associations with networks. This work discovers protein networks containing known and novel proteins and protein interactions associated with clinically relevant COPD phenotypes across race groups and cohorts.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

The association between aging-related monocyte transcriptional networks and comorbidity burden: the Multi-Ethnic Study of Atherosclerosis (MESA).

Translating our knowledge of the biological aging from animal models to humans may give rise to novel approaches of targeting multiple aging-related diseases simultaneously and increasing health span. Here, for the first time, we use transcriptomic signatures of monocytes to identify biological aging pathways underlying multiple aging-related diseases in humans. The ordinal logistic regression was used to cross-sectionally investigate transcriptomics of the comorbidity index in 1264 community-based Multi-Ethnic Study of Atherosclerosis (MESA) adults, 47% Caucasian, 32% Hispanic, 21% African American, and 51% female, aged 55-94 years. The comorbidity index was defined as the number of prevalent aging-related diseases including cardiovascular disease, type-2 diabetes, hypertension, cancer, dementia, chronic kidney disease, chronic obstructive pulmonary disease, and hip fracture. We identified 708 gene transcripts associated with the comorbidity index (FDR < 0.05) after adjusting for age, sex, ethnicity, and study site. In a weighted gene co-expression network analysis, as postulated, aging-related declines in apoptosis/autophagy (OR = 1.21 per SD increment, p = 0.0006) and ribosome/mitochondrion (OR = 0.90 per SD increment, p = 0.05) were positively associated with the comorbidity index. After adjusting for multiple comparisons, we identified 10 comorbidity-associated modules (FDR < 0.05), including the module of apoptosis/autophagy. There were three inter-correlated modules of these 10 involved in the complement subcomponent C1q, Fc gamma receptor I, and Fc gamma receptor III of the immune system, respectively. Aging-related upregulation of these three modules was positively associated with the comorbidity index. The odds of comorbidity increased with more of these modules acting together in a dose-response fashion. In conclusion, transcriptomic analysis of human immune cells may identify biomarker panels indicative of comprehensive biological mechanisms, especially immune signaling pathways, contributing to health aging.

Large scale proteomic studies create novel privacy considerations.

Privacy protection is a core principle of genomic but not proteomic research. We identified independent single nucleotide polymorphism (SNP) quantitative trait loci (pQTL) from COPDGene and Jackson Heart Study (JHS), calculated continuous protein level genotype probabilities, and then applied a naïve Bayesian approach to link SomaScan 1.3K proteomes to genomes for 2812 independent subjects from COPDGene, JHS, SubPopulations and InteRmediate Outcome Measures In COPD Study (SPIROMICS) and Multi-Ethnic Study of Atherosclerosis (MESA). We correctly linked 90-95% of proteomes to their correct genome and for 95-99% we identify the 1% most likely links. The linking accuracy in subjects with African ancestry was lower (~ 60%) unless training included diverse subjects. With larger profiling (SomaScan 5K) in the Atherosclerosis Risk Communities (ARIC) correct identification was > 99% even in mixed ancestry populations. We also linked proteomes-to-proteomes and used the proteome only to determine features such as sex, ancestry, and first-degree relatives. When serial proteomes are available, the linking algorithm can be used to identify and correct mislabeled samples. This work also demonstrates the importance of including diverse populations in omics research and that large proteomic datasets (> 1000 proteins) can be accurately linked to a specific genome through pQTL knowledge and should not be considered unidentifiable.

Protein-metabolite association studies identify novel proteomic determinants of metabolite levels in human plasma.

Although many novel gene-metabolite and gene-protein associations have been identified using high-throughput biochemical profiling, systematic studies that leverage human genetics to illuminate causal relationships between circulating proteins and metabolites are lacking. Here, we performed protein-metabolite association studies in 3,626 plasma samples from three human cohorts. We detected 171,800 significant protein-metabolite pairwise correlations between 1,265 proteins and 365 metabolites, including established relationships in metabolic and signaling pathways such as the protein thyroxine-binding globulin and the metabolite thyroxine, as well as thousands of new findings. In Mendelian randomization (MR) analyses, we identified putative causal protein-to-metabolite associations. We experimentally validated top MR associations in proof-of-concept plasma metabolomics studies in three murine knockout strains of key protein regulators. These analyses identified previously unrecognized associations between bioactive proteins and metabolites in human plasma. We provide publicly available data to be leveraged for studies in human metabolism and disease.

Correlations between complex human phenotypes vary by genetic background, gender, and environment.

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.