Whole genome sequence association with E-selectin levels reveals loss-of-function variant in African Americans

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 (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). 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.

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.

Inflammatory proteins are associated with mortality in a middle-aged diverse cohort

BACKGROUND

Recent data indicate a decline in overall longevity in the United States. Even prior to the COVID-19 pandemic, an increase in midlife mortality rates had been reported. Life expectancy disparities have persisted in the United States for racial and ethnic groups and for individuals living at low socioeconomic status. These continued trends in mortality indicate the importance of examining biomarkers of mortality at midlife in at-risk populations. Circulating levels of cytokines and inflammatory markers reflect systemic chronic inflammation, which is a well-known driver of many age-related diseases.

METHODS

In this study, we examined the relationship of nine different inflammatory proteins with mortality in a middle-aged socioeconomically diverse cohort of African-American and White men and women (n = 1122; mean age = 47.8 years).

RESULTS

We found significant differences in inflammatory-related protein serum levels between African-American and White middle-aged adults. E-selectin and fibrinogen were significantly higher in African-American adults. IFN-γ, TNF-α trimer, monocyte chemoattractant protein-1 (MCP-1), soluble receptor for advanced glycation end-products (sRAGE) and P-selectin were significantly higher in White participants compared to African-American participants. Higher levels of E-selectin, MCP-1 and P-selectin were associated with a higher mortality risk. Furthermore, there was a significant interaction between sex and IL-6 with mortality. IL-6 levels were associated with an increased risk of mortality, an association that was significantly greater in women than men. In addition, White participants with high levels of sRAGE had significantly higher survival probability than White participants with low levels of sRAGE, while African-American participants had similar survival probabilities across sRAGE levels.

CONCLUSIONS

These results suggest that circulating inflammatory markers can be utilized as indicators of midlife mortality risk in a socioeconomically diverse cohort of African-American and White individuals.

Multimodal measurement of glycocalyx degradation during coronary artery bypass grafting

Background

Glycocalyx shedding and subsequent endothelial dysfunction occur in many conditions, such as in sepsis, in critical illness, and during major surgery such as in coronary artery bypass grafting (CABG) where it has been shown to associate with organ dysfunction. Hitherto, there is no consensus about the golden standard in measuring glycocalyx properties in humans. The objective of this study was to compare different indices of glycocalyx shedding and dysfunction. To this end, we studied patients undergoing elective CABG surgery, which is a known cause of glycocalyx shedding.

Materials and methods

Sublingual glycocalyx thickness was measured in 23 patients by: 1) determining the perfused boundary region (PBR)-an inverse measure of glycocalyx thickness-by means of sidestream dark field imaging technique. This is stated double, 2) measuring plasma levels of the glycocalyx shedding products syndecan-1, hyaluronan, and heparan sulfate and 3) measuring plasma markers of impaired glycocalyx function and endothelial activation (Ang-2, Tie-2, E-selectin, and thrombomodulin). Measurements were performed directly after induction, directly after onset of cardiopulmonary bypass (CPB), and directly after cessation of CPB. We assessed changes over time as well as correlations between the various markers.

Results

The PBR increased from 1.81 ± 0.21 μm after induction of anesthesia to 2.27 ± 0.25 μm (p < 0.0001) directly after CPB was initiated and did not change further during CPB. A similar pattern was seen for syndecan-1, hyaluronan, heparan sulfate, Ang-2, Tie-2, and thrombomodulin. E-selectin levels also increased between induction and the start of CPB and increased further during CPB. The PBR correlated moderately with heparan sulfate, E-selectin, and thrombomodulin and weakly with Syndecan-1, hyaluronan, and Tie-2. Shedding markers syndecan-1 and hyaluronan correlated with all functional markers. Shedding marker heparan sulfate only correlated with Tie-2, thrombomodulin, and E-selectin. Thrombomodulin correlated with all shedding markers.

Conclusion

Our results show that glycocalyx thinning, illustrated by increased sublingual PBR and increased levels of shedding markers, is paralleled with impaired glycocalyx function and increased endothelial activation in CABG surgery with CPB. As correlations between different markers were limited, no single marker could be identified to represent the glycocalyx in its full complexity.

Genetic variants associated with circulating liver injury markers in Mexican Americans, a population at risk for non-alcoholic fatty liver disease

Objective: Mexican Americans are disproportionally affected by non-alcoholic fatty liver disease (NAFLD), liver fibrosis and hepatocellular carcinoma. Noninvasive means to identify those in this population at high risk for these diseases are urgently needed. Approach: The Cameron County Hispanic Cohort (CCHC) is a population-based cohort with high rates of obesity (51%), type 2 diabetes (28%) and NAFLD (49%). In a subgroup of 564 CCHC subjects, we evaluated 339 genetic variants previously reported to be associated with liver injury markers aspartate aminotransferase (AST) and alanine aminotransferase (ALT) in United Kingdom and Japanese cohorts. Results: Association was confirmed for 86 variants. Among them, 27 had higher effect allele frequency in the CCHC than in the United Kingdom and Japanese cohorts, and 16 had stronger associations with AST and ALT than rs738409 (PNPLA3). These included rs17710008 (MYCT1), rs2519093 (ABO), rs1801690 (APOH), rs10409243 (S1PR2), rs1800759 (LOC100507053) and rs2491441 (RGL1), which were also associated with steatosis and/or liver fibrosis measured by vibration-controlled transient elastography. Main contributors to advanced fibrosis risk were rs11240351 (CNTN2), rs1800759 (LOC100507053), rs738409 (PNPLA3) and rs1801690 (APOH), with advanced fibrosis detected in 37.5% of subjects with 3 of these 4 variants [AOR = 11.6 (95% CI) = 3.8-35.3]. AST- and ALT-associated variants implicated distinct pathways (ethanol and galactose degradation versus antigen presentation and B cell development). Finally, 8 variants, including rs62292950 (DNAJC13), were associated with gut microbiome changes. Conclusion: These genotype-phenotype findings may have utility in risk modeling and disease prevention in this high-risk population.

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.

Soluble Urokinase Plasminogen Activator Receptor: Genetic Variation and Cardiovascular Disease Risk in Black Adults

BACKGROUND

suPAR (Soluble urokinase plasminogen activator receptor) has emerged as an important biomarker of coagulation, inflammation, and cardiovascular disease (CVD) risk. The contribution of suPAR to CVD risk and its genetic influence in Black populations have not been evaluated.

METHODS

We measured suPAR in 3492 Black adults from the prospective, community-based JHS (Jackson Heart Study). Cross-sectional associations of suPAR with lifestyle and CVD risk factors were assessed, whole-genome sequence data were used to evaluate genetic associations of suPAR, and relationships of suPAR with incident CVD outcomes and overall mortality were estimated over follow-up.

RESULTS

In Cox models adjusted for traditional CVD risk factors, estimated glomerular filtration rate, and CRP (C-reactive protein), each 1-SD higher suPAR was associated with a 21% to 31% increased risk of incident coronary heart disease, heart failure, stroke, and mortality. In the genome-wide association study, 2 missense (rs399145 encoding p.Thr86Ala, rs4760 encoding p.Phe272Leu) and 2 noncoding regulatory variants (rs73935023 within an enhancer element and rs4251805 within the promoter) of PLAUR on chromosome 19 were each independently associated with suPAR and together explained 14% of suPAR phenotypic variation. The allele frequencies of each of the four suPAR-associated genetic variants differ considerably across African and European populations. We further show that PLAUR rs73935023 can alter transcriptional activity in vitro. We did not find any association between genetically determined suPAR and CVD in JHS or a larger electronic medical record-based analyses of Blacks or Whites.

CONCLUSIONS

Our results demonstrate the importance of ancestry-differentiated genetic variation on suPAR levels and indicate suPAR is a CVD biomarker in Black adults.

Genome-wide association study revealed novel candidate gene loci associated with soluble E-selectin levels in a Taiwanese population

BACKGROUND AND AIMS

Increase soluble E-selectin (sE-selectin) levels are associated with various inflammation and cardiometabolic disorders.

METHODS

This study aimed to investigate the genetic determinants of circulating sE-selectin levels by genome-wide association study (GWAS) in 4,525 Taiwan Biobank (TWB) participants and genotype-phenotype association analysis for sE-selectin level-determining alleles in over 80,000 TWB participants.

RESULTS

By GWAS, ABO, SELE, and FUT6 gene variants were identified as the determinants of sE-selectin levels, which reach genome-wide significance (maximum p = 3.25 × 10-271, 4.81 × 10-14, and 9.64 × 10-12, respectively). After further adjustment for the lead ABO rs2519093 genotypes, three novel gene loci, EVI5, FER and DMAC1, were associated with sE-selectin levels at p < 5 × 10-7. Three other previously reported gene loci, CELSR2, ST3GAL6-AS1, and HNF1A-AS1, also showed supportive evidence for the association with sE-selectin levels (maximum p < 0.0073). A multivariate analysis revealed age, body mass index, current smoking, hemoglobin A1C, hematocrit, leukocyte and platelet counts, serum alanine aminotransferase, triglycerides, and uric acid levels were independently associated with sE-selectin levels, in which the above ten gene loci contribute to 27.68% of the variance. For genotype-phenotype association analysis, a pleiotropic effect was demonstrated with genome-wide significant association between ABO gene variants and total and low-density-lipoprotein cholesterol levels, leukocyte counts and hematocrit.

CONCLUSIONS

Our data provide novel insight into the regulation of sE-selectin levels. These results may open new avenues in understanding the critical role of E-selectin on the pathogenesis of inflammatory and cardiometabolic disorders.

What the Jackson Heart Study Has Taught Us About Diabetes and Cardiovascular Disease in the African American Community: a 20-year Appreciation

PURPOSE OF REVIEW

The burden of cardiometabolic diseases such as cardiovascular disease (CVD) and type 2 diabetes (T2D) is pronounced among African Americans. Research has shown that behavioral, social, metabolic, psychosocial, and genetic risk factors of CVD and T2D are closely interwoven. Approximately 20 years ago, the Jackson Heart Study (JHS) was established to investigate this constellation of risk factors.

RECENT FINDINGS

Findings from neighborhood studies emphasize the importance of social cohesion and physical environment in the context CVD and T2D risk. Socioeconomic status factors such as income and education were significant predictors for CVD and T2D. Behavioral studies indicate that modifiable risk factors such as smoking, physical inactivity, lack of sleep, and poor nutrition are associated with CVD risk and all-cause mortality. Mental health also was found to be associated with CVD and T2D. Genetic influences are associated with disease etiology. This review summarizes the joint contributions of CVD and cardiometabolic risk factors in an African American population.

Association between ABO and Duffy blood types and circulating chemokines and cytokines

Blood group antigens are inherited traits that may play a role in immune and inflammatory processes. We investigated associations between blood groups and circulating inflammation-related molecules in 3537 non-Hispanic white participants selected from the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Whole-genome scans were used to infer blood types for 12 common antigen systems based on well-characterized single-nucleotide polymorphisms. Serum levels of 96 biomarkers were measured on multiplex fluorescent bead-based panels. We estimated marker associations with blood type using weighted linear or logistic regression models adjusted for age, sex, smoking status, and principal components of population substructure. Bonferroni correction was used to control for multiple comparisons, with two-sided p values < 0.05 considered statistically significant. Among the 1152 associations tested, 10 were statistically significant. Duffy blood type was associated with levels of CXCL6/GCP2, CXCL5/ENA78, CCL11/EOTAXIN, CXCL1/GRO, CCL2/MCP1, CCL13/MCP4, and CCL17/TARC, whereas ABO blood type was associated with levels of sVEGFR2, sVEGFR3, and sGP130. Post hoc pairwise t-tests showed that individuals with type Fy(a+b−) had the lowest mean levels of all Duffy-associated markers, while individuals with type A blood had the lowest mean levels of all ABO-associated markers. Additional work is warranted to explore potential clinical implications of these differences.

PSCAN: Spatial scan tests guided by protein structures improve complex disease gene discovery and signal variant detection

Germline disease-causing variants are generally more spatially clustered in protein 3-dimensional structures than benign variants. Motivated by this tendency, we develop a fast and powerful protein-structure-based scan (PSCAN) approach for evaluating gene-level associations with complex disease and detecting signal variants. We validate PSCAN's performance on synthetic data and two real data sets for lipid traits and Alzheimer's disease. Our results demonstrate that PSCAN performs competitively with existing gene-level tests while increasing power and identifying more specific signal variant sets. Furthermore, PSCAN enables generation of hypotheses about the molecular basis for the associations in the context of protein structures and functional domains.

ABO O gene frequency increase in the US might be causing increased maternal mortality

Maternal mortality rate has increased in the United States over the past 30 years from 16 deaths per 100,000 births to 28 deaths per 100,000 births while the rest of the world is experiencing declining rates. Increasing obesity and c-section rates in the US have been cited as contributing factors needing remediation, and because of the two to three fold difference in maternal mortality rates in non-Hispanic black women compared to white women, inequality and implicit racial bias has been targeted as well for remediation. Using an epidemiologic approach, a hypothesis here brought to bear is that US immigration policy changes over the past 50 years have brought changes in the gene pool that have caused increasing obstetric hemorrhage and other causes of maternal death. ABO gene frequencies have changed in the US during this time such that ABO O, a gene associated with hemorrhage and mortality in pregnancy, has increased in frequency in the US thus increasing population maternal mortality rate. Using mendelian randomization logic, noting the increase in ABO O gene in the US population over the past 30 years and the association of ABO O gene with both hemorrhage and lower longevity, the increase in frequency of the ABO O gene in the past 30 years in the US population might be causative of an increase in maternal mortality rate. Consequences of this hypothesis would include recognition of the role of ABO gene and thus ABO blood group in prediction of risk of obstetric hemorrhage. Thus those at risk on this basis would be under high surveillance and would have medications and treatment strategies readily available. While research on ABO gene and pregnancy has been done, much of the research is being done in countries other than the US, and given the increasing mortality in the US as well as the role that ABO gene may have in that, further research needs to be done in US populations to quantify risk for all adverse events in pregnancy related to ABO blood type including hemorrhage as well as inter-related causes including pre-eclampsia, cardiovascular disease, thromboembolic disease and infection.

Comparison of Proteomic Assessment Methods in Multiple Cohort Studies

Novel proteomics platforms, such as the aptamer-based SOMAscan platform, can quantify large numbers of proteins efficiently and cost-effectively and are rapidly growing in popularity. However, comparisons to conventional immunoassays remain underexplored, leaving investigators unsure when cross-assay comparisons are appropriate. The correlation of results from immunoassays with relative protein quantification is explored by SOMAscan. For 63 proteins assessed in two chronic obstructive pulmonary disease (COPD) cohorts, subpopulations and intermediate outcome measures in COPD Study (SPIROMICS), and COPDGene, using myriad rules based medicine multiplex immunoassays and SOMAscan, Spearman correlation coefficients range from -0.13 to 0.97, with a median correlation coefficient of ≈0.5 and consistent results across cohorts. A similar range is observed for immunoassays in the population-based Multi-Ethnic Study of Atherosclerosis and for other assays in COPDGene and SPIROMICS. Comparisons of relative quantification from the antibody-based Olink platform and SOMAscan in a small cohort of myocardial infarction patients also show a wide correlation range. Finally, cis pQTL data, mass spectrometry aptamer confirmation, and other publicly available data are integrated to assess relationships with observed correlations. Correlation between proteomics assays shows a wide range and should be carefully considered when comparing and meta-analyzing proteomics data across assays and studies.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.