Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.

Central adiposity, obesity during early adulthood, and pancreatic cancer mortality in a pooled analysis of cohort studies

BACKGROUND

Body mass index (BMI), a measure of obesity typically assessed in middle age or later, is known to be positively associated with pancreatic cancer. However, little evidence exists regarding the influence of central adiposity, a high BMI during early adulthood, and weight gain after early adulthood on pancreatic cancer risk.

DESIGN

We conducted a pooled analysis of individual-level data from 20 prospective cohort studies in the National Cancer Institute BMI and Mortality Cohort Consortium to examine the association of pancreatic cancer mortality with measures of central adiposity (e.g. waist circumference; n = 647 478; 1947 pancreatic cancer deaths), BMI during early adulthood (ages 18-21 years) and BMI change between early adulthood and cohort enrollment, mostly in middle age or later (n = 1 096 492; 3223 pancreatic cancer deaths). Multivariable hazard ratios (HRs) and 95% confidence intervals (CIs) were calculated using Cox proportional hazards regression models.

RESULTS

Higher waist-to-hip ratio (HR = 1.09, 95% CI 1.02-1.17 per 0.1 increment) and waist circumference (HR = 1.07, 95% CI 1.00-1.14 per 10 cm) were associated with increased risk of pancreatic cancer mortality, even when adjusted for BMI at baseline. BMI during early adulthood was associated with increased pancreatic cancer mortality (HR = 1.18, 95% CI 1.11-1.25 per 5 kg/m(2)), with increased risk observed in both overweight and obese individuals (compared with BMI of 21.0 to <23 kg/m(2), HR = 1.36, 95% CI 1.20-1.55 for BMI 25.0 < 27.5 kg/m(2), HR = 1.48, 95% CI 1.20-1.84 for BMI 27.5 to <30 kg/m(2), HR = 1.43, 95% CI 1.11-1.85 for BMI ≥30 kg/m(2)). BMI gain after early adulthood, adjusted for early adult BMI, was less strongly associated with pancreatic cancer mortality (HR = 1.05, 95% CI 1.01-1.10 per 5 kg/m(2)).

CONCLUSIONS

Our results support an association between pancreatic cancer mortality and central obesity, independent of BMI, and also suggest that being overweight or obese during early adulthood may be important in influencing pancreatic cancer mortality risk later in life.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Cross Cancer Genomic Investigation of Inflammation Pathway for Five Common Cancers: Lung, Ovary, Prostate, Breast, and Colorectal Cancer

BACKGROUND

Inflammation has been hypothesized to increase the risk of cancer development as an initiator or promoter, yet no large-scale study of inherited variation across cancer sites has been conducted.

METHODS

We conducted a cross-cancer genomic analysis for the inflammation pathway based on 48 genome-wide association studies within the National Cancer Institute GAME-ON Network across five common cancer sites, with a total of 64 591 cancer patients and 74 467 control patients. Subset-based meta-analysis was used to account for possible disease heterogeneity, and hierarchical modeling was employed to estimate the effect of the subcomponents within the inflammation pathway. The network was visualized by enrichment map. All statistical tests were two-sided.

RESULTS

We identified three pleiotropic loci within the inflammation pathway, including one novel locus in Ch12q24 encoding SH2B3 (rs3184504), which reached GWAS significance with a P value of 1.78 x 10(-8), and it showed an association with lung cancer (P = 2.01 x 10(-6)), colorectal cancer (GECCO P = 6.72x10(-6); CORECT P = 3.32x10(-5)), and breast cancer (P = .009). We also identified five key subpathway components with genetic variants that are relevant for the risk of these five cancer sites: inflammatory response for colorectal cancer (P = .006), inflammation related cell cycle gene for lung cancer (P = 1.35x10(-6)), and activation of immune response for ovarian cancer (P = .009). In addition, sequence variations in immune system development played a role in breast cancer etiology (P = .001) and innate immune response was involved in the risk of both colorectal (P = .022) and ovarian cancer (P = .003).

CONCLUSIONS

Genetic variations in inflammation and its related subpathway components are keys to the development of lung, colorectal, ovary, and breast cancer, including SH2B3, which is associated with lung, colorectal, and breast cancer.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Genetic architecture of colorectal cancer

Colorectal cancer (CRC) is a complex disease that develops as a consequence of both genetic and environmental risk factors. A small proportion (3-5%) of cases arise from hereditary syndromes predisposing to early onset CRC as a result of mutations in over a dozen well defined genes. In contrast, CRC is predominantly a late onset 'sporadic' disease, developing in individuals with no obvious hereditary syndrome. In recent years, genome wide association studies have discovered that over 40 genetic regions are associated with weak effects on sporadic CRC, and it has been estimated that increasingly large genome wide scans will identify many additional novel genetic regions. Subsequent experimental validations have identified the causally related variant(s) in a limited number of these genetic regions. Further biological insight could be obtained through ethnically diverse study populations, larger genetic sequencing studies and development of higher throughput functional experiments. Along with inherited variation, integration of the tumour genome may shed light on the carcinogenic processes in CRC. In addition to summarising the genetic architecture of CRC, this review discusses genetic factors that modify environmental predictors of CRC, as well as examples of how genetic insight has improved clinical surveillance, prevention and treatment strategies. In summary, substantial progress has been made in uncovering the genetic architecture of CRC, and continued research efforts are expected to identify additional genetic risk factors that further our biological understanding of this disease. Subsequently these new insights will lead to improved treatment and prevention of colorectal cancer.

Rare and Coding Region Genetic Variants Associated With Risk of Ischemic Stroke: The NHLBI Exome Sequence Project

IMPORTANCE

Stroke is the second leading cause of death and the third leading cause of years of life lost. Genetic factors contribute to stroke prevalence, and candidate gene and genome-wide association studies (GWAS) have identified variants associated with ischemic stroke risk. These variants often have small effects without obvious biological significance. Exome sequencing may discover predicted protein-altering variants with a potentially large effect on ischemic stroke risk.

OBJECTIVE

To investigate the contribution of rare and common genetic variants to ischemic stroke risk by targeting the protein-coding regions of the human genome.

DESIGN, SETTING, AND PARTICIPANTS

The National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP) analyzed approximately 6000 participants from numerous cohorts of European and African ancestry. For discovery, 365 cases of ischemic stroke (small-vessel and large-vessel subtypes) and 809 European ancestry controls were sequenced; for replication, 47 affected sibpairs concordant for stroke subtype and an African American case-control series were sequenced, with 1672 cases and 4509 European ancestry controls genotyped. The ESP's exome sequencing and genotyping started on January 1, 2010, and continued through June 30, 2012. Analyses were conducted on the full data set between July 12, 2012, and July 13, 2013.

MAIN OUTCOMES AND MEASURES

Discovery of new variants or genes contributing to ischemic stroke risk and subtype (primary analysis) and determination of support for protein-coding variants contributing to risk in previously published candidate genes (secondary analysis).

RESULTS

We identified 2 novel genes associated with an increased risk of ischemic stroke: a protein-coding variant in PDE4DIP (rs1778155; odds ratio, 2.15; P = 2.63 × 10(-8)) with an intracellular signal transduction mechanism and in ACOT4 (rs35724886; odds ratio, 2.04; P = 1.24 × 10(-7)) with a fatty acid metabolism; confirmation of PDE4DIP was observed in affected sibpair families with large-vessel stroke subtype and in African Americans. Replication of protein-coding variants in candidate genes was observed for 2 previously reported GWAS associations: ZFHX3 (cardioembolic stroke) and ABCA1 (large-vessel stroke).

CONCLUSIONS AND RELEVANCE

Exome sequencing discovered 2 novel genes and mechanisms, PDE4DIP and ACOT4, associated with increased risk for ischemic stroke. In addition, ZFHX3 and ABCA1 were discovered to have protein-coding variants associated with ischemic stroke. These results suggest that genetic variation in novel pathways contributes to ischemic stroke risk and serves as a target for prediction, prevention, and therapy.

Abstract 862: Association of non-steroidal anti-inflammatory drugs with colorectal cancer by subgroups in the Vitamins and Lifestyle (VITAL) Study

Background: There is substantial evidence that use of non-steroidal anti-inflammatory drugs (NSAIDs) reduces the risk of colorectal cancer (CRC), but no subgroup has been identified for which the chemoprevention effect outweighs the risk of side effects. Methods: We tested the interaction between NSAID use and multiple risk factors on CRC risk in the VITAL cohort. A total of 73,458 individuals aged 50-76 completed a questionnaire between 2000 and 2002, and 674 incidental colorectal cancer cases were identified through 2010. Results: In stratified analysis, high use of any type of NSAIDs (4+days/week for 4+ years) was statistically significantly associated with a lower risk of CRC across all subgroups stratified by sex, BMI, physical activity, smoking, alcohol intake, screening and dietary factors, with a suggestion of stronger association among men, obese individuals, and heavier drinkers; however, none of these tests for interaction reached statistical significance. The associations were almost identical for subjects with higher overall CRC risk scores (HR: 0.62; 95% CI: 0.49-0.79) and those with lower risk scores (HR: 0.61; 95% CI: 0.42-0.88). Differential effects by cancer subsites and stages were tested. NSAID use was associated with a greater risk reduction of proximal colon cancer vs distal (p for difference = 0.06) and distant stage vs local (p for difference = 0.04). Conclusion: The association between high use of NSAIDs and CRC risk does not differ significantly among subgroups. Impact: Our results suggest that NSAIDs have a generally beneficial role in colorectal cancer prevention, largely unmodified by other exposures.

Citation Format: Xiaoliang Wang, Ulrike Peters, John D. Potter, Emily White. Association of non-steroidal anti-inflammatory drugs with colorectal cancer by subgroups in the Vitamins and Lifestyle (VITAL) Study. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 862. doi:10.1158/1538-7445.AM2015-862

Analyses of 7,635 Patients with Colorectal Cancer Using Independent Training and Validation Cohorts Show That rs9929218 in CDH1 Is a Prognostic Marker of Survival

PURPOSE

Genome-wide association studies have identified numerous loci associated with colorectal cancer risk. Several of these have also been associated with patient survival, although none have been validated. Here, we used large independent training and validation cohorts to identify robust prognostic biomarkers for colorectal cancer.

EXPERIMENTAL DESIGN

In our training phase, we analyzed 20 colorectal cancer-risk SNPs from 14 genome-wide associated loci, for their effects on survival in 2,083 patients with advanced colorectal cancer. A Cox survival model was used, stratified for treatment, adjusted for known prognostic factors, and corrected for multiple testing. Three SNPs were subsequently analyzed in an independent validation cohort of 5,552 colorectal cancer patients. A validated SNP was analyzed by disease stage and response to treatment.

RESULTS

Three variants associated with survival in the training phase; however, only rs9929218 at 16q22 (intron 2 of CDH1, encoding E-cadherin) was significant in the validation phase. Patients homozygous for the minor allele (AA genotype) had worse survival (training phase HR, 1.43; 95% confidence intervals; CI, 1.20-1.71, P = 5.8 × 10(-5); validation phase HR, 1.18; 95% CI, 1.01-1.37, P = 3.2 × 10(-2); combined HR, 1.28; 95% CI, 1.14-1.43, P = 2.2 × 10(-5)). This effect was independent of known prognostic factors, and was significant amongst patients with stage IV disease (P = 2.7 × 10(-5)). rs9929218 was also associated with poor response to chemotherapy (P = 3.9 × 10(-4)).

CONCLUSIONS

We demonstrate the potential of common inherited genetic variants to inform patient outcome and show that rs9929218 identifies approximately 8% of colorectal cancer patients with poor prognosis. rs9929218 may affect CDH1 expression and E-cadherin plays a role in epithelial-to-mesenchymal transition providing a mechanism underlying its prognostic potential. Clin Cancer Res; 21(15); 3453-61. ©2015 AACR.

Abstract 4609: Risk loci in telomere structure and maintenance genes across five cancer types: GAME-ON Consortium

Telomeres are complex structures that cap chromosome ends, protecting them from degradation, double strand breaks and end-to-end fusions. Telomeres are maintained by the enzyme telomerase, which is made up of a reverse transcriptase encoded by TERT, and an RNA template encoded by TERC. The telomere structure itself is composed of proteins encoded by ACD, ACYP2, BICD1, DKC1, DCLRE1B, MPHOSPH6, NAF1, NOLA1, NOLA2, NOLA3, OBFC1, PIK3C3, POT1, RTEL1, TEP1, TERF1, TERF2, TERF21P, TINF2, TNKS/PINX1, TPP1 and ZNF208. Several single nucleotide polymorphisms (SNPs) in the TERT and adjoining CLPTM1L gene region are associated with multiple cancer types, and some are associated with both increased and decreased risks across different cancer types. We systematically characterized the patterns of association between variants in these 25 telomere structure and maintenance genes and risk across five cancer types in the Genetic Association and Mechanisms in Oncology (GAME-ON) consortium.

We performed a subset-based meta-analysis (ASSET) of 209,367 directly measured and imputed SNPs, one megabase up- and downstream of these genes, across genome-wide association studies of colorectal (5,100 cases, 4,831 controls), lung (12,160 cases, 16,838 controls), breast (15,748 cases, 18,084 controls), ovarian (4,369 cases, 9,123 controls) and prostate (14,160 cases, 12,724 controls) cancers. Correlations (r2) between SNPs were examined in Haploview using the 1000 Genomes Project CEU population.

A total of 87 TERT, 123 TERC and 26 DCLRE1B SNPs were associated with cancer risk at gene-level Bonferroni-corrected p-values of 4.2-7.8×10−6. Patterns of association were similar for prostate and colorectal cancers in DCLRE1B. Much stronger associations were observed in TERT and TERC, with 63 and 24 SNPs reaching genome-wide significance (p&lt;5.0×10−8), respectively. Of these, 9 TERT and 9 TERC SNPs were correlated at r2&lt;0.75. The most strongly associated SNPs in TERT (rs37004, p = 2.6×10−11) and TERC (rs76925190, p = 1.5×10−15) are in regions containing documented risk loci for lung and prostate cancers, respectively. TERT rs37004 (p = 1.2×10−13) and 6 of the other 9 SNPs in TERT were associated only with lung cancer risk. For the other 2 SNPs, there was a suggestion that patterns of risk were opposite for lung and prostate cancers. While TERC rs76925190 was strongly associated with prostate cancer risk (p = 5.4 x10−17), it was also suggestively associated with colorectal cancer risk (p&lt;0.005). This pattern was similar for another 2 of the 9 SNPs, while for the rest, associations were limited to prostate cancer.

The complex patterns of association in telomere structure and maintenance genes observed across cancer types may provide insight about the mechanisms through which telomere dysfunction in different tissues influences cancer risk.

Citation Format: Sara Karami, Younghun Han, Fredrick R. Schumacher, Zsofia Kote-Jarai, Sara Lindstrom, John S. Witte, Iona Cheng, Shenying Fang, Jiali Han, Peter Kraft, Fengju Song, Rayjean J. Hung, James McKay, Stephen J. Chanock, Mala Pande, Angela Risch, Hongbing Shen, Christopher A. Haiman, Lisa Boardman, Cornelia M. Ulrich, Graham Casey, Ulrike Peters, Nilanjan Chatterjee, Brandon Pierce, Wei Zheng, Christopher I. Amos, Jennifer A. Doherty. Risk loci in telomere structure and maintenance genes across five cancer types: GAME-ON Consortium. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4609. doi:10.1158/1538-7445.AM2015-4609

Abstract 854: Ovarian cancer risk factors by histologic subtypes: evidence for etiologic heterogeneity

Background: A subset of high grade serous carcinomas may arise from the fallopian tube, while some endometrioid and clear cell carcinomas may derive from endometrial tissue. Previous studies have suggested differences in ovarian cancer risk factors by histologic subtypes, but had limited sample sizes. In the Ovarian Cancer Cohort Consortium (OC3), we evaluated associations of reproductive, hormonal, demographic and lifestyle factors, and family history of cancer with ovarian cancer subtypes. Identification of potential differences in associations among subtypes is important for clarifying ovarian cancer etiology and for developing novel prevention and risk prediction approaches.

Methods: Among over 1.2 million women from 21 cohort studies, 4,347 ovarian cancers with histology information were identified during follow-up (3223 serous, 555 endometrioid, 316 mucinous, 253 clear cell). We used competing risks Cox proportional hazards regression to compute risk factor associations by histologic subtype. Models were stratified on study and year of birth and adjusted for age, parity and oral contraceptive use; subtype heterogeneity was evaluated by a likelihood ratio test. Unsupervised hierarchical clustering was used to evaluate patterns of risk factors by histology.

Results: Most risk factors showed significant heterogeneity across histologic subtypes. Higher parity was most strongly associated with lower risks of endometrioid (RR per child: 0.79; 95%CI: 0.74-0.85) and clear cell (RR: 0.69; 95%CI: 0.61-0.78) carcinomas (p-het&lt;0.0001). Age at menopause was positively and tubal ligation was inversely associated only with endometrioid and clear cell carcinomas (p-het = 0.02 and 0.003, respectively). Long-term menopausal hormone use (&gt;5 years) was associated with endometrioid carcinomas (RR: 2.23; 95% CI: 1.46-3.42) and serous carcinomas (RR: 1.66; 95% CI: 1.44-1.92), and inversely associated with clear cell carcinomas (RR: 0.43; 95% CI: 0.20-0.91; p-het = 0.001). Family history of breast cancer was associated with increased risk of serous carcinomas (RR:1.13; 95% CI:1.02-1.27) and endometrioid carcinomas (RR: 1.44; 95% CI: 1.12-1.87; p-het = 0.02). Smoking (per 20 pack years) showed a positive association with mucinous carcinomas (RR: 1.38; 95% CI: 1.09-1.75) and an inverse association with clear cell carcinomas (RR:0.62; 95% CI: 0.46-0.85; p-het = 0.001). In unsupervised hierarchical clustering, serous and mucinous carcinomas clustered in one group and endometrioid and clear cell carcinomas in the other.

Conclusion: Our results demonstrate heterogeneous associations of risk factors with ovarian cancer subtypes, supporting the hypothesis that the subtypes develop through different pathways. Most established risk factors were more strongly associated with non-serous carcinomas, suggesting that risk prediction may be more challenging for serous cancers, the most fatal subtype.

Citation Format: Nicolas A. Wentzensen, Elizabeth Poole, Alan A. Arslan, Alpa V. Patel, V Wendy Setiawan, Kala Visvanathan, Elisabete Weiderpass, Emily White, Hans-Olov Adami, Louise A. Brinton, Leslie Bernstein, Julie Buring, Lesley M. Butler, Saioa Chamosa, Tess V. Clendenen, Laure Dossus, Renee Fortner, Susan M. Gapstur, Mia M. Gaudet, Inger Torhild Gram, Patricia Hartge, Judith Hoffman-Bolton, Annika Idahl, Michael Jones, Rudolf Kaaks, Vivki Kirsh, Woon-Puay Koh, James V. Lacey, I-Min Lee, Eva Lundin, Melissa Merritt, Ulrike Peters, Jenny Poynter, Sabina Rinaldi, Kim Robien, Thomas Rohan, Dale P. Sandler, Leo J. Schouten, Louise Sjöholm, Sabina Sieri, Anthony Swerdlow, Anne Tjønneland, Britton Trabert, Lynne Wilkens, Alicja Wolk, Hannah P. Yang, Anne Zeleniuch-Jacquotte, Shelley S. Tworoger. Ovarian cancer risk factors by histologic subtypes: evidence for etiologic heterogeneity. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 854. doi:10.1158/1538-7445.AM2015-854

Serum 25-hydroxyvitamin D concentrations and risk of prostate cancer: results from the Prostate Cancer Prevention Trial

BACKGROUND

Epidemiologic studies have reported inconsistent associations of vitamin D and prostate cancer risk; however, few have adequately controlled for detection bias related to prostate-specific antigen (PSA) screening, and the results of many studies may be affected by occult prostate cancers among controls.

METHODS

Data for this nested case-control analysis (n = 1,695 cases/1,682 controls) are from the Prostate Cancer Prevention Trial. Baseline serum was analyzed for 25-hydroxyvitamin D [25(OH)D]. The presence or absence of cancer was subsequently determined by prostate biopsy. Polytomous logistic regression models were used to estimate associations of 25(OH)D with risk of total, Gleason 2-6, Gleason 7, and Gleason 8-10 prostate cancer. Results are presented for placebo and finasteride arms separately and combined.

RESULTS

There were no associations of serum 25(OH)D with total prostate cancer risk. For Gleason 2-6 cancers, results were inconsistent across treatment arms with a suggestion of increased risk in the placebo arm only; however, there was no dose-response relationship. For Gleason 8-10 prostate cancers, 25(OH)D concentrations were associated with a linear decrease in risk among combined treatment arms [quartile 4 vs. 1: OR, 0.55; 95% confidence interval (CI), 0.32-0.94; P(trend) = 0.04]. These findings were somewhat stronger among men ≥65 versus 55-64 years at baseline (quartile 4 vs. 1: OR, 0.40; 95% CI, 0.18-0.88 vs. OR, 0.73; 95% CI, 0.35-1.52, respectively; P(interaction) = 0.52).

CONCLUSIONS

Higher serum 25(OH)D may modestly increase risk of Gleason 2-6 disease and more substantially reduce risk of Gleason 8-10 prostate cancer.

IMPACT

Vitamin D may have different effects for different stages of prostate cancers.

Identification of a common variant with potential pleiotropic effect on risk of inflammatory bowel disease and colorectal cancer

Although genome-wide association studies (GWAS) have separately identified many genetic susceptibility loci for ulcerative colitis (UC), Crohn's disease (CD) and colorectal cancer (CRC), there has been no large-scale examination for pleiotropy, or shared genetic susceptibility, for these conditions. We used logistic regression modeling to examine the associations of 181 UC and CD susceptibility variants previously identified by GWAS with risk of CRC using data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. We also examined associations of significant variants with clinical and molecular characteristics in a subset of the studies. Among 11794 CRC cases and 14190 controls, rs11676348, the susceptibility single nucleotide polymorphism (SNP) for UC, was significantly associated with reduced risk of CRC (P = 7E-05). The multivariate-adjusted odds ratio of CRC with each copy of the T allele was 0.93 (95% CI 0.89-0.96). The association of the SNP with risk of CRC differed according to mucinous histological features (P heterogeneity = 0.008). In addition, the (T) allele was associated with lower risk of tumors with Crohn's-like reaction but not tumors without such immune infiltrate (P heterogeneity = 0.02) and microsatellite instability-high (MSI-high) but not microsatellite stable or MSI-low tumors (P heterogeneity = 0.03). The minor allele (T) in SNP rs11676348, located downstream from CXCR2 that has been implicated in CRC progression, is associated with a lower risk of CRC, particularly tumors with a mucinous component, Crohn's-like reaction and MSI-high. Our findings offer the promise of risk stratification of inflammatory bowel disease patients for complications such as CRC.

Powerful Set-Based Gene-Environment Interaction Testing Framework for Complex Diseases

Identification of gene-environment interaction (G × E) is important in understanding the etiology of complex diseases. Based on our previously developed Set Based gene EnviRonment InterAction test (SBERIA), in this paper we propose a powerful framework for enhanced set-based G × E testing (eSBERIA). The major challenge of signal aggregation within a set is how to tell signals from noise. eSBERIA tackles this challenge by adaptively aggregating the interaction signals within a set weighted by the strength of the marginal and correlation screening signals. eSBERIA then combines the screening-informed aggregate test with a variance component test to account for the residual signals. Additionally, we develop a case-only extension for eSBERIA (coSBERIA) and an existing set-based method, which boosts the power not only by exploiting the G-E independence assumption but also by avoiding the need to specify main effects for a large number of variants in the set. Through extensive simulation, we show that coSBERIA and eSBERIA are considerably more powerful than existing methods within the case-only and the case-control method categories across a wide range of scenarios. We conduct a genome-wide G × E search by applying our methods to Illumina HumanExome Beadchip data of 10,446 colorectal cancer cases and 10,191 controls and identify two novel interactions between nonsteroidal anti-inflammatory drugs (NSAIDs) and MINK1 and PTCHD3.

MiRNA-Related SNPs and Risk of Esophageal Adenocarcinoma and Barrett's Esophagus: Post Genome-Wide Association Analysis in the BEACON Consortium

Incidence of esophageal adenocarcinoma (EA) has increased substantially in recent decades. Multiple risk factors have been identified for EA and its precursor, Barrett’s esophagus (BE), such as reflux, European ancestry, male sex, obesity, and tobacco smoking, and several germline genetic variants were recently associated with disease risk. Using data from the Barrett’s and Esophageal Adenocarcinoma Consortium (BEACON) genome-wide association study (GWAS) of 2,515 EA cases, 3,295 BE cases, and 3,207 controls, we examined single nucleotide polymorphisms (SNPs) that potentially affect the biogenesis or biological activity of microRNAs (miRNAs), small non-coding RNAs implicated in post-transcriptional gene regulation, and deregulated in many cancers, including EA. Polymorphisms in three classes of genes were examined for association with risk of EA or BE: miRNA biogenesis genes (157 SNPs, 21 genes); miRNA gene loci (234 SNPs, 210 genes); and miRNA-targeted mRNAs (177 SNPs, 158 genes). Nominal associations (P<0.05) of 29 SNPs with EA risk, and 25 SNPs with BE risk, were observed. None remained significant after correction for multiple comparisons (FDR q>0.50), and we did not find evidence for interactions between variants analyzed and two risk factors for EA/BE (smoking and obesity). This analysis provides the most extensive assessment to date of miRNA-related SNPs in relation to risk of EA and BE. While common genetic variants within components of the miRNA biogenesis core pathway appear unlikely to modulate susceptibility to EA or BE, further studies may be warranted to examine potential associations between unassessed variants in miRNA genes and targets with disease risk.

A model to determine colorectal cancer risk using common genetic susceptibility loci

BACKGROUND & AIMS

Risk for colorectal cancer (CRC) can be greatly reduced through screening. To aid in the development of screening strategies, we refined models designed to determine risk of CRC by incorporating information from common genetic susceptibility loci.

METHODS

By using data collected from more than 12,000 participants in 6 studies performed from 1990 through 2011 in the United States and Germany, we developed risk determination models based on sex, age, family history, genetic risk score (number of risk alleles carried at 27 validated common CRC susceptibility loci), and history of endoscopic examinations. The model was validated using data collected from approximately 1800 participants in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial, conducted from 1993 through 2001 in the United States.

RESULTS

We identified a CRC genetic risk score that independently predicted which patients in the training set would develop CRC. Compared with determination of risk based only on family history, adding the genetic risk score increased the discriminatory accuracy from 0.51 to 0.59 (P = .0028) for men and from 0.52 to 0.56 (P = .14) for women. We calculated age- and sex-specific 10-year CRC absolute risk estimates based on the number of risk alleles, family history, and history of endoscopic examinations. A model that included a genetic risk score better determined the recommended starting age for screening in subjects with and without family histories of CRC. The starting age for high-risk men (family history of CRC and genetic risk score, 90%) was 42 years, and for low-risk men (no family history of CRC and genetic risk score, 10%) was 52 years. For men with no family history and a high genetic risk score (90%), the starting age would be 47 years; this is an intermediate value that is 5 years earlier than it would be for men with a genetic risk score of 10%. Similar trends were observed in women.

CONCLUSIONS

By incorporating information on CRC risk alleles, we created a model to determine the risk for CRC more accurately. This model might be used to develop screening and prevention strategies.

Mendelian randomization study of height and risk of colorectal cancer

BACKGROUND

For men and women, taller height is associated with increased risk of all cancers combined. For colorectal cancer (CRC), it is unclear whether the differential association of height by sex is real or is due to confounding or bias inherent in observational studies. We performed a Mendelian randomization study to examine the association between height and CRC risk.

METHODS

To minimize confounding and bias, we derived a weighted genetic risk score predicting height (using 696 genetic variants associated with height) in 10,226 CRC cases and 10,286 controls. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (95% CI) for associations between height, genetically predicted height and CRC.

RESULTS

Using conventional methods, increased height (per 10-cm increment) was associated with increased CRC risk (OR = 1.08, 95% CI = 1.02-1.15). In sex-specific analyses, height was associated with CRC risk for women (OR = 1.15, 95% CI = 1.05-1.26), but not men (OR = 0.98, 95% CI = 0.92-1.05). Consistent with these results, carrying greater numbers of (weighted) height-increasing alleles (per 1-unit increase) was associated with higher CRC risk for women and men combined (OR = 1.07, 95% CI = 1.01-1.14) and for women (OR = 1.09, 95% CI =  .01-1.19). There was weaker evidence of an association for men (OR = 1.05, 95% CI = 0.96-1.15).

CONCLUSION

We provide evidence for a causal association between height and CRC for women. The CRC-height association for men remains unclear and warrants further investigation in other large studies.

Further Confirmation of Germline Glioma Risk Variant rs78378222 in TP53 and Its Implication in Tumor Tissues via Integrative Analysis of TCGA Data

We confirmed strong association of rs78378222:A>C (per allele odds ratio [OR] = 3.14; P = 6.48 × 10(-11) ), a germline rare single-nucleotide polymorphism (SNP) in TP53, via imputation of a genome-wide association study of glioma (1,856 cases and 4,955 controls). We subsequently performed integrative analyses on the Cancer Genome Atlas (TCGA) data for GBM (glioblastoma multiforme) and LUAD (lung adenocarcinoma). Based on SNP data, we imputed genotypes for rs78378222 and selected individuals carrying rare risk allele (C). Using RNA sequencing data, we observed aberrant transcripts with ∼3 kb longer than normal for those individuals. Using exome sequencing data, we further showed that loss of haplotype carrying common protective allele (A) occurred somatically in GBM but not in LUAD. Our bioinformatic analysis suggests rare risk allele (C) disrupts mRNA termination, and an allelic loss of a genomic region harboring common protective allele (A) occurs during tumor initiation or progression for glioma.

Mendelian Randomization Study of Body Mass Index and Colorectal Cancer Risk

BACKGROUND

High body mass index (BMI) is consistently linked to increased risk of colorectal cancer for men, whereas the association is less clear for women. As risk estimates from observational studies may be biased and/or confounded, we conducted a Mendelian randomization study to estimate the causal association between BMI and colorectal cancer.

METHODS

We used data from 10,226 colorectal cancer cases and 10,286 controls of European ancestry. The Mendelian randomization analysis used a weighted genetic risk score, derived from 77 genome-wide association study-identified variants associated with higher BMI, as an instrumental variable (IV). We compared the IV odds ratio (IV-OR) with the OR obtained using a conventional covariate-adjusted analysis.

RESULTS

Individuals carrying greater numbers of BMI-increasing alleles had higher colorectal cancer risk [per weighted allele OR, 1.31; 95% confidence interval (CI), 1.10-1.57]. Our IV estimation results support the hypothesis that genetically influenced BMI is directly associated with risk for colorectal cancer (IV-OR per 5 kg/m(2), 1.50; 95% CI, 1.13-2.01). In the sex-specific IV analyses higher BMI was associated with higher risk of colorectal cancer among women (IV-OR per 5 kg/m(2), 1.82; 95% CI, 1.26-2.61). For men, genetically influenced BMI was not associated with colorectal cancer (IV-OR per 5 kg/m(2), 1.18; 95% CI, 0.73-1.92).

CONCLUSIONS

High BMI was associated with increased colorectal cancer risk for women. Whether abdominal obesity, rather than overall obesity, is a more important risk factor for men requires further investigation.

IMPACT

Overall, conventional epidemiologic and Mendelian randomization studies suggest a strong association between obesity and the risk of colorectal cancer.

Identification of two metallothioneins as novel inhalative coffee allergens cof a 2 and cof a 3

Background

Dust of green coffee beans is known to be a relevant cause for occupational allergic disorders in coffee industry workers. Recently, we described the first coffee allergen (Cof a 1) establishing an allergenic potential of green coffee dust.

Objective

Our aim was to identify allergenic components of green coffee in order to enhance inhalative coffee allergy diagnosis.

Methods

A Coffea arabica pJuFo cDNA phage display library was created and screened for IgE binding with sera from allergic coffee workers. Two further coffee allergens were identified by sequence analysis, expressed in E. coli, and evaluated by Western blots. The prevalence of sensitization to recombinant Cof a 1, Cof a 2, and Cof a 3 and to commercially available extract was investigated by ELISA (enzyme-linked immunosorbent assay) respectively CAP (capacity test) screening in 18 sera of symptomatic coffee workers.

Results

In addition to the previously described chitinase Cof a 1, two Coffea arabica cysteine-rich metallothioneins of 9 and 7 kDa were identified and included in the IUIS Allergen Nomenclature as Cof a 2 and Cof a 3. Serum IgE antibodies to at least one of the recombinant allergens were found in 8 out of 18 symptomatic coffee workers (44%). Only 2 of the analysed sera (11%) had reacted previously to the commercial allergy test.

Conclusions

In addition to the previously described Cof a 1 we have identified two further coffee proteins to be type I coffee allergens (Cof a 2 and Cof a 3) which may have a relevant potential for the specific diagnosis and/or therapy of coffee allergy.

Proceedings of the second international molecular pathological epidemiology (MPE) meeting

Disease classification system increasingly incorporates information on pathogenic mechanisms to predict clinical outcomes and response to therapy and intervention. Technological advancements to interrogate omics (genomics, epigenomics, transcriptomics, proteomics, metabolomics, metagenomics, interactomics, etc.) provide widely open opportunities in population-based research. Molecular pathological epidemiology (MPE) represents integrative science of molecular pathology and epidemiology. This unified paradigm requires multidisciplinary collaboration between pathology, epidemiology, biostatistics, bioinformatics, and computational biology. Integration of these fields enables better understanding of etiologic heterogeneity, disease continuum, causal inference, and the impact of environment, diet, lifestyle, host factors (including genetics and immunity), and their interactions on disease evolution. Hence, the Second International MPE Meeting was held in Boston in December 2014, with aims to: (1) develop conceptual and practical frameworks; (2) cultivate and expand opportunities; (3) address challenges; and (4) initiate the effort of specifying guidelines for MPE. The meeting mainly consisted of presentations of method developments and recent data in various malignant neoplasms and tumors (breast, prostate, ovarian and colorectal cancers, renal cell carcinoma, lymphoma, and leukemia), followed by open discussion sessions on challenges and future plans. In particular, we recognized need for efforts to further develop statistical methodologies. This meeting provided an unprecedented opportunity for interdisciplinary collaboration, consistent with the purposes of the Big Data to Knowledge, Genetic Associations and Mechanisms in Oncology, and Precision Medicine Initiative of the US National Institute of Health. The MPE meeting series can help advance transdisciplinary population science and optimize training and education systems for twenty-first century medicine and public health.

Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.