Somatic retrotransposition in human cancer revealed by whole-genome and exome sequencing

Comprehensive genomic profiles of small cell lung cancer

We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. In nearly all the tumours analysed we found bi-allelic inactivation of TP53 and RB1, sometimes by complex genomic rearrangements. Two tumours with wild-type RB1 had evidence of chromothripsis leading to overexpression of cyclin D1 (encoded by the CCND1 gene), revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumour suppressors TP53 and RB1 is obligatory in SCLC. We discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In rare cases, SCLC tumours exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signalling in a pre-clinical SCLC mouse model strikingly reduced the number of tumours and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

Mobile element insertions are frequent in oesophageal adenocarcinomas and can mislead paired-end sequencing analysis

BACKGROUND

Mobile elements are active in the human genome, both in the germline and cancers, where they can mutate driver genes.

RESULTS

While analysing whole genome paired-end sequencing of oesophageal adenocarcinomas to find genomic rearrangements, we identified three ways in which new mobile element insertions appear in the data, resembling translocation or insertion junctions: inserts where unique sequence has been transduced by an L1 (Long interspersed element 1) mobile element; novel inserts that are confidently, but often incorrectly, mapped by alignment software to L1s or polyA tracts in the reference sequence; and a combination of these two ways, where different sequences within one insert are mapped to different loci. We identified nine unique sequences that were transduced by neighbouring L1s, both L1s in the reference genome and L1s not present in the reference. Many of the resulting inserts were small fragments that include little or no recognisable mobile element sequence. We found 6 loci in the reference genome to which sequence reads from inserts were frequently mapped, probably erroneously, by alignment software: these were either L1 sequence or particularly long polyA runs. Inserts identified from such apparent rearrangement junctions averaged 16 inserts/tumour, range 0-153 insertions in 43 tumours. However, many inserts would not be detected by mapping the sequences to the reference genome, because they do not include sufficient mappable sequence. To estimate total somatic inserts we searched for polyA sequences that were not present in the matched normal or other normals from the same tumour batch, and were not associated with known polymorphisms. Samples of these candidate inserts were verified by sequencing across them or manual inspection of surrounding reads: at least 85 % were somatic and resembled L1-mediated events, most including L1Hs sequence. Approximately 100 such inserts were detected per tumour on average (range zero to approximately 700).

CONCLUSIONS

Somatic mobile elements insertions are abundant in these tumours, with over 75 % of cases having a number of novel inserts detected. The inserts create a variety of problems for the interpretation of paired-end sequencing data.

Altering Genomic Integrity: Heavy Metal Exposure Promotes Transposable Element-Mediated Damage

Maintenance of genomic integrity is critical for cellular homeostasis and survival. The active transposable elements (TEs) composed primarily of three mobile element lineages LINE-1, Alu, and SVA comprise approximately 30% of the mass of the human genome. For the past 2 decades, studies have shown that TEs significantly contribute to genetic instability and that TE-caused damages are associated with genetic diseases and cancer. Different environmental exposures, including several heavy metals, influence how TEs interact with its host genome increasing their negative impact. This mini-review provides some basic knowledge on TEs, their contribution to disease, and an overview of the current knowledge on how heavy metals influence TE-mediated damage.

Response of transposable elements to environmental stressors

Transposable elements (TEs) comprise a group of repetitive sequences that bring positive, negative, as well as neutral effects to the host organism. Earlier considered as "junk DNA," TEs are now well-accepted driving forces of evolution and critical regulators of the expression of genetic information. Their activity is regulated by epigenetic mechanisms, including methylation of DNA and histone modifications. The loss of epigenetic control over TEs, exhibited as loss of DNA methylation and decondensation of the chromatin structure, may result in TEs reactivation, initiation of their insertional mutagenesis (retrotransposition) and has been reported in numerous human diseases, including cancer. Accumulating evidence suggests that these alterations are not the simple consequences of the disease, but often may drive the pathogenesis, as they can be detected early during disease development. Knowledge derived from the in vitro, in vivo, and epidemiological studies, clearly demonstrates that exposure to ubiquitous environmental stressors, many of which are carcinogens or suspected carcinogens, are capable of causing alterations in methylation and expression of TEs and initiate retrotransposition events. Evidence summarized in this review suggests that TEs are the sensitive endpoints for detection of effects caused by such environmental stressors, as ionizing radiation (terrestrial, space, and UV-radiation), air pollution (including particulate matter [PM]-derived and gaseous), persistent organic pollutants, and metals. Furthermore, the significance of these effects is characterized by their early appearance, persistence and presence in both, target organs and peripheral blood. Altogether, these findings suggest that TEs may potentially be introduced into safety and risk assessment and serve as biomarkers of exposure to environmental stressors. Furthermore, TEs also show significant potential to become invaluable surrogate biomarkers in clinic and possible targets for therapeutic modalities for disease treatment and prevention.

SoloDel: a probabilistic model for detecting low-frequent somatic deletions from unmatched sequencing data

MOTIVATION

Finding somatic mutations from massively parallel sequencing data is becoming a standard process in genome-based biomedical studies. There are a number of robust methods developed for detecting somatic single nucleotide variations However, detection of somatic copy number alteration has been substantially less explored and remains vulnerable to frequently raised sampling issues: low frequency in cell population and absence of the matched control samples.

RESULTS

We developed a novel computational method SoloDel that accurately classifies low-frequent somatic deletions from germline ones with or without matched control samples. We first constructed a probabilistic, somatic mutation progression model that describes the occurrence and propagation of the event in the cellular lineage of the sample. We then built a Gaussian mixture model to represent the mixed population of somatic and germline deletions. Parameters of the mixture model could be estimated using the expectation-maximization algorithm with the observed distribution of read-depth ratios at the points of discordant-read based initial deletion calls. Combined with conventional structural variation caller, SoloDel greatly increased the accuracy in classifying somatic mutations. Even without control, SoloDel maintained a comparable performance in a wide range of mutated subpopulation size (10-70%). SoloDel could also successfully recall experimentally validated somatic deletions from previously reported neuropsychiatric whole-genome sequencing data.

AVAILABILITY AND IMPLEMENTATION

Java-based implementation of the method is available at http://sourceforge.net/projects/solodel/

CONTACT

swkim@yuhs.ac or dhlee@biosoft.kaist.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Atrx promotes heterochromatin formation at retrotransposons

More than 50% of mammalian genomes consist of retrotransposon sequences. Silencing of retrotransposons by heterochromatin is essential to ensure genomic stability and transcriptional integrity. Here, we identified a short sequence element in intracisternal A particle (IAP) retrotransposons that is sufficient to trigger heterochromatin formation. We used this sequence in a genome-wide shRNA screen and identified the chromatin remodeler Atrx as a novel regulator of IAP silencing. Atrx binds to IAP elements and is necessary for efficient heterochromatin formation. In addition, Atrx facilitates a robust and largely inaccessible heterochromatin structure as Atrx knockout cells display increased chromatin accessibility at retrotransposons and non-repetitive heterochromatic loci. In summary, we demonstrate a direct role of Atrx in the establishment and robust maintenance of heterochromatin.

Somatic structural variation and cancer

With the advent of next-generation sequencing technologies, we have witnessed a rapid pace of discovery of new patterns of somatic structural variation in cancer genomes, and an attempt to figure out their underlying mechanisms. Some of these mechanisms are associated with particular cancer types, and in some cases are the main cause of the structural mutations that drive the oncogenic process. This review provides an overview of the patterns of somatic structural variation and chromosomal structures that characterize cancer genomes, their causal mechanisms and their impact in oncogenesis.

Sequencing, identification and mapping of primed L1 elements (SIMPLE) reveals significant variation in full length L1 elements between individuals

BACKGROUND

There are over a half a million copies of L1 retroelements in the human genome which are responsible for as much as 0.5% of new human genetic diseases. Most new L1 inserts arise from young source elements that are polymorphic in the human genome. Highly active polymorphic "hot" L1 source elements have been shown to be capable of extremely high levels of mobilization and result in numerous instances of disease. Additionally, hot polymorphic L1s have been described to be highly active within numerous cancer genomes. These hot L1s result in mutagenesis by insertion of new L1 copies elsewhere in the genome, but also have been shown to generate additional full length L1 insertions which are also hot and able to further retrotranspose. Through this mechanism, hot L1s may amplify within a tumor and result in a continued cycle of mutagenesis.

RESULTS AND CONCLUSIONS

We have developed a method to detect full-length, polymorphic L1 elements using a targeted next generation sequencing approach, Sequencing Identification and Mapping of Primed L1 Elements (SIMPLE). SIMPLE has 94% sensitivity and detects nearly all full-length L1 elements in a genome. SIMPLE will allow researchers to identify hot mutagenic full-length L1s as potential drivers of genome instability. Using SIMPLE we find that the typical individual has approximately 100 non-reference, polymorphic L1 elements in their genome. These elements are at relatively low population frequencies relative to previously identified polymorphic L1 elements and demonstrate the tremendous diversity in potentially active L1 elements in the human population.

Identification of novel exonic mobile element insertions in epithelial ovarian cancers

Mobile elements comprise about half of the human genome. Three active mobile element families (L1, Alu, and SVA) possibly cause diseases such as cancer. We conducted mobile element insertion (MEI) profiling of 44 epithelial ovarian cancers using exome-sequencing data. We identified a total of 106 MEIs using the Mobster program, 8 of which were novel exonic MEIs.

A droplet digital PCR detection method for rare L1 insertions in tumors

Background

The active human mobile element, long interspersed element 1 (L1) currently populates human genomes in excess of 500,000 copies per haploid genome. Through its mobility via a process called target primed reverse transcription (TPRT), L1 mobilization has resulted in over 100 de novo cases of human disease and has recently been associated with various cancer types. Large advances in high-throughput sequencing (HTS) technology have allowed for an increased understanding of the role of L1 in human cancer; however, researchers are still limited by the ability to validate potentially rare L1 insertion events detected by HTS that may occur in only a small fraction of tumor cells. Additionally, HTS detection of rare events varies greatly as a function of read depth, and new tools for de novo element discovery are needed to fill in gaps created by HTS.

Results

We have employed droplet digital PCR (ddPCR) to detect rare L1 loci in mosaic human genomes. Our assay allows for the detection of L1 insertions as rare as one cell in every 10,000.

Conclusions

ddPCR represents a robust method to be used alongside HTS techniques for detecting, validating and quantitating rare L1 insertion events in tumors and other tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-014-0030-4) contains supplementary material, which is available to authorized users.

Wip1 phosphatase in breast cancer

Understanding the factors contributing to tumor initiation, progression and evolution is of paramount significance. Among them, wild-type p53-induced phosphatase 1 (Wip1) is emerging as an important oncogene by virtue of its negative control on several key tumor suppressor pathways. Originally discovered as a p53-regulated gene, Wip1 has been subsequently found amplified and more recently mutated in a significant fraction of human cancers including breast tumors. Recent development in the field further uncovered the utility of anti-Wip1-directed therapies in delaying tumor onset or in reducing the tumor burden. Furthermore, Wip1 could be an important factor that contributes to tumor heterogeneity, suggesting that its inhibition may decrease the rate of cancer evolution. These effects depend on several signaling pathways modulated by Wip1 phosphatase in a spatial and temporal manner. In this review we discuss the recent development in understanding how Wip1 contributes to tumorigenesis with its relevance to breast cancer.

euL1db: the European database of L1HS retrotransposon insertions in humans

Retrotransposons account for almost half of our genome. They are mobile genetics elements—also known as jumping genes—but only the L1HS subfamily of Long Interspersed Nuclear Elements (LINEs) has retained the ability to jump autonomously in modern humans. Their mobilization in germline—but also some somatic tissues—contributes to human genetic diversity and to diseases, such as cancer. Here, we present euL1db, the European database of L1HS retrotransposon insertions in humans (available at http://euL1db.unice.fr). euL1db provides a curated and comprehensive summary of L1HS insertion polymorphisms identified in healthy or pathological human samples and published in peer-reviewed journals. A key feature of euL1db is its sample-wise organization. Hence L1HS insertion polymorphisms are connected to samples, individuals, families and clinical conditions. The current version of euL1db centralizes results obtained in 32 studies. It contains >900 samples, >140 000 sample-wise insertions and almost 9000 distinct merged insertions. euL1db will help understanding the link between L1 retrotransposon insertion polymorphisms and phenotype or disease.