Gene length corrected trimmed mean of M-values (GeTMM) processing of RNA-seq data performs similarly in intersample analyses while improving intrasample comparisons.

BACKGROUND: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. RESULTS: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. CONCLUSIONS: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

Ankyrin repeat and zinc-finger domain-containing 1 mutations are associated with infantile-onset inflammatory bowel disease.

Infantile-onset inflammatory bowel disease (IO IBD) is an invalidating illness with an onset before 2 years of age and has a complex pathophysiology in which genetic factors are important. Homozygosity mapping and whole-exome sequencing in an IO IBD patient and subsequent sequencing of the candidate gene in 12 additional IO IBD patients revealed two patients with two mutated ankyrin repeat and zinc-finger domain-containing 1 (ANKZF1) alleles (homozygous ANKZF1 R585Q mutation and compound heterozygous ANKZF1 E152K and V32_Q87del mutations, respectively) and two patients with one mutated ANKZF1 allele. Although the function of ANKZF1 in mammals had not been previously evaluated, we show that ANKZF1 has an indispensable role in the mitochondrial response to cellular stress. ANKZF1 is located diffusely in the cytoplasm and translocates to the mitochondria upon cellular stress. ANKZF1 depletion reduces mitochondrial integrity and mitochondrial respiration under conditions of cellular stress. The ANKZF1 mutations identified in IO IBD patients with two mutated ANKZF1 alleles result in dysfunctional ANKZF1, as shown by an increased level of apoptosis in patients' lymphocytes, a decrease in mitochondrial respiration in patient fibroblasts with a homozygous ANKZF1 R585Q mutation, and an inability of ANKZF1 R585Q and E152K to rescue the phenotype of yeast deficient in Vms1, the yeast homologue of ANKZF1. These data indicate that loss-of-function mutations in ANKZF1 result in deregulation of mitochondrial integrity, and this may play a pathogenic role in the development of IO IBD.

Chromothripsis in healthy individuals affects multiple protein-coding genes and can result in severe congenital abnormalities in offspring.

Chromothripsis represents an extreme class of complex chromosome rearrangements (CCRs) with major effects on chromosomal architecture. Although recent studies have associated chromothripsis with congenital abnormalities, the incidence and pathogenic effects of this phenomenon require further investigation. Here, we analyzed the genomes of three families in which chromothripsis rearrangements were transmitted from a mother to her child. The chromothripsis in the mothers resulted in completely balanced rearrangements involving 8-23 breakpoint junctions across three to five chromosomes. Two mothers did not show any phenotypic abnormalities, although 3-13 protein-coding genes were affected by breakpoints. Unbalanced but stable transmission of a subset of the derivative chromosomes caused apparently de novo complex copy-number changes in two children. This resulted in gene-dosage changes, which are probably responsible for the severe congenital phenotypes of these two children. In contrast, the third child, who has a severe congenital disease, harbored all three chromothripsis chromosomes from his healthy mother, but one of the chromosomes acquired de novo rearrangements leading to copy-number changes. These results show that the human genome can tolerate extreme reshuffling of chromosomal architecture, including breakage of multiple protein-coding genes, without noticeable phenotypic effects. The presence of chromothripsis in healthy individuals affects reproduction and is expected to substantially increase the risk of miscarriages, abortions, and severe congenital disease.

Characteristics of de novo structural changes in the human genome.

Small insertions and deletions (indels) and large structural variations (SVs) are major contributors to human genetic diversity and disease. However, mutation rates and characteristics of de novo indels and SVs in the general population have remained largely unexplored. We report 332 validated de novo structural changes identified in whole genomes of 250 families, including complex indels, retrotransposon insertions, and interchromosomal events. These data indicate a mutation rate of 2.94 indels (1-20 bp) and 0.16 SVs (>20 bp) per generation. De novo structural changes affect on average 4.1 kbp of genomic sequence and 29 coding bases per generation, which is 91 and 52 times more nucleotides than de novo substitutions, respectively. This contrasts with the equal genomic footprint of inherited SVs and substitutions. An excess of structural changes originated on paternal haplotypes. Additionally, we observed a nonuniform distribution of de novo SVs across offspring. These results reveal the importance of different mutational mechanisms to changes in human genome structure across generations.

Genomic and transcriptomic plasticity in treatment-naive ovarian cancer.

Intra-tumor heterogeneity is a hallmark of many cancers and may lead to therapy resistance or interfere with personalized treatment strategies. Here, we combined topographic mapping of somatic breakpoints and transcriptional profiling to probe intra-tumor heterogeneity of treatment-naïve stage IIIC/IV epithelial ovarian cancer. We observed that most substantial differences in genomic rearrangement landscapes occurred between metastases in the omentum and peritoneum versus tumor sites in the ovaries. Several cancer genes such as NF1, CDKN2A, and FANCD2 were affected by lesion-specific breakpoints. Furthermore, the intra-tumor variability involved different mutational hallmarks including lesion-specific kataegis (local mutation shower coinciding with genomic breakpoints), rearrangement classes, and coding mutations. In one extreme case, we identified two independent TP53 mutations in ovary tumors and omentum/peritoneum metastases, respectively. Examination of gene expression dynamics revealed up-regulation of key cancer pathways including WNT, integrin, chemokine, and Hedgehog signaling in only subsets of tumor samples from the same patient. Finally, we took advantage of the multilevel tumor analysis to understand the effects of genomic breakpoints on qualitative and quantitative gene expression changes. We show that intra-tumor gene expression differences are caused by site-specific genomic alterations, including formation of in-frame fusion genes. These data highlight the plasticity of ovarian cancer genomes, which may contribute to their strong capacity to adapt to changing environmental conditions and give rise to the high rate of recurrent disease following standard treatment regimes.

Same-day genomic and epigenomic diagnosis of brain tumors using real-time nanopore sequencing.

Molecular classification of cancer has entered clinical routine to inform diagnosis, prognosis, and treatment decisions. At the same time, new tumor entities have been identified that cannot be defined histologically. For central nervous system tumors, the current World Health Organization classification explicitly demands molecular testing, e.g., for 1p/19q-codeletion or IDH mutations, to make an integrated histomolecular diagnosis. However, a plethora of sophisticated technologies is currently needed to assess different genomic and epigenomic alterations and turnaround times are in the range of weeks, which makes standardized and widespread implementation difficult and hinders timely decision making. Here, we explored the potential of a pocket-size nanopore sequencing device for multimodal and rapid molecular diagnostics of cancer. Low-pass whole genome sequencing was used to simultaneously generate copy number (CN) and methylation profiles from native tumor DNA in the same sequencing run. Single nucleotide variants in IDH1, IDH2, TP53, H3F3A, and the TERT promoter region were identified using deep amplicon sequencing. Nanopore sequencing yielded ~0.1X genome coverage within 6 h and resulting CN and epigenetic profiles correlated well with matched microarray data. Diagnostically relevant alterations, such as 1p/19q codeletion, and focal amplifications could be recapitulated. Using ad hoc random forests, we could perform supervised pan-cancer classification to distinguish gliomas, medulloblastomas, and brain metastases of different primary sites. Single nucleotide variants in IDH1, IDH2, and H3F3A were identified using deep amplicon sequencing within minutes of sequencing. Detection of TP53 and TERT promoter mutations shows that sequencing of entire genes and GC-rich regions is feasible. Nanopore sequencing allows same-day detection of structural variants, point mutations, and methylation profiling using a single device with negligible capital cost. It outperforms hybridization-based and current sequencing technologies with respect to time to diagnosis and required laboratory equipment and expertise, aiming to make precision medicine possible for every cancer patient, even in resource-restricted settings.

Clinical and molecular characterization of an infant with a tandem duplication and deletion of 19p13.

Copy number variations (CNVs) on the short arm of chromosome 19 are relatively rare. We present a patient with a tandem de novo 3.9 Mb duplication of 19p13.12p13.2 and an adjacent 288 kb deletion of 19p13.12. The CNVs were detected by genome wide SNP-array and confirmed by fluorescence in situ hybridization. Mate-pair sequencing revealed two breakpoint junctions leading to a germline tandem inverted duplication and an adjacent deletion. The patient had a major congenital heart defect and refractory edema leading to metabolic and endocrinological disturbances. Further complications occurred due to refractory chylothorax, severe inflammatory response syndrome, and repeating sepsis. After 2 months, the child died due to intractable respiratory failure. The phenotype of this patient was compared with reported patients with overlapping deletions or duplications. We conclude that the congenital heart defect, respiratory insufficiency, and abnormal neurologic examination are most likely due the contiguous gene deletion/duplication.

The Diverse Effects of Complex Chromosome Rearrangements and Chromothripsis in Cancer Development.

In recent years, enormous progress has been made with respect to the identification of somatic mutations that contribute to cancer development. Mutation types range from small substitutions to large structural genomic rearrangements, including complex reshuffling of the genome. Sets of mutations in individual cancer genomes may show specific signatures, which can be provoked by both exogenous and endogenous forces. One of the most remarkable mutation patterns observed in human cancers involve massive rearrangement of just a few chromosomal regions. This phenomenon has been termed chromothripsis and appears widespread in a multitude of cancer types. Chromothripsis provides a way for cancer to rapidly evolve through a one-off massive change in genome structure as opposed to a gradual process of mutation and selection. This chapter focuses on the origin, prevalence and impact of chromothripsis and related complex genomic rearrangements during cancer development.

Prevalence and clinical implications of chromothripsis in cancer genomes.

PURPOSE OF REVIEW: A variety of mutational mechanisms shape the landscape of somatic mutations in cancer genomes. Although the contribution of single nucleotide mutations is well studied, getting a hold of structural genomic rearrangements is more difficult due to their complexity and diversity in sizes and classes. Here, we discuss the incidence of complex genomic rearrangements and their impact on cancer development and progression. RECENT FINDINGS: Catastrophic genome rearrangements have recently been described in various cancer genomes. Such complex genomic rearrangements may be a result of local shattering of chromosomes followed by reassembly of DNA fragments, a process termed chromothripsis. In addition, DNA replication errors may lead to complex genomic rearrangements in cancer. Complex reshuffling of chromosomes can cause formation of gene fusions, disruption of tumor suppressors, and amplification of oncogenes. Furthermore, the occurrence of chromothripsis has been associated with poor prognosis in neuroblastoma, melanoma, and multiple myeloma. SUMMARY: Complex genomic rearrangements, such as chromothripsis, may affect cancer gene function and thereby have a major impact on cancer progression, prognosis, and therapy response.

MEFV mutations affecting pyrin amino acid 577 cause autosomal dominant autoinflammatory disease.

OBJECTIVES: Autoinflammatory disorders are disorders of the innate immune system. Standard genetic testing provided no correct diagnosis in a female patient from a non-consanguineous family of British descent with a colchicine-responsive autosomal dominant periodic fever syndrome. We aimed to unravel the genetic cause of the symptoms. METHODS: Whole exome sequencing was used to screen for novel sequence variants, which were validated by direct Sanger sequencing. Ex vivo stimulation with peripheral blood mononuclear cells was performed to study the functional consequences of the mutation. mRNA and cytokine levels were measured by quantitative PCR and ELISA, respectively. RESULTS: Whole exome sequencing revealed a novel missense sequence variant, not seen in around 6800 controls, mapping to exon 8 of the MEFV gene (c.1730C>A; p.T577N), co-segregating perfectly with disease in this family. Other mutations at the same amino acid (c.1730C>G; p.T577S and c.1729A>T; p.T577S) were found in a family of Turkish descent, with autosomal dominant inheritance of familial Mediterranean fever (FMF)-like phenotype, and a Dutch patient, respectively. Moreover, a mutation (c.1729A>G; p.T577A) was detected in two Dutch siblings, who had episodes of inflammation of varying severity not resembling FMF. Peripheral blood mononuclear cells from one patient of the index family showed increased basal interleukin 1ß mRNA levels and cytokine responses after lipopolysaccharide stimulation. Responses normalised with colchicine treatment. CONCLUSIONS: Heterozygous mutations at amino acid position 577 of pyrin can induce an autosomal dominant autoinflammatory syndrome. This suggests that T577, located in front of the C-terminal B30.2/SPRY domain, is crucial for pyrin function.

The Neo-Open Reading Frame Peptides That Comprise the Tumor Framome Are a Rich Source of Neoantigens for Cancer Immunotherapy.

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate the whole-genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOP) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of noncoding regions of the genome downstream of a rearrangement breakpoint, i.e., where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for the detection of NOPs.

Improving mammalian genome scaffolding using large insert mate-pair next-generation sequencing.

BACKGROUND: Paired-tag sequencing approaches are commonly used for the analysis of genome structure. However, mammalian genomes have a complex organization with a variety of repetitive elements that complicate comprehensive genome-wide analyses. RESULTS: Here, we systematically assessed the utility of paired-end and mate-pair (MP) next-generation sequencing libraries with insert sizes ranging from 170 bp to 25 kb, for genome coverage and for improving scaffolding of a mammalian genome (Rattus norvegicus). Despite a lower library complexity, large insert MP libraries (20 or 25 kb) provided very high physical genome coverage and were found to efficiently span repeat elements in the genome. Medium-sized (5, 8 or 15 kb) MP libraries were much more efficient for genome structure analysis than the more commonly used shorter insert paired-end and 3 kb MP libraries. Furthermore, the combination of medium- and large insert libraries resulted in a 3-fold increase in N50 in scaffolding processes. Finally, we show that our data can be used to evaluate and improve contig order and orientation in the current rat reference genome assembly. CONCLUSIONS: We conclude that applying combinations of mate-pair libraries with insert sizes that match the distributions of repetitive elements improves contig scaffolding and can contribute to the finishing of draft genomes.

Chromothripsis as a mechanism driving complex de novo structural rearrangements in the germline.

A variety of mutational mechanisms shape the dynamic architecture of human genomes and occasionally result in congenital defects and disease. Here, we used genome-wide long mate-pair sequencing to systematically screen for inherited and de novo structural variation in a trio including a child with severe congenital abnormalities. We identified 4321 inherited structural variants and 17 de novo rearrangements. We characterized the de novo structural changes to the base-pair level revealing a complex series of balanced inter- and intra-chromosomal rearrangements consisting of 12 breakpoints involving chromosomes 1, 4 and 10. Detailed inspection of breakpoint regions indicated that a series of simultaneous double-stranded DNA breaks caused local shattering of chromosomes. Fusion of the resulting chromosomal fragments involved non-homologous end joining, since junction points displayed limited or no homology and small insertions and deletions. The pattern of random joining of chromosomal fragments that we observe here strongly resembles the somatic rearrangement patterns--termed chromothripsis--that have recently been described in deranged cancer cells. We conclude that a similar mechanism may also drive the formation of de novo structural variation in the germline.

Development of a Personalized Tumor Neoantigen Based Vaccine Formulation (FRAME-001) for Use in a Phase II Trial for the Treatment of Advanced Non-Small Cell Lung Cancer.

Stage III-IV non-small cell lung cancer (NSCLC) is a devastating disease characterized by a poor prognosis. NSCLC tumors carry genetic mutations, which can lead to the expression of altered protein sequences. Peptides originating from mutated proteins and bound to MHC molecules on the tumor cell surface are referred to as neoantigens, as they are tumor-specific and not expressed in normal cells. Due to their tumor specificity, neoantigens have a strong potential to induce an anti-tumor immune response and have been investigated for development of personalized therapeutic cancer vaccines. The current study describes the development of a clinical grade neoantigen vaccine formulation (FRAME-001) intended as immunotherapy in advanced NSCLC in combination with the immune checkpoint inhibitor pembrolizumab. The detection of aberrant tumor-specific transcripts as well as an algorithm to select immunogenic neoantigen peptides are described. Subsequently, selected neoantigen peptides were synthesized with a high throughput synthesis platform and aseptically formulated under good manufacturing practice (GMP) conditions into four aqueous peptides mixtures that each contained six neoantigen peptides. A validated stability-indicating analytical method was developed in which we considered the personalized nature of the formulation. An extensive stability study performed either at -25 °C or -80 °C showed that the formulation was stable for up to 32 weeks. The formulation was mixed with the vaccine adjuvant Montanide ISA 51 VG, which yielded the final vaccine emulsion. The stability of the vaccine emulsion was demonstrated using microscopic examination, differential light scattering, and the water-drop test. The presented data show that FRAME-001 is a feasible personalized vaccine formulation for the treatment of stage III-IV NSCLC. The presented data may give guidance in the development of novel personalized therapeutic vaccines since this formulation strategy could be used for any cancer indication.

Distinct Genomic Profiles Are Associated with Treatment Response and Survival in Ovarian Cancer.

The majority of patients with ovarian cancer ultimately develop recurrent chemotherapy-resistant disease. Treatment stratification is mainly based on histological subtype and stage, prior response to platinum-based chemotherapy, and time to recurrent disease. Here, we integrated clinical treatment, treatment response, and survival data with whole-genome sequencing profiles of 132 solid tumor biopsies of metastatic epithelial ovarian cancer to explore genome-informed stratification opportunities. Samples from primary and recurrent disease harbored comparable numbers of single nucleotide variants and structural variants. Mutational signatures represented platinum exposure, homologous recombination deficiency, and aging. Unsupervised hierarchical clustering based on genomic input data identified specific ovarian cancer subgroups, characterized by homologous recombination deficiency, genome stability, and duplications. The clusters exhibited distinct response rates and survival probabilities which could thus potentially be used for genome-informed therapy stratification for more personalized ovarian cancer treatment.

A multi-platform reference for somatic structural variation detection.

Accurate detection of somatic structural variation (SV) in cancer genomes remains a challenging problem. This is in part due to the lack of high-quality, gold-standard datasets that enable the benchmarking of experimental approaches and bioinformatic analysis pipelines. Here, we performed somatic SV analysis of the paired melanoma and normal lymphoblastoid COLO829 cell lines using four different sequencing technologies. Based on the evidence from multiple technologies combined with extensive experimental validation, we compiled a comprehensive set of carefully curated and validated somatic SVs, comprising all SV types. We demonstrate the utility of this resource by determining the SV detection performance as a function of tumor purity and sequence depth, highlighting the importance of assessing these parameters in cancer genomics projects. The truth somatic SV dataset as well as the underlying raw multi-platform sequencing data are freely available and are an important resource for community somatic benchmarking efforts.

The diverse functions of microRNAs in animal development and disease.

MicroRNAs (miRNAs) control gene expression by translational inhibition and destabilization of mRNAs. While hundreds of miRNAs have been found, only a few have been studied in detail. miRNAs have been implicated in tissue morphogenesis, cellular processes like apoptosis, and major signaling pathways. Emerging evidence suggests a direct link between miRNAs and disease, and miRNA expression signatures are associated with various types of cancer. In addition, the gain and loss of miRNA target sites appears to be causal to some genetic disorders. Here, we discuss the current literature on the role of miRNAs in animal development and disease.

miR-451 regulates zebrafish erythroid maturation in vivo via its target gata2.

We demonstrate that in zebrafish, the microRNA miR-451 plays a crucial role in promoting erythroid maturation, in part via its target transcript gata2. Zebrafish miR-144 and miR-451 are processed from a single precursor transcript selectively expressed in erythrocytes. In contrast to other hematopoietic mutants, the zebrafish mutant meunier (mnr) showed intact erythroid specification but diminished miR-144/451 expression. Although erythropoiesis initiated normally in mnr, erythrocyte maturation was morphologically retarded. Morpholino knockdown of miR-451 increased erythrocyte immaturity in wild-type embryos, and miR-451 RNA duplexes partially rescued erythroid maturation in mnr, demonstrating a requirement and role for miR-451 in erythrocyte maturation. mnr provided a selectively miR-144/451-deficient background, facilitating studies to discern miRNA function and validate candidate targets. Among computer-predicted miR-451 targets potentially mediating these biologic effects, the pro-stem cell transcription factor gata2 was an attractive candidate. In vivo reporter assays validated the predicted miR-451/gata2-3'UTR interaction, gata2 down-regulation was delayed in miR-451-knockdown and mnr embryos, and gata2 knockdown partially restored erythroid maturation in mnr, collectively confirming gata2 down-regulation as pivotal for miR-451-driven erythroid maturation. These studies define a new genetic pathway promoting erythroid maturation (mnr/miR-451/gata2) and provide a rare example of partial rescue of a mutant phenotype solely by miRNA overexpression.

Accurate detection of circulating tumor DNA using nanopore consensus sequencing.

Levels of circulating tumor DNA (ctDNA) in liquid biopsies may serve as a sensitive biomarker for real-time, minimally-invasive tumor diagnostics and monitoring. However, detecting ctDNA is challenging, as much fewer than 5% of the cell-free DNA in the blood typically originates from the tumor. To detect lowly abundant ctDNA molecules based on somatic variants, extremely sensitive sequencing methods are required. Here, we describe a new technique, CyclomicsSeq, which is based on Oxford Nanopore sequencing of concatenated copies of a single DNA molecule. Consensus calling of the DNA copies increased the base-calling accuracy ~60×, enabling accurate detection of TP53 mutations at frequencies down to 0.02%. We demonstrate that a TP53-specific CyclomicsSeq assay can be successfully used to monitor tumor burden during treatment for head-and-neck cancer patients. CyclomicsSeq can be applied to any genomic locus and offers an accurate diagnostic liquid biopsy approach that can be implemented in clinical workflows.