Exon-based phylogenomics strengthens the phylogeny of Neotropical cichlids and identifies remaining conflicting clades (Cichliformes: Cichlidae: Cichlinae).

The phenotypic, geographic, and species diversity of cichlid fishes have made them a group of great interest for studying evolutionary processes. Here we present a targeted-exon next-generation sequencing approach for investigating the evolutionary relationships of cichlid fishes (Cichlidae), with focus on the Neotropical subfamily Cichlinae using a set of 923 primarily single-copy exons designed through mining of the Nile tilapia (Oreochromis niloticus) genome. Sequence capture and assembly were robust, leading to a complete dataset of 415 exons for 139 species (147 terminals) that consisted of 128 Neotropical species, six African taxa, and five Indo-Malagasy cichlids. Gene and species trees were calculated using alternative partitioning schemes and reconstruction methods. In general, all methods yielded similar topologies to previously hypothesized relationships within the Cichlinae and clarified several relationships that were previously poorly supported or in conflict. Additional work will be needed to fully resolve all aspects of Cichlinae phylogeny. Overall, this approach yielded a well-resolved phylogeny of Neotropical cichlids that will be of utility for future assessments of the evolutionary and ecological processes within this diverse group of fishes. Furthermore, the general methodology employed here of exon targeting and capture should be applicable to any group of organisms with the availability of a reference genome.

Genetic and genomic architecture of the evolution of resistance to antifungal drug combinations.

The evolution of drug resistance in fungal pathogens compromises the efficacy of the limited number of antifungal drugs. Drug combinations have emerged as a powerful strategy to enhance antifungal efficacy and abrogate drug resistance, but the impact on the evolution of drug resistance remains largely unexplored. Targeting the molecular chaperone Hsp90 or its downstream effector, the protein phosphatase calcineurin, abrogates resistance to the most widely deployed antifungals, the azoles, which inhibit ergosterol biosynthesis. Here, we evolved experimental populations of the model yeast Saccharomyces cerevisiae and the leading human fungal pathogen Candida albicans with azole and an inhibitor of Hsp90, geldanamycin, or calcineurin, FK506. To recapitulate a clinical context where Hsp90 or calcineurin inhibitors could be utilized in combination with azoles to render resistant pathogens responsive to treatment, the evolution experiment was initiated with strains that are resistant to azoles in a manner that depends on Hsp90 and calcineurin. Of the 290 lineages initiated, most went extinct, yet 14 evolved resistance to the drug combination. Drug target mutations that conferred resistance to geldanamycin or FK506 were identified and validated in five evolved lineages. Whole-genome sequencing identified mutations in a gene encoding a transcriptional activator of drug efflux pumps, PDR1, and a gene encoding a transcriptional repressor of ergosterol biosynthesis genes, MOT3, that transformed azole resistance of two lineages from dependent on calcineurin to independent of this regulator. Resistance also arose by mutation that truncated the catalytic subunit of calcineurin, and by mutation in LCB1, encoding a sphingolipid biosynthetic enzyme. Genome analysis revealed extensive aneuploidy in four of the C. albicans lineages. Thus, we identify molecular determinants of the transition of azole resistance from calcineurin dependence to independence and establish multiple mechanisms by which resistance to drug combinations evolves, providing a foundation for predicting and preventing the evolution of drug resistance.

Early Changes in Tumor-Naive Cell-Free Methylomes and Fragmentomes Predict Outcomes in Pembrolizumab-Treated Solid Tumors.

Early kinetics of circulating tumor DNA (ctDNA) in plasma predict response to pembrolizumab but typically requires sequencing of matched tumor tissue or fixed gene panels. We analyzed genome-wide methylation and fragment-length profiles using cell-free methylated DNA immunoprecipitation and sequencing (cfMeDIP-seq) in 204 plasma samples from 87 patients before and during treatment with pembrolizumab from a pan-cancer phase II investigator-initiated trial (INSPIRE). We trained a pan-cancer methylation signature using independent methylation array data from The Cancer Genome Atlas to quantify cancer-specific methylation (CSM) and fragment-length score (FLS) for each sample. CSM and FLS are strongly correlated with tumor-informed ctDNA levels. Early kinetics of CSM predict overall survival and progression-free survival, independently of tumor type, PD-L1, and tumor mutation burden. Early kinetics of FLS are associated with overall survival independently of CSM. Our tumor-naïve mutation-agnostic ctDNA approach integrating methylomics and fragmentomics could predict outcomes in patients treated with pembrolizumab. SIGNIFICANCE: Analysis of methylation and fragment length in plasma using cfMeDIP-seq provides a tumor-naive approach to measure ctDNA with results comparable with a tumor-informed bespoke ctDNA. Early kinetics within the first weeks of treatment in methylation and fragment quantity can predict outcomes with pembrolizumab in patients with various advanced solid tumors. This article is featured in Selected Articles from This Issue, p. 897.

Early Cancer Detection in Li-Fraumeni Syndrome with Cell-Free DNA.

People with Li-Fraumeni syndrome (LFS) harbor a germline pathogenic variant in the TP53 tumor suppressor gene, face a near 100% lifetime risk of cancer, and routinely undergo intensive surveillance protocols. Liquid biopsy has become an attractive tool for a range of clinical applications, including early cancer detection. Here, we provide a proof-of-principle for a multimodal liquid biopsy assay that integrates a targeted gene panel, shallow whole-genome, and cell-free methylated DNA immunoprecipitation sequencing for the early detection of cancer in a longitudinal cohort of 89 LFS patients. Multimodal analysis increased our detection rate in patients with an active cancer diagnosis over uni-modal analysis and was able to detect cancer-associated signal(s) in carriers prior to diagnosis with conventional screening (positive predictive value = 67.6%, negative predictive value = 96.5%). Although adoption of liquid biopsy into current surveillance will require further clinical validation, this study provides a framework for individuals with LFS. SIGNIFICANCE: By utilizing an integrated cell-free DNA approach, liquid biopsy shows earlier detection of cancer in patients with LFS compared with current clinical surveillance methods such as imaging. Liquid biopsy provides improved accessibility and sensitivity, complementing current clinical surveillance methods to provide better care for these patients. See related commentary by Latham et al., p. 23. This article is featured in Selected Articles from This Issue, p. 5.

Integrated, Longitudinal Analysis of Cell-free DNA in Uveal Melanoma.

Uveal melanomas are rare tumors arising from melanocytes that reside in the eye. Despite surgical or radiation treatment, approximately 50% of patients with uveal melanoma will progress to metastatic disease, most often to the liver. Cell-free DNA (cfDNA) sequencing is a promising technology due to the minimally invasive sample collection and ability to infer multiple aspects of tumor response. We analyzed 46 serial cfDNA samples from 11 patients with uveal melanoma over a 1-year period following enucleation or brachytherapy (n = ∼4/patient) using targeted panel, shallow whole genome, and cell-free methylated DNA immunoprecipitation sequencing. We found detection of relapse was highly variable using independent analyses (P = 0.06-0.46), whereas a logistic regression model integrating all cfDNA profiles significantly improved relapse detection (P = 0.02), with greatest power derived from fragmentomic profiles. This work provides support for the use of integrated analyses to improve the sensitivity of circulating tumor DNA detection using multi-modal cfDNA sequencing. Significance: Here, we demonstrate integrated, longitudinal cfDNA sequencing using multi-omic approaches is more effective than unimodal analysis. This approach supports the use of frequent blood testing using comprehensive genomic, fragmentomic, and epigenomic techniques.

Identifying Mechanisms of Resistance by Circulating Tumor DNA in EVOLVE, a Phase II Trial of Cediranib Plus Olaparib for Ovarian Cancer at Time of PARP Inhibitor Progression.

PURPOSE: To evaluate the use of blood cell-free DNA (cfDNA) to identify emerging mechanisms of resistance to PARP inhibitors (PARPi) in high-grade serous ovarian cancer (HGSOC). EXPERIMENTAL DESIGN: We used targeted sequencing (TS) to analyze 78 longitudinal cfDNA samples collected from 30 patients with HGSOC enrolled in a phase II clinical trial evaluating cediranib (VEGF inhibitor) plus olaparib (PARPi) after progression on PARPi alone. cfDNA was collected at baseline, before treatment cycle 2, and at end of treatment. These were compared with whole-exome sequencing (WES) of baseline tumor tissues. RESULTS: At baseline (time of initial PARPi progression), cfDNA tumor fractions were 0.2% to 67% (median, 3.25%), and patients with high ctDNA levels (>15%) had a higher tumor burden (sum of target lesions; P = 0.043). Across all timepoints, cfDNA detected 74.4% of mutations known from prior tumor WES, including three of five expected BRCA1/2 reversion mutations. In addition, cfDNA identified 10 novel mutations not detected by WES, including seven TP53 mutations annotated as pathogenic by ClinVar. cfDNA fragmentation analysis attributed five of these novel TP53 mutations to clonal hematopoiesis of indeterminate potential (CHIP). At baseline, samples with significant differences in mutant fragment size distribution had shorter time to progression (P = 0.001). CONCLUSIONS: Longitudinal testing of cfDNA by TS provides a noninvasive tool for detection of tumor-derived mutations and mechanisms of PARPi resistance that may aid in directing patients to appropriate therapeutic strategies. With cfDNA fragmentation analyses, CHIP was identified in several patients and warrants further investigation.

Genomic Landscape of Malignant Peripheral Nerve Sheath Tumor‒Like Melanoma.

Malignant peripheral nerve sheath tumor (MPNST)-like melanoma is a rare malignancy with overlapping characteristics of both neural sarcoma and melanoma. Although the genomics of cutaneous melanoma has been extensively studied, those of MPNST-like melanoma have not. To characterize the genomic landscape of MPNST-like melanoma, we performed a single-center, retrospective cohort study at a tertiary academic cancer center. Consecutive patients with a confirmed histologic diagnosis of MPNST-like melanoma were screened, and those whose tissues were locally available were included in this analysis. Archival tissue from six patients (eight samples) was submitted for whole-exome and transcriptome sequencing analysis. We compared these data with available genomic studies of cutaneous melanoma and MPNST. NF1 was altered (mutated, deleted, or amplified) in 67% of patients. Genes related to cell cycle regulation were frequently altered, with frequent deletion of ZNF331, which, to the best of our knowledge, has not been previously described in cutaneous melanoma. The serine protease inhibitor SERPINB4 was deleted in 100% of the patients. We show that MPNST-like melanoma presents overlapping genomic features with cutaneous melanoma and MPNST, but it is unique by the frequency of loss of function of ZNF331 and SERPINB4.

Applications of Circulating Tumor DNA in a Cohort of Phase I Solid Tumor Patients Treated With Immunotherapy.

Background: The correlation between blood-based tumor mutation burden (bTMB) and tissue-based tumor mutation burden(tTMB) has not been broadly tested in a multicancer cohort. Here, we assess the correlation between bTMB with tTMB in phase I trial patients treated with immunotherapy. As an exploratory analysis, we evaluated circulating tumor DNA (ctDNA) dynamics in responders. Methods: Patients treated with immunotherapy at the Princess Margaret phase I trials unit were enrolled. Pretreatment plasma ctDNA and matched normal blood controls were collected. Available archival tissue formalin-fixed paraffin-embedded (FFPE) samples were analyzed. A 425-gene panel was used to sequence both ctDNA and FFPE samples. Samples with TMB within the highest tertile were considered as high TMB. Results: Thirty-eight patients were accrued from 25 different trials, 86.8% of which involved an anti-PD-1/PD-L1 agent. Thirty patients (78.9%) had detectable mutations in ctDNA, of which the median (range) bTMB was 5 (1-53) mutations per megabase (mut/Mb). Of the 22 patients with available FFPE samples, mutations were detected in 21 (95.4%); the median (range) tTMB was 6 (2-124) mut/Mb. Among the 16 patients with detectable mutations in both FFPE and ctDNA, a statistically significant correlation between bTMB and tTMB was observed (ρ = 0.71; P = .002). High TMB was not associated with better survival. All 3 responders had a decrease in the variant allele frequency of mutations detected in ctDNA at a second timepoint relative to baseline, indicating a potential early marker of response. Conclusions: In this small series, bTMB correlated with tTMB. An on-treatment decrease in VAF of mutations detected in ctDNA at baseline was observed in responders. Larger studies to verify our findings are warranted.

Tumor genomic, transcriptomic, and immune profiling characterizes differential response to first-line platinum chemotherapy in high grade serous ovarian cancer.

BACKGROUND: In high grade serous ovarian cancer (HGSOC), there is a spectrum of sensitivity to first line platinum-based chemotherapy. This study molecularly characterizes HGSOC patients from two distinct groups of chemotherapy responders (good vs. poor). METHODS: Following primary debulking surgery and intravenous carboplatin/paclitaxel, women with stage III-IV HGSOC were grouped by response. Patients in the good response (GR) and poor response (PR) groups respectively had a progression-free intervals (PFI) of ≥12 and ≤6 months. Analysis of surgical specimens interrogated genomic and immunologic features using whole exome sequencing. RNA-sequencing detected gene expression outliers and inference of immune infiltrate, with validation by targeted NanoString arrays. PD-L1 expression was scored by immunohistochemistry (IHC). RESULTS: A total of 39 patient samples were analyzed (GR = 20; PR = 19). Median PFI for GR and PR patient cohorts was 32 and 3 months, respectively. GR tumors were enriched for loss-of-function BRCA2 mutations and had a significantly higher nonsynonymous mutation rate compared to PR tumors (p = 0.001). Samples from the PR cohort were characterized by mutations in MGA and RAD51B and trended towards a greater rate of amplification of PIK3CA, MECOM, and ATR in comparison to GR tumors. Gene expression analysis by NanoString correlated increased PARP4 with PR and increased PD-L1 and EMSY with GR. There was greater tumor immune cell infiltration and higher immune cell PD-L1 protein expression in the GR group. CONCLUSIONS: Our research demonstrates that tumors from HGSOC patients responding poorly to first line chemotherapy have a distinct molecular profile characterized by actionable drug targets including PARP4.

Antitumor immune effects of preoperative sitravatinib and nivolumab in oral cavity cancer: SNOW window-of-opportunity study.

BACKGROUND: Sitravatinib, a tyrosine kinase inhibitor that targets TYRO3, AXL, MERTK and the VEGF receptor family, is predicted to increase the M1 to M2-polarized tumor-associated macrophages ratio in the tumor microenvironment and have synergistic antitumor activity in combination with anti-programmed death-1/ligand-1 agents. SNOW is a window-of-opportunity study designed to evaluate the immune and molecular effects of preoperative sitravatinib and nivolumab in patients with oral cavity squamous cell carcinoma. METHODS: Patients with newly-diagnosed untreated T2-4a, N0-2 or T1 >1 cm-N2 oral cavity carcinomas were eligible. All patients received sitravatinib 120 mg daily from day 1 up to 48 hours pre-surgery and one dose of nivolumab 240 mg on day 15. Surgery was planned between day 23 and 30. Standard of care adjuvant radiotherapy was given based on clinical stage. Tumor photographs, fresh tumor biopsies and blood samples were collected at baseline, at day 15 after sitravatinib alone, and at surgery after sitravatinib-nivolumab combination. Tumor flow cytometry, multiplex immunofluorescence staining and single-cell RNA sequencing (scRNAseq) were performed on tumor biopsies to study changes in immune-cell populations. Tumor whole-exome sequencing and circulating tumor DNA and cell-free DNA were evaluated at each time point. RESULTS: Ten patients were included. Grade 3 toxicity occurred in one patient (hypertension); one patient required sitravatinib dose reduction, and one patient required discontinuation and surgery delay due to G2 thrombocytopenia. Nine patients had clinical-to-pathological downstaging, with one complete response. Independent pathological treatment response (PTR) assessment confirmed a complete PTR and two major PTRs. With a median follow-up of 21 months, all patients are alive with no recurrence. Circulating tumor DNA and cell-free DNA dynamics correlated with clinical and pathological response and distinguished two patient groups with different tumor biological behavior after sitravatinib alone (1A) versus sitravatinib-nivolumab (1B). Tumor immunophenotyping and scRNAseq analyses revealed differential changes in the expression of immune cell populations and sitravatinib-targeted and hypoxia-related genes in group 1A vs 1B patients. CONCLUSIONS: The SNOW study shows sitravatinib plus nivolumab is safe and leads to deep clinical and pathological responses in oral cavity carcinomas. Multi-omic biomarker analyses dissect the differential molecular effects of sitravatinib versus the sitravatinib-nivolumab and revealed patients with distinct tumor biology behavior. TRIAL REGISTRATION NUMBER: NCT03575598.

Personalized circulating tumor DNA analysis as a predictive biomarker in solid tumor patients treated with pembrolizumab.

Immune checkpoint blockade (ICB) provides clinical benefit to a subset of patients with cancer. However, existing biomarkers do not reliably predict treatment response across diverse cancer types. Limited data exist to show how serial circulating tumor DNA (ctDNA) testing may perform as a predictive biomarker in patients receiving ICB. We conducted a prospective phase II clinical trial to assess ctDNA in five distinct cohorts of patients with advanced solid tumors treated with pembrolizumab (NCT02644369). We applied bespoke ctDNA assays to 316 serial plasma samples obtained at baseline and every three cycles from 94 patients. Baseline ctDNA concentration correlated with progression-free survival, overall survival, clinical response and clinical benefit. This association became stronger when considering ctDNA kinetics during treatment. All 12 patients with ctDNA clearance during treatment were alive with median 25 months follow up. This study demonstrates the potential for broad clinical utility of ctDNA-based surveillance in patients treated with ICB.

Clinical impact of combined epigenetic and molecular analysis of pediatric low-grade gliomas.

BACKGROUND: Both genetic and methylation analysis have been shown to provide insight into the diagnosis and prognosis of many brain tumors. However, the implication of methylation profiling and its interaction with genetic alterations in pediatric low-grade gliomas (PLGGs) are unclear. METHODS: We performed a comprehensive analysis of PLGG with long-term clinical follow-up. In total 152 PLGGs were analyzed from a range of pathological subtypes, including 40 gangliogliomas. Complete molecular analysis was compared with genome-wide methylation data and outcome in all patients. For further analysis of specific PLGG groups, including BRAF p.V600E mutant gliomas, we compiled an additional cohort of clinically and genetically defined tumors from 3 large centers. RESULTS: Unsupervised hierarchical clustering revealed 5 novel subgroups of PLGG. These were dominated by nonneoplastic factors such as tumor location and lymphocytic infiltration. Midline PLGG clustered together while deep hemispheric lesions differed from lesions in the periphery. Mutations were distributed throughout these location-driven clusters of PLGG. A novel methylation cluster suggesting high lymphocyte infiltration was confirmed pathologically and exhibited worse progression-free survival compared with PLGG harboring similar molecular alterations (P = 0.008; multivariate analysis: P = 0.035). Although the current methylation classifier revealed low confidence in 44% of cases and failed to add information in most PLGG, it was helpful in reclassifying rare cases. The addition of histopathological and molecular information to specific methylation subgroups such as pleomorphic xanthoastrocytoma-like tumors could stratify these tumors into low and high risk (P = 0.0014). CONCLUSION: The PLGG methylome is affected by multiple nonneoplastic factors. Combined molecular and pathological analysis is key to provide additional information when methylation classification is used for PLGG in the clinical setting.

CpG Island microarray probe sequences derived from a physical library are representative of CpG Islands annotated on the human genome.

An effective tool for the global analysis of both DNA methylation status and protein-chromatin interactions is a microarray constructed with sequences containing regulatory elements. One type of array suited for this purpose takes advantage of the strong association between CpG Islands (CGIs) and gene regulatory regions. We have obtained 20,736 clones from a CGI Library and used these to construct CGI arrays. The utility of this library requires proper annotation and assessment of the clones, including CpG content, genomic origin and proximity to neighboring genes. Alignment of clone sequences to the human genome (UCSC hg17) identified 9595 distinct genomic loci; 64% were defined by a single clone while the remaining 36% were represented by multiple, redundant clones. Approximately 68% of the loci were located near a transcription start site. The distribution of these loci covered all 23 chromosomes, with 63% overlapping a bioinformatically identified CGI. The high representation of genomic CGI in this rich collection of clones supports the utilization of microarrays produced with this library for the study of global epigenetic mechanisms and protein-chromatin interactions. A browsable database is available on-line to facilitate exploration of the CGIs in this library and their association with annotated genes or promoter elements.

Long read nanopore sequencing for detection of HLA and CYP2D6 variants and haplotypes.

Haplotypes are often critical for the interpretation of genetic laboratory observations into medically actionable findings. Current massively parallel DNA sequencing technologies produce short sequence reads that are often unable to resolve haplotype information. Phasing short read data typically requires supplemental statistical phasing based on known haplotype structure in the population or parental genotypic data. Here we demonstrate that the MinION nanopore sequencer is capable of producing very long reads to resolve both variants and haplotypes of HLA-A, HLA-B and CYP2D6 genes important in determining patient drug response in sample NA12878 of CEPH/UTAH pedigree 1463, without the need for statistical phasing. Long read data from a single 24-hour nanopore sequencing run was used to reconstruct haplotypes, which were confirmed by HapMap data and statistically phased Complete Genomics and Sequenom genotypes. Our results demonstrate that nanopore sequencing is an emerging standalone technology with potential utility in a clinical environment to aid in medical decision-making.

A genome scale overexpression screen to reveal drug activity in human cells.

Target identification is a critical step in the lengthy and expensive process of drug development. Here, we describe a genome-wide screening platform that uses systematic overexpression of pooled human ORFs to understand drug mode-of-action and resistance mechanisms. We first calibrated our screen with the well-characterized drug methotrexate. We then identified new genes involved in the bioactivity of diverse drugs including antineoplastic agents and biologically active molecules. Finally, we focused on the transcription factor RHOXF2 whose overexpression conferred resistance to DNA damaging agents. This approach represents an orthogonal method for functional screening and, to our knowledge, has never been reported before.

Miniature short hairpin RNA screens to characterize antiproliferative drugs.

The application of new proteomics and genomics technologies support a view in which few drugs act solely by inhibiting a single cellular target. Indeed, drug activity is modulated by complex, often incompletely understood cellular mechanisms. Therefore, efforts to decipher mode of action through genetic perturbation such as RNAi typically yields "hits" that fall into several categories. Of particular interest to the present study, we aimed to characterize secondary activities of drugs on cells. Inhibiting a known target can result in clinically relevant synthetic phenotypes. In one scenario, drug perturbation could, for example, improperly activate a protein that normally inhibits a particular kinase. In other cases, additional, lower affinity targets can be inhibited as in the example of inhibition of c-Kit observed in Bcr-Abl-positive cells treated with Gleevec. Drug transport and metabolism also play an important role in the way any chemicals act within the cells. Finally, RNAi per se can also affect cell fitness by more general off-target effects, e.g., via the modulation of apoptosis or DNA damage repair. Regardless of the root cause of these unwanted effects, understanding the scope of a drug's activity and polypharmacology is essential for better understanding its mechanism(s) of action, and such information can guide development of improved therapies. We describe a rapid, cost-effective approach to characterize primary and secondary effects of small-molecules by using small-scale libraries of virally integrated short hairpin RNAs. We demonstrate this principle using a "minipool" composed of shRNAs that target the genes encoding the reported protein targets of approved drugs. Among the 28 known reported drug-target pairs, we successfully identify 40% of the targets described in the literature and uncover several unanticipated drug-target interactions based on drug-induced synthetic lethality. We provide a detailed protocol for performing such screens and for analyzing the data. This cost-effective approach to mammalian knockdown screens, combined with the increasing maturation of RNAi technology will expand the accessibility of similar approaches in academic settings.

Chromatin is an ancient innovation conserved between Archaea and Eukarya.

The eukaryotic nucleosome is the fundamental unit of chromatin, comprising a protein octamer that wraps ∼147 bp of DNA and has essential roles in DNA compaction, replication and gene expression. Nucleosomes and chromatin have historically been considered to be unique to eukaryotes, yet studies of select archaea have identified homologs of histone proteins that assemble into tetrameric nucleosomes. Here we report the first archaeal genome-wide nucleosome occupancy map, as observed in the halophile Haloferax volcanii. Nucleosome occupancy was compared with gene expression by compiling a comprehensive transcriptome of Hfx. volcanii. We found that archaeal transcripts possess hallmarks of eukaryotic chromatin structure: nucleosome-depleted regions at transcriptional start sites and conserved -1 and +1 promoter nucleosomes. Our observations demonstrate that histones and chromatin architecture evolved before the divergence of Archaea and Eukarya, suggesting that the fundamental role of chromatin in the regulation of gene expression is ancient.DOI:http://dx.doi.org/10.7554/eLife.00078.001.