Retroelement distributions in the human genome: variations associated with age and proximity to genes

Different Founding Effects Underlie Dominant Blue Eyes (DBE) in the Domestic Cat

During the last twenty years, minimal white spotting associated with blue eyes was selected by feline breeders to create the Altai, Topaz, and Celestial breeds. Additionally, certain breeders introduced this trait in their lineages of purebred cats. The trait has been called "dominant blue eyes (DBE)" and was confirmed to be autosomal dominant in all lineages. DBE was initially described in outbred cats from Kazakhstan and Russia and in two purebred lineages of British cats from Russia, as well as in Dutch Maine Coon cats, suggesting different founding effects. We have previously identified two variants in the Paired Box 3 (PAX3) gene associated with DBE in Maine Coon and Celestial cats; however, the presence of an underlying variant remains undetermined in other DBE breeding lines. Using a genome-wide association study, we identified a single region on chromosome C1 that was associated with DBE in British cats. Within that region, we identified PAX3 as the strongest candidate gene. Whole-genome sequencing of a DBE cat revealed an RD-114 retrovirus LTR (long terminal repeat) insertion within PAX3 intron 4 (namely NC_018730.3:g.206975776_206975777insN[433]) known to contain regulatory sequences. Using a panel of 117 DBE cats, we showed that this variant was fully associated with DBE in two British lineages, in Altai cats, and in some other DBE lineages. We propose that this NC_018730.3:g.206975776_206975777insN[433] variant represents the DBEALT (Altai Dominant Blue Eye) allele in the domestic cat. Finally, we genotyped DBE cats from 14 lineages for the three PAX3 variants and showed that they were not present in four lineages, confirming genetic heterogeneity of the DBE trait in the domestic cat.

The impact of differential transposition activities of autonomous and nonautonomous hAT transposable elements on genome architecture and gene expression in Caenorhabditis inopinata

Transposable elements are DNA sequences capable of moving within genomes and significantly influence genomic evolution. The nematode Caenorhabditis inopinata exhibits a much higher transposable element copy number than its sister species, Caenorhabditis elegans. In this study, we identified a novel autonomous transposable element belonging to the hAT superfamily from a spontaneous transposable element-insertion mutant in C. inopinata and named this transposon Ci-hAT1. Further bioinformatic analyses uncovered 3 additional autonomous hAT elements-Ci-hAT2, Ci-hAT3, and Ci-hAT4-along with over 1,000 copies of 2 nonautonomous miniature inverted-repeat transposable elements, mCi-hAT1 and mCi-hAT4, likely derived from Ci-hAT1 and Ci-hAT4 through internal deletion. We tracked at least 3 sequential transpositions of Ci-hAT1 over several years. However, the transposition rates of the other 3 autonomous hAT elements were lower, suggesting varying activity levels. Notably, the distribution patterns of the 2 miniature inverted-repeat transposable element families differed significantly: mCi-hAT1 was primarily located in the chromosome arms, a pattern observed in the transposable elements of other Caenorhabditis species, whereas mCi-hAT4 was more evenly distributed across chromosomes. Additionally, interspecific transcriptome analysis indicated that C. inopinata genes with upstream or intronic these miniature inverted-repeat transposable element insertions tend to be more highly expressed than their orthologous genes in C. elegans. These findings highlight the significant role of de-silenced transposable elements in driving the evolution of genomes and transcriptomes, leading to species-specific genetic diversity.

Endogenous feline leukemia virus long terminal repeat integration site diversity is highly variable in related and unrelated domestic cats

Endogenous retroviruses (ERV) are indicators of vertebrate evolutionary history and play important roles as homeostatic regulators. ERV long terminal repeat (LTR) elements may act as cis-activating promoters or trans-activating enhancer elements modifying gene transcription distant from LTR insertion sites. We previously documented that endogenous feline leukemia virus (FeLV)-LTR copy number variation in individual cats tracks inversely with susceptibility to virulent FeLV disease. To evaluate FeLV-LTR insertion characteristics, we assessed enFeLV-LTR integration site diversity in 20 cats from three genetically distinct populations using a baited linker-mediated PCR approach. We documented 765 individual integration sites unequally represented among individuals. Only three LTR integration sites were shared among all individuals, while 412 sites were unique to a single individual. When primary fibroblast cultures were challenged with exogenous FeLV, we found significantly increased expression of both exogenous and endogenous FeLV orthologs, supporting previous findings of potential exFeLV-enFeLV interactions; however, viral challenge did not elicit transcriptional changes in genes associated with the vast majority of integration sites. This study assesses FeLV-LTR integration sites in individual animals, providing unique transposome genotypes. Further, we document substantial individual variation in LTR integration site locations, even in a highly inbred population, and provide a framework for understanding potential endogenous retroviral element position influence on host gene transcription.

Recent insights into crosstalk between genetic parasites and their host genome

The bulk of higher order organismal genomes is comprised of transposable element (TE) copies, i.e. genetic parasites. The host-parasite relation is multi-faceted, varying across genomic region (genic versus intergenic), life-cycle stages, tissue-type and of course in health versus pathological state. The reach of functional genomics though, in investigating genotype-to-phenotype relations, has been limited when TEs are involved. The aim of this review is to highlight recent progress made in understanding how TE origin biochemical activity interacts with the central dogma stages of the host genome. Such interaction can also bring about modulation of the immune context and this could have important repercussions in disease state where immunity has a role to play. Thus, the review is to instigate ideas and action points around identifying evolutionary adaptations that the host genome and the genetic parasite have evolved and why they could be relevant.

The pig pangenome provides insights into the roles of coding structural variations in genetic diversity and adaptation

Structural variations have emerged as an important driving force for genome evolution and phenotypic variation in various organisms, yet their contributions to genetic diversity and adaptation in domesticated animals remain largely unknown. Here we constructed a pangenome based on 250 sequenced individuals from 32 pig breeds in Eurasia and systematically characterized coding sequence presence/absence variations (PAVs) within pigs. We identified 308.3-Mb nonreference sequences and 3438 novel genes absent from the current reference genome. Gene PAV analysis showed that 16.8% of the genes in the pangene catalog undergo PAV. A number of newly identified dispensable genes showed close associations with adaptation. For instance, several novel swine leukocyte antigen (SLA) genes discovered in nonreference sequences potentially participate in immune responses to productive and respiratory syndrome virus (PRRSV) infection. We delineated previously unidentified features of the pig mobilome that contained 490,480 transposable element insertion polymorphisms (TIPs) resulting from recent mobilization of 970 TE families, and investigated their population dynamics along with influences on population differentiation and gene expression. In addition, several candidate adaptive TE insertions were detected to be co-opted into genes responsible for responses to hypoxia, skeletal development, regulation of heart contraction, and neuronal cell development, likely contributing to local adaptation of Tibetan wild boars. These findings enhance our understanding on hidden layers of the genetic diversity in pigs and provide novel insights into the role of SVs in the evolutionary adaptation of mammals.

The Regulation and Immune Signature of Retrotransposons in Cancer

Advances in sequencing technologies and the bioinformatic analysis of big data facilitate the study of jumping genes' activity in the human genome in cancer from a broad perspective. Retrotransposons, which move from one genomic site to another by a copy-and-paste mechanism, are regulated by various molecular pathways that may be disrupted during tumorigenesis. Active retrotransposons can stimulate type I IFN responses. Although accumulated evidence suggests that retrotransposons can induce inflammation, the research investigating the exact mechanism of triggering these responses is ongoing. Understanding these mechanisms could improve the therapeutic management of cancer through the use of retrotransposon-induced inflammation as a tool to instigate immune responses to tumors.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Widespread somatic L1 retrotransposition in normal colorectal epithelium

Throughout an individual's lifetime, genomic alterations accumulate in somatic cells^1-11. However, the mutational landscape induced by retrotransposition of long interspersed nuclear element-1 (L1), a widespread mobile element in the human genome^12-14, is poorly understood in normal cells. Here we explored the whole-genome sequences of 899 single-cell clones established from three different cell types collected from 28 individuals. We identified 1,708 somatic L1 retrotransposition events that were enriched in colorectal epithelium and showed a positive relationship with age. Fingerprinting of source elements showed 34 retrotransposition-competent L1s. Multidimensional analysis demonstrated that (1) somatic L1 retrotranspositions occur from early embryogenesis at a substantial rate, (2) epigenetic on/off of a source element is preferentially determined in the early organogenesis stage, (3) retrotransposition-competent L1s with a lower population allele frequency have higher retrotransposition activity and (4) only a small fraction of L1 transcripts in the cytoplasm are finally retrotransposed in somatic cells. Analysis of matched cancers further suggested that somatic L1 retrotransposition rate is substantially increased during colorectal tumourigenesis. In summary, this study illustrates L1 retrotransposition-induced somatic mosaicism in normal cells and provides insights into the genomic and epigenomic regulation of transposable elements over the human lifetime.

Identification and Characterization of the HERV-K (HML-8) Group of Human Endogenous Retroviruses in the Genome

Human endogenous retroviruses (HERVs) can be vertically transmitted in a Mendelian fashion, are stably maintained in the human genome, and are estimated to constitute ∼8% of the genome. HERVs affect human physiology and pathology through their provirus-encoded protein or long terminal repeat (LTR) element effect. Characterization of the genomic distribution is an essential step to understanding the relationships between endogenous retrovirus expression and diseases. However, the poor characterization of human MMTV-like (HML)-8 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, the definition of an accurate and updated HERV-K HML-8 genomic map is urgently needed. In this study, we report the results of a comprehensive analysis of HERV-K HML-8 sequence presence and distribution within the human genome and hominoids, with a detailed description of the different structural and phylogenetic aspects characterizing the group. A total of 40 proviruses and 5 solo LTR elements for human were characterized, which included a detailed description of provirus structure, integration time, potentially regulated genes, transcription factor-binding sites, and primer-binding site features. Besides, 9 chimpanzee sequences, 8 gorilla sequences, and 10 orangutan sequences belonging to the HML-8 subgroup were identified. The integration time results showed that the HML-8 elements were integrated into the primate lineage around 35 and 42 million years ago (mya), during primates evolutionary speciation. Overall, the results clarified the composition of the HML-8 groups, providing an exhaustive background for subsequent functional studies.

Transposable element and host silencing activity in gigantic genomes

Transposable elements (TEs) and the silencing machinery of their hosts are engaged in a germline arms-race dynamic that shapes TE accumulation and, therefore, genome size. In animal species with extremely large genomes (>10 Gb), TE accumulation has been pushed to the extreme, prompting the question of whether TE silencing also deviates from typical conditions. To address this question, we characterize TE silencing via two pathways-the piRNA pathway and KRAB-ZFP transcriptional repression-in the male and female gonads of Ranodon sibiricus, a salamander species with a ∼21 Gb genome. We quantify 1) genomic TE diversity, 2) TE expression, and 3) small RNA expression and find a significant relationship between the expression of piRNAs and TEs they target for silencing in both ovaries and testes. We also quantified TE silencing pathway gene expression in R. sibiricus and 14 other vertebrates with genome sizes ranging from 1 to 130 Gb and find no association between pathway expression and genome size. Taken together, our results reveal that the gigantic R. sibiricus genome includes at least 19 putatively active TE superfamilies, all of which are targeted by the piRNA pathway in proportion to their expression levels, suggesting comprehensive piRNA-mediated silencing. Testes have higher TE expression than ovaries, suggesting that they may contribute more to the species' high genomic TE load. We posit that apparently conflicting interpretations of TE silencing and genomic gigantism in the literature, as well as the absence of a correlation between TE silencing pathway gene expression and genome size, can be reconciled by considering whether the TE community or the host is currently "on the attack" in the arms race dynamic.

Transposable Elements in Bats Show Differential Accumulation Patterns Determined by Class and Functionality

Bat genomes are characterized by a diverse transposable element (TE) repertoire. In particular, the genomes of members of the family Vespertilionidae contain both active retrotransposons and active DNA transposons. Each TE type is characterized by a distinct pattern of accumulation over the past ~40 million years. Each also exhibits its own target site preferences (sometimes shared with other TEs) that impact where they are likely to insert when mobilizing. Therefore, bats provide a great resource for understanding the diversity of TE insertion patterns. To gain insight into how these diverse TEs impact genome structure, we performed comparative spatial analyses between different TE classes and genomic features, including genic regions and CpG islands. Our results showed a depletion of all TEs in the coding sequence and revealed patterns of species- and element-specific attraction in the transcript. Trends of attraction in the distance tests also suggested significant TE activity in regions adjacent to genes. In particular, the enrichment of small, non-autonomous TE insertions in introns and near coding regions supports the hypothesis that the genomic distribution of TEs is the product of a balance of the TE insertion preference in open chromatin regions and the purifying selection against TEs within genes.

TEspeX: consensus-specific quantification of transposable element expression preventing biases from exonized fragments

SUMMARY

Transposable elements (TEs) play key roles in crucial biological pathways. Therefore, several tools enabling the quantification of their expression were recently developed. However, many of the existing tools lack the capability to distinguish between the transcription of autonomously expressed TEs and TE fragments embedded in canonical coding/non-coding non-TE transcripts. Consequently, an apparent change in the expression of a given TE may simply reflect the variation in the expression of the transcripts containing TE-derived sequences. To overcome this issue, we have developed TEspeX, a pipeline for the quantification of TE expression at the consensus level. TEspeX uses Illumina RNA-seq short reads to quantify TE expression avoiding counting reads deriving from inactive TE fragments embedded in canonical transcripts.

AVAILABILITY AND IMPLEMENTATION

The tool is implemented in python3, distributed under the GNU General Public License (GPL) and available on Github at https://github.com/fansalon/TEspeX (Zenodo URL: https://doi.org/10.5281/zenodo.6800331).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

p53 Binding Sites in Long Terminal Repeat 5Hs (LTR5Hs) of Human Endogenous Retrovirus K Family (HML-2 Subgroup) Play Important Roles in the Regulation of LTR5Hs Transcriptional Activity

The long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) are distributed throughout the human genome and provide favorable conditions to regulate the expression of their adjacent genes. HML-2 is the most biologically active subgroup of the HERV-K family, and expression of its members has been associated with many cancer types. The LTRs of HML-2 have been classified into three subgroups (LTR5A, LTR5B, and LTR5Hs) based on phylogenetic analyses. The current study aimed to explore the LTR transcriptional activity differences among the three subtypes and further explore the underlying factors. A total of 43 LTR5A elements, 62 LTR5B elements, and 194 LTR5Hs elements were selected. A phylogenetic tree showed that the LTR5Hs group was clearly separated from the LTR5A and LTR5B groups. A luciferase reporter assay indicated that LTR5Hs had the strongest promoter activity, followed by LTR5A and LTR5B. To investigate the underlying factors, LTR5Hs was divided into 4 sections, and the homologous fragments in LTR5B were replaced successively. Replacement of the third section (−263 to 0) significantly increased LTR5B activity. Subsequent mutation experiments revealed that the increased transcriptional activity was induced by the TATA box and the two p53 binding sites within the section. Further interference with TP53 significantly decreased LTR5Hs transcriptional activity. Chromatin immunoprecipitation (ChIP) and CUT&Tag experiments finally confirmed the direct binding of the p53 protein with the two LTR5Hs p53 binding sites. Overall, the two p53 binding sites in the third section (−263 to 0) of LTR5Hs were revealed to play critical roles in the difference in transcriptional activity among the three subtypes.

IMPORTANCE Human endogenous retroviruses (HERVs) were integrated into the human genome in ancient times and have been coevolving with the host. Since the Human Genome Project, HERVs have attracted increasing attention. Many studies have focused on their characterization, evolution, and biological function. In particular, the expression of HERV-K has been associated with many diseases, such as germ cell tumors, neurotoxicity, ovarian cancer, prostate cancer, and melanoma. Indeed, two HML-2-produced proteins, Np9 and Rec, are associated with certain cancers. However, their roles in these disease associations remain unclear. The current work focused on subgroup HML-2 of HERV-K, which is recognized as the most biologically active subgroup, and aimed to explore the mechanistic basis of transcriptional activity. The results revealed that p53 deeply determined the activity of HML-2 LTR5Hs. p53 is a rather important tumor suppressor protein. It can regulate the expression of genes related to cell cycle arrest, organic processes, and apoptosis in response to cellular stress and is critical for the control of homeostasis. Previous ChIP and expression studies of individual genes suggested that p53 sites in HERV LTRs may be part of the p53 transcription program and directly regulate p53 target genes in a species-specific manner. However, the exact function of p53 in the regulation of HERV LTR expression is largely elusive. Our results clearly demonstrated the interaction between LTR5Hs of HML-2 and p53. They are of great significance for the future comprehensive study of the physiological and pathological functions of LTRs of HERVs.

ExoPLOT: Representation of alternative splicing in human tissues and developmental stages with transposed element (TE) involvement

Advances in RNA high-throughput sequencing and large-scale functional assays yield new insights into the multifaceted activities of transposed elements (TE) and many other previously undiscovered sequence elements. Currently, no tool for easy access, analysis, quantification, and visualization of alternatively spliced exons across multiple tissues or developmental stages is available. Also, analysis pipelines demand computational skills or hardware requirements, which often are hard to meet by wet-lab scientists. We developed ExoPLOT to enable simplified access to massive RNA high throughput sequencing datasets to facilitate the analysis of alternative splicing across many biological samples. To demonstrate the functonality of ExoPLOT, we analyzed the contributon of exonized TEs to human coding sequences (CDS). mRNA splice variants containing the TE-derived exon were quantified and compared to expression levels of TE-free splice variants. For analysis, we utilized 313 human cerebrum, cerebellum, heart, kidney, liver, ovary, and testis transcriptomes, representing various pre- and postnatal developmental stages. ExoPLOT visualizes the relative expression levels of alternative transcripts, e.g., caused by the insertion of new TE-derived exons, across different developmental stages of and among multiple tissues. This tool also provides a unique link between evolution and function during exonization (gain of a new exon) and exaptation (recruitment/co-optation) of a new exon. As input for analysis, we derived a database of 1151 repeat-masked, exonized TEs, representing all prominent families of transposons in the human genome and the collection of human consensus coding sequences (CCDS). ExoPLOT screened preprocessed RNA high-throughput sequencing datasets from seven human tissues to quantify and visualize the dynamics in RNA splicing for these 1151 TE-derived exons during the entire human organ development. In addition, we successfully mapped and analyzed 993 recently described exonized sequences from the human frontal cortex onto these 313 transcriptome libraries. ExoPLOT's approach to preprocessing RNA deep sequencing datasets facilitates alternative splicing analysis and significantly reduces processing times. In addition, ExoPLOT's design allows studying alternative RNA isoforms other than TE-derived in a customized - coordinate-based manner and is available at http://retrogenomics3.uni-muenster.de:3838/exz-plot-d/.

Zebrafish transposable elements show extensive diversification in age, genomic distribution, and developmental expression

There is considerable interest in understanding the effect of transposable elements (TEs) on embryonic development. Studies in humans and mice are limited by the difficulty of working with mammalian embryos and by the relative scarcity of active TEs in these organisms. The zebrafish is an outstanding model for the study of vertebrate development, and over half of its genome consists of diverse TEs. However, zebrafish TEs remain poorly characterized. Here we describe the demography and genomic distribution of zebrafish TEs and their expression throughout embryogenesis using bulk and single-cell RNA sequencing data. These results reveal a highly dynamic genomic ecosystem comprising nearly 2000 distinct TE families, which vary in copy number by four orders of magnitude and span a wide range of ages. Longer retroelements tend to be retained in intergenic regions, whereas short interspersed nuclear elements (SINEs) and DNA transposons are more frequently found nearby or within genes. Locus-specific mapping of TE expression reveals extensive TE transcription during development. Although two-thirds of TE transcripts are likely driven by nearby gene promoters, we still observe stage- and tissue-specific expression patterns in self-regulated TEs. Long terminal repeat (LTR) retroelements are most transcriptionally active immediately following zygotic genome activation, whereas DNA transposons are enriched among transcripts expressed in later stages of development. Single-cell analysis reveals several endogenous retroviruses expressed in specific somatic cell lineages. Overall, our study provides a valuable resource for using zebrafish as a model to study the impact of TEs on vertebrate development.

Comprehensive identification and characterization of the HERV-K (HML-9) group in the human genome

Background

Human endogenous retroviruses (HERVs) result from ancestral infections caused by exogenous retroviruses that became incorporated into the germline DNA and evolutionarily fixed in the human genome. HERVs can be transmitted vertically in a Mendelian fashion and be stably maintained in the human genome, of which they are estimated to comprise approximately 8%. HERV-K (HML1-10) transcription has been confirmed to be associated with a variety of diseases, such as breast cancer, lung cancer, prostate cancer, melanoma, rheumatoid arthritis, and amyotrophic lateral sclerosis. However, the poor characterization of HML-9 prevents a detailed understanding of the regulation of the expression of this family in humans and its impact on the host genome. In light of this, a precise and updated HERV-K HML-9 genomic map is urgently needed to better evaluate the role of these elements in human health.

Results

We report a comprehensive analysis of the presence and distribution of HERV-K HML-9 elements within the human genome, with a detailed characterization of the structural and phylogenetic properties of the group. A total of 23 proviruses and 47 solo LTR elements were characterized, with a detailed description of the provirus structure, integration time, potential regulated genes, transcription factor binding sites (TFBS), and primer binding site (PBS) features. The integration time results showed that the HML-9 elements found in the human genome integrated into the primate lineage between 17.5 and 48.5 million years ago (mya).

Conclusion

The results provide a clear characterization of HML-9 and a comprehensive background for subsequent functional studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12977-022-00596-2.

Comprehensive In Silico Analysis of Retrotransposon Insertions within the Survival Motor Neuron Genes Involved in Spinal Muscular Atrophy

Transposable elements (TEs) are interspersed repetitive and mobile DNA sequences within the genome. Better tools for evaluating TE-derived sequences have provided insights into the contribution of TEs to human development and disease. Spinal muscular atrophy (SMA) is an autosomal recessive motor neuron disease that is caused by deletions or mutations in the Survival Motor Neuron 1 (SMN1) gene but retention of its nearly perfect orthologue SMN2. Both genes are highly enriched in TEs. To establish a link between TEs and SMA, we conducted a comprehensive, in silico analysis of TE insertions within the SMN1/2 loci of SMA, carrier and healthy genomes. We found an Alu insertion in the promoter region and one L1 element in the 3'UTR that may play an important role in alternative promoter as well as in alternative transcriptional termination. Additionally, several intronic Alu repeats may influence alternative splicing via RNA circularization and causes the presence of new alternative exons. These Alu repeats present throughout the genes are also prone to recombination events that could lead to SMN1 exons deletions and, ultimately, SMA. TE characterization of the SMA genomic region could provide for a better understanding of the implications of TEs on human disease and genomic evolution.

On the Base Composition of Transposable Elements

Transposable elements exhibit a base composition that is often different from the genomic average and from hosts’ genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.

Characterizing mobile element insertions in 5675 genomes

Mobile element insertions (MEIs) are a major class of structural variants (SVs) and have been linked to many human genetic disorders, including hemophilia, neurofibromatosis, and various cancers. However, human MEI resources from large-scale genome sequencing are still lacking compared to those for SNPs and SVs. Here, we report a comprehensive map of 36 699 non-reference MEIs constructed from 5675 genomes, comprising 2998 Chinese samples (∼26.2×, NyuWa) and 2677 samples from the 1000 Genomes Project (∼7.4×, 1KGP). We discovered that LINE-1 insertions were highly enriched in centromere regions, implying the role of chromosome context in retroelement insertion. After functional annotation, we estimated that MEIs are responsible for about 9.3% of all protein-truncating events per genome. Finally, we built a companion database named HMEID for public use. This resource represents the latest and largest genomewide study on MEIs and will have broad utility for exploration of human MEI findings.

Clustering of Sex-Biased Genes and Transposable Elements in the Genome of the Medaka Fish Oryzias latipes

Although genes with similar expression patterns are sometimes found in the same genomic regions, almost nothing is known about the relative organization in genomes of genes and transposable elements (TEs), which might influence each other at the regulatory level. In this study, we used transcriptomic data from male and female gonads of the Japanese medaka Oryzias latipes to define sexually biased genes and TEs and analyze their relative genomic localization. We identified 20,588 genes expressed in the adult gonads of O. latipes. Around 39% of these genes are differentially expressed between male and female gonads. We further analyzed the expression of TEs using the program SQuIRE and showed that more TE copies are overexpressed in testis than in ovaries (36% vs. 10%, respectively). We then developed a method to detect genomic regions enriched in testis- or ovary-biased genes. This revealed that sex-biased genes and TEs are not randomly distributed in the genome and a part of them form clusters with the same expression bias. We also found a correlation of expression between TE copies and their closest genes, which increases with decreasing intervening distance. Such a genomic organization suggests either that TEs hijack the regulatory sequences of neighboring sexual genes, allowing their expression in germ line cells and consequently new insertions to be transmitted to the next generation, or that TEs are involved in the regulation of sexual genes, and might therefore through their mobility participate in the rewiring of sex regulatory networks.

Hallmarks of Retroelement Expression in T-Cells Treated With HDAC Inhibitors

A wide spectrum of drugs have been assessed as latency reversal agents (LRA) to reactivate HIV-1 from cellular reservoirs and aid in viral eradication strategies. Histone deacetylase inhibitors (HDACi) have been studied in vitro and in vivo as potential candidates for HIV-1 latency reversion. Suberoylanilide hydroxamic acid (SAHA) and romidepsin (RMD) are two HDACi able to reverse HIV latency, however studies of potential off-target effects on retroelement expression have been limited. Retroelements constitute a large portion of the human genome, and some are considered “fossil viruses” as they constitute remnants of ancient exogenous retroviruses infections. Retroelements are reactivated during certain disease conditions like cancer or during HIV-1 infection. In this study, we analyzed differential expression of retroelements using publicly available RNA-seq datasets (GSE102187 and GSE114883) obtained from uninfected CD4+, and HIV-1 latently infected CD4+ T-cells treated with HDACi (SAHA and RMD). We found a total of 712 and 1,380 differentially expressed retroelements in HIV-1 latently infected cells following a 24-h SAHA and RMD treatment, respectively. Furthermore, we found that 531 retroelement sequences (HERVs and L1) were differentially expressed under both HDACi treatments, while 1,030 HERV/L1 were exclusively regulated by each drug. Despite differences in specific HERV loci expression, the overall pattern at the HERV family level was similar for both treatments. We detected differential expression of full-length HERV families including HERV-K, HERV-W and HERV-H. Furthermore, we analyzed the link between differentially expressed retroelements and nearby immune genes. TRAF2 (TNF receptor) and GBP5 (inflammasome activator) were upregulated in HDACi treated samples and their expression was correlated with nearby HERV (MERV101_9q34.3) and L1 (L1FLnI_1p22.2k, L1FLnI_1p22.2j, L1FLnI_1p22.2i). Our findings suggest that HDACi have an off-target effect on the expression of retroelements and on the expression of immune associated genes in treated CD4+ T-cells. Furthermore, our data highlights the importance of exploring the interaction between HIV-1 and retroelement expression in LRA treated samples to understand their role and impact on “shock and kill” strategies and their potential use as reservoir biomarkers.

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

One-carbon metabolism and global DNA methylation in mothers of individuals with Down syndrome

Down syndrome (DS) is the most common chromosomal disorder, resulting from the failure of normal chromosome 21 segregation. Studies have suggested that impairments within the one-carbon metabolic pathway can be of relevance for the global genome instability observed in mothers of individuals with DS. Based on the association between global DNA hypomethylation, genome instability, and impairments within the one-carbon metabolic pathway, the present study aimed to identify possible predictors, within the one-carbon metabolism, of global DNA methylation, measured by methylation patterns of LINE-1 and Alu repetitive sequences, in mothers of individuals with DS and mothers of individuals without the syndrome. In addition, we investigated one-carbon genetic polymorphisms and metabolites as maternal predisposing factors for the occurrence of trisomy 21 in children. Eighty-three samples of mothers of children with DS with karyotypically confirmed free trisomy 21 (case group) and 84 of mothers who had at least one child without DS or any other aneuploidy were included in the study. Pyrosequencing assays were performed to access global methylation. The results showed that group affiliation (case or control), betaine-homocysteine methyltransferase (BHMT) G742A and transcobalamin 2 (TCN2) C776G polymorphisms, and folate concentration were identified as predictors of global Alu DNA methylation values. In addition, thymidylate synthase (TYMS) 28-bp repeats 2R/3R or 3R/3R genotypes are independent maternal predisposing factors for having a child with DS. This study adds evidence that supports the association of impairments in the one-carbon metabolism, global DNA methylation, and the possibility of having a child with DS.

Transposable element sequence fragments incorporated into coding and noncoding transcripts modulate the transcriptome of human pluripotent stem cells

Transposable elements (TEs) occupy nearly 40% of mammalian genomes and, whilst most are fragmentary and no longer capable of transposition, they can nevertheless contribute to cell function. TEs within genes transcribed by RNA polymerase II can be copied as parts of primary transcripts; however, their full contribution to mature transcript sequences remains unresolved. Here, using long and short read (LR and SR) RNA sequencing data, we show that 26% of coding and 65% of noncoding transcripts in human pluripotent stem cells (hPSCs) contain TE-derived sequences. Different TE families are incorporated into RNAs in unique patterns, with consequences to transcript structure and function. The presence of TE sequences within a transcript is correlated with TE-type specific changes in its subcellular distribution, alterations in steady-state levels and half-life, and differential association with RNA Binding Proteins (RBPs). We identify hPSC-specific incorporation of endogenous retroviruses (ERVs) and LINE:L1 into protein-coding mRNAs, which generate TE sequence-derived peptides. Finally, single cell RNA-seq reveals that hPSCs express ERV-containing transcripts, whilst differentiating subpopulations lack ERVs and express SINE and LINE-containing transcripts. Overall, our comprehensive analysis demonstrates that the incorporation of TE sequences into the RNAs of hPSCs is more widespread and has a greater impact than previously appreciated.

The Simons Genome Diversity Project: A Global Analysis of Mobile Element Diversity

Ongoing retrotransposition of Alu, LINE-1, and SINE–VNTR–Alu elements generates diversity and variation among human populations. Previous analyses investigating the population genetics of mobile element insertions (MEIs) have been limited by population ascertainment bias or by relatively small numbers of populations and low sequencing coverage. Here, we use 296 individuals representing 142 global populations from the Simons Genome Diversity Project (SGDP) to discover and characterize MEI diversity from deeply sequenced whole-genome data. We report 5,742 MEIs not originally reported by the 1000 Genomes Project and show that high sampling diversity leads to a 4- to 7-fold increase in MEI discovery rates over the original 1000 Genomes Project data. As a result of negative selection, nonreference polymorphic MEIs are underrepresented within genes, and MEIs within genes are often found in the transcriptional orientation opposite that of the gene. Globally, 80% of Alu subfamilies predate the expansion of modern humans from Africa. Polymorphic MEIs show heterozygosity gradients that decrease from Africa to Eurasia to the Americas, and the number of MEIs found uniquely in a single individual are also distributed in this general pattern. The maximum fraction of MEI diversity partitioned among the seven major SGDP population groups (F_ST) is 7.4%, similar to, but slightly lower than, previous estimates and likely attributable to the diverse sampling strategy of the SGDP. Finally, we utilize these MEIs to extrapolate the primary Native American shared ancestry component to back to Asia and provide new evidence from genome-wide identical-by-descent genetic markers that add additional support for a southeastern Siberian origin for most Native Americans.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

The worldwide invasion of Drosophila suzukii is accompanied by a large increase of transposable element load and a small number of putatively adaptive insertions

Transposable Elements (TEs) are ubiquitous and mobile repeated sequences. They are major determinants of host fitness. Here, we characterized the TE content of the spotted wing fly Drosophila suzukii. Using a recently improved genome assembly, we reconstructed TE sequences de novo, and found that TEs occupy 47% of the genome and are mostly located in gene poor regions. The majority of TE insertions segregate at low frequencies, indicating a recent and probably ongoing TE activity. To explore TE dynamics in the context of biological invasions, we studied variation of TE abundance in genomic data from 16 invasive and six native populations of D. suzukii. We found a large increase of the TE load in invasive populations correlated with a reduced Watterson estimate of genetic diversity θ̂w a proxy of effective population size. We did not find any correlation between TE contents and bioclimatic variables, indicating a minor effect of environmentally induced TE activity. A genome-wide association study revealed that ca. 2,000 genomic regions are associated with TE abundance. We did not find, however, any evidence in such regions of an enrichment for genes known to interact with TE activity (e.g. transcription factor encoding genes or genes of the piRNA pathway). Finally, the study of TE insertion frequencies revealed 15 putatively adaptive TE insertions, six of them being likely associated with the recent invasion history of the species.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

Endogenous retroviruses in the origins and treatment of cancer

Endogenous retroviruses (ERVs) are emerging as promising therapeutic targets in cancer. As remnants of ancient retroviral infections, ERV-derived regulatory elements coordinate expression from gene networks, including those underpinning embryogenesis and immune cell function. ERV activation can promote an interferon response, a phenomenon termed viral mimicry. Although ERV expression is associated with cancer, and provisionally with autoimmune and neurodegenerative diseases, ERV-mediated inflammation is being explored as a way to sensitize tumors to immunotherapy. Here we review ERV co-option in development and innate immunity, the aberrant contribution of ERVs to tumorigenesis, and the wider biomedical potential of therapies directed at ERVs.

Expression of Human Endogenous Retroviruses in Systemic Lupus Erythematosus: Multiomic Integration With Gene Expression

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies predominantly to nuclear material. Many aspects of disease pathology are mediated by the deposition of nucleic acid containing immune complexes, which also induce the type 1interferon response, a characteristic feature of SLE. Notably, SLE is remarkably heterogeneous, with a variety of organs involved in different individuals, who also show variation in disease severity related to their ancestries. Here, we probed one potential contribution to disease heterogeneity as well as a possible source of immunoreactive nucleic acids by exploring the expression of human endogenous retroviruses (HERVs). We investigated the expression of HERVs in SLE and their potential relationship to SLE features and the expression of biochemical pathways, including the interferon gene signature (IGS). Towards this goal, we analyzed available and new RNA-Seq data from two independent whole blood studies using Telescope. We identified 481 locus specific HERV encoding regions that are differentially expressed between case and control individuals with only 14% overlap of differentially expressed HERVs between these two datasets. We identified significant differences between differentially expressed HERVs and non-differentially expressed HERVs between the two datasets. We also characterized the host differentially expressed genes and tested their association with the differentially expressed HERVs. We found that differentially expressed HERVs were significantly more physically proximal to host differentially expressed genes than non-differentially expressed HERVs. Finally, we capitalized on locus specific resolution of HERV mapping to identify key molecular pathways impacted by differential HERV expression in people with SLE.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes

One life-threatening outcome of cardiovascular disease is myocardial infarction, where cardiomyocytes are deprived of oxygen. To study inter-individual differences in response to hypoxia, we established an in vitro model of induced pluripotent stem cell-derived cardiomyocytes from 15 individuals. We measured gene expression levels, chromatin accessibility, and methylation levels in four culturing conditions that correspond to normoxia, hypoxia, and short- or long-term re-oxygenation. We characterized thousands of gene regulatory changes as the cells transition between conditions. Using available genotypes, we identified 1,573 genes with a cis expression quantitative locus (eQTL) in at least one condition, as well as 367 dynamic eQTLs, which are classified as eQTLs in at least one, but not in all conditions. A subset of genes with dynamic eQTLs is associated with complex traits and disease. Our data demonstrate how dynamic genetic effects on gene expression, which are likely relevant for disease, can be uncovered under stress.

Asymmetron: a toolkit for the identification of strand asymmetry patterns in biological sequences

DNA strand asymmetries can have a major effect on several biological functions, including replication, transcription and transcription factor binding. As such, DNA strand asymmetries and mutational strand bias can provide information about biological function. However, a versatile tool to explore this does not exist. Here, we present Asymmetron, a user-friendly computational tool that performs statistical analysis and visualizations for the evaluation of strand asymmetries. Asymmetron takes as input DNA features provided with strand annotation and outputs strand asymmetries for consecutive occurrences of a single DNA feature or between pairs of features. We illustrate the use of Asymmetron by identifying transcriptional and replicative strand asymmetries of germline structural variant breakpoints. We also show that the orientation of the binding sites of 45% of human transcription factors analyzed have a significant DNA strand bias in transcribed regions, that is also corroborated in ChIP-seq analyses, and is likely associated with transcription. In summary, we provide a novel tool to assess DNA strand asymmetries and show how it can be used to derive new insights across a variety of biological disciplines.

Transcript assembly improves expression quantification of transposable elements in single-cell RNA-seq data

Transposable elements (TEs) are an integral part of the host transcriptome. TE-containing noncoding RNAs (ncRNAs) show considerable tissue specificity and play important roles during development, including stem cell maintenance and cell differentiation. Recent advances in single-cell RNA-seq (scRNA-seq) revolutionized cell type-specific gene expression analysis. However, effective scRNA-seq quantification tools tailored for TEs are lacking, limiting our ability to dissect TE expression dynamics at single-cell resolution. To address this issue, we established a TE expression quantification pipeline that is compatible with scRNA-seq data generated across multiple technology platforms. We constructed TE-containing ncRNA references using bulk RNA-seq data and showed that quantifying TE expression at the transcript level effectively reduces noise. As proof of principle, we applied this strategy to mouse embryonic stem cells and successfully captured the expression profile of endogenous retroviruses in single cells. We further expanded our analysis to scRNA-seq data from early stages of mouse embryogenesis. Our results illustrated the dynamic TE expression at preimplantation stages and revealed 146 TE-containing ncRNA transcripts with substantial tissue specificity during gastrulation and early organogenesis.

A Field Guide to Eukaryotic Transposable Elements

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes. Through diverse invasion strategies, TEs have come to occupy a substantial fraction of nearly all eukaryotic genomes, and they represent a major source of genetic variation and novelty. Here we review the defining features of each major group of eukaryotic TEs and explore their evolutionary origins and relationships. We discuss how the unique biology of different TEs influences their propagation and distribution within and across genomes. Environmental and genetic factors acting at the level of the host species further modulate the activity, diversification, and fate of TEs, producing the dramatic variation in TE content observed across eukaryotes. We argue that cataloging TE diversity and dissecting the idiosyncratic behavior of individual elements are crucial to expanding our comprehension of their impact on the biology of genomes and the evolution of species.

Unveiling Human Non-Random Genome Editing Mechanisms Activated in Response to Chronic Environmental Changes: I. Where Might These Mechanisms Come from and What Might They Have Led To?

This article challenges the notion of the randomness of mutations in eukaryotic cells by unveiling stress-induced human non-random genome editing mechanisms. To account for the existence of such mechanisms, I have developed molecular concepts of the cell environment and cell environmental stressors and, making use of a large quantity of published data, hypothesised the origin of some crucial biological leaps along the evolutionary path of life on Earth under the pressure of natural selection, in particular, (1) virus-cell mating as a primordial form of sexual recombination and symbiosis; (2) Lamarckian CRISPR-Cas systems; (3) eukaryotic gene development; (4) antiviral activity of retrotransposon-guided mutagenic enzymes; and finally, (5) the exaptation of antiviral mutagenic mechanisms to stress-induced genome editing mechanisms directed at "hyper-transcribed" endogenous genes. Genes transcribed at their maximum rate (hyper-transcribed), yet still unable to meet new chronic environmental demands generated by "pollution", are inadequate and generate more and more intronic retrotransposon transcripts. In this scenario, RNA-guided mutagenic enzymes (e.g., Apolipoprotein B mRNA editing catalytic polypeptide-like enzymes, APOBECs), which have been shown to bind to retrotransposon RNA-repetitive sequences, would be surgically targeted by intronic retrotransposons on opened chromatin regions of the same "hyper-transcribed" genes. RNA-guided mutagenic enzymes may therefore "Lamarkianly" generate single nucleotide polymorphisms (SNP) and gene copy number variations (CNV), as well as transposon transposition and chromosomal translocations in the restricted areas of hyper-functional and inadequate genes, leaving intact the rest of the genome. CNV and SNP of hyper-transcribed genes may allow cells to surgically explore a new fitness scenario, which increases their adaptability to stressful environmental conditions. Like the mechanisms of immunoglobulin somatic hypermutation, non-random genome editing mechanisms may generate several cell mutants, and those codifying for the most environmentally adequate proteins would have a survival advantage and would therefore be Darwinianly selected. Non-random genome editing mechanisms represent tools of evolvability leading to organismal adaptation including transgenerational non-Mendelian gene transmission or to death of environmentally inadequate genomes. They are a link between environmental changes and biological novelty and plasticity, finally providing a molecular basis to reconcile gene-centred and "ecological" views of evolution.

Transposon clusters as substrates for aberrant splice-site activation

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5ʹ splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5ʹ splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

DNA methylation patterns of LINE-1 and Alu for pre-symptomatic dementia in type 2 diabetes

The identification of early markers of dementia is important for higher-risk populations such as those with type 2 diabetes (T2D). Retrotransposons, including long interspersed nuclear element 1 (LINE-1) and Alu, comprise ~40% of the human genome. Although dysregulation of these retrotransposons can induce aberrant gene regulation and genomic instability, their role in the development of pre-symptomatic dementia (PSD) among T2D patients is unknown. Here, we examined locus-specific changes in LINE-1 and Alu methylation in PSD and the potential to offset these changes via supplementation with folate and vitamin B₁₂. We interrogated DNA methylation patterns corresponding to 22,352 probes for LINE-1 and Alu elements using publicly-available Illumina Infinium 450K methylation datasets from i) an 18-month prospective study in 28 T2D patients (GSE62003) and ii) an intervention study in which 44 individuals were supplemented with folic acid (400 μg/day) and vitamin B₁₂ (500 μg/day) over two years (GSE74548). We identified 714 differentially methylated positions (DMP) mapping to retrotransposons in T2D patients who developed PSD in comparison to those who did not (P_FDR < 0.05), comprised of 2.4% (228 probes) of all LINE-1 probes and 3.8% (486 probes) of all Alu probes. These loci were enriched in genes with functions related to Alzheimer's disease and cognitive decline, including GNB5, GNG7 and PKN3 (p < 0.05). In older individuals supplemented with folate/vitamin B₁₂, 85 (11.9%) PSD retrotransposon loci showed significant changes in methylation (p < 0.05): participants with the MTHFR CC genotype predominantly showed hypermethylation at these loci, while hypomethylation was observed more frequently in those with the TT genotype. In T2D patients, LINE-1 and Alu elements are differentially methylated in PSD in a locus-specific manner and may offer clinical utility in monitoring risk of dementia. Further work is required to examine the potential for dietary supplementation in lowering the risk of PSD.

Variation in base composition underlies functional and evolutionary divergence in non-LTR retrotransposons

Background

Non-LTR retrotransposons often exhibit base composition that is markedly different from the nucleotide content of their host’s gene. For instance, the mammalian L1 element is AT-rich with a strong A bias on the positive strand, which results in a reduced transcription. It is plausible that the A-richness of mammalian L1 is a self-regulatory mechanism reflecting a trade-off between transposition efficiency and the deleterious effect of L1 on its host. We examined if the A-richness of L1 is a general feature of non-LTR retrotransposons or if different clades of elements have evolved different nucleotide content. We also investigated if elements belonging to the same clade evolved towards different base composition in different genomes or if elements from different clades evolved towards similar base composition in the same genome.

Results

We found that non-LTR retrotransposons differ in base composition among clades within the same host but also that elements belonging to the same clade differ in base composition among hosts. We showed that nucleotide content remains constant within the same host over extended period of evolutionary time, despite mutational patterns that should drive nucleotide content away from the observed base composition.

Conclusions

Our results suggest that base composition is evolving under selection and may be reflective of the long-term co-evolution between non-LTR retrotransposons and their host. Finally, the coexistence of elements with drastically different base composition suggests that these elements may be using different strategies to persist and multiply in the genome of their host.

Identification and characterization of ERV-W-like sequences in Platyrrhini species provides new insights into the evolutionary history of ERV-W in primates

Background

Endogenous Retroviruses (ERVs) constitute approximately 8% of every human genome and are relics of ancestral infections that affected the germ line cells. The ERV-W group contributed to primate physiology by providing an envelope protein (Syncytin-1) that has been adopted for placenta development in hominoids. Expression of Human ERV-W (HERV-W) sequences is investigated for a pathological role in various human diseases.

Results

We previously characterized ERV-W group genomic sequences in human and non-human Catarrhini species. We now investigated ERV-W-like sequences in the parvorder Platyrrhini, especially regarding two species with complete genome assemblies, namely marmoset (Callithrix jacchus) and squirrel monkey (Saimiri boliviensis). We identified in both species proviral sequences, annotated as ERV1–1 in respective genome assemblies, sharing high sequence similarities with Catarrhini ERV-W. A total of 130 relatively intact proviruses from the genomes of marmoset and squirrel monkey were characterized regarding their structural and evolutionarily relationships with Catarrhini ERV-W elements. Platyrrhini ERV-W sequences share several structural features with Catarrhini ERV-W elements and are closely related phylogenetically with the latter as well as with other ERV-W-related gammaretrovirus-like ERVs. The ERV-W group colonized Platyrrhini primates of both Callitrichidae and Atelidae lineages, with provirus formations having occurred mostly between 25 and 15 mya. Two LTR subgroups were associated with monophyletic proviral bodies. A pre-gag region appears to be a sequence feature common to the ERV-W group: it harbors a putative intron sequence that is missing in some ERV-W loci, holding a putative ORF as well. The presence of a long pre-gag portion was confirmed among all gammaretroviral ERV analyzed, suggesting a role in the latter biology. It is noteworthy that, contrary to Catarrhini ERV-W, there was no evidence of L1-mediated mobilization for Platyrrhini ERV-W sequences.

Conclusions

Our data establish that ERV-W is not exclusive to Catarrhini primates but colonized both parvorders of Simiiformes, providing further insight into the evolution of ERV-W and the colonization of primate genomes.

Retrotransposons shuttling genetic and epigenetic information from the nuclear to the mitochondrial compartment: Do they play a pathogenetic role in scleroderma?

Endogenous retroelements are a class of ancient defective viral insertions contained in the genome of host cells, where they account for up to 40% of all DNA. Centuries of co-existence in host genome have led to the development of immunotolerance to endogenous retroelements, most of which are defective and unable to replicate or transcribe functional proteins. However, given their capacity to move across the nuclear and mitochondrial genome and recombine, they could mix phenotypes and give rise to infections that may trigger innate and adaptive immune responses by sensing receptors capable of recognising foreign nucleic acids and proteins. It has recently been suggested that they play a role in the pathogenesis of autoimmune diseases on the grounds of their partial reactivation or the epigenetic control of host gene transcription. A number of studies have confirmed their contribution to the development of rheumatoid arthritis, multiple sclerosis and systemic lupus erythematosus, but there is still a lack of data concerning systemic sclerosis (SSc). Their role in the pathogenesis of SSc can be hypothesised on the basis of mitochondrial and nuclear chromatinic damage, and hyper-activation of the immune pathway involved in antiviral defense. SSc is characterised by genetic and immunological evidence of a viral infection but, as no viral agent has yet been isolated from SSc patients, the hypothesis that partial reactivation of endogenous retroviruses may trigger the disease cannot be excluded and deserves further investigation.

Endogenous retroviral insertions drive non-canonical imprinting in extra-embryonic tissues

BACKGROUND

Genomic imprinting is an epigenetic phenomenon that allows a subset of genes to be expressed mono-allelically based on the parent of origin and is typically regulated by differential DNA methylation inherited from gametes. Imprinting is pervasive in murine extra-embryonic lineages, and uniquely, the imprinting of several genes has been found to be conferred non-canonically through maternally inherited repressive histone modification H3K27me3. However, the underlying regulatory mechanisms of non-canonical imprinting in post-implantation development remain unexplored.

RESULTS

We identify imprinted regions in post-implantation epiblast and extra-embryonic ectoderm (ExE) by assaying allelic histone modifications (H3K4me3, H3K36me3, H3K27me3), gene expression, and DNA methylation in reciprocal C57BL/6 and CAST hybrid embryos. We distinguish loci with DNA methylation-dependent (canonical) and independent (non-canonical) imprinting by assaying hybrid embryos with ablated maternally inherited DNA methylation. We find that non-canonical imprints are localized to endogenous retrovirus-K (ERVK) long terminal repeats (LTRs), which act as imprinted promoters specifically in extra-embryonic lineages. Transcribed ERVK LTRs are CpG-rich and located in close proximity to gene promoters, and imprinting status is determined by their epigenetic patterning in the oocyte. Finally, we show that oocyte-derived H3K27me3 associated with non-canonical imprints is not maintained beyond pre-implantation development at these elements and is replaced by secondary imprinted DNA methylation on the maternal allele in post-implantation ExE, while being completely silenced by bi-allelic DNA methylation in the epiblast.

CONCLUSIONS

This study reveals distinct epigenetic mechanisms regulating non-canonical imprinted gene expression between embryonic and extra-embryonic development and identifies an integral role for ERVK LTR repetitive elements.

SQuIRE reveals locus-specific regulation of interspersed repeat expression

Transposable elements (TEs) are interspersed repeat sequences that make up much of the human genome. Their expression has been implicated in development and disease. However, TE-derived RNA-seq reads are difficult to quantify. Past approaches have excluded these reads or aggregated RNA expression to subfamilies shared by similar TE copies, sacrificing quantitative accuracy or the genomic context necessary to understand the basis of TE transcription. As a result, the effects of TEs on gene expression and associated phenotypes are not well understood. Here, we present Software for Quantifying Interspersed Repeat Expression (SQuIRE), the first RNA-seq analysis pipeline that provides a quantitative and locus-specific picture of TE expression (https://github.com/wyang17/SQuIRE). SQuIRE is an accurate and user-friendly tool that can be used for a variety of species. We applied SQuIRE to RNA-seq from normal mouse tissues and a Drosophila model of amyotrophic lateral sclerosis. In both model organisms, we recapitulated previously reported TE subfamily expression levels and revealed locus-specific TE expression. We also identified differences in TE transcription patterns relating to transcript type, gene expression and RNA splicing that would be lost with other approaches using subfamily-level analyses. Altogether, our findings illustrate the importance of studying TE transcription with locus-level resolution.

MIR sequences recruit zinc finger protein ZNF768 to expressed genes

Mammalian-wide interspersed repeats (MIRs) are retrotransposed elements of mammalian genomes. Here, we report the specific binding of zinc finger protein ZNF768 to the sequence motif GCTGTGTG (N20) CCTCTCTG in the core region of MIRs. ZNF768 binding is preferentially associated with euchromatin and promoter regions of genes. Binding was observed for genes expressed in a cell type-specific manner in human B cell line Raji and osteosarcoma U2OS cells. Mass spectrometric analysis revealed binding of ZNF768 to Elongator components Elp1, Elp2 and Elp3 and other nuclear factors. The N-terminus of ZNF768 contains a heptad repeat array structurally related to the C-terminal domain (CTD) of RNA polymerase II. This array evolved in placental animals but not marsupials and monotreme species, displays species-specific length variations, and possibly fulfills CTD related functions in gene regulation. We propose that the evolution of MIRs and ZNF768 has extended the repertoire of gene regulatory mechanisms in mammals and that ZNF768 binding is associated with cell type-specific gene expression.

Retrotransposons evolution and impact on lncRNA and protein coding genes in pigs

BACKGROUND

Retrotransposons are the major determinants of genome sizes and they have shaped both genes and genomes in mammalian organisms, but their overall activity, diversity, and evolution dynamics, particularly their impact on protein coding and lncRNA genes in pigs remain largely unknown.

RESULTS

In the present study, we performed de novo detection of retrotransposons in pigs by using multiple pipelines, four distinct families of pig-specific L1 s classified into 51 distinct subfamilies and representing four evolution models and three expansion waves of pig-specific SINEs represented by three distinct families were identified. ERVs were classified into 18 families and found two most "modern" subfamilies in the pig genome. The transposition activity of pig L1 was verified by experiment, the sense and antisense promoter activities of young L1 5'UTRs and ERV LTRs and expression profiles of young retrotransposons in multiple tissues and cell lines were also validated. Furthermore, retrotransposons had an extensive impact on lncRNA and protein coding genes at both the genomic and transcriptomic levels. Most protein coding and lncRNA (> 80%) genes contained retrotransposon insertions, and about half of protein coding genes (44.30%) and one-fourth (24.13%) of lncRNA genes contained the youngest retrotransposon insertions. Nearly half of protein coding genes (43.78%) could generate chimeric transcripts with retrotransposons. Significant distribution bias of retrotransposon composition, location, and orientation in lncRNA and protein coding genes, and their transcripts, were observed.

CONCLUSIONS

In the current study, we characterized the classification and evolution profile of retrotransposons in pigs, experimentally proved the transposition activity of the young pig L1 subfamily, characterized the sense and antisense expression profiles and promoter activities of young retrotransposons, and investigated their impact on lncRNA and protein coding genes by defining the mobilome landscapes at the genomic and transcriptomic levels. These findings help provide a better understanding of retrotransposon evolution in mammal and their impact on the genome and transcriptome.

Mouse germ line mutations due to retrotransposon insertions

Transposable element (TE) insertions are responsible for a significant fraction of spontaneous germ line mutations reported in inbred mouse strains. This major contribution of TEs to the mutational landscape in mouse contrasts with the situation in human, where their relative contribution as germ line insertional mutagens is much lower. In this focussed review, we provide comprehensive lists of TE-induced mouse mutations, discuss the different TE types involved in these insertional mutations and elaborate on particularly interesting cases. We also discuss differences and similarities between the mutational role of TEs in mice and humans.

The Landscape of L1 Retrotransposons in the Human Genome Is Shaped by Pre-insertion Sequence Biases and Post-insertion Selection

L1 retrotransposons are transposable elements and major contributors of genetic variation in humans. Where L1 integrates into the genome can directly impact human evolution and disease. Here, we experimentally induced L1 retrotransposition in cells and mapped integration sites at nucleotide resolution. At local scales, L1 integration is mostly restricted by genome sequence biases and the specificity of the L1 machinery. At regional scales, L1 shows a broad capacity for integration into all chromatin states, in contrast to other known mobile genetic elements. However, integration is influenced by the replication timing of target regions, suggesting a link to host DNA replication. The distribution of new L1 integrations differs from those of preexisting L1 copies, which are significantly reshaped by natural selection. Our findings reveal that the L1 machinery has evolved to efficiently target all genomic regions and underline a predominant role for post-integrative processes on the distribution of endogenous L1 elements.