Plasma proteome association with coronary heart disease and carotid intima media thickness: results from the KORA F4 study.

BACKGROUND AND AIMS: Atherosclerosis is the main cause of stroke and coronary heart disease (CHD), both leading mortality causes worldwide. Proteomics, as a high-throughput method, could provide helpful insights into the pathological mechanisms underlying atherosclerosis. In this study, we characterized the associations of plasma protein levels with CHD and with carotid intima-media thickness (CIMT), as a surrogate measure of atherosclerosis. METHODS: The discovery phase included 1000 participants from the KORA F4 study, whose plasma protein levels were quantified using the aptamer-based SOMAscan proteomics platform. We evaluated the associations of plasma protein levels with CHD using logistic regression, and with CIMT using linear regression. For both outcomes we applied two models: an age-sex adjusted model, and a model additionally adjusted for body mass index, smoking status, physical activity, diabetes status, hypertension status, low density lipoprotein, high density lipoprotein, and triglyceride levels (fully-adjusted model). The replication phase included a matched case-control sample from the independent KORA F3 study, using ELISA-based measurements of galectin-4. Pathway analysis was performed with nominally associated proteins (p-value < 0.05) from the fully-adjusted model. RESULTS: In the KORA F4 sample, after Bonferroni correction, we found CHD to be associated with five proteins using the age-sex adjusted model: galectin-4 (LGALS4), renin (REN), cathepsin H (CTSH), and coagulation factors X and Xa (F10). The fully-adjusted model yielded only the positive association of galectin-4 (OR = 1.58, 95% CI = 1.30-1.93), which was successfully replicated in the KORA F3 sample (OR = 1.40, 95% CI = 1.09-1.88). For CIMT, we found four proteins to be associated using the age-sex adjusted model namely: cytoplasmic protein NCK1 (NCK1), insulin-like growth factor-binding protein 2 (IGFBP2), growth hormone receptor (GHR), and GDNF family receptor alpha-1 (GFRA1). After assessing the fully-adjusted model, only NCK1 remained significant (ß = 0.017, p-value = 1.39e-06). Upstream regulators of galectin-4 and NCK1 identified from pathway analysis were predicted to be involved in inflammation pathways. CONCLUSIONS: Our proteome-wide association study identified galectin-4 to be associated with CHD and NCK1 to be associated with CIMT. Inflammatory pathways underlying the identified associations highlight the importance of inflammation in the development and progression of CHD.

Epigenome-wide association study of dietary fatty acid intake.

BACKGROUND: Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS). RESULTS: DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively). CONCLUSIONS: Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

Epidemiology of cervical cancer in elderly women: Analysis of incidence, treatment, and survival using German registry data.

BACKGROUND: Cervical cancer (CC) screening is generally recommended until age 65. The incidence of CC could be underestimated, particularly in older women, due to a lack of hysterectomy correction. Furthermore, elderly women (≥65 years) are more often diagnosed with late-stage disease and have worse outcomes than younger patients. This study aims to provide an in-depth overview of CC in Germany. METHODS: Incidence rates of CC (ICD-10 C53) were determined using data from the German Centre of Cancer Registry data (ZfKD) of six federal state registries. Incidence was corrected by using hysterectomy prevalence values from a real-world study. The distribution of treatment modalities (surgery, chemotherapy, radiation therapy) was assessed. Relative survival was calculated using the period approach (2011-2015). Survival was stratified by tumor (T) stage and histological type. RESULTS: In total, 14,528 CC cases were included, 27.6% of which occurred in elderly women. Cumulative (2001-2015) age-standardized incidence rates were 12.5 per 100,000 women without hysterectomy correction and 15.5 per 100,000 women after hysterectomy correction (+24% relative change). A lower proportion of elderly women were treated, especially in advanced tumor stages. Younger women (20-64 years) had a higher 5-year relative survival compared to elderly women: 76.7% versus 46.9%, respectively. Survival was worse with increasing stage and for glandular histological subgroups, particularly among elderly women. CONCLUSIONS: CC incidence in elderly women is underestimated and survival is lower compared to younger women in Germany. Due to the high disease burden in elderly women, screening and treatment strategies need to be improved.

Mendelian randomisation study of smoking exposure in relation to breast cancer risk.

BACKGROUND: Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk. METHODS: We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy. RESULTS: Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect. CONCLUSION: Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Whole-genome sequencing analysis of the cardiometabolic proteome.

The human proteome is a crucial intermediate between complex diseases and their genetic and environmental components, and an important source of drug development targets and biomarkers. Here, we comprehensively assess the genetic architecture of 257 circulating protein biomarkers of cardiometabolic relevance through high-depth (22.5×) whole-genome sequencing (WGS) in 1328 individuals. We discover 131 independent sequence variant associations (P < 7.45 × 10-11) across the allele frequency spectrum, all of which replicate in an independent cohort (n = 1605, 18.4x WGS). We identify for the first time replicating evidence for rare-variant cis-acting protein quantitative trait loci for five genes, involving both coding and noncoding variation. We construct and validate polygenic scores that explain up to 45% of protein level variation. We find causal links between protein levels and disease risk, identifying high-value biomarkers and drug development targets.

Genetic Variants in the Regulatory T cell-Related Pathway and Colorectal Cancer Prognosis.

BACKGROUND: High numbers of lymphocytes in tumor tissue, including T regulatory cells (Treg), have been associated with better colorectal cancer survival. Tregs, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and therefore variants in genes related to Treg differentiation and function could be associated with colorectal cancer prognosis. METHODS: In a prospective German cohort of 3,593 colorectal cancer patients, we assessed the association of 771 single-nucleotide polymorphisms (SNP) in 58 Treg-related genes with overall and colorectal cancer-specific survival using Cox regression models. Effect modification by microsatellite instability (MSI) status was also investigated because tumors with MSI show greater lymphocytic infiltration and have been associated with better prognosis. Replication of significant results was attempted in 2,047 colorectal cancer patients of the International Survival Analysis in Colorectal Cancer Consortium (ISACC). RESULTS: A significant association of the TGFBR3 SNP rs7524066 with more favorable colorectal cancer-specific survival [hazard ratio (HR) per minor allele: 0.83; 95% confidence interval (CI), 0.74-0.94; P value: 0.0033] was replicated in ISACC (HR: 0.82; 95% CI, 0.68-0.98; P value: 0.03). Suggestive evidence for association was found with two IL7 SNPs, rs16906568 and rs7845577. Thirteen SNPs with differential associations with overall survival according to MSI in the discovery analysis were not confirmed. CONCLUSIONS: Common genetic variation in the Treg pathway implicating genes such as TGFBR3 and IL7 was shown to be associated with prognosis of colorectal cancer patients. IMPACT: The implicated genes warrant further investigation.

Genetic Predictors of Circulating 25-Hydroxyvitamin D and Prognosis after Colorectal Cancer.

BACKGROUND: Low serum 25-hydroxyvitamin D [25(OH)D] concentrations in patients with colorectal cancer have been consistently associated with higher mortality in observational studies. It is unclear whether low 25(OH)D levels directly influence colorectal cancer mortality. To minimize bias, we use genetic variants associated with vitamin D levels to evaluate the association with overall and colorectal cancer-specific survival. METHODS: Six genetic variants have been robustly identified to be associated with 25(OH)D levels in genome-wide association studies. On the basis of data from the International Survival Analysis in Colorectal Cancer Consortium, the individual genetic variants and a weighted genetic risk score were tested for association with overall and colorectal cancer-specific survival using Cox proportional hazards models in 7,657 patients with stage I to IV colorectal cancer, of whom 2,438 died from any cause and 1,648 died from colorectal cancer. RESULTS: The 25(OH)D decreasing allele of SNP rs2282679 (GC gene, encodes group-specific component/vitamin D-binding protein) was associated with poorer colorectal cancer-specific survival, although not significant after multiple-testing correction. None of the other five SNPs showed an association. The genetic risk score showed nonsignificant associations with increased overall [HR = 1.54; confidence interval (CI), 0.86-2.78] and colorectal cancer-specific mortality (HR = 1.76; 95% CI, 0.86-3.58). A significant increased risk of overall mortality was observed in women (HR = 3.26; 95% CI, 1.45-7.33; P heterogeneity = 0.01) and normal-weight individuals (HR = 4.14; 95% CI, 1.50-11.43, P heterogeneity = 0.02). CONCLUSIONS: Our results provided little evidence for an association of genetic predisposition of lower vitamin D levels with increased overall or colorectal cancer-specific survival, although power might have been an issue. IMPACT: Further studies are warranted to investigate the association in specific subgroups.

Strengthening Causal Inference for Complex Disease Using Molecular Quantitative Trait Loci.

Large genome-wide association studies (GWAS) have identified loci that are associated with complex traits and diseases, but index variants are often not causal and reside in non-coding regions of the genome. To gain a better understanding of the relevant biological mechanisms, intermediate traits such as gene expression and protein levels are increasingly being investigated because these are likely mediators between genetic variants and disease outcome. Genetic variants associated with intermediate traits, termed molecular quantitative trait loci (molQTLs), can then be used as instrumental variables in a Mendelian randomization (MR) approach to identify the causal features and mechanisms of complex traits. Challenges such as pleiotropy and the non-specificity of molQTLs remain, and further approaches and methods need to be developed.

DNA methylation profiling to explore colorectal tumor differences according to menopausal hormone therapy use in women.

Aim: Use of menopausal hormone therapy (MHT) has been associated with a reduced risk for colorectal cancer, but mechanisms underlying this relationship are not well understood. In the colon, MHT appears to act through estrogen receptor ß (ERß) which may influence DNA methylation by binding to DNA. Using genome-wide methylation profiling data, we aimed to identify genes that may be differentially methylated according to MHT use. Materials & methods: DNA methylation was measured using Illumina HumanMethylation450k arrays in two independent tumor sample sets of colorectal cancer patients. Differential methylation was determined using R/limma. Results: In the discovery analysis, two CpG sites showed differential DNA methylation according to MHT use, both were not replicated. In stratified analyses, 342 CpG sites were associated with current MHT use only in ERß-positive tumors. Conclusion: The suggestive findings of differential methylation according to current MHT use in ERß-positive tumors warrant further investigation.

Genome-wide DNA methylation differences according to oestrogen receptor beta status in colorectal cancer.

Involvement of sex hormones in colorectal cancer (CRC) development has been linked to oestrogen receptor ß (ERß). Expression of ERß is found reduced in tumour tissue and inversely related to mortality. However, mechanisms are not well understood. Our study aimed to detect differentially methylated genes associated with ERß expression, which could point to mechanisms by which ERß could influence risk and prognosis of CRC. Epigenome-wide DNA methylation profiling was performed using Illumina HumanMethylation450k BeadChip arrays in two independent tumour sample sets of CRC patients recruited in 2003-2010 by the German DACHS study (discovery cohort n = 917, replication cohort n = 907). ERß expression was measured using immunohistochemistry and scored as negative, moderate and high. Differentially methylated CpG sites and genomic regions were determined using limma in the R-package RnBeads. For the comparison of tumours with moderate/high ERß versus negative expression, differentially methylated CpG sites were identified but not confirmed by replication. Comparing tumours of high with tumours of negative ERß expression revealed 2,904 differentially methylated CpG sites of which 403 were replicated (FDR adjusted p-value<0.05). Replicated CpGs were annotated to genes such as CD36, HK1 or LRP5. A survival analysis indicates that 30 of the replicated CpGs are also associated with overall survival (FDR-adjusted p-value<0.05). The regional analysis identified 60 differentially methylated promotor regions. The epigenome-wide analysis identified both novel genes as well as genes already implicated in CRC. Follow-up mechanistic studies to better understand the regulatory role of ERß could inform potential targets for improving treatment or prevention of CRC.

Mendelian randomisation study of age at menarche and age at menopause and the risk of colorectal cancer.

BACKGROUND: Substantial evidence supports an association between use of menopausal hormone therapy and decreased colorectal cancer (CRC) risk, indicating a role of exogenous sex hormones in CRC development. However, findings on endogenous oestrogen exposure and CRC are inconsistent. METHODS: We used a Mendelian randomisation approach to test for a causal effect of age at menarche and age at menopause as surrogates for endogenous oestrogen exposure on CRC risk. Weighted genetic risk scores based on 358 single-nucleotide polymorphisms associated with age at menarche and 51 single-nucleotide polymorphisms associated with age at menopause were used to estimate the association with CRC risk using logistic regression in 12,944 women diagnosed with CRC and 10,741 women without CRC from three consortia. Sensitivity analyses were conducted to address pleiotropy and possible confounding by body mass index. RESULTS: Genetic risk scores for age at menarche (odds ratio per year 0.98, 95% confidence interval: 0.95-1.02) and age at menopause (odds ratio 0.98, 95% confidence interval: 0.94-1.01) were not significantly associated with CRC risk. The sensitivity analyses yielded similar results. CONCLUSIONS: Our study does not support a causal relationship between genetic risk scores for age at menarche and age at menopause and CRC risk.