A LINE-1-encoded reverse transcriptase-dependent regulatory mechanism is active in embryogenesis and tumorigenesis

Transposon insertion profiling by sequencing (TIPseq) identifies novel LINE-1 insertions in human sperm

PURPOSE

Long interspersed nuclear element-1 (LINE-1 or L1) comprises 17% of the human genome. Retrotransposons may perturb gene integrity or alter gene expression by altering regulatory regions in the genome. The germline employs a number of mechanisms, including cytosine methylation, to repress retrotransposon transcription throughout most of life. Demethylation during germ cell and early embryo development de-represses retrotransposons. Intriguingly, de novo genetic variation appearing in sperm has been implicated in a number of disorders in offspring, including autism spectrum disorder, schizophrenia, and bipolar disorder. We hypothesize that human sperm exhibit de novo retrotransposition and employ a new sequencing method, single cell transposon insertion profiling by sequencing (scTIPseq) to map them in small amounts of human sperm.

METHODS

Cross-sectional case-control study of sperm samples (n=10 men; ages 32-55 years old) from consenting men undergoing IVF at NYU Langone Fertility Center. scTIPseq identified novel LINE-1 insertions in individual sperm and TIPseqHunter, a custom bioinformatics pipeline, compared the architecture of sperm LINE-1 to known LINE-1 insertions from the European database of Human specific LINE-1 (L1Hs) retrotransposon insertions (euL1db).

RESULTS

scTIPseq identified 17 novel insertions in sperm. New insertions were mainly intergenic or intronic. Only one sample did not exhibit new insertions. The location or number of novel insertions did not differ by paternal age.

CONCLUSION

This study for the first time reports novel LINE-1 insertions in human sperm, demonstrating the feasibility of scTIPseq, and identifies new contributors to genetic diversity in the human germ line.

The epigenetic basis of evolution

An increasing body of data are revealing key roles of epigenetics in evolutionary processes. The scope of this manuscript is to assemble in a coherent frame experimental evidence supporting a role of epigenetic factors and networks, active during embryogenesis, in orchestrating variation-inducing phenomena underlying evolution, seen as a global process. This process unfolds over two crucial levels: i) a flow of RNA-based information - predominantly small regulatory RNAs released from somatic cells exposed to environmental stimuli - taken up by spermatozoa and delivered to oocytes at fertilization and ii) the highly permissive and variation-prone environments offered by zygotes and totipotent early embryos. Totipotent embryos provide a variety of biological tools favouring the emergence of evolutionarily significant phenotypic novelties driven by RNA information. Under this light, neither random genomic mutations, nor the sieving role of natural selection are required, as the sperm-delivered RNA cargo conveys specific information and acts as "phenotypic-inducer" of defined environmentally acquired traits.

Integrated analysis of methylation profiles and transcriptome of MDV-infected chicken spleens reveal hypomethylation of CD4 and HMGB1 genes might promote MD tumorigenesis

Marek's disease (MD) is a lymphoproliferative neoplastic disease caused by Marek's disease virus (MDV). Previous studies have showed that DNA methylation was involved in MD development, but systematic studies are still lacking. Herein, we performed whole genome bisulfite sequencing (WGBS) and RNA-seq in MDV-infected tumorous spleens (IN), noninfected spleens (NoIN), and survivor (SUR) spleens of chickens to identify the genes playing important roles in MD tumor transformation. We generated the first genome-wide DNA methylation profile of MDV-infected, noninfected, and survivor chickens. Combined the WGBS and RNA-Seq, we found that the expression of 25% differential expression genes (DEGs) were significantly correlated with methylation of CpG sites in their gene bodies or promoters. Further, we focused on the DEGs with differentially methylated regions (DMRs) on genes' body and promoter, and it showed the expression of 60% DEGs were significantly correlated with methylation of CpG sites in DMRs. Finally, we identified 8 genes, including CD4, CTLA4, DTL, HMGB1, LGMN, NUP210, RAD52, and ZAP70, and their expression was negatively correlated with methylation of DMRs in their promoters in both IN vs. NoIN and IN vs. SUR. These 8 genes showed specifically high expression in IN groups and clustered in module turquoise analyzed by WGCNA. Out of 8 genes, CD4 and HMGB1 were drop in QTLs associated with MD resistance. Thus, we overexpressed the 2 genes to simulate their high expression in the IN group and found they significantly promoted MDCC-MSB-1 cell proliferation, which revealed they might play promoting roles in MD tumorigenesis in IN due to their high expression induced by hypomethylation.

An Epigenetic LINE-1-Based Mechanism in Cancer

In the last fifty years, large efforts have been deployed in basic research, clinical oncology, and clinical trials, yielding an enormous amount of information regarding the molecular mechanisms of cancer and the design of effective therapies. The knowledge that has accumulated underpins the complexity, multifactoriality, and heterogeneity of cancer, disclosing novel landscapes in cancer biology with a key role of genome plasticity. Here, we propose that cancer onset and progression are determined by a stress-responsive epigenetic mechanism, resulting from the convergence of upregulation of LINE-1 (long interspersed nuclear element 1), the largest family of human retrotransposons, genome damage, nuclear lamina fragmentation, chromatin remodeling, genome reprogramming, and autophagy activation. The upregulated expression of LINE-1 retrotransposons and their protein products plays a key role in these processes, yielding an increased plasticity of the nuclear architecture with the ensuing reprogramming of global gene expression, including the reactivation of embryonic transcription profiles. Cancer phenotypes would thus emerge as a consequence of the unscheduled reactivation of embryonic gene expression patterns in an inappropriate context, triggering de-differentiation and aberrant proliferation in differentiated cells. Depending on the intensity of the stressing stimuli and the level of LINE-1 response, diverse degrees of malignity would be generated.

Association between cancer genes and germ layer specificity

Cancer signaling pathways defining cell fates are related to differentiation. During the developmental process, three germ layers (endoderm, mesoderm, and ectoderm) are formed during embryonic development that differentiate into organs via the epigenetic regulation of specific genes. To examine the relationship, the specificities of cancer gene mutations that depend on the germ layers are studied. The major organs affected by cancer were determined based on statistics from the National Cancer Information Center of Korea, and were grouped according to their germ layer origins. Then, the gene mutation frequencies were evaluated to identify any bias based on the differentiation group using the Catalogue of Somatic Mutations in Cancer (COSMIC) database. The chi-square test showed that the p-value of 152 of 166 genes was less than 0.05, and 151 genes showed p-values of less than 0.05 even after adjusting for the false discovery rate (FDR). The germ layer-specific genes were evaluated using visualization based on basic statistics, and the results matched the top ranking genes depending on organs in the COSMIC database.The current study confirmed the germ layer specificity of major cancer genes. The germ layer specificity of mutated driver genes is possibly important in cancer treatments because each mutated gene may react differently depending on the germ layer of origin. By understanding the mechanism of gene mutation in the development and progression of cancer in the context of cell-fate pathways, a more effective therapeutic strategy for cancer can be established.

A proposed Information-Based modality for the treatment of cancer

Treatment modalities for cancer involve physical manipulations such as surgery, immunology, radiation, chemotherapy or gene editing. This is a proposal for an information-based modality. This modality does not change the internal state of the cancer cell directly - instead, the cancer cell is manipulated by giving it information to instruct the cell to perform an action. This modality is based on a theory of Structure Encoding in DNA, where information about body part structure controls the epigenetic state of cells in the process of development from pluripotent cells to fully differentiated cells. It has been noted that cancer is often due to errors in morphogenetic differentiation accompanied by associated epigenetic processes. This implies a model of cancer called the Epigenetic Differentiation Model. A major feature of the Structure Encoding Theory is that the characteristics of the differentiated cell are affected by inter-cellular information passed in the tissue microenvironment, which specifies the exact location of a cell in a body part structure. This is done by exosomes that carry fragments of long non-coding RNA and transposons, which convey structure information. In the normal process of epigenetic differentiation, the information passed may lead to apoptosis due to the constraints of a particular body part structure. The proposed treatment involves determining what structure information is being passed in a particular tumor, then adding artificial exosomes that overwhelm the current information with commands for the cells to go into apoptosis.

Zidovudine inhibits telomere elongation, increases the transposable element LINE-1 copy number and compromises mouse embryo development

PURPOSE

Millions of pregnant, HIV-infected women take reverse transcriptase inhibitors, such as zidovudine (azidothymidine or AZT), during pregnancy. Reverse transcription plays important roles in early development, including regulation of telomere length (TL) and activity of transposable elements (TE). So we evaluated the effects of AZT on embryo development, TL, and copy number of an active TE, Long Interspersed Nuclear Element 1 (LINE-1), during early development in a murine model.

DESIGN

Experimental study.

METHODS

In vivo fertilized mouse zygotes from B6C3F1/B6D2F1 mice were cultured for 48 h in KSOM with no AZT (n = 45), AZT 1 μM (n = 46) or AZT 10 μM (n = 48). TL was measured by single-cell quantitative PCR (SC-pqPCR) and LINE-1 copy number by qPCR. The percentage of morulas at 48 h, TL and LINE-1 copy number were compared among groups.

RESULTS

Exposure to AZT 1 μM or 10 μM significantly impairs early embryo development. TL elongates from oocyte to control embryos. TL in AZT 1 μM embryos is shorter than in control embryos. LINE-1 copy number is significantly lower in oocytes than control embryos. AZT 1 μM increases LINE-1 copy number compared to oocytes controls, and AZT 10 μM embryos.

CONCLUSION

AZT at concentrations approaching those used to prevent perinatal HIV transmission compromises mouse embryo development, prevents telomere elongation and increases LINE-1 copy number after 48 h treatment. The impact of these effects on the trajectory of aging of children exposed to AZT early during development deserves further investigation.

Eukaryotic Genomes Show Strong Evolutionary Conservation of k-mer Composition and Correlation Contributions between Introns and Intergenic Regions

Several strongly conserved DNA sequence patterns in and between introns and intergenic regions (IIRs) consisting of short tandem repeats (STRs) with repeat lengths <3 bp have already been described in the kingdom of Animalia. In this work, we expanded the search and analysis of conserved DNA sequence patterns to a wider range of eukaryotic genomes. Our aims were to confirm the conservation of these patterns, to support the hypothesis on their functional constraints and/or the identification of unknown patterns. We pairwise compared genomic DNA sequences of genes, exons, CDS, introns and intergenic regions of 34 Embryophyta (land plants), 30 Protista and 29 Fungi using established k-mer-based (alignment-free) comparison methods. Additionally, the results were compared with values derived for Animalia in former studies. We confirmed strong correlations between the sequence structures of IIRs spanning over the entire domain of Eukaryotes. We found that the high correlations within introns, intergenic regions and between the two are a result of conserved abundancies of STRs with repeat units ≤2 bp (e.g., (AT)n). For some sequence patterns and their inverse complementary sequences, we found a violation of equal distribution on complementary DNA strands in a subset of genomes. Looking at mismatches within the identified STR patterns, we found specific preferences for certain nucleotides stable over all four phylogenetic kingdoms. We conclude that all of these conserved patterns between IIRs indicate a shared function of these sequence structures related to STRs.

MicroRNAs' role in the environment-related non-communicable diseases and link to multidrug resistance, regulation, or alteration

The discovery of microRNAs (miRNAs) 20 years ago has advocated a new era of "small molecular genetics." About 2000 miRNAs are present that regulate one third of the genome. MiRNA dysregulated expression arising as a response to our environment insult or stress or changes may contribute to several diseases, namely non-communicable diseases, including tumor growth. Their presence in body fluids, reflecting level alteration in various cancers, merit circulating miRNAs as the "next-generation biomarkers" for early-stage tumor diagnosis and/or prognosis. Herein, we performed a comprehensive literature search focusing on the origin, biosynthesis, and role of miRNAs and summarized the foremost studies centering on miR value as non-invasive biomarkers in different environment-related non-communicable diseases, including various cancer types. Moreover, during chemotherapy, many miRNAs were linked to multidrug resistance, via modulating numerous, environment triggered or not, biological processes and/or pathways that will be highlighted as well.

Control of LINE-1 Expression Maintains Genome Integrity in Germline and Early Embryo Development

Maintenance of genome integrity in the germline and in preimplantation embryos is crucial for mammalian development. Epigenetic remodeling during primordial germ cell (PGC) and preimplantation embryo development may contribute to genomic instability in these cells, since DNA methylation is an important mechanism to silence retrotransposons. Long interspersed elements 1 (LINE-1 or L1) are the most common autonomous retrotransposons in mammals, corresponding to approximately 17% of the human genome. Retrotransposition events are more frequent in germ cells and in early stages of embryo development compared with somatic cells. It has been shown that L1 activation and expression occurs in germline and is essential for preimplantation development. In this review, we focus on the role of L1 retrotransposon in mouse and human germline and early embryo development and discuss the possible relationship between L1 expression and genomic instability during these stages. Although several studies have addressed L1 expression at different stages of development, the developmental consequences of this expression remain poorly understood. Future research is still needed to highlight the relationship between L1 retrotransposition events and genomic instability during germline and early embryo development.

Differential genetic mutations of ectoderm, mesoderm, and endoderm-derived tumors in TCGA database

Background

In terms of biological behavior, gene regulation, or signaling pathways, there is a certain similarity between tumorigenesis and embryonic development of humans. Three germ layer structure exhibits the distinct ability to form specific tissues and organs.

Methods

The present study set out to investigate the genetic mutation characteristics of germ layer differentiation-related genes using the tumor cases of the cancer genome atlas (TCGA) database.

Results

These tumor samples were divided into three groups, including the ectoderm, mesoderm, and endoderm. Children cases less than 9 years old accounted for a larger proportion for the cases in the ectoderm and mesoderm groups; whereas the middle-aged and elderly individuals (from 50 to 89 years old) were more susceptible to tumors of endoderm. There was a better prognosis for the cases of mesoderm, especially the male with the race of White, compared with the other groups. A missense mutation was frequently detected for the cases of ectoderm and endoderm, while deletion mutation was common for that of mesoderm. We could not identify the ectoderm, mesoderm, or endoderm-specific mutated genes or variants with high mutation frequency. However, there was a relatively higher mutation incidence of endoderm markers (GATA6, FOXA2, GATA4, AFP) in the endoderm group, compared with the groups of ectoderm and mesoderm. Additionally, four members (SMO, GLI1, GLI2, GLI3) within the Hedgehog signaling pathway genes showed a relatively higher mutation rate in the endoderm group than the other two groups.

Conclusions

TCGA tumors of ectoderm, mesoderm, and endoderm groups exhibit the distinct subject distribution, survival status, and genomic alteration characteristics. The synergistic mutation effect of specific genes closely related to embryonic development may contribute to the tumorigenesis of tissues or organs derived from the specific germ layers. This study provides a novel reference for exploring the functional connection between embryogenesis and tumorigenesis.

RNA element discovery from germ cell to blastocyst

Recent studies have shown that tissue-specific transcriptomes contain multiple types of RNAs that are transcribed from intronic and intergenic sequences. The current study presents a tool for the discovery of transcribed, unannotated sequence elements from RNA-seq libraries. This RNA Element (RE) discovery algorithm (REDa) was applied to a spectrum of tissues and cells representing germline, embryonic, and somatic tissues and examined as a function of differentiation through the first set of cell divisions of human development. This highlighted extensive transcription throughout the genome, yielding previously unidentified human spermatogenic RNAs. Both exonic and novel X-chromosome REs were subject to robust meiotic sex chromosome inactivation, although an extensive de-repression occurred in the post-meiotic stages of spermatogenesis. Surprisingly, 2.4% of the 10,395 X chromosome exonic REs were present in mature sperm. Transcribed genomic repetitive sequences, including simple centromeric repeats, HERVE and HSAT1, were also shown to be associated with RE expression during spermatogenesis. These results suggest that pervasive intergenic repetitive sequence expression during human spermatogenesis may play a role in regulating chromatin dynamics. Repetitive REs switching repeat classes during differentiation upon fertilization and embryonic genome activation was evident.

Conservation of k-mer Composition and Correlation Contribution between Introns and Intergenic Regions of Animalia Genomes

In this study, we pairwise-compared multiple genome regions, including genes, exons, coding DNA sequences (CDS), introns, and intergenic regions of 39 Animalia genomes, including Deuterostomia (27 species) and Protostomia (12 species), by applying established k-mer-based (alignment-free) comparison methods. We found strong correlations between the sequence structure of introns and intergenic regions, individual organisms, and within wider phylogenetical ranges, indicating the conservation of certain structures over the full range of analyzed organisms. We analyzed these sequence structures by quantifying the contribution of different sets of DNA words to the average correlation value by decomposing the correlation coefficients with respect to these word sets. We found that the conserved structures within introns, intergenic regions, and between the two were mainly a result of conserved tandem repeats with repeat units ≤ 2 bp (e.g., (AT)n), while other conserved sequence structures, such as those found between exons and CDS, were dominated by tandem repeats with repeat unit sizes of 3 bp in length and more complex DNA word patterns. We conclude that the conservation between intron and intergenic regions indicates a shared function of these sequence structures. Also, the similar differences in conserved structures with known origin, especially to the conservation between exons and CDS resulting from DNA codons, indicate that k-mer composition-based functional properties of introns and intergenic regions may differ from those of exons and CDS.

The "evolutionary field" hypothesis. Non-Mendelian transgenerational inheritance mediates diversification and evolution

Epigenetics is increasingly regarded as a potential contributing factor to evolution. Building on apparently unrelated results, here I propose that RNA-containing nanovesicles, predominantly small regulatory RNAs, are released from somatic tissues in the bloodstream, cross the Weismann barrier, reach the epididymis, and are eventually taken up by spermatozoa; henceforth the information is delivered to oocytes at fertilization. In the model, a LINE-1-encoded reverse transcriptase activity, present in spermatozoa and early embryos, plays a key role in amplifying and propagating these RNAs as extrachromosomal structures. It may be conceived that, over generations, the cumulative effects of sperm-delivered RNAs would cross a critical threshold and overcome the buffering capacity of embryos. As a whole, the process can promote the generation of an information-containing platform that drives the reshaping of the embryonic epigenetic landscape with the potential to generate ontogenic changes and redirect the evolutionary trajectory. Over time, evolutionary significant, stably acquired variations could be generated through the process. The interplay between these elements defines the concept of "evolutionary field", a self-consistent, comprehensive information-containing platform and a source of discontinuous evolutionary novelty.

Two genetic codes: Repetitive syntax for active non-coding RNAs; non-repetitive syntax for the DNA archives

Current knowledge of the RNA world indicates 2 different genetic codes being present throughout the living world. In contrast to non-coding RNAs that are built of repetitive nucleotide syntax, the sequences that serve as templates for proteins share-as main characteristics-a non-repetitive syntax. Whereas non-coding RNAs build groups that serve as regulatory tools in nearly all genetic processes, the coding sections represent the evolutionarily successful function of the genetic information storage medium. This indicates that the differences in their syntax structure are coherent with the differences of the functions they represent. Interestingly, these 2 genetic codes resemble the function of all natural languages, i.e., the repetitive non-coding sequences serve as appropriate tool for organization, coordination and regulation of group behavior, and the non-repetitive coding sequences are for conservation of instrumental constructions, plans, blueprints for complex protein-body architecture. This differentiation may help to better understand RNA group behavioral motifs.

Role of microRNAs in translation regulation and cancer

MicroRNAs (miRNAs) are pervasively expressed and regulate most biological functions. They function by modulating transcriptional and translational programs and therefore they orchestrate both physiological and pathological processes, such as development, cell differentiation, proliferation, apoptosis and tumor growth. miRNAs work as small guide molecules in RNA silencing, by negatively regulating the expression of several genes both at mRNA and protein level, by degrading their mRNA target and/or by silencing translation. One of the most recent advances in the field is the comprehension of their role in oncogenesis. The number of miRNA genes is increasing and an alteration in the level of miRNAs is involved in the initiation, progression and metastases formation of several tumors. Some tumor types show a distinct miRNA signature that distinguishes them from normal tissues and from other cancer types. Genetic and biochemical evidence supports the essential role of miRNAs in tumor development. Although the abnormal expression of miRNAs in cancer cells is a widely accepted phenomenon, the cause of this dysregulation is still unknown. Here, we discuss the biogenesis of miRNAs, focusing on the mechanisms by which they regulate protein synthesis. In addition we debate on their role in cancer, highlighting their potential to become therapeutic targets.

The Reverse Transcriptase Encoded by LINE-1 Retrotransposons in the Genesis, Progression, and Therapy of Cancer

In higher eukaryotic genomes, Long Interspersed Nuclear Element 1 (LINE-1) retrotransposons represent a large family of repeated genomic elements. They transpose using a reverse transcriptase (RT), which they encode as part of the ORF2p product. RT inhibition in cancer cells, either via RNA interference-dependent silencing of active LINE-1 elements, or using RT inhibitory drugs, reduces cancer cell proliferation, promotes their differentiation and antagonizes tumor progression in animal models. Indeed, the non-nucleoside RT inhibitor efavirenz has recently been tested in a phase II clinical trial with metastatic prostate cancer patients. An in-depth analysis of ORF2p in a mouse model of breast cancer showed ORF2p to be precociously expressed in precancerous lesions and highly abundant in advanced cancer stages, while being barely detectable in normal breast tissue, providing a rationale for the finding that RT-expressing tumors are therapeutically sensitive to RT inhibitors. We summarize mechanistic and gene profiling studies indicating that abundant LINE-1-derived RT can “sequester” RNA substrates for reverse transcription in tumor cells, entailing the formation of RNA:DNA hybrid molecules and impairing the overall production of regulatory miRNAs, with a global impact on the cell transcriptome. Based on these data, LINE-1-ORF2 encoded RT has a tumor-promoting potential that is exerted at an epigenetic level. We propose a model whereby LINE1-RT drives a previously unrecognized global regulatory process, the deregulation of which drives cell transformation and tumorigenesis with possible implications for cancer cell heterogeneity.

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.