Citations

Abstract Massively parallel next generation sequencing data (NGS) has proved to be useful in identification of rare Mendelian disorders. Efforts are on-going to use NGS for discovery of novel rare variants in common disease etiology. Distinguishing rare private mutations from causal variants remains a major challenge in complex disorders. With decreasing cost, investigators are generating NGS data for common cancers. However, power and other considerations suggest that a consortia approach of pooling data is more likely to succeed where individual efforts fail. Combining data within such consortia brings its own level of challenges. We enumerate these challenges and suggest recommendations based on our experience of combining and harmonizing NGS data on 358 breast cancer samples from four centers; (City of Hope n=8, Mayo Clinic n=221, MSKCC n=96 and University of Pennsylvania n=33). Major concerns were ability of individual datasets to discover known and novel variations, depth of coverage across the exome, and quality of variants. Using pre- and post-processed data, we demonstrate an increased power for detecting variants when using a joint-calling method, quality control (QC) and quality assurance pre- and post-merge, as well as the advantages of combined annotations and filtering procedures. To generate a filtered list of high confident calls, we subset the data on Depth per Allele (DP), Genotyping Quality (GQ), Allele balance, Phred likelihood score (PL) as well as using frequency filters on the alternate allele both within the data and comparing against convenience control summaries. We describe general principles and QC measures to be adopted when NGS data is used from diverse sources, specifically applicable to large consortia. Note: This abstract was not presented at the meeting. Citation Format: Joseph Vijai, Steven Hart, Tinu Thomas, Bradley Wubbenhorst, Lucia Guidugli, Kasmintan Schrader, Kara Maxwell, Lauren Jacobs, Danylo Villano, Robert Klein, Steven Lipkin, Susan Neuhausen, Jeffrey Weitzel, David Altshuler, Fergus Couch, Kenneth Offit, Katherine Nathanson. Harmonization of next generation sequencing data within consortia for gene discovery in familial breast 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 2378. doi:10.1158/1538-7445.AM2014-2378

Abstract BRCA1 and BRCA2 (BRCA) testing is uninformative in approximately 80% of families with clinical features of inherited susceptibility to breast cancer. Early identification of the individuals at risk in BRCAx families may lead to enhanced screening and prevention strategies, and potentially improved overall survival as has been seen for carriers of BRCA1/2 mutations. A proportion of these families are likely explained by mutations in other high penetrance and moderate penetrance breast cancer susceptibility genes. Although multi-gene panels are now commercially available for clinical testing, only limited information is available on the spectrum of mutations in breast cancer susceptibility genes in specific subgroups of high-risk breast cancer patients. We performed whole exome sequencing in 278 BRCA1/2 negative individuals with high-risk familial breast cancer, defined as a proband and at least two first or second degree relatives with breast cancer. Samples included 232 individuals from 97 families and 49 unrelated individuals, for a total of 146 independent cases. Data were analyzed for identification of all variants in 18 known high and moderate penetrance breast cancer susceptibility genes, in addition to BRCA1 and BRCA2. Two of the 97 families and one of the 49 singletons (2% of independent cases) were found to have deleterious mutations in high penetrance genes, namely TP53 (2) and CDH1 (1). In addition, nine of 146 independent cases (6%) were found to have deleterious mutations in the other predisposition genes. Mutations in ATM and CHEK2 accounted for seven of the ten moderate penetrance mutations identified; additionally one family had a mutation in BLM, one in MRE11A, and one in PALB2. Furthermore, 15 of 146 independent cases (10%) were found to have variants of unknown significance (VUS) in one of eight genes. There were no deleterious mutations or VUS in eight of the 18 genes. Our data show that significant genetic heterogeneity exists in BRCAx families. Large-scale collaborative efforts will be required to attain sufficient power to understand how to appropriately apply these results clinically in cancer risk evaluation. Citation Format: Kara N. Maxwell, Lucia Guidugli, Kasmintan Schrader, Steven Hart, Vijai Joseph, Tinu Thomas, Xianshu Wang, Bradley Wubbenhorst, Robert Klein, Susan M. Domchek, Csilla Szabo, Susan Neuhausen, Jeffrey Weitzel, Katherine L. Nathanson, Kenneth Offit, Fergus Couch. High and moderate penetrance germline mutations in a number of genes are responsible for a small proportion of familial breast cancer risk in BRCAx families. [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 1291. doi:10.1158/1538-7445.AM2014-1291

Abstract INTRODUCTION Patient derived xenografts (PDX) may better reflect individual patient (pt) tumor biology; however, the feasibility of collecting PDX from percutaneous tumor biopsies (PTB) in the neoadjuvant setting is unknown. Furthermore, drug response phenotypes observed in PDX have not been prospectively compared to the corresponding pt clinical outcomes. METHODS The Breast Cancer Genome Guided Therapy Study (BEAUTY) is a prospective Mayo study of pts with high-risk breast cancer treated with neoadjuvant weekly paclitaxel (T) +/- trastuzumab followed by anthracycline based chemotherapy. PTB (at baseline) and residual surgical tissue (after all chemotherapy) are obtained for next generation sequencing (NGS) and PDX. Tumor biopsies (1-2 cores from 14 gauge needle) were implanted with Matrigel <1 hour of collection in the flanks of NOD-SCID or NSG mice. Low dose estradiol was supplemented in the drinking water. Primary outgrowth rate was defined as PDX tumor volume >50 mm3. Take rate was defined as development of at least 1 stably transplantable xenograft line/pt. To determine whether clinical T response assessed by MRI corresponded with in vivo T response, pretreatment PDX from 5 pts were injected into NOD-SCID mice (20 mice per pt PDX) and when tumors reached 100-200mm3, mice were randomized to no treatment vs T (20 mg/kg, ip. every 3-4 days). Two of these 5 patients had a MRI response defined as >30% decrease in longest lesion. RESULTS Pretreatment PTB from 81 unique pts were implanted in 251 mice (2-4 mice/pt). PDX outgrowth rates were 33.3% (27/81 pts) and 22 stable PDX were established (overall take rate 27.2%). Take rates were as follows: triple negative breast cancer (46%; 13/28); HER2 (27%; 6/22), Luminal B (13%; 3/22), and luminal A (0%; 0/9). Residual surgical tumor (after all treatment) from 17 pts was injected into 85 mice (average 5 mice/pt) and the initial outgrowth rate was 23% (4/17) with 3 stably transplantable lines established. PDX, derived from pretreatment PTB of 5 pts (2 responders and 3 non-responders), were assessed for in vivo T response. The size of the T treated group was significantly smaller than the no treatment group for the PDX derived from the 2 clinical responders, with complete disappearance of tumor by 18 days. In contrast, the PDX derived from the 3 clinical non-responders had no evidence for T response. CONCLUSIONS We have demonstrated the feasibility of using PTB to establish PDX in a prospective neoadjuvant clinical study and have demonstrated similar T drug response phenotypes in in the PDX as seen in the corresponding pt. These data suggest that PDX generated prospectively may be useful for biomarker validation and the development and individualization of new drug therapy. Funded by the Mayo Clinic Center for Individualized Medicine and the Mayo Clinic Cancer Center Citation Format: Jia Yu, Ping Yin, Bowen Gao, Jason P. Sinnwell, Ann M. Moyer, Daniel W. Visscher, Amy L. Conners, Travis J. Dockter, Krishna R. Kalari, Xiaojia Tang, Kevin J. Thompson, Hugues Sicotte, Douglas W. Mahoney, Steven N. Hart, Peter T. Vedell, Poulami Barman, Katie N. Jones, Sarah A. McLaughlin, John A. Copland, Alvaro Moreno Aspitia, Donald W. Northfelt, Richard J. Gray, Vera J. Suman, Jeanette E. Eckel Passow, Eric D. Wieben, James N. Ingle, Zhenkun Lou, Gianrico Farrugia, Richard Weinshilboum, Matthew P. Goetz, Judy C. Boughey, Liewei Wang. Feasibility of using percutaneous tumor biopsies from a prospective neoadjuvant breast cancer study to develop patient derived xenografts and assess in vivo chemotherapy sensitivity. [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 1195. doi:10.1158/1538-7445.AM2014-1195

Purpose Recent advances in DNA sequencing have led to the development of breast cancer susceptibility gene panels for germline genetic testing of patients. We assessed the frequency of mutations in 17 predisposition genes, including BRCA1 and BRCA2, in a large cohort of patients with triple-negative breast cancer (TNBC) unselected for family history of breast or ovarian cancer to determine the utility of germline genetic testing for those with TNBC. Patients and Methods Patients with TNBC (N = 1,824) unselected for family history of breast or ovarian cancer were recruited through 12 studies, and germline DNA was sequenced to identify mutations. Results Deleterious mutations were identified in 14.6% of all patients. Of these, 11.2% had mutations in the BRCA1 (8.5%) and BRCA2 (2.7%) genes. Deleterious mutations in 15 other predisposition genes were detected in 3.7% of patients, with the majority observed in genes involved in homologous recombination, including PALB2 (1.2%) and BARD1, RAD51D, RAD51C, and BRIP1 (0.3% to 0.5%). Patients with TNBC with mutations were diagnosed at an earlier age (P < .001) and had higher-grade tumors (P = .01) than those without mutations. Conclusion Deleterious mutations in predisposition genes are present at high frequency in patients with TNBC unselected for family history of cancer. Mutation prevalence estimates suggest that patients with TNBC, regardless of age at diagnosis or family history of cancer, should be considered for germline genetic testing of BRCA1 and BRCA2. Although mutations in other predisposition genes are observed among patients with TNBC, better cancer risk estimates are needed before these mutations are used for clinical risk assessment in relatives.

Abstract BACKGROUND Next generation sequencing (NGS)-based assays continue to redefine the field of genetic testing. Owing to the complexity of the data, bioinformatics has become a necessary component in any laboratory implementing a clinical NGS test. CONTENT The computational components of an NGS-based work flow can be conceptualized as primary, secondary, and tertiary analytics. Each of these components addresses a necessary step in the transformation of raw data into clinically actionable knowledge. Understanding the basic concepts of these analysis steps is important in assessing and addressing the informatics needs of a molecular diagnostics laboratory. Equally critical is a familiarity with the regulatory requirements addressing the bioinformatics analyses. These and other topics are covered in this review article. SUMMARY Bioinformatics has become an important component in clinical laboratories generating, analyzing, maintaining, and interpreting data from molecular genetics testing. Given the rapid adoption of NGS-based clinical testing, service providers must develop informatics work flows that adhere to the rigor of clinical laboratory standards, yet are flexible to changes as the chemistry and software for analyzing sequencing data mature.

Significance The major portion of hereditary breast cancer still remains unexplained, and many susceptibility loci are yet to be found. Exome sequencing of 24 high-risk familial BRCA1/2 -negative breast cancer patients and further genotyping of a large sample set of breast/ovarian cancer cases and controls was used to discover previously unidentified susceptibility alleles and genes. A significant association of a FANCM nonsense mutation with breast cancer, especially triple-negative breast cancer, identifies FANCM as a breast cancer susceptibility gene. Identification of such risk alleles is expected to improve cancer risk assessment for breast cancer patients and families, and may lead to improvements in the prevention, early diagnosis, and treatment of cancer.

Abstract Motivation: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. Results: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. Availability and implementation: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/ , and R package at https://github.com/topsoil/patternCNV/ . Contact: Asmann.Yan@mayo.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract Motivation: The Biological Reference Repository (BioR) is a toolkit for annotating variants. BioR stores public and user-specific annotation sources in indexed JSON-encoded flat files (catalogs). The BioR toolkit provides the functionality to combine and retrieve annotation from these catalogs via the command-line interface. Several catalogs from commonly used annotation sources and instructions for creating user-specific catalogs are provided. Commands from the toolkit can be combined with other UNIX commands for advanced annotation processing. We also provide instructions for the development of custom annotation pipelines. Availability and implementation: The package is implemented in Java and makes use of external tools written in Java and Perl. The toolkit can be executed on Mac OS X 10.5 and above or any Linux distribution. The BioR application, quickstart, and user guide documents and many biological examples are available at http://bioinformaticstools.mayo.edu . Contact: Kocher.JeanPierre@mayo.edu Supplementary information: Supplementary data are available at Bioinformatics online.

Abstract Ovarian cancer is a clinically and molecularly heterogeneous disease. The driving forces behind this variability are unknown. Here, we report wide variation in the expression of the DNA cytosine deaminase APOBEC3B, with elevated expression in the majority of ovarian cancer cell lines (three SDs above the mean of normal ovarian surface epithelial cells) and high-grade primary ovarian cancers. APOBEC3B is active in the nucleus of several ovarian cancer cell lines and elicits a biochemical preference for deamination of cytosines in 5′-TC dinucleotides. Importantly, examination of whole-genome sequence from 16 ovarian cancers reveals that APOBEC3B expression correlates with total mutation load as well as elevated levels of transversion mutations. In particular, high APOBEC3B expression correlates with C-to-A and C-to-G transversion mutations within 5′-TC dinucleotide motifs in early-stage high-grade serous ovarian cancer genomes, suggesting that APOBEC3B-catalyzed genomic uracil lesions are further processed by downstream DNA "repair" enzymes including error-prone translesion polymerases. These data identify a potential role for APOBEC3B in serous ovarian cancer genomic instability. Cancer Res; 73(24); 7222–31. ©2013 AACR.

Abstract Background: Given the high technical reproducibility and orders of magnitude greater resolution than microarrays, next-generation sequencing of mRNA (RNA-Seq) is quickly becoming the de facto standard for measuring levels of gene expression in biological experiments. Two important questions must be taken into consideration when designing a particular experiment, namely, 1) how deep does one need to sequence? and, 2) how many biological replicates are necessary to observe a significant change in expression? Results: Based on the gene expression distributions from 127 RNA-Seq experiments, we find evidence that 91% ± 4% of all annotated genes are sequenced at a frequency of 0.1 times per million bases mapped, regardless of sample source. Based on this observation, and combining this information with other parameters such as biological variation and technical variation that we empirically estimate from our large datasets, we developed a model to estimate the statistical power needed to identify differentially expressed genes from RNA-Seq experiments. Conclusions: Our results provide a needed reference for ensuring RNA-Seq gene expression studies are conducted with the optimally sample size, power, and sequencing depth. We also make available both R code and an Excel worksheet for investigators to calculate for their own experiments.