Association of exome sequences with plasma C-reactive protein levels in >9000 participants

Genetic variant rs1205 is associated with COVID-19 outcomes: The Strong Heart Study and Strong Heart Family Study

BACKGROUND

Although COVID-19 infection has been associated with a number of clinical and environmental risk factors, host genetic variation has also been associated with the incidence and morbidity of infection. The CRP gene codes for a critical component of the innate immune system and CRP variants have been reported associated with infectious disease and vaccination outcomes. We investigated possible associations between COVID-19 outcome and a limited number of candidate gene variants including rs1205.

METHODOLOGY/PRINCIPAL FINDINGS

The Strong Heart and Strong Heart Family studies have accumulated detailed genetic, cardiovascular risk and event data in geographically dispersed American Indian communities since 1988. Genotypic data and 91 COVID-19 adjudicated deaths or hospitalizations from 2/1/20 through 3/1/23 were identified among 3,780 participants in two subsets. Among 21 candidate variants including genes in the interferon response pathway, APOE, TMPRSS2, TLR3, the HLA complex and the ABO blood group, only rs1205, a 3' untranslated region variant in the CRP gene, showed nominally significant association in T-dominant model analyses (odds ratio 1.859, 95%CI 1.001-3.453, p = 0.049) after adjustment for age, sex, center, body mass index, and a history of cardiovascular disease. Within the younger subset, association with the rs1205 T-Dom genotype was stronger, both in the same adjusted logistic model and in the SOLAR analysis also adjusting for other genetic relatedness.

CONCLUSION

A T-dominant genotype of rs1205 in the CRP gene is associated with COVID-19 death or hospitalization, even after adjustment for relevant clinical factors and potential participant relatedness. Additional study of other populations and genetic variants of this gene are warranted.

Rare functional variants in the CRP and G6PC genes modify the relationship between obesity and serum C-reactive protein in white British population

BACKGROUND

C-reactive protein (CRP) is a sensitive biomarker of inflammation with moderate heritability. The role of rare functional genetic variants in relation to serum CRP is understudied. We aimed to examine gene mutation burden of protein-altering (PA) and loss-of-function (LOF) variants in association with serum CRP, and to further explore the clinical relevance.

METHODS

We included 161,430 unrelated participants of European ancestry from the UK Biobank. Of the rare (minor allele frequency <0.1%) and functional variants, 1,776,249 PA and 266,226 LOF variants were identified. Gene-based burden tests, linear regressions, and logistic regressions were performed to identify the candidate mutations at the gene and variant levels, to estimate the potential interaction effect between the identified PA mutation and obesity, and to evaluate the relative risk of 16 CRP-associated diseases.

RESULTS

At the gene level, PA mutation burdens of the CRP (β = -0.685, p = 2.87e-28) and G6PC genes (β = 0.203, p = 1.50e-06) were associated with reduced and increased serum CRP concentration, respectively. At the variant level, seven PA alleles in the CRP gene decreased serum CRP, of which the per-allele effects were approximately three to seven times greater than that of a common variant in the same locus. The effects of obesity and central obesity on serum CRP concentration were smaller among the PA mutation carriers in the CRP (pinteraction  = 0.008) and G6PC gene (pinteraction  = 0.034) compared to the corresponding non-carriers.

CONCLUSION

PA mutation burdens in the CRP and G6PC genes are strongly associated with decreased serum CRP concentrations. As serum CRP and obesity are important predictors of cardiovascular risks in clinics, our observations suggest taking rare genetic factors into consideration might improve the delivery of precision medicine.

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.

Whole-exome sequence analysis of anthropometric traits illustrates challenges in identifying effects of rare genetic variants

Anthropometric traits, measuring body size and shape, are highly heritable and significant clinical risk factors for cardiometabolic disorders. These traits have been extensively studied in genome-wide association studies (GWASs), with hundreds of genome-wide significant loci identified. We performed a whole-exome sequence analysis of the genetics of height, body mass index (BMI) and waist/hip ratio (WHR). We meta-analyzed single-variant and gene-based associations of whole-exome sequence variation with height, BMI, and WHR in up to 22,004 individuals, and we assessed replication of our findings in up to 16,418 individuals from 10 independent cohorts from Trans-Omics for Precision Medicine (TOPMed). We identified four trait associations with single-nucleotide variants (SNVs; two for height and two for BMI) and replicated the LECT2 gene association with height. Our expression quantitative trait locus (eQTL) analysis within previously reported GWAS loci implicated CEP63 and RFT1 as potential functional genes for known height loci. We further assessed enrichment of SNVs, which were monogenic or syndromic variants within loci associated with our three traits. This led to the significant enrichment results for height, whereas we observed no Bonferroni-corrected significance for all SNVs. With a sample size of ∼20,000 whole-exome sequences in our discovery dataset, our findings demonstrate the importance of genomic sequencing in genetic association studies, yet they also illustrate the challenges in identifying effects of rare genetic variants.

Molecular Biology Networks and Key Gene Regulators for Inflammatory Biomarkers Shared by Breast Cancer Development: Multi-Omics Systems Analysis

As key inflammatory biomarkers C-reactive protein (CRP) and interleukin-6 (IL6) play an important role in the pathogenesis of non-inflammatory diseases, including specific cancers, such as breast cancer (BC). Previous genome-wide association studies (GWASs) have neither explained the large proportion of genetic heritability nor provided comprehensive understanding of the underlying regulatory mechanisms. We adopted an integrative genomic network approach by incorporating our previous GWAS data for CRP and IL6 with multi-omics datasets, such as whole-blood expression quantitative loci, molecular biologic pathways, and gene regulatory networks to capture the full range of genetic functionalities associated with CRP/IL6 and tissue-specific key drivers (KDs) in gene subnetworks. We applied another systematic genomics approach for BC development to detect shared gene sets in enriched subnetworks across BC and CRP/IL6. We detected the topmost significant common pathways across CRP/IL6 (e.g., immune regulatory; chemokines and their receptors; interferon γ, JAK-STAT, and ERBB4 signaling), several of which overlapped with BC pathways. Further, in gene–gene interaction networks enriched by those topmost pathways, we identified KDs—both well-established (e.g., JAK1/2/3, STAT3) and novel (e.g., CXCR3, CD3D, CD3G, STAT6)—in a tissue-specific manner, for mechanisms shared in regulating CRP/IL6 and BC risk. Our study may provide robust, comprehensive insights into the mechanisms of CRP/IL6 regulation and highlight potential novel genetic targets as preventive and therapeutic strategies for associated disorders, such as BC.

Genetically determined elevated C-reactive protein associated with primary colorectal cancer risk: Mendelian randomization with lifestyle interactions

Systemic inflammation-related etiologic pathways via inflammatory cytokines in the development of colorectal cancer (CRC) have not been convincingly determined and may be confounded by lifestyle factors or reverse causality. We investigated the genetically predicted C-reactive protein (CRP) phenotype in the potential causal pathway of primary CRC risk in postmenopausal women in a Mendelian randomization (MR) framework. We employed individual-level data of the Women's Health Initiative Database for Genotypes and Phenotypes Study, which consists of 5 genome-wide association (GWA) studies, including 10,142 women, 737 of whom developed primary CRC. We examined 61 GWA single-nucleotide polymorphisms (SNPs) associated with CRP by using weighted/penalized MR weighted-medians and MR gene-environment interactions that allow some relaxation of the strict variable requirements and attenuate the heterogeneous estimates of outlying SNPs. In lifestyle-stratification analyses, genetically determined CRP exhibited its effects on the decreased CRC risk in non-viscerally obese and high-fat diet subgroups. In contrast, genetically driven CRP was associated with an increased risk for CRC in women who smoked ≥ 15 cigarettes/day, with significant interaction of the gene-smoking relationship. Further, a substantially increased risk of CRC induced by CRP was observed in relatively short-term users (< 5 years) of estrogen (E)-only and also longer-term users (5 to > 10 years) of E plus progestin. Our findings may provide novel evidence on immune-related etiologic pathways connected to CRC risk and suggest the possible use of CRP as a CRC-predictive biomarker in women with particular behaviors and CRP marker-informed interventions to reduce CRC risk.

Genetically Predicted C-Reactive Protein Associated With Postmenopausal Breast Cancer Risk: Interrelation With Estrogen and Cancer Molecular Subtypes Using Mendelian Randomization

Background

Immune-related etiologic pathways that influence breast cancer risk are incompletely understood and may be confounded by lifestyles or reverse causality. Using a Mendelian randomization (MR) approach, we investigated the potential causal relationship between genetically elevated C-reactive protein (CRP) concentrations and primary invasive breast cancer risk in postmenopausal women.

Methods

We used individual-level data obtained from 10,179 women, including 537 who developed breast cancer, from the Women’s Health Initiative Database for Genotypes and Phenotypes Study, which consists of five genome-wide association (GWA) studies. We examined 61 GWA single-nucleotide polymorphisms (SNPs) previously associated with CRP. We employed weighted/penalized weighted–medians and MR gene–environment interactions that allow instruments’ invalidity to some extent and attenuate the heterogeneous estimates of outlying SNPs.

Results

In lifestyle-stratification analyses, genetically elevated CRP decreased risk for breast cancer in exogenous estrogen-only, estrogen + progestin, and past oral contraceptive (OC) users, but only among relatively short-term users (<5 years). Estrogen-only users for ≥5 years had more profound CRP-decreased breast cancer risk in dose–response fashion, whereas past OC users for ≥5 years had CRP-increased cancer risk. Also, genetically predicted CRP was strongly associated with increased risk for hormone-receptor positive or human epidermal growth factor receptor-2 negative breast cancer.

Conclusions

Our findings may provide novel evidence on the immune-related molecular pathways linking to breast cancer risk and suggest potential clinical use of CRP to predict the specific cancer subtypes. Our findings suggest potential interventions targeting CRP–inflammatory markers to reduce breast cancer risk.

Pro-inflammatory cytokine polymorphisms and interactions with dietary alcohol and estrogen, risk factors for invasive breast cancer using a post genome-wide analysis for gene-gene and gene-lifestyle interaction

Molecular and genetic immune-related pathways connected to breast cancer and lifestyles in postmenopausal women are not fully characterized. In this study, we explored the role of pro-inflammatory cytokines such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways at the genome-wide level. With single-nucleotide polymorphisms (SNPs) in the biomarkers and lifestyles together, we further constructed risk profiles to improve predictability for breast cancer. Our earlier genome-wide association gene-environment interaction study used large cohort data from the Women's Health Initiative Database for Genotypes and Phenotypes Study and identified 88 SNPs associated with CRP and IL-6. For this study, we added an additional 68 SNPs from previous GWA studies, and together with 48 selected lifestyles, evaluated for the association with breast cancer risk via a 2-stage multimodal random survival forest and generalized multifactor dimensionality reduction methods. Overall and in obesity strata (by body mass index, waist, waist-to-hip ratio, exercise, and dietary fat intake), we identified the most predictive genetic and lifestyle variables. Two SNPs (SALL1 rs10521222 and HLA-DQA1 rs9271608) and lifestyles, including alcohol intake, lifetime cumulative exposure to estrogen, and overall and visceral obesity, are the most common and strongest predictive markers for breast cancer across the analyses. The risk profile that combined those variables presented their synergistic effect on the increased breast cancer risk in a gene-lifestyle dose-dependent manner. Our study may contribute to improved predictability for breast cancer and suggest potential interventions for the women with the risk genotypes and lifestyles to reduce their breast cancer risk.

Inflammation and cardiovascular disease: From mechanisms to therapeutics

Inflammation constitutes a complex, highly conserved cascade of molecular and cellular events. Inflammation has been labeled as “the fire within,” is highly regulated, and is critical to host defense and tissue repair. In general, inflammation is beneficial and has evolved to promote survival. However, inflammation can also be maladaptive when chronically activated and sustained, leading to progressive tissue injury and reduced survival. Examples of a maladaptive response include rheumatologic disease and atherosclerosis. Despite evidence gathered by Virchow over 100 years ago showing that inflammatory white cells play a role in atherogenesis, atherosclerosis was until recently viewed as a disease of passive cholesterol accumulation in the subendothelial space. This view has been supplanted by considerable basic scientific and clinical evidence demonstrating that every step of atherogenesis, from the development of endothelial cell dysfunction to foam cell formation, plaque formation and progression, and ultimately plaque rupture stemming from architectural instability, is driven by the cytokines, interleukins, and cellular constituents of the inflammatory response. Herein we provide an overview of the role of inflammation in atherosclerotic cardiovascular disease, discuss the predictive value of various biomarkers involved in inflammation, and summarize recent clinical trials that evaluated the capacity of various pharmacologic interventions to attenuate the intensity of inflammation and impact risk for acute cardiovascular events.

Pro-inflammatory cytokine polymorphisms in ONECUT2 and HNF4A and primary colorectal carcinoma: a post genome-wide gene-lifestyle interaction study

Immune-related molecular and genetic pathways that are connected to colorectal cancer (CRC) and lifestyles in postmenopausal women are incompletely characterized. In this study, we examined the role of pro-inflammatory biomarkers such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways. Through selection of the best predictive single-nucleotide polymorphisms (SNPs) and lifestyles, our goal was to improve the prediction accuracy and ability for CRC risk. Using large cohort data of postmenopausal women from the Women's Health Initiative Database for Genotypes and Phenotypes Study, we previously conducted a genome-wide association (GWA) for a CRP and IL-6 gene-behavioral interaction study. For the present study, we added GWA-SNPs from outside GWA studies, resulting in a total of 152 SNPs. Together with 41 selected lifestyles, we performed a 2-stage multimodal random survival forest analysis with generalized multifactor dimensionality reduction approach to construct CRC risk profiles. Overall and in obesity strata (by body mass index, waist circumference, waist-to-hip ratio, exercise, and dietary fat intake), we identified the best predictive genetic markers in inflammatory cytokines and lifestyles. Across the strata, 2 SNPs (ONECUT2 rs4092465 and HNF4A rs1800961) and 1 lifestyle factor (relatively short-term past use of oral contraceptives) were the most common and strongest predictive markers for CRC risk. The risk profile that combined those variables exhibited synergistically increased risk for CRC; this pattern appeared more strongly in obese and inactive subgroups. Our results may contribute to improved predictability for CRC and suggest genetically targeted lifestyle interventions for women carrying the inflammatory-risk genotypes, reducing CRC risk.

Burden of rare exome sequence variants in PROC gene is associated with venous thromboembolism: a population-based study

BACKGROUND

Rare coding mutations underlying deficiencies of antithrombin and proteins C and S contribute to familial venous thromboembolism (VTE). It is uncertain whether rare variants play a role in the etiology of VTE in the general population.

OBJECTIVES

We conducted a deep whole-exome sequencing (WES) study to investigate the associations between rare coding variants and the risk of VTE in two population-based prospective cohorts.

PATIENTS/METHODS

Whole-exome sequencing was performed in the Longitudinal Investigation of Thromboembolism Etiology (LITE), which combines the Atherosclerosis Risk in Communities (ARIC) study (316 incident VTE events among 3159 African Americans [AAs] and 458 incident VTEs among 7772 European Americans [EAs]) and the Cardiovascular Healthy Study (CHS; 60 incident VTEs among 1751 EAs). We performed gene-based tests of rare variants (allele frequency < 1%, exome-wide significance P < 1.47 × 10-6 ) separately in each study and ancestry group, and meta-analyzed the results for the EAs in ARIC and CHS.

RESULTS

In the meta-analysis of EAs, we identified one gene, PROC, in which the burden of rare, coding variants was significantly associated with increased risk of VTE (HR = 5.42 [3.11, 9.42] for carriers versus non-carriers, P = 2.27 × 10-9 ). In ARIC EAs, carriers of the PROC rare variants had on average 0.75 standard deviation (SD) lower concentrations of plasma protein C and 0.28 SD higher D-dimer (P < .05) than non-carriers. Adjustment for low protein C status did not eliminate the association of PROC burden with VTE. In AAs, rare coding PROC variants were not associated with VTE.

CONCLUSIONS

Rare coding variants in PROC contribute to increased VTE risk in EAs in this general population sample.

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.

Genetic analysis of hsCRP in American Indians: The Strong Heart Family Study

Background

Increased serum levels of C-reactive protein (CRP), an important component of the innate immune response, are associated with increased risk of cardiovascular disease (CVD). Multiple single nucleotide polymorphisms (SNP) have been identified which are associated with CRP levels, and Mendelian randomization studies have shown a positive association between SNPs increasing CRP expression and risk of colon cancer (but thus far not CVD). The effects of individual genetic variants often interact with the genetic background of a population and hence we sought to resolve the genetic determinants of serum CRP in a number of American Indian populations.

Methods

The Strong Heart Family Study (SHFS) has serum CRP measurements from 2428 tribal members, recruited as large families from three regions of the United States. Microsatellite markers and MetaboChip defined SNP genotypes were incorporated into variance components, decomposition-based linkage and association analyses.

Results

CRP levels exhibited significant heritability (h² = 0.33 ± 0.05, p<1.3 X 10⁻²⁰). A locus on chromosome (chr) 6, near marker D6S281 (approximately at 169.6 Mb, GRCh38/hg38) showed suggestive linkage (LOD = 1.9) to CRP levels. No individual SNPs were found associated with CRP levels after Bonferroni adjustment for multiple testing (threshold <7.77 x 10⁻⁷), however, we found nominal associations, many of which replicate previous findings at the CRP, HNF1A and 7 other loci. In addition, we report association of 46 SNPs located at 7 novel loci on chromosomes 2, 5, 6(2 loci), 9, 10 and 17, with an average of 15.3 Kb between SNPs and all with p-values less than 7.2 X 10⁻⁴.

Conclusion

In agreement with evidence from other populations, these data show CRP serum levels are under considerable genetic influence; and include loci, such as near CRP and other genes, that replicate results from other ethnic groups. These findings also suggest possible novel loci on chr 6 and other chromosomes that warrant further investigation.

A Rare Variant Nonparametric Linkage Method for Nuclear and Extended Pedigrees with Application to Late-Onset Alzheimer Disease via WGS Data

To analyze family-based whole-genome sequence (WGS) data for complex traits, we developed a rare variant (RV) non-parametric linkage (NPL) analysis method, which has advantages over association methods. The RV-NPL differs from the NPL in that RVs are analyzed, and allele sharing among affected relative-pairs is estimated only for minor alleles. Analyzing families can increase power because causal variants with familial aggregation usually have larger effect sizes than those underlying sporadic diseases. Differing from association analysis, for NPL only affected individuals are analyzed, which can increase power, since unaffected family members can be susceptibility variant carriers. RV-NPL is robust to population substructure and admixture, inclusion of nonpathogenic variants, as well as allelic and locus heterogeneity and can readily be applied outside of coding regions. In contrast to analyzing common variants using NPL, where loci localize to large genomic regions (e.g., >50 Mb), mapped regions are well defined for RV-NPL. Using simulation studies, we demonstrate that RV-NPL is substantially more powerful than applying traditional NPL methods to analyze RVs. The RV-NPL was applied to analyze 107 late-onset Alzheimer disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with WGS data, and statistically significant linkage (LOD ≥ 3.8) was found with RVs in PSMF1 and PTPN21 which have been shown to be involved in LOAD etiology. Additionally, nominally significant linkage was observed with RVs in ABCA7, ACE, EPHA1, and SORL1, genes that were previously reported to be associated with LOAD. RV-NPL is an ideal method to elucidate the genetic etiology of complex familial diseases.

C-reactive protein gene 1846C>T polymorphism is associated with increased risk and clinical features of lung cancer: a case-control study

Chronic inflammation plays an important role in lung carcinogenesis. Recently, several studies investigated the association of C-reactive protein (CRP) gene 1846C>T polymorphism and lung cancer (LC) risk, but with conflicting findings. In the present study, we conducted this case-control study with 408 LC patients and 472 healthy controls in a Chinese Han population. Genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLR) method. Our data found that CRP gene 1846C>T polymorphism increased the risk of LC. Subgroup analyses obtained significant associations among the groups of males, ≥50 years old, smoking, and non-drinkers. Bioinformatics analysis showed that the expression levels of CRP in LC tissues were significantly increased compared with normal tissues. Additionally, the present study found CRP mRNA high expression was associated with worse survival in LC patients. Furthermore, our data indicated that TT genotype of 1846C>T polymorphism was associated with a larger size of tumor and was related with lymphatic metastasis in LC patients. In conclusion, the present study suggests that CRP gene 1846C>T polymorphism is associated with increased risk of LC. CRP gene 1846C>T polymorphism may be a potential marker for the diagnosis of LC.

Low-dose aspirin protective effects are correlated with deregulation of HNF factor expression in the preeclamptic placentas from mice and humans

Aspirin (acetyl-salicylic acid) is one of the most ancient drugs of the human pharmacopeia. Nonetheless, its action at low doses is not well understood at the molecular level. One of the applications of low-dose aspirin treatment is the prevention of preeclampsia (PE) in patients at risk. Foeto-placental overexpression of the STOX1A transcription factor in mice triggers PE symptoms. Transcriptomic analysis of the placentas, showed that aspirin massively down-regulates genes of the coagulation and complement cascade, as well as genes involved in lipid transport. The genes modified by aspirin treatment are not the ones that are modified by STOX1 overexpression, suggesting that aspirin could act downstream, symptomatically on the preeclamptic disease. Bioinformatics analysis of the promoters of the deregulated genes showed that they are strongly enriched in HNF transcription factors-binding sites, in accordance with existing literature showing their roles as regulators of coagulation. Two of these transcription factors, Hnf1β and Hnf4α are found down-regulated by aspirin treatment. In parallel, we show that in human patient placentas, aspirin-induced deregulations of genes of the coagulation cascade are also observed. Finally, the expression of Hnf1β target sequences (Kif12, F2, Hnf4α promoters and a synthetic concatemer of the Hnf1β-binding site) were investigated by transfection in trophoblast cell models, with or without aspirin treatment and with or without STOX1A overexpression. In this model we observed that STOX1A and aspirin tended to synergize in the down-regulation of Hnf1β target genes in trophoblasts.

Decreased plasma C-reactive protein levels in APOE ε4 allele carriers

Objective

Apolipoprotein E (APOE) ε4 allele is a well‐established risk factor in Alzheimer's disease (AD). Here, we assessed the effects of APOE polymorphism on cardiovascular, metabolic, and inflammation‐related parameters in population‐based cohorts.

Methods

Association of cardiovascular, metabolic, and inflammation‐related parameters with the APOE polymorphism in a large Finnish Metabolic Syndrome in Men (METSIM) cohort and Finnish Geriatric Intervention study to prevent cognitive impairment and disability (FINGER) were investigated. Brain‐specific effects were addressed in postmortem brain samples.

Results

Individuals carrying the APOE ε4 allele displayed significantly elevated serum/plasma LDL cholesterol and apolipoprotein B levels. APOE ε3ε4 and ε4ε4 significantly associated with lower levels of plasma high‐sensitivity C‐reactive protein (hs‐CRP). Plasma amyloid‐β 42 (Aβ42) and reduced hs‐CRP levels showed an association independently of the APOE status.

Interpretation

These data suggest that the APOE ε4 allele associates with lower levels of hs‐CRP in individuals without dementia. Moreover, Aβ42 may encompass anti‐inflammatory effects reflected by reduced hs‐CRP levels.

Discovery, fine-mapping, and conditional analyses of genetic variants associated with C-reactive protein in multiethnic populations using the Metabochip in the Population Architecture using Genomics and Epidemiology (PAGE) study

C-reactive protein (CRP) is a circulating biomarker indicative of systemic inflammation. We aimed to evaluate genetic associations with CRP levels among non-European-ancestry populations through discovery, fine-mapping and conditional analyses. A total of 30 503 non-European-ancestry participants from 6 studies participating in the Population Architecture using Genomics and Epidemiology study had serum high-sensitivity CRP measurements and ∼200 000 single nucleotide polymorphisms (SNPs) genotyped on the Metabochip. We evaluated the association between each SNP and log-transformed CRP levels using multivariate linear regression, with additive genetic models adjusted for age, sex, the first four principal components of genetic ancestry, and study-specific factors. Differential linkage disequilibrium patterns between race/ethnicity groups were used to fine-map regions associated with CRP levels. Conditional analyses evaluated for multiple independent signals within genetic regions. One hundred and sixty-three unique variants in 12 loci in overall or race/ethnicity-stratified Metabochip-wide scans reached a Bonferroni-corrected P-value <2.5E-7. Three loci have no (HACL1, OLFML2B) or only limited (PLA2G6) previous associations with CRP levels. Six loci had different top hits in race/ethnicity-specific versus overall analyses. Fine-mapping refined the signal in six loci, particularly in HNF1A. Conditional analyses provided evidence for secondary signals in LEPR, IL1RN and HNF1A, and for multiple independent signals in CRP and APOE. We identified novel variants and loci associated with CRP levels, generalized known CRP associations to a multiethnic study population, refined association signals at several loci and found evidence for multiple independent signals at several well-known loci. This study demonstrates the benefit of conducting inclusive genetic association studies in large multiethnic populations.

Population-based differences in immune system response contribute to an increased risk of schizophrenia in African migrants?

Among the highest incidences of schizophrenia is the one documented in second-generation migrants of African descent in the Western countries. Interestingly, people of African and European ancestry demonstrate significant genetic-based differences in immune system regulation and response. As a result, the pro-inflammatory phenotype is more pronounced in people of African descent than it is in Europeans. At the same time, the role of the immune system in the etiology of schizophrenia is gaining increased recognition. Here, we propose that the population-specific genetic variation within the immune system interacts with unfavourable environments to contribute to a higher risk of schizophrenia in people of African ancestry. Our hypothesis is supported by recent findings from two separate fields of research-population genetics and psychoneuroimmunology. Moreover, we highlight the need to include African populations in genetic studies of schizophrenia.

Cardiovascular Health in African Americans: A Scientific Statement From the American Heart Association

BACKGROUND AND PURPOSE

Population-wide reductions in cardiovascular disease incidence and mortality have not been shared equally by African Americans. The burden of cardiovascular disease in the African American community remains high and is a primary cause of disparities in life expectancy between African Americans and whites. The objectives of the present scientific statement are to describe cardiovascular health in African Americans and to highlight unique considerations for disease prevention and management.

METHOD

The primary sources of information were identified with PubMed/Medline and online sources from the Centers for Disease Control and Prevention.

RESULTS

The higher prevalence of traditional cardiovascular risk factors (eg, hypertension, diabetes mellitus, obesity, and atherosclerotic cardiovascular risk) underlies the relatively earlier age of onset of cardiovascular diseases among African Americans. Hypertension in particular is highly prevalent among African Americans and contributes directly to the notable disparities in stroke, heart failure, and peripheral artery disease among African Americans. Despite the availability of effective pharmacotherapies and indications for some tailored pharmacotherapies for African Americans (eg, heart failure medications), disease management is less effective among African Americans, yielding higher mortality. Explanations for these persistent disparities in cardiovascular disease are multifactorial and span from the individual level to the social environment.

CONCLUSIONS

The strategies needed to promote equity in the cardiovascular health of African Americans require input from a broad set of stakeholders, including clinicians and researchers from across multiple disciplines.

Genome-wide association study of bronchopulmonary dysplasia: a potential role for variants near the CRP gene

Bronchopulmonary dysplasia (BPD), the main consequence of prematurity, has a significant heritability, but little is known about predisposing genes. The aim of this study was to identify gene loci predisposing infants to BPD. The initial genome-wide association study (GWAS) included 174 Finnish preterm infants of gestational age 24–30 weeks. Thereafter, the most promising single-nucleotide polymorphisms (SNPs) associated with BPD were genotyped in both Finnish (n = 555) and non-Finnish (n = 388) replication cohorts. Finally, plasma CRP levels from the first week of life and the risk of BPD were assessed. SNP rs11265269, flanking the CRP gene, showed the strongest signal in GWAS (odds ratio [OR] 3.2, p = 3.4 × 10⁻⁶). This association was nominally replicated in Finnish and French African populations. A number of other SNPs in the CRP region, including rs3093059, had nominal associations with BPD. During the first week of life the elevated plasma levels of CRP predicted the risk of BPD (OR 3.4, p = 2.9 × 10^–4) and the SNP rs3093059 associated nominally with plasma CRP levels. Finally, SNP rs11265269 was identified as a risk factor of BPD (OR 1.8, p = 5.3 × 10⁻⁵), independently of the robust antenatal risk factors. As such, in BPD, a potential role for variants near CRP gene is proposed.

SEQSpark: A Complete Analysis Tool for Large-Scale Rare Variant Association Studies Using Whole-Genome and Exome Sequence Data

Massively parallel sequencing technologies provide great opportunities for discovering rare susceptibility variants involved in complex disease etiology via large-scale imputation and exome and whole-genome sequence-based association studies. Due to modest effect sizes, large sample sizes of tens to hundreds of thousands of individuals are required for adequately powered studies. Current analytical tools are obsolete when it comes to handling these large datasets. To facilitate the analysis of large-scale sequence-based studies, we developed SEQSpark which implements parallel processing based on Spark to increase the speed and efficiency of performing data quality control, annotation, and association analysis. To demonstrate the versatility and speed of SEQSpark, we analyzed whole-genome sequence data from the UK10K, testing for associations with waist-to-hip ratios. The analysis, which was completed in 1.5 hr, included loading data, annotation, principal component analysis, and single variant and rare variant aggregate association analysis of >9 million variants. For rare variant aggregate analysis, an exome-wide significant association (p < 2.5 × 10^-6) was observed with CCDC62 (SKAT-O [p = 6.89 × 10^-7], combined multivariate collapsing [p = 1.48 × 10^-6], and burden of rare variants [p = 1.48 × 10^-6]). SEQSpark was also used to analyze 50,000 simulated exomes and it required 1.75 hr for the analysis of a quantitative trait using several rare variant aggregate association methods. Additionally, the performance of SEQSpark was compared to Variant Association Tools and PLINK/SEQ. SEQSpark was always faster and in some situations computation was reduced to a hundredth of the time. SEQSpark will empower large sequence-based epidemiological studies to quickly elucidate genetic variation involved in the etiology of complex traits.

Practical Approaches for Whole-Genome Sequence Analysis of Heart- and Blood-Related Traits

Whole-genome sequencing (WGS) allows for a comprehensive view of the sequence of the human genome. We present and apply integrated methodologic steps for interrogating WGS data to characterize the genetic architecture of 10 heart- and blood-related traits in a sample of 1,860 African Americans. In order to evaluate the contribution of regulatory and non-protein coding regions of the genome, we conducted aggregate tests of rare variation across the entire genomic landscape using a sliding window, complemented by an annotation-based assessment of the genome using predefined regulatory elements and within the first intron of all genes. These tests were performed treating all variants equally as well as with individual variants weighted by a measure of predicted functional consequence. Significant findings were assessed in 1,705 individuals of European ancestry. After these steps, we identified and replicated components of the genomic landscape significantly associated with heart- and blood-related traits. For two traits, lipoprotein(a) levels and neutrophil count, aggregate tests of low-frequency and rare variation were significantly associated across multiple motifs. For a third trait, cardiac troponin T, investigation of regulatory domains identified a locus on chromosome 9. These practical approaches for WGS analysis led to the identification of informative genomic regions and also showed that defined non-coding regions, such as first introns of genes and regulatory domains, are associated with important risk factor phenotypes. This study illustrates the tractable nature of WGS data and outlines an approach for characterizing the genetic architecture of complex traits.

Genetics of Coronary Artery Disease

Genetic factors contribute importantly to the risk of coronary artery disease (CAD), and in the past decade, there has been major progress in this area. The tools applied include genome-wide association studies encompassing >200,000 individuals complemented by bioinformatic approaches, including 1000 Genomes imputation, expression quantitative trait locus analyses, and interrogation of Encyclopedia of DNA Elements, Roadmap, and other data sets. close to 60 common SNPs (minor allele frequency>0.05) associated with CAD risk and reaching genome-wide significance (P<5 × 10(-8)) have been identified. Furthermore, a total of 202 independent signals in 109 loci have achieved a false discovery rate (q<0.05) and together explain 28% of the estimated heritability of CAD. These data have been used successfully to create genetic risk scores that can improve risk prediction beyond conventional risk factors and identify those individuals who will benefit most from statin therapy. Such information also has important applications in clinical medicine and drug discovery by using a Mendelian randomization approach to interrogate the causal nature of many factors found to associate with CAD risk in epidemiological studies. In contrast to genome-wide association studies, whole-exome sequencing has provided valuable information directly relevant to genes with known roles in plasma lipoprotein metabolism but has, thus far, failed to identify other rare coding variants linked to CAD. Overall, recent studies have led to a broader understanding of the genetic architecture of CAD and demonstrate that it largely derives from the cumulative effect of multiple common risk alleles individually of small effect size rather than rare variants with large effects on CAD risk. Despite this success, there has been limited progress in understanding the function of the novel loci; the majority of which are in noncoding regions of the genome.

The effect of phenotypic outliers and non-normality on rare-variant association testing

Rare-variant association studies (RVAS) have made important contributions to human complex trait genetics. These studies rely on specialized statistical methods for analyzing rare-variant associations, both individually and in aggregate. We investigated the impact that phenotypic outliers and non-normality have on the performance of rare-variant association testing procedures. Ignoring outliers or non-normality can significantly inflate Type I error rates. We found that rank-based inverse normal transformation (INT) and trait winsorisation were both effective at maintaining Type I error control without sacrificing power in the presence of outliers. INT was the optimal method for non-normally distributed traits. For RVAS of quantitative traits with outliers or non-normality, we recommend using INT to transform phenotypic values before association testing.

Rare Exome Sequence Variants in CLCN6 Reduce Blood Pressure Levels and Hypertension Risk

BACKGROUND

Rare genetic variants influence blood pressure (BP).

METHODS AND RESULTS

Whole-exome sequencing was performed on DNA samples from 17 956 individuals of European ancestry and African ancestry (14 497, first-stage discovery and 3459, second-stage discovery) to examine the effect of rare variants on hypertension and 4 BP traits: systolic BP, diastolic BP, pulse pressure, and mean arterial pressure. Tests of ≈170 000 common variants (minor allele frequency, ≥1%; statistical significance, P≤2.9×10(-7)) and gene-based tests of rare variants (minor allele frequency, <1%; ≈17 000 genes; statistical significance, P≤1.5×10(-6)) were evaluated for each trait and ancestry, followed by multiethnic meta-analyses. In the first-stage discovery, rare coding variants (splicing, stop-gain, stop-loss, nonsynonymous variants, or indels) in CLCN6 were associated with lower diastolic BP (cumulative minor allele frequency, 1.3%; β=-3.20; P=4.1×10(-6)) and were independent of a nearby common variant (rs17367504) previously associated with BP. CLCN6 rare variants were also associated with lower systolic BP (β=-4.11; P=2.8×10(-4)), mean arterial pressure (β=-3.50; P=8.9×10(-6)), and reduced hypertension risk (odds ratio, 0.72; P=0.017). Meta-analysis of the 2-stage discovery samples showed that CLCN6 was associated with lower diastolic BP at exome-wide significance (cumulative minor allele frequency, 1.1%; β=-3.30; P=5.0×10(-7)).

CONCLUSIONS

These findings implicate the effect of rare coding variants in CLCN6 in BP variation and offer new insights into BP regulation.

Genetic evidence for common pathways in human age-related diseases

Aging is the single largest risk factor for chronic disease. Studies in model organisms have identified conserved pathways that modulate aging rate and the onset and progression of multiple age-related diseases, suggesting that common pathways of aging may influence age-related diseases in humans as well. To determine whether there is genetic evidence supporting the notion of common pathways underlying age-related diseases, we analyzed the genes and pathways found to be associated with five major categories of age-related disease using a total of 410 genomewide association studies (GWAS). While only a small number of genes are shared among all five disease categories, those found in at least three of the five major age-related disease categories are highly enriched for apoliprotein metabolism genes. We found that a more substantial number of gene ontology (GO) terms are shared among the 5 age-related disease categories and shared GO terms include canonical aging pathways identified in model organisms, such as nutrient-sensing signaling, translation, proteostasis, stress responses, and genome maintenance. Taking advantage of the vast amount of genetic data from the GWAS, our findings provide the first direct evidence that conserved pathways of aging simultaneously influence multiple age-related diseases in humans as has been demonstrated in model organisms.

Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers

We report ∼17.6 million genetic variants from whole-genome sequencing of 2,120 Sardinians; 22% are absent from previous sequencing-based compilations and are enriched for predicted functional consequences. Furthermore, ∼76,000 variants common in our sample (frequency >5%) are rare elsewhere (<0.5% in the 1000 Genomes Project). We assessed the impact of these variants on circulating lipid levels and five inflammatory biomarkers. We observe 14 signals, including 2 major new loci, for lipid levels and 19 signals, including 2 new loci, for inflammatory markers. The new associations would have been missed in analyses based on 1000 Genomes Project data, underlining the advantages of large-scale sequencing in this founder population.

Collaborative science in the next-generation sequencing era: a viewpoint on how to combine exome sequencing data across sites to identify novel disease susceptibility genes

The purpose of this article is to inform readers about technical challenges that we encountered when assembling exome sequencing data from the 'Simplifying Complex Exomes' (SIMPLEXO) consortium-whose mandate is the discovery of novel genes predisposing to breast and ovarian cancers. Our motivation is to share these obstacles-and our solutions to them-as a means of communicating important technical details that should be discussed early in projects involving massively parallel sequencing.

Discovery of coding genetic variants influencing diabetes-related serum biomarkers and their impact on risk of type 2 diabetes

CONTEXT

Type 2 diabetes (T2D) prevalence is spiraling globally, and knowledge of its pathophysiological signatures is crucial for a better understanding and treatment of the disease.

OBJECTIVE

We aimed to discover underlying coding genetic variants influencing fasting serum levels of nine biomarkers associated with T2D: adiponectin, C-reactive protein, ferritin, heat shock 70-kDa protein 1B, IGF binding protein 1 and IGF binding protein 2, IL-18, IL-2 receptor-α, and leptin.

DESIGN AND PARTICIPANTS

A population-based sample of 6215 adult Danes was genotyped for 16 340 coding single-nucleotide polymorphisms and were tested for association with each biomarker. Identified loci were tested for association with T2D through a large-scale meta-analysis involving up to 17 024 T2D cases and up to 64 186 controls.

RESULTS

We discovered 11 associations between single-nucleotide polymorphisms and five distinct biomarkers at a study-wide P < 3.4 × 10(-7). Nine associations were novel: IL18: BIRC6, RAD17, MARVELD2; ferritin: F5; IGF binding protein 1: SERPING1, KLKB, GCKR, CELSR2, and heat shock 70-kDa protein 1B: CFH. Three of the identified loci (CELSR2, HNF1A, and GCKR) were significantly associated with T2D, of which the association with the CELSR2 locus has not been shown previously.

CONCLUSION

The identified loci influence processes related to insulin signaling, cell communication, immune function, apoptosis, DNA repair, and oxidative stress, all of which could provide a rationale for novel diabetes therapeutic strategies.