Mobile DNA in Health and Disease

Evidence of the pathogenic HERV-W envelope expression in T lymphocytes in association with the respiratory outcome of COVID-19 patients

Background

Despite an impressive effort in clinical research, no standard therapeutic approach for coronavirus disease 2019 (COVID-19) patients has been established, highlighting the need to identify early biomarkers for predicting disease progression and new therapeutic interventions for patient management. The present study aimed to evaluate the involvement of the human endogenous retrovirus -W envelope (HERV-W ENV) in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection considering recent findings that HERVs are activated in response to infectious agents and lead to various immunopathological effects. We analysed HERV-W ENV expression in blood cells of COVID-19 patients in correlation with clinical characteristics and have discussed its potential role in the outcome of the disease.

Methods

We analysed HERV-W ENV expression in blood samples of COVID-19 patients and healthy donors by flow cytometry and quantitative reverse transcriptase PCR analysis, and evaluated its correlation with clinical signs, inflammatory markers, cytokine expression, and disease progression.

Findings

HERV-W ENV was highly expressed in the leukocytes of COVID-19 patients but not in those of healthy donors. Its expression correlated with the markers of T-cell differentiation and exhaustion and blood cytokine levels. The percentage of HERV-W ENV-positive lymphocytes correlated with inflammatory markers and pneumonia severity in COVID-19 patients. Notably, HERV-W ENV expression reflects the respiratory outcome of patients during hospitalization.

Interpretation

Given the known immuno- and neuro-pathogenicity of HERV-W ENV protein, it could promote certain pathogenic features of COVID-19 and therefore serve as a biomarker to predict clinical progression of disease and open to further studies for therapeutic intervention.

Funding

Information available at the end of the manuscript.

De novo structural mutation rates and gamete-of-origin biases revealed through genome sequencing of 2,396 families

Each human genome includes de novo mutations that arose during gametogenesis. While these germline mutations represent a fundamental source of new genetic diversity, they can also create deleterious alleles that impact fitness. Whereas the rate and patterns of point mutations in the human germline are now well understood, far less is known about the frequency and features that impact de novo structural variants (dnSVs). We report a family-based study of germline mutations among 9,599 human genomes from 33 multigenerational CEPH-Utah families and 2,384 families from the Simons Foundation Autism Research Initiative. We find that de novo structural mutations detected by alignment-based, short-read WGS occur at an overall rate of at least 0.160 events per genome in unaffected individuals, and we observe a significantly higher rate (0.206 per genome) in ASD-affected individuals. In both probands and unaffected samples, nearly 73% of de novo structural mutations arose in paternal gametes, and we predict most de novo structural mutations to be caused by mutational mechanisms that do not require sequence homology. After multiple testing correction, we did not observe a statistically significant correlation between parental age and the rate of de novo structural variation in offspring. These results highlight that a spectrum of mutational mechanisms contribute to germline structural mutations and that these mechanisms most likely have markedly different rates and selective pressures than those leading to point mutations.

Intercellular viral spread and intracellular transposition of Drosophila gypsy

It has become increasingly clear that retrotransposons (RTEs) are more widely expressed in somatic tissues than previously appreciated. RTE expression has been implicated in a myriad of biological processes ranging from normal development and aging, to age related diseases such as cancer and neurodegeneration. Long Terminal Repeat (LTR)-RTEs are evolutionary ancestors to, and share many features with, exogenous retroviruses. In fact, many organisms contain endogenous retroviruses (ERVs) derived from exogenous retroviruses that integrated into the germ line. These ERVs are inherited in Mendelian fashion like RTEs, and some retain the ability to transmit between cells like viruses, while others develop the ability to act as RTEs. The process of evolutionary transition between LTR-RTE and retroviruses is thought to involve multiple steps by which the element loses or gains the ability to transmit copies between cells versus the ability to replicate intracellularly. But, typically, these two modes of transmission are incompatible because they require assembly in different sub-cellular compartments. Like murine IAP/IAP-E elements, the gypsy family of retroelements in arthropods appear to sit along this evolutionary transition. Indeed, there is some evidence that gypsy may exhibit retroviral properties. Given that gypsy elements have been found to actively mobilize in neurons and glial cells during normal aging and in models of neurodegeneration, this raises the question of whether gypsy replication in somatic cells occurs via intracellular retrotransposition, intercellular viral spread, or some combination of the two. These modes of replication in somatic tissues would have quite different biological implications. Here, we demonstrate that Drosophila gypsy is capable of both cell-associated and cell-free viral transmission between cultured S2 cells of somatic origin. Further, we demonstrate that the ability of gypsy to move between cells is dependent upon a functional copy of its viral envelope protein. This argues that the gypsy element has transitioned from an RTE into a functional endogenous retrovirus with the acquisition of its envelope gene. On the other hand, we also find that intracellular retrotransposition of the same genomic copy of gypsy can occur in the absence of the Env protein. Thus, gypsy exhibits both intracellular retrotransposition and intercellular viral transmission as modes of replicating its genome.

Human Endogenous Retrovirus Type W Envelope from Multiple Sclerosis Demyelinating Lesions Shows Unique Solubility and Antigenic Characteristics

In multiple sclerosis (MS), human endogenous retrovirus W family (HERV-W) envelope protein, pHERV-W ENV, limits remyelination and induces microglia-mediated neurodegeneration. To better understand its role, we examined the soluble pHERV-W antigen from MS brain lesions detected by specific antibodies. Physico-chemical and antigenic characteristics confirmed differences between pHERV-W ENV and syncytin-1. pHERV-W ENV monomers and trimers remained associated with membranes, while hexamers self-assembled from monomers into a soluble macrostructure involving sulfatides in MS brain. Extracellular hexamers are stabilized by internal hydrophobic bonds and external hydrophilic moieties. HERV-W studies in MS also suggest that this diffusible antigen may correspond to a previously described high-molecular-weight neurotoxic factor secreted by MS B-cells and thus represents a major agonist in MS pathogenesis. Adapted methods are now needed to identify encoding HERV provirus(es) in affected cells DNA. The properties and origin of MS brain pHERV-W ENV soluble antigen will allow a better understanding of the role of HERVs in MS pathogenesis. The present results anyhow pave the way to an accurate detection of the different forms of pHERV-W ENV antigen with appropriate conditions that remained unseen until now.

The "missing heritability"-Problem in psychiatry: Is the interaction of genetics, epigenetics and transposable elements a potential solution?

Psychiatric disorders exhibit an enormous burden on the health care systems worldwide accounting for around one-third of years lost due to disability among adults. Their etiology is largely unknown and diagnostic classification is based on symptomatology and course of illness and not on objective biomarkers. Most psychiatric disorders are moderately to highly heritable. However, it is still unknown what mechanisms may explain the discrepancy between heritability estimates and the present data from genetic analysis. In addition to genetic differences also epigenetic modifications are considered as potentially relevant in the transfer of susceptibility to psychiatric diseases. Though, whether or not epigenetic alterations can be inherited for many generations is highly controversial. In the present article, we will critically summarize both the genetic findings and the results from epigenetic analyses, including also those of noncoding RNAs. We will argue that one possible solution to the "missing heritability" problem in psychiatry is a potential role of retrotransposons, the exploration of which is presently only in its beginnings.

A genome-wide strategy to identify causes and consequences of retrotransposon expression finds activation by BRCA1 in ovarian cancer

It is challenging to identify the causes and consequences of retrotransposon expression in human disease due to the hundreds of active genomic copies and their poor conservation across species. We profiled genomic insertions of retrotransposons in ovarian cancer. In addition, in ovarian and breast cancer we analyzed RNAs exhibiting Bayesian correlation with retrotransposon RNA to identify causes and consequences of retrotransposon expression. This strategy finds divergent inflammatory responses associated with retrotransposon expression in ovarian and breast cancer and identifies new factors inducing expression of endogenous retrotransposons including anti-viral responses and the common tumor suppressor BRCA1. In cell lines, mouse ovarian epithelial cells and patient-derived tumor spheroids, BRCA1 promotes accumulation of retrotransposon RNA. BRCA1 promotes transcription of active families of retrotransposons and their insertion into the genome. Intriguingly, elevated retrotransposon expression predicts survival in ovarian cancer patients. Retrotransposons are part of a complex regulatory network in ovarian cancer including BRCA1 that contributes to patient survival. The described strategy can be used to identify the regulators and impacts of retrotransposons in various contexts of biology and disease in humans.

HBV-Integration Studies in the Clinic: Role in the Natural History of Infection

Hepatitis B virus (HBV) infection is a major global health problem causing acute and chronic liver disease that can lead to liver cirrhosis and hepatocellular carcinoma (HCC). HBV covalently closed circular DNA (cccDNA) is essential for viral replication and the establishment of a persistent infection. Integrated HBV DNA represents another stable form of viral DNA regularly observed in the livers of infected patients. HBV DNA integration into the host genome occurs early after HBV infection. It is a common occurrence during the HBV life cycle, and it has been detected in all the phases of chronic infection. HBV DNA integration has long been considered to be the main contributor to liver tumorigenesis. The recent development of highly sensitive detection methods and research models has led to the clarification of some molecular and pathogenic aspects of HBV integration. Though HBV integration does not lead to replication-competent transcripts, it can act as a stable source of viral RNA and proteins, which may contribute in determining HBV-specific T-cell exhaustion and favoring virus persistence. The relationship between HBV DNA integration and the immune response in the liver microenvironment might be closely related to the development and progression of HBV-related diseases. While many new antiviral agents aimed at cccDNA elimination or silencing have been developed, integrated HBV DNA remains a difficult therapeutic challenge.

Sense-oriented AluYRa1 elements provide a lineage-specific transcription environment for polyadenylation

Transposable elements cause alternative splicing (AS) in different ways, contributing to transcript diversification. Alternative polyadenylation (APA), one of the AS events, is related to the generation of mRNA isoforms in 70% of human genes. In this study, we tried to investigate AluYRa1s located at the terminal region of cynomolgus monkey genes, utilizing both computational analysis and molecular experimentation. We found that ten genes had AluYRa1 at their 3' end, and nine of these AluYRa1s were sense-oriented. Furthermore, in seven genes, AluYRa1s were expected to have a similar consensus sequence for polyadenylation cleavage. Additional computational analysis using the annotation files from the UCSC database showed that AluYRa1 was more involved in polyadenylation than in open reading frame exon splicing. To examine the extent of AluYRa1 involvement in polyadenylation, RNA-seq data from 30 normal cynomolgus monkeys were analyzed using TAPAS, a recently devised software that detects all the promising polyadenylation sites including APA sites. We observed that approximately 74% of possible polyadenylation sites in the analyzed genes were provided by sense-oriented AluYRa1. In conclusion, AluYRa1 is an Old-World monkey-specific TE, and its sense-oriented insertion at the 3'UTR region tends to provide a favorable environment for polyadenylation, diversifying gene transcripts.

Cytoplasmic synthesis of endogenous Alu complementary DNA via reverse transcription and implications in age-related macular degeneration

Alu retroelements propagate via retrotransposition by hijacking long interspersed nuclear element-1 (L1) reverse transcriptase (RT) and endonuclease activities. Reverse transcription of Alu RNA into complementary DNA (cDNA) is presumed to occur exclusively in the nucleus at the genomic integration site. Whether Alu cDNA is synthesized independently of genomic integration is unknown. Alu RNA promotes retinal pigmented epithelium (RPE) death in geographic atrophy, an untreatable type of age-related macular degeneration. We report that Alu RNA-induced RPE degeneration is mediated via cytoplasmic L1-reverse-transcribed Alu cDNA independently of retrotransposition. Alu RNA did not induce cDNA production or RPE degeneration in L1-inhibited animals or human cells. Alu reverse transcription can be initiated in the cytoplasm via self-priming of Alu RNA. In four health insurance databases, use of nucleoside RT inhibitors was associated with reduced risk of developing atrophic macular degeneration (pooled adjusted hazard ratio, 0.616; 95% confidence interval, 0.493-0.770), thus identifying inhibitors of this Alu replication cycle shunt as potential therapies for a major cause of blindness.

Role of Transposable Elements in Gene Regulation in the Human Genome

Transposable elements (TEs), also known as mobile elements (MEs), are interspersed repeats that constitute a major fraction of the genomes of higher organisms. As one of their important functional impacts on gene function and genome evolution, TEs participate in regulating the expression of genes nearby and even far away at transcriptional and post-transcriptional levels. There are two known principal ways by which TEs regulate the expression of genes. First, TEs provide cis-regulatory sequences in the genome with their intrinsic regulatory properties for their own expression, making them potential factors for regulating the expression of the host genes. TE-derived cis-regulatory sites are found in promoter and enhancer elements, providing binding sites for a wide range of trans-acting factors. Second, TEs encode for regulatory RNAs with their sequences showed to be present in a substantial fraction of miRNAs and long non-coding RNAs (lncRNAs), indicating the TE origin of these RNAs. Furthermore, TEs sequences were found to be critical for regulatory functions of these RNAs, including binding to the target mRNA. TEs thus provide crucial regulatory roles by being part of cis-regulatory and regulatory RNA sequences. Moreover, both TE-derived cis-regulatory sequences and TE-derived regulatory RNAs have been implicated in providing evolutionary novelty to gene regulation. These TE-derived regulatory mechanisms also tend to function in a tissue-specific fashion. In this review, we aim to comprehensively cover the studies regarding these two aspects of TE-mediated gene regulation, mainly focusing on the mechanisms, contribution of different types of TEs, differential roles among tissue types, and lineage-specificity, based on data mostly in humans.

Sea Urchin as a Universal Model for Studies of Gene Networks

The purple sea urchin Strongylocentrotus purpuratus has been used for over 150 years as a model organism in developmental biology. Using this model species, scientists have been able to describe, in detail, the mechanisms of cell cycle control and cell adhesion, fertilization, calcium signaling, cell differentiation, and death. Massive parallel sequencing of the sea urchin genome enabled the deciphering of the main components of gene regulatory networks during the activation of embryonic signaling pathways. This knowledge helped to extrapolate aberrations in somatic cells that may lead to diseases, including cancer in humans. Furthermore, since many, if not all, developmental signaling pathways were shown to be controlled by non-coding RNAs (ncRNAs), the sea urchin organism represents an attractive experimental model. In this review, we discuss the main discoveries in the genetics, genomics, and transcriptomics of sea urchins during embryogenesis with the main focus on the role of ncRNAs. This information may be useful for comparative studies between different organisms, and may help identify new regulatory networks controlled by ncRNAs.

A multi-omics view of the complex mechanism of vascular calcification

Vascular calcification is a high incidence and high risk disease with increasing morbidity and high mortality, which is considered the consequence of smooth muscle cell transdifferentiation initiating the mechanism of accumulation of hydroxyl calcium phosphate. Vascular calcification is also thought to be strongly associated with poor outcomes in diabetes and chronic kidney disease. Numerous studies have been accomplished; however, the specific mechanism of the disease remains unclear. Development of the genome project enhanced the understanding of life science and has entered the post-genomic era resulting in a variety of omics techniques used in studies and a large amount of available data; thus, a new perspective on data analysis has been revealed. Omics has a broader perspective and is thus advantageous over a single pathway analysis in the study of complex vascular calcification mechanisms. This paper reviews in detail various omics studies including genomics, proteomics, transcriptomics, metabolomics and multiple group studies on vascular calcification. Advances and deficiencies in the use of omics to study vascular calcification are presented in a comprehensive view. We also review the methodology of the omics studies and omics data analysis and processing. In addition, the methodology and data processing presented here can be applied to other areas. An omics landscape perspective across the boundaries between genomics, transcriptomics, proteomics and metabolomics is used to examine the mechanisms of vascular calcification. The perspective combined with various technologies also provides a direction for the subsequent exploration of clinical significance.

SeqURE - a new copy-capture based method for sequencing of unknown Retroposition events

Background

Retroelements (REs) occupy a significant part of all eukaryotic genomes including humans. The majority of retroelements in the human genome are inactive and unable to retrotranspose. Dozens of active copies are repressed in most normal tissues by various cellular mechanisms. These copies can become active in normal germline and brain tissues or in cancer, leading to new retroposition events. The consequences of such events and their role in normal cell functioning and carcinogenesis are not yet fully understood. If new insertions occur in a small portion of cells they can be found only with the use of specific methods based on RE enrichment and high-throughput sequencing. The downside of the high sensitivity of such methods is the presence of various artifacts imitating real insertions, which in many cases cannot be validated due to lack of the initial template DNA. For this reason, adequate assessment of rare (< 1%) subclonal cancer specific RE insertions is complicated.

Results

Here we describe a new copy-capture technique which we implemented in a method called SeqURE for Sequencing Unknown of Retroposition Events that allows for efficient and reliable identification of new genomic RE insertions. The method is based on the capture of copies of target molecules (copy-capture), selective amplification and sequencing of genomic regions adjacent to active RE insertions from both sides. Importantly, the template genomic DNA remains intact and can be used for validation experiments. In addition, we applied a novel system for testing method sensitivity and precisely showed the ability of the developed method to reliably detect insertions present in 1 out of 100 cells and a substantial portion of insertions present in 1 out of 1000 cells. Using advantages of the method we showed the absence of somatic Alu insertions in colorectal cancer samples bearing tumor-specific L1HS insertions.

Conclusions

This study presents the first description and implementation of the copy-capture technique and provides the first methodological basis for the quantitative assessment of RE insertions present in a small portion of cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-020-00228-6.

SARS-CoV-2 RNA reverse-transcribed and integrated into the human genome

Prolonged SARS-CoV-2 RNA shedding and recurrence of PCR-positive tests have been widely reported in patients after recovery, yet these patients most commonly are non-infectious. Here we investigated the possibility that SARS-CoV-2 RNAs can be reverse-transcribed and integrated into the human genome and that transcription of the integrated sequences might account for PCR-positive tests. In support of this hypothesis, we found chimeric transcripts consisting of viral fused to cellular sequences in published data sets of SARS-CoV-2 infected cultured cells and primary cells of patients, consistent with the transcription of viral sequences integrated into the genome. To experimentally corroborate the possibility of viral retro-integration, we describe evidence that SARS-CoV-2 RNAs can be reverse transcribed in human cells by reverse transcriptase (RT) from LINE-1 elements or by HIV-1 RT, and that these DNA sequences can be integrated into the cell genome and subsequently be transcribed. Human endogenous LINE-1 expression was induced upon SARS-CoV-2 infection or by cytokine exposure in cultured cells, suggesting a molecular mechanism for SARS-CoV-2 retro-integration in patients. This novel feature of SARS-CoV-2 infection may explain why patients can continue to produce viral RNA after recovery and suggests a new aspect of RNA virus replication.

SVA retrotransposon insertion in exon of MMR genes results in aberrant RNA splicing and causes Lynch syndrome

Lynch syndrome is an autosomal dominant hereditary cancer syndrome in which many cancers develop, the main one being colorectal cancer. Germline pathogenic variants in one of four mismatch repair (MMR) genes are known to be causative of this disease. Accurate diagnosis using genetic testing can greatly benefit the health of those affected. Recently, owing to the improvement of sequence techniques, complicated variants affecting the functions of MMR genes were discovered. In this study, we analyzed insertions of a retrotransposon-like sequence in exon 5 of the MSH6 gene and exon 3 of the MSH2 gene found in Japanese families suspected of having Lynch syndrome. Both of these insertions induced aberrant splicing, and these variants were successfully identified by mRNA sequencing or visual observation of mapping results, although a standard DNA-seq analysis pipeline failed to detect them. The insertion sequences were ~2.5 kbp in length and were found to have the structure of an SVA retrotransposon (SVA). One SVA sequence was not present in the hg19 or hg38 reference genome, but was in a Japanese-specific reference sequence (JRGv2). Our study illustrates the difficulties of identifying SVA insertions in disease genes, and that the possibility of polymorphic insertions should be considered when analyzing mobile elements.

Frequency and methylation status of selected retrotransposition competent L1 loci in amyotrophic lateral sclerosis

Long interspersed element-1 (LINE-1/L1) is the only autonomous transposable element in the human genome that currently mobilises in both germline and somatic tissues. Recent studies have identified correlations between altered retrotransposon expression and the fatal neurodegenerative disease amyotrophic lateral sclerosis (ALS) in a subset of patients. The risk of an individual developing ALS is dependent on an interaction of genetic variants and subsequent modifiers during life. These modifiers could include environmental factors, which can lead to epigenetic and genomic changes, such as somatic mutations, occurring in the neuronal cells that degenerate as the disease develops. There are more than 1 million L1 copies in the human genome today, but only 80-100 L1 loci in the reference genome are considered to be retrotransposition-competent (RC) and an even smaller number of these RC-L1s loci are highly active. We hypothesise that RC-L1s could affect normal cellular function through their mutagenic potential conferred by their ability to retrotranspose in neuronal cells and through DNA damage caused by the endonuclease activity of the L1-encoded ORF2 protein. To investigate whether either an increase in the genomic burden of RC-L1s or epigenetic changes to RC-L1s altering their expression, could play a role in disease development, we chose a set of seven well characterised genomic RC-L1 loci that were reported earlier to be highly active in a cellular L1 retrotransposition reporter assay or serve as major source elements for germline and/or somatic retrotransposition events. Analysis of the insertion allele frequency of five polymorphic RC-L1s, out of the set of seven, for their presence or absence, did not identify an increased number individually or when combined in individuals with the disease. However, we did identify reduced levels of methylation of RC-L1s in the motor cortex of those individuals with both familial and sporadic ALS compared to control brains. The changes to the regulation of the loci encompassing these RC-L1s demonstrated tissue specificity and could be related to the disease process.

Nanopore Sequencing Enables Comprehensive Transposable Element Epigenomic Profiling

Transposable elements (TEs) drive genome evolution and are a notable source of pathogenesis, including cancer. While CpG methylation regulates TE activity, the locus-specific methylation landscape of mobile human TEs has to date proven largely inaccessible. Here, we apply new computational tools and long-read nanopore sequencing to directly infer CpG methylation of novel and extant TE insertions in hippocampus, heart, and liver, as well as paired tumor and non-tumor liver. As opposed to an indiscriminate stochastic process, we find pronounced demethylation of young long interspersed element 1 (LINE-1) retrotransposons in cancer, often distinct to the adjacent genome and other TEs. SINE-VNTR-Alu (SVA) retrotransposons, including their internal tandem repeat-associated CpG island, are near-universally methylated. We encounter allele-specific TE methylation and demethylation of aberrantly expressed young LINE-1s in normal tissues. Finally, we recover the complete sequences of tumor-specific LINE-1 insertions and their retrotransposition hallmarks, demonstrating how long-read sequencing can simultaneously survey the epigenome and detect somatic TE mobilization.

Functional Long Non-coding RNAs Evolve from Junk Transcripts

Transcriptome studies reveal pervasive transcription of complex genomes, such as those of mammals. Despite popular arguments for functionality of most, if not all, of these transcripts, genome-wide analysis of selective constraints indicates that most of the produced RNA are junk. However, junk is not garbage. On the contrary, junk transcripts provide the raw material for the evolution of diverse long non-coding (lnc) RNAs by non-adaptive mechanisms, such as constructive neutral evolution. The generation of many novel functional entities, such as lncRNAs, that fuels organismal complexity does not seem to be driven by strong positive selection. Rather, the weak selection regime that dominates the evolution of most multicellular eukaryotes provides ample material for functional innovation with relatively little adaptation involved.

Replication protein A: a multifunctional protein with roles in DNA replication, repair and beyond

Single-stranded DNA (ssDNA) forms continuously during DNA replication and is an important intermediate during recombination-mediated repair of damaged DNA. Replication protein A (RPA) is the major eukaryotic ssDNA-binding protein. As such, RPA protects the transiently formed ssDNA from nucleolytic degradation and serves as a physical platform for the recruitment of DNA damage response factors. Prominent and well-studied RPA-interacting partners are the tumor suppressor protein p53, the RAD51 recombinase and the ATR-interacting proteins ATRIP and ETAA1. RPA interactions are also documented with the helicases BLM, WRN and SMARCAL1/HARP, as well as the nucleotide excision repair proteins XPA, XPG and XPF–ERCC1. Besides its well-studied roles in DNA replication (restart) and repair, accumulating evidence shows that RPA is engaged in DNA activities in a broader biological context, including nucleosome assembly on nascent chromatin, regulation of gene expression, telomere maintenance and numerous other aspects of nucleic acid metabolism. In addition, novel RPA inhibitors show promising effects in cancer treatment, as single agents or in combination with chemotherapeutics. Since the biochemical properties of RPA and its roles in DNA repair have been extensively reviewed, here we focus on recent discoveries describing several non-canonical functions.

Transposon expression in the Drosophila brain is driven by neighboring genes and diversifies the neural transcriptome

Somatic transposon expression in neural tissue is commonly considered as a measure of mobilization and has therefore been linked to neuropathology and organismal individuality. We combined genome sequencing data with single-cell mRNA sequencing of the same inbred fly strain to map transposon expression in the Drosophila midbrain and found that transposon expression patterns are highly stereotyped. Every detected transposon is resident in at least one cellular gene with a matching expression pattern. Bulk RNA sequencing from fly heads of the same strain revealed that coexpression is a physical link in the form of abundant chimeric transposon-gene mRNAs. We identified 264 genes where transposons introduce cryptic splice sites into the nascent transcript and thereby significantly expand the neural transcript repertoire. Some genes exclusively produce chimeric mRNAs with transposon sequence; on average, 11.6% of the mRNAs produced from a given gene are chimeric. Conversely, most transposon-containing transcripts are chimeric, which suggests that somatic expression of these transposons is largely driven by cellular genes. We propose that chimeric mRNAs produced by alternative splicing into polymorphic transposons, rather than transposon mobilization, may contribute to functional differences between individual cells and animals.

Repurposing anti-inflammasome NRTIs for improving insulin sensitivity and reducing type 2 diabetes development

Innate immune signaling through the NLRP3 inflammasome is activated by multiple diabetes-related stressors, but whether targeting the inflammasome is beneficial for diabetes is still unclear. Nucleoside reverse-transcriptase inhibitors (NRTI), drugs approved to treat HIV-1 and hepatitis B infections, also block inflammasome activation. Here, we show, by analyzing five health insurance databases, that the adjusted risk of incident diabetes is 33% lower in patients with NRTI exposure among 128,861 patients with HIV-1 or hepatitis B (adjusted hazard ratio for NRTI exposure, 0.673; 95% confidence interval, 0.638 to 0.710; P < 0.0001; 95% prediction interval, 0.618 to 0.734). Meanwhile, an NRTI, lamivudine, improves insulin sensitivity and reduces inflammasome activation in diabetic and insulin resistance-induced human cells, as well as in mice fed with high-fat chow; mechanistically, inflammasome-activating short interspersed nuclear element (SINE) transcripts are elevated, whereas SINE-catabolizing DICER1 is reduced, in diabetic cells and mice. These data suggest the possibility of repurposing an approved class of drugs for prevention of diabetes.

Mosaicism in Human Health and Disease

Mosaicism refers to the occurrence of two or more genomes in an individual derived from a single zygote. Germline mosaicism is a mutation that is limited to the gonads and can be transmitted to offspring. Somatic mosaicism is a postzygotic mutation that occurs in the soma, and it may occur at any developmental stage or in adult tissues. Mosaic variation may be classified in six ways: (a) germline or somatic origin, (b) class of DNA mutation (ranging in scale from single base pairs to multiple chromosomes), (c) developmental context, (d) body location(s), (e) functional consequence (including deleterious, neutral, or advantageous), and (f) additional sources of mosaicism, including mitochondrial heteroplasmy, exogenous DNA sources such as vectors, and epigenetic changes such as imprinting and X-chromosome inactivation. Technological advances, including single-cell and other next-generation sequencing, have facilitated improved sensitivity and specificity to detect mosaicism in a variety of biological contexts.

Promoter methylation, transcription, and retrotransposition of LINE-1 in colorectal adenomas and adenocarcinomas

BACKGROUND

The methylation of the CpG islands of the LINE-1 promoter is a tight control mechanism on the function of mobile elements. However, simultaneous quantification of promoter methylation and transcription of LINE-1 has not been performed in progressive stages of colorectal cancer. In addition, the insertion of mobile elements in the genome of advanced adenoma stage, a precancerous stage before colorectal carcinoma has not been emphasized. In this study, we quantify promoter methylation and transcripts of LINE-1 in three stages of colorectal non-advanced adenoma, advanced adenoma, and adenocarcinoma. In addition, we analyze the insertion of LINE-1, Alu, and SVA elements in the genome of patient tumors with colorectal advanced adenomas.

METHODS

LINE-1 hypomethylation status was evaluated by absolute quantitative analysis of methylated alleles (AQAMA) assay. To quantify the level of transcripts for LINE-1, quantitative RT-PCR was performed. To find mobile element insertions, the advanced adenoma tissue samples were subjected to whole genome sequencing and MELT analysis.

RESULTS

We found that the LINE-1 promoter methylation in advanced adenoma and adenocarcinoma was significantly lower than that in non-advanced adenomas. Accordingly, the copy number of LINE-1 transcripts in advanced adenoma was significantly higher than that in non-advanced adenomas, and in adenocarcinomas was significantly higher than that in the advanced adenomas. Whole-genome sequencing analysis of colorectal advanced adenomas revealed that at this stage polymorphic insertions of LINE-1, Alu, and SVA comprise approximately 16%, 51%, and 74% of total insertions, respectively.

CONCLUSIONS

Our correlative analysis showing a decreased methylation of LINE-1 promoter accompanied by the higher level of LINE-1 transcription, and polymorphic genomic insertions in advanced adenoma, suggests that the early and advanced polyp stages may host very important pathogenic processes concluding to cancer.

Early-onset breast cancer in a woman with a germline mobile element insertion resulting in BRCA2 disruption: a case report

Mobile element insertions (MEIs) contribute to genomic diversity, but they can be responsible for human disease in some cases. Initial clinical testing (BRCA1, BRCA2 and PALB2) in a 40-year-old female with unilateral breast cancer did not detect any pathogenic variants. Subsequent reanalysis for MEIs detected a novel likely pathogenic insertion of the retrotransposon element (RE) c.7894_7895insSVA in BRCA2. This case highlights the importance of bioinformatic pipeline optimization for the detection of MEIs in genes associated with hereditary cancer, as early detection can significantly impact clinical management.

What impact does oocyte vitrification have on epigenetics and gene expression?

Children conceived by assisted reproductive technologies (ART) have a moderate risk for a number of adverse events and conditions. The question whether this additional risk is associated with specific procedures used in ART or whether it is related to the intrinsic biological factors associated with infertility remains unresolved. One of the main hypotheses is that laboratory procedures could have an effect on the epigenome of gametes and embryos. This suspicion is linked to the fact that ART procedures occur precisely during the period when there are major changes in the organization of the epigenome. Oocyte freezing protocols are generally considered safe; however, some evidence suggests that vitrification may be associated with modifications of the epigenetic marks. In this manuscript, after describing the main changes that occur during epigenetic reprogramming, we will provide current information regarding the impact of oocyte vitrification on epigenetic regulation and the consequences on gene expression, both in animals and humans. Overall, the literature suggests that epigenetic and transcriptomic profiles are sensitive to the stress induced by oocyte vitrification, and it also underlines the need to improve our knowledge in this field.

New Understanding of the Relevant Role of LINE-1 Retrotransposition in Human Disease and Immune Modulation

Long interspersed nuclear element-1 (LINE-1) retrotransposition is a major hallmark of cancer accompanied by global chromosomal instability, genomic instability, and genetic heterogeneity and has become one indicator for the occurrence, development, and poor prognosis of many diseases. LINE-1 also modulates the immune system and affects the immune microenvironment in a variety of ways. Aberrant expression of LINE-1 retrotransposon can provide strong stimuli for an innate immune response, activate the immune system, and induce autoimmunity and inflammation. Therefore, inhibition the activity of LINE-1 has become a potential treatment strategy for various diseases. In this review, we discussed the components and regulatory mechanisms involved with LINE-1, its correlations with disease and immunity, and multiple inhibitors of LINE-1, providing a new understanding of LINE-1.

Polymorphic mobile element insertions contribute to gene expression and alternative splicing in human tissues

BACKGROUND

Mobile elements are a major source of structural variants in the human genome, and some mobile elements can regulate gene expression and transcript splicing. However, the impact of polymorphic mobile element insertions (pMEIs) on gene expression and splicing in diverse human tissues has not been thoroughly studied. The multi-tissue gene expression and whole genome sequencing data generated by the Genotype-Tissue Expression (GTEx) project provide a great opportunity to systematically evaluate the role of pMEIs in regulating gene expression in human tissues.

RESULTS

Using the GTEx whole genome sequencing data, we identify 20,545 high-quality pMEIs from 639 individuals. Coupling pMEI genotypes with gene expression profiles, we identify pMEI-associated expression quantitative trait loci (eQTLs) and splicing quantitative trait loci (sQTLs) in 48 tissues. Using joint analyses of pMEIs and other genomic variants, pMEIs are predicted to be the potential causal variant for 3522 eQTLs and 3717 sQTLs. The pMEI-associated eQTLs and sQTLs show a high level of tissue specificity, and these pMEIs are enriched in the proximity of affected genes and in regulatory elements. Using reporter assays, we confirm that several pMEIs associated with eQTLs and sQTLs can alter gene expression levels and isoform proportions, respectively.

CONCLUSION

Overall, our study shows that pMEIs are associated with thousands of gene expression and splicing variations, indicating that pMEIs could have a significant role in regulating tissue-specific gene expression and transcript splicing. Detailed mechanisms for the role of pMEIs in gene regulation in different tissues will be an important direction for future studies.

Human TRIM5α senses and restricts LINE-1 elements

Mobile genetic elements have significantly shaped our genomic landscape. LINE-1 retroelements are the only autonomously active elements left in the human genome. Since new insertions can have detrimental consequences, cells need to efficiently control LINE-1 retrotransposition. Here, we demonstrate that the intrinsic immune factor TRIM5α senses and restricts LINE-1 retroelements. Previously, rhesus TRIM5α has been shown to efficiently block HIV-1 replication, while human TRIM5α was found to be less active. Surprisingly, we found that both human and rhesus TRIM5α efficiently repress human LINE-1 retrotransposition. TRIM5α interacts with LINE-1 ribonucleoprotein complexes in the cytoplasm, which is essential for restriction. In line with its postulated role as pattern recognition receptor, we show that TRIM5α also induces innate immune signaling upon interaction with LINE-1 ribonucleoprotein complexes. The signaling events activate the transcription factors AP-1 and NF-κB, leading to the down-regulation of LINE-1 promoter activity. Together, our findings identify LINE-1 as important target of human TRIM5α, which restricts and senses LINE-1 via two distinct mechanisms. Our results corroborate TRIM5α as pattern recognition receptor and shed light on its previously undescribed activity against mobile genetic elements, such as LINE-1, to protect the integrity of our genome.

ERVcaller: identifying polymorphic endogenous retrovirus and other transposable element insertions using whole-genome sequencing data

MOTIVATION

Approximately 8% of the human genome is derived from endogenous retroviruses (ERVs). In recent years, an increasing number of human diseases have been found to be associated with ERVs. However, it remains challenging to accurately detect the full spectrum of polymorphic (unfixed) ERVs using whole-genome sequencing (WGS) data.

RESULTS

We designed a new tool, ERVcaller, to detect and genotype transposable element (TE) insertions, including ERVs, in the human genome. We evaluated ERVcaller using both simulated and real benchmark WGS datasets. Compared to existing tools, ERVcaller consistently obtained both the highest sensitivity and precision for detecting simulated ERV and other TE insertions derived from real polymorphic TE sequences. For the WGS data from the 1000 Genomes Project, ERVcaller detected the largest number of TE insertions per sample based on consensus TE loci. By analyzing the experimentally verified TE insertions, ERVcaller had 94.0% TE detection sensitivity and 96.6% genotyping accuracy. Polymerase chain reaction and Sanger sequencing in a small sample set verified 86.7% of examined insertion statuses and 100% of examined genotypes. In conclusion, ERVcaller is capable of detecting and genotyping TE insertions using WGS data with both high sensitivity and precision. This tool can be applied broadly to other species.

AVAILABILITY AND IMPLEMENTATION

http://www.uvm.edu/genomics/software/ERVcaller.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Measuring and interpreting transposable element expression

Transposable elements (TEs) are insertional mutagens that contribute greatly to the plasticity of eukaryotic genomes, influencing the evolution and adaptation of species as well as physiology or disease in individuals. Measuring TE expression helps to understand not only when and where TE mobilization can occur but also how this process alters gene expression, chromatin accessibility or cellular signalling pathways. Although genome-wide gene expression assays such as RNA sequencing include transposon-derived transcripts, most computational analytical tools discard or misinterpret TE-derived reads. Emerging approaches are improving the identification of expressed TE loci and helping to discriminate TE transcripts that permit TE mobilization from chimeric gene-TE transcripts or pervasive transcription. Here we review the main challenges associated with the detection of TE expression, including mappability, insertional and internal sequence polymorphisms, and the diversity of the TE transcriptional landscape, as well as the different experimental and computational strategies to solve them.

Piwi reduction in the aged niche eliminates germline stem cells via Toll-GSK3 signaling

Transposons are known to participate in tissue aging, but their effects on aged stem cells remain unclear. Here, we report that in the Drosophila ovarian germline stem cell (GSC) niche, aging-related reductions in expression of Piwi (a transposon silencer) derepress retrotransposons and cause GSC loss. Suppression of Piwi expression in the young niche mimics the aged niche, causing retrotransposon depression and coincident activation of Toll-mediated signaling, which promotes Glycogen synthase kinase 3 activity to degrade β-catenin. Disruption of β-catenin-E-cadherin-mediated GSC anchorage then results in GSC loss. Knocking down gypsy (a highly active retrotransposon) or toll, or inhibiting reverse transcription in the piwi-deficient niche, suppresses GSK3 activity and β-catenin degradation, restoring GSC-niche attachment. This retrotransposon-mediated impairment of aged stem cell maintenance may have relevance in many tissues, and could represent a viable therapeutic target for aging-related tissue degeneration.

Stress, Adaptation, and the Deep Genome: Why Transposons Matter

Stress is a common, if often unpredictable life event. It can be defined from an evolutionary perspective as a force an organism perceives it must adapt to. Thus stress is a useful tool to study adaptation and the adaptive capacity of organisms. The deep genome, long neglected as a pile of “junk” has emerged as a source of regulatory DNA and RNA as well as a potential stockpile of adaptive capacity at the organismal and species levels. Recent work on the regulation of transposable elements (TEs), the principle constituents of the deep genome, by stress has shown that these elements are responsive to host stress and other environmental cues. Further, we have shown that some are likely directly regulated by the glucocorticoid receptor (GR), one of the two major vertebrate stress steroid receptors in a fashion that appears adaptive. On the basis of this and other emerging evidence I argue that the deep genome may represent an adaptive toolkit for organisms to respond to their environments at both individual and evolutionary scales. This argues that genomes may be adapted for what Waddington called “trait adaptability” rather than being purely passive objects of natural selection and single nucleotide level mutation.

PiRNAs Rise to Rescue Koalas

PIWI-interacting small RNAs (piRNAs) establish sequence-specific adaptive restriction of resident genomic parasites to guard genome integrity. In this issue of Cell, Yu, Koppetsch, et al. describe an innate piRNA-response that specifically fragments the viral RNA genome in the germline of recently invaded koalas. This first line of defense might ensure survival until adaptive immunity develops.

Alu-Mediated MEN1 Gene Deletion and Loss of Heterozygosity in a Patient with Multiple Endocrine Neoplasia Type 1

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant disorder caused by mutations of the tumor suppressor gene MEN1. Most of the germline MEN1 gene mutations have been small mutations, and the whole gene deletion is rarely observed. In the present study, we revealed Alu retrotransposon-mediated de novo germline deletion of the whole MEN1 gene and somatic copy-neutral loss of heterozygosity (LOH) in a patient with MEN1. The patient is a 39-year-old woman who was referred to our department for the management of prolactinoma. She was also diagnosed with primary hyperparathyroidism and suspected of MEN1. Although nucleotide sequencing did not detect any MEN1 gene mutations, multiplex ligation-dependent probe amplification (MLPA) revealed a large germline deletion of the MEN1 gene. Subsequent quantitative polymerase chain reaction (qPCR)–based copy number mapping showed a monoallelic loss of approximately 18.5-kilobase region containing the whole MEN1 gene. Intriguingly, the 2 breakpoints were flanked by Alu repetitive elements, suggesting the contribution of Alu/Alu-mediated rearrangements (AAMR) to the whole MEN1 gene deletion. Furthermore, copy number mapping using MLPA and qPCR in combination with single nucleotide polymorphism analysis revealed copy-neutral LOH as a somatic event for parathyroid tumorigenesis. In conclusion, copy number mapping revealed a novel combination of Alu/Alu-mediated de novo germline deletion of the MEN1 gene and somatic copy-neutral LOH as a cytogenetic basis for the MEN1 pathogenesis. Moreover, subsequent in silico analysis highlighted the possible predisposition of the MEN1 gene to Alu retrotransposon-mediated genomic deletion.

Established and Novel Mechanisms Leading to de novo Genomic Rearrangements in the Human Germline

During gametogenesis, the human genome can acquire various de novo rearrangements. Most constitutional genomic rearrangements are created through 1 of the 4 well-known mechanisms, i.e., nonallelic homologous recombination, erroneous repair after double-strand DNA breaks, replication errors, and retrotransposition. However, recent studies have identified 2 types of extremely complex rearrangements that cannot be simply explained by these mechanisms. The first type consists of chaotic structural changes in 1 or a few chromosomes that result from “chromoanagenesis (an umbrella term that covers chromothripsis, chromoanasynthesis, and chromoplexy).” The other type is large independent rearrangements in multiple chromosomes indicative of “transient multifocal genomic crisis.” Germline chromoanagenesis (chromothripsis) likely occurs predominantly during spermatogenesis or postzygotic embryogenesis, while multifocal genomic crisis appears to be limited to a specific time window during oogenesis and early embryogenesis or during spermatogenesis. This review article introduces the current understanding of the molecular basis of de novo rearrangements in the germline.

Genomic Indexing by Somatic Gene Recombination of mRNA/ncRNA - Does It Play a Role in Genomic Mosaicism, Memory Formation, and Alzheimer's Disease?

Recent evidence indicates that genomic individuality of neurons, characterized by DNA-content variation, is a common if not universal phenomenon in the human brain that occurs naturally but can also show aberrancies that have been linked to the pathomechanism of Alzheimer’s disease and related neurodegenerative disorders. Etiologically, this genomic mosaic has been suggested to arise from defects of cell cycle regulation that may occur either during brain development or in the mature brain after terminal differentiation of neurons. Here, we aim to draw attention towards another mechanism that can give rise to genomic individuality of neurons, with far-reaching consequences. This mechanism has its origin in the transcriptome rather than in replication defects of the genome, i.e., somatic gene recombination of RNA. We continue to develop the concept that somatic gene recombination of RNA provides a physiological process that, through integration of intronless mRNA/ncRNA into the genome, allows a particular functional state at the level of the individual neuron to be indexed. By insertion of defined RNAs in a somatic recombination process, the presence of specific mRNA transcripts within a definite temporal context can be “frozen” and can serve as an index that can be recalled at any later point in time. This allows information related to a specific neuronal state of differentiation and/or activity relevant to a memory trace to be fixed. We suggest that this process is used throughout the lifetime of each neuron and might have both advantageous and deleterious consequences.

Integrative and Analytical Review of the 5-Alpha-Reductase Type 2 Deficiency Worldwide

Introduction

The conversion of testosterone into dihydrotestosterone is catalyzed by the 5α-reductase type 2 enzyme which plays a crucial role in the external genitalia virilization. It is encoded by the SRD5A2 gene. Allelic variants in this gene cause a 46,XY DSD with no genotype-phenotype relationship. It was firstly reported in the early 70s from isolated clusters. Since then, several cases have been reported. Putting together, it will expand the knowledge on the molecular bases of androgen milieu.

Methods

We searched for SRD5A2 allelic variants (AV) in the literature (PubMed, Embase, MEDLINE) and websites (ensembl, HGMD, ClinVar). Only cases with AV in both alleles, either in homozygous or compound heterozygous were included. The included cases were analyzed according to ethnicity, exon, domain, aminoacid (aa) conservation, age at diagnosis, sex assignment, gender reassignment, external genitalia virilization and functional studies. External genitalia virilization was scored using Sinnecker scale. Conservation analysis was carried out using the CONSURF platform. For categorical variables, we used X2 test and Cramer's V. Continuous variables were analyzed by t test or ANOVA. Concordance was estimated by Kappa.

Results

We identified 434 cases of 5ARD2 deficiencies from 44 countries. Most came from Turkey (23%), China (17%), Italy (9%), and Brazil (7%). Sixty-nine percent were assigned as female. There were 70% of homozygous allelic variants and 30% compound heterozygous. Most were missense variants (76%). However, small indels (11%), splicing (5%) and large deletions (4%) were all reported. They were distributed along with all exons with exon 1 (33%) and exon 4 (25%) predominance. Allelic variants in the exon 4 (NADPH-binding domain) resulted in lower virilization (p<0.0001). The codons 55, 65, 196, 235 and 246 are hotspots making up 25% of all allelic variants. Most of them (76%) were located at conserved aa. However, allelic variants at non-conserved aa were more frequently indels (28% vs 6%; p<0.01). The overall rate of gender change from female to male ranged from 16% to 70%. The lowest rate of gender change from female to male occurred in Turkey and the highest in Brazil. External genitalia virilization was similar between those who changed and those who kept their assigned gender. The gender change rate was significantly different across the countries (V=0.44; p<0.001) even with similar virilization scores.

Conclusion

5ARD2 deficiency has a worldwide distribution. Allelic variants at the NADPH-ligand region cause lower virilization. Genitalia virilization influenced sex assignment but not gender change which was influenced by cultural aspects across the countries. Molecular diagnosis influenced on sex assignment, favoring male sex assignment in newborns with 5α-reductase type 2 deficiency.

Long INterspersed element-1 mobility as a sensor of environmental stresses

Long INterspersed element (LINE-1, L1) retrotransposons are the most abundant transposable elements in the human genome, constituting approximately 17%. They move by a "copy-paste" mechanism, involving reverse transcription of an RNA intermediate and insertion of its cDNA copy at a new site in the genome. L1 retrotransposition (L1-RTP) can cause insertional mutations, alter gene expression, transduce exons, and induce epigenetic dysregulation. L1-RTP is generally repressed; however, a number of observations collected over about 15 years revealed that it can occur in response to environmental stresses. Moreover, emerging evidence indicates that L1-RTP can play a role in the onset of several neurological and oncological diseases in humans. In recent years, great attention has been paid to the exposome paradigm, which proposes that health effects of an environmental factor should be evaluated considering both cumulative environmental exposures and the endogenous processes resulting from the biological response. L1-RTP could be an endogenous process considered for this application. Here, we summarize the current understanding of environmental factors that can affect the retrotransposition of human L1 elements. Evidence indicates that L1-RTP alteration is triggered by numerous and various environmental stressors, such as chemical agents (heavy metals, carcinogens, oxidants, and drugs), physical agents (ionizing and non-ionizing radiations), and experiential factors (voluntary exercise, social isolation, maternal care, and environmental light/dark cycles). These data come from in vitro studies on cell lines and in vivo studies on transgenic animals: future investigations should be focused on physiologically relevant models to gain a better understanding of this topic.

Contribution of unfixed transposable element insertions to human regulatory variation

Thousands of unfixed transposable element (TE) insertions segregate in the human population, but little is known about their impact on genome function. Recently, a few studies associated unfixed TE insertions to mRNA levels of adjacent genes, but the biological significance of these associations, their replicability across cell types and the mechanisms by which they may regulate genes remain largely unknown. Here, we performed a TE-expression QTL analysis of 444 lymphoblastoid cell lines (LCL) and 289 induced pluripotent stem cells using a newly developed set of genotypes for 2743 polymorphic TE insertions. We identified 211 and 176 TE-eQTL acting in cis in each respective cell type. Approximately 18% were shared across cell types with strongly correlated effects. Furthermore, analysis of chromatin accessibility QTL in a subset of the LCL suggests that unfixed TEs often modulate the activity of enhancers and other distal regulatory DNA elements, which tend to lose accessibility when a TE inserts within them. We also document a case of an unfixed TE likely influencing gene expression at the post-transcriptional level. Our study points to broad and diverse cis-regulatory effects of unfixed TEs in the human population and underscores their plausible contribution to phenotypic variation. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

Jump around: transposons in and out of the laboratory

Since Barbara McClintock’s groundbreaking discovery of mobile DNA sequences some 70 years ago, transposable elements have come to be recognized as important mutagenic agents impacting genome composition, genome evolution, and human health. Transposable elements are a major constituent of prokaryotic and eukaryotic genomes, and the transposition mechanisms enabling transposon proliferation over evolutionary time remain engaging topics for study, suggesting complex interactions with the host, both antagonistic and mutualistic. The impact of transposition is profound, as over 100 human heritable diseases have been attributed to transposon insertions. Transposition can be highly mutagenic, perturbing genome integrity and gene expression in a wide range of organisms. This mutagenic potential has been exploited in the laboratory, where transposons have long been utilized for phenotypic screening and the generation of defined mutant libraries. More recently, barcoding applications and methods for RNA-directed transposition are being used towards new phenotypic screens and studies relevant for gene therapy. Thus, transposable elements are significant in affecting biology both in vivo and in the laboratory, and this review will survey advances in understanding the biological role of transposons and relevant laboratory applications of these powerful molecular tools.

Integrated Mobile Element Scanning (ME-Scan) method for identifying multiple types of polymorphic mobile element insertions

Background

Mobile elements are ubiquitous components of mammalian genomes and constitute more than half of the human genome. Polymorphic mobile element insertions (pMEIs) are a major source of human genomic variation and are gaining research interest because of their involvement in gene expression regulation, genome integrity, and disease.

Results

Building on our previous Mobile Element Scanning (ME-Scan) protocols, we developed an integrated ME-Scan protocol to identify three major active families of human mobile elements, AluYb, L1HS, and SVA. This approach selectively amplifies insertion sites of currently active retrotransposons for Illumina sequencing. By pooling the libraries together, we can identify pMEIs from all three mobile element families in one sequencing run. To demonstrate the utility of the new ME-Scan protocol, we sequenced 12 human parent-offspring trios. Our results showed high sensitivity (> 90%) and accuracy (> 95%) of the protocol for identifying pMEIs in the human genome. In addition, we also tested the feasibility of identifying somatic insertions using the protocol.

Conclusions

The integrated ME-Scan protocol is a cost-effective way to identify novel pMEIs in the human genome. In addition, by developing the protocol to detect three mobile element families, we demonstrate the flexibility of the ME-Scan protocol. We present instructions for the library design, a sequencing protocol, and a computational pipeline for downstream analyses as a complete framework that will allow researchers to easily adapt the ME-Scan protocol to their own projects in other genomes.

LINE-1 retrotransposon encoded ORF1p expression and promoter methylation in oral squamous cell carcinoma: a pilot study

Oral squamous cell carcinoma (OSCC) is highly predominant in India due to excessive use of tobacco. Here we investigated Long INterpersed Element 1 (LINE or L1) retrotransposon activity in OSCC samples in the same population. There are almost 500,000 copies of L1 occupied around 30%  of the human genome. Although most of them are inactive, around 150-200 copies are actively jumping in a human genome. L1 encodes two proteins designated as ORF1p and ORF2p and expression of both proteins are critical for the process of retrotransposition. Here we have analyzed L1 ORF1p expression in a small cohort (n = 15) of paired cancer-normal tissues obtained from operated oral cancer patients. Immunohistochemistry (IHC) with the human ORF1 antibody showed the presence of ORF1p in almost 60%  cancer samples, and few of them also showed aberrant p53 expression.  Investigating L1 promoter methylation status, showed certain trends towards hypomethylation of the L1 promoter in cancer tissues compared to its normal counterpart. Our data raise the possibility that L1ORF1p expression might have some role in the onset and progression of this particular type of cancer.

Short interspersed nuclear element (SINE)-mediated post-transcriptional effects on human and mouse gene expression: SINE-UP for active duty

Primate-specific Alu short interspersed nuclear elements (SINEs) and rodent-specific B and ID (B/ID) SINEs are non-autonomous and generally non-coding retrotransposons that have been copied and pasted into the respective genomes so as to constitute what is estimated to be a remarkable 13% and 8% of those genomes. In the context of messenger RNAs (mRNAs), those residing within 3'-untranslated regions (3'UTRs) can influence mRNA export from the nucleus to the cytoplasm, mRNA translation and/or mRNA decay via proteins with which they associate either individually or base-paired in cis or in trans with a partially complementary SINE. Each of these influences impinges on the primary function of mRNA, which is to serve as a template for protein synthesis. This review describes how human cells have used 3'UTR Alu elements to mediate post-transcriptional gene regulation and also describes examples of convergent evolution between human and mouse 3'UTR SINEs. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Human transposon insertion profiling by sequencing (TIPseq) to map LINE-1 insertions in single cells

Long interspersed element-1 (LINE-1, L1) sequences, which comprise about 17% of human genome, are the product of one of the most active types of mobile DNAs in modern humans. LINE-1 insertion alleles can cause inherited and de novo genetic diseases, and LINE-1-encoded proteins are highly expressed in some cancers. Genome-wide LINE-1 mapping in single cells could be useful for defining somatic and germline retrotransposition rates, and for enabling studies to characterize tumour heterogeneity, relate insertions to transcriptional and epigenetic effects at the cellular level, or describe cellular phylogenies in development. Our laboratories have reported a genome-wide LINE-1 insertion site mapping method for bulk DNA, named transposon insertion profiling by sequencing (TIPseq). There have been significant barriers applying LINE-1 mapping to single cells, owing to the chimeric artefacts and features of repetitive sequences. Here, we optimize a modified TIPseq protocol and show its utility for LINE-1 mapping in single lymphoblastoid cells. Results from single-cell TIPseq experiments compare well to known LINE-1 insertions found by whole-genome sequencing and TIPseq on bulk DNA. Among the several approaches we tested, whole-genome amplification by multiple displacement amplification followed by restriction enzyme digestion, vectorette ligation and LINE-1-targeted PCR had the best assay performance. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

The Developmental Gene Hypothesis for Punctuated Equilibrium: Combined Roles of Developmental Regulatory Genes and Transposable Elements

Theories of the genetics underlying punctuated equilibrium (PE) have been vague to date. Here the developmental gene hypothesis is proposed, which states that: 1) developmental regulatory (DevReg) genes are responsible for the orchestration of metazoan morphogenesis and their extreme conservation and mutation intolerance generates the equilibrium or stasis present throughout much of the fossil record and 2) the accumulation of regulatory elements and recombination within these same genes-often derived from transposable elements-drives punctuated bursts of morphological divergence and speciation across metazoa. This two-part hypothesis helps to explain the features that characterize PE, providing a theoretical genetic basis for the once-controversial theory. Also see the video abstract here https://youtu.be/C-fu-ks5yDs.

Maximizing the ovarian reserve in mice by evading LINE-1 genotoxicity

Female reproductive success critically depends on the size and quality of a finite ovarian reserve. Paradoxically, mammals eliminate up to 80% of the initial oocyte pool through the enigmatic process of fetal oocyte attrition (FOA). Here, we interrogate the striking correlation of FOA with retrotransposon LINE-1 (L1) expression in mice to understand how L1 activity influences FOA and its biological relevance. We report that L1 activity triggers FOA through DNA damage-driven apoptosis and the complement system of immunity. We demonstrate this by combined inhibition of L1 reverse transcriptase activity and the Chk2-dependent DNA damage checkpoint to prevent FOA. Remarkably, reverse transcriptase inhibitor AZT-treated Chk2 mutant oocytes that evade FOA initially accumulate, but subsequently resolve, L1-instigated genotoxic threats independent of piRNAs and differentiate, resulting in an increased functional ovarian reserve. We conclude that FOA serves as quality control for oocyte genome integrity, and is not obligatory for oogenesis nor fertility.

LINE-1 ORF2p expression is nearly imperceptible in human cancers

BACKGROUND

Long interspersed element-1 (LINE-1, L1) is the major driver of mobile DNA activity in modern humans. When expressed, LINE-1 loci produce bicistronic transcripts encoding two proteins essential for retrotransposition, ORF1p and ORF2p. Many types of human cancers are characterized by L1 promoter hypomethylation, L1 transcription, L1 ORF1p protein expression, and somatic L1 retrotransposition. ORF2p encodes the endonuclease and reverse transcriptase activities required for L1 retrotransposition. Its expression is poorly characterized in human tissues and cell lines.

RESULTS

We report mass spectrometry-based tumor proteome profiling studies wherein ORF2p eludes detection. To test whether ORF2p could be detected with specific reagents, we developed and validated five rabbit monoclonal antibodies with immunoreactivity for specific epitopes on the protein. These reagents readily detect ectopic ORF2p expressed from bicistronic L1 constructs. However, endogenous ORF2p is not detected in human tumor samples or cell lines by western blot, immunoprecipitation, or immunohistochemistry despite high levels of ORF1p expression. Moreover, we report endogenous ORF1p-associated interactomes, affinity isolated from colorectal cancers, wherein we similarly fail to detect ORF2p. These samples include primary tumors harboring hundreds of somatically acquired L1 insertions. The new data are available via ProteomeXchange with identifier PXD013743.

CONCLUSIONS

Although somatic retrotransposition provides unequivocal genetic evidence for the expression of ORF2p in human cancers, we are unable to directly measure its presence using several standard methods. Experimental systems have previously indicated an unequal stoichiometry between ORF1p and ORF2p, but in vivo, the expression of these two proteins may be more strikingly uncoupled. These findings are consistent with observations that ORF2p is not tolerable for cell growth.