Endogenous retroviral LTRs as promoters for human genes: a critical assessment

Construction of the porcine genome mobile element variations and investigation of its role in population diversity and gene expression

BACKGROUND

Mobile element variants (MEVs) have a significant and complex impact on genomic diversity and phenotypic traits. However, the quantity, distribution, and relationship with gene expression and complex traits of MEVs in the pig genome remain poorly understood.

RESULTS

We constructed the most comprehensive porcine MEV library based on high-depth whole genome sequencing (WGS) data from 747 pigs across 59 breeds worldwide. This database identified a total of 147,993 polymorphic MEVs, including 121,099 short interspersed nuclear elements (SINEs), 26,053 long interspersed nuclear elements (LINEs), 802 long terminal repeats (LTRs), and 39 other transposons, among which 54% are newly discovered. We found that MEVs are unevenly distributed across the genome and are strongly influenced by negative selection effects. Importantly, we identified 514, 530, and 584 candidate MEVs associated with population differentiation, domestication, and breed formation, respectively. For example, a significantly differentiated MEV is located in the ATRX intron between Asian and European pigs, whereas ATRX is also differentially expressed between Asian and European pigs in muscle tissue. In addition, we identified 4,169 expressed MEVs (eMEVs) significantly associated with gene expression and 6,914 splicing MEVs (sMEVs) associated with gene splicing based on RNA-seq data from 266 porcine liver tissues. These eMEVs and sMEVs explain 6.24% and 9.47%, respectively, of the observed cis-heritability and highlight the important role of MEVs in the regulation of gene expression. Finally, we provide a high-quality SNP-MEV reference haplotype panel to impute MEV genotypes from genome-wide SNPs. Notably, we identified a candidate MEV significantly associated with total teat number, demonstrating the functionality of this reference panel.

CONCLUSIONS

The present investigation demonstrated the importance of MEVs in pigs in terms of population diversity, gene expression and phenotypic traits, which may provide useful resources and theoretical support for pig genetics and breeding.

Extensive binding of uncharacterized human transcription factors to genomic dark matter

Most of the human genome is thought to be non-functional, and includes large segments often referred to as "dark matter" DNA. The genome also encodes hundreds of putative and poorly characterized transcription factors (TFs). We determined genomic binding locations of 166 uncharacterized human TFs in living cells. Nearly half of them associated strongly with known regulatory regions such as promoters and enhancers, often at conserved motif matches and co-localizing with each other. Surprisingly, the other half often associated with genomic dark matter, at largely unique sites, via intrinsic sequence recognition. Dozens of these, which we term "Dark TFs", mainly bind within regions of closed chromatin. Dark TF binding sites are enriched for transposable elements, and are rarely under purifying selection. Some Dark TFs are KZNFs, which contain the repressive KRAB domain, but many are not: the Dark TFs also include known or potential pioneer TFs. Compiled literature information supports that the Dark TFs exert diverse functions ranging from early development to tumor suppression. Thus, our results sheds light on a large fraction of previously uncharacterized human TFs and their unappreciated activities within the dark matter genome.

Shared and unique patterns of autonomous human endogenous retrovirus loci transcriptomes in CD14 + monocytes from individuals with physical trauma or infection with COVID-19

Since previous studies have suggested that the RNAs of human endogenous retrovirus (HERV) might be involved in regulating innate immunity, it is important to investigate the HERV transcriptome patterns in innate immune cell types such as CD14 + monocytes. Using single cell RNA-seq datasets from resting or stimulated PBMCs mapped to 3,220 known discrete autonomous proviral HERV loci, we found individual-specific variation in HERV transcriptomes between HERV loci in CD14 + monocytes. Analysis of paired datasets from the same individual that were cultured in vitro with LPS or without (i.e. control) revealed 36 HERV loci in CD14 + monocytes that were detected only after activation. To extend our analysis to in vivo activated CD14 + monocytes, we used two scRNA-seq datasets from studies that had demonstrated activation of circulating CD14 + monocytes in patients with physical trauma or patients hospitalized with COVID-19 infections. For direct comparison between the trauma and COVID-19 datasets, we first analyzed 1.625 billion sequence reads from a composite pangenome control of 21 normal individuals. Comparison of the sequence read depth of HERV loci in the trauma or COVID-19 samples to the pangenome control revealed that 39 loci in the COVID-19 and 11 HERV loci in the trauma samples were significantly different (Mann-Whitney U test), with 9 HERV loci shared between the COVID-19 and trauma datasets. The capacity to compare HERV loci transcriptome patterns in innate immune cells, like CD14 + monocytes, across different pathological conditions will lead to greater understanding of the physiological role of HERV expression in health and disease.

Potential Contributions of Human Endogenous Retroviruses in Innate Immune Memory

The phenomenon wherein innate immune cells adopt long-term inflammatory phenotypes following the first stimuli is named trained immunity and can improve host defense against infections. Transcriptional and epigenetic reprogramming are critical mechanisms of trained immunity; however, the regulatory networks are not entirely clear at present. The human endogenous retroviruses (HERVs) provide large amounts of transcriptional regulators in the regulatory pathways. In this study, we analyzed published large omics data to explore the roles of such "dark matter" of the human genome in trained and tolerant macrophages. We collected 80 RNA sequencing data and 62 sequencing data to detect histone modifications and active regulatory regions from nine published studies on trained and tolerant macrophages. By analyzing the characteristics of transcription and epigenetic modification of HERVs, as well as their association with gene expression, we found that 15.3% of HERVs were transcribed nonrandomly from noncoding regions and enriched in specific HERV families and specific chromosomes, such as chromosomes 11, 15, 17, and 19, and they were highly related with the expression of adjacent genes. We found that 295 differentially expressed HERVs are located in 50-kbp flanking regions of 142 differentially expressed genes. We found epigenetic changes of these HERVs and that overlap with predicted enhancers and identified 35 enhancer-like HERVs. The related genes were highly involved in the activation and inflammatory responses, such as the TLR pathway. Other pathways including phosphoinositide signaling and transport of folate and K+ might be also related with trained immunity, which require further study. These results demonstrated that HERVs might play important roles in trained immunity.

Functional Bidirectionality of ERV-Derived Long Non-Coding RNAs in Humans

Human endogenous retroviruses (HERVs) are widely recognized as the result of exogenous retroviruses infecting the ancestral germline, stabilizing integration and vertical transmission during human genetic evolution. To date, endogenous retroviruses (ERVs) appear to have been selected for human physiological functions with the loss of retrotransposable capabilities. ERV elements were previously regarded as junk DNA for a long time. Since then, the aberrant activation and expression of ERVs have been observed in the development of many kinds of human diseases, and their role has been explored in a variety of human disorders such as cancer. The results show that specific ERV elements play respective crucial roles. Among them, long non-coding RNAs (lncRNAs) transcribed from specific long-terminal repeat regions of ERVs are often key factors. lncRNAs are over 200 nucleotides in size and typically bind to DNA, RNA, and proteins to perform biological functions. Dysregulated lncRNAs have been implicated in a variety of diseases. In particular, studies have shown that the aberrant expression of some ERV-derived lncRNAs has a tumor-suppressive or oncogenic effect, displaying significant functional bidirectionality. Therefore, theses lncRNAs have a promising future as novel biomarkers and therapeutic targets to explore the concise relationship between ERVs and cancers. In this review, we first summarize the role of ERV-derived lncRNAs in physiological regulation, mainly including immunomodulation, the maintenance of pluripotency, and erythropoiesis. In addition, pathological regulation examples of their aberrant activation and expression leading to carcinogenesis are highlighted, and specific mechanisms of occurrence are discussed.

Novel function of U7 snRNA in the repression of HERV1/LTR12s and lincRNAs in human cells

U7 snRNA is part of the U7 snRNP complex, required for the 3' end processing of replication-dependent histone pre-mRNAs in S phase of the cell cycle. Here, we show that U7 snRNA plays another function in inhibiting the expression of a subset of long terminal repeats of human endogenous retroviruses (HERV1/LTR12s) and LTR12-containing long intergenic noncoding RNAs (lincRNAs), both bearing sequence motifs that perfectly match the 5' end of U7 snRNA. We demonstrate that U7 snRNA inhibits LTR12 and lincRNA transcription and propose a mechanism in which U7 snRNA hampers the binding/activity of the NF-Y transcription factor to CCAAT motifs within LTR12 elements. Thereby, U7 snRNA plays a protective role in maintaining the silencing of deleterious genetic elements in selected types of cells.

Epigenetic modifications of gonadotropin receptors can regulate follicular development

The spatiotemporal transcription of follicle-stimulating hormone receptor (FSHR) and luteinizing hormone/human chorionic gonadotropin receptor (LHCGR) are crucial events for follicular development. However, their regulatory mechanisms are unclear. DNA methylation and histone acetylation are the main epigenetic modifications, and play important roles in transcriptional expression, which regulate cell responses including cell proliferation, senescence and apoptosis. This review will discuss the dynamic epigenetic modifications of FSHR and LHCGR that occur during the process of follicular development and their response to gonadotropins. In addition, some alteration patterns that occur during these epigenetic modifications, as well as their retrospect retrotransposons, which regulate the gene expression levels of FSHR and LHCGR will be discussed.

Dissecting transposable elements and endogenous retroviruses upregulation by HDAC inhibitors in leiomyosarcoma cells: Implications for the interferon response

Transposable elements (TEs) are of interest as immunomodulators for cancer therapies. TEs can fold into dsRNAs that trigger the interferon response. Here, we investigated the effect of different HDAC inhibitors (HDACIs) on the expression of TEs in leiomyosarcoma cells. Our data show that endogenous retroviruses (ERVs), especially ERV1 elements, are upregulated after treatment with HDAC1/2/3-specific inhibitors. Surprisingly, the interferon response was not activated. We observed an increase in A-to-I editing of upregulated ERV1. This could have an impact on the stability of dsRNAs and the activation of the interferon response. We also found that H3K27ac levels are increased in the LTR12 subfamilies, which could be regulatory elements controlling the expression of proapoptotic genes such as TNFRSF10B. In summary, we provide a detailed characterization of TEs modulation in response to HDACIs and suggest the use of HDACIs in combination with ADAR inhibitors to induce cell death and support immunotherapy in cancer.

Transcription of Endogenous Retroviruses: Broad and Precise Mechanisms of Control

Endogenous retroviruses (ERVs) are the remnants of retroviral germline infections and are highly abundant in the genomes of vertebrates. At one time considered to be nothing more than inert 'junk' within genomes, ERVs have been tolerated within host genomes over vast timescales, and their study continues to reveal complex co-evolutionary histories within their respective host species. For example, multiple instances have been characterized of ERVs having been 'borrowed' for normal physiology, from single copies to ones involved in various regulatory networks such as innate immunity and during early development. Within the cell, the accessibility of ERVs is normally tightly controlled by epigenetic mechanisms such as DNA methylation or histone modifications. However, these silencing mechanisms of ERVs are reversible, and epigenetic alterations to the chromatin landscape can thus lead to their aberrant expression, as is observed in abnormal cellular environments such as in tumors. In this review, we focus on ERV transcriptional control and draw parallels and distinctions concerning the loss of regulation in disease, as well as their precise regulation in early development.

The regulatory potential of transposable elements in maize

Since their initial discovery in maize, transposable elements (TEs) have emerged as being integral to the evolution of maize, accounting for 80% of its genome. However, the repetitive nature of TEs has hindered our understanding of their regulatory potential. Here, we demonstrate that long-read chromatin fiber sequencing (Fiber-seq) permits the comprehensive annotation of the regulatory potential of maize TEs. We uncover that only 94 LTR retrotransposons contain the functional epigenetic architecture required for mobilization within maize leaves. This epigenetic architecture degenerates with evolutionary age, resulting in solo TE enhancers being preferentially marked by simultaneous hyper-CpG methylation and chromatin accessibility, an architecture markedly divergent from canonical enhancers. We find that TEs shape maize gene regulation by creating novel promoters within the TE itself as well as through TE-mediated gene amplification. Lastly, we uncover a pervasive epigenetic code directing TEs to specific loci, including that locus that sparked McClintock's discovery of TEs.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Endogenous retroviral solo-LTRs in human genome

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology.

Nanopore-based DNA long-read sequencing analysis of the aged human brain

Aging disrupts cellular processes such as DNA repair and epigenetic control, leading to a gradual buildup of genomic alterations that can have detrimental effects in post-mitotic cells. Genomic alterations in regions of the genome that are rich in repetitive sequences, often termed "dark loci," are difficult to resolve using traditional sequencing approaches. New long-read technologies offer promising avenues for exploration of previously inaccessible regions of the genome. Using nanopore-based long-read whole-genome sequencing of DNA extracted from aged 18 human brains, we identify previously unreported structural variants and methylation patterns within repetitive DNA, focusing on transposable elements ("jumping genes") as crucial sources of variation, particularly in dark loci. Our analyses reveal potential somatic insertion variants and provides DNA methylation frequencies for many retrotransposon families. We further demonstrate the utility of this technology for the study of these challenging genomic regions in brains affected by Alzheimer's disease and identify significant differences in DNA methylation in pathologically normal brains versus those affected by Alzheimer's disease. Highlighting the power of this approach, we discover specific polymorphic retrotransposons with altered DNA methylation patterns. These retrotransposon loci have the potential to contribute to pathology, warranting further investigation in Alzheimer's disease research. Taken together, our study provides the first long-read DNA sequencing-based analysis of retrotransposon sequences, structural variants, and DNA methylation in the aging brain affected with Alzheimer's disease neuropathology.

Possibilities of skin coat color-dependent risks and risk factors of squamous cell carcinoma and deafness of domestic cats inferred via RNA-seq data

The transcriptome data of skin cells from domestic cats with brown, orange, and white coats were analyzed using a public database to investigate the possible relationship between coat color-related gene expression and squamous cell carcinoma risk, as well as the mechanism of deafness in white cats. We found that the ratio of the expression level of genes suppressing squamous cell carcinoma to that of genes promoting squamous cell carcinoma might be considerably lower than the theoretical estimation in skin cells with orange and white coats in white-spotted cat. We also found the possibility of the frequent production of KIT lacking the first exon (d1KIT) in skin cells with white coats, and d1KIT production exhibited a substantial negative correlation with the expression of SOX10, which is essential for melanocyte formation and adjustment of hearing function. Additionally, the production of d1KIT was expected to be due to the insulating activity of the feline endogenous retrovirus 1 (FERV1) LTR in the first intron of KIT by its CTCF binding sequence repeat. These results contribute to basic veterinary research to understand the relationship between cat skin coat and disease risk, as well as the underlying mechanism.

Expression and high levels of insertional polymorphism of an endogenous gammaretrovirus lineage in dogs

Despite the absence of a confirmed exogenously replicating retrovirus in Canis lupus familiaris (C. familiaris), past retroviral infections are evident in the genomes of living animals via the presence of endogenous retroviruses (ERVs). Although gammaretrovirus-like transcripts and enzyme activities were previously reported to be present in canine leukemias and lymphomas, those findings were not further explored. Initial analysis of the C. familiaris reference genome revealed a minor subset of one ERV lineage, classified as CfERV-Fc1(a), or Fc1(a) here, with features characteristic of recent integration, including the presence of ORFs and identical or nearly identical LTRs. Our previous analysis of whole genome sequence data belonging to extant Canidae revealed a burst of past infections in Canis ancestors resulting in numerous young, polymorphic, and highly intact loci now segregating in dogs. Here, we demonstrate the expression of full-length Fc1(a) proviruses in tissues collected from healthy animals and from animals with cancer. We observed significantly higher expression in samples of dogs with various cancer diagnoses when compared to samples from healthy dogs. Genotyping of insertionally polymorphic Fc1(a) loci identified candidate expressed proviruses and delineated distributions over sample groups. Collectively, the data show that Fc1(a) proviruses retain biological activity in the domestic dog and provides a means to examine potential genetic links with disease states in this species.

Human Endogenous Retrovirus-K (HML-2)-Related Genetic Variation: Human Genome Diversity and Disease

Human endogenous retroviruses (HERVs) comprise a significant portion of the human genome, making up roughly 8%, a notable comparison to the 2-3% represented by coding sequences. Numerous studies have underscored the critical role and importance of HERVs, highlighting their diverse and extensive influence on the evolution of the human genome and establishing their complex correlation with various diseases. Among HERVs, the HERV-K (HML-2) subfamily has recently attracted significant attention, integrating into the human genome after the divergence between humans and chimpanzees. Its insertion in the human genome has received considerable attention due to its structural and functional characteristics and the time of insertion. Originating from ancient exogenous retroviruses, these elements succeeded in infecting germ cells, enabling vertical transmission and existing as proviruses within the genome. Remarkably, these sequences have retained the capacity to form complete viral sequences, exhibiting activity in transcription and translation. The HERV-K (HML-2) subfamily is the subject of active debate about its potential positive or negative effects on human genome evolution and various pathologies. This review summarizes the variation, regulation, and diseases in human genome evolution arising from the influence of HERV-K (HML-2).

Progressive Exaptation of Endogenous Retroviruses in Placental Evolution in Cattle

Viviparity is made possible by the placenta, a structure acquired relatively recently in the evolutionary history of eutherian mammals. Compared to oviparity, it increases the survival rate of the fetus, owing to the eutherian placenta. Questions such as "How was the placenta acquired?" and "Why is there diversity in placental morphology among mammalian species?" remain largely unsolved. Our present understanding of the molecules regulating placental development remains unclear, owing in no small part to the persistent obscurity surrounding the molecular mechanisms underlying placental acquisition. Numerous genes associated with the development of eutherian placental morphology likely evolved to function at the fetal-maternal interface in conjunction with those participating in embryogenesis. Therefore, identifying these genes, how they were acquired, and how they came to be expressed specifically at the fetal-maternal interface will shed light on some crucial molecular mechanisms underlying placental evolution. Exhaustive studies support the hypothesis that endogenous retroviruses (ERVs) could be evolutional driving forces for trophoblast cell fusion and placental structure in mammalian placentas including those of the bovine species. This review focuses on bovine ERVs (BERVs) and their expression and function in the placenta.

Unveiling the hidden players: The crucial role of transposable elements in the placenta and their potential contribution to pre-eclampsia

The human placenta is a vital connection between maternal and fetal tissues, allowing for the exchange of molecules and modulation of immune interactions during pregnancy. Interestingly, some of the placenta's unique functions can be attributed to transposable elements (TEs), which are DNA sequences that have mobilised into the genome. Co-option throughout mammalian evolution has led to the generation of TE-derived regulators and TE-derived genes, some of which are expressed in the placenta but silenced in somatic tissues. TE genes encompass both TE-derived genes with a repeat element in the coding region and TE-derived regulatory regions such as alternative promoters and enhancers. Placental-specific TE genes are known to contribute to the placenta's unique functions, and interestingly, they are also expressed in some cancers and share similar functions. There is evidence to support that aberrant activity of TE genes may contribute to placental pathologies, cancer and autoimmunity. In this review, we highlight the crucial roles of TE genes in placental function, and how their dysregulation may lead to pre-eclampsia, a common and dangerous placental condition. We provide a summary of the functional TE genes in the placenta to offer insight into their significance in normal and abnormal human development. Ultimately, this review highlights an opportunity for future research to investigate the potential dysregulation of TE genes in the development of placental pathologies such as pre-eclampsia. Further understanding of TE genes and their role in the placenta could lead to significant improvements in maternal and fetal health.

Genome-wide identification of endogenous retrovirus elements and their active transcription in mink genome

Mammalian endogenous retroviruses (ERVs) are ancient retroviruses that have been integrated into genomes. ERVs were believed to be inactive until the discovery of ERV transcription in the mouse genome. However, the transcription level and function of ERV elements in mammalian genomes are not well understood. In this study, we performed the first genome-wide scanning of ERV loci in the American mink (Neogale vison) genome (NeoERV) followed by transcriptomic analysis to detect actively transcribed NeoERV elements. A total of 365,791 NeoERV loci were identified, and161,205 (44%) of these loci were found to be actively transcribed based on transcriptomic data from three types of tissues (amygdala, trachea and lung). More than one third of the actively transcribed NeoERV loci were tissue-specific. Furthermore, some of the active loci were associated with host gene transcription, and the level of NeoERV transcription was positively correlated with that of host genes, specifically when active loci were located in overlapped gene regions. An in-depth analysis of the envelope protein coding env gene showed that, in general, its transcription level was higher than that of NeoERVs, which is believed to be associated with host immunity.

An old model with new insights: endogenous retroviruses drive the evolvement toward ASD susceptibility and hijack transcription machinery during development

The BTBR T+Itpr3tf/J (BTBR/J) strain is one of the most valid models of idiopathic autism, serving as a potent forward genetics tool to dissect the complexity of autism. We found that a sister strain with an intact corpus callosum, BTBR TF/ArtRbrc (BTBR/R), showed more prominent autism core symptoms but moderate ultrasonic communication/normal hippocampus-dependent memory, which may mimic autism in the high functioning spectrum. Intriguingly, disturbed epigenetic silencing mechanism leads to hyperactive endogenous retrovirus (ERV), a mobile genetic element of ancient retroviral infection, which increases de novo copy number variation (CNV) formation in the two BTBR strains. This feature makes the BTBR strain a still evolving multiple-loci model toward higher ASD susceptibility. Furthermore, active ERV, analogous to virus infection, evades the integrated stress response (ISR) of host defense and hijacks the transcriptional machinery during embryonic development in the BTBR strains. These results suggest dual roles of ERV in the pathogenesis of ASD, driving host genome evolution at a long-term scale and managing cellular pathways in response to viral infection, which has immediate effects on embryonic development. The wild-type Draxin expression in BTBR/R also makes this substrain a more precise model to investigate the core etiology of autism without the interference of impaired forebrain bundles as in BTBR/J.

Regulation of human trophoblast gene expression by endogenous retroviruses

The placenta is a fast-evolving organ with large morphological and histological differences across eutherians, but the genetic changes driving placental evolution have not been fully elucidated. Transposable elements, through their capacity to quickly generate genetic variation and affect host gene regulation, may have helped to define species-specific trophoblast gene expression programs. Here we assess the contribution of transposable elements to human trophoblast gene expression as enhancers or promoters. Using epigenomic data from primary human trophoblast and trophoblast stem-cell lines, we identified multiple endogenous retrovirus families with regulatory potential that lie close to genes with preferential expression in trophoblast. These largely primate-specific elements are associated with inter-species gene expression differences and are bound by transcription factors with key roles in placental development. Using genetic editing, we demonstrate that several elements act as transcriptional enhancers of important placental genes, such as CSF1R and PSG5. We also identify an LTR10A element that regulates ENG expression, affecting secretion of soluble endoglin, with potential implications for preeclampsia. Our data show that transposons have made important contributions to human trophoblast gene regulation, and suggest that their activity may affect pregnancy outcomes.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone H4 or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

HERV-K (HML-2) insertion polymorphisms in the 8q24.13 region and their potential etiological associations with acute myeloid leukemia

Human endogenous retroviruses (HERVs) are LTR retrotransposons that are present in the human genome. Among them, members of the HERV-K (HML-2) group are suspected to play a role in the development of different types of cancer, including lung, ovarian, and prostate cancer, as well as leukemia. Acute myeloid leukemia (AML) is an important disease that causes 1% of cancer deaths in the United States and has a survival rate of 28.7%. Here, we describe a method for assessing the statistical association between HERV-K (HML-2) transposable element insertion polymorphisms (or TIPs) and AML, using whole-genome sequencing and read mapping using TIP_finder software. Our results suggest that 101 polymorphisms involving HERV-K (HML-2) elements were correlated with AML, with a percentage between 44.4 to 56.6%, most of which (70) were located in the region from 8q24.13 to 8q24.21. Moreover, it was found that the TRIB1, LRATD2, POU5F1B, MYC, PCAT1, PVT1, and CCDC26 genes could be displaced or fragmented by TIPs. Furthermore, a general method was devised to facilitate analysis of the correlation between transposable element insertions and specific diseases. Finally, although the relationship between HERV-K (HML-2) TIPs and AML remains unclear, the data reported in this study indicate a statistical correlation, as supported by the χ² test with p-values < 0.05.

Evolvability and emergence of tumor heterogeneity as a space-time function

The loss of control of cell proliferation, apoptosis regulation and contact inhibition leads to tumor development. While benign tumors are restricted to their primary space, i.e. where these tumors first originate, the metastatic tumors not only disseminate- facilitated by hypoxia-driven neovascularization- to distant secondary sites but also show substantial changes in metabolism, tissue architectures, gene expression profiles and immune phenotypes. All these alterations result in radio-, chemo- and immune-resistance rendering these metastatic tumor cells refractory to therapy. Since the beginning of the transformation, these factors- which influence each other- are incorporated to the developing and metastasizing tumor. As a result, the complexities in the heterogeneity of tumor progressively increase. This space-time function in the heterogeneity of tumors is generated by various conditions and factors at the genetic as well as microenvironmental levels, for example, endogenous retroviruses, methylation and epigenetic dysregulation that may be etiology-specific, cancer associated inflammation, remodeling of the extracellular matrix and mesenchymal cell shifted functions. On the one hand, these factors may cause de-differentiation of the tumor cells leading to cancer stem cells that contribute to radio-, chemo- and immune-resistance and recurrence of tumors. On the other hand, they may also enhance the heterogeneity under specific microenvironment-driven proliferation. In this editorial, we intend to underline the importance of heterogeneity in cancer progress, its evaluation and its use in correlation with the tumor evolution in a specific patient as a field of research for achieving precise patient-tailored treatments and amelioration of diagnostic (monitoring) tools and prognostic capacity.

T3E: a tool for characterising the epigenetic profile of transposable elements using ChIP-seq data

BACKGROUND

Despite the advent of Chromatin Immunoprecipitation Sequencing (ChIP-seq) having revolutionised our understanding of the mammalian genome's regulatory landscape, many challenges remain. In particular, because of their repetitive nature, the sequencing reads derived from transposable elements (TEs) pose a real bioinformatics challenge, to the point that standard analysis pipelines typically ignore reads whose genomic origin cannot be unambiguously ascertained.

RESULTS

We show that discarding ambiguously mapping reads may lead to a systematic underestimation of the number of reads associated with young TE families/subfamilies. We also provide evidence suggesting that the strategy of randomly permuting the location of the read mappings (or the TEs) that is often used to compute the background for enrichment calculations at TE families/subfamilies can result in both false positive and negative enrichments. To address these problems, we present the Transposable Element Enrichment Estimator (T3E), a tool that makes use of ChIP-seq data to characterise the epigenetic profile of associated TE families/subfamilies. T3E weights the number of read mappings assigned to the individual TE copies of a family/subfamily by the overall number of genomic loci to which the corresponding reads map, and this is done at the single nucleotide level. In addition, T3E computes ChIP-seq enrichment relative to a background estimated based on the distribution of the read mappings in the input control DNA. We demonstrated the capabilities of T3E on 23 different ChIP-seq libraries. T3E identified enrichments that were consistent with previous studies. Furthermore, T3E detected context-specific enrichments that are likely to pinpoint unexplored TE families/subfamilies with individual TE copies that have been frequently exapted as cis-regulatory elements during the evolution of mammalian regulatory networks.

CONCLUSIONS

T3E is a novel open-source computational tool (available for use at: https://github.com/michelleapaz/T3E ) that overcomes some of the pitfalls associated with the analysis of ChIP-seq data arising from the repetitive mammalian genome and provides a framework to shed light on the epigenetics of entire TE families/subfamilies.

Primate-specific transposable elements shape transcriptional networks during human development

The human genome contains more than 4.5 million inserts derived from transposable elements (TEs), the result of recurrent waves of invasion and internal propagation throughout evolution. For new TE copies to be inherited, they must become integrated in the genome of the germline or pre-implantation embryo, which requires that their source TE be expressed at these stages. Accordingly, many TEs harbor DNA binding sites for the pluripotency factors OCT4, NANOG, SOX2, and KLFs and are transiently expressed during embryonic genome activation. Here, we describe how many primate-restricted TEs have additional binding sites for lineage-specific transcription factors driving their expression during human gastrulation and later steps of fetal development. These TE integrants serve as lineage-specific enhancers fostering the transcription, amongst other targets, of KRAB-zinc finger proteins (KZFPs) of comparable evolutionary age, which in turn corral the activity of TE-embedded regulatory sequences in a similarly lineage-restricted fashion. Thus, TEs and their KZFP controllers play broad roles in shaping transcriptional networks during early human development.

Alternative and aberrant splicing of human endogenous retroviruses in cancer. What about head and neck? —A mini review

Human endogenous retroviruses (HERVs) are transcribed in many cancer types, including head and neck cancer. Because of accumulating mutations at proviral loci over evolutionary time, HERVs are functionally defective and cannot complete their viral life cycle. Despite that, HERV transcripts, including full-length viral RNAs and viral RNAs spliced as expected at the conventional viral splice sites, can be detected in particular conditions, such as cancer. Interestingly, non-viral–related transcription, including aberrant, non-conventionally spliced RNAs, has been reported as well. The role of HERV transcription in cancer and its contribution to oncogenesis or progression are still debated. Nonetheless, HERVs may constitute a suitable cancer biomarker or a target for therapy. Thus, ongoing research aims both to clarify the basic mechanisms underlying HERV transcription in cancer and to exploit its potential toward clinical application. In this mini-review, we summarize the current knowledge, the most recent findings, and the future perspectives of research on HERV transcription and splicing, with particular focus on head and neck cancer.

Comparative Methylome Analysis Reveals Epigenetic Signatures Associated with Growth and Shell Color in the Pacific Oyster, Crassostrea gigas

Fast growth is one of the most important breeding goals for all economic species such as the Pacific oyster (Crassostrea gigas), an aquaculture mollusk with top global production. Although the genetic basis and molecular mechanisms of growth-related traits have been widely investigated in the oyster, the role of DNA methylation involved in growth regulation remains largely unexplored. In this study, we performed a comparative DNA methylome analysis of two selectively bred C. gigas strains with contrasted phenotypes in growth and shell color based on whole-genome bisulfite sequencing (WGBS). Genome-wide profiling of DNA methylation at the single-base resolution revealed that DNA methylations were widely spread across the genome with obvious hotspots, coinciding with the distribution of genes and repetitive elements. Higher methylation levels were observed within genic regions compared with intergenic and promoter regions. Comparative analysis of DNA methylation allowed the identification of 339,604 differentially methylated CpG sites (DMCs) clustering in 27,600 differentially methylated regions (DMRs). Functional annotation analysis identified 11,033 genes from DMRs which were enriched in biological processes including cytoskeleton system, cell cycle, signal transduction, and protein biosynthesis. Integrative analysis of methylome and transcriptome profiles revealed a positive correlation between gene expression and DNA methylation within gene-body regions. Protein-protein interaction (PPI) analysis of differentially expressed and methylated genes allowed for the identification of integrin beta-6 (homolog of human ITGB3) as a hub modulator of the PI3K/Akt signaling pathway that was involved in various growth-related processes. This work provided insights into epigenetic regulation of growth in oysters and will be valuable resources for studying DNA methylation in invertebrates.

Endogenous Retroviral Sequences Behave as Putative Enhancers Controlling Gene Expression through HP1-Regulated Long-Range Chromatin Interactions

About half of the mammalian genome is constituted of repeated elements, among which endogenous retroviruses (ERVs) are known to influence gene expression and cancer development. The HP1 (Heterochromatin Protein 1) proteins are known to be essential for heterochromatin establishment and function and its loss in hepatocytes leads to the reactivation of specific ERVs and to liver tumorigenesis. Here, by studying two ERVs located upstream of genes upregulated upon loss of HP1, Mbd1 and Trim24, we show that these HP1-dependent ERVs behave as either alternative promoters or as putative enhancers forming a loop with promoters of endogenous genes depending on the genomic context and HP1 expression level. These ERVs are characterised by a specific HP1-independent enrichment in heterochromatin-associated marks H3K9me3 and H4K20me3 as well as in the enhancer-specific mark H3K4me1, a combination that might represent a bookmark of putative ERV-derived enhancers. These ERVs are further enriched in a HP1-dependent manner in H3K27me3, suggesting a critical role of this mark together with HP1 in the silencing of the ERVs, as well as for the repression of the associated genes. Altogether, these results lead to the identification of a new regulatory hub involving the HP1-dependent formation of a physical loop between specific ERVs and endogenous genes.

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.

Roles of transposable elements in the regulation of mammalian transcription

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

Scan of the endogenous retrovirus sequences across the swine genome and survey of their copy number variation and sequence diversity among various Chinese and Western pig breeds

In pig-to-human xenotransplantation, the transmission risk of porcine endogenous retroviruses (PERVs) is of great concern. However, the distribution of PERVs in pig genomes, their genetic variation among Eurasian pigs, and their evolutionary history remain unclear. We scanned PERVs in the current pig reference genome (assembly Build 11.1), and identified 36 long complete or near-complete PERVs (lcPERVs) and 23 short incomplete PERVs (siPERVs). Besides three known PERVs (PERV-A, -B, and -C), four novel types (PERV-JX1, -JX2, -JX3, and -JX4) were detected in this study. According to evolutionary analyses, the newly discovered PERVs were more ancient, and PERV-Bs probably experienced a bottleneck ~0.5 million years ago (Ma). By analyzing 63 high-quality porcine whole-genome resequencing data, we found that the PERV copy numbers in Chinese pigs were lower (32.0±4.0) than in Western pigs (49.1±6.5). Additionally, the PERV sequence diversity was lower in Chinese pigs than in Western pigs. Regarding the lcPERV copy numbers, PERV-A and -JX2 in Western pigs were higher than in Chinese pigs. Notably, Bama Xiang (BMX) pigs had the lowest PERV copy number (27.8±5.1), and a BMX individual had no PERV-C and the lowest PERV copy number (23), suggesting that BMX pigs were more suitable for screening and/or modification as xenograft donors. Furthermore, we identified 451 PERV transposon insertion polymorphisms (TIPs), of which 86 were shared by all 10 Chinese and Western pig breeds. Our findings provide systematic insights into the genomic distribution, variation, evolution, and possible biological function of PERVs.

Endogenous retroviruses transcriptomes in response to four avian pathogenic microorganisms infection in chicken

The impact of Endogenous retroviruses (ERVs) on chicken disease is not well understood. Here, we systematically identified 436 relatively complete ChERVs from the chicken genome. Subsequently, ChERV transcriptomes were analyzed in chicken after subgroup J avian leukosis virus (ALV-J), avian influenza virus (AIV), Marek's disease virus (MDV) and avian pathogenic Escherichia coli (APEC) infection. We found that about 50%-68% of ChERVs were transcriptionally active in infected and uninfected-samples, although the abundance of most ChERVs is relatively low. Moreover, compared to uninfected-samples, 49, 18, 66 and 17 ChERVs were significantly differentially expressed in ALV-J, AIV, MDV and APEC infected-samples, respectively. These findings may be of significance for understanding the role and function of ChERVs to response the pathogenic microorganism infection.

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

Positive Selection and Enhancer Evolution Shaped Lifespan and Body Mass in Great Apes

Within primates, the great apes are outliers both in terms of body size and lifespan, since they include the largest and longest-lived species in the order. Yet, the molecular bases underlying such features are poorly understood. Here, we leveraged an integrated approach to investigate multiple sources of molecular variation across primates, focusing on over 10,000 genes, including approximately 1,500 previously associated with lifespan, and additional approximately 9,000 for which an association with longevity has never been suggested. We analyzed dN/dS rates, positive selection, gene expression (RNA-seq), and gene regulation (ChIP-seq). By analyzing the correlation between dN/dS, maximum lifespan, and body mass, we identified 276 genes whose rate of evolution positively correlates with maximum lifespan in primates. Further, we identified five genes, important for tumor suppression, adaptive immunity, metastasis, and inflammation, under positive selection exclusively in the great ape lineage. RNA-seq data, generated from the liver of six species representing all the primate lineages, revealed that 8% of approximately 1,500 genes previously associated with longevity are differentially expressed in apes relative to other primates. Importantly, by integrating RNA-seq with ChIP-seq for H3K27ac (which marks active enhancers), we show that the differentially expressed longevity genes are significantly more likely than expected to be located near a novel "ape-specific" enhancer. Moreover, these particular ape-specific enhancers are enriched for young transposable elements, and specifically SINE-Vntr-Alus. In summary, we demonstrate that multiple evolutionary forces have contributed to the evolution of lifespan and body size in primates.

Mammalian tumor-like organs. 1. The role of tumor-like normal organs and atypical tumor organs in the evolution of development (carcino-evo-devo)

Background

Earlier I hypothesized that hereditary tumors might participate in the evolution of multicellular organisms. I formulated the hypothesis of evolution by tumor neofunctionalization, which suggested that the evolutionary role of hereditary tumors might consist in supplying evolving multicellular organisms with extra cell masses for the expression of evolutionarily novel genes and the origin of new cell types, tissues, and organs. A new theory—the carcino-evo-devo theory—has been developed based on this hypothesis.

Main text

My lab has confirmed several non-trivial predictions of this theory. Another non-trivial prediction is that evolutionarily new organs if they originated from hereditary tumors or tumor-like structures, should recapitulate some tumor features in their development. This paper reviews the tumor-like features of evolutionarily novel organs. It turns out that evolutionarily new organs such as the eutherian placenta, mammary gland, prostate, the infantile human brain, and hoods of goldfishes indeed have many features of tumors. I suggested calling normal organs, which have many tumor features, the tumor-like organs.

Conclusion

Tumor-like organs might originate from hereditary atypical tumor organs and represent the part of carcino-evo-devo relationships, i.e., coevolution of normal and neoplastic development. During subsequent evolution, tumor-like organs may lose the features of tumors and the high incidence of cancer and become normal organs without (or with almost no) tumor features.

Emerging roles for endogenous retroviruses in immune epigenetic regulation

In recent years, there has been significant progress toward understanding the transcriptional networks underlying mammalian immune responses, fueled by advances in regulatory genomic technologies. Epigenomic studies profiling immune cells have generated detailed genome-wide maps of regulatory elements that will be key to deciphering the regulatory networks underlying cellular immune responses and autoimmune disorders. Unbiased analyses of these genomic maps have uncovered endogenous retroviruses as an unexpected ally in the regulation of human immune systems. Despite their parasitic origins, studies are finding an increasing number of examples of retroviral sequences having been co-opted for beneficial immune function and regulation by the host cell. Here, we review how endogenous retroviruses have given rise to numerous regulatory elements that shape the epigenetic landscape of host immune responses. We will discuss the implications of these elements on the function, dysfunction, and evolution of innate immunity.

Long Terminal Repeat Retrotransposon Afut4 Promotes Azole Resistance of Aspergillus fumigatus by Enhancing the Expression of sac1 Gene

Aspergillus fumigatus causes a series of invasive diseases, including the high-mortality invasive aspergillosis, and has been a serious global health threat because of its increased resistance to the first-line clinical triazoles. We analyzed the whole-genome sequence of 15 A. fumigatus strains from China and found that long terminal repeat retrotransposons (LTR-RTs), including Afut1, Afut2, Afut3, and Afut4, are most common and have the largest total nucleotide length among all transposable elements in A. fumigatus. Deleting one of the most enriched Afut4_977-sac1 in azole-resistant strains decreased azole resistance and downregulated its nearby gene, sac1, but it did not significantly affect the expression of genes of the ergosterol synthesis pathway. We then discovered that 5'LTR of Afut4_977-sac1 had promoter activity and enhanced the adjacent sac1 gene expression. We found that sac1 is important to A. fumigatus, and the upregulated sac1 caused elevated resistance of A. fumigatus to azoles. Finally, we showed that Afut4_977-sac1 has an evolution pattern similar to that of the whole genome of azole-resistant strains due to azoles; phylogenetic analysis of both the whole genome and Afut4_977-sac1 suggests that the insertion of Afut4_977-sac1 might have preceded the emergence of azole-resistant strains. Taking these data together, we found that the Afut4_977-sac1 LTR-RT might be involved in the regulation of azole resistance of A. fumigatus by upregulating its nearby sac1 gene.

Embryonic LTR retrotransposons supply promoter modules to somatic tissues

Long terminal repeat (LTR) retrotransposons are widely distributed across the human genome. They have accumulated through retroviral integration into germline DNA and are latent genetic modules. Active LTR promoters are observed in germline cells; however, little is known about the mechanisms underlying their active transcription in somatic tissues. Here, by integrating our previous transcriptome data set with publicly available data sets, we show that the LTR families MLT2A1 and MLT2A2 are primarily expressed in human four-cell and eight-cell embryos and are also activated in some adult somatic tissues, particularly pineal gland. Three MLT2A elements function as the promoters and first exons of the protein-coding genes ABCE1, COL5A1, and GALNT13 specifically in the pineal gland of humans but not in that of macaques, suggesting that the exaptation of these LTRs as promoters occurred during recent primate evolution. This analysis provides insight into the possible transition from germline insertion to somatic expression of LTR retrotransposons.

The Long Terminal Repeats of ERV6 Are Activated in Pre-Implantation Embryos of Cynomolgus Monkey

Precise gene regulation is critical during embryo development. Long terminal repeat elements (LTRs) of endogenous retroviruses (ERVs) are dynamically expressed in blastocysts of mammalian embryos. However, the expression pattern of LTRs in monkey blastocyst is still unknown. By single-cell RNA-sequencing (seq) data of cynomolgus monkeys, we found that LTRs of several ERV families, including MacERV6, MacERV3, MacERV2, MacERVK1, and MacERVK2, were highly expressed in pre-implantation embryo cells including epiblast (EPI), trophectoderm (TrB), and primitive endoderm (PrE), but were depleted in post-implantation. We knocked down MacERV6-LTR1a in cynomolgus monkeys with a short hairpin RNA (shRNA) strategy to examine the potential function of MacERV6-LTR1a in the early development of monkey embryos. The silence of MacERV6-LTR1a mainly postpones the differentiation of TrB, EPI, and PrE cells in embryos at day 7 compared to control. Moreover, we confirmed MacERV6-LTR1a could recruit Estrogen Related Receptor Beta (ESRRB), which plays an important role in the maintenance of self-renewal and pluripotency of embryonic and trophoblast stem cells through different signaling pathways including FGF and Wnt signaling pathways. In summary, these results suggest that MacERV6-LTR1a is involved in gene regulation of the pre-implantation embryo of the cynomolgus monkeys.

Recombination marks the evolutionary dynamics of a recently endogenized retrovirus

All vertebrate genomes have been colonized by retroviruses along their evolutionary trajectory. Although endogenous retroviruses (ERVs) can contribute important physiological functions to contemporary hosts, such benefits are attributed to long-term coevolution of ERV and host because germline infections are rare and expansion is slow, and because the host effectively silences them. The genomes of several outbred species including mule deer (Odocoileus hemionus) are currently being colonized by ERVs, which provides an opportunity to study ERV dynamics at a time when few are fixed. We previously established the locus-specific distribution of cervid ERV (CrERV) in populations of mule deer. In this study, we determine the molecular evolutionary processes acting on CrERV at each locus in the context of phylogenetic origin, genome location, and population prevalence. A mule deer genome was de novo assembled from short- and long-insert mate pair reads and CrERV sequence generated at each locus. We report that CrERV composition and diversity have recently measurably increased by horizontal acquisition of a new retrovirus lineage. This new lineage has further expanded CrERV burden and CrERV genomic diversity by activating and recombining with existing CrERV. Resulting interlineage recombinants then endogenize and subsequently expand. CrERV loci are significantly closer to genes than expected if integration were random and gene proximity might explain the recent expansion of one recombinant CrERV lineage. Thus, in mule deer, retroviral colonization is a dynamic period in the molecular evolution of CrERV that also provides a burst of genomic diversity to the host population.

Transposable elements and their KZFP controllers are drivers of transcriptional innovation in the developing human brain

Transposable elements (TEs) account for more than 50% of the human genome and many have been co-opted throughout evolution to provide regulatory functions for gene expression networks. Several lines of evidence suggest that these networks are fine-tuned by the largest family of TE controllers, the KRAB-containing zinc finger proteins (KZFPs). One tissue permissive for TE transcriptional activation (termed "transposcription") is the adult human brain, however comprehensive studies on the extent of this process and its potential contribution to human brain development are lacking. To elucidate the spatiotemporal transposcriptome of the developing human brain, we have analyzed two independent RNA-seq data sets encompassing 16 brain regions from eight weeks postconception into adulthood. We reveal a distinct KZFP:TE transcriptional profile defining the late prenatal to early postnatal transition, and the spatiotemporal and cell type-specific activation of TE-derived alternative promoters driving the expression of neurogenesis-associated genes. Long-read sequencing confirmed these TE-driven isoforms as significant contributors to neurogenic transcripts. We also show experimentally that a co-opted antisense L2 element drives temporal protein relocalization away from the endoplasmic reticulum, suggestive of novel TE dependent protein function in primate evolution. This work highlights the widespread dynamic nature of the spatiotemporal KZFP:TE transcriptome and its importance throughout TE mediated genome innovation and neurotypical human brain development. To facilitate interactive exploration of these spatiotemporal gene and TE expression dynamics, we provide the "Brain TExplorer" web application freely accessible for the community.

RepExpress: A Novel Pipeline for the Quantification and Characterization of Transposable Element Expression from RNA-seq Data

Transposable elements (TEs) are key regulators of both development and disease; however, their repetitive nature presents substantial computational challenges to their analysis. Due to a lack of computational tools and suitable analysis frameworks, TE expression is often not quantified at the locus level. Therefore, we have developed RepExpress, a novel pipeline that enables locus-level TE quantification and characterization. RepExpress enables the characterization of TE expression in a genomic context, and is the first tool focusing on the identification of tissue-specific TE-derived and TE-regulated genes. RepExpress identifies expressed TEs overlapping with annotated genomic features and enables tissue-specific profiles of TE-derived genes. TEs that are expressed with no overlap with any known genomic features are characterized by the closest downstream genomic feature enabling identification of novel TE-gene regulatory relationships. RepExpress takes standard RNA-seq data as input and performs genomic alignment optimized for TEs. Our novel pipeline quantifies expression of both TEs and genes using featureCounts and Stringtie, respectively. RepExpress then filters expressed repeats and characterizes their genomic context, enabling the identification of TEs that overlap with genes, or that may be influencing gene expression. Here, we describe RepExpress, and provide a step-by-step protocol detailing its workflow. We also discuss other TE analysis tools and their applicability to addressing different biological questions. © 2021 Wiley Periodicals LLC. Basic Protocol: RepExpress workflow.

TFs for TEs: the transcription factor repertoire of mammalian transposable elements

In this review, Hermant and Torres-Padilla summarize and discuss the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific transposable element families or subfamilies.

Transposable elements (TEs) are genetic elements capable of changing position within the genome. Although their mobilization can constitute a threat to genome integrity, nearly half of modern mammalian genomes are composed of remnants of TE insertions. The first critical step for a successful transposition cycle is the generation of a full-length transcript. TEs have evolved cis-regulatory elements enabling them to recruit host-encoded factors driving their own, selfish transcription. TEs are generally transcriptionally silenced in somatic cells, and the mechanisms underlying their repression have been extensively studied. However, during germline formation, preimplantation development, and tumorigenesis, specific TE families are highly expressed. Understanding the molecular players at stake in these contexts is of utmost importance to establish the mechanisms regulating TEs, as well as the importance of their transcription to the biology of the host. Here, we review the transcription factors known to be involved in the sequence-specific recognition and transcriptional activation of specific TE families or subfamilies. We discuss the diversity of TE regulatory elements within mammalian genomes and highlight the importance of TE mobilization in the dispersal of transcription factor-binding sites over the course of evolution.

Endogenous retroviruses drive lineage-specific regulatory evolution across primate and rodent placentae

In mammals, the placenta mediates maternal-fetal nutrient and waste exchange and acts in an immunomodulatory way to facilitate maternal-fetal tolerance. The placenta is highly diverse across mammalian species, yet the molecular mechanisms that distinguish the placenta of human from other mammals are not fully understood. Using an interspecies transcriptomic comparison of human, macaque, and mouse late-gestation placentae, we identified hundreds of genes with lineage-specific expression-including dozens that are placentally-enriched and potentially related to pregnancy. We further annotated the enhancers for different human tissues using epigenomic data and demonstrate that the placenta and chorion are unique in that their enhancers display the least conservation. We identified numerous lineage-specific human placental enhancers and found they highly overlap with specific families of endogenous retroviruses (ERVs), including MER21A, MER41A/B and MER39B that were previously linked to immune response and placental function. Among these ERV families, we further demonstrate that MER41A/B insertions create dozens of lineage-specific Serum Response Factor (SRF) binding loci in human, including one adjacent to FBN2, a placenta-specific gene with increased expression in humans that produces the peptide hormone placensin to stimulate glucose secretion and trophoblast invasion. Overall, our results demonstrate the prevalence of lineage-specific placental enhancers which are frequently associated with ERV insertions and likely facilitate the lineage-specific evolution of the mammalian placenta.

Expression of Retroelements in Cervical Cancer and Their Interplay with HPV Infection and Host Gene Expression

Retroelements are expressed in diverse types of cancer and are related to tumorigenesis and to cancer progression. We characterized the expression of retroelements in cervical cancer and explored their interplay with HPV infection and their association with expression of neighboring genes. Forty biopsies of invasive cervical carcinoma (squamous cell carcinomas and adenocarcinomas) with genotyped HPV were selected and analyzed for human endogenous retrovirus (HERV) and long interspersed nuclear element 1 (L1) expression through RNA-seq data. We found 8060 retroelements expressed in the samples and a negative correlation of DNA methyltransferase 1 expression with the two most expressed L1 elements. A total of 103 retroelements were found differentially expressed between tumor histological types and between HPV types, including several HERV families (HERV-K, HERV-H, HERV-E, HERV-I and HERV-L). The comparison between HPV mono- and co-infections showed the highest proportion of differentially expressed L1 elements. The location of retroelements affected neighboring gene expression, such as shown for the interleukin-20 gene family. Three HERVs and seven L1 were located close to this gene family and two L1 showed a positive association with IL20RB expression. This study describes the expression of retroelements in cervical cancer and shows their association with HPV status and host gene expression.

A Human Endogenous Bornavirus-Like Nucleoprotein Encodes a Mitochondrial Protein Associated with Cell Viability

Endogenous retroviruses (ERVs) are sequences in animal genomes that originated from ancient retrovirus infections; they provide genetic novelty in hosts by being coopted as functional genes or elements during evolution. Recently, we demonstrated that endogenous elements from not only from retroviruses but also nonretroviral RNA viruses are a possible source of functional genes in host animals. The remnants of ancient bornavirus infections, called endogenous bornavirus-like elements (EBLs), are present in the genomes of a wide variety of vertebrate species, and some express functional products in host cells. Previous studies have predicted that the human EBL locus derived from bornavirus nucleoprotein, termed hsEBLN-2, expresses mRNA encoding a protein, suggesting that hsEBLN-2 has acquired a cellular function during evolution. However, the detailed function of the hsEBLN-2-derived product remains to be elucidated. In this study, we show that the hsEBLN-2-derived protein E2 acts as a mitochondrial protein that interacts with mitochondrial host factors associated with apoptosis, such as HAX-1. We also demonstrate that knockdown of hsEBLN-2-derived RNA increased the levels of PARP and caspase-3 cleavage and markedly decreased cell viability. In contrast, overexpression of E2 enhanced cell viability, as well as the intracellular stability of HAX-1, under stress conditions. Our results suggest that hsEBLN-2 has been coopted as a host gene, the product of which is involved in cell viability by interacting with mitochondrial proteins. IMPORTANCE Our genomes contain molecular fossils of ancient viruses, called endogenous virus elements (EVEs). Mounting evidence suggests that EVEs derived from nonretroviral RNA viruses have acquired functions in host cells during evolution. Previous studies have revealed that a locus encoding a bornavirus-derived EVE, hsEBLN-2, which was generated approximately 43 million years ago in a human ancestor, may be linked to the development of some tumors. However, the function of hsEBLN-2 has not been determined. In this study, we found that the E2 protein, an expression product of hsEBLN-2, interacts with apoptosis-related host proteins as a mitochondrial protein and affects cell viability. This study suggests that nonretroviral RNA viral EVEs have been coopted by hosts with more diverse functions than previously thought, showing a pivotal role for RNA virus infection in evolution.

Functions and Regulation of Endogenous Retrovirus Elements during Zygotic Genome Activation: Implications for Improving Somatic Cell Nuclear Transfer Efficiency

Endogenous retroviruses (ERVs), previously viewed as deleterious relics of ancestral retrovirus infections, are silenced in the vast majority of cells to minimize the risk of retrotransposition. Counterintuitively, bursts of ERV transcription usually occur during maternal-to-zygotic transition (MZT) in preimplantation embryos; this is regarded as a major landmark event in the zygotic genome activation (ZGA) process, indicating that ERVs play an active part in ZGA. Evolutionarily, the interaction between ERVs and hosts is mutually beneficial. The endogenization of retrovirus sequences rewires the gene regulatory network during ZGA, and ERV repression may lower germline fitness. Unfortunately, owing to various limitations of somatic cell nuclear transfer (SCNT) technology, both developmental arrest and ZGA abnormalities occur in a high percentage of cloned embryos, accompanied by ERV silencing, which may be caused by the activation failure of upstream ERV inducers. In this review, we discuss the functions and regulation of ERVs during the ZGA process and the feasibility of temporal control over ERVs in cloned embryos via exogenous double homeobox (DUX). We hypothesize that further accurate characterization of the ERV-rewired gene regulatory network during ZGA may provide a novel perspective on the development of preimplantation embryos.

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.