Transposable elements: an abundant and natural source of regulatory sequences for host genes

Transposon-mediated genic rearrangements underlie variation in small RNA pathways

Transposable elements (TEs) can alter host gene structure and expression, whereas host organisms develop mechanisms to repress TE activities. In the nematode Caenorhabditis elegans, a small interfering RNA pathway dependent on the helicase ERI-6/7 primarily silences retrotransposons and recent genes of likely viral origin. By studying gene expression variation among wild C. elegans strains, we found that structural variants and transposon remnants likely underlie expression variation in eri-6/7 and the pathway targets. We further found that multiple insertions of the DNA transposons, Polintons, reshuffled the eri-6/7 locus and induced inversion of eri-6 in some wild strains. In the inverted configuration, gene function was previously shown to be repaired by unusual trans-splicing mediated by direct repeats. We identified that these direct repeats originated from terminal inverted repeats of Polintons. Our findings highlight the role of host-transposon interactions in driving rapid host genome diversification among natural populations and shed light on evolutionary novelty in genes and splicing mechanisms.

Multi-ancestry GWAS reveals loci linked to human variation in LINE-1- and Alu-copy numbers

Long INterspersed Element-1 (LINE-1; L1) and Alu are two families of transposable elements (TEs) occupying ~17% and ~11% of the human genome, respectively. Though only a small fraction of L1 copies is able to produce the machinery to mobilize autonomously, Alu elements and degenerate L1 copies can hijack their functional machinery and mobilize in trans. The expression and subsequent copy number expansion of L1 and Alu can exert pathological effects on their hosts, promoting genome instability, inflammation, and cell cycle alterations. These features have made L1 and Alu promising focus subjects in studies of aging and aging diseases where they can become active. However, the mechanisms regulating variation in their expression and copy number remain incompletely characterized. Moreover, the relevance of known mechanisms to diverse human populations remains unclear, as mechanisms are often characterized in isogenic cell culture models. To address these gaps, we leveraged genomic data from the 1000 Genomes Project to carry out a trans-ethnic GWAS of L1 and Alu insertion global singletons. These singletons are rare insertions observed only once in a population, potentially reflecting recently acquired L1 and Alu integrants or structural variants, and which we used as proxies for L1/Alu-associated copy number variation. Our computational approach identified single nucleotide variants in genomic regions containing genes with potential and known TE regulatory properties, and it enriched for single nucleotide variants in regions containing known regulators of L1 expression. Moreover, we identified many reference TE copies and polymorphic structural variants that were associated with L1/Alu singletons, suggesting their potential contribution to TE copy number variation through transposition-dependent or transposition-independent mechanisms. Finally, a transcriptional analysis of lymphoblastoid cells highlighted potential cell cycle alterations in a subset of samples harboring L1/Alu singletons. Collectively, our results (i) suggest that known TE regulatory mechanisms may also play regulatory roles in diverse human populations, (ii) expand the list of genic and repetitive genomic loci implicated in TE copy number variation, and (iii) reinforce the links between TEs and disease.

Regulatory transposable elements in the encyclopedia of DNA elements

Transposable elements (TEs) comprise ~50% of our genome, but knowledge of how TEs affect genome evolution remains incomplete. Leveraging ENCODE4 data, we provide the most comprehensive study to date of TE contributions to the regulatory genome. We find 236,181 (~25%) human candidate cis-regulatory elements (cCREs) are TE-derived, with over 90% lineage-specific since the human-mouse split, accounting for 8-36% of lineage-specific cCREs. Except for SINEs, cCRE-associated transcription factor (TF) motifs in TEs are derived from ancestral TE sequence more than expected by chance. We show that TEs may adopt similar regulatory activities of elements near their integration site. Since human-mouse divergence, TEs have contributed 3-56% of TF binding site turnover events across 30 examined TFs. Finally, TE-derived cCREs are similar to non-TE cCREs in terms of MPRA activity and GWAS variant enrichment. Overall, our results substantiate the notion that TEs have played an important role in shaping the human regulatory genome.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Mutant IDH1 inhibition induces dsDNA sensing to activate tumor immunity

Isocitrate dehydrogenase 1 (IDH1) is the most commonly mutated metabolic gene across human cancers. Mutant IDH1 (mIDH1) generates the oncometabolite (R)-2-hydroxyglutarate, disrupting enzymes involved in epigenetics and other processes. A hallmark of IDH1-mutant solid tumors is T cell exclusion, whereas mIDH1 inhibition in preclinical models restores antitumor immunity. Here, we define a cell-autonomous mechanism of mIDH1-driven immune evasion. IDH1-mutant solid tumors show selective hypermethylation and silencing of the cytoplasmic double-stranded DNA (dsDNA) sensor CGAS, compromising innate immune signaling. mIDH1 inhibition restores DNA demethylation, derepressing CGAS and transposable element (TE) subclasses. dsDNA produced by TE-reverse transcriptase (TE-RT) activates cGAS, triggering viral mimicry and stimulating antitumor immunity. In summary, we demonstrate that mIDH1 epigenetically suppresses innate immunity and link endogenous RT activity to the mechanism of action of a US Food and Drug Administration-approved oncology drug.

Schistosoma haematobium and Schistosoma bovis first generation hybrids undergo gene expressions changes consistent with species compatibility and heterosis

When two species hybridize, the two parental genomes are brought together and some alleles might interact for the first time. To date, the extent of the transcriptomic changes in first hybrid generations, along with their functional outcome constitute an important knowledge gap, especially in parasite species. Here we explored the molecular and functional outcomes of hybridization in first-generation hybrids between the blood fluke parasites Schistosoma haematobium and S. bovis. Through a transcriptomic approach, we measured gene expression in both parental species and hybrids. We described and quantified expression profiles encountered in hybrids along with the main biological processes impacted. Up to 7,100 genes fell into a particular hybrid expression profile (intermediate between the parental expression levels, over-expressed, under-expressed, or expressed like one of the parental lines). Most of these genes were different depending on the direction of the parental cross (S. bovis mother and S. haematobium father or the reverse) and depending on the sex. For a given sex and cross direction, the vast majority of genes were hence unassigned to a hybrid expression profile: either they were differentially expressed genes but not typical of any hybrid expression profiles or they were not differentially expressed neither between hybrids and parental lines nor between parental lines. The most prevalent profile of gene expression in hybrids was the intermediate one (24% of investigated genes). These results suggest that transcriptomic compatibility between S. haematobium and S. bovis remains quite high. We also found support for an over-dominance model (over- and under-expressed genes in hybrids compared to parental lines) potentially associated with heterosis. In females in particular, processes such as reproductive processes, metabolism and cell interactions as well as signaling pathways were indeed affected. Our study hence provides new insight on the biology of Schistosoma hybrids with evidences supporting compatibility and heterosis.

Characterization of the CsCENH3 protein and centromeric DNA profiles reveal the structures of centromeres in cucumber

Centromeres in eukaryotes mediate the accurate segregation of chromosomes during cell division. They serve as essential functional units of chromosomes and play a core role in the process of genome evolution. Centromeres are composed of satellite repeats and highly repetitive centromeric retrotransposons (CRs), which vary greatly even among closely related species. Cucumber (Cucumis sativus) is a globally cultivated and economically important vegetable and the only species in the Cucumis genus with seven pairs of chromosomes. Therefore, studying the centromeres of the Cucumis subgenus may yield valuable insights into its genome structure and evolution. Using chromatin immunoprecipitation (ChIP) techniques, we isolated centromeric DNA from cucumber reference line 9930. Our investigation into cucumber centromeres uncovered the centromeric satellite sequence, designated as CentCs, and the prevalence of Ty1/Copia long terminal repeat retrotransposons. In addition, active genes were identified in the CsCENH3 nucleosome regions with low transcription levels. To the best of our knowledge, this is the first time that characterization of centromeres has been achieved in cucumber. Meanwhile, our results on the distribution of CentCs and CsCRs in the subgenus Cucumis indicate that the content of centromeric repeats in the wild variants was significantly reduced compared with the cultivated cucumber. The results provide evidence for centromeric DNA amplification that occurred during the domestication process from wild to cultivated cucumber. Furthermore, these findings may offer new information for enhancing our understanding of phylogenetic relationships in the Cucumis genus.

An eQTL-based approach reveals candidate regulators of LINE-1 RNA levels in lymphoblastoid cells

Long interspersed element 1 (LINE-1; L1) are a family of transposons that occupy ~17% of the human genome. Though a small number of L1 copies remain capable of autonomous transposition, the overwhelming majority of copies are degenerate and immobile. Nevertheless, both mobile and immobile L1s can exert pleiotropic effects (promoting genome instability, inflammation, or cellular senescence) on their hosts, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSD17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for regulators of transposon RNA levels.

Molecular Mimicry of Transposable Elements in Plants

Transposable elements (TEs) are mobile DNA elements that are particularly abundant in the plant genomes. They have long been considered as junk DNA; however, a growing body of evidence suggests that TE insertions promote genetic diversity that is essential for the adaptive evolution of a species. Thus far, studies have mainly investigated the cis-acting regulatory roles of TEs generated by their insertions nearby or within the host genes. However, the trans-acting effects of TE-derived RNA and DNA remained obscure to date. TEs contain various regulatory elements within their sequences that can accommodate the binding of specific RNAs and proteins. Recently, it was suggested that some of these cellular regulators are shared between TEs and the host genes, and the competition for the common host factors underlies the fine-tuned developmental reprogramming. In this review, we will highlight and discuss the latest discoveries on the biological functions of plant TEs, with a particular focus on their competitive binding with specific developmental regulators.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Transposable Element Expression Profiles in Premalignant Pigment Cell Lesions and Melanoma of Xiphophorus

Transposable elements (TEs) are characterized by their ability to change their genomic position. Through insertion or recombination leading to deletions and other chromosomal aberrations, they can cause genetic instability. The extent to which they thereby exert regulatory influence on cellular functions is unclear. To better characterize TEs in processes such as carcinogenesis, we used the well-established Xiphophorus melanoma model. By transcriptome sequencing, we show that an increasing total number in transposons correlates with progression of malignancy in melanoma samples from Xiphophorus interspecific hybrids. Further, by comparing the presence of TEs in the parental genomes of Xiphophorus maculatus and Xiphophorus hellerii, we could show that even in closely related species, genomic location and spectrum of TEs are considerably different.

piRNA-Guided Transposon Silencing and Response to Stress in Drosophila Germline

Transposons are integral genome constituents that can be domesticated for host functions, but they also represent a significant threat to genome stability. Transposon silencing is especially critical in the germline, which is dedicated to transmitting inherited genetic material. The small Piwi-interacting RNAs (piRNAs) have a deeply conserved function in transposon silencing in the germline. piRNA biogenesis and function are particularly well understood in Drosophila melanogaster, but some fundamental mechanisms remain elusive and there is growing evidence that the pathway is regulated in response to genotoxic and environmental stress. Here, we review transposon regulation by piRNAs and the piRNA pathway regulation in response to stress, focusing on the Drosophila female germline.

Look4LTRs: a Long terminal repeat retrotransposon detection tool capable of cross species studies and discovering recently nested repeats

Plant genomes include large numbers of transposable elements. One particular type of these elements is flanked by two Long Terminal Repeats (LTRs) and can translocate using RNA. Such elements are known as LTR-retrotransposons; they are the most abundant type of transposons in plant genomes. They have many important functions involving gene regulation and the rise of new genes and pseudo genes in response to severe stress. Additionally, LTR-retrotransposons have several applications in biotechnology. Due to the abundance and the importance of LTR-retrotransposons, multiple computational tools have been developed for their detection. However, none of these tools take advantages of the availability of related genomes; they process one chromosome at a time. Further, recently nested LTR-retrotransposons (multiple elements of the same family are inserted into each other) cannot be annotated accurately - or cannot be annotated at all - by the currently available tools. Motivated to overcome these two limitations, we built Look4LTRs, which can annotate LTR-retrotransposons in multiple related genomes simultaneously and discover recently nested elements. The methodology of Look4LTRs depends on techniques imported from the signal-processing field, graph algorithms, and machine learning with a minimal use of alignment algorithms. Four plant genomes were used in developing Look4LTRs and eight plant genomes for evaluating it in contrast to three related tools. Look4LTRs is the fastest while maintaining better or comparable F1 scores (the harmonic average of recall and precision) to those obtained by the other tools. Our results demonstrate the added benefit of annotating LTR-retrotransposons in multiple related genomes simultaneously and the ability to discover recently nested elements. Expert human manual examination of six elements - not included in the ground truth - revealed that three elements belong to known families and two elements are likely from new families. With respect to examining recently nested LTR-retrotransposons, three out of five were confirmed to be valid elements. Look4LTRs - with its speed, accuracy, and novel features - represents a true advancement in the annotation of LTR-retrotransposons, opening the door to many studies focused on understanding their functions in plants.

Precise genome-editing in human diseases: mechanisms, strategies and applications

Precise genome-editing platforms are versatile tools for generating specific, site-directed DNA insertions, deletions, and substitutions. The continuous enhancement of these tools has led to a revolution in the life sciences, which promises to deliver novel therapies for genetic disease. Precise genome-editing can be traced back to the 1950s with the discovery of DNA's double-helix and, after 70 years of development, has evolved from crude in vitro applications to a wide range of sophisticated capabilities, including in vivo applications. Nonetheless, precise genome-editing faces constraints such as modest efficiency, delivery challenges, and off-target effects. In this review, we explore precise genome-editing, with a focus on introduction of the landmark events in its history, various platforms, delivery systems, and applications. First, we discuss the landmark events in the history of precise genome-editing. Second, we describe the current state of precise genome-editing strategies and explain how these techniques offer unprecedented precision and versatility for modifying the human genome. Third, we introduce the current delivery systems used to deploy precise genome-editing components through DNA, RNA, and RNPs. Finally, we summarize the current applications of precise genome-editing in labeling endogenous genes, screening genetic variants, molecular recording, generating disease models, and gene therapy, including ex vivo therapy and in vivo therapy, and discuss potential future advances.

Transposable elements as scaffold/matrix attachment regions: shaping organization and functions in genomes

The hierarchical structure of eukaryotic genomes has regulatory layers, one of them being epigenetic "indexing" of the genome that leads to cell-type-specific patterns of gene expression. By establishing loops and defining chromatin domains, cells can achieve coordinated control over multi-locus segments of the genome. This is thought to be achieved using scaffold/matrix attachment regions (S/MARs) that establish structural and functional loops and topologically associating domains (TADs) that define a self-interacting region of the genome. Large-scale genome-wide mapping of S/MARs has begun to uncover these aspects of genome organization. A recent genome-wide study showed the association of transposable elements (TEs) with a significant fraction of S/MARs, suggesting that the multitude of TE-derived repeats constitute a class of anchorage sites of chromatin loops to nuclear architecture. In this study, we provide an insight that TE-driven dispersal of S/MARs has the potential to restructure the chromosomes by creating novel loops and domains. The combination of TEs and S/MARs, as elements that can hop through the genome along with regulatory capabilities, may provide an active mechanism of genome evolution leading to the emergence of novel features in biological systems. The significance is that a genome-wide study mapping developmental S/MARs reveals an intriguing link between these elements and TEs. This article highlights the potential of the TE-S/MAR combination to drive evolution by restructuring and shaping the genome.

Towards targeting transposable elements for cancer therapy

Transposable elements (TEs) represent almost half of the human genome. Historically deemed 'junk DNA', recent technological advancements have stimulated a wave of research into the functional impact of TEs on gene-regulatory networks in evolution and development, as well as in diseases including cancer. The genetic and epigenetic evolution of cancer involves the exploitation of TEs, whereby TEs contribute directly to cancer-specific gene activities. This Review provides a perspective on the role of TEs in cancer as being a 'double-edged sword', both promoting cancer evolution and representing a vulnerability that could be exploited in cancer therapy. We discuss how TEs affect transcriptome regulation and other cellular processes in cancer. We highlight the potential of TEs as therapeutic targets for cancer. We also summarize technical hurdles in the characterization of TEs with genomic assays. Last, we outline open questions and exciting future research avenues.

Transposable elements mediate genetic effects altering the expression of nearby genes in colorectal cancer

Transposable elements (TEs) are prevalent repeats in the human genome, play a significant role in the regulome, and their disruption can contribute to tumorigenesis. However, TE influence on gene expression in cancer remains unclear. Here, we analyze 275 normal colon and 276 colorectal cancer samples from the SYSCOL cohort, discovering 10,231 and 5,199 TE-expression quantitative trait loci (eQTLs) in normal and tumor tissues, respectively, of which 376 are colorectal cancer specific eQTLs, likely due to methylation changes. Tumor-specific TE-eQTLs show greater enrichment of transcription factors, compared to shared TE-eQTLs suggesting specific regulation of their expression in tumor. Bayesian networks reveal 1,766 TEs as mediators of genetic effects, altering the expression of 1,558 genes, including 55 known cancer driver genes and show that tumor-specific TE-eQTLs trigger the driver capability of TEs. These insights expand our knowledge of cancer drivers, deepening our understanding of tumorigenesis and presenting potential avenues for therapeutic interventions.

Advances in Genetic Reprogramming: Prospects from Developmental Biology to Regenerative Medicine

The foundations of cell reprogramming were laid by Yamanaka and co-workers, who showed that somatic cells can be reprogrammed into pluripotent cells (induced pluripotency). Since this discovery, the field of regenerative medicine has seen advancements. For example, because they can differentiate into multiple cell types, pluripotent stem cells are considered vital components in regenerative medicine aimed at the functional restoration of damaged tissue. Despite years of research, both replacement and restoration of failed organs/ tissues have remained elusive scientific feats. However, with the inception of cell engineering and nuclear reprogramming, useful solutions have been identified to counter the need for compatible and sustainable organs. By combining the science underlying genetic engineering and nuclear reprogramming with regenerative medicine, scientists have engineered cells to make gene and stem cell therapies applicable and effective. These approaches have enabled the targeting of various pathways to reprogramme cells, i.e., make them behave in beneficial ways in a patient-specific manner. Technological advancements have clearly supported the concept and realization of regenerative medicine. Genetic engineering is used for tissue engineering and nuclear reprogramming and has led to advances in regenerative medicine. Targeted therapies and replacement of traumatized , damaged, or aged organs can be realized through genetic engineering. Furthermore, the success of these therapies has been validated through thousands of clinical trials. Scientists are currently evaluating induced tissue-specific stem cells (iTSCs), which may lead to tumour-free applications of pluripotency induction. In this review, we present state-of-the-art genetic engineering that has been used in regenerative medicine. We also focus on ways that genetic engineering and nuclear reprogramming have transformed regenerative medicine and have become unique therapeutic niches.

Comprehensive identification of onco-exaptation events in bladder cancer cell lines revealed L1PA2-SYT1 as a prognosis-relevant event

Transposable elements (TEs) can provide ectopic promoters to drive the expression of oncogenes in cancer, a mechanism known as onco-exaptation. Onco-exaptation events have been extensively identified in various cancers, with bladder cancer showing a high frequency of onco-exaptation events (77%). However, the effect of most of these events in bladder cancer remains unclear. This study identified 44 onco-exaptation events in 44 bladder cancer cell lines in 137 RNA-seq datasets from six publicly available cohorts, with L1PA2 contributing the most events. L1PA2-SYT1, L1PA2-MET, and L1PA2-XCL1 had the highest frequency not only in cell lines but also in TCGA-BLCA samples. L1PA2-SYT1 showed significant tumor specificity and was found to be activated by CpG island demethylation in its promoter. The upregulation of L1PA2-SYT1 enhances the in vitro invasion of bladder cancer and is an independent risk factor for patient's overall survival, suggesting L1PA2-SYT1 being an important event that promotes the development of bladder cancer.

m6A RNA demethylase AtALKBH9B promotes mobilization of a heat-activated long terminal repeat retrotransposon in Arabidopsis

Transposons are mobile and ubiquitous DNA molecules that can cause vast genomic alterations. In plants, it is well documented that transposon mobilization is strongly repressed by DNA methylation; however, its regulation at the posttranscriptional level remains relatively uninvestigated. Here, we suggest that transposon RNA is marked by m6A RNA methylation and can be localized in stress granules (SGs). Intriguingly, SG-localized AtALKBH9B selectively demethylates a heat-activated retroelement, Onsen, and thereby releases it from spatial confinement, allowing for its mobilization. In addition, we show evidence that m6A RNA methylation contributes to transpositional suppression by inhibiting virus-like particle assembly and extrachromosomal DNA production. In summary, this study unveils a previously unknown role for m6A in the suppression of transposon mobility and provides insight into how transposons counteract the m6A-mediated repression mechanism by hitchhiking the RNA demethylase of the host.

Advances in transposable elements: from mechanisms to applications in mammalian genomics

It has been 70 years since Barbara McClintock discovered transposable elements (TE), and the mechanistic studies and functional applications of transposable elements have been at the forefront of life science research. As an essential part of the genome, TEs have been discovered in most species of prokaryotes and eukaryotes, and the relative proportion of the total genetic sequence they comprise gradually increases with the expansion of the genome. In humans, TEs account for about 40% of the genome and are deeply involved in gene regulation, chromosome structure maintenance, inflammatory response, and the etiology of genetic and non-genetic diseases. In-depth functional studies of TEs in mammalian cells and the human body have led to a greater understanding of these fundamental biological processes. At the same time, as a potent mutagen and efficient genome editing tool, TEs have been transformed into biological tools critical for developing new techniques. By controlling the random insertion of TEs into the genome to change the phenotype in cells and model organisms, critical proteins of many diseases have been systematically identified. Exploiting the TE's highly efficient in vitro insertion activity has driven the development of cutting-edge sequencing technologies. Recently, a new technology combining CRISPR with TEs was reported, which provides a novel targeted insertion system to both academia and industry. We suggest that interrogating biological processes that generally depend on the actions of TEs with TEs-derived genetic tools is a very efficient strategy. For example, excessive activation of TEs is an essential factor in the occurrence of cancer in humans. As potent mutagens, TEs have also been used to unravel the key regulatory elements and mechanisms of carcinogenesis. Through this review, we aim to effectively combine the traditional views of TEs with recent research progress, systematically link the mechanistic discoveries of TEs with the technological developments of TE-based tools, and provide a comprehensive approach and understanding for researchers in different fields.

An eQTL-based Approach Reveals Candidate Regulators of LINE-1 RNA Levels in Lymphoblastoid Cells

Long interspersed element 1 (L1) are a family of autonomous, actively mobile transposons that occupy ~17% of the human genome. A number of pleiotropic effects induced by L1 (promoting genome instability, inflammation, or cellular senescence) have been observed, and L1's contributions to aging and aging diseases is an area of active research. However, because of the cell type-specific nature of transposon control, the catalogue of L1 regulators remains incomplete. Here, we employ an eQTL approach leveraging transcriptomic and genomic data from the GEUVADIS and 1000Genomes projects to computationally identify new candidate regulators of L1 RNA levels in lymphoblastoid cell lines. To cement the role of candidate genes in L1 regulation, we experimentally modulate the levels of top candidates in vitro, including IL16, STARD5, HSDB17B12, and RNF5, and assess changes in TE family expression by Gene Set Enrichment Analysis (GSEA). Remarkably, we observe subtle but widespread upregulation of TE family expression following IL16 and STARD5 overexpression. Moreover, a short-term 24-hour exposure to recombinant human IL16 was sufficient to transiently induce subtle, but widespread, upregulation of L1 subfamilies. Finally, we find that many L1 expression-associated genetic variants are co-associated with aging traits across genome-wide association study databases. Our results expand the catalogue of genes implicated in L1 RNA control and further suggest that L1-derived RNA contributes to aging processes. Given the ever-increasing availability of paired genomic and transcriptomic data, we anticipate this new approach to be a starting point for more comprehensive computational scans for transposon transcriptional regulators.

Human endogenous retroviruses of the HERV-K (HML-2) family are expressed in the brain of healthy individuals and modify the composition of the brain-infiltrating immune cells

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the human genome. RNA expression of individual HERVs has frequently been observed in various pathologic conditions, but some activity can also be seen in healthy individuals, e.g. in the blood. To quantitate the basal expression levels of HERVs in the brain, we now used high-throughput sequencing-based metagenomic analysis to characterize the expression profiles of the HERV-K (HML-2) family proviruses in different brain regions of healthy brain tissue. To this end, RNA-seq data from the Genotype-Tissue Expression (GTEx) project was used. The GTEx project is a public resource to study tissue-specific gene expression and regulation, consisting of a large selection of sequenced samples from different tissues. The GTEx data used in this study consisted of 378 samples taken from 13 brain regions from 55 individuals. The data demonstrated that out of 99 intact proviruses in the family 58 were expressed, but the expression profiles were highly divergent and there were no significant differences in the expression profiles between the various anatomic regions of the brain. It is known that the brain contains a variety of infiltrating immune cells, which are probably of great importance both in the normal defense mechanisms as well as in the various pathogenic processes. Digital cytometry (CIBERSORTx) was used to quantify the proportions of the infiltrating immune cells in the same brain samples. Six most abundant (>5 % of the total population) cell types were observed to be CD4 memory resting T cells, M0 macrophages, plasma cells, CD8 T cells, CD4 memory activated T cells, and monocytes. Analysis of the correlations between the individual HERVs and infiltrating cell types indicated that a cluster of 6 HERVs had a notable correlation signature between T cell type infiltrating cell proportions and HERV RNA expression intensity. The correlations between inflammatory type infiltrating cells were negative or weak. Taken together, these data indicate that the expression of HERVs is associated with a T cell type immunity.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Genome-wide enhancer identification by massively parallel reporter assay in Arabidopsis

Enhancers are critical cis-regulatory elements controlling gene expression during cell development and differentiation. However, genome-wide enhancer characterization has been challenging due to the lack of a well-defined relationship between enhancers and genes. Function-based methods are the gold standard for determining the biological function of cis-regulatory elements; however, these methods have not been widely applied to plants. Here, we applied a massively parallel reporter assay on Arabidopsis to measure enhancer activities across the genome. We identified 4327 enhancers with various combinations of epigenetic modifications distinctively different from animal enhancers. Furthermore, we showed that enhancers differ from promoters in their preference for transcription factors. Although some enhancers are not conserved and overlap with transposable elements forming clusters, enhancers are generally conserved across thousand Arabidopsis accessions, suggesting they are selected under evolution pressure and could play critical roles in the regulation of important genes. Moreover, comparison analysis reveals that enhancers identified by different strategies do not overlap, suggesting these methods are complementary in nature. In sum, we systematically investigated the features of enhancers identified by functional assay in A. thaliana, which lays the foundation for further investigation into enhancers' functional mechanisms in plants.

Genomic object detection: An improved approach for transposable elements detection and classification using convolutional neural networks

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO.

A Quantitative, Genome-Wide Analysis in Drosophila Reveals Transposable Elements' Influence on Gene Expression Is Species-Specific

Transposable elements (TEs) are parasite DNA sequences that are able to move and multiply along the chromosomes of all genomes. They can be controlled by the host through the targeting of silencing epigenetic marks, which may affect the chromatin structure of neighboring sequences, including genes. In this study, we used transcriptomic and epigenomic high-throughput data produced from ovarian samples of several Drosophila melanogaster and Drosophila simulans wild-type strains, in order to finely quantify the influence of TE insertions on gene RNA levels and histone marks (H3K9me3 and H3K4me3). Our results reveal a stronger epigenetic effect of TEs on ortholog genes in D. simulans compared with D. melanogaster. At the same time, we uncover a larger contribution of TEs to gene H3K9me3 variance within genomes in D. melanogaster, which is evidenced by a stronger correlation of TE numbers around genes with the levels of this chromatin mark in D. melanogaster. Overall, this work contributes to the understanding of species-specific influence of TEs within genomes. It provides a new light on the considerable natural variability provided by TEs, which may be associated with contrasted adaptive and evolutionary potentials.

The transcriptional landscape of endogenous retroelements delineates esophageal adenocarcinoma subtypes

Most cancer types exhibit aberrant transcriptional activity, including derepression of retrotransposable elements (RTEs). However, the degree, specificity and potential consequences of RTE transcriptional activation may differ substantially among cancer types and subtypes. Representing one extreme of the spectrum, we characterize the transcriptional activity of RTEs in cohorts of esophageal adenocarcinoma (EAC) and its precursor Barrett's esophagus (BE) from the OCCAMS (Oesophageal Cancer Clinical and Molecular Stratification) consortium, and from TCGA (The Cancer Genome Atlas). We found exceptionally high RTE inclusion in the EAC transcriptome, driven primarily by transcription of genes incorporating intronic or adjacent RTEs, rather than by autonomous RTE transcription. Nevertheless, numerous chimeric transcripts straddling RTEs and genes, and transcripts from stand-alone RTEs, particularly KLF5- and SOX9-controlled HERVH proviruses, were overexpressed specifically in EAC. Notably, incomplete mRNA splicing and EAC-characteristic intronic RTE inclusion was mirrored by relative loss of the respective fully-spliced, functional mRNA isoforms, consistent with compromised cellular fitness. Defective RNA splicing was linked with strong transcriptional activation of a HERVH provirus on Chr Xp22.32 and defined EAC subtypes with distinct molecular features and prognosis. Our study defines distinguishable RTE transcriptional profiles of EAC, reflecting distinct underlying processes and prognosis, thus providing a framework for targeted studies.

Transposable elements as essential elements in the control of gene expression

Interspersed repetitions called transposable elements (TEs), commonly referred to as mobile elements, make up a significant portion of the genomes of higher animals. TEs contribute in controlling the expression of genes locally and even far away at the transcriptional and post-transcriptional levels, which is one of their significant functional effects on gene function and genome evolution. There are different mechanisms through which TEs control the expression of genes. First, TEs offer cis-regulatory regions in the genome with their inherent regulatory features for their own expression, making them potential factors for controlling the expression of the host genes. Promoter and enhancer elements contain cis-regulatory sites generated from TE, which function as binding sites for a variety of trans-acting factors. Second, a significant portion of miRNAs and long non-coding RNAs (lncRNAs) have been shown to have TEs that encode for regulatory RNAs, revealing the TE origin of these RNAs. Furthermore, it was shown that TE sequences are essential for these RNAs' regulatory actions, which include binding to the target mRNA. By being a member of cis-regulatory and regulatory RNA sequences, TEs therefore play essential regulatory roles. Additionally, it has been suggested that TE-derived regulatory RNAs and cis-regulatory regions both contribute to the evolutionary novelty of gene regulation. Additionally, these regulatory systems arising from TE frequently have tissue-specific functions. The objective of this review is to discuss TE-mediated gene regulation, with a particular emphasis on the processes, contributions of various TE types, differential roles of various tissue types, based mostly on recent studies on humans.

Transposons are important contributors to gene expression variability under selection in rice populations

Transposable elements (TEs) are an important source of genome variability. Here, we analyze their contribution to gene expression variability in rice by performing a TE insertion polymorphism expression quantitative trait locus mapping using expression data from 208 varieties from the Oryza sativa ssp. indica and O. sativa ssp. japonica subspecies. Our data show that TE insertions are associated with changes of expression of many genes known to be targets of rice domestication and breeding. An important fraction of these insertions were already present in the rice wild ancestors, and have been differentially selected in indica and japonica rice populations. Taken together, our results show that small changes of expression in signal transduction genes induced by TE insertions accompany the domestication and adaptation of rice populations.

Locus specific endogenous retroviral expression associated with Alzheimer's disease

Introduction

Human endogenous retroviruses (HERVs) are transcriptionally-active remnants of ancient retroviral infections that may play a role in Alzheimer's disease.

Methods

We combined two, publicly available RNA-Seq datasets with a third, novel dataset for a total cohort of 103 patients with Alzheimer's disease and 45 healthy controls. We use telescope to perform HERV quantification for these samples and simultaneously perform gene expression analysis.

Results

We identify differentially expressed genes and differentially expressed HERVs in Alzheimer's disease patients. Differentially expressed HERVs are scattered throughout the genome; many of them are members of the HERV-K superfamily. A number of HERVs are correlated with the expression of dysregulated genes in Alzheimer's and are physically proximal to genes which drive disease pathways.

Discussion

Dysregulated expression of ancient retroviral insertions in the human genome are present in Alzheimer's disease and show localization patterns that may explain how these elements drive pathogenic gene expression.

Single-cell analysis reveals transcriptomic and epigenomic impacts on the maternal-fetal interface following SARS-CoV-2 infection

During pregnancy the maternal-fetal interface plays vital roles in fetal development. Its disruption is frequently found in pregnancy complications. Recent studies show increased incidences of adverse pregnancy outcomes in patients with COVID-19; however, the mechanism remains unclear. Here we analysed the molecular impacts of SARS-CoV-2 infection on the maternal-fetal interface. Generating bulk and single-nucleus transcriptomic and epigenomic profiles from patients with COVID-19 and control samples, we discovered aberrant immune activation and angiogenesis patterns in distinct cells from patients. Surprisingly, retrotransposons were also dysregulated in specific cell types. Notably, reduced enhancer activities of LTR8B elements were functionally linked to the downregulation of pregnancy-specific glycoprotein genes in syncytiotrophoblasts. Our findings revealed that SARS-CoV-2 infection induced substantial changes to the epigenome and transcriptome at the maternal-fetal interface, which may be associated with pregnancy complications.

KRAB-zinc-finger proteins regulate endogenous retroviruses to sculpt germline transcriptomes and genome evolution

As transposable elements (TEs) coevolved with the host genome, the host genome exploited TEs as functional regulatory elements. What remains largely unknown are how the activity of TEs, namely, endogenous retroviruses (ERVs), are regulated and how TEs evolved in the germline. Here we show that KRAB domain-containing zinc-finger proteins (KZFPs), which are highly expressed in mitotically dividing spermatogonia, bind to suppressed ERVs that function following entry into meiosis as active enhancers. These features are observed for independently evolved KZFPs and ERVs in mice and humans, i.e., are evolutionarily conserved in mammals. Further, we show that meiotic sex chromosome inactivation (MSCI) antagonizes the coevolution of KZFPs and ERVs in mammals. Our study uncovers a mechanism by which KZFPs regulate ERVs to sculpt germline transcriptomes. We propose that epigenetic programming in the mammalian germline during the mitosis-to-meiosis transition facilitates coevolution of KZFPs and TEs on autosomes and is antagonized by MSCI.

SINE Insertion in the Pig Carbonic Anhydrase 5B (CA5B) Gene Is Associated with Changes in Gene Expression and Phenotypic Variation

Transposons are genetic elements that are present in mammalian genomes and occupy a large proportion of the pig genome, with retrotransposons being the most abundant. In a previous study, it was found that a SINE retrotransposon was inserted in the 1^st intron of the CA5B gene in pigs, and the present study aimed to investigate the SINE insertion polymorphism in this gene in different pig breeds. Polymerase chain reaction (PCR) was used to confirm the polymorphism in 11 pig breeds and wild boars), and it was found that there was moderate polymorphism information content in 9 of the breeds. Further investigation in cell experiments revealed that the 330 bp SINE insertion in the RIP-CA5B site promoted expression activity in the weak promoter region of this site. Additionally, an enhancer verification vector experiment showed that the 330 bp SINE sequence acted as an enhancer on the core promoter region upstream of the CA5B gene region. The expression of CA5B in adipose tissue (back fat and leaf fat) in individuals with the (SINE^+/+) genotype was significantly higher than those with (SINE^+/-) and (SINE^-/-) genotypes. The association analysis revealed that the (SINE^+/+) genotype was significantly associated with a higher back fat thickness than the (SINE^-/-) genotype. Moreover, it was observed that the insertion of SINE at the RIP-CA5B site carried ATTT repeats, and three types of (ATTT) repeats were identified among different individuals/breeds (i.e., (ATTT)₄, (ATTT)₆ and (ATTT)₉). Overall, the study provides insights into the genetic basis of adipose tissue development in pigs and highlights the role of a SINE insertion in the CA5B gene in this process.

Long interspersed nuclear element 1 and B1/Alu repeats blueprint genome compartmentalization

The organization of the genome into euchromatin and heterochromatin has been known for almost 100 years [1]. More than 50% of mammalian genomes contain repetitive sequences [2,3]. Recently, a functional link between the genome and its folding has been identified [4,5]. Homotypic clustering of long interspersed nuclear element 1 (LINE1 or L1) and B1/Alu retrotransposons forms grossly exclusive nuclear domains that characterize and predict heterochromatin and euchromatin, respectively. The spatial segregation of L1 and B1/Alu-rich compartments is conserved in mammalian cells and can be rebuilt during the cell cycle and established de novo in early embryogenesis. Inhibition of L1 RNA drastically weakened homotypic repeat contacts and compartmental segregation, indicating that L1 plays a more significant role than just being a compartmental marker. This simple and inclusive genetic coding model of L1 and B1/Alu in shaping the macroscopic structure of the genome provides a plausible explanation for the remarkable conservation and robustness of its folding in mammalian cells. It also proposes a conserved core structure on which subsequent dynamic regulation takes place.

Gene editing tools for mycoplasmas: references and future directions for efficient genome manipulation

Mycoplasmas are successful pathogens that cause debilitating diseases in humans and various animal hosts. Despite the exceptionally streamlined genomes, mycoplasmas have evolved specific mechanisms to access essential nutrients from host cells. The paucity of genetic tools to manipulate mycoplasma genomes has impeded studies of the virulence factors of pathogenic species and mechanisms to access nutrients. This review summarizes several strategies for editing of mycoplasma genomes, including homologous recombination, transposons, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas system, and synthetic biology. In addition, the mechanisms and features of different tools are discussed to provide references and future directions for efficient manipulation of mycoplasma genomes.

The Immunological Conundrum of Endogenous Retroelements

Our defenses against infection rely on the ability of the immune system to distinguish invading pathogens from self. This task is exceptionally challenging, if not seemingly impossible, in the case of retroviruses that have integrated almost seamlessly into the host. This review examines the limits of innate and adaptive immune responses elicited by endogenous retroviruses and other retroelements, the targets of immune recognition, and the consequences for host health and disease. Contrary to theoretical expectation, endogenous retroelements retain substantial immunogenicity, which manifests most profoundly when their epigenetic repression is compromised, contributing to autoinflammatory and autoimmune disease and age-related inflammation. Nevertheless, recent evidence suggests that regulated immune reactivity to endogenous retroelements is integral to immune system development and function, underpinning cancer immunosurveillance, resistance to infection, and responses to the microbiota. Elucidation of the interaction points with endogenous retroelements will therefore deepen our understanding of immune system function and contribution to disease.

Evolution of piggyBac Transposons in Apoidea

In this study, we investigated the presence of piggyBac (PB) transposons in 44 bee genomes from the Apoidea order, which is a superfamily within the Hymenoptera, which includes a large number of bee species crucial for pollination. We annotated the PB transposons in these 44 bee genomes and examined their evolution profiles, including structural characteristics, distribution, diversity, activity, and abundance. The mined PB transposons were divided into three clades, with uneven distribution in each genus of PB transposons in Apoidea. The complete PB transposons we discovered are around 2.23-3.52 kb in length and encode transposases of approximately 580 aa, with terminal inverted repeats (TIRs) of about 14 bp and 4 bp (TTAA) target-site duplications. Long TIRs (200 bp, 201 bp, and 493 bp) were also detected in some species of bees. The DDD domains of the three transposon types were more conserved, while the other protein domains were less conserved. Generally, most PB transposons showed low abundance in the genomes of Apoidea. Divergent evolution dynamics of PB were observed in the genomes of Apoidea. PB transposons in some identified species were relatively young, whiles others were older and with some either active or inactive. In addition, multiple invasions of PB were also detected in some genomes of Apoidea. Our findings highlight the contribution of PB transposons to genomic variation in these species and suggest their potential as candidates for future gene transfer tools.

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.

RNA-mediated heterochromatin formation at repetitive elements in mammals

The failure to repress transcription of repetitive genomic elements can lead to catastrophic genome instability and is associated with various human diseases. As such, multiple parallel mechanisms cooperate to ensure repression and heterochromatinization of these elements, especially during germline development and early embryogenesis. A vital question in the field is how specificity in establishing heterochromatin at repetitive elements is achieved. Apart from trans-acting protein factors, recent evidence points to a role of different RNA species in targeting repressive histone marks and DNA methylation to these sites in mammals. Here, we review recent discoveries on this topic and predominantly focus on the role of RNA methylation, piRNAs, and other localized satellite RNAs.

Recent transposable element bursts are associated with the proximity to genes in a fungal plant pathogen

The activity of transposable elements (TEs) contributes significantly to pathogen genome evolution. TEs often destabilize genome integrity but may also confer adaptive variation in pathogenicity or resistance traits. De-repression of epigenetically silenced TEs often initiates bursts of transposition activity that may be counteracted by purifying selection and genome defenses. However, how these forces interact to determine the expansion routes of TEs within a pathogen species remains largely unknown. Here, we analyzed a set of 19 telomere-to-telomere genomes of the fungal wheat pathogen Zymoseptoria tritici. Phylogenetic reconstruction and ancestral state estimates of individual TE families revealed that TEs have undergone distinct activation and repression periods resulting in highly uneven copy numbers between genomes of the same species. Most TEs are clustered in gene poor niches, indicating strong purifying selection against insertions near coding sequences, or as a consequence of insertion site preferences. TE families with high copy numbers have low sequence divergence and strong signatures of defense mechanisms (i.e., RIP). In contrast, small non-autonomous TEs (i.e., MITEs) are less impacted by defense mechanisms and are often located in close proximity to genes. Individual TE families have experienced multiple distinct burst events that generated many nearly identical copies. We found that a Copia element burst was initiated from recent copies inserted substantially closer to genes compared to older copies. Overall, TE bursts tended to initiate from copies in GC-rich niches that escaped inactivation by genomic defenses. Our work shows how specific genomic environments features provide triggers for TE proliferation in pathogen genomes.

Contribution of Retrotransposons to the Pathogenesis of Type 1 Diabetes and Challenges in Analysis Methods

Type 1 diabetes (T1D) is one of the most common chronic diseases of the endocrine system, associated with several life-threatening comorbidities. While the etiopathogenesis of T1D remains elusive, a combination of genetic susceptibility and environmental factors, such as microbial infections, are thought to be involved in the development of the disease. The prime model for studying the genetic component of T1D predisposition encompasses polymorphisms within the HLA (human leukocyte antigen) region responsible for the specificity of antigen presentation to lymphocytes. Apart from polymorphisms, genomic reorganization caused by repeat elements and endogenous viral elements (EVEs) might be involved in T1D predisposition. Such elements are human endogenous retroviruses (HERVs) and non-long terminal repeat (non-LTR) retrotransposons, including long and short interspersed nuclear elements (LINEs and SINEs). In line with their parasitic origin and selfish behaviour, retrotransposon-imposed gene regulation is a major source of genetic variation and instability in the human genome, and may represent the missing link between genetic susceptibility and environmental factors long thought to contribute to T1D onset. Autoreactive immune cell subtypes with differentially expressed retrotransposons can be identified with single-cell transcriptomics, and personalized assembled genomes can be constructed, which can then serve as a reference for predicting retrotransposon integration/restriction sites. Here we review what is known to date about retrotransposons, we discuss the involvement of viruses and retrotransposons in T1D predisposition, and finally we consider challenges in retrotransposons analysis methods.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Characterisation of retrotransposon insertion polymorphisms in whole genome sequencing data from individuals with amyotrophic lateral sclerosis

The genetics of an individual is a crucial factor in understanding the risk of developing the neurodegenerative disease amyotrophic lateral sclerosis (ALS). There is still a large proportion of the heritability of ALS, particularly in sporadic cases, to be understood. Among others, active transposable elements drive inter-individual variability, and in humans long interspersed element 1 (LINE1, L1), Alu and SINE-VNTR-Alu (SVA) retrotransposons are a source of polymorphic insertions in the population. We undertook a pilot study to characterise the landscape of non-reference retrotransposon insertion polymorphisms (non-ref RIPs) in 15 control and 15 ALS individuals' whole genomes from Project MinE, an international project to identify potential genetic causes of ALS. The combination of two bioinformatics tools (mobile element locator tool (MELT) and TEBreak) identified on average 1250 Alu, 232 L1 and 77 SVA non-ref RIPs per genome across the 30 analysed. Further PCR validation of individual polymorphic retrotransposon insertions showed a similar level of accuracy for MELT and TEBreak. Our preliminary study did not identify a specific RIP or a significant difference in the total number of non-ref RIPs in ALS compared to control genomes. The use of multiple bioinformatic tools improved the accuracy of non-ref RIP detection and our study highlights the potential importance of studying these elements further in ALS.

The single-cell expression profile of transposable elements and transcription factors in human early biparental and uniparental embryonic development

Transposable elements (TEs) and transcription factors (TFs) are involved in the precise regulation of gene expression during the preimplantation stage. Activation of TEs is a key event for mammalian embryonic genome activation and preimplantation early embryonic development. TFs are involved in the regulation of drastic changes in gene expression patterns, but an inventory of the interplay between TEs and TFs during normal/abnormal human embryonic development is still lacking. Here we used single-cell RNA sequencing data generated from biparental and uniparental embryos to perform an integrative analysis of TE and TF expression. Our results showed that endogenous retroviruses (ERVs) are mainly expressed during the minor embryonic genome activation (EGA) process of early embryos, while Alu is gradually expressed in the middle and later stages. Some important ERVs (e.g., LTR5_Hs, MLT2A1) and Alu TEs are expressed at significantly lower levels in androgenic embryos. Integrative analysis revealed that the expression of the transcription factors CTCF and POU5F1 is correlated with the differential expression of ERV TEs. Comparative coexpression network analysis further showed distinct expression levels of important TFs (e.g., LEUTX and ZSCAN5A) in dizygotic embryos vs. parthenogenetic and androgenic embryos. This systematic investigation of TE and TF expression in human early embryonic development by single-cell RNA sequencing provides valuable insights into mammalian embryonic development.

Genome-wide chromatin accessibility analysis unveils open chromatin convergent evolution during polyploidization in cotton

Allopolyploidization, resulting in divergent genomes in the same cell, is believed to trigger a “genome shock”, leading to broad genetic and epigenetic changes. However, little is understood about chromatin and gene-expression dynamics as underlying driving forces during allopolyploidization. Here, we examined the genome-wide DNase I-hypersensitive site (DHS) and its variations in domesticated allotetraploid cotton ( Gossypium hirsutum and Gossypium barbadense , AADD) and its extant AA ( Gossypium arboreum ) and DD ( Gossypium raimondii ) progenitors. We observed distinct DHS distributions between G. arboreum and G. raimondii . In contrast, the DHSs of the two subgenomes of G. hirsutum and G. barbadense showed a convergent distribution. This convergent distribution of DHS was also present in the wild allotetraploids Gossypium darwinii and G. hirsutum var. yucatanense , but absent from a resynthesized hybrid of G. arboreum and G. raimondii , suggesting that it may be a common feature in polyploids, and not a consequence of domestication after polyploidization. We revealed that putative cis -regulatory elements (CREs) derived from polyploidization-related DHSs were dominated by several families, including Dof, ERF48, and BPC1. Strikingly, 56.6% of polyploidization-related DHSs were derived from transposable elements (TEs). Moreover, we observed positive correlations between DHS accessibility and the histone marks H3K4me3, H3K27me3, H3K36me3, H3K27ac, and H3K9ac, indicating that coordinated interplay among histone modifications, TEs, and CREs drives the DHS landscape dynamics under polyploidization. Collectively, these findings advance our understanding of the regulatory architecture in plants and underscore the complexity of regulome evolution during polyploidization.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae

The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.

Transcriptional and reverse transcriptional regulation of host genes by human endogenous retroviruses in cancers

Human endogenous retroviruses (HERVs) originated from ancient retroviral infections of germline cells millions of years ago and have evolved as part of the host genome. HERVs not only retain the capacity as retroelements but also regulate host genes. The expansion of HERVs involves transcription by RNA polymerase II, reverse transcription, and re-integration into the host genome. Fast progress in deep sequencing and functional analysis has revealed the importance of domesticated copies of HERVs, including their regulatory sequences, transcripts, and proteins in normal cells. However, evidence also suggests the involvement of HERVs in the development and progression of many types of cancer. Here we summarize the current state of knowledge about the expression of HERVs, transcriptional regulation of host genes by HERVs, and the functions of HERVs in reverse transcription and gene editing with their reverse transcriptase.