Genomic data in the All of Us Research Program

Comprehensively mapping the genetic basis of human disease across diverse individuals is a long-standing goal for the field of human genetics1-4. The All of Us Research Program is a longitudinal cohort study aiming to enrol a diverse group of at least one million individuals across the USA to accelerate biomedical research and improve human health5,6. Here we describe the programme's genomics data release of 245,388 clinical-grade genome sequences. This resource is unique in its diversity as 77% of participants are from communities that are historically under-represented in biomedical research and 46% are individuals from under-represented racial and ethnic minorities. All of Us identified more than 1 billion genetic variants, including more than 275 million previously unreported genetic variants, more than 3.9 million of which had coding consequences. Leveraging linkage between genomic data and the longitudinal electronic health record, we evaluated 3,724 genetic variants associated with 117 diseases and found high replication rates across both participants of European ancestry and participants of African ancestry. Summary-level data are publicly available, and individual-level data can be accessed by researchers through the All of Us Researcher Workbench using a unique data passport model with a median time from initial researcher registration to data access of 29 hours. We anticipate that this diverse dataset will advance the promise of genomic medicine for all.

Researching COVID to Enhance Recovery (RECOVER) adult study protocol: Rationale, objectives, and design

IMPORTANCE

SARS-CoV-2 infection can result in ongoing, relapsing, or new symptoms or other health effects after the acute phase of infection; termed post-acute sequelae of SARS-CoV-2 infection (PASC), or long COVID. The characteristics, prevalence, trajectory and mechanisms of PASC are ill-defined. The objectives of the Researching COVID to Enhance Recovery (RECOVER) Multi-site Observational Study of PASC in Adults (RECOVER-Adult) are to: (1) characterize PASC prevalence; (2) characterize the symptoms, organ dysfunction, natural history, and distinct phenotypes of PASC; (3) identify demographic, social and clinical risk factors for PASC onset and recovery; and (4) define the biological mechanisms underlying PASC pathogenesis.

METHODS

RECOVER-Adult is a combined prospective/retrospective cohort currently planned to enroll 14,880 adults aged ≥18 years. Eligible participants either must meet WHO criteria for suspected, probable, or confirmed infection; or must have evidence of no prior infection. Recruitment occurs at 86 sites in 33 U.S. states, Washington, DC and Puerto Rico, via facility- and community-based outreach. Participants complete quarterly questionnaires about symptoms, social determinants, vaccination status, and interim SARS-CoV-2 infections. In addition, participants contribute biospecimens and undergo physical and laboratory examinations at approximately 0, 90 and 180 days from infection or negative test date, and yearly thereafter. Some participants undergo additional testing based on specific criteria or random sampling. Patient representatives provide input on all study processes. The primary study outcome is onset of PASC, measured by signs and symptoms. A paradigm for identifying PASC cases will be defined and updated using supervised and unsupervised learning approaches with cross-validation. Logistic regression and proportional hazards regression will be conducted to investigate associations between risk factors, onset, and resolution of PASC symptoms.

DISCUSSION

RECOVER-Adult is the first national, prospective, longitudinal cohort of PASC among US adults. Results of this study are intended to inform public health, spur clinical trials, and expand treatment options.

REGISTRATION

NCT05172024.

Pragmatic evaluation of events and benefits of lipid lowering in older adults (PREVENTABLE): Trial design and rationale

Whether initiation of statins could increase survival free of dementia and disability in adults aged ≥75 years is unknown. PREVENTABLE, a double-blind, placebo-controlled randomized pragmatic clinical trial, will compare high-intensity statin therapy (atorvastatin 40 mg) with placebo in 20,000 community-dwelling adults aged ≥75 years without cardiovascular disease, disability, or dementia at baseline. Exclusion criteria include statin use in the prior year or for >5 years and inability to take a statin. Potential participants are identified using computable phenotypes derived from the electronic health record and local referrals from the community. Participants will undergo baseline cognitive testing, with physical testing and a blinded lipid panel if feasible. Cognitive testing and disability screening will be conducted annually. Multiple data sources will be queried for cardiovascular events, dementia, and disability; survival is site-reported and supplemented by a National Death Index search. The primary outcome is survival free of new dementia or persisting disability. Co-secondary outcomes are a composite of cardiovascular death, hospitalization for unstable angina or myocardial infarction, heart failure, stroke, or coronary revascularization; and a composite of mild cognitive impairment or dementia. Ancillary studies will offer mechanistic insights into the effects of statins on key outcomes. Biorepository samples are obtained and stored for future study. These results will inform the benefit of statins for increasing survival free of dementia and disability among older adults. This is a pioneering pragmatic study testing important questions with low participant burden to align with the needs of the growing population of older adults.

Researching COVID to enhance recovery (RECOVER) pediatric study protocol: Rationale, objectives and design

Importance

The prevalence, pathophysiology, and long-term outcomes of COVID-19 (post-acute sequelae of SARS-CoV-2 [PASC] or "Long COVID") in children and young adults remain unknown. Studies must address the urgent need to define PASC, its mechanisms, and potential treatment targets in children and young adults.

Observations

We describe the protocol for the Pediatric Observational Cohort Study of the NIH's RE searching COV ID to E nhance R ecovery (RECOVER) Initiative. RECOVER-Pediatrics is an observational meta-cohort study of caregiver-child pairs (birth through 17 years) and young adults (18 through 25 years), recruited from more than 100 sites across the US. This report focuses on two of five cohorts that comprise RECOVER-Pediatrics: 1) a de novo RECOVER prospective cohort of children and young adults with and without previous or current infection; and 2) an extant cohort derived from the Adolescent Brain Cognitive Development (ABCD) study ( n =10,000). The de novo cohort incorporates three tiers of data collection: 1) remote baseline assessments (Tier 1, n=6000); 2) longitudinal follow-up for up to 4 years (Tier 2, n=6000); and 3) a subset of participants, primarily the most severely affected by PASC, who will undergo deep phenotyping to explore PASC pathophysiology (Tier 3, n=600). Youth enrolled in the ABCD study participate in Tier 1. The pediatric protocol was developed as a collaborative partnership of investigators, patients, researchers, clinicians, community partners, and federal partners, intentionally promoting inclusivity and diversity. The protocol is adaptive to facilitate responses to emerging science.

Conclusions and Relevance

RECOVER-Pediatrics seeks to characterize the clinical course, underlying mechanisms, and long-term effects of PASC from birth through 25 years old. RECOVER-Pediatrics is designed to elucidate the epidemiology, four-year clinical course, and sociodemographic correlates of pediatric PASC. The data and biosamples will allow examination of mechanistic hypotheses and biomarkers, thus providing insights into potential therapeutic interventions.

Clinical Trialsgov Identifier

Clinical Trial Registration: http://www.clinicaltrials.gov . Unique identifier: NCT05172011.

Antibodies to SARS-CoV-2 in All of Us Research Program Participants, January 2-March 18, 2020

Background

With limited severe acute respiratory syndrome coronavirus (SARS-CoV-2) testing capacity in the United States at the start of the epidemic (January–March 2020), testing was focused on symptomatic patients with a travel history throughout February, obscuring the picture of SARS-CoV-2 seeding and community transmission. We sought to identify individuals with SARS-CoV-2 antibodies in the early weeks of the US epidemic.

Methods

All of Us study participants in all 50 US states provided blood specimens during study visits from 2 January to 18 March 2020. Participants were considered seropositive if they tested positive for SARS-CoV-2 immunoglobulin G (IgG) antibodies with the Abbott Architect SARS-CoV-2 IgG enzyme-linked immunosorbent assay (ELISA) and the EUROIMMUN SARS-CoV-2 ELISA in a sequential testing algorithm. The sensitivity and specificity of these ELISAs and the net sensitivity and specificity of the sequential testing algorithm were estimated, along with 95% confidence intervals (CIs).

Results

The estimated sensitivities of the Abbott and EUROIMMUN assays were 100% (107 of 107 [95% CI: 96.6%–100%]) and 90.7% (97 of 107 [83.5%–95.4%]), respectively, and the estimated specificities were 99.5% (995 of 1000 [98.8%–99.8%]) and 99.7% (997 of 1000 [99.1%–99.9%]), respectively. The net sensitivity and specificity of our sequential testing algorithm were 90.7% (97 of 107 [95% CI: 83.5%–95.4%]) and 100.0% (1000 of 1000 [99.6%–100%]), respectively. Of the 24 079 study participants with blood specimens from 2 January to 18 March 2020, 9 were seropositive, 7 before the first confirmed case in the states of Illinois, Massachusetts, Wisconsin, Pennsylvania, and Mississippi.

Conclusions

Our findings identified SARS-CoV-2 infections weeks before the first recognized cases in 5 US states.

Mini-Review of Laboratory Operations in Biobanking: Building Biobanking Resources for Translational Research

Biobanks have become integral to improving population health. We are in a new era in medicine as patients, health professionals, and researchers increasingly collaborate to gain new knowledge and explore new paradigms for diagnosing and treating disease. Many large-scale biobanking efforts are underway worldwide at the institutional, national, and even international level. When linked with subject data from questionnaires and medical records, biobanks serve as valuable resources in translational research. A biobank must have high quality samples that meet researcher's needs. Biobank laboratory operations require an enormous amount of support—from lab and storage space, information technology expertise, and a laboratory management information system to logistics for sample movement, quality management systems, and appropriate facilities. A paramount metric of success for a biobank is the concept of every biospecimen coming to the repository belongs to a participant who has something to contribute to research for a healthier future. This article will discuss the importance of biorepository operations, specific to the collection and storage of participants materials. Specific focus will be given to maintaining the quality of samples, along with the various levels of support biorepositories need to fulfill their purpose and ensure the integrity of each specimen is maintained.

Impact of Diverse Data Sources on Computational Phenotyping

Electronic health records (EHRs) are widely adopted with a great potential to serve as a rich, integrated source of phenotype information. Computational phenotyping, which extracts phenotypes from EHR data automatically, can accelerate the adoption and utilization of phenotype-driven efforts to advance scientific discovery and improve healthcare delivery. A list of computational phenotyping algorithms has been published but data fragmentation, i.e., incomplete data within one single data source, has been raised as an inherent limitation of computational phenotyping. In this study, we investigated the impact of diverse data sources on two published computational phenotyping algorithms, rheumatoid arthritis (RA) and type 2 diabetes mellitus (T2DM), using Mayo EHRs and Rochester Epidemiology Project (REP) which links medical records from multiple health care systems. Results showed that both RA (less prevalent) and T2DM (more prevalent) case selections were markedly impacted by data fragmentation, with positive predictive value (PPV) of 91.4 and 92.4%, false-negative rate (FNR) of 26.6 and 14% in Mayo data, respectively, PPV of 97.2 and 98.3%, FNR of 5.2 and 3.3% in REP. T2DM controls also contain biases, with PPV of 91.2% and FNR of 1.2% for Mayo. We further elaborated underlying reasons impacting the performance.

The impact of sample processing on inflammatory markers in serum: Lessons learned

Objectives: To investigate the effect of sample handling on inflammatory cytokines in serum and highlight challenges with using samples pre-collected from biobanks for biomarker research.Methods: Cytokine concentrations (IL-1β, IL-2, IL-6, IL-8, IL-10, TNFα, and IFNγ) were measured in serum samples of 205 patients with bipoldar disorder (BD) from the Mayo Clinic Bipolar Disorder Biobank and 205 non-psychiatric controls from the Mayo Clinic Biobank. As cytokine concentrations varied by recruitment site, post-hoc models were used to test the effect of clinical variables and pre-processing time on cytokines. To evaluate the effect of pre-processing time experimentally, cytokines were assayed in serum and plasma from 6 healthy volunteers processed at different time points.Results: Cytokine levels were significantly higher in the BD group. However, both cytokine levels and pre-processing times differed by recruitment site, and post-hoc analyses revealed that pre-processing time was significantly associated with several cytokines. An experiment using samples from healthy volunteers confirmed that concentrations for most cytokines increased with longer pre-processing times.Conclusions: Delays in processing influence cytokine concentrations in blood samples. Given the increasing use of biobanks in research, this study highlights the need to carefully evaluate sample collection and handling methods when designing biomarker studies.

Characteristics Associated With Recruitment and Re-contact in Mayo Clinic Biobank

Objective: To better understand the characteristics associated with a participant's willingness to consent to the Mayo Clinic Biobank (MCB) and examine factors associated with willingness to participate in follow-up studies embedded within MCB that require re-contact and participant approval. Participants and Methods: Consent rates were compared across patient demographics to the MCB. Rates of participation to follow-up studies were also compared across demographics and request types. Results: Among 272,102 Mayo Clinic patients invited to the MCB, 48,314 (19%) consented across the three recruitment sites within 90 days of initial invitation. A significant age by gender interaction was identified, showing young males consent at a lower rate than young females and older males consent at a higher rate than older females. Over the recruitment time frame of 2009-2015, there was a significant decrease in consent rates (decline of 2.5%/year). Of the 57,041 consented MCB participants, 33,487 participants (59%) have been invited to participate in follow-up studies via re-contact. Follow-up studies of the MCB may require participants to provide additional samples, complete questionnaires, and/or release their identity to a research team. MCB participants have been invited to enroll in a median of two studies (IQR: 1-3). Seventy-one percent of participants consented to at least one follow-up study, with individual follow-up study consent rates ranging from 14 to 87% depending on study type, with a median consent rate of 61% (IQR: 47-70%). Studies requesting return of a questionnaire had the highest participation rates. White participants, older participants, and participants with some college or a degree were significantly more likely to participate to follow-up studies, while there was no association with gender. Conclusion: Consent rates among younger and non-white patients were lower than in older, white patients. However, we also found that participation rates among those already enrolled in the biobank were much higher than those seen in new recruitment efforts, external to an existing biobank. We thus demonstrate an important way that biobanks can advance precision medicine goals: through provision of populations from which studies can draw participants for future studies.

Association of mitochondrial DNA copy number with self-rated health status

Purpose

In aging adults, mitochondrial dysfunction may be an important contributor. We evaluated the association between mitochondrial DNA (mtDNA) copy number, which is a biomarker for mitochondrial function, and self-rated health status.

Patients and methods

We conducted a cross-sectional study of patients enrolled within the Mayo Clinic Biobank. We utilized the questionnaire and sequence data from 944 patients. We examined the association between mtDNA copy number and self-rated health status with 3 collapsed categories for the latter variable (excellent/very good, good, and fair/poor). For analysis, we used proportional odds models after log-transforming mtDNA copy number, and we adjusted for age and sex.

Results

We found the median age at enrollment was 61 years (25th–75th percentile: 51–71), and 64% reported excellent or very good health, 31% reported good health, and 6% reported fair/poor health. Overall, the median mtDNA copy number was 88.9 (25th–75th percentile: 77.6–101.1). Higher mtDNA copy number was found for subjects reporting better self-rated health status after adjusting for age, sex, and comorbidity burden (OR =2.3 [95% CI: 1.2–4.5] for having better self-rated health for a one-unit increase in log-transformed mtDNA copy number).

Conclusion

We found that a higher mtDNA copy number is associated with better self-rated health status after adjustment for age, sex, and comorbidity burden. The current study implies that mtDNA copy number may serve as a biomarker for self-reported health. Further studies, potentially including cohort studies, may be required.

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.

PPM1D Mosaic Truncating Variants in Ovarian Cancer Cases May Be Treatment-Related Somatic Mutations

Mosaic truncating mutations in the protein phosphatase, Mg(2+)/Mn(2+)-dependent, 1D (PPM1D) gene have recently been reported with a statistically significantly greater frequency in lymphocyte DNA from ovarian cancer case patients compared with unaffected control patients. Using massively parallel sequencing (MPS) we identified truncating PPM1D mutations in 12 of 3236 epithelial ovarian cancer (EOC) case patients (0.37%) but in only one of 3431 unaffected control patients (0.03%) (P = .001). All statistical tests were two-sided. A combination of Sanger sequencing, pyrosequencing, and MPS data suggested that 12 of the 13 mutations were mosaic. All mutations were identified in post-chemotherapy treatment blood samples from case patients (n = 1827) (average 1234 days post-treatment in carriers) rather than from cases collected pretreatment (less than 14 days after diagnosis, n = 1384) (P = .002). These data suggest that PPM1D variants in EOC cases are primarily somatic mosaic mutations caused by treatment and are not associated with germline predisposition to EOC.

Germline Mutations in the BRIP1, BARD1, PALB2, and NBN Genes in Women With Ovarian Cancer

BACKGROUND

Epithelial ovarian cancer (EOC) is the most lethal gynecological malignancy, responsible for 13 000 deaths per year in the United States. Risk prediction based on identifying germline mutations in ovarian cancer susceptibility genes could have a clinically significant impact on reducing disease mortality.

METHODS

Next generation sequencing was used to identify germline mutations in the coding regions of four candidate susceptibility genes-BRIP1, BARD1, PALB2 and NBN-in 3236 invasive EOC case patients and 3431 control patients of European origin, and in 2000 unaffected high-risk women from a clinical screening trial of ovarian cancer (UKFOCSS). For each gene, we estimated the prevalence and EOC risks and evaluated associations between germline variant status and clinical and epidemiological risk factor information. All statistical tests were two-sided.

RESULTS

We found an increased frequency of deleterious mutations in BRIP1 in case patients (0.9%) and in the UKFOCSS participants (0.6%) compared with control patients (0.09%) (P = 1 x 10(-4) and 8 x 10(-4), respectively), but no differences for BARD1 (P = .39), NBN1 ( P = .61), or PALB2 (P = .08). There was also a difference in the frequency of rare missense variants in BRIP1 between case patients and control patients (P = 5.5 x 10(-4)). The relative risks associated with BRIP1 mutations were 11.22 for invasive EOC (95% confidence interval [CI] = 3.22 to 34.10, P = 1 x 10(-4)) and 14.09 for high-grade serous disease (95% CI = 4.04 to 45.02, P = 2 x 10(-5)). Segregation analysis in families estimated the average relative risks in BRIP1 mutation carriers compared with the general population to be 3.41 (95% CI = 2.12 to 5.54, P = 7×10(-7)).

CONCLUSIONS

Deleterious germline mutations in BRIP1 are associated with a moderate increase in EOC risk. These data have clinical implications for risk prediction and prevention approaches for ovarian cancer and emphasize the critical need for risk estimates based on very large sample sizes before genes of moderate penetrance have clinical utility in cancer prevention.

Contribution of Germline Mutations in the RAD51B, RAD51C, and RAD51D Genes to Ovarian Cancer in the Population

PURPOSE

The aim of this study was to estimate the contribution of deleterious mutations in the RAD51B, RAD51C, and RAD51D genes to invasive epithelial ovarian cancer (EOC) in the population and in a screening trial of individuals at high risk of ovarian cancer.

PATIENTS AND METHODS

The coding sequence and splice site boundaries of the three RAD51 genes were sequenced and analyzed in germline DNA from a case-control study of 3,429 patients with invasive EOC and 2,772 controls as well as in 2,000 unaffected women who were BRCA1/BRCA2 negative from the United Kingdom Familial Ovarian Cancer Screening Study (UK_FOCSS) after quality-control analysis.

RESULTS

In the case-control study, we identified predicted deleterious mutations in 28 EOC cases (0.82%) compared with three controls (0.11%; P < .001). Mutations in EOC cases were more frequent in RAD51C (14 occurrences, 0.41%) and RAD51D (12 occurrences, 0.35%) than in RAD51B (two occurrences, 0.06%). RAD51C mutations were associated with an odds ratio of 5.2 (95% CI, 1.1 to 24; P = .035), and RAD51D mutations conferred an odds ratio of 12 (95% CI, 1.5 to 90; P = .019). We identified 13 RAD51 mutations (0.65%) in unaffected UK_FOCSS participants (RAD51C, n = 7; RAD51D, n = 5; and RAD51B, n = 1), which was a significantly greater rate than in controls (P < .001); furthermore, RAD51 mutation carriers were more likely than noncarriers to have a family history of ovarian cancer (P < .001).

CONCLUSION

These results confirm that RAD51C and RAD51D are moderate ovarian cancer susceptibility genes and suggest that they confer levels of risk of EOC that may warrant their use alongside BRCA1 and BRCA2 in routine clinical genetic testing.

Abstract 260: Integrative genomic analysis identifies epigenetic marks that mediate genetic risk for epithelial ovarian cancer

Both genetic and epigenetic factors influence the development and progression of epithelial ovarian cancer (EOC). However, there is an incomplete understanding of the interrelationship between these factors and whether they interact to impact disease risk. Given that the development of ovarian cancer is complex, we aimed to identify DNA methylation marks that are candidates for mediating ovarian cancer genetic risk. We used 214 cases and 214 age-matched controls from the Mayo Clinic Ovarian Cancer Study. Pretreatment, blood-derived DNA was profiled for genome-wide methylation (Illumina Infinium HumanMethylation27 BeadArray) and single nucleotide polymorphisms (SNPs, Illumina Infinium HD Human610-Quad BeadArray). We employed a three-step filtering procedure, followed by the Causal Inference Test (CIT), to distinguish CpG sites that mediate genetic risk, from those that are consequential or independently acted on by genotype. Controlling for the estimated distribution of immune cells and other key covariates, we identified 1,993 out of 25,926 (7.7%) CpGs that were significantly differentially methylated between cases and controls (FDR, q &lt; 0.05). The relationship between methylation and case-control status among these 1,993 CpGs was found to be highly consistent with the results of an independent study. Implementation of the CIT test revealed 13 CpG/SNP pairs, comprising 13 unique CpGs and 17 unique SNPs, which represent potential methylation-mediated relationships between genotype and EOC risk. These findings provide additional insight into EOC etiology and may serve as novel biomarkers for EOC susceptibility.

Citation Format: Devin C. Koestler, Prabhakar Chalise, Mine S. Cicek, Julie M. Cunningham, Sebastian Armasu, Melissa C. Larson, Jeremy Chien, Matthew Block, Kimberly R. Kalli, Thomas A. Sellers, Ellen L. Goode, Brooke L. Fridley. Integrative genomic analysis identifies epigenetic marks that mediate genetic risk for epithelial ovarian cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 260. doi:10.1158/1538-7445.AM2014-260

The contribution of deleterious germline mutations in BRCA1, BRCA2 and the mismatch repair genes to ovarian cancer in the population

The aim of this study was to estimate the contribution of deleterious mutations in BRCA1, BRCA2, MLH1, MSH2, MSH6 and PMS2 to invasive epithelial ovarian cancer (EOC) in the population. The coding sequence and splice site boundaries of all six genes were amplified in germline DNA from 2240 invasive EOC cases and 1535 controls. Barcoded fragment libraries were sequenced using the Illumina GAII or HiSeq and sequence data for each subject de-multiplexed prior to interpretation. GATK and Annovar were used for variant detection and annotation. After quality control 2222 cases (99.2%) and 1528 controls (99.5%) were included in the final analysis. We identified 193 EOC cases (8.7%) carrying a deleterious mutation in at least one gene compared with 10 controls (0.65%). Mutations were most frequent in BRCA1 and BRCA2, with 84 EOC cases (3.8%) carrying a BRCA1 mutation and 94 EOC cases (4.2%) carrying a BRCA2 mutation. The combined BRCA1 and BRCA2 mutation prevalence was 11% in high-grade serous disease. Seventeen EOC cases carried a mutation in a mismatch repair gene, including 10 MSH6 mutation carriers (0.45%) and 4 MSH2 mutation carriers (0.18%). At least 1 in 10 women with high-grade serous EOC has a BRCA1 or BRCA2 mutation. The development of next generation sequencing technologies enables rapid mutation screening for multiple susceptibility genes at once, suggesting that routine clinical testing of all incidence cases should be considered.

Genome-wide investigation of regional blood-based DNA methylation adjusted for complete blood counts implicates BNC2 in ovarian cancer

Due to its potential as a biomarker for early cancer detection, blood-based DNA methylation (DNAm) is of interest in cancer research. Specifically, highly predictive mechanisms for early detection of epithelial ovarian cancer (EOC) are desired, so previous studies have compared DNAm between EOC cases and controls. However, case-control studies are confounded by the distribution of white blood cell types through an immune response induced by the cancer. Rather than determining the distribution of the cell types manually or investigating isolated cell types, an alternative approach involves the use of complete blood count (CBC), which is routinely collected. In the analysis of an EOC case-control study of DNAm, we incorporate CBC measures to adjust for this confounding and compare DNAm between 242 EOC cases and 181 age-matched controls (assayed on the Illumina Infinium HumanMethylation27 or HumanMethylation450 Beadchips), at both the individual CpG and CpG island levels. We found that adjustment for leukocyte distribution using CBC measurements dramatically reduced confounding, with 62 single CpG sites found to be associated with EOC status after adjustment (P < 5E-8). Additionally, regional DNAm was assessed by applying principal components analysis to CpG islands. The top associated CpG island (P = 7E-6) was located in the promoter/transcription start site of the human basonuclin 2 gene (BNC2), a known susceptibility gene for EOC risk identified through GWAS. Follow-up studies are necessary to establish the role of BNC2 in blood-based DNA and EOC, including prospective studies to validate this region as a potential biomarker and predictor of EOC susceptibility.

Methylation of leukocyte DNA and ovarian cancer: relationships with disease status and outcome

Background

Genome-wide interrogation of DNA methylation (DNAm) in blood-derived leukocytes has become feasible with the advent of CpG genotyping arrays. In epithelial ovarian cancer (EOC), one report found substantial DNAm differences between cases and controls; however, many of these disease-associated CpGs were attributed to differences in white blood cell type distributions.

Methods

We examined blood-based DNAm in 336 EOC cases and 398 controls; we included only high-quality CpG loci that did not show evidence of association with white blood cell type distributions to evaluate association with case status and overall survival.

Results

Of 13,816 CpGs, no significant associations were observed with survival, although eight CpGs associated with survival at p < 10^-3, including methylation within a CpG island located in the promoter region of GABRE (p = 5.38 x 10^-5, HR = 0.95). In contrast, 53 CpG methylation sites were significantly associated with EOC risk (p <5 x10^-6). The top association was observed for the methylation probe cg04834572 located approximately 315 kb upstream of DUSP13 (p = 1.6 x10^-14). Other disease-associated CpGs included those near or within HHIP (cg14580567; p =5.6x10^-11), HDAC3 (cg10414058; p = 6.3x10^-12), and SCR (cg05498681; p = 4.8x10^-7).

Conclusions

We have identified several CpGs in leukocytes that are differentially methylated by case-control status. Since a retrospective study design was used, we cannot differentiate whether DNAm was etiologic or resulting from EOC; thus, prospective studies of EOC-associated loci are the critical next step.

Tumor hypomethylation at 6p21.3 associates with longer time to recurrence of high-grade serous epithelial ovarian cancer

To reveal biologic mechanisms underlying clinical outcome of high-grade serous (HGS) epithelial ovarian carcinomas (EOC), we evaluated the association between tumor epigenetic changes and time to recurrence (TTR). We assessed methylation at approximately 450,000 genome-wide CpGs in tumors of 337 Mayo Clinic (Rochester, MN) patients. Semi-supervised clustering of discovery (n=168) and validation (n=169) sets was used to determine clinically relevant methylation classes. Clustering identified two methylation classes based on 60 informative CpGs, which differed in TTR in the validation set [R vs. L class, P=2.9×10(-3), HR=0.52; 95% confidence interval (CI), 0.34-0.80]. Follow-up analyses considered genome-wide tumor mRNA expression (n=104) and CD8 T-cell infiltration (n=89) in patient subsets. Hypomethylation of CpGs located in 6p21.3 in the R class associated with cis upregulation of genes enriched in immune response processes (TAP1, PSMB8, PSMB9, HLA-DQB1, HLA-DQB2, HLA-DMA, and HLA-DOA), increased CD8 T-cell tumor infiltration (P=7.6×10(-5)), and trans-regulation of genes in immune-related pathways (P=1.6×10(-32)). This is the most comprehensive assessment of clinical outcomes with regard to epithelial ovarian carcinoma tumor methylation to date. Collectively, these results suggest that an epigenetically mediated immune response is a predictor of recurrence and, possibly, treatment response for HGS EOC.

Integrative genomic analysis identifies epigenetic marks that mediate genetic risk for epithelial ovarian cancer

BACKGROUND

Both genetic and epigenetic factors influence the development and progression of epithelial ovarian cancer (EOC). However, there is an incomplete understanding of the interrelationship between these factors and the extent to which they interact to impact disease risk. In the present study, we aimed to gain insight into this relationship by identifying DNA methylation marks that are candidate mediators of ovarian cancer genetic risk.

METHODS

We used 214 cases and 214 age-matched controls from the Mayo Clinic Ovarian Cancer Study. Pretreatment, blood-derived DNA was profiled for genome-wide methylation (Illumina Infinium HumanMethylation27 BeadArray) and single nucleotide polymorphisms (SNPs, Illumina Infinium HD Human610-Quad BeadArray). The Causal Inference Test (CIT) was implemented to distinguish CpG sites that mediate genetic risk, from those that are consequential or independently acted on by genotype.

RESULTS

Controlling for the estimated distribution of immune cells and other key covariates, our initial epigenome-wide association analysis revealed 1,993 significantly differentially methylated CpGs that between cases and controls (FDR, q < 0.05). The relationship between methylation and case-control status for these 1,993 CpGs was found to be highly consistent with the results of previously published, independent study that consisted of peripheral blood DNA methylation signatures in 131 pretreatment cases and 274 controls. Implementation of the CIT test revealed 17 CpG/SNP pairs, comprising 13 unique CpGs and 17 unique SNPs, which represent potential methylation-mediated relationships between genotype and EOC risk. Of these 13 CpGs, several are associated with immune related genes and genes that have been previously shown to exhibit altered expression in the context of cancer.

CONCLUSIONS

These findings provide additional insight into EOC etiology and may serve as novel biomarkers for EOC susceptibility.

Identification of novel variants in colorectal cancer families by high-throughput exome sequencing

BACKGROUND

Colorectal cancer (CRC) in densely affected families without Lynch Syndrome may be due to mutations in undiscovered genetic loci. Familial linkage analyses have yielded disparate results; the use of exome sequencing in coding regions may identify novel segregating variants.

METHODS

We completed exome sequencing on 40 affected cases from 16 multicase pedigrees to identify novel loci. Variants shared among all sequenced cases within each family were identified and filtered to exclude common variants and single-nucleotide variants (SNV) predicted to be benign.

RESULTS

We identified 32 nonsense or splice-site SNVs, 375 missense SNVs, 1,394 synonymous or noncoding SNVs, and 50 indels in the 16 families. Of particular interest are two validated and replicated missense variants in CENPE and KIF23, which are both located within previously reported CRC linkage regions, on chromosomes 1 and 15, respectively.

CONCLUSIONS

Whole-exome sequencing identified DNA variants in multiple genes. Additional sequencing of these genes in additional samples will further elucidate the role of variants in these regions in CRC susceptibility.

IMPACT

Exome sequencing of familial CRC cases can identify novel rare variants that may influence disease risk.

Epigenome-wide ovarian cancer analysis identifies a methylation profile differentiating clear-cell histology with epigenetic silencing of the HERG K+ channel

Ovarian cancer remains the leading cause of death in women with gynecologic malignancies, despite surgical advances and the development of more effective chemotherapeutics. As increasing evidence indicates that clear-cell ovarian cancer may have unique pathogenesis, further understanding of molecular features may enable us to begin to understand the underlying biology and histology-specific information for improved outcomes. To study epigenetics in clear-cell ovarian cancer, fresh frozen tumor DNA (n = 485) was assayed on Illumina Infinium HumanMethylation450 BeadChips. We identified a clear-cell ovarian cancer tumor methylation profile (n = 163) which we validated in two independent replication sets (set 1, n = 163; set 2, n = 159), highlighting 22 CpG loci associated with nine genes (VWA1, FOXP1, FGFRL1, LINC00340, KCNH2, ANK1, ATXN2, NDRG21 and SLC16A11). Nearly all of the differentially methylated CpGs showed a propensity toward hypermethylation among clear-cell cases. Several loci methylation inversely correlated with tumor gene expression, most notably KCNH2 (HERG, a potassium channel) (P = 9.5 × 10(-7)), indicating epigenetic silencing. In addition, a predicted methylation class mainly represented by the clear-cell cases (20 clear cell out of 23 cases) had improved survival time. Although these analyses included only 30 clear-cell carcinomas, results suggest that loss of expression of KCNH2 (HERG) by methylation could be a good prognostic marker, given that overexpression of the potassium (K(+)) channel Eag family members promotes increased proliferation and results in poor prognosis. Validation in a bigger cohort of clear-cell tumors of the ovary is warranted.

GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer

Genome-wide association studies (GWAS) have identified four susceptibility loci for epithelial ovarian cancer (EOC), with another two suggestive loci reaching near genome-wide significance. We pooled data from a GWAS conducted in North America with another GWAS from the UK. We selected the top 24,551 SNPs for inclusion on the iCOGS custom genotyping array. We performed follow-up genotyping in 18,174 individuals with EOC (cases) and 26,134 controls from 43 studies from the Ovarian Cancer Association Consortium. We validated the two loci at 3q25 and 17q21 that were previously found to have associations close to genome-wide significance and identified three loci newly associated with risk: two loci associated with all EOC subtypes at 8q21 (rs11782652, P = 5.5 × 10(-9)) and 10p12 (rs1243180, P = 1.8 × 10(-8)) and another locus specific to the serous subtype at 17q12 (rs757210, P = 8.1 × 10(-10)). An integrated molecular analysis of genes and regulatory regions at these loci provided evidence for functional mechanisms underlying susceptibility and implicated CHMP4C in the pathogenesis of ovarian cancer.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Survival prediction based on inherited gene variation analysis

There is a significant variation of outcome among ovarian cases. Clinical features such as age, stage, comorbidities, or degree of debulking are known prognostic factors for the disease. However, additional variation remains unexplained, some of which may be due to inherited factors. Here, we describe identification of survival-associated inherited variants in ovarian cancer that can enhance our current prognostic capabilities.

Identification and molecular characterization of a new ovarian cancer susceptibility locus at 17q21.31

Epithelial ovarian cancer (EOC) has a heritable component that remains to be fully characterized. Most identified common susceptibility variants lie in non-protein-coding sequences. We hypothesized that variants in the 3' untranslated region at putative microRNA (miRNA)-binding sites represent functional targets that influence EOC susceptibility. Here, we evaluate the association between 767 miRNA-related single-nucleotide polymorphisms (miRSNPs) and EOC risk in 18,174 EOC cases and 26,134 controls from 43 studies genotyped through the Collaborative Oncological Gene-environment Study. We identify several miRSNPs associated with invasive serous EOC risk (odds ratio=1.12, P=10(-8)) mapping to an inversion polymorphism at 17q21.31. Additional genotyping of non-miRSNPs at 17q21.31 reveals stronger signals outside the inversion (P=10(-10)). Variation at 17q21.31 is associated with neurological diseases, and our collaboration is the first to report an association with EOC susceptibility. An integrated molecular analysis in this region provides evidence for ARHGAP27 and PLEKHM1 as candidate EOC susceptibility genes.

Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Colorectal cancer linkage on chromosomes 4q21, 8q13, 12q24, and 15q22

A substantial proportion of familial colorectal cancer (CRC) is not a consequence of known susceptibility loci, such as mismatch repair (MMR) genes, supporting the existence of additional loci. To identify novel CRC loci, we conducted a genome-wide linkage scan in 356 white families with no evidence of defective MMR (i.e., no loss of tumor expression of MMR proteins, no microsatellite instability (MSI)-high tumors, or no evidence of linkage to MMR genes). Families were ascertained via the Colon Cancer Family Registry multi-site NCI-supported consortium (Colon CFR), the City of Hope Comprehensive Cancer Center, and Memorial University of Newfoundland. A total of 1,612 individuals (average 5.0 per family including 2.2 affected) were genotyped using genome-wide single nucleotide polymorphism linkage arrays; parametric and non-parametric linkage analysis used MERLIN in a priori-defined family groups. Five lod scores greater than 3.0 were observed assuming heterogeneity. The greatest were among families with mean age of diagnosis less than 50 years at 4q21.1 (dominant HLOD = 4.51, α = 0.84, 145.40 cM, rs10518142) and among all families at 12q24.32 (dominant HLOD = 3.60, α = 0.48, 285.15 cM, rs952093). Among families with four or more affected individuals and among clinic-based families, a common peak was observed at 15q22.31 (101.40 cM, rs1477798; dominant HLOD = 3.07, α = 0.29; dominant HLOD = 3.03, α = 0.32, respectively). Analysis of families with only two affected individuals yielded a peak at 8q13.2 (recessive HLOD = 3.02, α = 0.51, 132.52 cM, rs1319036). These previously unreported linkage peaks demonstrate the continued utility of family-based data in complex traits and suggest that new CRC risk alleles remain to be elucidated.

Abstract 5008: Genome-wide methylation and ovarian cancer survival

Background Ovarian cancer is the fifth most common cause of cancer death in women in the US; approximately 15,000 women will die from this disease in 2011 (ACS, 2011). The 5 year survival rate for stage IV disease is only 18%. Only a handful of independent significant prognostic factors are known, including tumor stage at surgery, size of residual tumor at the end of primary surgery, and age at diagnosis. Understanding how the variation of cytosine methylation in the context of cytosine-guanine dinucleotides (CpGs) affects complex disease etiology and survival is the current focus of many ongoing research studies. We conducted an epigenome-wide study of methylation in ovarian cancer cases to evaluate whether variation in methylation of CpG sites in genomic DNA may be associated with ovarian cancer survival. Methods Subjects included cases from an ongoing case-control study. The Mayo Clinic Ovarian Cancer Study is a clinic-based study that recruits ovarian cancer cases within 12 months of diagnosis. Pre-treatment, germline DNA from lymphocytes from 137 cases was bisulfite-modified, amplified and genotyped on the Illumina 450K methylation chip and included in these analyses. Quality control (QC) analyses include the elimination of all methylation probes with a single nucleotide polymorphism (SNP) or insertion/deletion (in/del) at the same location; removal of non-specific probes; and estimation of bisulfite conversion rates using embedded control probes. Planned analyses include assessment of clinical covariates (age, stage, extent of surgical debulking, grade, presence of ascites) to determine which are associated with overall survival in our data. Additional analyses will assess the age- and multi-variate-adjusted association of each methylation probe with overall survival using Cox proportional hazards regression models with each methylation probe modeled as a continuous variable. Data will be summarized using Manhattan plots. These analyses will be repeated using ‘time to recurrence’ as the outcome variable, with recurrence defined as time to ovarian cancer recurrence/progression or death. Results Bioinformatics analyses suggest that 36,659 probes were located at the site of a SNP or in/del and should be removed from future analyses. Another 448 probes had a detection rate of p&lt; 0.01 and 6630 probes had signals for the negative control samples; all were removed from analyses. A total of 43,802 probes were removed. In addition, ten cases were excluded due to either poor bisulfite conversion (N=9) or extremely low raw methylation values across all probes (N=1). This leaves a total of 127 cases for analyses. Results from planned analyses will be presented, including estimation of the association of methylation probes with both overall survival and time to recurrence. Conclusions Data to be presented will describe whether there is any evidence that methylation of genomic DNA in ovarian cancer cases is associated with outcome of ovarian cancer.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5008. doi:1538-7445.AM2012-5008

Abstract 2603: Subtype-specific ovarian cancer risk associated with SNPs in IL1A

As pathways of the immune system play an important role in multiple malignancies, we assessed the association between single-nucleotide polymorphisms (SNPs) in inflammation-related genes and risk of ovarian cancer in a multi-site case-control study. Using a custom array, we genotyped 930 epithelial ovarian cancer cases and 1,037 controls for SNPs in ALOX12, ALOX15, ALOX5, CCL11, CCL2, CCL3, CCR3, CRP, CXCL16, IL10, IL15RA, IL18, IL1A, IL1B, IL1RN, IL4R, IL6, IL6R, IL7R, IL8RA, IL8RB, IL9, NOS3, PTGS1, PTGS2, TLR2, and TNF. After analysis by logistic regression, SNPs with p&lt;0.10 were evaluated in an independent data set of 3,143 cases and 2,102 controls from the Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) collaboration. Combined analysis revealed association with interleukin 1, alpha (IL1A) SNPs (rs17561 and rs4848300, r2=0.990, p=0.005) which varied by histological subtype (heterogeneity p=0.03). For example, IL1A rs17561, which results in a missense change in amino acid at position 114 and correlates with multiple inflammatory phenotypes, was inversely associated with risk of clear cell, mucinous, and endometrioid subtypes, but not with the serous subtype (Table 1). Additional data on IL1A rs17561 in an even larger collection of ovarian cancer cases and controls are pending and will be analyzed soon. We also found that genotype at rs1864414 in ALOX5 (arachidonate 5-lipoxygenase) was associated with decreased risk combined across subtypes (OR 0.86, 95% CI 0.78-0.95; p=0.002). Thus, inherited variation in IL1A and ALOX5 appears to have a role in ovarian cancer risk, though it is limited to rarer subtypes for IL1A. These data add to existing literature on the importance of inflammation in tumorigenesis and growing evidence of subtype-specific traits in ovarian cancer and may provide important clues on prevention.

Abstract 2619: Relevance of mismatch repair deficiency in ovarian cancer

Ovarian cancer is the leading cause of death in gynecologic malignancies in the United States with most patients presenting with late stage disease; patients receive the same treatment regardless of histology or tumor characteristics. The mismatch repair (MMR) molecular pathway is well-defined and known to be involved in up to 20% of ovarian cancers; however, association of deficient v proficient MMR status (dMMR v pMMR status) with risk factors and outcome is not well understood. We are conducting extensive molecular testing of the MMR genes MLH1, MSH2, MSH3, MSH6, and PMS2 on Mayo Clinic invasive epithelial ovarian cancer cases. First, the prevalence of dMMR tumors will be estimated in ovarian cancer cases (n=496) based on loss of expression of MLH1, MSH2, MSH3, MSH6, and PMS2 determined by immunohistochemistry staining of primary antibodies using negative and positive controls. Staining of tissue microarrays (TMAs) is complete; triplicate cores will be scored independently for presence or loss of protein expression by a board-certified pathologist. Second, for those dMMR cases, the possible cause of MMR deficiency will be evaluated considering germline and somatic mutations. The loss of expression of each protein will be characterized as germline if a deleterious mutation is present in peripheral blood or somatic if the pathogenic mutation is only present in tumor tissue. In cases with MLH1 loss somatic hypermethylation will also be considered. Analysis of germline mutation data will include all 496 cases and somatic mutations data will include a subset of cases tested with frozen tumor DNA (n=292) using Fluidigm access array targeted resequencing. Third, the clinical features of dMMR disease will be compared to those of pMMR; in particular, age at diagnosis, histopathology, and disease outcome. Among cases in this analysis, the majority have high-grade serous disease with a mean age at diagnosis of 61 years, about half had recurrence or progression (n=243) of which two-thirds (n=167) are deceased. Finally, known and suspected ovarian cancer risk factors (family history, smoking, body mass index, parity, use of hormone therapy, and genetic variation in susceptibility loci at 2q31, 3q25, 8q24, 9p22, 17q21, and 19p13) will be analyzed by MMR phenotype. Risk factor questionnaire data is available on 462 cases (93%), as well as 737 population controls with germline genome-wide genotype data. About half of the cases (n=231) report a 1st or 2nd degree family member with breast or ovarian cancer. Given the limited knowledge about ovarian cancer etiology and the poor prognosis of advanced disease, it is very important to study risk factors with this molecular subtype. The ovarian cancer dMMR phenotype may have distinct clinicopathological features that could guide risk stratification and be useful in diagnostics and therapeutics. This is the largest study characterizing MMR in an ovarian cancer population to date.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2619. doi:1538-7445.AM2012-2619

Ovarian cancer risk associated with inherited inflammation-related variants

The importance of inflammation pathways to the development of many human cancers prompted us to examine the associations between single-nucleotide polymorphisms (SNP) in inflammation-related genes and risk of ovarian cancer. In a multisite case-control study, we genotyped SNPs in a large panel of inflammatory genes in 930 epithelial ovarian cancer cases and 1,037 controls using a custom array and analyzed by logistic regression. SNPs with P < 0.10 were evaluated among 3,143 cases and 2,102 controls from the Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) post-GWAS collaboration. Combined analysis revealed association with SNPs rs17561 and rs4848300 in the interleukin gene IL1A which varied by histologic subtype (P(heterogeneity) = 0.03). For example, IL1A rs17561, which correlates with numerous inflammatory phenotypes, was associated with decreased risk of clear cell, mucinous, and endometrioid subtype, but not with the most common serous subtype. Genotype at rs1864414 in the arachidonate 5-lipoxygenase ALOX5 was also associated with decreased risk. Thus, inherited variation in IL1A and ALOX5 seems to affect ovarian cancer risk which, for IL1A, is limited to rarer subtypes. Given the importance of inflammation in tumorigenesis and growing evidence of subtype-specific features in ovarian cancer, functional investigations will be important to help clarify the importance of inherited variation related to inflammation in ovarian carcinogenesis.

Germline PKHD1 mutations are protective against colorectal cancer

The autosomal recessive polycystic kidney disease (ARPKD) gene, PKHD1, has been implicated in the genesis or growth of colorectal adenocarcinoma, as a high level of somatic mutations was found in colorectal tumor tissue. To determine whether carriers of a single PKHD1 mutation are at increased risk of colorectal carcinoma, we assessed the prevalence of the commonest European mutation, T36M. First, we assayed a European cohort of ARPKD patients and found T36M was responsible for 13.1% of mutations. We then investigated two European cohorts with colorectal adenocarcinoma versus two control cohorts of similar age and gender. Screening for the most common PKHD1 mutation, T36M, we detected 15:3,603 (0.42%) controls versus 1:3,767 (0.027%) colorectal cancer individuals, indicating that heterozygous PKHD1 mutations are not a risk factor and are protective (p=0.0002). We also show that the carriage rate for PKHD1 mutations in the European population is higher than previous accepted at 3.2% (1:31 genomes).

Functional and clinical significance of variants localized to 8q24 in colon cancer

Multiple genome-wide association studies have identified several susceptibility variants for colon cancer at 8q24. However, the functional roles of these variants have yet to be elucidated. Here, we evaluated the potential role of these markers in tumor progression and examined association with commonly observed structural abnormalities in this region, c-MYC amplification and chromosome fragility at FRA8C and FRA8D. We first replicated the previously reported association by testing 1,178 cases and 1,009 clinic-based controls with eight markers localized to three specific regions at 8q24. We observed significant associations with colon cancer risk with markers rs13254738 (ordinal odds ratio, 0.82; 95% confidence interval, 0.072-0.94; P(trend) = 0.0037) and rs6983267 (ordinal odds ratio, 1.17; 95% confidence interval, 1.03-1.32, P(trend) = 0.013). Survival analysis was done using a separate set of 460 cases to evaluate the clinical significance of these markers. Overall, univariate analysis did not detect survival differences for any of the markers. We also tested a subset of the 460 cases (n = 380) for structural abnormalities at or near the c-MYC locus using fluorescence in situ hybridization analysis. Furthermore, we evaluated a small number of cases homozygous for the rs6983267 alleles to test for differences in fragile site induction. None of the 8q markers correlated with amplification at the c-MYC locus as detected by fluorescence in situ hybridization, and no clear pattern of breakage was observed at the FRA8C and FRA8D sites. In this study, we confirm the association for several single nucleotide polymorphisms at 8q24 in colon cancer but have not detected any structural role relating to c-MYC amplification or chromosomal fragility. Finally, these risk alleles do not seem to be associated with survival.

Carboxypeptidase 4 gene variants and early-onset intermediate-to-high risk prostate cancer

Background

Carboxypeptidase 4 (CPA4) is a zinc-dependent metallocarboxypeptidase on chromosome 7q32 in a region linked to prostate cancer aggressiveness. CPA4 is involved in the histone hyperacetylation pathway and may modulate the function of peptides that affect the growth and regulation of prostate epithelial cells. We examined the association between genetic variation in CPA4 and intermediate-to-high risk prostate cancer.

Methods

We studied 1012 men (506 cases and 506 controls) from Cleveland, Ohio. All cases had Gleason ≥ 7, clinical stage ≥ T2c, or PSA ≥ 10 ng/mL at diagnosis. Six CPA4 single-nucleotide polymorphisms were genotyped, and evaluated for their relation to prostate cancer. We also evaluated whether CPA4 variants influence risk of disease among men diagnosed at an earlier age (< 66 years).

Results

The nonsynonymous coding SNP (rs2171492, Cys303Gly) in CPA4 was associated with an increased risk of aggressive prostate cancer among younger patients (< 66 years). Specifically, men carrying the TT genotype had an approximately two-fold increased risk for being diagnosed with intermediate-to-high risk disease (Odds Ratio = 1.83, p = 0.04). In the overall population (all ages) none of the CPA4 SNPs demonstrated a statistically significant association with prostate cancer.

Conclusion

Coding variation in CPA4 may confer increased risk of intermediate-to-high risk prostate cancer among younger patients. Further work is needed to identify the functional aspects of this variation and understand its biological effects on prostate cancer. Such work may translate into more precise screening of higher risk individuals as well as guiding clinicians and patients toward earlier and more definitive treatment modalities in patients genetically identified as higher risk.

BRMS1 contributes to the negative regulation of uPA gene expression through recruitment of HDAC1 to the NF-kappaB binding site of the uPA promoter

The BRMS1 metastasis suppressor was recently shown to negatively regulate NF-kappaB signaling and down regulate NF-kappaB-dependent uPA expression. Here we confirm that BRMS1 expression correlates with reduced NF-kappaB DNA binding activity in independently derived human melanoma C8161.9 cells stably expressing BRMS1. We show that knockdown of BRMS1 expression in these cells using small interfering RNA (siRNA) leads to the reactivation of NF-kappaB DNA binding activity and re-expression of uPA. Further, we confirm that BRMS1 expression does not alter IKKbeta kinase activity suggesting that BRMS1-dependent uPA regulation does not occur through inhibition of the classical upstream activators of NF-kappaB. BRMS1 has been implicated as a corepressor of HDAC1 and consistent with this, we show that BRMS1 promotes HDAC1 recruitment to the NF-kappaB binding site of the uPA promoter and is associated with reduced H3 acetylation. We also confirm that BRMS1 expression stimulates disassociation of p65 from the NF-kappaB binding site of the uPA promoter consistent with its reduced DNA binding activity. These data suggest that BRMS1 recruits HDAC1 to the NF-kappaB binding site of the uPA promoter, modulates histone acetylation of p65 on the uPA promoter, leading to reduced NF-kappaB binding activity on its consensus sequence, and reduced transactivation of uPA expression.

Association of vitamin D receptor gene variants, adiposity and colon cancer

Vitamin D receptor (VDR) gene variants have been variably associated with risk of colon cancer in epidemiologic studies. We sought to further clarify the relationship between colon cancer and three single-nucleotide polymorphisms (SNPs) in the VDR gene (Cdx-2, FokI and TaqI) in a population-based case-control study of 250 incident cases and 246 controls. Colon cancer cases were more frequently homozygous for the Cdx-2 A allele (9.2 versus 4.1%, P = 0.06). Cdx-2 AA homozygotes were at increased risk with an unadjusted odds ratio (OR) of 2.47 [95% confidence interval (CI): 1.13-5.37, P = 0.022]; adjustment for age, sex, body mass index (BMI), non-steroidal anti-inflammatory use and family history of colorectal cancer yielded an OR of 2.27 (CI: 0.95-5.41, P = 0.065). Carriers of the FokI TT genotype were also at increased risk with an adjusted OR of 1.87 (CI: 1.03-3.38, P = 0.038). Haplotype analyses showed significant increased colon cancer risk for carriers of the Cdx-2-FokI A-T haplotype and the FokI-TaqI T-G haplotype. The three-SNP Cdx-2-FokI-TaqI (A-T-G) haplotype showed a similar association with an adjusted OR of 3.63 (CI: 1.01-13.07). A strong positive association was observed for the Cdx-2 variant among individuals with low BMI or low waist circumference. Our results suggest that genetic variation at the VDR locus, in particular Cdx-2 and FokI SNPs, may influence colon cancer risk and these associations may be modified by adiposity.

Evaluation of genetic variations in the androgen and estrogen metabolic pathways as risk factors for sporadic and familial prostate cancer

Previous studies suggest that enzymes involved in the androgen metabolic pathway are susceptibility factors for prostate cancer. Estrogen metabolites functioning as genotoxins have also been proposed as risk factors. In this study, we systematically tested the hypothesis that common genetic variations for those enzymes involved in the androgen and estrogen metabolic pathways increase risk for sporadic and familial prostate cancer. From these two pathways, 46 polymorphisms (34 single nucleotide polymorphisms, 10 short tandem repeat polymorphisms, and 2 null alleles) in 25 genes were tested for possible associations. Those genes tested included PRL, LHB, CYP11A1, HSD3B1, HSD3B2, HSD17B2, CYP17, SRD5A2, AKR1C3, UGT2B15, AR, SHBG, and KLK3 from the androgen pathway and CYP19, HSD17B1, CYP1A1, CYP1A2, CYP1B1, COMT, GSTP1, GSTT1, GSTM1, NQO1, ESR1, and ESR2 from the estrogen pathway. A case-control study design was used with two sets of cases: familial cases with a strong prostate cancer family history (n = 438 from 178 families) and sporadic cases with a negative prostate cancer family history (n = 499). The controls (n = 493) were derived from a population-based collection. Our results provide suggestive findings for an association with either familial or sporadic prostate cancer with polymorphisms in four genes: AKR1C3, HSD17B1, NQO1, and GSTT1. Additional suggestive findings for an association with clinical variables (disease stage, grade, and/or node status) were observed for single nucleotide polymorphisms in eight genes: HSD3B2, SRD5A2, SHBG, ESR1, CYP1A1, CYP1B1, GSTT1, and NQO1. However, none of the findings were statistically significant after appropriate corrections for multiple comparisons. Given that the point estimates for the odds ratio for each of these polymorphisms are <2.0, much larger sample sizes will be required for confirmation.

Vitamin D receptor genotypes/haplotypes and prostate cancer risk

The vitamin D receptor (VDR) gene has been associated with prostate cancer, although previous results are somewhat equivocal. To further study this, we did a family-based case-control study (N = 918) of the association between prostate cancer and six common VDR variants: Cdx2, FokI, BsmI, ApaI, TaqI, and the poly-A microsatellite. Looking at each variant alone, only FokI and ApaI were associated with disease. The FokI FF genotype was inversely associated with prostate cancer among men with less advanced disease (i.e., Gleason score <7 and tumor stage <T2c), where the odds ratio OR was 0.56 [95% confidence interval (95% CI), 0.31-1.01; P = 0.05]. ApaI, carrying one or two copies of the A allele, exhibited a weak inverse association with disease (OR, 0.64; 95% CI, 0.39-1.03; P = 0.06); this association was strengthened in Caucasian men with more advanced disease (OR, 0.44; 95% CI, 0.21-0.93; P = 0.03). We observed inverse associations between disease and the four-locus FBAt haplotype (OR, 0.48; 95% CI, 0.30-0.76; P = 0.002) and the fbaT haplotype (OR, 0.60; 95% CI, 0.38-0.95; P = 0.03; i.e., in comparison with the FbaT haplotype). These were stronger among men with more advanced disease: for FBAt, the OR was 0.31 (95% CI, 0.16-0.61; P = 0.0008); for fbaT, the OR was 0.32 (95% CI, 0.16-0.64; P = 0.001). These observations support a role for VDR variants in prostate cancer risk but suggest that any potential causal variant(s) may reside on the haplotypes reported here. This would help explain the somewhat equivocal results for VDR genotype-level associations with prostate cancer.

Association of testis derived transcript gene variants and prostate cancer risk

PURPOSE

The testis derived transcript gene has been suggested as a tumor suppressor gene for prostate cancer at 7q31. To investigate this concept we evaluated the effects of 7 tagging single nucleotide polymorphisms that comprehensively captured the common genetic variants in TES on aggressive prostate cancer in a case-control study.

MATERIALS AND METHODS

A total of 506 cases diagnosed with aggressive prostate cancer, and an equal number of age, institute and ethnicity matched controls, were recruited from the major medical institutions in Cleveland, Ohio. A logistic regression model was used to evaluate the association between SNPs/multimarker haplotypes and prostate cancer.

RESULTS

When looking at all study subjects and white men only, no statistically significant associations were observed between any variants and more aggressive disease. However, 3 variants showed inverse associations with disease in black men (178), including 2 intronic SNPs (rs2402056, rs1004109) and 1 SNP close to the 3' untranslated region (rs4730721) with ORs of 0.57 (95% CI 0.36-0.90, under an additive mode of inheritance), 0.57 (95% CI 0.36-0.91, under an additive mode of inheritance) and 0.45 (95% CI 0.21-0.98, under a dominant mode of inheritance), respectively. Variants rs2402056 and rs1004109 are in tight linkage disequilibrium (r2=0.8) and the reconstructed haplotype did not provide any additional evidence for association than their genotype level results.

CONCLUSIONS

Our findings suggest that the variants in TES, or in a nearby gene, may be associated with prostate cancer in black men.

Polymorphisms in polycyclic aromatic hydrocarbon metabolism and conjugation genes, interactions with smoking and prostate cancer risk

The relationship between cigarette smoking and prostate cancer remains unclear. Any potential association may depend on the individuals' ability to metabolize and detoxify cigarette carcinogens--such as polycyclic aromatic hydrocarbons. To investigate this, we studied the association between prostate cancer and smoking, as well as the main and modifying effects of functional polymorphisms in genes that metabolize polycyclic aromatic hydrocarbons (CYP1A1 Ile(462)Val, microsomal epoxide hydrolase His(139)Arg) and detoxify reactive derivatives (GSTM1 null deletion, GSTT1 null deletion, GSTP1 Ile(105)Val and Ala(114)Val) using a family-based case-control design (439 prostate cancer cases and 479 brother controls). Within the entire study population, there were no main effects for smoking or any of the polymorphisms. However, the nondeleted GSTM1 allele was inversely associated with prostate cancer [odds ratio (OR), 0.50; 95% confidence interval (95% CI), 0.26-0.94] among men with less aggressive disease (Gleason score < 7 and clinical tumor stage < T2c) and positively associated (OR, 1.68; 95% CI, 1.01-2.79) with prostate cancer in men with more aggressive disease (Gleason score > or = 7 or clinical tumor stage > or = T2c). We also found a statistically significant negative multiplicative interaction between the GSTM1 nondeleted allele and heavy smoking (> 20 pack-years) in the total study population (P = 0.01) and in Caucasians (P = 0.01). Among Caucasians, heavy smoking increased prostate cancer risk nearly 2-fold in those with the GSTM1 null genotype (OR, 1.73; 95% CI, 0.99-3.05) but this increased risk was not observed in heavy smokers who carried the GSTM1 nondeleted allele (OR, 0.95; 95% CI, 0.53-1.71). Our results highlight the importance of considering genetic modifiers of carcinogens when evaluating smoking in prostate cancer.

Polymorphisms in estrogen bioactivation, detoxification and oxidative DNA base excision repair genes and prostate cancer risk

To date, the potential impact of hormones on prostate cancer has predominantly focused on receptor-mediated events. However, catechol estrogens, if not inactivated by catechol-O-methyltransferase (COMT), can generate large quantities of reactive oxygen species (ROS). ROS may cause a spectrum of damage including oxidative DNA base lesions, which can lead to irreversible mutation(s) if they are not repaired by base excision repair (BER) systems. hOGG1 is a key enzyme in short patch BER because it recognizes and performs initial excision of the most common form of oxidative DNA base damage, 8-hydroxyguanine (8-oxo-dG). To investigate potential non-receptor-mediated estrogen effects, we evaluated the association between COMT Val158Met and hOGG1 Ser326Cys polymorphisms and prostate cancer in a family-based case-control study (439 prostate cancer cases, 479 brother controls). We observed no noteworthy associations between these polymorphisms and prostate cancer risk in the total study population. However, among men with more aggressive prostate cancer, the hOGG1 326 Cys/Cys genotype was inversely associated with disease (OR=0.30; 95% CI=0.09-0.98). Combining the lower activity CYP1B1 432 Leu/Leu or Leu/Val genotypes (which may decrease the level of catechol estrogens and ROS generated) with the hOGG1 326 Cys/Cys genotype and the XRCC1 399 Arg/Arg or Arg/Gln genotypes (which may enhance BER) resulted in an even further reduced risk in Caucasians with more aggressive disease (OR=0.09; 95% CI=0.01-0.56). Including the high-activity COMT 158Val allele to this combination also lowered aggressive prostate cancer risk but the effect was not as strong (OR=0.20; 95% CI=0.05-0.88). The decreased risk we observed with the hOGG1 326 Cys/Cys genotype confirms an earlier report and the further reduced risk found with the CYP1B1 (432 Leu/Leu or Leu/Val)-hOGG1 (326 Cys/Cys)-XRCC1 (Arg/Arg or Arg/Gln) genotype combination may lend new insights to the importance of ROS generated from non-receptor-mediated estrogenic mechanisms in more aggressive prostate cancer.

Podocalyxin variants and risk of prostate cancer and tumor aggressiveness

We previously reported linkage of a prostate cancer tumor aggressiveness locus to chromosome 7q32-q33, a region also associated with a high frequency of allelic imbalance in prostate tumors. The smallest region of allelic imbalance contains the podocalyxin-like (PODXL) gene, which we evaluate here as a candidate prostate cancer aggressiveness gene mapping to 7q32-q33. DNA from probands of linked families was examined for germ-line mutations in PODXL. A variable in-frame deletion, four missense variants and two nonsense variants were identified in linked men. Variants that affected amino acid sequence were further evaluated for association with risk of prostate cancer and tumor aggressiveness in a family-based case-control population (439 cases and 479 sibling controls). The presence of any single in-frame deletion was positively associated with prostate cancer [odds ratio (OR)=2.14, 95% confidence interval (95%CI)=1.09-4.20, P=0.03] and the presence of two copies of any deletion further increased risk (OR=2.58, 95%CI=1.23-5.45, P=0.01). This finding was strengthened when stratifying among men with more aggressive disease (high grade or stage): OR=3.04 for one deletion (95%CI=1.01-9.15) and OR=4.42 for two deletions (95%CI=1.32-14.85, P=0.02). A weak positive association was also observed between prostate cancer risk and PODXL variant 340A (in linkage disequilibrium with another variant, 587T) (OR=1.48, 95%CI=1.02-2.14, P=0.04). These results implicate PODXL as a candidate prostate cancer tumor aggressiveness gene mapping to chromosome 7q32-q33.

Role of Androgen Metabolism Genes CYP1B1, PSA/KLK3, and CYP11 in Prostate Cancer Risk and Aggressiveness

Candidate genes involved with androgen metabolism have been hypothesized to affect the risk of prostate cancer. To further investigate this, we evaluated the relationship between prostate cancer and multiple potentially functional polymorphisms in three genes involved in androgen metabolism: CYP1B1 (two single nucleotide polymorphisms: 355G/T and 4326C/G), prostate-specific antigen (PSA/KLK3 (three single nucleotide polymorphisms: -158A/G, -4643G/A, and -5412C/T), and CYP11alpha [(tttta)(n) repeat], using a moderately large (n = 918) sibling-based case-control population. When looking at all subjects combined, no association was observed between any polymorphism-or their haplotypes-and prostate cancer risk. However, among men with more aggressive prostate cancer, the CYP1B1 355G/T variant was positively associated with disease: carrying one or two T alleles gave odds ratios (OR) of 1.90 [95% confidence interval (95% CI), 1.09-3.31; P = 0.02] and 3.73 (95% CI, 1.39-10.0; P = 0.009), respectively. Similarly, carrying the CYP1B1 355T-4326C haplotype was positively associated with prostate cancer among men with high aggressive disease (P = 0.01). In addition, the PSA -158G/-158G genotype was positively associated with prostate cancer among men with less aggressive disease (OR, 2.71; 95% CI, 1.06-6.94; P = 0.04). Our findings suggest that CYP1B1 and PSA variants may affect the risk of prostate cancer and tumor aggressiveness.

Relationship between Methylenetetrahydrofolate Reductase C677T and A1298C Genotypes and Haplotypes and Prostate Cancer Risk and Aggressiveness

Previous reports indicate that polymorphisms in the MTHFR gene play a role in cancer development, but their potential impact on prostate cancer has not been well studied. Here, we evaluate the association between two MTHFR polymorphisms, C677T and A1298C, and prostate cancer risk and aggressiveness in a moderately large family-based case-control study (439 cases and 479 sibling controls). Among all study subjects, we observed no association between the C677T variant and prostate cancer but a slight positive association between the A1298C variant and risk of this disease [odds ratio (OR) 1.41, 95% confidence interval (CI) 0.96-2.06; P = 0.08]. When stratifying the study population by disease aggressiveness at diagnosis, the C677T variant was positively associated with risk among men with less advanced disease (OR 1.86, 95% CI 1.00-3.46; P = 0.05). In contrast, when looking at men with more advanced disease, the C677T variant was inversely associated with risk (OR 0.51, 95% CI 0.32-0.82; P = 0.01), whereas the A1298C variant was positively associated with risk (OR 1.79, 95% CI 1.06-3.02; P = 0.03). Furthermore, the 677T-1298A haplotype was positively associated with prostate cancer among men with less advanced disease (OR 1.84, 95% CI 1.07-3.16; P = 0.03) and inversely associated with risk of more advanced disease (OR 0.47, 95% CI 0.29-0.76; P = 0.002). Our findings suggest that 677T and 1298A, or another variant on their haplotype, may be associated with a reduced risk of progression to more advanced prostate cancer.

Relationship between methylenetetrahydrofolate reductase C677T and A1298C genotypes and haplotypes and prostate cancer risk and aggressiveness

Previous reports indicate that polymorphisms in the MTHFR gene play a role in cancer development, but their potential impact on prostate cancer has not been well studied. Here, we evaluate the association between two MTHFR polymorphisms, C677T and A1298C, and prostate cancer risk and aggressiveness in a moderately large family-based case-control study (439 cases and 479 sibling controls). Among all study subjects, we observed no association between the C677T variant and prostate cancer but a slight positive association between the A1298C variant and risk of this disease [odds ratio (OR) 1.41, 95% confidence interval (CI) 0.96-2.06; P = 0.08]. When stratifying the study population by disease aggressiveness at diagnosis, the C677T variant was positively associated with risk among men with less advanced disease (OR 1.86, 95% CI 1.00-3.46; P = 0.05). In contrast, when looking at men with more advanced disease, the C677T variant was inversely associated with risk (OR 0.51, 95% CI 0.32-0.82; P = 0.01), whereas the A1298C variant was positively associated with risk (OR 1.79, 95% CI 1.06-3.02; P = 0.03). Furthermore, the 677T-1298A haplotype was positively associated with prostate cancer among men with less advanced disease (OR 1.84, 95% CI 1.07-3.16; P = 0.03) and inversely associated with risk of more advanced disease (OR 0.47, 95% CI 0.29-0.76; P = 0.002). Our findings suggest that 677T and 1298A, or another variant on their haplotype, may be associated with a reduced risk of progression to more advanced prostate cancer.

Association of prostate cancer risk and aggressiveness to androgen pathway genes: SRD5A2, CYP17, and the AR

BACKGROUND

The prostate is an androgen-regulated organ and polymorphisms in genes involved in testosterone synthesis, in particular, SRD5A2 (A49T and V89L variants), CYP17 (MspAI variant), and the AR (CAG, GGC repeats), represent candidate risk factors for prostate cancer incidence and aggressiveness.

METHODS

We evaluated the relationship between these five polymorphisms and prostate cancer risk in a family-based case-control study (N = 920). Cases were diagnosed at major medical institutions in Cleveland Ohio, and Detroit Michigan, and their unaffected brothers were used as controls. Associations were investigated with regard to prostate cancer risk, and clinical characteristics at diagnosis (i.e., tumor stage/grade, age, family history).

RESULTS

The SRD5A2 V89L variant was associated with an increased risk of prostate cancer (OR = 1.56, P = 0.02). This association was driven primarily by men diagnosed at an earlier age (OR = 2.35, P = 0.001), or with more aggressive disease (OR = 1.63, P = 0.06). None of the other variants exhibited noteworthy associations with disease.

CONCLUSIONS

These findings suggest that the SRD5A2 V89L variant may influence risk of developing prostate cancer, especially among men with a younger age of diagnosis or more aggressive disease.