Search

Abstract Background Atlantic salmon (Salmo salar L.) is an environmentally and economically important organism and its gene content is reasonably well characterized. From a transcriptional standpoint, it is important to characterize the changes in gene expression over the course of unperturbed early development, from fertilization through to the parr stage. Findings S. salar samples were taken at 17 time points from 2 to 89 days post fertilization. Total RNA was extracted and cRNA was synthesized and hybridized to a newly developed 44K oligo salmonid microarray platform. Quantified results were subjected to preliminary data analysis and submitted to NCBI's Gene Expression Omnibus (GEO). Data can be found under the GEO accession number GSE25938. http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25938 Conclusions Throughout the entire period of development, several thousand genes were found to be differentially regulated. This work represents the trancriptional characterization of a very large geneset that will be extremely valuable in further examination of the transcriptional changes in Atlantic salmon during the first few months of development. The expression profiles can help to annotate salmon genes in addition to being used as references against any number of experimental variables to which developing salmonids might be subjected.

Abstract Background The volume of viral genomic sequence data continues to increase rapidly. This is especially true for the smaller RNA viruses, which are relatively easy to sequence in large numbers. The data volumes cause a number of significant problems for research applications that require large multiple alignments of essentially complete genomes, which are of the order of 10 kb. Findings We present a simple strategy to enable the creation of large quasi-multiple sequence alignments from pairwise alignment data. This process is suitable for large, closely related sequences such as the polyproteins of dengue viruses, which need the insertion of very few indels. Conclusion The quasi-multiple sequence alignments generated by KISSa are sufficiently accurate to support tree-based genome selection for interactive bioinformatics analysis tools. The speed of this process is critical to providing an interactive experience for the user.

Abstract Background Many sequence-searching tools have limiting factors for their use. For example, they may be platform specific, enforce restrictive size limits and sequences to be searched, or only allow searches of one of DNA or protein. Findings We present an easy-to-use, fast, platform-independent tool to search for amino acid or nucleotide patterns within one or many protein or nucleic acid sequences. The user can choose to search for regular expressions or perform a fuzzy search in which a particular number of errors is accepted during matching of a sequence. Positions of mismatches in fuzzy searches are displayed graphically the user. Conclusion SeqS provides an improved feature set and functions as a stand-alone tool or could be integrated into other bioinformatics platforms.

This proof-of-concept evaluation demonstrates that next-generation sequencing-based liquid biopsy can detect genomic alterations in the blood of cats with cancer and the absence of such alterations in the blood of presumably cancer-free cats. Two cats with cytologically confirmed lymphoma and nine presumably cancer-free cats were included in this analysis. Whole blood was collected from each subject and samples were subjected to DNA extraction, library preparation, and next-generation sequencing. Both cancer-diagnosed subjects had somatic copy number variants (a "cancer signal") identified in cell-free DNA, suggesting the current presence of cancer in these subjects. All nine presumably cancer-free subjects had unremarkable genomic profiles, suggesting the absence of cancer in these subjects. Liquid biopsy using next-generation sequencing of cell-free DNA allows for blood-based detection of cancer-associated genomic alterations in cats. Such technology has the potential to offer considerable utility in veterinary medicine, particularly for the non-invasive prioritization of small cell intestinal lymphoma versus inflammatory bowel disease in cats with gastrointestinal signs. This study lays the foundation for future studies to fully validate this type of testing for use in clinical practice., Competing Interests: CR-P, PN, MM, FM, TT, LM, MK, BF, SG, JR, DG, SH, JT, DT, AF, KK were employed by PetDx. This study received funding from PetDx. The funder had the following involvement with the study: study design, data collection and analysis, decision to publish, and preparation of the manuscript., (Copyright © 2024 Ruiz-Perez, Nakashe, Marshall, Marass, Tang, McLennan, Kroll, Flesner, Gray, Rafalko, Grosu, Hicks, Tynan, Tsui, Flory and Kruglyak.)

Background: Liquid biopsy is a minimally-invasive method of sampling bodily fluids, capable of revealing evidence of cancer. The distribution of cell-free DNA (cfDNA) fragment lengths has been shown to differ between healthy subjects and cancer patients, whereby the distributional shift correlates with the sample's tumour content. These fragmentomic data have not yet been utilised to directly quantify the proportion of tumour-derived cfDNA in a liquid biopsy., Results: We used statistical learning to predict tumour content from Fourier and wavelet transforms of cfDNA length distributions in samples from 118 cancer patients. The model was validated on an independent dilution series of patient plasma., Conclusions: This proof of concept suggests that our fragmentomic methodology could be useful for predicting tumour content in liquid biopsies., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Age-related somatic genomic alterations in hematopoietic cell lines have been well characterized in humans; however, this phenomenon has not been well studied in other species. Next-generation sequencing-based liquid biopsy testing for cancer detection was recently developed for dogs and has been used to study the genomic profiles of blood samples from thousands of canine patients since 2021. In this study, 4870 client-owned dogs with and without a diagnosis or suspicion of cancer underwent liquid biopsy testing by this method. Copy number variants detected exclusively in genomic DNA derived from white blood cells (WBC gDNA-specific CNVs) were observed in 126 dogs (2.6%; 95% CI: 2.2-3.1); these copy number variants were absent from matched plasma cell-free DNA, and from tumor tissue in dogs with concurrent cancer. These findings were more common in older dogs and were persistent in WBC gDNA in over 70% of patients, with little to no change in the amplitude of the signal across longitudinal samples. Many of these alterations were observed at recurrent locations in the genome across subjects; the most common finding was a partial loss on CFA25, typically accompanied by a partial gain on the same chromosome. These early findings suggest that age-related somatic alterations may be present at an appreciable frequency in the general canine population. Further research is needed to determine the clinical significance of these findings.

Objective: To review ordering patterns, positivity rates, and outcome data for a subset of consecutive samples submitted for a commercially available, blood-based multicancer early-detection liquid biopsy test for dogs using next-generation sequencing at 1 laboratory., Sample: 1,500 consecutively submitted blood samples from client-owned dogs with and without clinical suspicion and/or history of cancer for prospective liquid biopsy testing between December 28, 2021, and June 28, 2022., Procedures: We performed a retrospective observational study, reviewing data from 1,500 consecutive clinical samples submitted for liquid biopsy testing. Outcome data were obtained via medical record review, direct communication with the referring clinic, and/or a patient outcome survey through October 16, 2022., Results: Sixty-four percent (910/1,419) of reportable samples were submitted for cancer screening, 26% (366/1,419) for aid in diagnosis, and 10% (143/1,419) for other indications. The positivity rate was 25.4% (93/366) in aid-in-diagnosis patients and 4.5% (41/910) in screening patients. Outcome data were available for 33% (465/1,401) of patients, and outcomes were classifiable for 428 patients. The relative observed sensitivity was 61.5% (67/109) and specificity was 97.5% (311/319). The positive predictive value was 75.0% (21/28) for screening patients and 97.7% (43/44) for aid-in-diagnosis patients, and the time to diagnostic resolution following a positive result was < 2 weeks in most cases., Clinical Relevance: Liquid biopsy using next-generation sequencing represents a novel tool for noninvasive detection of cancer in dogs. Real-world clinical performance meets or exceeds expectations established in the test's clinical validation study.

Cancer is the leading cause of death in dogs, yet there are no established screening paradigms for early detection. Liquid biopsy methods that interrogate cancer-derived genomic alterations in cell-free DNA in blood are being adopted for multi-cancer early detection in human medicine and are now available for veterinary use. The CANcer Detection in Dogs (CANDiD) study is an international, multi-center clinical study designed to validate the performance of a novel multi-cancer early detection "liquid biopsy" test developed for noninvasive detection and characterization of cancer in dogs using next-generation sequencing (NGS) of blood-derived DNA; study results are reported here. In total, 1,358 cancer-diagnosed and presumably cancer-free dogs were enrolled in the study, representing the range of breeds, weights, ages, and cancer types seen in routine clinical practice; 1,100 subjects met inclusion criteria for analysis and were used in the validation of the test. Overall, the liquid biopsy test demonstrated a 54.7% (95% CI: 49.3-60.0%) sensitivity and a 98.5% (95% CI: 97.0-99.3%) specificity. For three of the most aggressive canine cancers (lymphoma, hemangiosarcoma, osteosarcoma), the detection rate was 85.4% (95% CI: 78.4-90.9%); and for eight of the most common canine cancers (lymphoma, hemangiosarcoma, osteosarcoma, soft tissue sarcoma, mast cell tumor, mammary gland carcinoma, anal sac adenocarcinoma, malignant melanoma), the detection rate was 61.9% (95% CI: 55.3-68.1%). The test detected cancer signal in patients representing 30 distinct cancer types and provided a Cancer Signal Origin prediction for a subset of patients with hematological malignancies. Furthermore, the test accurately detected cancer signal in four presumably cancer-free subjects before the onset of clinical signs, further supporting the utility of liquid biopsy as an early detection test. Taken together, these findings demonstrate that NGS-based liquid biopsy can offer a novel option for noninvasive multi-cancer detection in dogs., Competing Interests: AF, KMK, JAT, LMM, JMR, PCF, GEH, FM, PN, CAR-P, DMF, TJ, RM-P, KW, ALM-W, DSG, AP, KM, SCH, KML, LH, JC, IC, and DWYT are employed by or affiliated with PetDx; and receive compensation from PetDx and/or hold vested or unvested equity in PetDx. Drs. Flesner and Leibman are members of the PetDx clinical advisory board and receive compensation and equity. Dr. Borgatti has ownership interest (including patents) in a patent entitled “Reduction of EGFR therapeutic toxicity” filed by the University of Minnesota Office of Technology Commercialization. Drs. Cahalane, Phillips, and Aiken are investors in PetDx. Drs. Jensen, van den Boom, Dharajiya, and Diaz are advisors for PetDx and hold vested or unvested equity in PetDx. Drs. Chorny, Kruglyak, Grosu, Marass, Ruiz-Perez and Tsui hold pending patent applications related to technology described in this work. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Background: The analysis of liquid biopsies, e.g., circulating tumor cells (CTCs) is an appealing diagnostic concept for targeted therapy selection. In this proof-of-concept study, we aimed to perform multiparametric analyses of CTCs to select targeted therapies for metastatic breast cancer patients., Methods: First, CTCs of five metastatic breast cancer patients were analyzed by whole exome sequencing (WES). Based on the results, one patient was selected and monitored by longitudinal and multiparametric liquid biopsy analyses over more than three years, including WES, RNA profiling, and in vitro drug testing of CTCs., Results: Mutations addressable by targeted therapies were detected in all patients, including mutations that were not detected in biopsies of the primary tumor. For the index patient, the clonal evolution of the tumor cells was retraced and resistance mechanisms were identified. The AKT1 E17K mutation was uncovered as the driver of the metastatic process. Drug testing on the patient's CTCs confirmed the efficacy of drugs targeting the AKT1 pathway. During a targeted therapy chosen based on the CTC characterization and including the mTOR inhibitor everolimus, CTC numbers dropped by 97.3% and the disease remained stable as determined by computer tomography/magnetic resonance imaging., Conclusion: These results illustrate the strength of a multiparametric CTC analysis to choose and validate targeted therapies to optimize cancer treatment in the future. Furthermore, from a scientific point of view, such studies promote the understanding of the biology of CTCs during different treatment regimens.

Motivation: The high resolution of single-cell DNA sequencing (scDNA-seq) offers great potential to resolve intratumor heterogeneity (ITH) by distinguishing clonal populations based on their mutation profiles. However, the increasing size of scDNA-seq datasets and technical limitations, such as high error rates and a large proportion of missing values, complicate this task and limit the applicability of existing methods., Results: Here, we introduce BnpC, a novel non-parametric method to cluster individual cells into clones and infer their genotypes based on their noisy mutation profiles. We benchmarked our method comprehensively against state-of-the-art methods on simulated data using various data sizes, and applied it to three cancer scDNA-seq datasets. On simulated data, BnpC compared favorably against current methods in terms of accuracy, runtime and scalability. Its inferred genotypes were the most accurate, especially on highly heterogeneous data, and it was the only method able to run and produce results on datasets with 5000 cells. On tumor scDNA-seq data, BnpC was able to identify clonal populations missed by the original cluster analysis but supported by Supplementary Experimental Data. With ever growing scDNA-seq datasets, scalable and accurate methods such as BnpC will become increasingly relevant, not only to resolve ITH but also as a preprocessing step to reduce data size., Availability and Implementation: BnpC is freely available under MIT license at https://github.com/cbg-ethz/BnpC., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press.)

Purpose: The enucleation rate for retinoblastoma has dropped from over 95% to under 10% in the past 10 years as a result of improvements in therapy. This reduces access to tumor tissue for molecular profiling, especially in unilateral retinoblastoma, and hinders the confirmation of somatic RB1 mutations necessary for genetic counseling. Plasma cell-free DNA (cfDNA) has provided a platform for noninvasive molecular profiling in cancer, but its applicability in low tumor burden retinoblastoma has not been shown. We analyzed cfDNA collected from 10 patients with available tumor tissue to determine whether sufficient tumorderived cfDNA is shed in plasma from retinoblastoma tumors to enable noninvasive RB1 mutation detection., Methods: Tumor tissue was collected from eye enucleations in 10 patients diagnosed with advanced intra-ocular unilateral retinoblastoma, three of which went on to develop metastatic disease. Tumor RB1 mutation status was determined using an FDA-cleared tumor sequencing assay, MSK-IMPACT. Plasma samples were collected before eye enucleation and analyzed with a customized panel targeting all exons of RB1., Results: Tumor-guided genotyping detected 10 of the 13 expected somatic RB1 mutations in plasma cfDNA in 8 of 10 patients (average variant allele frequency 3.78%). Without referring to RB1 status in the tumor, de novo mutation calling identified 7 of the 13 expected RB1 mutations (in 6 of 10 patients) with high confidence., Conclusion: Plasma cfDNA can detect somatic RB1 mutations in patients with unilateral retinoblastoma. Since intraocular biopsies are avoided in these patients because of concern about spreading tumor, cfDNA can potentially offer a noninvasive platform to guide clinical decisions about treatment, follow-up schemes, and risk of metastasis., (© 2020 The Authors. Cancer Medicine published by John Wiley & Sons Ltd.)

High-grade serous ovarian carcinoma (HGSOC) remains the deadliest form of epithelial ovarian cancer and despite major efforts little improvement in overall survival has been achieved. Identification of recurring "driver" genetic lesions has the potential to enable design of novel therapies for cancer. Here, we report on a study to find such new therapeutic targets for HGSOC using exome-capture sequencing approach targeting all kinase genes in 127 patient samples. Consistent with previous reports, the most frequently mutated gene was TP53 (97% mutation frequency) followed by BRCA1 (10% mutation frequency). The average mutation frequency of the kinase genes mutated from our panel was 1.5%. Intriguingly, after BRCA1, JAK3 was the most frequently mutated gene (4% mutation frequency). We tested the transforming properties of JAK3 mutants using the Ba/F3 cell-based in vitro functional assay and identified a novel gain-of-function mutation in the kinase domain of JAK3 (p.T1022I). Importantly, p.T1022I JAK3 mutants displayed higher sensitivity to the JAK3-selective inhibitor Tofacitinib compared to controls. For independent validation, we re-sequenced the entire JAK3 coding sequence using tagged amplicon sequencing (TAm-Seq) in 463 HGSOCs resulting in an overall somatic mutation frequency of 1%. TAm-Seq screening of CDK12 in the same population revealed a 7% mutation frequency. Our data confirms that the frequency of mutations in kinase genes in HGSOC is low and provides accurate estimates for the frequency of JAK3 and CDK12 mutations in a large well characterized cohort. Although p.T1022I JAK3 mutations are rare, our functional validation shows that if detected they should be considered as potentially actionable for therapy. The observation of CDK12 mutations in 7% of HGSOC cases provides a strong rationale for routine somatic testing, although more functional and clinical characterization is required to understand which nonsynonymous mutations alterations are associated with homologous recombination deficiency., Competing Interests: The authors have declared that no competing interests exist.

Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of >10 5 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to parts per million. This resulted in median area under the curve (AUC) values of 0.98 in advanced cancers and 0.80 in early-stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that INVAR may be applied without requiring personalized sequencing panels so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Next-generation sequencing of DNA and RNA obtained from liquid biopsies of cancer patients may reveal important insights into disease progression and metastasis formation, and it holds the promise to enable new methods for noninvasive screening and clinical decision support. However, implementing liquid biopsy sequencing protocols is challenged by capturing circulating tumor cells or cell-free tumor DNA from blood samples, by amplifying genomic DNA and RNA in a reliable and unbiased manner, and by extracting biologically meaningful signals from the noisy sequencing data. In this chapter, we discuss computational methods for the analysis of DNA and RNA sequencing data obtained from liquid biopsies, addressing these challenges.

The factors responsible for the low detection rate of cell-free tumor DNA (ctDNA) in the plasma of patients with glioblastoma (GBM) are currently unknown. In this study, we measured circulating nucleic acids in patient-derived orthotopically implanted xenograft (PDOX) models of GBM ( n = 64) and show that tumor size and cell proliferation, but not the integrity of the blood-brain barrier or cell death, affect the release of ctDNA in treatment-naïve GBM PDOX. Analysis of fragment length profiles by shallow genome-wide sequencing (<0.2× coverage) of host (rat) and tumor (human) circulating DNA identified a peak at 145 bp in the human DNA fragments, indicating a difference in the origin or processing of the ctDNA. The concentration of ctDNA correlated with cell death only after treatment with temozolomide and radiotherapy. Digital PCR detection of plasma tumor mitochondrial DNA (tmtDNA), an alternative to detection of nuclear ctDNA, improved plasma DNA detection rate (82% vs. 24%) and allowed detection in cerebrospinal fluid and urine. Mitochondrial mutations are prevalent across all cancers and can be detected with high sensitivity, at low cost, and without prior knowledge of tumor mutations via capture-panel sequencing. Coupled with the observation that mitochondrial copy number increases in glioma, these data suggest analyzing tmtDNA as a more sensitive method to detect and monitor tumor burden in cancer, specifically in GBM, where current methods have largely failed. SIGNIFICANCE: These findings show that detection of tumor mitochondrial DNA is more sensitive than circulating tumor DNA analysis to detect and monitor tumor burden in patient-derived orthotopic xenografts of glioblastoma., (©2018 American Association for Cancer Research.)

Glioma is difficult to detect or characterize using current liquid biopsy approaches. Detection of cell-free tumor DNA (cftDNA) in cerebrospinal fluid (CSF) has been proposed as an alternative to detection in plasma. We used shallow whole-genome sequencing (sWGS, at a coverage of < 0.4×) of cell-free DNA from the CSF of 13 patients with primary glioma to determine somatic copy number alterations and DNA fragmentation patterns. This allowed us to determine the presence of cftDNA in CSF without any prior knowledge of point mutations present in the tumor. We also showed that the fragmentation pattern of cell-free DNA in CSF is different from that in plasma. This low-cost screening method provides information on the tumor genome and can be used to target those patients with high levels of cftDNA for further larger-scale sequencing, such as by whole-exome and whole-genome sequencing., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Existing methods to improve detection of circulating tumor DNA (ctDNA) have focused on genomic alterations but have rarely considered the biological properties of plasma cell-free DNA (cfDNA). We hypothesized that differences in fragment lengths of circulating DNA could be exploited to enhance sensitivity for detecting the presence of ctDNA and for noninvasive genomic analysis of cancer. We surveyed ctDNA fragment sizes in 344 plasma samples from 200 patients with cancer using low-pass whole-genome sequencing (0.4×). To establish the size distribution of mutant ctDNA, tumor-guided personalized deep sequencing was performed in 19 patients. We detected enrichment of ctDNA in fragment sizes between 90 and 150 bp and developed methods for in vitro and in silico size selection of these fragments. Selecting fragments between 90 and 150 bp improved detection of tumor DNA, with more than twofold median enrichment in >95% of cases and more than fourfold enrichment in >10% of cases. Analysis of size-selected cfDNA identified clinically actionable mutations and copy number alterations that were otherwise not detected. Identification of plasma samples from patients with advanced cancer was improved by predictive models integrating fragment length and copy number analysis of cfDNA, with area under the curve (AUC) >0.99 compared to AUC <0.80 without fragmentation features. Increased identification of cfDNA from patients with glioma, renal, and pancreatic cancer was achieved with AUC > 0.91 compared to AUC < 0.5 without fragmentation features. Fragment size analysis and selective sequencing of specific fragment sizes can boost ctDNA detection and could complement or provide an alternative to deeper sequencing of cfDNA., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Circulating tumor DNA (ctDNA) offers new opportunities for noninvasive cancer management. Detecting ctDNA in plasma is challenging because it constitutes only a minor fraction of the total cell-free DNA (cfDNA). Pre-analytical factors affect cfDNA levels contributed from leukocyte lysis, hence the ability to detect low-frequency mutant alleles. This study investigates the effects of the delay in processing, storage temperatures, different blood collection tubes, centrifugation protocols, and sample shipment on cfDNA levels. Peripheral blood (n = 231) from cancer patients (n = 62) were collected into K 3 EDTA or Cell-free DNA BCT tubes and analyzed by digital PCR, targeted amplicon, or shallow whole-genome sequencing. To assess pre-analytic effects, plasma was processed under different conditions after 0, 6, 24, 48, 96 hours, and 1 week at room temperature or 4°C, or using different centrifugation protocols. Digital PCR showed that cfDNA levels increased gradually with time in K 3 EDTA tubes, but were stable in BCT tubes. K 3 EDTA samples stored at 4°C showed less variation than room temperature storage, but levels were elevated compared with BCT. A second centrifugation at 3000 × g gave similar cfDNA yields compared with higher-speed centrifugation. Next-generation sequencing showed negligible differences in background error or copy number changes between K 3 EDTA and BCT, or following shipment in BCT. This study provides insights into the effects of sample processing on ctDNA analysis., (Copyright © 2018 American Society for Investigative Pathology and the Association for Molecular Pathology. Published by Elsevier Inc. All rights reserved.)

Tumour heterogeneity leads to the development of multiple resistance mechanisms during targeted therapies. Identifying the dominant driver(s) is critical for treatment decision. We studied the relative dynamics of multiple oncogenic drivers in longitudinal plasma of 50 EGFR -mutant non-small-cell lung cancer patients receiving gefitinib and hydroxychloroquine. We performed digital PCR and targeted sequencing on samples from all patients and shallow whole-genome sequencing on samples from three patients who underwent histological transformation to small-cell lung cancer. In 43 patients with known EGFR mutations from tumour, we identified them accurately in plasma of 41 patients (95%, 41/43). We also found additional mutations, including EGFR T790M (31/50, 62%), TP53 (23/50, 46%), PIK3CA (7/50, 14%) and PTEN (4/50, 8%). Patients with both TP53 and EGFR mutations before treatment had worse overall survival than those with only EGFR Patients who progressed without T790M had worse PFS during TKI continuation and developed alternative alterations, including small-cell lung cancer-associated copy number changes and TP53 mutations, that tracked subsequent treatment responses. Longitudinal plasma analysis can help identify dominant resistance mechanisms, including non-druggable genetic information that may guide clinical management., (© 2018 The Authors. Published under the terms of the CC BY 4.0 license.)

Background: An increasing number of precision oncology programmes are being launched world-wide. To support this development, we present the Cancer Variant Explorer (CVE), an R package with an interactive Shiny web browser interface., Results: Leveraging Oncotator and the Drug Gene Interaction Database, CVE offers exploration of variants within single or multiple tumour exomes to identify drivers, resistance mechanisms and to assess druggability. We present example applications including the analysis of an individual patient and a cohort-wide study, and provide a first extension of CVE by adding a tumour-specific co-expression network., Conclusions: The CVE package allows interactive variant prioritisation to expedite the analysis of cancer sequencing studies. Our framework also includes the prioritisation of druggable targets, allows exploratory analysis of tissue specific networks and is extendable for specific applications by virtue of its modular design. We encourage the use of CVE within translational research studies and molecular tumour boards. The CVE package is available via Bioconductor ( http://bioconductor.org/packages/CVE/ ).

TP53 mutations are ubiquitous in high-grade serous ovarian carcinomas (HGSOC), and the presence of TP53 mutation discriminates between high and low-grade serous carcinomas and is now an important biomarker for clinical trials targeting mutant p53. p53 immunohistochemistry (IHC) is widely used as a surrogate for TP53 mutation but its accuracy has not been established. The objective of this study was to test whether improved methods for p53 IHC could reliably predict TP53 mutations independently identified by next generation sequencing (NGS). Four clinical p53 IHC assays and tagged-amplicon NGS for TP53 were performed on 171 HGSOC and 80 endometrioid carcinomas (EC). p53 expression was scored as overexpression (OE), complete absence (CA), cytoplasmic (CY) or wild type (WT). p53 IHC was evaluated as a binary classifier where any abnormal staining predicted deleterious TP53 mutation and as a ternary classifier where OE, CA or WT staining predicted gain-of-function (GOF or nonsynonymous), loss-of-function (LOF including stopgain, indel, splicing) or no detectable TP53 mutations (NDM), respectively. Deleterious TP53 mutations were detected in 169/171 (99%) HGSOC and 7/80 (8.8%) EC. The overall accuracy for the best performing IHC assay for binary and ternary prediction was 0.94 and 0.91 respectively, which improved to 0.97 (sensitivity 0.96, specificity 1.00) and 0.95 after secondary analysis of discordant cases. The sensitivity for predicting LOF mutations was lower at 0.76 because p53 IHC detected mutant p53 protein in 13 HGSOC with LOF mutations. CY staining associated with LOF was seen in 4 (2.3%) of HGSOC. Optimized p53 IHC can approach 100% specificity for the presence of TP53 mutation and its high negative predictive value is clinically useful as it can exclude the possibility of a low-grade serous tumour. 4.1% of HGSOC cases have detectable WT staining while harboring a TP53 LOF mutation, which limits sensitivity for binary prediction of mutation to 96%.

Circulating tumour DNA analysis can be used to track tumour burden and analyse cancer genomes non-invasively but the extent to which it represents metastatic heterogeneity is unknown. Here we follow a patient with metastatic ER-positive and HER2-positive breast cancer receiving two lines of targeted therapy over 3 years. We characterize genomic architecture and infer clonal evolution in eight tumour biopsies and nine plasma samples collected over 1,193 days of clinical follow-up using exome and targeted amplicon sequencing. Mutation levels in the plasma samples reflect the clonal hierarchy inferred from sequencing of tumour biopsies. Serial changes in circulating levels of sub-clonal private mutations correlate with different treatment responses between metastatic sites. This comparison of biopsy and plasma samples in a single patient with metastatic breast cancer shows that circulating tumour DNA can allow real-time sampling of multifocal clonal evolution.

Tumour heterogeneity in primary prostate cancer is a well-established phenomenon. However, how the subclonal diversity of tumours changes during metastasis and progression to lethality is poorly understood. Here we reveal the precise direction of metastatic spread across four lethal prostate cancer patients using whole-genome and ultra-deep targeted sequencing of longitudinally collected primary and metastatic tumours. We find one case of metastatic spread to the surgical bed causing local recurrence, and another case of cross-metastatic site seeding combining with dynamic remoulding of subclonal mixtures in response to therapy. By ultra-deep sequencing end-stage blood, we detect both metastatic and primary tumour clones, even years after removal of the prostate. Analysis of mutations associated with metastasis reveals an enrichment of TP53 mutations, and additional sequencing of metastases from 19 patients demonstrates that acquisition of TP53 mutations is linked with the expansion of subclones with metastatic potential which we can detect in the blood.

VHL is mutated in the majority of patients with clear cell renal cell carcinoma (ccRCC), with conflicting clinical relevance. Recent studies have identified recurrent mutations in histone modifying and chromatin remodeling genes, including BAP1, PBRM1, SETD2, KDM6A, and JARID1c. Current evidence suggests that BAP1 mutations are associated with aggressive disease. The clinical significance of the remaining genes is unknown. In this study, targeted sequencing of VHL and JARID1c (entire genes) and coding regions of BAP1, PBRM1, SETD2, and KDM6A was performed on 132 ccRCCs and matched normal tissues. Associations between mutations and clinical and pathological outcomes were interrogated. Inactivation of VHL (coding mutation or promoter methylation) was seen in 75% of ccRCCs. Somatic noncoding VHL alterations were identified in 29% of ccRCCs and may be associated with improved overall survival. BAP1 (11%), PBRM1 (33%), SETD2 (16%), JARID1c (4%), and KDM6A (3%) mutations were identified. BAP1-mutated tumors were associated with metastatic disease at presentation (P = 0.023), advanced clinical stage (P = 0.042) and a trend towards shorter recurrence free survival (P = 0.059) when compared with tumors exclusively mutated for PBRM1. Our results support those of recent publications pointing towards a role for BAP1 and PBRM1 mutations in risk stratifying ccRCCs. Further investigation of noncoding alterations in VHL is warranted., (Copyright © 2013 Wiley Periodicals, Inc.)

At least 5% of individuals with hypertension have adrenal aldosterone-producing adenomas (APAs). Gain-of-function mutations in KCNJ5 and apparent loss-of-function mutations in ATP1A1 and ATP2A3 were reported to occur in APAs. We find that KCNJ5 mutations are common in APAs resembling cortisol-secreting cells of the adrenal zona fasciculata but are absent in a subset of APAs resembling the aldosterone-secreting cells of the adrenal zona glomerulosa. We performed exome sequencing of ten zona glomerulosa-like APAs and identified nine with somatic mutations in either ATP1A1, encoding the Na(+)/K(+) ATPase α1 subunit, or CACNA1D, encoding Cav1.3. The ATP1A1 mutations all caused inward leak currents under physiological conditions, and the CACNA1D mutations induced a shift of voltage-dependent gating to more negative voltages, suppressed inactivation or increased currents. Many APAs with these mutations were <1 cm in diameter and had been overlooked on conventional adrenal imaging. Recognition of the distinct genotype and phenotype for this subset of APAs could facilitate diagnosis.

Cancers acquire resistance to systemic treatment as a result of clonal evolution and selection. Repeat biopsies to study genomic evolution as a result of therapy are difficult, invasive and may be confounded by intra-tumour heterogeneity. Recent studies have shown that genomic alterations in solid cancers can be characterized by massively parallel sequencing of circulating cell-free tumour DNA released from cancer cells into plasma, representing a non-invasive liquid biopsy. Here we report sequencing of cancer exomes in serial plasma samples to track genomic evolution of metastatic cancers in response to therapy. Six patients with advanced breast, ovarian and lung cancers were followed over 1-2 years. For each case, exome sequencing was performed on 2-5 plasma samples (19 in total) spanning multiple courses of treatment, at selected time points when the allele fraction of tumour mutations in plasma was high, allowing improved sensitivity. For two cases, synchronous biopsies were also analysed, confirming genome-wide representation of the tumour genome in plasma. Quantification of allele fractions in plasma identified increased representation of mutant alleles in association with emergence of therapy resistance. These included an activating mutation in PIK3CA (phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha) following treatment with paclitaxel; a truncating mutation in RB1 (retinoblastoma 1) following treatment with cisplatin; a truncating mutation in MED1 (mediator complex subunit 1) following treatment with tamoxifen and trastuzumab, and following subsequent treatment with lapatinib, a splicing mutation in GAS6 (growth arrest-specific 6) in the same patient; and a resistance-conferring mutation in EGFR (epidermal growth factor receptor; T790M) following treatment with gefitinib. These results establish proof of principle that exome-wide analysis of circulating tumour DNA could complement current invasive biopsy approaches to identify mutations associated with acquired drug resistance in advanced cancers. Serial analysis of cancer genomes in plasma constitutes a new paradigm for the study of clonal evolution in human cancers.

Background: Atlantic salmon (Salmo salar L.) is an environmentally and economically important organism and its gene content is reasonably well characterized. From a transcriptional standpoint, it is important to characterize the changes in gene expression over the course of unperturbed early development, from fertilization through to the parr stage., Findings: S. salar samples were taken at 17 time points from 2 to 89 days post fertilization. Total RNA was extracted and cRNA was synthesized and hybridized to a newly developed 44K oligo salmonid microarray platform. Quantified results were subjected to preliminary data analysis and submitted to NCBI's Gene Expression Omnibus (GEO). Data can be found under the GEO accession number GSE25938. http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25938, Conclusions: Throughout the entire period of development, several thousand genes were found to be differentially regulated. This work represents the trancriptional characterization of a very large geneset that will be extremely valuable in further examination of the transcriptional changes in Atlantic salmon during the first few months of development. The expression profiles can help to annotate salmon genes in addition to being used as references against any number of experimental variables to which developing salmonids might be subjected.

We describe the use of Java Pattern Finder (JaPaFi) to identify short (<100 nt) highly conserved sequences in a series of poxvirus genomes. The algorithm utilizes pattern matching to identify approximate matches appearing at least once in each member of a set of genomes; a key feature is that the genomes do not need to be aligned. The user simply specifies the genomes to search, minimum length of sequences to find and the maximum number of mismatches and indels allowed. Many of the most highly conserved segments contain poxvirus promoter elements.

Background: The volume of viral genomic sequence data continues to increase rapidly. This is especially true for the smaller RNA viruses, which are relatively easy to sequence in large numbers. The data volumes cause a number of significant problems for research applications that require large multiple alignments of essentially complete genomes, which are of the order of 10 kb., Findings: We present a simple strategy to enable the creation of large quasi-multiple sequence alignments from pairwise alignment data. This process is suitable for large, closely related sequences such as the polyproteins of dengue viruses, which need the insertion of very few indels., Conclusion: The quasi-multiple sequence alignments generated by KISSa are sufficiently accurate to support tree-based genome selection for interactive bioinformatics analysis tools. The speed of this process is critical to providing an interactive experience for the user.

Background: Many sequence-searching tools have limiting factors for their use. For example, they may be platform specific, enforce restrictive size limits and sequences to be searched, or only allow searches of one of DNA or protein., Findings: We present an easy-to-use, fast, platform-independent tool to search for amino acid or nucleotide patterns within one or many protein or nucleic acid sequences. The user can choose to search for regular expressions or perform a fuzzy search in which a particular number of errors is accepted during matching of a sequence. Positions of mismatches in fuzzy searches are displayed graphically the user., Conclusion: SeqS provides an improved feature set and functions as a stand-alone tool or could be integrated into other bioinformatics platforms.