The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758.

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

Improved structural variant interpretation for hereditary cancer susceptibility using long-read sequencing.

PURPOSE: Structural variants (SVs) may be an underestimated cause of hereditary cancer syndromes given the current limitations of short-read next-generation sequencing. Here we investigated the utility of long-read sequencing in resolving germline SVs in cancer susceptibility genes detected through short-read genome sequencing. METHODS: Known or suspected deleterious germline SVs were identified using Illumina genome sequencing across a cohort of 669 advanced cancer patients with paired tumor genome and transcriptome sequencing. Candidate SVs were subsequently assessed by Oxford Nanopore long-read sequencing. RESULTS: Nanopore sequencing confirmed eight simple pathogenic or likely pathogenic SVs, resolving three additional variants whose impact could not be fully elucidated through short-read sequencing. A recurrent sequencing artifact on chromosome 16p13 and one complex rearrangement on chromosome 5q35 were subsequently classified as likely benign, obviating the need for further clinical assessment. Variant configuration was further resolved in one case with a complex pathogenic rearrangement affecting TSC2. CONCLUSION: Our findings demonstrate that long-read sequencing can improve the validation, resolution, and classification of germline SVs. This has important implications for return of results, cascade carrier testing, cancer screening, and prophylactic interventions.

Sample Tracking Using Unique Sequence Controls.

Sample tracking and identity are essential when processing multiple samples in parallel. Sequencing applications often involve high sample numbers, and the data are frequently used in a clinical setting. As such, a simple and accurate intrinsic sample tracking process through a sequencing pipeline is essential. Various solutions have been implemented to verify sample identity, including variant detection at the start and end of the pipeline using arrays or genotyping, bioinformatic comparisons, and optical barcoding of samples. None of these approaches are optimal. To establish a more effective approach using genetic barcoding, we developed a panel of unique DNA sequences cloned into a common vector. A unique DNA sequence is added to the sample when it is first received and can be detected by PCR and/or sequencing at any stage of the process. The control sequences are approximately 200 bases long with low identity to any sequence in the National Center for Biotechnology Information nonredundant database (<30 bases) and contain no long homopolymer (>7) stretches. When a spiked next-generation sequencing library is sequenced, sequence reads derived from this control sequence are generated along with the standard sequencing run and are used to confirm sample identity and determine cross-contamination levels. This approach is used in our targeted clinical diagnostic whole-genome and RNA-sequencing pipelines and is an inexpensive, flexible, and platform-agnostic solution.

The Genome of the Steller Sea Lion (Eumetopias jubatus).

The Steller sea lion is the largest member of the Otariidae family and is found in the coastal waters of the northern Pacific Rim. Here, we present the Steller sea lion genome, determined through DNA sequencing approaches that utilized microfluidic partitioning library construction, as well as nanopore technologies. These methods constructed a highly contiguous assembly with a scaffold N50 length of over 14 megabases, a contig N50 length of over 242 kilobases and a total length of 2.404 gigabases. As a measure of completeness, 95.1% of 4104 highly conserved mammalian genes were found to be complete within the assembly. Further annotation identified 19,668 protein coding genes. The assembled genome sequence and underlying sequence data can be found at the National Center for Biotechnology Information (NCBI) under the BioProject accession number PRJNA475770.

The Genome of the North American Brown Bear or Grizzly: Ursus arctos ssp. horribilis.

The grizzly bear (Ursus arctos ssp. horribilis) represents the largest population of brown bears in North America. Its genome was sequenced using a microfluidic partitioning library construction technique, and these data were supplemented with sequencing from a nanopore-based long read platform. The final assembly was 2.33 Gb with a scaffold N50 of 36.7 Mb, and the genome is of comparable size to that of its close relative the polar bear (2.30 Gb). An analysis using 4104 highly conserved mammalian genes indicated that 96.1% were found to be complete within the assembly. An automated annotation of the genome identified 19,848 protein coding genes. Our study shows that the combination of the two sequencing modalities that we used is sufficient for the construction of highly contiguous reference quality mammalian genomes. The assembled genome sequence and the supporting raw sequence reads are available from the NCBI (National Center for Biotechnology Information) under the bioproject identifier PRJNA493656, and the assembly described in this paper is version QXTK01000000.

Interfaces of Malignant and Immunologic Clonal Dynamics in Ovarian Cancer.

High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.

The Genome of the Beluga Whale (Delphinapterus leucas).

The beluga whale is a cetacean that inhabits arctic and subarctic regions, and is the only living member of the genus Delphinapterus. The genome of the beluga whale was determined using DNA sequencing approaches that employed both microfluidic partitioning library and non-partitioned library construction. The former allowed for the construction of a highly contiguous assembly with a scaffold N50 length of over 19 Mbp and total reconstruction of 2.32 Gbp. To aid our understanding of the functional elements, transcriptome data was also derived from brain, duodenum, heart, lung, spleen, and liver tissue. Assembled sequence and all of the underlying sequence data are available at the National Center for Biotechnology Information (NCBI) under the Bioproject accession number PRJNA360851A.

Spatial heterogeneity in medulloblastoma.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways.

Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.

Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

MEN1 mutations in Hürthle cell (oncocytic) thyroid carcinoma.

CONTEXT AND OBJECTIVE: Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation, and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management. DESIGN AND PATIENTS: Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line, and five Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in the multiple endocrine neoplasia I (MEN1) gene. RESULTS: Here we report the identification of MEN1 loss-of-function mutations in 4% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumors. CONCLUSION: Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of MEN1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis, and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

Clonal analysis via barcoding reveals diverse growth and differentiation of transplanted mouse and human mammary stem cells.

Cellular barcoding offers a powerful approach to characterize the growth and differentiation activity of large numbers of cotransplanted stem cells. Here, we describe a lentiviral genomic-barcoding and analysis strategy and its use to compare the clonal outputs of transplants of purified mouse and human basal mammary epithelial cells. We found that both sources of transplanted cells produced many bilineage mammary epithelial clones in primary recipients, although primary clones containing only one detectable mammary lineage were also common. Interestingly, regardless of the species of origin, many clones evident in secondary recipients were not detected in the primary hosts, and others that were changed from appearing luminal-restricted to appearing bilineage. This barcoding methodology has thus revealed conservation between mice and humans of a previously unknown diversity in the growth and differentiation activities of their basal mammary epithelial cells stimulated to grow in transplanted hosts.

Distinct evolutionary trajectories of primary high-grade serous ovarian cancers revealed through spatial mutational profiling.

High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to the emergence of treatment-resistant subclones. We sought to measure the degree of genomic diversity within primary, untreated HGSCs to examine the natural state of tumour evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on 31 spatially and temporally separated HGSC tumour specimens (six patients), including ovarian masses, distant metastases and fallopian tube lesions. We found widespread intratumoural variation in mutation, copy number and gene expression profiles, with key driver alterations in genes present in only a subset of samples (eg PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range 10.2-91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole-genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology, with common aetiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating that diversity arises at early stages of tumourigenesis. Our results revealed that HGSCs exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug-resistance mechanisms.

Complete genomic landscape of a recurring sporadic parathyroid carcinoma.

Parathyroid carcinoma is a rare endocrine malignancy with an estimated incidence of less than 1 per million population. Excessive secretion of parathyroid hormone, extremely high serum calcium level, and the deleterious effects of hypercalcaemia are the clinical manifestations of the disease. Up to 60% of patients develop multiple disease recurrences and although long-term survival is possible with palliative surgery, permanent remission is rarely achieved. Molecular drivers of sporadic parathyroid carcinoma have remained largely unknown. Previous studies, mostly based on familial cases of the disease, suggested potential roles for the tumour suppressor MEN1 and proto-oncogene RET in benign parathyroid tumourigenesis, while the tumour suppressor HRPT2 and proto-oncogene CCND1 may also act as drivers in parathyroid cancer. Here, we report the complete genomic analysis of a sporadic and recurring parathyroid carcinoma. Mutational landscapes of the primary and recurrent tumour specimens were analysed using high-throughput sequencing technologies. Such molecular profiling allowed for identification of somatic mutations never previously identified in this malignancy. These included single nucleotide point mutations in well-characterized cancer genes such as mTOR, MLL2, CDKN2C, and PIK3CA. Comparison of acquired mutations in patient-matched primary and recurrent tumours revealed loss of PIK3CA activating mutation during the evolution of the tumour from the primary to the recurrence. Structural variations leading to gene fusions and regions of copy loss and gain were identified at a single-base resolution. Loss of the short arm of chromosome 1, along with somatic missense and truncating mutations in CDKN2C and THRAP3, respectively, provides new evidence for the potential role of these genes as tumour suppressors in parathyroid cancer. The key somatic mutations identified in this study can serve as novel diagnostic markers as well as therapeutic targets.

The E3 ubiquitin ligase UBR5 is recurrently mutated in mantle cell lymphoma.

We have recently reported the application of RNAseq to mantle cell lymphoma (MCL) transcriptomes revealing recurrent mutations in NOTCH1. Here we describe the targeted resequencing of 18 genes mutated in this discovery cohort using a larger cohort of MCL tumors. In addition to frequent mutations in ATM, CCND1, TP53, and NOTCH1, mutations were also observed recurrently in MEF2B, TRAF2, and TET2. Interestingly, the third most frequently mutated gene was UBR5, a gene encoding a 2799aa protein, with multiple functions, including E3 ligase activity based on a conserved cysteine residue at the C-terminus. Nonsynonymous mutations were detected in 18% (18/102) of tumors, with 61% of the mutations resulting in frameshifts in, or around, exon 58, predicted to result in the loss of this conserved cysteine residue. The recurrence and clustering of deleterious mutations implicate UBR5 mutations as a critical pathogenic event in a subgroup of MCL.

Clonal evolution of high-grade serous ovarian carcinoma from primary to recurrent disease.

High-grade serous carcinoma (HGSC) is the most common and fatal form of ovarian cancer. While most tumours are highly sensitive to cytoreductive surgery and platinum- and taxane-based chemotherapy, the majority of patients experience recurrence of treatment-resistant tumours. The clonal origin and mutational adaptations associated with recurrent disease are poorly understood. We performed whole exome sequencing on tumour cells harvested from ascites at three time points (primary, first recurrence, and second recurrence) for three HGSC patients receiving standard treatment. Somatic point mutations and small insertions and deletions were identified by comparison to constitutional DNA. The clonal structure and evolution of tumours were inferred from patterns of mutant allele frequencies. TP53 mutations were predominant in all patients at all time points, consistent with the known founder role of this gene. Tumours from all three patients also harboured mutations associated with cell cycle checkpoint function and Golgi vesicle trafficking. There was convergence of germline and somatic variants within the DNA repair, ECM, cell cycle control, and Golgi vesicle pathways. The vast majority of somatic variants found in recurrent tumours were present in primary tumours. Our findings highlight both known and novel pathways that are commonly mutated in HGSC. Moreover, they provide the first evidence at single nucleotide resolution that recurrent HGSC arises from multiple clones present in the primary tumour with negligible accumulation of new mutations during standard treatment.

Genetic alterations activating kinase and cytokine receptor signaling in high-risk acute lymphoblastic leukemia.

Genomic profiling has identified a subtype of high-risk B-progenitor acute lymphoblastic leukemia (B-ALL) with alteration of IKZF1, a gene expression profile similar to BCR-ABL1-positive ALL and poor outcome (Ph-like ALL). The genetic alterations that activate kinase signaling in Ph-like ALL are poorly understood. We performed transcriptome and whole genome sequencing on 15 cases of Ph-like ALL and identified rearrangements involving ABL1, JAK2, PDGFRB, CRLF2, and EPOR, activating mutations of IL7R and FLT3, and deletion of SH2B3, which encodes the JAK2-negative regulator LNK. Importantly, several of these alterations induce transformation that is attenuated with tyrosine kinase inhibitors, suggesting the treatment outcome of these patients may be improved with targeted therapy.

Mutation discovery in regions of segmental cancer genome amplifications with CoNAn-SNV: a mixture model for next generation sequencing of tumors.

Next generation sequencing has now enabled a cost-effective enumeration of the full mutational complement of a tumor genome-in particular single nucleotide variants (SNVs). Most current computational and statistical models for analyzing next generation sequencing data, however, do not account for cancer-specific biological properties, including somatic segmental copy number alterations (CNAs)-which require special treatment of the data. Here we present CoNAn-SNV (Copy Number Annotated SNV): a novel algorithm for the inference of single nucleotide variants (SNVs) that overlap copy number alterations. The method is based on modelling the notion that genomic regions of segmental duplication and amplification induce an extended genotype space where a subset of genotypes will exhibit heavily skewed allelic distributions in SNVs (and therefore render them undetectable by methods that assume diploidy). We introduce the concept of modelling allelic counts from sequencing data using a panel of Binomial mixture models where the number of mixtures for a given locus in the genome is informed by a discrete copy number state given as input. We applied CoNAn-SNV to a previously published whole genome shotgun data set obtained from a lobular breast cancer and show that it is able to discover 21 experimentally revalidated somatic non-synonymous mutations in a lobular breast cancer genome that were not detected using copy number insensitive SNV detection algorithms. Importantly, ROC analysis shows that the increased sensitivity of CoNAn-SNV does not result in disproportionate loss of specificity. This was also supported by analysis of a recently published lymphoma genome with a relatively quiescent karyotype, where CoNAn-SNV showed similar results to other callers except in regions of copy number gain where increased sensitivity was conferred. Our results indicate that in genomically unstable tumors, copy number annotation for SNV detection will be critical to fully characterize the mutational landscape of cancer genomes.

Use of mutation profiles to refine the classification of endometrial carcinomas.

The classification of endometrial carcinomas is based on pathological assessment of tumour cell type; the different cell types (endometrioid, serous, carcinosarcoma, mixed, undifferentiated, and clear cell) are associated with distinct molecular alterations. This current classification system for high-grade subtypes, in particular the distinction between high-grade endometrioid (EEC-3) and serous carcinomas (ESC), is limited in its reproducibility and prognostic abilities. Therefore, a search for specific molecular classifiers to improve endometrial carcinoma subclassification is warranted. We performed target enrichment sequencing on 393 endometrial carcinomas from two large cohorts, sequencing exons from the following nine genes: ARID1A, PPP2R1A, PTEN, PIK3CA, KRAS, CTNNB1, TP53, BRAF, and PPP2R5C. Based on this gene panel, each endometrial carcinoma subtype shows a distinct mutation profile. EEC-3s have significantly different frequencies of PTEN and TP53 mutations when compared to low-grade endometrioid carcinomas. ESCs and EEC-3s are distinct subtypes with significantly different frequencies of mutations in PTEN, ARID1A, PPP2R1A, TP53, and CTNNB1. From the mutation profiles, we were able to identify subtype outliers, ie cases diagnosed morphologically as one subtype but with a mutation profile suggestive of a different subtype. Careful review of these diagnostically challenging cases suggested that the original morphological classification was incorrect in most instances. The molecular profile of carcinosarcomas suggests two distinct mutation profiles for these tumours: endometrioid-type (PTEN, PIK3CA, ARID1A, KRAS mutations) and serous-type (TP53 and PPP2R1A mutations). While this nine-gene panel does not allow for a purely molecularly based classification of endometrial carcinoma, it may prove useful as an adjunct to morphological classification and serve as an aid in the classification of problematic cases. If used in practice, it may lead to improved diagnostic reproducibility and may also serve to stratify patients for targeted therapeutics.

The clonal and mutational evolution spectrum of primary triple-negative breast cancers.

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.