Real-world Trends, Rural-urban Differences, and Socioeconomic Disparities in Utilization of Narrow versus Broad Next-generation Sequencing Panels.

Advances in genetic technology have led to the increasing use of genomic panels in precision oncology practice, with panels ranging from a couple to hundreds of genes. However, the clinical utilization and utility of oncology genomic panels, especially among vulnerable populations, is unclear. We examined the association of panel size with socioeconomic status and clinical trial matching. We retrospectively identified 9,886 eligible adult subjects in the Mayo Clinic Health System who underwent genomic testing between January 1, 2016 and June 30, 2020. Patient data were retrieved from structured and unstructured data sources of institutional collections, including cancer registries, clinical data warehouses, and clinical notes. Socioeconomic surrogates were approximated using the Area Deprivation Index (ADI) corresponding to primary residence addresses. Logistic regression was performed to analyze relationships between ADI or rural/urban status and (i) use of genomic test by panel size; (ii) clinical trial matching status. Compared with patients from the most affluent areas, patients had a lower odds of receiving a panel test (vs. a single-gene test) if from areas of higher socioeconomic deprivation [OR (95% confidence interval (CI): 0.71 (0.61-0.83), P < 0.01] or a rural area [OR (95% CI): 0.85 (0.76-0.96), P < 0.01]. Patients in areas of higher socioeconomic deprivation were less likely to be matched to clinical trials if receiving medium panel tests [(OR) (95% CI): 0.69 (0.49-0.97), P = 0.03]; however, there was no difference among patients receiving large panel tests (P > 0.05) and rural patients were almost 2x greater odds of being matched if receiving a large panel test [(OR) (95% CI): 1.76 (1.21-2.55), P < 0.01]. SIGNIFICANCE: We identified socioeconomic and rurality disparities in the use of genomic tests and trial matching by panel size, which may have implications for equal access to targeted therapies. The lack of association between large panel tests and clinical trial matching by socioeconomic status, suggests a potential health equity impact, while removing barriers in access to large panels for rural patients may improve access to trials. However, further research is needed.

Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways.

Rare familial spontaneous coronary artery dissection (SCAD) kindreds implicate genetic disease predisposition and provide a unique opportunity for candidate gene discovery. Whole-genome sequencing was performed in fifteen probands with non-syndromic SCAD who had a relative with SCAD, eight of whom had a second relative with extra-coronary arteriopathy. Co-segregating variants and associated genes were prioritized by quantitative variant, gene, and disease-level metrics. Curated public databases were queried for functional relationships among encoded proteins. Fifty-four heterozygous coding variants in thirteen families co-segregated with disease and fulfilled primary filters of rarity, gene variation constraint, and predicted-deleterious protein effect. Secondary filters yielded 11 prioritized candidate genes in 12 families, with high arterial tissue expression (n = 7), high-confidence protein-level interactions with genes associated with SCAD previously (n = 10), and/or previous associations with connective tissue disorders and aortopathies (n = 3) or other vascular phenotypes in mice or humans (n = 11). High-confidence associations were identified among 10 familial SCAD candidate-gene-encoded proteins. A collagen-encoding gene was identified in five families, two with distinct variants in COL4A2. Familial SCAD is genetically heterogeneous, yet perturbations of extracellular matrix, cytoskeletal, and cell-cell adhesion proteins implicate common disease-susceptibility pathways. Incomplete penetrance and variable expression suggest genetic or environmental modifiers.

Mitochondrial MICOS complex genes, implicated in hypoplastic left heart syndrome, maintain cardiac contractility and actomyosin integrity.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease (CHD) with a likely oligogenic etiology, but our understanding of the genetic complexities and pathogenic mechanisms leading to HLHS is limited. We performed whole genome sequencing (WGS) on 183 HLHS patient-parent trios to identify candidate genes, which were functionally tested in the Drosophila heart model. Bioinformatic analysis of WGS data from an index family of a HLHS proband born to consanguineous parents prioritized 9 candidate genes with rare, predicted damaging homozygous variants. Of them, cardiac-specific knockdown (KD) of mitochondrial MICOS complex subunit dCHCHD3/6 resulted in drastically compromised heart contractility, diminished levels of sarcomeric actin and myosin, reduced cardiac ATP levels, and mitochondrial fission-fusion defects. These defects were similar to those inflicted by cardiac KD of ATP synthase subunits of the electron transport chain (ETC), consistent with the MICOS complex's role in maintaining cristae morphology and ETC assembly. Five additional HLHS probands harbored rare, predicted damaging variants in CHCHD3 or CHCHD6. Hypothesizing an oligogenic basis for HLHS, we tested 60 additional prioritized candidate genes from these patients for genetic interactions with CHCHD3/6 in sensitized fly hearts. Moderate KD of CHCHD3/6 in combination with Cdk12 (activator of RNA polymerase II), RNF149 (goliath, E3 ubiquitin ligase), or SPTBN1 (ß-Spectrin, scaffolding protein) caused synergistic heart defects, suggesting the likely involvement of diverse pathways in HLHS. Further elucidation of novel candidate genes and genetic interactions of potentially disease-contributing pathways is expected to lead to a better understanding of HLHS and other CHDs.

Genome-Wide Association and Inheritance-Based Analyses Implicate Unconventional Myosin Genes in Hypoplastic Left Heart Syndrome.

BACKGROUND: Deciphering hypoplastic left heart syndrome (HLHS) pathogenesis is confounded by its genetic heterogeneity and oligogenic underpinnings. METHODS: Whole genome sequences were analyzed by 3 independent strategies to identify HLHS gene candidates, ranked by variant, gene, and disease-level metrics. RESULTS: First, a genome-wide association study of 174 cases and 853 controls revealed suggestive association with a MYO18B intron 33 variant (rs2269628-G; frequency=0.55 versus 0.39; OR, 1.97 [95% CI, 1.54-2.52]; P=6.70×10-8). Second, transmission disequilibrium testing of 161 HLHS proband-parent trios revealed overrepresentation of a MYO18B intron 42 variant (rs73154186-A; frequency=0.05; OR, 24 [95% CI, 3.2-177.4]; P=4.23×10-6). Third, rare, predicted-damaging variants were filtered in 2 multiplex families. In 141H, 2 fifth-degree relatives with HLHS shared a paternally-inherited MYO5A missense variant (p.Arg801Trp; frequency=0.00003; combined annotation-dependent depletion score=29), each with a maternally-inherited or de novo candidate modifier variant in a MYO5A-interacting conventional myosin. In 442H, a HLHS proband was compound heterozygous for MYO15A variants-a maternally-inherited pathogenic stop-gain variant co-segregating with tetralogy of Fallot and bicuspid aortic valve in maternal relatives (p.Tyr2819Ter; frequency=0.00003) and a paternally-inherited intronic variant altering a canonical transcription factor binding site (rs1277068603; frequency=0.00001; position weight matrix score=0.98). CONCLUSIONS: Collectively, these findings suggest that common and rare alleles within unconventional myosin genes are associated with HLHS susceptibility. The identified candidate MYO18B regulates cardiac sarcomerogenesis, supporting the hypothesis of intrinsic myogenic perturbation in arrested left heart development.

CELSR1 Risk Alleles in Familial Bicuspid Aortic Valve and Hypoplastic Left Heart Syndrome.

BACKGROUND: Whole-genome sequencing in families enables deciphering of congenital heart disease causes. A shared genetic basis for familial bicuspid aortic valve (BAV) and hypoplastic left heart syndrome (HLHS) was postulated. METHODS: Whole-genome sequencing was performed in affected members of 6 multiplex BAV families, an HLHS cohort of 197 probands and 546 relatives, and 813 controls. Data were filtered for rare, predicted-damaging variants that cosegregated with familial BAV and disrupted genes associated with congenital heart disease in humans and mice. Candidate genes were further prioritized by rare variant burden testing in HLHS cases versus controls. Modifier variants in HLHS proband-parent trios were sought to account for the severe developmental phenotype. RESULTS: In 5 BAV families, missense variants in 6 ontologically diverse genes for structural (SPTBN1, PAXIP1, and FBLN1) and signaling (CELSR1, PLXND1, and NOS3) proteins fulfilled filtering metrics. CELSR1, encoding cadherin epidermal growth factor laminin G seven-pass G-type receptor, was identified as a candidate gene in 2 families and was the only gene demonstrating rare variant enrichment in HLHS probands (P=0.003575). HLHS-associated CELSR1 variants included 16 missense, one splice site, and 3 noncoding variants predicted to disrupt canonical transcription factor binding sites, most of which were inherited from a parent without congenital heart disease. Filtering whole-genome sequencing data for rare, predicted-damaging variants inherited from the other parent revealed 2 cases of CELSR1 compound heterozygosity, one case of CELSR1-CELSR3 synergistic heterozygosity, and 4 cases of CELSR1-MYO15A digenic heterozygosity. CONCLUSIONS: CELSR1 is a susceptibility gene for familial BAV and HLHS, further implicating planar cell polarity pathway perturbation in congenital heart disease.

Frequent POLE-driven hypermutation in ovarian endometrioid cancer revealed by mutational signatures in RNA sequencing.

BACKGROUND: DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT > TAT and TCG > TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC. METHODS: Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n = 195), The Cancer Genome Atlas (TCGA) OC tumors (n = 419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n = 84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10. RESULTS: In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P = 5e-04). These six tumors had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P = .008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations. CONCLUSIONS: From the largest cohort of RNA-seq from endometrioid OC to date (n = 53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

Patient-specific genomics and cross-species functional analysis implicate LRP2 in hypoplastic left heart syndrome.

Congenital heart diseases (CHDs), including hypoplastic left heart syndrome (HLHS), are genetically complex and poorly understood. Here, a multidisciplinary platform was established to functionally evaluate novel CHD gene candidates, based on whole-genome and iPSC RNA sequencing of a HLHS family-trio. Filtering for rare variants and altered expression in proband iPSCs prioritized 10 candidates. siRNA/RNAi-mediated knockdown in healthy human iPSC-derived cardiomyocytes (hiPSC-CM) and in developing Drosophila and zebrafish hearts revealed that LDL receptor-related protein LRP2 is required for cardiomyocyte proliferation and differentiation. Consistent with hypoplastic heart defects, compared to patents the proband's iPSC-CMs exhibited reduced proliferation. Interestingly, rare, predicted-damaging LRP2 variants were enriched in a HLHS cohort; however, understanding their contribution to HLHS requires further investigation. Collectively, we have established a multi-species high-throughput platform to rapidly evaluate candidate genes and their interactions during heart development, which are crucial first steps toward deciphering oligogenic underpinnings of CHDs, including hypoplastic left hearts.

An expanded variant list and assembly annotation identifies multiple novel coding and noncoding genes for prostate cancer risk using a normal prostate tissue eQTL data set.

Prostate cancer (PrCa) is highly heritable; 284 variants have been identified to date that are associated with increased prostate cancer risk, yet few genes contributing to its development are known. Expression quantitative trait loci (eQTL) studies link variants with affected genes, helping to determine how these variants might regulate gene expression and may influence prostate cancer risk. In the current study, we performed eQTL analysis on 471 normal prostate epithelium samples and 249 PrCa-risk variants in 196 risk loci, utilizing RNA sequencing transcriptome data based on ENSEMBL gene definition and genome-wide variant data. We identified a total of 213 genes associated with known PrCa-risk variants, including 141 protein-coding genes, 16 lncRNAs, and 56 other non-coding RNA species with differential expression. Compared to our previous analysis, where RefSeq was used for gene annotation, we identified an additional 130 expressed genes associated with known PrCa-risk variants. We detected an eQTL signal for more than half (n = 102, 52%) of the 196 loci tested; 52 (51%) of which were a Group 1 signal, indicating high linkage disequilibrium (LD) between the peak eQTL variant and the PrCa-risk variant (r2>0.5) and may help explain how risk variants influence the development of prostate cancer.

Molecular signatures of X chromosome inactivation and associations with clinical outcomes in epithelial ovarian cancer.

X chromosome inactivation (XCI) is a key epigenetic gene expression regulatory process, which may play a role in women's cancer. In particular tissues, some genes are known to escape XCI, yet patterns of XCI in ovarian cancer (OC) and their clinical associations are largely unknown. To examine XCI in OC, we integrated germline genotype with tumor copy number, gene expression and DNA methylation information from 99 OC patients. Approximately 10% of genes showed different XCI status (either escaping or being subject to XCI) compared with the studies of other tissues. Many of these genes are known oncogenes or tumor suppressors (e.g. DDX3X, TRAPPC2 and TCEANC). We also observed strong association between cis promoter DNA methylation and allele-specific expression imbalance (P = 2.0 × 10-10). Cluster analyses of the integrated data identified two molecular subgroups of OC patients representing those with regulated (N = 47) and dysregulated (N = 52) XCI. This XCI cluster membership was associated with expression of X inactive specific transcript (P = 0.002), a known driver of XCI, as well as age, grade, stage, tumor histology and extent of residual disease following surgical debulking. Patients with dysregulated XCI (N = 52) had shorter time to recurrence (HR = 2.34, P = 0.001) and overall survival time (HR = 1.87, P = 0.02) than those with regulated XCI, although results were attenuated after covariate adjustment. Similar findings were observed when restricted to high-grade serous tumors. We found evidence of a unique OC XCI profile, suggesting that XCI may play an important role in OC biology. Additional studies to examine somatic changes with paired tumor-normal tissue are needed.

An integrative approach to assess X-chromosome inactivation using allele-specific expression with applications to epithelial ovarian cancer.

X-chromosome inactivation (XCI) epigenetically silences transcription of an X chromosome in females; patterns of XCI are thought to be aberrant in women's cancers, but are understudied due to statistical challenges. We develop a two-stage statistical framework to assess skewed XCI and evaluate gene-level patterns of XCI for an individual sample by integration of RNA sequence, copy number alteration, and genotype data. Our method relies on allele-specific expression (ASE) to directly measure XCI and does not rely on male samples or paired normal tissue for comparison. We model ASE using a two-component mixture of beta distributions, allowing estimation for a given sample of the degree of skewness (based on a composite likelihood ratio test) and the posterior probability that a given gene escapes XCI (using a Bayesian beta-binomial mixture model). To illustrate the utility of our approach, we applied these methods to data from tumors of ovarian cancer patients. Among 99 patients, 45 tumors were informative for analysis and showed evidence of XCI skewed toward a particular parental chromosome. For 397 X-linked genes, we observed tumor XCI patterns largely consistent with previously identified consensus states based on multiple normal tissue types. However, 37 genes differed in XCI state between ovarian tumors and the consensus state; 17 genes aberrantly escaped XCI in ovarian tumors (including many oncogenes), whereas 20 genes were unexpectedly inactivated in ovarian tumors (including many tumor suppressor genes). These results provide evidence of the importance of XCI in ovarian cancer and demonstrate the utility of our two-stage analysis.

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

Importance: Cytotoxic CD8+ tumor-infiltrating lymphocytes (TILs) participate in immune control of epithelial ovarian cancer; however, little is known about prognostic patterns of CD8+ TILs by histotype and in relation to other clinical factors. Objective: To define the prognostic role of CD8+ TILs in epithelial ovarian cancer. Design, Setting, and Participants: This was a multicenter observational, prospective survival cohort study of the Ovarian Tumor Tissue Analysis Consortium. More than 5500 patients, including 3196 with high-grade serous ovarian carcinomas (HGSOCs), were followed prospectively for over 24 650 person-years. Exposures: Following immunohistochemical analysis, CD8+ TILs were identified within the epithelial components of tumor islets. Patients were grouped based on the estimated number of CD8+ TILs per high-powered field: negative (none), low (1-2), moderate (3-19), and high (≥20). CD8+ TILs in a subset of patients were also assessed in a quantitative, uncategorized manner, and the functional form of associations with survival was assessed using penalized B-splines. Main Outcomes and Measures: Overall survival time. Results: The final sample included 5577 women; mean age at diagnosis was 58.4 years (median, 58.2 years). Among the 5 major invasive histotypes, HGSOCs showed the most infiltration. CD8+ TILs in HGSOCs were significantly associated with longer overall survival; median survival was 2.8 years for patients with no CD8+ TILs and 3.0 years, 3.8 years, and 5.1 years for patients with low, moderate, or high levels of CD8+ TILs, respectively (P value for trend = 4.2 × 10−16). A survival benefit was also observed among women with endometrioid and mucinous carcinomas, but not for those with the other histotypes. Among HGSOCs, CD8+ TILs were favorable regardless of extent of residual disease following cytoreduction, known standard treatment, and germline BRCA1 pathogenic mutation, but were not prognostic for BRCA2 mutation carriers. Evaluation of uncategorized CD8+ TIL counts showed a near-log-linear functional form. Conclusions and Relevance: This study demonstrates the histotype-specific nature of immune infiltration and provides definitive evidence for a dose-response relationship between CD8+ TILs and HGSOC survival. That the extent of infiltration is prognostic, not merely its presence or absence, suggests that understanding factors that drive infiltration will be the key to unraveling outcome heterogeneity in this cancer.

Targeted sequencing of 36 known or putative colorectal cancer susceptibility genes.

BACKGROUND: Mutations in several genes predispose to colorectal cancer. Genetic testing for hereditary colorectal cancer syndromes was previously limited to single gene tests; thus, only a very limited number of genes were tested, and rarely those infrequently mutated in colorectal cancer. Next-generation sequencing technologies have made it possible to sequencing panels of genes known and suspected to influence colorectal cancer susceptibility. METHODS: Targeted sequencing of 36 known or putative CRC susceptibility genes was conducted for 1231 CRC cases from five subsets: (1) Familial Colorectal Cancer Type X (n = 153); (2) CRC unselected by tumor immunohistochemical or microsatellite stability testing (n = 548); (3) young onset (age <50 years) (n = 333); (4) proficient mismatch repair (MMR) in cases diagnosed at ≥50 years (n = 68); and (5) deficient MMR CRCs with no germline mutations in MLH1, MSH2, MSH6, or PMS2 (n = 129). Ninety-three unaffected controls were also sequenced. RESULTS: Overall, 29 nonsense, 43 frame-shift, 13 splice site, six initiator codon variants, one stop codon, 12 exonic deletions, 658 missense, and 17 indels were identified. Missense variants were reviewed by genetic counselors to determine pathogenicity; 13 were pathogenic, 61 were not pathogenic, and 584 were variants of uncertain significance. Overall, we identified 92 cases with pathogenic mutations in APC,MLH1,MSH2,MSH6, or multiple pathogenic MUTYH mutations (7.5%). Four cases with intact MMR protein expression by immunohistochemistry carried pathogenic MMR mutations. CONCLUSIONS: Results across case subsets may help prioritize genes for inclusion in clinical gene panel tests and underscore the issue of variants of uncertain significance both in well-characterized genes and those for which limited experience has accumulated.

Mapping Complex Traits in a Diversity Outbred F1 Mouse Population Identifies Germline Modifiers of Metastasis in Human Prostate Cancer.

It is unclear how standing genetic variation affects the prognosis of prostate cancer patients. To provide one controlled answer to this problem, we crossed a dominant, penetrant mouse model of prostate cancer to Diversity Outbred mice, a collection of animals that carries over 40 million SNPs. Integration of disease phenotype and SNP variation data in 493 F1 males identified a metastasis modifier locus on Chromosome 8 (LOD = 8.42); further analysis identified the genes Rwdd4, Cenpu, and Casp3 as functional effectors of this locus. Accordingly, analysis of over 5,300 prostate cancer patient samples revealed correlations between the presence of genetic variants at these loci, their expression levels, cancer aggressiveness, and patient survival. We also observed that ectopic overexpression of RWDD4 and CENPU increased the aggressiveness of two human prostate cancer cell lines. In aggregate, our approach demonstrates how well-characterized genetic variation in mice can be harnessed in conjunction with systems genetics approaches to identify and characterize germline modifiers of human disease processes.

Germline miRNA DNA variants and the risk of colorectal cancer by subtype.

MicroRNAs (miRNAs) regulate up to one-third of all protein-coding genes including genes relevant to cancer. Variants within miRNAs have been reported to be associated with prognosis, survival, response to chemotherapy across cancer types, in vitro parameters of cell growth, and altered risks for development of cancer. Five miRNA variants have been reported to be associated with risk for development of colorectal cancer (CRC). In this study, we evaluated germline genetic variation in 1,123 miRNAs in 899 individuals with CRCs categorized by clinical subtypes and in 204 controls. The role of common miRNA variation in CRC was investigated using single variant and miRNA-level association tests. Twenty-nine miRNAs and 30 variants exhibited some marginal association with CRC in at least one subtype of CRC. Previously reported associations were not confirmed (n = 4) or could not be evaluated (n = 1). The variants noted for the CRCs with deficient mismatch repair showed little overlap with the variants noted for CRCs with proficient mismatch repair, consistent with our evolving understanding of the distinct biology underlying these two groups. © 2016 The Authors Genes, Chromosomes & Cancer Published by Wiley Periodicals, Inc.

The inflammatory microenvironment in epithelial ovarian cancer: a role for TLR4 and MyD88 and related proteins.

The tumor-associated inflammatory microenvironment may play a pivotal role in epithelial ovarian cancer (EOC) carcinogenesis and outcomes, but a detailed profile in patient-derived tumors is needed. Here, we investigated the expression of TLR4- and MyD88-associated markers in tumors from over 500 EOC patients using immunohistochemical staining. We demonstrate that high expression of TLR4 and MyD88 predicts poorer overall survival in patients with EOC; most likely, this is due to their association with serous histology and features of high tumor burden and aggressiveness, including stage, grade, and ascites at surgery. Combined TLR4 and MyD88 expression appears to serve as an independent risk factor for shortened survival time, even after covariate adjustment (both moderate HR 1.1 [95 % CI 0.7-1.8], both strong HR 2.1 [95 % CI 1.1-3.8], both weak as referent; p = 0.027). We reveal that in EOC tissues with elevated expression of both TLR4 and MyD88 and activated NF-κB signaling pathway, expression of hsp60, hsp70, beta 2 defensin, and HMGB1 are also enriched. In total, these results suggest that activation of TLR4/MyD88/NF-κB signaling by endogenous ligands may contribute to an inflammatory microenvironment that drives a more aggressive phenotype with poorer clinical outcome in EOC patients.

A targeted genetic association study of epithelial ovarian cancer susceptibility.

BACKGROUND: Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci. RESULTS: At nine of eleven previously published EOC susceptibility regions (2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13), novel variants were identified that were more strongly associated with risk than previously reported variants. Beyond known susceptibility regions, no variants were found to be associated with EOC risk at genome-wide statistical significance (p <5x10(-8)), nor were any significant after Bonferroni correction for 17,000 variants (p< 3x10-6). METHODS: A customized genotyping array was used to assess over 17,000 variants in coding, non-coding, regulatory, and known susceptibility regions in 4,973 EOC cases and 5,640 controls from 13 independent studies. Susceptibility for EOC overall and for select histotypes was evaluated using logistic regression adjusted for age, study site, and population substructure. CONCLUSION: Given the novel variants identified within the 2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13 regions, larger follow-up genotyping studies, using imputation where necessary, are needed for fine-mapping and confirmation of low frequency variants that fall below statistical significance.

Investigation of Exomic Variants Associated with Overall Survival in Ovarian Cancer.

BACKGROUND: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). METHODS: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). RESULTS: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E-6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E-6; Pcorrected = 0.01). CONCLUSIONS: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. IMPACT: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446-54. ©2016 AACR.

Large-scale evaluation of common variation in regulatory T cell-related genes and ovarian cancer outcome.

The presence of regulatory T cells (Treg) in solid tumors is known to play a role in patient survival in ovarian cancer and other malignancies. We assessed inherited genetic variations via 749 tag single-nucleotide polymorphisms (SNP) in 25 Treg-associated genes (CD28, CTLA4, FOXP3, IDO1, IL10, IL10RA, IL15, 1L17RA, IL23A, IL23R, IL2RA, IL6, IL6R, IL8, LGALS1, LGALS9, MAP3K8, STAT5A, STAT5B, TGFB1, TGFB2, TGFB3, TGFBR1, TGRBR2, and TGFBR3) in relation to ovarian cancer survival. We analyzed genotype and overall survival in 10,084 women with invasive epithelial ovarian cancer, including 5,248 high-grade serous, 1,452 endometrioid, 795 clear cell, and 661 mucinous carcinoma cases of European descent across 28 studies from the Ovarian Cancer Association Consortium (OCAC). The strongest associations were found for endometrioid carcinoma and IL2RA SNPs rs11256497 [HR, 1.42; 95% confidence interval (CI), 1.22-1.64; P = 5.7 × 10(-6)], rs791587 (HR, 1.36; 95% CI, 1.17-1.57; P = 6.2 × 10(-5)), rs2476491 (HR, = 1.40; 95% CI, 1.19-1.64; P = 5.6 × 10(-5)), and rs10795763 (HR, 1.35; 95% CI, 1.17-1.57; P = 7.9 × 10(-5)), and for clear cell carcinoma and CTLA4 SNP rs231775 (HR, 0.67; 95% CI, 0.54-0.82; P = 9.3 × 10(-5)) after adjustment for age, study site, population stratification, stage, grade, and oral contraceptive use. The rs231775 allele associated with improved survival in our study also results in an amino acid change in CTLA4 and previously has been reported to be associated with autoimmune conditions. Thus, we found evidence that SNPs in genes related to Tregs seem to play a role in ovarian cancer survival, particularly in patients with clear cell and endometrioid epithelial ovarian cancer.

Analysis of over 10,000 Cases finds no association between previously reported candidate polymorphisms and ovarian cancer outcome.

BACKGROUND: Ovarian cancer is a leading cause of cancer-related death among women. In an effort to understand contributors to disease outcome, we evaluated single-nucleotide polymorphisms (SNP) previously associated with ovarian cancer recurrence or survival, specifically in angiogenesis, inflammation, mitosis, and drug disposition genes. METHODS: Twenty-seven SNPs in VHL, HGF, IL18, PRKACB, ABCB1, CYP2C8, ERCC2, and ERCC1 previously associated with ovarian cancer outcome were genotyped in 10,084 invasive cases from 28 studies from the Ovarian Cancer Association Consortium with over 37,000-observed person-years and 4,478 deaths. Cox proportional hazards models were used to examine the association between candidate SNPs and ovarian cancer recurrence or survival with and without adjustment for key covariates. RESULTS: We observed no association between genotype and ovarian cancer recurrence or survival for any of the SNPs examined. CONCLUSIONS: These results refute prior associations between these SNPs and ovarian cancer outcome and underscore the importance of maximally powered genetic association studies. IMPACT: These variants should not be used in prognostic models. Alternate approaches to uncovering inherited prognostic factors, if they exist, are needed.