Adult height is associated with increased risk of ovarian cancer: a Mendelian randomisation study

BACKGROUND

Observational studies suggest greater height is associated with increased ovarian cancer risk, but cannot exclude bias and/or confounding as explanations for this. Mendelian randomisation (MR) can provide evidence which may be less prone to bias.

METHODS

We pooled data from 39 Ovarian Cancer Association Consortium studies (16,395 cases; 23,003 controls). We applied two-stage predictor-substitution MR, using a weighted genetic risk score combining 609 single-nucleotide polymorphisms. Study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted height and risk were pooled using random-effects meta-analysis.

RESULTS

Greater genetically predicted height was associated with increased ovarian cancer risk overall (pooled-OR (pOR) = 1.06; 95% CI: 1.01-1.11 per 5 cm increase in height), and separately for invasive (pOR = 1.06; 95% CI: 1.01-1.11) and borderline (pOR = 1.15; 95% CI: 1.02-1.29) tumours.

CONCLUSIONS

Women with a genetic propensity to being taller have increased risk of ovarian cancer. This suggests genes influencing height are involved in pathways promoting ovarian carcinogenesis.

MyD88 and TLR4 Expression in Epithelial Ovarian Cancer

OBJECTIVE

To evaluate myeloid differentiation primary response gene 88 (MyD88) and Toll-like receptor 4 (TLR4) expression in relation to clinical features of epithelial ovarian cancer, histologic subtypes, and overall survival.

PATIENTS AND METHODS

We conducted centralized immunohistochemical staining, semi-quantitative scoring, and survival analysis in 5263 patients participating in the Ovarian Tumor Tissue Analysis consortium. Patients were diagnosed between January 1, 1978, and December 31, 2014, including 2865 high-grade serous ovarian carcinomas (HGSOCs), with more than 12,000 person-years of follow-up time. Tissue microarrays were stained for MyD88 and TLR4, and staining intensity was classified using a 2-tiered system for each marker (weak vs strong).

RESULTS

Expression of MyD88 and TLR4 was similar in all histotypes except clear cell ovarian cancer, which showed reduced expression compared with other histotypes (P<.001 for both). In HGSOC, strong MyD88 expression was modestly associated with shortened overall survival (hazard ratio [HR], 1.13; 95% CI, 1.01-1.26; P=.04) but was also associated with advanced stage (P<.001). The expression of TLR4 was not associated with survival. In low-grade serous ovarian cancer (LGSOC), strong expression of both MyD88 and TLR4 was associated with favorable survival (HR [95% CI], 0.49 [0.29-0.84] and 0.44 [0.21-0.89], respectively; P=.009 and P=.02, respectively).

CONCLUSION

Results are consistent with an association between strong MyD88 staining and advanced stage and poorer survival in HGSOC and demonstrate correlation between strong MyD88 and TLR4 staining and improved survival in LGSOC, highlighting the biological differences between the 2 serous histotypes.

Circulating sex hormones in relation to anthropometric, sociodemographic and behavioural factors in an international dataset of 12,300 men

INTRODUCTION

Sex hormones have been implicated in the etiology of a number of diseases. To better understand disease etiology and the mechanisms of disease-risk factor associations, this analysis aimed to investigate the associations of anthropometric, sociodemographic and behavioural factors with a range of circulating sex hormones and sex hormone-binding globulin.

METHODS

Statistical analyses of individual participant data from 12,330 male controls aged 25-85 years from 25 studies involved in the Endogenous Hormones Nutritional Biomarkers and Prostate Cancer Collaborative Group. Analysis of variance was used to estimate geometric means adjusted for study and relevant covariates.

RESULTS

Older age was associated with higher concentrations of sex hormone-binding globulin and dihydrotestosterone and lower concentrations of dehydroepiandrosterone sulfate, free testosterone, androstenedione, androstanediol glucuronide and free estradiol. Higher body mass index was associated with higher concentrations of free estradiol, androstanediol glucuronide, estradiol and estrone and lower concentrations of dihydrotestosterone, testosterone, sex hormone-binding globulin, free testosterone, androstenedione and dehydroepiandrosterone sulfate. Taller height was associated with lower concentrations of androstenedione, testosterone, free testosterone and sex hormone-binding globulin and higher concentrations of androstanediol glucuronide. Current smoking was associated with higher concentrations of androstenedione, sex hormone-binding globulin and testosterone. Alcohol consumption was associated with higher concentrations of dehydroepiandrosterone sulfate, androstenedione and androstanediol glucuronide. East Asians had lower concentrations of androstanediol glucuronide and African Americans had higher concentrations of estrogens. Education and marital status were modestly associated with a small number of hormones.

CONCLUSION

Circulating sex hormones in men are strongly associated with age and body mass index, and to a lesser extent with smoking status and alcohol consumption.

Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer

Importance

Cytotoxic CD8+ tumor-infiltrating lymphocytes (TILs) participate in immune control of epithelial ovarian cancer; however, little is known about prognostic patterns of CD8+ TILs by histotype and in relation to other clinical factors.

Objective

To define the prognostic role of CD8+ TILs in epithelial ovarian cancer.

Design, Setting, and Participants

This was a multicenter observational, prospective survival cohort study of the Ovarian Tumor Tissue Analysis Consortium. More than 5500 patients, including 3196 with high-grade serous ovarian carcinomas (HGSOCs), were followed prospectively for over 24 650 person-years.

Exposures

Following immunohistochemical analysis, CD8+ TILs were identified within the epithelial components of tumor islets. Patients were grouped based on the estimated number of CD8+ TILs per high-powered field: negative (none), low (1-2), moderate (3-19), and high (≥20). CD8+ TILs in a subset of patients were also assessed in a quantitative, uncategorized manner, and the functional form of associations with survival was assessed using penalized B-splines.

Main Outcomes and Measures

Overall survival time.

Results

The final sample included 5577 women; mean age at diagnosis was 58.4 years (median, 58.2 years). Among the 5 major invasive histotypes, HGSOCs showed the most infiltration. CD8+ TILs in HGSOCs were significantly associated with longer overall survival; median survival was 2.8 years for patients with no CD8+ TILs and 3.0 years, 3.8 years, and 5.1 years for patients with low, moderate, or high levels of CD8+ TILs, respectively (P value for trend = 4.2 × 10−16). A survival benefit was also observed among women with endometrioid and mucinous carcinomas, but not for those with the other histotypes. Among HGSOCs, CD8+ TILs were favorable regardless of extent of residual disease following cytoreduction, known standard treatment, and germline BRCA1 pathogenic mutation, but were not prognostic for BRCA2 mutation carriers. Evaluation of uncategorized CD8+ TIL counts showed a near-log-linear functional form.

Conclusions and Relevance

This study demonstrates the histotype-specific nature of immune infiltration and provides definitive evidence for a dose-response relationship between CD8+ TILs and HGSOC survival. That the extent of infiltration is prognostic, not merely its presence or absence, suggests that understanding factors that drive infiltration will be the key to unraveling outcome heterogeneity in this cancer.

Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer

Importance

Cytotoxic CD8+ tumor-infiltrating lymphocytes (TILs) participate in immune control of epithelial ovarian cancer; however, little is known about prognostic patterns of CD8+ TILs by histotype and in relation to other clinical factors.

Objective

To define the prognostic role of CD8+ TILs in epithelial ovarian cancer.

Design, Setting, and Participants

This was a multicenter observational, prospective survival cohort study of the Ovarian Tumor Tissue Analysis Consortium. More than 5500 patients, including 3196 with high-grade serous ovarian carcinomas (HGSOCs), were followed prospectively for over 24 650 person-years.

Exposures

Following immunohistochemical analysis, CD8+ TILs were identified within the epithelial components of tumor islets. Patients were grouped based on the estimated number of CD8+ TILs per high-powered field: negative (none), low (1-2), moderate (3-19), and high (≥20). CD8+ TILs in a subset of patients were also assessed in a quantitative, uncategorized manner, and the functional form of associations with survival was assessed using penalized B-splines.

Main Outcomes and Measures

Overall survival time.

Results

The final sample included 5577 women; mean age at diagnosis was 58.4 years (median, 58.2 years). Among the 5 major invasive histotypes, HGSOCs showed the most infiltration. CD8+ TILs in HGSOCs were significantly associated with longer overall survival; median survival was 2.8 years for patients with no CD8+ TILs and 3.0 years, 3.8 years, and 5.1 years for patients with low, moderate, or high levels of CD8+ TILs, respectively (P value for trend = 4.2 × 10−16). A survival benefit was also observed among women with endometrioid and mucinous carcinomas, but not for those with the other histotypes. Among HGSOCs, CD8+ TILs were favorable regardless of extent of residual disease following cytoreduction, known standard treatment, and germline BRCA1 pathogenic mutation, but were not prognostic for BRCA2 mutation carriers. Evaluation of uncategorized CD8+ TIL counts showed a near-log-linear functional form.

Conclusions and Relevance

This study demonstrates the histotype-specific nature of immune infiltration and provides definitive evidence for a dose-response relationship between CD8+ TILs and HGSOC survival. That the extent of infiltration is prognostic, not merely its presence or absence, suggests that understanding factors that drive infiltration will be the key to unraveling outcome heterogeneity in this cancer.

Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer

Importance

Cytotoxic CD8+ tumor-infiltrating lymphocytes (TILs) participate in immune control of epithelial ovarian cancer; however, little is known about prognostic patterns of CD8+ TILs by histotype and in relation to other clinical factors.

Objective

To define the prognostic role of CD8+ TILs in epithelial ovarian cancer.

Design, Setting, and Participants

This was a multicenter observational, prospective survival cohort study of the Ovarian Tumor Tissue Analysis Consortium. More than 5500 patients, including 3196 with high-grade serous ovarian carcinomas (HGSOCs), were followed prospectively for over 24 650 person-years.

Exposures

Following immunohistochemical analysis, CD8+ TILs were identified within the epithelial components of tumor islets. Patients were grouped based on the estimated number of CD8+ TILs per high-powered field: negative (none), low (1-2), moderate (3-19), and high (≥20). CD8+ TILs in a subset of patients were also assessed in a quantitative, uncategorized manner, and the functional form of associations with survival was assessed using penalized B-splines.

Main Outcomes and Measures

Overall survival time.

Results

The final sample included 5577 women; mean age at diagnosis was 58.4 years (median, 58.2 years). Among the 5 major invasive histotypes, HGSOCs showed the most infiltration. CD8+ TILs in HGSOCs were significantly associated with longer overall survival; median survival was 2.8 years for patients with no CD8+ TILs and 3.0 years, 3.8 years, and 5.1 years for patients with low, moderate, or high levels of CD8+ TILs, respectively (P value for trend = 4.2 × 10−16). A survival benefit was also observed among women with endometrioid and mucinous carcinomas, but not for those with the other histotypes. Among HGSOCs, CD8+ TILs were favorable regardless of extent of residual disease following cytoreduction, known standard treatment, and germline BRCA1 pathogenic mutation, but were not prognostic for BRCA2 mutation carriers. Evaluation of uncategorized CD8+ TIL counts showed a near-log-linear functional form.

Conclusions and Relevance

This study demonstrates the histotype-specific nature of immune infiltration and provides definitive evidence for a dose-response relationship between CD8+ TILs and HGSOC survival. That the extent of infiltration is prognostic, not merely its presence or absence, suggests that understanding factors that drive infiltration will be the key to unraveling outcome heterogeneity in this cancer.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

Susceptibility Loci-Associated Cutaneous Squamous Cell Carcinoma Invasiveness

Genome-wide association studies have identified genetic loci associated with cutaneous squamous cell carcinoma (cSCC) risk, but single-nucleotide polymorphism associations with cSCC invasiveness have not been investigated. We examined associations between cSCC invasiveness and 23 reported single-nucleotide polymorphisms among 67,833 non-Hispanic white subjects. Additionally, we performed a genome-wide scan and identified one SNP with significantly different frequencies in 5,724 subjects with at least one invasive tumor and 1,943 subjects with in situ tumors only. We then compared genotype frequencies among the invasive and in situ groups with those of 60,166 control subjects. The genome-wide scan identified that the T allele in single-nucleotide polymorphism rs41269979 in the class II human leukocyte antigen region was more frequent in the invasive than the in situ group (P = 4.93 × 10^-8). Single-nucleotide polymorphisms in five of the 23 previously associated loci showed odds ratio heterogeneity between the in situ and invasive groups: rs447510 in HLA-DQA1 (P_het = 2.93 × 10^-3), rs12203592 in IRF4 (P_het = 3.94 × 10^-4), rs1805007 in MC1R (P_het = 7.71 × 10^-3), and two SNPs in DEF8 (rs4268748, P_het = 1.09 × 10^-4 and rs8063761, P_het = 1.40 × 10^-4). These findings may provide new insight into the genetic basis of cSCC invasiveness and may help identify individuals at higher risk for developing clinically aggressive cSCC.

Evaluating disease prediction models using a cohort whose covariate distribution differs from that of the target population

Personal predictive models for disease development play important roles in chronic disease prevention. The performance of these models is evaluated by applying them to the baseline covariates of participants in external cohort studies, with model predictions compared to subjects' subsequent disease incidence. However, the covariate distribution among participants in a validation cohort may differ from that of the population for which the model will be used. Since estimates of predictive model performance depend on the distribution of covariates among the subjects to which it is applied, such differences can cause misleading estimates of model performance in the target population. We propose a method for addressing this problem by weighting the cohort subjects to make their covariate distribution better match that of the target population. Simulations show that the method provides accurate estimates of model performance in the target population, while un-weighted estimates may not. We illustrate the method by applying it to evaluate an ovarian cancer prediction model targeted to US women, using cohort data from participants in the California Teachers Study. The methods can be implemented using open-source code for public use as the R-package RMAP (Risk Model Assessment Package) available at http://stanford.edu/~ggong/rmap/ .

Germline whole exome sequencing and large-scale replication identifies FANCM as a likely high grade serous ovarian cancer susceptibility gene

We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8×10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 - were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 - 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.

Age at Menarche and Late Adolescent Adiposity Associated with Mammographic Density on Processed Digital Mammograms in 24,840 Women

Background: High mammographic density is strongly associated with increased breast cancer risk. Some, but not all, risk factors for breast cancer are also associated with higher mammographic density.Methods: The study cohort (N = 24,840) was drawn from the Research Program in Genes, Environment and Health of Kaiser Permanente Northern California and included non-Hispanic white females ages 40 to 74 years with a full-field digital mammogram (FFDM). Percent density (PD) and dense area (DA) were measured by a radiological technologist using Cumulus. The association of age at menarche and late adolescent body mass index (BMI) with PD and DA were modeled using linear regression adjusted for confounders.Results: Age at menarche and late adolescent BMI were negatively correlated. Age at menarche was positively associated with PD (P value for trend <0.0001) and DA (P value for trend <0.0001) in fully adjusted models. Compared with the reference category of ages 12 to 13 years at menarche, menarche at age >16 years was associated with an increase in PD of 1.47% (95% CI, 0.69-2.25) and an increase in DA of 1.59 cm² (95% CI, 0.48-2.70). Late adolescent BMI was inversely associated with PD (P < 0.0001) and DA (P < 0.0001) in fully adjusted models.Conclusions: Age at menarche and late adolescent BMI are both associated with Cumulus measures of mammographic density on processed FFDM images.Impact: Age at menarche and late adolescent BMI may act through different pathways. The long-term effects of age at menarche on cancer risk may be mediated through factors besides mammographic density. Cancer Epidemiol Biomarkers Prev; 26(9); 1450-8. ©2017 AACR.

Pubertal development in girls by breast cancer family history: the LEGACY girls cohort

Background

Pubertal milestones, such as onset of breast development and menstruation, play an important role in breast cancer etiology. It is unclear if these milestones are different in girls with a first- or second-degree breast cancer family history (BCFH).

Methods

In the LEGACY Girls Study (n = 1040), we examined whether three mother/guardian-reported pubertal milestones (having reached Tanner Stage 2 or higher (T2+) for breast and pubic hair development, and having started menstruation) differed by BCFH. We also examined whether associations between body size and race/ethnicity and pubertal milestones were modified by BCFH. We used mother/guardian reports as the primary measure of pubertal milestones, but also conducted sensitivity analyses using clinical Tanner measurements available for a subcohort (n = 204). We analyzed cross-sectional baseline data with logistic regression models for the entire cohort, and longitudinal data with Weibull survival models for the subcohort of girls that were aged 5–7 years at baseline (n = 258).

Results

BCFH was modestly, but not statistically significantly, associated with Breast T2+ (odds ratio (OR) = 1.36, 95% confidence interval (CI) = 0.88–2.10), with a stronger association seen in the subcohort of girls with clinical breast Tanner staging (OR = 2.20, 95% CI = 0.91–5.32). In a longitudinal analysis of girls who were aged 5–7 years at baseline, BCFH was associated with a 50% increased rate of having early breast development (hazard ratio (HR) = 1.49, 95% CI = 1.0–2.21). This association increased to twofold in girls who were not overweight at baseline (HR = 2.04, 95% CI = 1.29–3.21). BCFH was not associated with pubic hair development and post-menarche status. The median interval between onset of breast development and menarche was longer for BCFH+ than BCFH– girls (2.3 versus 1.7 years), suggesting a slower developmental tempo for BCFH+ girls. Associations between pubertal milestones and body size and race/ethnicity were similar in girls with or without a BCFH. For example, weight was positively associated with Breast T2+ in both girls with (OR = 1.06 per 1 kg, 95% CI = 1.03–1.10) and without (OR = 1.14 per 1 kg, 95% CI = 1.04–1.24) a BCFH.

Conclusions

These results suggest that BCFH may be related to earlier breast development and slower pubertal tempo independent of body size and race/ethnicity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0849-y) contains supplementary material, which is available to authorized users.

Genetic modifiers of CHEK2*1100delC-associated breast cancer risk

PURPOSE

CHEK2*1100delC is a founder variant in European populations that confers a two- to threefold increased risk of breast cancer (BC). Epidemiologic and family studies have suggested that the risk associated with CHEK2*1100delC is modified by other genetic factors in a multiplicative fashion. We have investigated this empirically using data from the Breast Cancer Association Consortium (BCAC).

METHODS

Using genotype data from 39,139 (624 1100delC carriers) BC patients and 40,063 (224) healthy controls from 32 BCAC studies, we analyzed the combined risk effects of CHEK2*1100delC and 77 common variants in terms of a polygenic risk score (PRS) and pairwise interaction.

RESULTS

The PRS conferred odds ratios (OR) of 1.59 (95% CI: 1.21-2.09) per standard deviation for BC for CHEK2*1100delC carriers and 1.58 (1.55-1.62) for noncarriers. No evidence of deviation from the multiplicative model was found. The OR for the highest quintile of the PRS was 2.03 (0.86-4.78) for CHEK2*1100delC carriers, placing them in the high risk category according to UK NICE guidelines. The OR for the lowest quintile was 0.52 (0.16-1.74), indicating a lifetime risk close to the population average.

CONCLUSION

Our results confirm the multiplicative nature of risk effects conferred by CHEK2*1100delC and the common susceptibility variants. Furthermore, the PRS could identify carriers at a high lifetime risk for clinical actions.Genet Med advance online publication 06 October 2016.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

No Evidence That Genetic Variation in the Myeloid-Derived Suppressor Cell Pathway Influences Ovarian Cancer Survival

Background: The precise mechanism by which the immune system is adversely affected in cancer patients remains poorly understood, but the accumulation of immunosuppressive/protumorigenic myeloid-derived suppressor cells (MDSCs) is thought to be a prominent mechanism contributing to immunologic tolerance of malignant cells in epithelial ovarian cancer (EOC). To this end, we hypothesized genetic variation in MDSC pathway genes would be associated with survival after EOC diagnoses.Methods: We measured the hazard of death due to EOC within 10 years of diagnosis, overall and by invasive subtype, attributable to SNPs in 24 genes relevant in the MDSC pathway in 10,751 women diagnosed with invasive EOC. Versatile Gene-based Association Study and the admixture likelihood method were used to test gene and pathway associations with survival.Results: We did not identify individual SNPs that were significantly associated with survival after correction for multiple testing (P < 3.5 × 10-5), nor did we identify significant associations between the MDSC pathway overall, or the 24 individual genes and EOC survival.Conclusions: In this well-powered analysis, we observed no evidence that inherited variations in MDSC-associated SNPs, individual genes, or the collective genetic pathway contributed to EOC survival outcomes.Impact: Common inherited variation in genes relevant to MDSCs was not associated with survival in women diagnosed with invasive EOC. Cancer Epidemiol Biomarkers Prev; 26(3); 420-4. ©2016 AACR.

gsSKAT: Rapid gene set analysis and multiple testing correction for rare-variant association studies using weighted linear kernels

Next-generation sequencing technologies have afforded unprecedented characterization of low-frequency and rare genetic variation. Due to low power for single-variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel-machine regression and adaptive testing methods for aggregative rare-variant association testing have been demonstrated to be powerful approaches for pathway-level analysis, although these methods tend to be computationally intensive at high-variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare-variant analysis using component gene-level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family-wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case-control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open-source R code for public use to facilitate easy application of our methods to existing rare-variant analysis results.

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

BACKGROUND

Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis.

METHODS

All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals).

RESULTS

The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10-5 (including six with P<5 × 10-8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation.

CONCLUSIONS

Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC.

Cutaneous squamous cell cancer (cSCC) risk and the human leukocyte antigen (HLA) system

Cutaneous squamous cell carcinoma (cSCC) is the second most common cancer among Caucasians in the United States, with rising incidence over the past decade. Treatment for non-melanoma skin cancer, including cSCC, in the United States was estimated to cost $4.8 billion in 2014. Thus, an understanding of cSCC pathogenesis could have important public health implications. Immune function impacts cSCC risk, given that cSCC incidence rates are substantially higher in patients with compromised immune systems. We report a systematic review of published associations between cSCC risk and the human leukocyte antigen (HLA) system. This review includes studies that analyze germline class I and class II HLA allelic variation as well as HLA cell-surface protein expression levels associated with cSCC risk. We propose biological mechanisms for these HLA-cSCC associations based on known mechanisms of HLA involvement in other diseases. The review suggests that immunity regulates the development of cSCC and that HLA-cSCC associations differ between immunocompetent and immunosuppressed patients. This difference may reflect the presence of viral co-factors that affect tumorigenesis in immunosuppressed patients. Finally, we highlight limitations in the literature on HLA-cSCC associations, and suggest directions for future research aimed at understanding, preventing and treating cSCC.

Predicting Prostate Cancer Recurrence After Radical Prostatectomy

BACKGROUND

Prostate cancer prognosis is variable, and management decisions involve balancing patients' risks of recurrence and recurrence-free death. Moreover, the roles of body mass index (BMI) and race in risk of recurrence are controversial [1,2]. To address these issues, we developed and cross-validated RAPS (Risks After Prostate Surgery), a personal prediction model for biochemical recurrence (BCR) within 10 years of radical prostatectomy (RP) that includes BMI and race as possible predictors, and recurrence-free death as a competing risk.

METHODS

RAPS uses a patient's risk factors at surgery to assign him a recurrence probability based on statistical learning methods applied to a cohort of 1,276 patients undergoing RP at the University of Pennsylvania. We compared the performance of RAPS to that of an existing model with respect to calibration (by comparing observed and predicted outcomes), and discrimination (using the area under the receiver operating characteristic curve (AUC)).

RESULTS

RAPS' cross-validated BCR predictions provided better calibration than those of an existing model that underestimated patients' risks. Discrimination was similar for the two models, with BCR AUCs of 0.793, 95% confidence interval (0.766-0.820) for RAPS, and 0.780 (0.745-0.815) for the existing model. RAPS' most important BCR predictors were tumor grade, preoperative prostate-specific antigen (PSA) level and BMI; race was less important [3]. RAPS' predictions can be obtained online at https://predict.shinyapps.io/raps.

CONCLUSION

RAPS' cross-validated BCR predictions were better calibrated than those of an existing model, and BMI information contributed substantially to these predictions. RAPS predictions for recurrence-free death were limited by lack of co-morbidity data; however the model provides a simple framework for extension to include such data. Its use and extension should facilitate decision strategies for post-RP prostate cancer management. Prostate 77:291-298, 2017. © 2016 Wiley Periodicals, Inc.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Does radiotherapy still have a role in unresected biliary tract cancer?

The benefits of radiotherapy for inoperable biliary tract cancer remain unclear due to the lack of randomized data. We evaluated the impact of radiotherapy on survival in elderly patients using the SEER‐Medicare database. Patients in the SEER‐Medicare database with inoperable biliary tract tumors diagnosed between 1998 and 2011 were included. We used multivariate logistic regression to evaluate factors associated with treatment selection, and multivariate Cox regression and propensity score matching to evaluate treatment selection in relation to subsequent survival. Of the 2343 patients included, 451 (19%) received radiotherapy within 4 months of diagnosis. The use of radiotherapy declined over time, and was influenced by receipt of chemotherapy and patient age, race, marital status, poverty status, and tumor stage and type. Median survival was 9.3 (95% CI 8.7–9.7) months among patients who did not receive radiation and 10.0 (95% CI 9.1–11.3) months among those who received radiation, conditional on having survived 4 months. In patients who received chemotherapy (n = 1053), receipt of radiation was associated with improved survival, with an adjusted hazard ratio of 0.82 (95% 0.70–0.97, P = 0.02). In patients who did not receive chemotherapy (n = 1290), receipt of radiation was not associated with improved survival, with an adjusted hazard ratio of 1.09 (95% 0.91–1.30, P = 0.34). Propensity‐scored matched analyses showed similar results. Despite the survival benefit associated with the addition of radiotherapy to chemotherapy, the use of radiation for unresectable biliary tract cancers has declined over time.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with risk of clear cell ovarian cancer

BACKGROUND

Regulatory T (Treg) cells, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and, thus, variants in genes encoding Treg cell immune molecules could be associated with ovarian cancer.

METHODS

In a population of 15,596 epithelial ovarian cancer (EOC) cases and 23,236 controls, we measured genetic associations of 1,351 SNPs in Treg cell pathway genes with odds of ovarian cancer and tested pathway and gene-level associations, overall and by histotype, for the 25 genes, using the admixture likelihood (AML) method. The most significant single SNP associations were tested for correlation with expression levels in 44 ovarian cancer patients.

RESULTS

The most significant global associations for all genes in the pathway were seen in endometrioid ( p = 0.082) and clear cell ( p = 0.083), with the most significant gene level association seen with TGFBR2 ( p = 0.001) and clear cell EOC. Gene associations with histotypes at p < 0.05 included: IL12 ( p = 0.005 and p = 0.008, serous and high-grade serous, respectively), IL8RA ( p = 0.035, endometrioid and mucinous), LGALS1 ( p = 0.03, mucinous), STAT5B ( p = 0.022, clear cell), TGFBR1 ( p = 0.021 endometrioid) and TGFBR2 ( p = 0.017 and p = 0.025, endometrioid and mucinous, respectively).

CONCLUSIONS

Common inherited gene variation in Treg cell pathways shows some evidence of germline genetic contribution to odds of EOC that varies by histologic subtype and may be associated with mRNA expression of immune-complex receptor in EOC patients.