1
|
Nieuwoudt C, Farooq FB, Brooks-Wilson A, Bureau A, Graham J. Statistics to prioritize rare variants in family-based sequencing studies with disease subtypes. Genet Epidemiol 2024. [PMID: 38940260 DOI: 10.1002/gepi.22579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 03/26/2024] [Accepted: 06/13/2024] [Indexed: 06/29/2024]
Abstract
Family-based sequencing studies are increasingly used to find rare genetic variants of high risk for disease traits with familial clustering. In some studies, families with multiple disease subtypes are collected and the exomes of affected relatives are sequenced for shared rare variants (RVs). Since different families can harbor different causal variants and each family harbors many RVs, tests to detect causal variants can have low power in this study design. Our goal is rather to prioritize shared variants for further investigation by, for example, pathway analyses or functional studies. The transmission-disequilibrium test prioritizes variants based on departures from Mendelian transmission in parent-child trios. Extending this idea to families, we propose methods to prioritize RVs shared in affected relatives with two disease subtypes, with one subtype more heritable than the other. Global approaches condition on a variant being observed in the study and assume a known probability of carrying a causal variant. In contrast, local approaches condition on a variant being observed in specific families to eliminate the carrier probability. Our simulation results indicate that global approaches are robust to misspecification of the carrier probability and prioritize more effectively than local approaches even when the carrier probability is misspecified.
Collapse
|
2
|
Saner FAM, Takahashi K, Budden T, Pandey A, Ariyaratne D, Zwimpfer TA, Meagher NS, Fereday S, Twomey L, Pishas KI, Hoang T, Bolithon A, Traficante N, Alsop K, Christie EL, Kang EY, Nelson GS, Ghatage P, Lee CH, Riggan MJ, Alsop J, Beckmann MW, Boros J, Brand AH, Brooks-Wilson A, Carney ME, Coulson P, Courtney-Brooks M, Cushing-Haugen KL, Cybulski C, El-Bahrawy M, Elishaev E, Erber R, Gayther SA, Gentry-Maharaj A, Gilks CB, Harnett PR, Harris HR, Hartmann A, Hein A, Hendley J, Hernandez BY, Jakubowska A, Jimenez-Linan M, Jones ME, Kaufmann SH, Kennedy CJ, Kluz T, Koziak JM, Kristjansdottir B, Le ND, Lener M, Lester J, Lubiński J, Mateoiu C, Orsulic S, Ruebner M, Schoemaker MJ, Shah M, Sharma R, Sherman ME, Shvetsov YB, Soong TR, Steed H, Sukumvanich P, Talhouk A, Taylor SE, Vierkant RA, Wang C, Widschwendter M, Wilkens LR, Winham SJ, Anglesio MS, Berchuck A, Brenton JD, Campbell I, Cook LS, Doherty JA, Fasching PA, Fortner RT, Goodman MT, Gronwald J, Huntsman DG, Karlan BY, Kelemen LE, Menon U, Modugno F, Pharoah PDP, Schildkraut JM, Sundfeldt K, Swerdlow AJ, Goode EL, DeFazio A, Köbel M, Ramus SJ, Bowtell DDL, Garsed DW. Concurrent RB1 Loss and BRCA-Deficiency Predicts Enhanced Immunological Response and Long-Term Survival in Tubo-Ovarian High-Grade Serous Carcinoma. Clin Cancer Res 2024:745770. [PMID: 38837893 DOI: 10.1158/1078-0432.ccr-23-3552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Revised: 02/08/2024] [Accepted: 05/31/2024] [Indexed: 06/07/2024]
Abstract
PURPOSE To evaluate RB1 expression and survival across ovarian carcinoma histotypes, and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC). EXPERIMENTAL DESIGN RB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 primary HGSC to characterize tumors with concurrent BRCA-deficiency and RB1 loss. RESULTS RB1 loss was associated with longer OS in HGSC, but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared to patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA-deficiency correlated with transcriptional markers of enhanced interferon response, cell-cycle deregulation, and reduced epithelial-mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1. CONCLUSIONS Co-occurrence of RB1 loss and BRCA-deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.
Collapse
|
3
|
Vishnyakova O, Song X, Rockwood K, Elliott LT, Brooks-Wilson A. Physiological phenotypes have optimal values relevant to healthy aging: sweet spots deduced from the Canadian Longitudinal Study on Aging. GeroScience 2024; 46:1589-1605. [PMID: 37688655 PMCID: PMC10828371 DOI: 10.1007/s11357-023-00895-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Accepted: 07/27/2023] [Indexed: 09/11/2023] Open
Abstract
Previous observations on a group of exceptionally healthy "Super-Seniors" showed a lower variance of multiple physiological measures relevant for health than did a less healthy group of the same age. The finding was interpreted as the healthier individuals having physiological measurement values closer to an optimal level, or "sweet spot." Here, we tested the generalizability of the sweet-spot hypothesis in a larger community sample, comparing differences in the variance between healthier and less healthy groups. We apply this method to the Canadian Longitudinal Study on Aging (CLSA) comprehensive cohort of 30,097 participants aged 45 to 85 years with deep phenotype data. Data from both sexes and four age ranges were analyzed. Five instruments were used to represent different aspects of health, physical, and cognitive functioning. We tested 231 phenotypic measures for lower variance in the most healthy vs. least healthy quartile of each sex and age group, as classified by the five instruments. Segmented regression was used to determine sex-specific optimal values. One hundred forty-two physiological measures (61%) showed lower variance in the healthiest than in the least healthy group, in at least one sex and age group. The difference in variance was most significant for hemoglobin A1c and was also significant for many body composition measurements, but not for bone mineral density. Ninety-four phenotypes showed a nonmonotonic relationship with health, consistent with the idea of a sweet spot; for these, we determined optimal values and 95% confidence intervals that were generally narrower than the ranges of current clinical reference intervals. These findings for sweet spot discovery validate the proposed approach for identifying traits important for healthy aging.
Collapse
|
4
|
Saner FAM, Takahashi K, Budden T, Pandey A, Ariyaratne D, Zwimpfer TA, Meagher NS, Fereday S, Twomey L, Pishas KI, Hoang T, Bolithon A, Traficante N, Alsop K, Christie EL, Kang EY, Nelson GS, Ghatage P, Lee CH, Riggan MJ, Alsop J, Beckmann MW, Boros J, Brand AH, Brooks-Wilson A, Carney ME, Coulson P, Courtney-Brooks M, Cushing-Haugen KL, Cybulski C, El-Bahrawy MA, Elishaev E, Erber R, Gayther SA, Gentry-Maharaj A, Blake Gilks C, Harnett PR, Harris HR, Hartmann A, Hein A, Hendley J, Hernandez BY, Jakubowska A, Jimenez-Linan M, Jones ME, Kaufmann SH, Kennedy CJ, Kluz T, Koziak JM, Kristjansdottir B, Le ND, Lener M, Lester J, Lubiński J, Mateoiu C, Orsulic S, Ruebner M, Schoemaker MJ, Shah M, Sharma R, Sherman ME, Shvetsov YB, Singh N, Rinda Soong T, Steed H, Sukumvanich P, Talhouk A, Taylor SE, Vierkant RA, Wang C, Widschwendter M, Wilkens LR, Winham SJ, Anglesio MS, Berchuck A, Brenton JD, Campbell I, Cook LS, Doherty JA, Fasching PA, Fortner RT, Goodman MT, Gronwald J, Huntsman DG, Karlan BY, Kelemen LE, Menon U, Modugno F, Pharoah PD, Schildkraut JM, Sundfeldt K, Swerdlow AJ, Goode EL, DeFazio A, Köbel M, Ramus SJ, Bowtell DDL, Garsed DW. Concurrent RB1 loss and BRCA-deficiency predicts enhanced immunological response and long-term survival in tubo-ovarian high-grade serous carcinoma. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.11.09.23298321. [PMID: 37986741 PMCID: PMC10659507 DOI: 10.1101/2023.11.09.23298321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2023]
Abstract
Background Somatic loss of the tumour suppressor RB1 is a common event in tubo-ovarian high-grade serous carcinoma (HGSC), which frequently co-occurs with alterations in homologous recombination DNA repair genes including BRCA1 and BRCA2 (BRCA). We examined whether tumour expression of RB1 was associated with survival across ovarian cancer histotypes (HGSC, endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous carcinoma (LGSC)), and how co-occurrence of germline BRCA pathogenic variants and RB1 loss influences long-term survival in a large series of HGSC. Patients and methods RB1 protein expression patterns were classified by immunohistochemistry in epithelial ovarian carcinomas of 7436 patients from 20 studies participating in the Ovarian Tumor Tissue Analysis consortium and assessed for associations with overall survival (OS), accounting for patient age at diagnosis and FIGO stage. We examined RB1 expression and germline BRCA status in a subset of 1134 HGSC, and related genotype to survival, tumour infiltrating CD8+ lymphocyte counts and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cell lines with and without BRCA1 mutations to model co-loss with treatment response. We also performed genomic analyses on 126 primary HGSC to explore the molecular characteristics of concurrent homologous recombination deficiency and RB1 loss. Results RB1 protein loss was most frequent in HGSC (16.4%) and was highly correlated with RB1 mRNA expression. RB1 loss was associated with longer OS in HGSC (hazard ratio [HR] 0.74, 95% confidence interval [CI] 0.66-0.83, P = 6.8 ×10-7), but with poorer prognosis in ENOC (HR 2.17, 95% CI 1.17-4.03, P = 0.0140). Germline BRCA mutations and RB1 loss co-occurred in HGSC (P < 0.0001). Patients with both RB1 loss and germline BRCA mutations had a superior OS (HR 0.38, 95% CI 0.25-0.58, P = 5.2 ×10-6) compared to patients with either alteration alone, and their median OS was three times longer than non-carriers whose tumours retained RB1 expression (9.3 years vs. 3.1 years). Enhanced sensitivity to cisplatin (P < 0.01) and paclitaxel (P < 0.05) was seen in BRCA1 mutated cell lines with RB1 knockout. Among 126 patients with whole-genome and transcriptome sequence data, combined RB1 loss and genomic evidence of homologous recombination deficiency was correlated with transcriptional markers of enhanced interferon response, cell cycle deregulation, and reduced epithelial-mesenchymal transition in primary HGSC. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1. Conclusions Co-occurrence of RB1 loss and BRCA mutation was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.
Collapse
|
5
|
Middha P, Wang X, Behrens S, Bolla MK, Wang Q, Dennis J, Michailidou K, Ahearn TU, Andrulis IL, Anton-Culver H, Arndt V, Aronson KJ, Auer PL, Augustinsson A, Baert T, Freeman LEB, Becher H, Beckmann MW, Benitez J, Bojesen SE, Brauch H, Brenner H, Brooks-Wilson A, Campa D, Canzian F, Carracedo A, Castelao JE, Chanock SJ, Chenevix-Trench G, Cordina-Duverger E, Couch FJ, Cox A, Cross SS, Czene K, Dossus L, Dugué PA, Eliassen AH, Eriksson M, Evans DG, Fasching PA, Figueroa JD, Fletcher O, Flyger H, Gabrielson M, Gago-Dominguez M, Giles GG, González-Neira A, Grassmann F, Grundy A, Guénel P, Haiman CA, Håkansson N, Hall P, Hamann U, Hankinson SE, Harkness EF, Holleczek B, Hoppe R, Hopper JL, Houlston RS, Howell A, Hunter DJ, Ingvar C, Isaksson K, Jernström H, John EM, Jones ME, Kaaks R, Keeman R, Kitahara CM, Ko YD, Koutros S, Kurian AW, Lacey JV, Lambrechts D, Larson NL, Larsson S, Le Marchand L, Lejbkowicz F, Li S, Linet M, Lissowska J, Martinez ME, Maurer T, Mulligan AM, Mulot C, Murphy RA, Newman WG, Nielsen SF, Nordestgaard BG, Norman A, O'Brien KM, Olson JE, Patel AV, Prentice R, Rees-Punia E, Rennert G, Rhenius V, Ruddy KJ, Sandler DP, Scott CG, Shah M, Shu XO, Smeets A, Southey MC, Stone J, Tamimi RM, Taylor JA, Teras LR, Tomczyk K, Troester MA, Truong T, Vachon CM, Wang SS, Weinberg CR, Wildiers H, Willett W, Winham SJ, Wolk A, Yang XR, Zamora MP, Zheng W, Ziogas A, Dunning AM, Pharoah PDP, García-Closas M, Schmidt MK, Kraft P, Milne RL, Lindström S, Easton DF, Chang-Claude J. A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry. Breast Cancer Res 2023; 25:93. [PMID: 37559094 PMCID: PMC10411002 DOI: 10.1186/s13058-023-01691-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 07/27/2023] [Indexed: 08/11/2023] Open
Abstract
BACKGROUND Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer. METHODS Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs. RESULTS Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94). CONCLUSIONS Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.
Collapse
|
6
|
Fu Z, Brooks MM, Irvin S, Jordan S, Aben KKH, Anton-Culver H, Bandera EV, Beckmann MW, Berchuck A, Brooks-Wilson A, Chang-Claude J, Cook LS, Cramer DW, Cushing-Haugen KL, Doherty JA, Ekici AB, Fasching PA, Fortner RT, Gayther SA, Gentry-Maharaj A, Giles GG, Goode EL, Goodman MT, Harris HR, Hein A, Kaaks R, Kiemeney LA, Köbel M, Kotsopoulos J, Le ND, Lee AW, Matsuo K, McGuire V, McLaughlin JR, Menon U, Milne RL, Moysich KB, Pearce CL, Pike MC, Qin B, Ramus SJ, Riggan MJ, Rothstein JH, Schildkraut JM, Sieh W, Sutphen R, Terry KL, Thompson PJ, Titus L, van Altena AM, White E, Whittemore AS, Wu AH, Zheng W, Ziogas A, Taylor SE, Tang L, Songer T, Wentzensen N, Webb PM, Risch HA, Modugno F. Lifetime ovulatory years and risk of epithelial ovarian cancer: a multinational pooled analysis. J Natl Cancer Inst 2023; 115:539-551. [PMID: 36688720 PMCID: PMC10165492 DOI: 10.1093/jnci/djad011] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Revised: 11/10/2022] [Accepted: 01/06/2023] [Indexed: 01/24/2023] Open
Abstract
BACKGROUND The role of ovulation in epithelial ovarian cancer (EOC) is supported by the consistent protective effects of parity and oral contraceptive use. Whether these factors protect through anovulation alone remains unclear. We explored the association between lifetime ovulatory years (LOY) and EOC. METHODS LOY was calculated using 12 algorithms. Odds ratios (ORs) and 95% confidence intervals (CIs) estimated the association between LOY or LOY components and EOC among 26 204 control participants and 21 267 case patients from 25 studies. To assess whether LOY components act through ovulation suppression alone, we compared beta coefficients obtained from regression models with expected estimates assuming 1 year of ovulation suppression has the same effect regardless of source. RESULTS LOY was associated with increased EOC risk (OR per year increase = 1.014, 95% CI = 1.009 to 1.020 to OR per year increase = 1.044, 95% CI = 1.041 to 1.048). Individual LOY components, except age at menarche, also associated with EOC. The estimated model coefficient for oral contraceptive use and pregnancies were 4.45 times and 12- to 15-fold greater than expected, respectively. LOY was associated with high-grade serous, low-grade serous, endometrioid, and clear cell histotypes (ORs per year increase = 1.054, 1.040, 1.065, and 1.098, respectively) but not mucinous tumors. Estimated coefficients of LOY components were close to expected estimates for high-grade serous but larger than expected for low-grade serous, endometrioid, and clear cell histotypes. CONCLUSIONS LOY is positively associated with nonmucinous EOC. Differences between estimated and expected model coefficients for LOY components suggest factors beyond ovulation underlie the associations between LOY components and EOC in general and for non-HGSOC.
Collapse
|
7
|
Weir A, Kang EY, Meagher NS, Nelson GS, Ghatage P, Lee CH, Riggan MJ, Gentry-Maharaj A, Ryan A, Singh N, Widschwendter M, Alsop J, Anglesio MS, Beckmann MW, Berger J, Bisinotto C, Boros J, Brand AH, Brenton JD, Brooks-Wilson A, Carney ME, Cunningham JM, Cushing-Haugen KL, Cybulski C, Elishaev E, Erber R, Fereday S, Fischer A, Paz-Ares L, Gayarre J, Gilks BC, Grube M, Harnett PR, Harris HR, Hartmann A, Hein A, Hendley J, Hernandez BY, Heublein S, Huang Y, Huzarski T, Jakubowska A, Jimenez-Linan M, Kennedy CJ, Kommoss FKF, Koziak JM, Kraemer B, Le ND, Lesnock J, Lester J, Lubiński J, Menkiszak J, Ney B, Olawaiye A, Orsulic S, Osorio A, Robles-Díaz L, Ruebner M, Shah M, Sharma R, Shvetsov YB, Steed H, Talhouk A, Taylor SE, Traficante N, Vierkant RA, Wang C, Wilkens LR, Winham SJ, Benitez J, Berchuck A, Bowtell DD, Candido Dos Reis FJ, Cook LS, DeFazio A, Doherty JA, Fasching PA, García MJ, Goode EL, Goodman MT, Gronwald J, Huntsman DG, Karlan BY, Kommoss S, Modugno F, Schildkraut JM, Sinn HP, Staebler A, Kelemen LE, Ford CE, Menon U, Pharoah PDP, Köbel M, Ramus SJ. Increased FOXJ1 protein expression is associated with improved overall survival in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis Consortium Study. Br J Cancer 2023; 128:137-147. [PMID: 36323878 DOI: 10.1038/s41416-022-02014-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2022] [Revised: 09/01/2022] [Accepted: 10/04/2022] [Indexed: 11/07/2022] Open
Abstract
BACKGROUND Recently, we showed a >60% difference in 5-year survival for patients with tubo-ovarian high-grade serous carcinoma (HGSC) when stratified by a 101-gene mRNA expression prognostic signature. Given the varied patient outcomes, this study aimed to translate prognostic mRNA markers into protein expression assays by immunohistochemistry and validate their survival association in HGSC. METHODS Two prognostic genes, FOXJ1 and GMNN, were selected based on high-quality antibodies, correlation with protein expression and variation in immunohistochemical scores in a preliminary cohort (n = 134 and n = 80, respectively). Six thousand four hundred and thirty-four (FOXJ1) and 5470 (GMNN) formalin-fixed, paraffin-embedded ovarian neoplasms (4634 and 4185 HGSC, respectively) represented on tissue microarrays from the Ovarian Tumor Tissue Analysis consortium underwent immunohistochemical staining and scoring, then univariate and multivariate survival analysis. RESULTS Consistent with mRNA, FOXJ1 protein expression exhibited a linear, increasing association with improved overall survival in HGSC patients. Women with >50% expression had the most favourable outcomes (HR = 0.78, 95% CI 0.67-0.91, p < 0.0001). GMNN protein expression was not significantly associated with overall HSGC patient survival. However, HGSCs with >35% GMNN expression showed a trend for better outcomes, though this was not significant. CONCLUSION We provide foundational evidence for the prognostic value of FOXJ1 in HGSC, validating the prior mRNA-based prognostic association by immunohistochemistry.
Collapse
|
8
|
Meagher NS, Hamilton P, Milne K, Thornton S, Harris B, Weir A, Alsop J, Bisinoto C, Brenton JD, Brooks-Wilson A, Chiu DS, Cushing-Haugen KL, Fereday S, Garsed DW, Gayther SA, Gentry-Maharaj A, Gilks B, Jimenez-Linan M, Kennedy CJ, Le ND, Piskorz AM, Riggan MJ, Shah M, Singh N, Talhouk A, Widschwendter M, Bowtell DDL, Candido Dos Reis FJ, Cook LS, Fortner RT, García MJ, Harris HR, Huntsman DG, Karnezis AN, Köbel M, Menon U, Pharoah PDP, Doherty JA, Anglesio MS, Pike MC, Pearce CL, Friedlander ML, DeFazio A, Nelson BH, Ramus SJ. Profiling the immune landscape in mucinous ovarian carcinoma. Gynecol Oncol 2023; 168:23-31. [PMID: 36368129 PMCID: PMC10374276 DOI: 10.1016/j.ygyno.2022.10.022] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 10/21/2022] [Accepted: 10/26/2022] [Indexed: 11/09/2022]
Abstract
OBJECTIVE Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer, with low response rates to standard chemotherapy, and very poor survival for patients diagnosed at advanced stage. There is a limited understanding of the MOC immune landscape, and consequently whether immune checkpoint inhibitors could be considered for a subset of patients. METHODS We performed multicolor immunohistochemistry (IHC) and immunofluorescence (IF) on tissue microarrays in a cohort of 126 MOC patients. Cell densities were calculated in the epithelial and stromal components for tumor-associated macrophages (CD68+/PD-L1+, CD68+/PD-L1-), T cells (CD3+/CD8-, CD3+/CD8+), putative T-regulatory cells (Tregs, FOXP3+), B cells (CD20+/CD79A+), plasma cells (CD20-/CD79a+), and PD-L1+ and PD-1+ cells, and compared these values with clinical factors. Univariate and multivariable Cox Proportional Hazards assessed overall survival. Unsupervised k-means clustering identified patient subsets with common patterns of immune cell infiltration. RESULTS Mean densities of PD1+ cells, PD-L1- macrophages, CD4+ and CD8+ T cells, and FOXP3+ Tregs were higher in the stroma compared to the epithelium. Tumors from advanced (Stage III/IV) MOC had greater epithelial infiltration of PD-L1- macrophages, and fewer PD-L1+ macrophages compared with Stage I/II cancers (p = 0.004 and p = 0.014 respectively). Patients with high epithelial density of FOXP3+ cells, CD8+/FOXP3+ cells, or PD-L1- macrophages, had poorer survival, and high epithelial CD79a + plasma cells conferred better survival, all upon univariate analysis only. Clustering showed that most MOC (86%) had an immune depleted (cold) phenotype, with only a small proportion (11/76,14%) considered immune inflamed (hot) based on T cell and PD-L1 infiltrates. CONCLUSION In summary, MOCs are mostly immunogenically 'cold', suggesting they may have limited response to current immunotherapies.
Collapse
|
9
|
Yu Y, Chen S, Jones SJ, Hoque R, Vishnyakova O, Brooks-Wilson A, McNeney B. Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes. Hum Hered 2022; 87:000525650. [PMID: 35767963 DOI: 10.1159/000525650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Accepted: 04/11/2022] [Indexed: 11/19/2022] Open
Abstract
INTRODUCTION Increasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero. For a given m, log-F-penalized logistic regression may be easily implemented using data augmentation and standard software. METHOD We propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F-penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation. RESULTS We evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F-penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls. DISCUSSION/CONCLUSION We have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach.
Collapse
|
10
|
Wang SS, Vajdic CM, Linet MS, Slager SL, Voutsinas J, Nieters A, Casabonne D, Cerhan JR, Cozen W, Alarcón G, Martínez-Maza O, Brown EE, Bracci PM, Turner J, Hjalgrim H, Bhatti P, Zhang Y, Birmann BM, Flowers CR, Paltiel O, Holly EA, Kane E, Weisenburger DD, Maynadié M, Cocco P, Foretova L, Breen EC, Lan Q, Brooks-Wilson A, De Roos AJ, Smith MT, Roman E, Boffetta P, Kricker A, Zheng T, Skibola CF, Clavel J, Monnereau A, Chanock SJ, Rothman N, Benavente Y, Hartge P, Smedby KE. B-Cell NHL Subtype Risk Associated with Autoimmune Conditions and PRS. Cancer Epidemiol Biomarkers Prev 2022; 31:1103-1110. [PMID: 35244686 PMCID: PMC9081255 DOI: 10.1158/1055-9965.epi-21-0875] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 12/02/2021] [Accepted: 02/16/2022] [Indexed: 12/14/2022] Open
Abstract
BACKGROUND A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non-Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes. METHODS In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell-mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression. RESULTS We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08-1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59-2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell-mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction. CONCLUSIONS Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions. IMPACT Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.
Collapse
|
11
|
Wang X, Kapoor PM, Auer PL, Dennis J, Dunning AM, Wang Q, Lush M, Michailidou K, Bolla MK, Aronson KJ, Murphy RA, Brooks-Wilson A, Lee DG, Cordina-Duverger E, Guénel P, Truong T, Mulot C, Teras LR, Patel AV, Dossus L, Kaaks R, Hoppe R, Lo WY, Brüning T, Hamann U, Czene K, Gabrielson M, Hall P, Eriksson M, Jung A, Becher H, Couch FJ, Larson NL, Olson JE, Ruddy KJ, Giles GG, MacInnis RJ, Southey MC, Le Marchand L, Wilkens LR, Haiman CA, Olsson H, Augustinsson A, Krüger U, Wagner P, Scott C, Winham SJ, Vachon CM, Perou CM, Olshan AF, Troester MA, Hunter DJ, Eliassen HA, Tamimi RM, Brantley K, Andrulis IL, Figueroa J, Chanock SJ, Ahearn TU, García-Closas M, Evans GD, Newman WG, van Veen EM, Howell A, Wolk A, Håkansson N, Anton-Culver H, Ziogas A, Jones ME, Orr N, Schoemaker MJ, Swerdlow AJ, Kitahara CM, Linet M, Prentice RL, Easton DF, Milne RL, Kraft P, Chang-Claude J, Lindström S. Genome-wide interaction analysis of menopausal hormone therapy use and breast cancer risk among 62,370 women. Sci Rep 2022; 12:6199. [PMID: 35418701 PMCID: PMC9007944 DOI: 10.1038/s41598-022-10121-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2021] [Accepted: 03/11/2022] [Indexed: 02/02/2023] Open
Abstract
Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10-8 as genome-wide significant, and p-values < 1 × 10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 105. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10-7), which showed genome-wide significant interaction (p-value = 3.8 × 10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT-breast cancer risk association.
Collapse
|
12
|
Chinda B, Tran KH, Doesburg S, Siu W, Medvedev G, Liang SS, Brooks-Wilson A, Song X. Functional MRI evaluation of cognitive effects of carotid stenosis revascularization. Brain Behav 2022; 12:e2512. [PMID: 35233977 PMCID: PMC9014987 DOI: 10.1002/brb3.2512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 12/10/2021] [Accepted: 01/07/2022] [Indexed: 11/12/2022] Open
Abstract
INTRODUCTION Severe internal carotid stenosis, if left untreated, can pose serious risks for ischemic stroke and cognitive impairments. The effects of revascularization on any aspects of cognition, however, are not well understood, as conflicting results are reported, which have mainly been centered on paper-based cognitive analyses. Here, we summarized and evaluated the publications to date of functional MRI (fMRI) studies that examined the mechanisms of functional brain activation and connectivity as a way to reflect cognitive effects of revascularization on patients with carotid stenosis. METHODS A PubMed and Google Scholar (covering the relevant literature until November 1, 2021) search yielded eight original studies of the research line, including seven resting-state and one task-based fMRI reports. RESULTS Findings demonstrated treatment-related alterations in fMRI signal intensity and symmetry level, regional fMRI activation pattern, and functional brain network connectivity. The functional brain changes were associated largely with improvement in cognitive function assessed using standard cognitive test scores. CONCLUSIONS These findings support the contribution of fMRI to the understanding of brain functional activation and connectivity changes revealing cognitive effects of revascularization in the management of severe carotid stenosis. The review also highlighted the importance of reproducibility through enhancing experimental designs and cognitive task applications with future research for potential clinical translation.
Collapse
|
13
|
Ahearn TU, Zhang H, Michailidou K, Milne RL, Bolla MK, Dennis J, Dunning AM, Lush M, Wang Q, Andrulis IL, Anton-Culver H, Arndt V, Aronson KJ, Auer PL, Augustinsson A, Baten A, Becher H, Behrens S, Benitez J, Bermisheva M, Blomqvist C, Bojesen SE, Bonanni B, Børresen-Dale AL, Brauch H, Brenner H, Brooks-Wilson A, Brüning T, Burwinkel B, Buys SS, Canzian F, Castelao JE, Chang-Claude J, Chanock SJ, Chenevix-Trench G, Clarke CL, Collée JM, Cox A, Cross SS, Czene K, Daly MB, Devilee P, Dörk T, Dwek M, Eccles DM, Evans DG, Fasching PA, Figueroa J, Floris G, Gago-Dominguez M, Gapstur SM, García-Sáenz JA, Gaudet MM, Giles GG, Goldberg MS, González-Neira A, Alnæs GIG, Grip M, Guénel P, Haiman CA, Hall P, Hamann U, Harkness EF, Heemskerk-Gerritsen BAM, Holleczek B, Hollestelle A, Hooning MJ, Hoover RN, Hopper JL, Howell A, Jakimovska M, Jakubowska A, John EM, Jones ME, Jung A, Kaaks R, Kauppila S, Keeman R, Khusnutdinova E, Kitahara CM, Ko YD, Koutros S, Kristensen VN, Krüger U, Kubelka-Sabit K, Kurian AW, Kyriacou K, Lambrechts D, Lee DG, Lindblom A, Linet M, Lissowska J, Llaneza A, Lo WY, MacInnis RJ, Mannermaa A, Manoochehri M, Margolin S, Martinez ME, McLean C, Meindl A, Menon U, Nevanlinna H, Newman WG, Nodora J, Offit K, Olsson H, Orr N, Park-Simon TW, Patel AV, Peto J, Pita G, Plaseska-Karanfilska D, Prentice R, Punie K, Pylkäs K, Radice P, Rennert G, Romero A, Rüdiger T, Saloustros E, Sampson S, Sandler DP, Sawyer EJ, Schmutzler RK, Schoemaker MJ, Schöttker B, Sherman ME, Shu XO, Smichkoska S, Southey MC, Spinelli JJ, Swerdlow AJ, Tamimi RM, Tapper WJ, Taylor JA, Teras LR, Terry MB, Torres D, Troester MA, Vachon CM, van Deurzen CHM, van Veen EM, Wagner P, Weinberg CR, Wendt C, Wesseling J, Winqvist R, Wolk A, Yang XR, Zheng W, Couch FJ, Simard J, Kraft P, Easton DF, Pharoah PDP, Schmidt MK, García-Closas M, Chatterjee N. Common variants in breast cancer risk loci predispose to distinct tumor subtypes. Breast Cancer Res 2022; 24:2. [PMID: 34983606 PMCID: PMC8725568 DOI: 10.1186/s13058-021-01484-x] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Accepted: 11/02/2021] [Indexed: 12/14/2022] Open
Abstract
BACKGROUND Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. METHODS Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. RESULTS Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions. CONCLUSION This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.
Collapse
|
14
|
Jones SJ, Brooks-Wilson A. Anticipation in multiple-case lymphoid cancer families after controlling for ascertainment biases. Leuk Lymphoma 2021; 62:3147-3151. [PMID: 34612777 DOI: 10.1080/10428194.2021.1948026] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
15
|
Martins FC, Couturier DL, Paterson A, Karnezis AN, Chow C, Nazeran TM, Odunsi A, Gentry-Maharaj A, Vrvilo A, Hein A, Talhouk A, Osorio A, Hartkopf AD, Brooks-Wilson A, DeFazio A, Fischer A, Hartmann A, Hernandez BY, McCauley BM, Karpinskyj C, de Sousa CB, Høgdall C, Tiezzi DG, Herpel E, Taran FA, Modugno F, Keeney G, Nelson G, Steed H, Song H, Luk H, Benitez J, Alsop J, Koziak JM, Lester J, Rothstein JH, de Andrade JM, Lundvall L, Paz-Ares L, Robles-Díaz L, Wilkens LR, Garcia MJ, Intermaggio MP, Alcaraz ML, Brett MA, Beckmann MW, Jimenez-Linan M, Anglesio M, Carney ME, Schneider M, Traficante N, Pejovic N, Singh N, Le N, Sinn P, Ghatage P, Erber R, Edwards R, Vierkant R, Ness RB, Leung S, Orsulic S, Brucker SY, Kaufmann SH, Fereday S, Gayther S, Winham SJ, Kommoss S, Pejovic T, Longacre TA, McGuire V, Rhenius V, Sieh W, Shvetsov YB, Whittemore AS, Staebler A, Karlan BY, Rodriguez-Antona C, Bowtell DD, Goode EL, Høgdall E, Candido Dos Reis FJ, Gronwald J, Chang-Claude J, Moysich KB, Kelemen LE, Cook LS, Goodman MT, Fasching PA, Crawford R, Deen S, Menon U, Huntsman DG, Köbel M, Ramus SJ, Pharoah PDP, Brenton JD. Clinical and pathological associations of PTEN expression in ovarian cancer: a multicentre study from the Ovarian Tumour Tissue Analysis Consortium. Br J Cancer 2020; 123:793-802. [PMID: 32555365 PMCID: PMC7463007 DOI: 10.1038/s41416-020-0900-0] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2020] [Revised: 03/14/2020] [Accepted: 04/29/2020] [Indexed: 01/30/2023] Open
Abstract
BACKGROUND PTEN loss is a putative driver in histotypes of ovarian cancer (high-grade serous (HGSOC), endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous (LGSOC)). We aimed to characterise PTEN expression as a biomarker in epithelial ovarian cancer in a large population-based study. METHODS Tumours from 5400 patients from a multicentre observational, prospective cohort study of the Ovarian Tumour Tissue Analysis Consortium were used to evaluate associations between immunohistochemical PTEN patterns and overall survival time, age, stage, grade, residual tumour, CD8+ tumour-infiltrating lymphocytes (TIL) counts, expression of oestrogen receptor (ER), progesterone receptor (PR) and androgen receptor (AR) by means of Cox proportional hazard models and generalised Cochran-Mantel-Haenszel tests. RESULTS Downregulation of cytoplasmic PTEN expression was most frequent in ENOC (most frequently in younger patients; p value = 0.0001) and CCOC and was associated with longer overall survival in HGSOC (hazard ratio: 0.78, 95% CI: 0.65-0.94, p value = 0.022). PTEN expression was associated with ER, PR and AR expression (p values: 0.0008, 0.062 and 0.0002, respectively) in HGSOC and with lower CD8 counts in CCOC (p value < 0.0001). Heterogeneous expression of PTEN was more prevalent in advanced HGSOC (p value = 0.019) and associated with higher CD8 counts (p value = 0.0016). CONCLUSIONS PTEN loss is a frequent driver in ovarian carcinoma associating distinctly with expression of hormonal receptors and CD8+ TIL counts in HGSOC and CCOC histotypes.
Collapse
MESH Headings
- Adenocarcinoma, Clear Cell/enzymology
- Adenocarcinoma, Clear Cell/mortality
- Adenocarcinoma, Clear Cell/pathology
- Age Factors
- Biomarkers, Tumor/biosynthesis
- Biomarkers, Tumor/genetics
- Carcinoma, Ovarian Epithelial/enzymology
- Carcinoma, Ovarian Epithelial/genetics
- Carcinoma, Ovarian Epithelial/mortality
- Carcinoma, Ovarian Epithelial/pathology
- Cohort Studies
- Down-Regulation
- Female
- Gene Knockout Techniques
- Humans
- Middle Aged
- Neoplasm Staging
- Ovarian Neoplasms/enzymology
- Ovarian Neoplasms/genetics
- Ovarian Neoplasms/mortality
- Ovarian Neoplasms/pathology
- PTEN Phosphohydrolase/biosynthesis
- PTEN Phosphohydrolase/deficiency
- PTEN Phosphohydrolase/genetics
- Prospective Studies
- Receptors, Androgen/biosynthesis
- Receptors, Estrogen/biosynthesis
- Receptors, Progesterone/biosynthesis
- Tissue Array Analysis
- Tumor Suppressor Proteins/biosynthesis
- Tumor Suppressor Proteins/deficiency
Collapse
|
16
|
Pemov A, Pathak A, Jones SJ, Dewan R, Merberg J, Karra S, Kim J, Arons E, Ravichandran S, Luke BT, Suman S, Yeager M, Dyer MJS, Lynch HT, Greene MH, Caporaso NE, Kreitman RJ, Goldin LR, Spinelli JJ, Brooks-Wilson A, McMaster ML, Stewart DR. In search of genetic factors predisposing to familial hairy cell leukemia (HCL): exome-sequencing of four multiplex HCL pedigrees. Leukemia 2020; 34:1934-1938. [PMID: 31992839 PMCID: PMC8499084 DOI: 10.1038/s41375-019-0702-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Revised: 12/02/2019] [Accepted: 12/24/2019] [Indexed: 11/09/2022]
|
17
|
Nieuwoudt C, Brooks-Wilson A, Graham J. SimRVSequences: an R package to simulate genetic sequence data for pedigrees. Bioinformatics 2020; 36:2295-2297. [PMID: 31764964 DOI: 10.1093/bioinformatics/btz881] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Revised: 11/12/2019] [Accepted: 11/22/2019] [Indexed: 11/12/2022] Open
Abstract
SUMMARY We present the R package SimRVSequences to simulate sequence data for pedigrees. SimRVSequences allows for simulations of large numbers of single-nucleotide variants (SNVs) and scales well with increasing numbers of pedigrees. Users provide a sample of pedigrees and SNV data from a sample of unrelated individuals. AVAILABILITY AND IMPLEMENTATION SimRVSequences is publicly-available on CRAN https://cran.r-project.org/web/packages/SimRVSequences/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
|
18
|
Gu Q, Sable CM, Brooks-Wilson A, Murphy RA. Dietary patterns in the healthy oldest old in the healthy aging study and the Canadian longitudinal study of aging: a cohort study. BMC Geriatr 2020; 20:106. [PMID: 32178631 PMCID: PMC7077120 DOI: 10.1186/s12877-020-01507-w] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Accepted: 03/06/2020] [Indexed: 11/14/2022] Open
Abstract
Background Very few people live to eighty-five years and older (the ‘oldest old’), and even fewer live to this age without developing chronic diseases. It is important to understand the relationship, if any, of modifiable factors such as diet on healthy aging. However, there are few studies of diet among healthy oldest old, especially in North American populations. We aimed to characterize dietary patterns among ‘super-seniors’ (SS) within the Canadian Healthy Aging Study. Methods 122 SS aged 85 years or older and free of cancer, cardiovascular or pulmonary disease, dementia and diabetes were recruited. Comparisons were made to 12,626 participants aged 65–86 in the Canadian Longitudinal Study of Aging who completed the same 36-item food frequency questionnaire that queried consumption over the prior 12 months of nutrients and foods thought to be important for aging. Dietary patterns were identified with principal component analysis. The odds of being a SS were determined for quartiles of each dietary pattern with logistic regression. Results Two dietary patterns were identified; a western diet characterized by french fries, red meat, processed meat and a nutrient-rich diet which included fruits, vegetables, whole grains, nuts and seeds among other healthy food choices. Higher scores for both dietary patterns were associated with increased odds of being a SS, however, only the western dietary pattern remained associated with adjustment for covariates (Quartile 4: OR = 3.21, 95% CI 1.91–5.51). Conclusions Our finding adds to the limited evidence on dietary intake among the healthiest oldest old but it is unclear whether assocations reflect generational differences between groups or possible contributions to longevity.
Collapse
|
19
|
Fachal L, Aschard H, Beesley J, Barnes DR, Allen J, Kar S, Pooley KA, Dennis J, Michailidou K, Turman C, Soucy P, Lemaçon A, Lush M, Tyrer JP, Ghoussaini M, Moradi Marjaneh M, Jiang X, Agata S, Aittomäki K, Alonso MR, Andrulis IL, Anton-Culver H, Antonenkova NN, Arason A, Arndt V, Aronson KJ, Arun BK, Auber B, Auer PL, Azzollini J, Balmaña J, Barkardottir RB, Barrowdale D, Beeghly-Fadiel A, Benitez J, Bermisheva M, Białkowska K, Blanco AM, Blomqvist C, Blot W, Bogdanova NV, Bojesen SE, Bolla MK, Bonanni B, Borg A, Bosse K, Brauch H, Brenner H, Briceno I, Brock IW, Brooks-Wilson A, Brüning T, Burwinkel B, Buys SS, Cai Q, Caldés T, Caligo MA, Camp NJ, Campbell I, Canzian F, Carroll JS, Carter BD, Castelao JE, Chiquette J, Christiansen H, Chung WK, Claes KBM, Clarke CL, Collée JM, Cornelissen S, Couch FJ, Cox A, Cross SS, Cybulski C, Czene K, Daly MB, de la Hoya M, Devilee P, Diez O, Ding YC, Dite GS, Domchek SM, Dörk T, Dos-Santos-Silva I, Droit A, Dubois S, Dumont M, Duran M, Durcan L, Dwek M, Eccles DM, Engel C, Eriksson M, Evans DG, Fasching PA, Fletcher O, Floris G, Flyger H, Foretova L, Foulkes WD, Friedman E, Fritschi L, Frost D, Gabrielson M, Gago-Dominguez M, Gambino G, Ganz PA, Gapstur SM, Garber J, García-Sáenz JA, Gaudet MM, Georgoulias V, Giles GG, Glendon G, Godwin AK, Goldberg MS, Goldgar DE, González-Neira A, Tibiletti MG, Greene MH, Grip M, Gronwald J, Grundy A, Guénel P, Hahnen E, Haiman CA, Håkansson N, Hall P, Hamann U, Harrington PA, Hartikainen JM, Hartman M, He W, Healey CS, Heemskerk-Gerritsen BAM, Heyworth J, Hillemanns P, Hogervorst FBL, Hollestelle A, Hooning MJ, Hopper JL, Howell A, Huang G, Hulick PJ, Imyanitov EN, Isaacs C, Iwasaki M, Jager A, Jakimovska M, Jakubowska A, James PA, Janavicius R, Jankowitz RC, John EM, Johnson N, Jones ME, Jukkola-Vuorinen A, Jung A, Kaaks R, Kang D, Kapoor PM, Karlan BY, Keeman R, Kerin MJ, Khusnutdinova E, Kiiski JI, Kirk J, Kitahara CM, Ko YD, Konstantopoulou I, Kosma VM, Koutros S, Kubelka-Sabit K, Kwong A, Kyriacou K, Laitman Y, Lambrechts D, Lee E, Leslie G, Lester J, Lesueur F, Lindblom A, Lo WY, Long J, Lophatananon A, Loud JT, Lubiński J, MacInnis RJ, Maishman T, Makalic E, Mannermaa A, Manoochehri M, Manoukian S, Margolin S, Martinez ME, Matsuo K, Maurer T, Mavroudis D, Mayes R, McGuffog L, McLean C, Mebirouk N, Meindl A, Miller A, Miller N, Montagna M, Moreno F, Muir K, Mulligan AM, Muñoz-Garzon VM, Muranen TA, Narod SA, Nassir R, Nathanson KL, Neuhausen SL, Nevanlinna H, Neven P, Nielsen FC, Nikitina-Zake L, Norman A, Offit K, Olah E, Olopade OI, Olsson H, Orr N, Osorio A, Pankratz VS, Papp J, Park SK, Park-Simon TW, Parsons MT, Paul J, Pedersen IS, Peissel B, Peshkin B, Peterlongo P, Peto J, Plaseska-Karanfilska D, Prajzendanc K, Prentice R, Presneau N, Prokofyeva D, Pujana MA, Pylkäs K, Radice P, Ramus SJ, Rantala J, Rau-Murthy R, Rennert G, Risch HA, Robson M, Romero A, Rossing M, Saloustros E, Sánchez-Herrero E, Sandler DP, Santamariña M, Saunders C, Sawyer EJ, Scheuner MT, Schmidt DF, Schmutzler RK, Schneeweiss A, Schoemaker MJ, Schöttker B, Schürmann P, Scott C, Scott RJ, Senter L, Seynaeve CM, Shah M, Sharma P, Shen CY, Shu XO, Singer CF, Slavin TP, Smichkoska S, Southey MC, Spinelli JJ, Spurdle AB, Stone J, Stoppa-Lyonnet D, Sutter C, Swerdlow AJ, Tamimi RM, Tan YY, Tapper WJ, Taylor JA, Teixeira MR, Tengström M, Teo SH, Terry MB, Teulé A, Thomassen M, Thull DL, Tischkowitz M, Toland AE, Tollenaar RAEM, Tomlinson I, Torres D, Torres-Mejía G, Troester MA, Truong T, Tung N, Tzardi M, Ulmer HU, Vachon CM, van Asperen CJ, van der Kolk LE, van Rensburg EJ, Vega A, Viel A, Vijai J, Vogel MJ, Wang Q, Wappenschmidt B, Weinberg CR, Weitzel JN, Wendt C, Wildiers H, Winqvist R, Wolk A, Wu AH, Yannoukakos D, Zhang Y, Zheng W, Hunter D, Pharoah PDP, Chang-Claude J, García-Closas M, Schmidt MK, Milne RL, Kristensen VN, French JD, Edwards SL, Antoniou AC, Chenevix-Trench G, Simard J, Easton DF, Kraft P, Dunning AM. Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes. Nat Genet 2020; 52:56-73. [PMID: 31911677 PMCID: PMC6974400 DOI: 10.1038/s41588-019-0537-1] [Citation(s) in RCA: 97] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Accepted: 10/24/2019] [Indexed: 02/08/2023]
Abstract
Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.
Collapse
|
20
|
Kim S, Wang M, Tyrer JP, Jensen A, Wiensch A, Liu G, Lee AW, Ness RB, Salvatore M, Tworoger SS, Whittemore AS, Anton-Culver H, Sieh W, Olson SH, Berchuck A, Goode EL, Goodman MT, Doherty JA, Chenevix-Trench G, Rossing MA, Webb PM, Giles GG, Terry KL, Ziogas A, Fortner RT, Menon U, Gayther SA, Wu AH, Song H, Brooks-Wilson A, Bandera EV, Cook LS, Cramer DW, Milne RL, Winham SJ, Kjaer SK, Modugno F, Thompson PJ, Chang-Claude J, Harris HR, Schildkraut JM, Le ND, Wentzensen N, Trabert B, Høgdall E, Huntsman D, Pike MC, Pharoah PD, Pearce CL, Mukherjee B. A comprehensive gene-environment interaction analysis in Ovarian Cancer using genome-wide significant common variants. Int J Cancer 2019; 144:2192-2205. [PMID: 30499236 PMCID: PMC6399057 DOI: 10.1002/ijc.32029] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2018] [Accepted: 10/24/2018] [Indexed: 12/18/2022]
Abstract
As a follow-up to genome-wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well-known ovarian cancer risk factors and their interactions with 28 genome-wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case-control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10-4 ). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46-0.60) compared to 0.71 (95%CI = 0.66-0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1-5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene-environment analysis in ovarian cancer.
Collapse
|
21
|
Mavaddat N, Michailidou K, Dennis J, Lush M, Fachal L, Lee A, Tyrer JP, Chen TH, Wang Q, Bolla MK, Yang X, Adank MA, Ahearn T, Aittomäki K, Allen J, Andrulis IL, Anton-Culver H, Antonenkova NN, Arndt V, Aronson KJ, Auer PL, Auvinen P, Barrdahl M, Beane Freeman LE, Beckmann MW, Behrens S, Benitez J, Bermisheva M, Bernstein L, Blomqvist C, Bogdanova NV, Bojesen SE, Bonanni B, Børresen-Dale AL, Brauch H, Bremer M, Brenner H, Brentnall A, Brock IW, Brooks-Wilson A, Brucker SY, Brüning T, Burwinkel B, Campa D, Carter BD, Castelao JE, Chanock SJ, Chlebowski R, Christiansen H, Clarke CL, Collée JM, Cordina-Duverger E, Cornelissen S, Couch FJ, Cox A, Cross SS, Czene K, Daly MB, Devilee P, Dörk T, Dos-Santos-Silva I, Dumont M, Durcan L, Dwek M, Eccles DM, Ekici AB, Eliassen AH, Ellberg C, Engel C, Eriksson M, Evans DG, Fasching PA, Figueroa J, Fletcher O, Flyger H, Försti A, Fritschi L, Gabrielson M, Gago-Dominguez M, Gapstur SM, García-Sáenz JA, Gaudet MM, Georgoulias V, Giles GG, Gilyazova IR, Glendon G, Goldberg MS, Goldgar DE, González-Neira A, Grenaker Alnæs GI, Grip M, Gronwald J, Grundy A, Guénel P, Haeberle L, Hahnen E, Haiman CA, Håkansson N, Hamann U, Hankinson SE, Harkness EF, Hart SN, He W, Hein A, Heyworth J, Hillemanns P, Hollestelle A, Hooning MJ, Hoover RN, Hopper JL, Howell A, Huang G, Humphreys K, Hunter DJ, Jakimovska M, Jakubowska A, Janni W, John EM, Johnson N, Jones ME, Jukkola-Vuorinen A, Jung A, Kaaks R, Kaczmarek K, Kataja V, Keeman R, Kerin MJ, Khusnutdinova E, Kiiski JI, Knight JA, Ko YD, Kosma VM, Koutros S, Kristensen VN, Krüger U, Kühl T, Lambrechts D, Le Marchand L, Lee E, Lejbkowicz F, Lilyquist J, Lindblom A, Lindström S, Lissowska J, Lo WY, Loibl S, Long J, Lubiński J, Lux MP, MacInnis RJ, Maishman T, Makalic E, Maleva Kostovska I, Mannermaa A, Manoukian S, Margolin S, Martens JWM, Martinez ME, Mavroudis D, McLean C, Meindl A, Menon U, Middha P, Miller N, Moreno F, Mulligan AM, Mulot C, Muñoz-Garzon VM, Neuhausen SL, Nevanlinna H, Neven P, Newman WG, Nielsen SF, Nordestgaard BG, Norman A, Offit K, Olson JE, Olsson H, Orr N, Pankratz VS, Park-Simon TW, Perez JIA, Pérez-Barrios C, Peterlongo P, Peto J, Pinchev M, Plaseska-Karanfilska D, Polley EC, Prentice R, Presneau N, Prokofyeva D, Purrington K, Pylkäs K, Rack B, Radice P, Rau-Murthy R, Rennert G, Rennert HS, Rhenius V, Robson M, Romero A, Ruddy KJ, Ruebner M, Saloustros E, Sandler DP, Sawyer EJ, Schmidt DF, Schmutzler RK, Schneeweiss A, Schoemaker MJ, Schumacher F, Schürmann P, Schwentner L, Scott C, Scott RJ, Seynaeve C, Shah M, Sherman ME, Shrubsole MJ, Shu XO, Slager S, Smeets A, Sohn C, Soucy P, Southey MC, Spinelli JJ, Stegmaier C, Stone J, Swerdlow AJ, Tamimi RM, Tapper WJ, Taylor JA, Terry MB, Thöne K, Tollenaar RAEM, Tomlinson I, Truong T, Tzardi M, Ulmer HU, Untch M, Vachon CM, van Veen EM, Vijai J, Weinberg CR, Wendt C, Whittemore AS, Wildiers H, Willett W, Winqvist R, Wolk A, Yang XR, Yannoukakos D, Zhang Y, Zheng W, Ziogas A, Dunning AM, Thompson DJ, Chenevix-Trench G, Chang-Claude J, Schmidt MK, Hall P, Milne RL, Pharoah PDP, Antoniou AC, Chatterjee N, Kraft P, García-Closas M, Simard J, Easton DF. Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes. Am J Hum Genet 2019; 104:21-34. [PMID: 30554720 PMCID: PMC6323553 DOI: 10.1016/j.ajhg.2018.11.002] [Citation(s) in RCA: 560] [Impact Index Per Article: 112.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2018] [Accepted: 11/03/2018] [Indexed: 12/29/2022] Open
Abstract
Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.
Collapse
|
22
|
Kelemen LE, Earp M, Fridley BL, Chenevix-Trench G, Fasching PA, Beckmann MW, Ekici AB, Hein A, Lambrechts D, Lambrechts S, Van Nieuwenhuysen E, Vergote I, Rossing MA, Doherty JA, Chang-Claude J, Behrens S, Moysich KB, Cannioto R, Lele S, Odunsi K, Goodman MT, Shvetsov YB, Thompson PJ, Wilkens LR, Dörk T, Antonenkova N, Bogdanova N, Hillemanns P, Runnebaum IB, du Bois A, Harter P, Heitz F, Schwaab I, Butzow R, Pelttari LM, Nevanlinna H, Modugno F, Edwards RP, Kelley JL, Ness RB, Karlan BY, Lester J, Orsulic S, Walsh C, Kjaer SK, Jensen A, Cunningham JM, Vierkant RA, Giles GG, Bruinsma F, Southey MC, Hildebrandt MA, Liang D, Lu K, Wu X, Sellers TA, Levine DA, Schildkraut JM, Iversen ES, Terry KL, Cramer DW, Tworoger SS, Poole EM, Bandera EV, Olson SH, Orlow I, Vestrheim Thomsen LC, Bjorge L, Krakstad C, Tangen IL, Kiemeney LA, Aben KK, Massuger LF, van Altena AM, Pejovic T, Bean Y, Kellar M, Cook LS, Le ND, Brooks-Wilson A, Gronwald J, Cybulski C, Jakubowska A, Lubiński J, Wentzensen N, Brinton LA, Lissowska J, Hogdall E, Engelholm SA, Hogdall C, Lundvall L, Nedergaard L, Pharoah PD, Dicks E, Song H, Tyrer JP, McNeish I, Siddiqui N, Carty K, Glasspool R, Paul J, Campbell IG, Eccles D, Whittemore AS, McGuire V, Rothstein JH, Sieh W, Narod SA, Phelan CM, McLaughlin JR, Risch HA, Anton-Culver H, Ziogas A, Menon U, Gayther SA, Gentry-Maharaj A, Ramus SJ, Wu AH, Pearce CL, Lee AW, Pike MC, Kupryjanczyk J, Podgorska A, Plisiecka-Halasa J, Sawicki W, Goode EL, Berchuck A. rs495139 in the TYMS-ENOSF1 Region and Risk of Ovarian Carcinoma of Mucinous Histology. Int J Mol Sci 2018; 19:E2473. [PMID: 30134598 PMCID: PMC6163881 DOI: 10.3390/ijms19092473] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2018] [Revised: 08/05/2018] [Accepted: 08/09/2018] [Indexed: 12/20/2022] Open
Abstract
Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3' gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97⁻1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03⁻1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10-28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.
Collapse
|
23
|
Wang SS, Carrington M, Berndt SI, Slager SL, Bracci PM, Voutsinas J, Cerhan JR, Smedby KE, Hjalgrim H, Vijai J, Morton LM, Vermeulen R, Paltiel O, Vajdic CM, Linet MS, Nieters A, de Sanjose S, Cozen W, Brown EE, Turner J, Spinelli JJ, Zheng T, Birmann BM, Flowers CR, Becker N, Holly EA, Kane E, Weisenburger D, Maynadie M, Cocco P, Albanes D, Weinstein SJ, Teras LR, Diver WR, Lax SJ, Travis RC, Kaaks R, Riboli E, Benavente Y, Brennan P, McKay J, Delfau-Larue MH, Link BK, Magnani C, Ennas MG, Latte G, Feldman AL, Doo NW, Giles GG, Southey MC, Milne RL, Offit K, Musinsky J, Arslan AA, Purdue MP, Adami HO, Melbye M, Glimelius B, Conde L, Camp NJ, Glenn M, Curtin K, Clavel J, Monnereau A, Cox DG, Ghesquières H, Salles G, Bofetta P, Foretova L, Staines A, Davis S, Severson RK, Lan Q, Brooks-Wilson A, Smith MT, Roman E, Kricker A, Zhang Y, Kraft P, Chanock SJ, Rothman N, Hartge P, Skibola CF. HLA Class I and II Diversity Contributes to the Etiologic Heterogeneity of Non-Hodgkin Lymphoma Subtypes. Cancer Res 2018; 78:4086-4096. [PMID: 29735552 PMCID: PMC6065509 DOI: 10.1158/0008-5472.can-17-2900] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2017] [Revised: 02/07/2018] [Accepted: 04/24/2018] [Indexed: 01/07/2023]
Abstract
A growing number of loci within the human leukocyte antigen (HLA) region have been implicated in non-Hodgkin lymphoma (NHL) etiology. Here, we test a complementary hypothesis of "heterozygote advantage" regarding the role of HLA and NHL, whereby HLA diversity is beneficial and homozygous HLA loci are associated with increased disease risk. HLA alleles at class I and II loci were imputed from genome-wide association studies (GWAS) using SNP2HLA for 3,617 diffuse large B-cell lymphomas (DLBCL), 2,686 follicular lymphomas (FL), 2,878 chronic lymphocytic leukemia/small lymphocytic lymphomas (CLL/SLL), 741 marginal zone lymphomas (MZL), and 8,753 controls of European descent. Both DLBCL and MZL risk were elevated with homozygosity at class I HLA-B and -C loci (OR DLBCL = 1.31, 95% CI = 1.06-1.60; OR MZL = 1.45, 95% CI = 1.12-1.89) and class II HLA-DRB1 locus (OR DLBCL = 2.10, 95% CI = 1.24-3.55; OR MZL = 2.10, 95% CI = 0.99-4.45). Increased FL risk was observed with the overall increase in number of homozygous HLA class II loci (P trend < 0.0001, FDR = 0.0005). These results support a role for HLA zygosity in NHL etiology and suggests that distinct immune pathways may underly the etiology of the different NHL subtypes.Significance: HLA gene diversity reduces risk for non-Hodgkin lymphoma. Cancer Res; 78(14); 4086-96. ©2018 AACR.
Collapse
|
24
|
Earp M, Tyrer JP, Winham SJ, Lin HY, Chornokur G, Dennis J, Aben KKH, Anton‐Culver H, Antonenkova N, Bandera EV, Bean YT, Beckmann MW, Bjorge L, Bogdanova N, Brinton LA, Brooks-Wilson A, Bruinsma F, Bunker CH, Butzow R, Campbell IG, Carty K, Chang-Claude J, Cook LS, Cramer DW, Cunningham JM, Cybulski C, Dansonka-Mieszkowska A, Despierre E, Doherty JA, Dörk T, du Bois A, Dürst M, Easton DF, Eccles DM, Edwards RP, Ekici AB, Fasching PA, Fridley BL, Gentry-Maharaj A, Giles GG, Glasspool R, Goodman MT, Gronwald J, Harter P, Hein A, Heitz F, Hildebrandt MAT, Hillemanns P, Hogdall CK, Høgdall E, Hosono S, Iversen ES, Jakubowska A, Jensen A, Ji BT, Jung AY, Karlan BY, Kellar M, Kiemeney LA, Kiong Lim B, Kjaer SK, Krakstad C, Kupryjanczyk J, Lambrechts D, Lambrechts S, Le ND, Lele S, Lester J, Levine DA, Li Z, Liang D, Lissowska J, Lu K, Lubinski J, Lundvall L, Massuger LFAG, Matsuo K, McGuire V, McLaughlin JR, McNeish I, Menon U, Milne RL, Modugno F, Moysich KB, Ness RB, Nevanlinna H, Odunsi K, Olson SH, Orlow I, Orsulic S, Paul J, Pejovic T, Pelttari LM, Permuth JB, Pike MC, Poole EM, Rosen B, Rossing MA, Rothstein JH, Runnebaum IB, Rzepecka IK, Schernhammer E, Schwaab I, Shu XO, Shvetsov YB, Siddiqui N, Sieh W, Song H, Southey MC, Spiewankiewicz B, Sucheston-Campbell L, Tangen IL, Teo SH, Terry KL, Thompson PJ, Thomsen L, Tworoger SS, van Altena AM, Vergote I, Vestrheim Thomsen LC, Vierkant RA, Walsh CS, Wang-Gohrke S, Wentzensen N, Whittemore AS, Wicklund KG, Wilkens LR, Woo YL, Wu AH, Wu X, Xiang YB, Yang H, Zheng W, Ziogas A, Lee AW, Pearce CL, Berchuck A, Schildkraut JM, Ramus SJ, Monteiro ANA, Narod SA, Sellers TA, Gayther SA, Kelemen LE, Chenevix-Trench G, Risch HA, Pharoah PDP, Goode EL, Phelan CM. Variants in genes encoding small GTPases and association with epithelial ovarian cancer susceptibility. PLoS One 2018; 13:e0197561. [PMID: 29979793 PMCID: PMC6034790 DOI: 10.1371/journal.pone.0197561] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2017] [Accepted: 05/06/2018] [Indexed: 11/29/2022] Open
Abstract
Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.
Collapse
|
25
|
Jones S, Uvarov O, Salema D, Connors J, Brooks-Wilson A. PO-064 Early birth order and increased risk of lymphoid cancers and allergies in lymphoid cancer families. ESMO Open 2018. [DOI: 10.1136/esmoopen-2018-eacr25.594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022] Open
|