Winham, Stacey J.

Abstract Background: Mammographic breast density declines during menopause. We assessed changes in volumetric breast density across the menopausal transition and factors that influence these changes. Methods: Women without a history of breast cancer, who had full field digital mammograms during both pre- and postmenopausal periods, at least 2 years apart, were sampled from four facilities within the San Francisco Mammography Registry from 2007 to 2013. Dense breast volume (DV) was assessed using Volpara on mammograms across the time period. Annualized change in DV from pre- to postmenopause was estimated using linear mixed models adjusted for covariates and per-woman random effects. Multiplicative interactions were evaluated between premenopausal risk factors and time to determine whether these covariates modified the annualized changes. Results: Among the 2,586 eligible women, 1,802 had one premenopausal and one postmenopausal mammogram, 628 had an additional perimenopausal mammogram, and 156 had two perimenopausal mammograms. Women experienced an annualized decrease in DV [−2.2 cm3 (95% confidence interval, −2.7 to −1.7)] over the menopausal transition. Declines were greater among women with a premenopausal DV above the median (54 cm3) versus below (DV, −3.5 cm3 vs. −1.0 cm3; P < 0.0001). Other breast cancer risk factors, including race, body mass index, family history, alcohol, and postmenopausal hormone therapy, had no effect on change in DV over the menopausal transition. Conclusions: High premenopausal DV was a strong predictor of greater reductions in DV across the menopausal transition. Impact: We found that few factors other than premenopausal density influence changes in DV across the menopausal transition, limiting targeted prevention efforts.

Abstract Background Archived formalin fixed paraffin embedded (FFPE) samples are valuable clinical resources to examine clinically relevant morphology features and also to study genetic changes. However, DNA quality and quantity of FFPE samples are often sub-optimal, and resulting NGS-based genetics variant detections are prone to false positives. Evaluations of wet-lab and bioinformatics approaches are needed to optimize variant detection from FFPE samples. Results As a pilot study, we designed within-subject triplicate samples of DNA derived from paired FFPE and fresh frozen breast tissues to highlight FFPE-specific artifacts. For FFPE samples, we tested two FFPE DNA extraction methods to determine impact of wet-lab procedures on variant calling: QIAGEN QIAamp DNA Mini Kit ("QA"), and QIAGEN GeneRead DNA FFPE Kit ("QGR"). We also used negative-control (NA12891) and positive control samples (Horizon Discovery Reference Standard FFPE). All DNA sample libraries were prepared for NGS according to the QIAseq Human Breast Cancer Targeted DNA Panel protocol and sequenced on the HiSeq 4000. Variant calling and filtering were performed using QIAGEN Gene Globe Data Portal. Detailed variant concordance comparisons and mutational signature analysis were performed to investigate effects of FFPE samples compared to paired fresh frozen samples, along with different library preparations. In this study, we found that five times or more variants were called with FFPE samples, compared to their paired fresh-frozen tissue samples even after applying molecular barcoding error-correction and default bioinformatics filtering recommended by the vendor. We also found that QGR as an optimized FFPE-DNA extraction approach leads to much fewer discordant variants between paired fresh frozen and FFPE samples. Approximately 92% of the uniquely called FFPE variants were of low allelic frequency range (<5%), and collectively shared a “C>T|G>A” mutational signature known to be representative of FFPE artifacts resulting fromcytosine deamination. Based on control samples and FFPE-frozen replicates, we derived an effective filtering strategy with associated empirical false-discovery estimates. Conclusions Through this study, we demonstrated feasibility of calling and filtering genetic variants from FFPE tissue samples using a combined strategy with molecular barcodes, optimized DNA extraction, and bioinformatics methods incorporating genomics context such as mutational signature and variant allelic frequency.

Transparent reporting is essential for the critical evaluation of studies. However, the reporting of statistical methods for studies in the biomedical sciences is often limited. This systematic review examines the quality of reporting for two statistical tests, t-tests and ANOVA, for papers published in a selection of physiology journals in June 2017. Of the 328 original research articles examined, 277 (84.5%) included an ANOVA or t-test or both. However, papers in our sample were routinely missing essential information about both types of tests: 213 papers (95% of the papers that used ANOVA) did not contain the information needed to determine what type of ANOVA was performed, and 26.7% of papers did not specify what post-hoc test was performed. Most papers also omitted the information needed to verify ANOVA results. Essential information about t-tests was also missing in many papers. We conclude by discussing measures that could be taken to improve the quality of reporting.

AbstractX 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 rl 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.

As a follow‐up to genome‐wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well‐known ovarian cancer risk factors and their interactions with 28 genome‐wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case–control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10−4). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46–0.60) compared to 0.71 (95%CI = 0.66–0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1–5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene–environment analysis in ovarian cancer.

Abstract Introduction: Mitochondria are intracellular organelles that play an important role in energy metabolism, reactive oxygen species generation, and apoptosis. They possess their own haploid genome distinct from nuclear DNA, and mutations in mitochondrial DNA (mtDNA) have previously been linked to cancer risk, disease progression, and treatment response for multiple tumor types. However, the spectrum of mtDNA variation in epithelial ovarian carcinoma (EOC) and associations with clinical characteristics have not been fully explored. Methods: We extracted mtDNA from fresh frozen EOC tumors of women diagnosed with primary invasive epithelial ovarian, fallopian tube, or peritoneal carcinoma, and performed targeted sequencing using the Illumina MiSeq Reporter mtDNA workflow. Homoplasmic and heteroplasmic SNVs and indels were called using an in-house bioinformatics pipeline. We estimated mtDNA haplogroup using Haplogrep and compared to results from available germline mtDNA SNP genotyping. We derived regional weighted variant burden scores based on mtDNA functional domains, accounting for underlying haplogroup, heteroplasmic fraction, and predicted functional impact. To assess potential differences in variant burden by histotype, we compared burden scores using Kruskal-Wallis testing. Among high-grade serous tubo-ovarian cases, we explored associations with overall survival (OS) using Cox proportional hazards regression. After quality control filtering, 389 mitochondrial genomes were analyzed, including those from 320 tubo-ovarian high-grade serous tumors, 42 endometrioid tumors, and 18 clear cell tumors. Results: Accounting for predicted haplogroup, we identified a median count of two global private heteroplasmies per tumor (range: 0-9). Elevated mutation rates were identified in the regulatory control region as well as among the various mt-tRNAs (&gt;3 per 10kb per sample) relative to the remainder of the mtDNA. The majority of recurrent variants across tumors were also identified in the mtDNA control region. We did not observe an association with overall mtDNA variant burden and patient age at diagnosis. Regional mtDNA burden scores revealed a significant association with histotype (P&lt;0.001) for the oxidative phosphorylation complex IV genes, with higher burden in tubo-ovarian high-grade serous and clear-cell tumors relative to endometrioid tumors. For patients with tubo-ovarian high-grade serous tumors, increased variant burden in MT-CO1 (cytochrome C oxidase I) was associated with improved OS (HR = 0.32; 95% CI = [0.14,0.70], LRT P = 0.001) after adjustment for age, stage, and debulking status. Conclusions: We demonstrate a wide range of tumor mtDNA variant patterns in the largest known collection of EOC patients. This work provides the first extensive evaluation of mtDNA variation in these diseases and highlights the potential prognostic relevance of mitochondrial dysfunction. Citation Format: Nicholas B. Larson, Stacey J. Winham, Chen Wang, Jared M. Evans, Sebastian M. Armasu, Bryan M. McCauley, Sean J. Yoder, Lan Min Zhang, Catherine M. Phelan, Ellen L. Goode. Targeted deep sequencing of mitochondrial DNA in epithelial ovarian carcinomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 420.

Abstract Purpose: Volumetric breast density and body mass index (BMI) are increasingly used for breast cancer risk stratification. We examine if the effect of volumetric breast density on breast cancer risk increases with increasing BMI. Methods: Participants were from two case-control studies within mammography practices, the San Francisco Mammography Registry and the Mayo Clinic Rochester, Breast Screening Practice. Breast cancers were included if diagnosed between January 2007 and 2015 and with a screening full-field digital mammogram at least 6 months prior to diagnosis; the earliest mammogram within 5 years of diagnosis was selected. Up to three controls were matched to each case on age, race, date, mammography machine, and state of residence. Volumetric percent density (VPD), dense volume (DV) and non-dense volume (NDV) were measured using VolparaTM software. Breast cancer risk was assessed using logistic regression stratified by menopause status and adjusted for matching factors, family history of breast cancer, parity/age at first birth, and postmenopausal hormone therapy. Models for DV were additionally adjusted for NDV, and NDV models for DV. Multiplicative interactions were fit between BMI categories [&lt;25 (underweight), 25-30 (normal weight) and &gt;30 (overweight/obese) kg/m2] and VPD, DV, and NDV, and trend tests assessed for increasing odds ratios (OR) with increasing BMI. Results: A total of 781 premenopausal and 1850 postmenopausal breast cancers and 3535 controls were included in the analysis. Median age of premenopausal women was 45 years (IQR: 6.1) and 63.3 years (IQR: 14) for postmenopausal women. Cases vs. controls had greater VPD and DV for both premenopausal (VPD: 14.9% vs. 12.0%, DV: 74.1 cm3 vs. 64.4 cm3) and postmenopausal women (VPD: 6.8% vs. 6.1%, DV: 53.4 cm3 vs. 48.0 cm3)(all p's&lt;0.001). Trends between increasing BMI and VPD were evident for both pre (ptrend=0.0007) and postmenopausal (ptrend=0.0005) women. Among premenopausal women, the odds ratio (OR) for breast cancer associated with a 10% increase in VPD was of 1.39, 2.19 and 2.88 for BMI &lt;25, 25-30 and &gt;30 kg/m2 (p-trend=0.0007), respectively. For DV, OR's were 1.39, 1.33 and 1.51 for a 1 SD increase in DV, respectively, though the interaction with DV was not significant (ptrend=0.68). Among postmenopausal women, a 10% increase in VPD was associated with OR's of 1.35, 2.03, 3.6 for BMI &lt;25, 25-29, &gt;30-kg/m2 (ptrend =0.0001), respectively, and 1.31, 1.34 and 1.65 for a 1 SD increase in DV (ptrend =0.01), respectively. Associations between NDV and breast cancer risk did not differ by BMI category for premenopausal (ptrend =0.52) or postmenopausal (ptrend =0.07) women. Conclusions: The effect of VPD on breast cancer risk is strongest in overweight/obese women. As volumetric breast density and BMI are commonly used in clinical risk stratification, these differences in risk have high clinical relevance for informing prevention decisions. Citation Format: Natalie J. Engmann, Christopher G. Scott, Matthew R. Jensen, Stacey J. Winham, Lin Ma, Kathleen R. Brandt, Amir Mahmoudzadeh, Serghei Malkov, Dana Whaley, Carrie Hruska, Fang Fang Wu, Diana L. Miglioretti, Aaron D. Norman, John Heine, John Shepherd, Vernon S. Pankratz, Celine M. Vachon, Karla Kerlikowske. Overweight and obese women with high volumetric breast density at high breast cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3226.

Abstract Tubo-ovarian high-grade serous cancer (HGSC) is the most common histotype of epithelial ovarian cancer (70%-80%) and the leading cause of death from gynecological malignancy. Novel monoclonal antibody and adoptive cell transfer immunotherapies have found success in cancer in the past five years; however success in HGSC treatment has been elusive thus far. We have recently shown that CD8+ T-lymphocytes infiltrating the tumor proper (TILs) associate in a dose response manner with longer patient survival time regardless of age, stage or residual disease. As understanding host factors associated with this clinically relevant tumor phenotype will contribute to the development of targeted therapies, we sought to identify germline genetic factors associated with CD8+ TIL levels in HGSC. We performed a germline genome-wide association study (GWAS) of CD8+ TILs using 1204 HGSC patients of European ancestry from twelve studies with centralized genotyping, immunohistochemistry, and harmonization of clinical data. Germline genotypes were assayed using the Illumina Infinium OncoArray Beadchip and were imputed to over 11.4 million variants based on the 1000 Genomes Project reference panel. CD8 immunohistochemistry was performed on tissue microarrays using the Leica Bond Rx stainer and scored into four levels of CD8+ TILs per x400 magnification of a 0.55-mm diameter field for each tumor hotspot: none (no CD8+ TILs), low (1-2 CD8+ TILs), moderate (3-19 CD8+ TILs), high (20 or more CD8+ TILs). GWAS used multinomial regression modelling of these four CD8+ TIL levels adjusted for study site, age, and European ancestry principal components. On the average, CD8+ TIL levels were 18%, 19%, 44% and 18% for none, low, moderate and high levels, respectively. As expected, CD8+ TIL levels associated with overall survival time (p &lt; 2.6x10-8), and CD8+ TIL levels and age at diagnosis were similarly distributed across study sites. We identified a region on chromosome 22 with multiple, correlated variants (MAF &gt; 20%) associated with CD8+ TIL levels at p-value &lt; 2.0x10-5. The most statistically significant variant was an indel (rs557925408) located in the an intron of the MYO18B gene (Myosin XVIIIB, p-value &lt; 8.7x10-6). Each copy of this variant associated with a reduced extent of CD8+ TIL infiltration (OR low v none = 0.92; OR moderate v none = 0.61). An independent set of 1200 OncoArray genotyped HGSC cases are undergoing CD8+ TIL scoring in December 2017 and will be included in final genome-wide analysis. This will provide improved statistical power and enable consideration of additional genetic models (e.g., linear trend tests). Comprehensive integrative analysis of germline, tumor, and clinical features provides a model for clinical molecular epidemiology studies and will be key to increasing our understanding of this important HGSC immunophenotype. Citation Format: Yanina Natanzon, Martin Köbel, Stacey J. Winham, Sebastian M. Armasu, Bryan M. McCauley, Robert A. Vierkant, Julie M. Cunningham, David Bowtell, Ian G. Campbell, Jenny Chang-Claude, Anna deFazio, Peter A. Fasching, Mark T. Goodman, Beth Y. Karlan, Francesmary Modugno, Kirsten B. Moysich, Roberta B. Ness, Weiva Sieh, Paul D. Pharoah, Susan J. Ramus, Ellen L. Goode. Tumor-infiltrating CD8-positive T-lymphocytes in tubo-ovarian high-grade serous cancer are associated with multiple germline variants in 22q12.1 in a genome-wide association analysis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 225.