Human transposable elements in Repbase: genomic footprints from fish to humans

REPrise: de novo interspersed repeat detection using inexact seeding

BACKGROUND

Interspersed repeats occupy a large part of many eukaryotic genomes, and thus their accurate annotation is essential for various genome analyses. Database-free de novo repeat detection approaches are powerful for annotating genomes that lack well-curated repeat databases. However, existing tools do not yet have sufficient repeat detection performance.

RESULTS

In this study, we developed REPrise, a de novo interspersed repeat detection software program based on a seed-and-extension method. Although the algorithm of REPrise is similar to that of RepeatScout, which is currently the de facto standard tool, we incorporated three unique techniques into REPrise: inexact seeding, affine gap scoring and loose masking. Analyses of rice and simulation genome datasets showed that REPrise outperformed RepeatScout in terms of sensitivity, especially when the repeat sequences contained many mutations. Furthermore, when applied to the complete human genome dataset T2T-CHM13, REPrise demonstrated the potential to detect novel repeat sequence families.

CONCLUSION

REPrise can detect interspersed repeats with high sensitivity even in long genomes. Our software enhances repeat annotation in diverse genomic studies, contributing to a deeper understanding of genomic structures.

Aberrant expression of human endogenous retrovirus K9-derived elements is associated with better clinical outcome of acute myelocytic leukemia

BACKGROUND

Acute myelocytic leukemia (AML) is a common hematological malignancy in adults. Although several risk stratifications based on cytogenetic and molecular abnormalities are available to guide the indications for allogeneic hematopoietic cell transplantation (allo-HCT), determining optimal treatment strategies for AML remains challenging. In this study, using transcriptome datasets, we investigated the association between event-free survival (EFS) in intensively treated AML patients and the aberrant expression of endogenous viral element (EVE)-derived open reading frames (ORFs), which have been reported to be associated with the pathophysiology of various malignancies and have the potential to serve as neoantigens in specific cancers.

RESULTS

The expression levels of human endogenous retrovirus family K9 (HERVK9) ORFs were associated with EFS, independent of conventional risk stratification. Furthermore, AML cells with higher levels of HERVK9 expression exhibited enhanced antigen processing and presentation, along with increased expression of genes associated with adaptive immune responses and apoptosis, indicating that aberrant HERVK9 expression may initiate an anti-neoplastic immune response via increased antigen presentation.

CONCLUSIONS

HERVK9 expression may have serve as a crucial prognostic indicator that could aid in determining the indications for upfront allo-HCT in AML patients.

Further varieties of ancient endogenous retrovirus in human DNA

A retrovirus inserts its genome into the DNA of a cell, occasionally a germ-line cell that gives rise to descendants of the host organism: it is then called an endogenous retrovirus (ERV). The human genome contains relics from many kinds of ancient ERV. Some relics contributed new genes and regulatory elements. This study finds further kinds of ancient ERV, in the thoroughly-studied human genome version hg38: ERV-Hako, ERV-Saru, ERV-Hou, ERV-Han, and ERV-Goku. It also finds many relics of ERV-V, previously known from just two copies on chromosome 19 with placental genes. It finds a type of ERV flanked by MER41E long terminal repeats (LTRs), with surprisingly little similarity to the known MER41 ERV. ERV-Hako has subtypes that contain sequence from host genes SUSD6 and SPHKAP: the SUSD6 variant was transferred between catarrhine and platyrrhine primates. A retrovirus uses tRNA to prime reverse transcription: Hako is the only human ERV relic that used tRNA-Trp (tryptophan, symbol W), and HERV-W is misnamed because it used tRNA-Arg, based on the Genomic tRNA Database. One ERV-Saru LTR is the previously-described enhancer of AIM2 in innate immunity. This study contributes to understanding primate ERV history, but also shows that related ERVs can have drastic differences, challenging the goal of clearly annotating all ERV relics in genomes.

Inverted Alu repeats in loop-out exon skipping across hominoid evolution

Background

Changes in RNA splicing over the course of evolution have profoundly diversified the functional landscape of the human genome. While DNA sequences proximal to intron-exon junctions are known to be critical for RNA splicing, the impact of distal intronic sequences remains underexplored. Emerging evidence suggests that inverted pairs of intronic Alu elements can promote exon skipping by forming RNA stem-loop structures. However, their prevalence and influence throughout evolution remain unknown.

Results

Here, we present a systematic analysis of inverted Alu pairs across the human genome to assess their impact on exon skipping through predicted RNA stem-loop formation and their relevance to hominoid evolution. We found that inverted Alu pairs, particularly pairs of AluY-AluSx1 and AluSz-AluSx, are enriched in the flanking regions of skippable exons genome-wide and are predicted to form stable stem-loop structures. Exons defined by weak 3' acceptor and strong 5' donor splice sites appear especially prone to this skipping mechanism. Through comparative genome analysis across nine primate species, we identified 67,126 hominoid-specific Alu insertions, primarily from AluY and AluS subfamilies, which form inverted pairs enriched across skippable exons in genes of ubiquitination-related pathways. Experimental validation of exon skipping among several hominoid-specific inverted Alu pairs further reinforced their potential evolutionary significance.

Conclusion

This work extends our current knowledge of the roles of RNA secondary structure formed by inverted Alu pairs and details a newly emerging mechanism through which transposable elements have contributed to genomic innovation across hominoid evolution at the transcriptomic level.

Massively parallel reporter assay investigates shared genetic variants of eight psychiatric disorders

A meta-genome-wide association study across eight psychiatric disorders has highlighted the genetic architecture of pleiotropy in major psychiatric disorders. However, mechanisms underlying pleiotropic effects of the associated variants remain to be explored. We conducted a massively parallel reporter assay to decode the regulatory logic of variants with pleiotropic and disorder-specific effects. Pleiotropic variants differ from disorder-specific variants by exhibiting chromatin accessibility that extends across diverse cell types in the neuronal lineage and by altering motifs for transcription factors with higher connectivity in protein-protein interaction networks. We mapped pleiotropic and disorder-specific variants to putative target genes using functional genomics approaches and CRISPR perturbation. In vivo CRISPR perturbation of a pleiotropic and a disorder-specific gene suggests that pleiotropy may involve the regulation of genes expressed broadly across neuronal cell types and with higher network connectivity.

Flipons enable genomes to learn by intermediating the exchange of energy for information

Recent findings have confirmed the long-held belief that alternative DNA conformations encoded by genetic elements called flipons have important biological roles. Many of these alternative structures are formed by sequences originally spread throughout the human genome by endogenous retroelements (ERE) that captured 50% of the territory before being disarmed. Only 2.6% of the remaining DNA codes for proteins. Other organisms have instead streamlined their genomes by eliminating invasive retroelements and other repeat elements. The question arises, why retain any ERE at all? A new synthesis suggests that flipons enable genomes to learn and programme the context-specific readout of information by altering the transcripts produced. The exchange of energy for information is mediated through changes in DNA topology. Here I provide a formulation for how genomes learn and describe the underlying p-bit algorithm through which flipons are tuned. The framework suggests new strategies for the therapeutic reprogramming of cells.

Chromosome-scale genome assembly reveals how repeat elements shape non-coding RNA landscapes active during newt limb regeneration

Newts have large genomes harboring many repeat elements. How these elements shape the genome and relate to newts' unique regeneration ability remains unknown. We present here the chromosome-scale assembly of the 20.3 Gb genome of the Iberian ribbed newt, Pleurodeles waltl, with a hitherto unprecedented contiguity and completeness among giant genomes. Utilizing this assembly, we demonstrate conserved synteny as well as genetic rearrangements, such as in the major histocompatibility complex locus. We provide evidence suggesting that intronic repeat elements drive newt-specific circular RNA (circRNA) biogenesis and show their regeneration-specific expression. We also present a comprehensive in-depth annotation and chromosomal mapping of microRNAs, highlighting genomic expansion profiles as well as a distinct regulatory pattern in the regenerating limb. These data reveal links between repeat elements, non-coding RNAs, and adult regeneration and provide key resources for addressing developmental, regenerative, and evolutionary principles.

Transposable element activity captures human pluripotent cell states

Human pluripotent stem cells (hPSCs) exist in multiple, transcriptionally distinct states and serve as powerful models for studying human development. Despite their significance, the molecular determinants and pathways governing these pluripotent states remain incompletely understood. Here, we demonstrate that transposable elements act as sensitive indicators of distinct pluripotent cell states. We engineered hPSCs with fluorescent reporters to capture the temporal expression dynamics of two state-specific transposable elements, LTR5_Hs, and MER51B. This dual reporter system enables real-time monitoring and isolation of stem cells transitioning from naïve to primed pluripotency and further towards differentiation, serving as a more accurate readout of pluripotency states compared to conventional systems. Unexpectedly, we identified a rare, metastable cell population within primed hPSCs, marked by transcripts related to preimplantation embryo development and which is associated with a DNA damage response. Moreover, our system establishes the chromatin factor NSD1 and the RNA-binding protein FUS as potent molecular safeguards of primed pluripotency. Our study introduces a novel system for investigating cellular potency and provides key insights into the regulation of embryonic development.

HERV reactivation by adenovirus infection is associated with viral immune regulation

Human endogenous retroviruses (HERVs), which are normally silenced by methylation or mutation, can be reactivated by a variety of environmental factors, including infection with exogenous viruses. In this work, we investigated the transcriptional activity of HERVs following infection of human liver cells (HepaRG) with human adenovirus C serotype 5 (HAdV-C5). HAdV-C5 infection results in reactivation of several HERV groups as well as differentially expressed genes. Interestingly, in HAdV-C5 infection, upregulated genes that were in close chromosomal proximity to upregulated HERV loci were associated with influencing viral carcinogenesis and inflammatory signaling. We also identified an FBXO17 transcript encoding an intronic ERVK9-11 sense sequence upon HAdV-C5 infection. FBXO17 has previously been described as an important factor in the regulation of the interferon response. This suggests that specific HERV groups may have the potential to trigger gene networks and influence viral immune responses.

Genome-Wide Tool for Sensitive de novo Identification and Visualisation of Interspersed and Tandem Repeats

Genomic repeats are functionally ubiquitous structural units found in all genomes. Studying these repeats of different origins is essential for understanding the evolution and adaptation of a given organism. These repeating patterns have manifold signatures and structures with varying degrees of homology, making their identification challenging. To address this challenge, we developed a new algorithm and software that can rapidly and accurately detect any repeated sequences de novo with varying degrees of homology in genomic sequences in interspersed or clustered repeats. Numerous forms of repeated sequences and complex patterns can be identified, even for complex sequence variants and implicit or mixed types of repeat blocks. Direct and inverted-repeat elements, perfect and imperfect microsatellite repeats, and any short or long tandem repeat belonging to a wide range of higher-order repeat structures of telomeres or large satellite sequences can be detected. By combining precision and versatility, our tool contributes significantly to elucidating the intricate landscape of genomic repeats.

Cancer cells restrict immunogenicity of retrotransposon expression via distinct mechanisms

To thrive, cancer cells must navigate acute inflammatory signaling accompanying oncogenic transformation, such as via overexpression of repeat elements. We examined the relationship between immunostimulatory repeat expression, tumor evolution, and the tumor-immune microenvironment. Integration of multimodal data from a cohort of pancreatic ductal adenocarcinoma (PDAC) patients revealed expression of specific Alu repeats predicted to form double-stranded RNAs (dsRNAs) and trigger retinoic-acid-inducible gene I (RIG-I)-like-receptor (RLR)-associated type-I interferon (IFN) signaling. Such Alu-derived dsRNAs also anti-correlated with pro-tumorigenic macrophage infiltration in late stage tumors. We defined two complementary pathways whereby PDAC may adapt to such anti-tumorigenic signaling. In mutant TP53 tumors, ORF1p from long interspersed nuclear element (LINE)-1 preferentially binds Alus and decreases their expression, whereas adenosine deaminases acting on RNA 1 (ADAR1) editing primarily reduces dsRNA formation in wild-type TP53 tumors. Depletion of either LINE-1 ORF1p or ADAR1 reduced tumor growth in vitro. The fact that tumors utilize multiple pathways to mitigate immunostimulatory repeats implies the stress from their expression is a fundamental phenomenon to which PDAC, and likely other tumors, adapt.

Landscape of evolutionary arms races between transposable elements and KRAB-ZFP family

Transposable elements (TEs) are mobile parasitic sequences that have expanded within the host genome. It has been hypothesized that host organisms have expanded the Krüppel-associated box-containing zinc finger proteins (KRAB-ZFPs), which epigenetically suppress TEs, to counteract disorderly TE transpositions. This process is referred to as the evolutionary arms race. However, the extent to which this evolutionary arms race occurred across various TE families remains unclear. In the present study, we systematically explored the evolutionary arms race between TE families and human KRAB-ZFPs using public ChIP-seq data. We discovered and characterized new instances of evolutionary arms races with KRAB-ZFPs in endogenous retroviruses. Furthermore, we found that the regulatory landscape shaped by this arms race contributed to the gene regulatory networks. In summary, our results provide insight into the impact of the evolutionary arms race on TE families, the KRAB-ZFP family, and host gene regulatory networks.

Evolution of Repetitive Elements, Their Roles in Homeostasis and Human Disease, and Potential Therapeutic Applications

Repeating sequences of DNA, or repetitive elements (REs), are common features across both prokaryotic and eukaryotic genomes. Unlike many of their protein-coding counterparts, the functions of REs in host cells remained largely unknown and have often been overlooked. While there is still more to learn about their functions, REs are now recognized to play significant roles in both beneficial and pathological processes in their hosts at the cellular and organismal levels. Therefore, in this review, we discuss the various types of REs and review what is known about their evolution. In addition, we aim to classify general mechanisms by which REs promote processes that are variously beneficial and harmful to host cells/organisms. Finally, we address the emerging role of REs in cancer, aging, and neurological disorders and provide insights into how RE modulation could provide new therapeutic benefits for these specific conditions.

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.

Dysregulation of Transposon Transcription Profiles in Cancer Cells Resembles That of Embryonic Stem Cells

Transposable elements (TEs) comprise a substantial portion of the mammalian genome, with potential implications for both embryonic development and cancer. This study aimed to characterize the expression profiles of TEs in embryonic stem cells (ESCs), cancer cell lines, tumor tissues, and the tumor microenvironment (TME). We observed similarities in TE expression profiles between cancer cells and ESCs, suggesting potential parallels in regulatory mechanisms. Notably, four TE RNAs (HERVH, LTR7, HERV-Fc1, HERV-Fc2) exhibited significant downregulation across cancer cell lines and tumor tissues compared to ESCs, highlighting potential roles in pluripotency regulation. The strong up-regulation of the latter two TEs (HERV-Fc1, HERV-Fc2) in ESCs has not been previously demonstrated and may be a first indication of their role in the regulation of pluripotency. Conversely, tandemly repeated sequences (MSR1, CER, ALR) showed up-regulation in cancer contexts. Moreover, a difference in TE expression was observed between the TME and the tumor bulk transcriptome, with distinct dysregulated TE profiles. Some TME-specific TEs were absent in normal tissues, predominantly belonging to LTR and L1 retrotransposon families. These findings not only shed light on the regulatory roles of TEs in both embryonic development and cancer but also suggest novel targets for anti-cancer therapy. Understanding the interplay between cancer cells and the TME at the TE level may pave the way for further research into therapeutic interventions.

Recombinant origin and interspecies transmission of a HERV-K(HML-2)-related primate retrovirus with a novel RNA transport element

HERV-K(HML-2), the youngest clade of human endogenous retroviruses (HERVs), includes many intact or nearly intact proviruses, but no replication competent HML-2 proviruses have been identified in humans. HML-2-related proviruses are present in other primates, including rhesus macaques, but the extent and timing of HML-2 activity in macaques remains unclear. We have identified 145 HML-2-like proviruses in rhesus macaques, including a clade of young, rhesus-specific insertions. Age estimates, intact open reading frames, and insertional polymorphism of these insertions are consistent with recent or ongoing infectious activity in macaques. 106 of the proviruses form a clade characterized by an ~750 bp sequence between env and the 3' long terminal repeat (LTR), derived from an ancient recombination with a HERV-K(HML-8)-related virus. This clade is found in Old World monkeys (OWM), but not great apes, suggesting it originated after the ape/OWM split. We identified similar proviruses in white-cheeked gibbons; the gibbon insertions cluster within the OWM recombinant clade, suggesting interspecies transmission from OWM to gibbons. The LTRs of the youngest proviruses have deletions in U3, which disrupt the Rec Response Element (RcRE), required for nuclear export of unspliced viral RNA. We show that the HML-8-derived region functions as a Rec-independent constitutive transport element (CTE), indicating the ancestral Rec-RcRE export system was replaced by a CTE mechanism.

Endogenous retroviruses are dysregulated in ALS

Amyotrophic lateral sclerosis (ALS) is a universally fatal neurodegenerative disease with no cure. Human endogenous retroviruses (HERVs) have been implicated in its pathogenesis but their relevance to ALS is not fully understood. We examined bulk RNA-seq data from almost 2,000 ALS and unaffected control samples derived from the cortex and spinal cord. Using different methods of feature selection, including differential expression analysis and machine learning, we discovered that transcription of HERV-K loci 1q22 and 8p23.1 were significantly upregulated in the spinal cord of individuals with ALS. Additionally, we identified a subset of ALS patients with upregulated HERV-K expression in the cortex and spinal cord. We also found the expression of HERV-K loci 19q11 and 8p23.1 was correlated with protein coding genes previously implicated in ALS and dysregulated in ALS patients in this study. These results clarify the association of HERV-K and ALS and highlight specific genes in the pathobiology of late-stage ALS.

Conserved and divergent DNA recognition specificities and functions of R2 retrotransposon N-terminal domains

R2 non-long terminal repeat (non-LTR) retrotransposons are among the most extensively distributed mobile genetic elements in multicellular eukaryotes and show promise for applications in transgene supplementation of the human genome. They insert new gene copies into a conserved site in 28S ribosomal DNA with exquisite specificity. R2 clades are defined by the number of zinc fingers (ZFs) at the N terminus of the retrotransposon-encoded protein, postulated to additively confer DNA site specificity. Here, we illuminate general principles of DNA recognition by R2 N-terminal domains across and between clades, with extensive, specific recognition requiring only one or two compact domains. DNA-binding and protection assays demonstrate broadly shared as well as clade-specific DNA interactions. Gene insertion assays in cells identify the N-terminal domains sufficient for target-site insertion and reveal roles in second-strand cleavage or synthesis for clade-specific ZFs. Our results have implications for understanding evolutionary diversification of non-LTR retrotransposon insertion mechanisms and the design of retrotransposon-based gene therapies.

Long-read sequencing improves the genetic diagnosis of retinitis pigmentosa by identifying an Alu retrotransposon insertion in the EYS gene

BACKGROUND

Biallelic variants in EYS are the major cause of autosomal recessive retinitis pigmentosa (arRP) in certain populations, a clinically and genetically heterogeneous disease that may lead to legal blindness. EYS is one of the largest genes (~ 2 Mb) expressed in the retina, in which structural variants (SVs) represent a common cause of disease. However, their identification using short-read sequencing (SRS) is not always feasible. Here, we conducted targeted long-read sequencing (T-LRS) using adaptive sampling of EYS on the MinION sequencing platform (Oxford Nanopore Technologies) to definitively diagnose an arRP family, whose affected individuals (n = 3) carried the heterozygous pathogenic deletion of exons 32-33 in the EYS gene. As this was a recurrent variant identified in three additional families in our cohort, we also aimed to characterize the known deletion at the nucleotide level to assess a possible founder effect.

RESULTS

T-LRS in family A unveiled a heterozygous AluYa5 insertion in the coding exon 43 of EYS (chr6(GRCh37):g.64430524_64430525ins352), which segregated with the disease in compound heterozygosity with the previously identified deletion. Visual inspection of previous SRS alignments using IGV revealed several reads containing soft-clipped bases, accompanied by a slight drop in coverage at the Alu insertion site. This prompted us to develop a simplified program using grep command to investigate the recurrence of this variant in our cohort from SRS data. Moreover, LRS also allowed the characterization of the CNV as a ~ 56.4kb deletion spanning exons 32-33 of EYS (chr6(GRCh37):g.64764235_64820592del). The results of further characterization by Sanger sequencing and linkage analysis in the four families were consistent with a founder variant.

CONCLUSIONS

To our knowledge, this is the first report of a mobile element insertion into the coding sequence of EYS, as a likely cause of arRP in a family. Our study highlights the value of LRS technology in characterizing and identifying hidden pathogenic SVs, such as retrotransposon insertions, whose contribution to the etiopathogenesis of rare diseases may be underestimated.

Endogenous retroviruses in development and health

Endogenous retroviruses (ERVs) are evolutionary remnants of retroviral infections in which the viral genome became embedded as a dormant regulatory element within the host germline. When ERVs become activated, they comprehensively rewire genomic regulatory networks of the host and facilitate critical developmental events, such as preimplantation development and placentation, in a manner specific to species, developmental stage, and tissues. However, accumulating evidence suggests that aberrant ERV transcription compromises genome stability and has been implicated in cellular senescence and various pathogenic processes, underscoring the significance of host genomic surveillance mechanisms. Here, we revisit the prominent functions of ERVs in early development and highlight their emerging roles in mammalian post-implantation development and organogenesis. We also discuss their implications for aging and pathological processes such as microbial infection, immune response. Furthermore, we discuss recent advances in stem-cell-based models, single-cell omics, and genome editing technologies, which serve as beacons illuminating the versatile nature of ERVs in mammalian development and health.

Transposable Elements: Emerging Therapeutic Targets in Neurodegenerative Diseases

Neurodegenerative diseases, such as Alzheimer's disease (AD), Parkinson's disease (PD), and amyotrophic lateral sclerosis (ALS), are characterized by the progressive loss of neuronal function and structure. While several genetic and environmental factors have been implicated in the pathogenesis of these disorders, emerging evidence suggests that transposable elements (TEs), once considered "junk DNA," play a significant role in their development and progression. TEs are mobile genetic elements capable of moving within the genome, and their dysregulation has been associated with genomic instability, altered gene expression, and neuroinflammation. This review provides an overview of TEs, including long interspersed nuclear elements (LINEs), short interspersed nuclear elements (SINEs), and endogenous retroviruses (ERVs), mechanisms of repression and derepression, and their potential impact on neurodegeneration. The evidence linking TEs to AD, PD, and ALS by shedding light on the complex interactions between TEs and neurodegeneration has been discussed. Furthermore, the therapeutic potential of targeting TEs in neurodegenerative diseases has been explored. Understanding the role of TEs in neurodegeneration holds promise for developing novel therapeutic strategies aimed at mitigating disease progression and preserving neuronal health.

Statistical learning quantifies transposable element-mediated cis-regulation

BACKGROUND

Transposable elements (TEs) have colonized the genomes of most metazoans, and many TE-embedded sequences function as cis-regulatory elements (CREs) for genes involved in a wide range of biological processes from early embryogenesis to innate immune responses. Because of their repetitive nature, TEs have the potential to form CRE platforms enabling the coordinated and genome-wide regulation of protein-coding genes by only a handful of trans-acting transcription factors (TFs).

RESULTS

Here, we directly test this hypothesis through mathematical modeling and demonstrate that differences in expression at protein-coding genes alone are sufficient to estimate the magnitude and significance of TE-contributed cis-regulatory activities, even in contexts where TE-derived transcription fails to do so. We leverage hundreds of overexpression experiments and estimate that, overall, gene expression is influenced by TE-embedded CREs situated within approximately 500 kb of promoters. Focusing on the cis-regulatory potential of TEs within the gene regulatory network of human embryonic stem cells, we find that pluripotency-specific and evolutionarily young TE subfamilies can be reactivated by TFs involved in post-implantation embryogenesis. Finally, we show that TE subfamilies can be split into truly regulatorily active versus inactive fractions based on additional information such as matched epigenomic data, observing that TF binding may better predict TE cis-regulatory activity than differences in histone marks.

CONCLUSION

Our results suggest that TE-embedded CREs contribute to gene regulation during and beyond gastrulation. On a methodological level, we provide a statistical tool that infers TE-dependent cis-regulation from RNA-seq data alone, thus facilitating the study of TEs in the next-generation sequencing era.

Annotation of transposable elements in the transcriptome of the Neotropical brown stink bug Euschistus heros and its chromosomal distribution

Transposable elements (TEs) are DNA sequences capable of moving within the genome. Their distribution is very dynamic among organisms, and despite advances, there are still gaps in the understanding of the diversity and evolution of TEs in many insect species. In the case of Euschistus heros, considered the main stink bug in the soybean crop in Brazil, little is known about the participation of these elements. Therefore, the objective of the current work was to identify the different groups of transposable elements present in the E. heros transcriptome, evidencing their chromosomal distribution. Through RNA-Seq and de novo assembly, 60,009 transcripts were obtained, which were annotated locally via Blastn against specific databases. Of the 367 transcripts identified as TEs, 202 belong to Class II, with emphasis on the TIR order. Among Class I elements or retrotransposons, most were characterized as LINE. Phylogenetic analyses were performed with the protein domains, evidencing differences between Tc1-mariner sequences, which may be related to possible horizontal transfer events. The transposable elements that stood out in the transcriptome were selected for fluorescent in situ hybridization. DNA transposon probes hAT, Helitron, and Tc1-mariner showed mostly scattered signals, with the presence of some blocks. Retrotransposon probes Copia, Gypsy, Jockey, and RTE showed a more pulverized hybridization pattern, with the presence of small interstitial and/or terminal blocks. Studies like this one, integrating functional genomics and molecular cytogenetic tools, are essential to expanding knowledge about transcriptionally active mobile elements, and their behavior in the chromosomes.

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

Genetic features and genomic targets of human KRAB-zinc finger proteins

Krüppel-associated box (KRAB) domain-containing zinc finger proteins (KZFPs) are one of the largest groups of transcription factors encoded by tetrapods, with 378 members in human alone. KZFP genes are often grouped in clusters reflecting amplification by gene and segment duplication since the gene family first emerged more than 400 million years ago. Previous work has revealed that many KZFPs recognize transposable element (TE)-embedded sequences as genomic targets, and that KZFPs facilitate the co-option of the regulatory potential of TEs for the benefit of the host. Here, we present a comprehensive survey of the genetic features and genomic targets of human KZFPs, notably completing past analyses by adding data on close to a hundred family members. General principles emerge from our study of the TE-KZFP regulatory system, which point to multipronged evolutionary mechanisms underlaid by highly complex and combinatorial modes of action with strong influences on human speciation.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry

Overexpression of repetitive elements is an emerging hallmark of human cancers ¹ . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system ^2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

Improved DNA-Versus-Protein Homology Search for Protein Fossils

Protein fossils, i.e., noncoding DNA descended from coding DNA, arise frequently from transposable elements (TEs), decayed genes, and viral integrations. They can reveal, and mislead about, evolutionary history and relationships. They have been detected by comparing DNA to protein sequences, but current methods are not optimized for this task. We describe a powerful DNA-protein homology search method. We use a 64×21 substitution matrix, which is fitted to sequence data, automatically learning the genetic code. We detect subtly homologous regions by considering alternative possible alignments between them, and calculate significance (probability of occurring by chance between random sequences). Our method detects TE protein fossils much more sensitively than blastx, and faster. Of the  ∼ 7 major categories of eukaryotic TE, three were long thought absent in mammals: we find two of them in the human genome, polinton and DIRS/Ngaro. This method increases our power to find ancient fossils, and perhaps to detect non-standard genetic codes. The alternative-alignments and significance paradigm is not specific to DNA-protein comparison, and could benefit homology search generally. This is an extended version of a conference paper (Yao & Frith, 2021).

Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon

BACKGROUND

Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE).

RESULTS

We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association.

CONCLUSIONS

Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

Regulation of endogenous retrovirus-derived regulatory elements by GATA2/3 and MSX2 in human trophoblast stem cells

The placenta is an organ with extraordinary phenotypic diversity in eutherian mammals. Recent evidence suggests that numerous human placental enhancers are evolved from lineage-specific insertions of endogenous retroviruses (ERVs), yet the transcription factors (TFs) underlying their regulation remain largely elusive. Here, by first focusing on MER41, a primate-specific ERV family previously linked to placenta and innate immunity, we uncover the binding motifs of multiple crucial trophoblast TFs (GATA2/3, MSX2, GRHL2) in addition to innate immunity TFs STAT1 and IRF1. Integration of ChIP-seq data confirms the binding of GATA2/3, MSX2, and their related factors on the majority of MER41-derived enhancers in human trophoblast stem cells (TSCs). MER41-derived enhancers that are constitutively active in human TSCs are distinct from those activated upon interferon stimulation, which is determined by the binding of relevant TFs and their subfamily compositions. We further demonstrate that GATA2/3 and MSX2 have prevalent binding to numerous other ERV families - indicating their broad impact on ERV-derived enhancers. Functionally, the derepression of many syncytiotrophoblast genes after MSX2 knockdown is likely to be mediated by regulatory elements derived from ERVs - suggesting ERVs are also important for mediating transcriptional repression. Overall, this study characterizes the regulation of ERV-derived regulatory elements by GATA2/3, MSX2, and their cofactors in human TSCs, and provides mechanistic insights into the importance of ERVs in human trophoblast regulatory network.

Targets of histone H3 lysine 9 methyltransferases

Histone H3 lysine 9 di- and trimethylation are well-established marks of constitutively silenced heterochromatin domains found at repetitive DNA elements including pericentromeres, telomeres, and transposons. Loss of heterochromatin at these sites causes genomic instability in the form of aberrant DNA repair, chromosome segregation defects, replication stress, and transposition. H3K9 di- and trimethylation also regulate cell type-specific gene expression during development and form a barrier to cellular reprogramming. However, the role of H3K9 methyltransferases extends beyond histone methylation. There is a growing list of non-histone targets of H3K9 methyltransferases including transcription factors, steroid hormone receptors, histone modifying enzymes, and other chromatin regulatory proteins. Additionally, two classes of H3K9 methyltransferases modulate their own function through automethylation. Here we summarize the structure and function of mammalian H3K9 methyltransferases, their roles in genome regulation and constitutive heterochromatin, as well as the current repertoire of non-histone methylation targets including cases of automethylation.

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.

A Systems Biology Approach on the Regulatory Footprint of Human Endogenous Retroviruses (HERVs)

Human endogenous retroviruses (HERVs) are a family of endogenous retroviruses that comprise the ~8.93% of the human genome sequence, with a high proportion being human specific. The recent expansion of repeated HERV sequences has offered a framework for genetic and epigenetic innovation. In the current report, a systematic approach is implemented to catalogue regulatory elements within HERVs, as a roadmap to potential functions of HERV sequences in gene networks. ENCODE Project has offered a wealth of epigenetic data based on omics technologies. I analyzed the presence of HERV sequences on consensus cis-regulatory elements (cCREs) from ENCODE data. On the one side, HERVs are in 1 out of 9 cCREs (>100.000 cCREs in total), dispersed within the genome and present in cis-regulatory regions of ~81% of human genes, as calculated following gene enrichment analysis. On the other side, promoter-associated HERV cCREs are present adjacent to (in a 200 bp window) the transcription start sites of 256 human genes. Regulatory network production, followed by centrality analysis led to the discovery of 90 core genes containing HERV-associated promoters. Pathway analysis on the core network genes and their immediate neighbors revealed a regulatory footprint that, among others, is associated with inflammation, chemokine signaling and response to viral infection. Collectively, these results support the concept that the expansion of regulatory sequences derived from HERVs is critical for epigenetic innovation that may have wired together genes into novel transcriptional networks with critical roles in cellular physiology and pathology.

Molecular diversity and phenotypic pleiotropy of ancient genomic regulatory loci derived from human endogenous retrovirus type H (HERVH) promoter LTR7 and HERVK promoter LTR5_Hs and their contemporary impacts on pathophysiology of Modern Humans

Timelines of population-level effects of viruses on humans varied from the evolutionary scale of million years to contemporary spread of viral infections. Correspondingly, these events are exemplified by: (i) emergence of human endogenous retroviruses (HERVs) from ancient germline infections leading to stable integration of viral genomes into human chromosomes; and (ii) wide-spread viral infections reaching a global pandemic state such as the COVID-19 pandemic. Despite significant efforts, understanding of HERV’s roles in governance of genomic regulatory networks, their impacts on primate evolution and development of human-specific physiological and pathological phenotypic traits remains limited. Remarkably, present analyses revealed that expression of a dominant majority of genes (1696 of 1944 genes; 87%) constituting high-confidence down-steam regulatory targets of defined HERV loci was significantly altered in cells infected with the SARS-CoV-2 coronavirus, a pathogen causing the global COVID-19 pandemic. This study focused on defined sub-sets of DNA sequences derived from HERVs that are expressed at specific stages of human preimplantation embryogenesis and exert regulatory actions essential for self-renewal and pluripotency. Evolutionary histories of LTR7/HERVH and LTR5_Hs/HERVK were charted based on evidence of the earliest presence and expansion of highly conserved (HC) LTR sequences. Sequence conservation analyses of most recent releases 17 primate species’ genomes revealed that LTR7/HERVH have entered germlines of primates in Africa after the separation of the New World Monkey lineage, while LTR5_Hs/HERVK successfully colonized primates’ germlines after the segregation of Gibbons’ species. Subsequently, both LTR7 and LTR5_Hs undergo a marked ~ fourfold–fivefold expansion in genomes of Great Apes. Timelines of quantitative expansion of both LTR7 and LTR5_Hs loci during evolution of Great Apes appear to replicate the consensus evolutionary sequence of increasing cognitive and behavioral complexities of non-human primates, which seems particularly striking for LTR7 loci and 11 distinct LTR7 subfamilies. Consistent with previous reports, identified in this study, 351 human-specific (HS) insertions of LTR7 (175 loci) and LTR5_Hs (176 loci) regulatory sequences have been linked to genes implicated in establishment and maintenance of naïve and primed pluripotent states and preimplantation embryogenesis phenotypes. Unexpectedly, HS-LTRs manifest regulatory connectivity to genes encoding markers of 12 distinct cells’ populations of fetal gonads, as well as genes implicated in physiology and pathology of human spermatogenesis, including Y-linked spermatogenic failure, oligo- and azoospermia. Granular interrogations of genes linked with 11 distinct LTR7 subfamilies revealed that mammalian offspring survival (MOS) genes seem to remain one of consistent regulatory targets throughout ~ 30 MYA of the divergent evolution of LTR7 loci. Differential GSEA of MOS versus non-MOS genes identified clearly discernable dominant enrichment patterns of phenotypic traits affected by MOS genes linked with LTR7 (562 MOS genes) and LTR5_Hs (126 MOS genes) regulatory loci across the large panel of genomics and proteomics databases reflecting a broad spectrum of human physiological and pathological traits. GSEA of LTR7-linked MOS genes identified more than 2200 significantly enriched records of human common and rare diseases and gene signatures of 466 significantly enriched records of Human Phenotype Ontology traits, including Autosomal Dominant (92 genes) and Autosomal Recessive (93 genes) Inheritance. LTR7 regulatory elements appear linked with genes implicated in functional and morphological features of central nervous system, including synaptic transmission and protein–protein interactions at synapses, as well as gene signatures differentially regulated in cells of distinct neurodevelopmental stages and morphologically diverse cell types residing and functioning in human brain. These include Neural Stem/Precursor cells, Radial Glia cells, Bergman Glia cells, Pyramidal cells, Tanycytes, Immature neurons, Interneurons, Trigeminal neurons, GABAergic neurons, and Glutamatergic neurons. GSEA of LTR7-linked genes identified significantly enriched gene sets encoding markers of more than 80 specialized types of neurons and markers of 521 human brain regions, most prominently, subiculum and dentate gyrus. Identification and characterization of 1944 genes comprising high-confidence down-steam regulatory targets of LTR7 and/or LTR5_Hs loci validated and extended these observations by documenting marked enrichments for genes implicated in neoplasm metastasis, intellectual disability, autism, multiple cancer types, Alzheimer’s, schizophrenia, and other brain disorders. Overall, genes representing down-stream regulatory targets of ancient retroviral LTRs exert the apparently cooperative and exceedingly broad phenotypic impacts on human physiology and pathology. This is exemplified by altered expression of 93% high-confidence LTR targets in cells infected by contemporary viruses, revealing a convergence of virus-inflicted aberrations on genomic regulatory circuitry governed by ancient retroviral LTR elements and interference with human cells’ differentiation programs.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

Influenza A Virus Infection Reactivates Human Endogenous Retroviruses Associated with Modulation of Antiviral Immunity

Human endogenous retrovirus (HERVs), normally silenced by methylation or mutations, can be reactivated by multiple environmental factors, including infections with exogenous viruses. In this work, we investigated the transcriptional activity of HERVs in human A549 cells infected by two wild-type (PR8M, SC35M) and one mutated (SC35MΔNS1) strains of Influenza A virus (IAVs). We found that the majority of differentially expressed HERVs (DEHERVS) and genes (DEGs) were up-regulated in the infected cells, with the most significantly enriched biological processes associated with the genes differentially expressed exclusively in SC35MΔNS1 being linked to the immune system. Most DEHERVs in PR8M and SC35M are mammalian apparent LTR retrotransposons, while in SC35MΔNS1, more HERV loci from the HERVW9 group were differentially expressed. Furthermore, up-regulated pairs of HERVs and genes in close chromosomal proximity to each other tended to be associated with immune responses, which implies that specific HERV groups might have the potential to trigger specific gene networks and influence host immunological pathways.

Role of Transposable Elements in Genome Stability: Implications for Health and Disease

Most living organisms have in their genome a sizable proportion of DNA sequences capable of mobilization; these sequences are commonly referred to as transposons, transposable elements (TEs), or jumping genes. Although long thought to have no biological significance, advances in DNA sequencing and analytical technologies have enabled precise characterization of TEs and confirmed their ubiquitous presence across all forms of life. These findings have ignited intense debates over their biological significance. The available evidence now supports the notion that TEs exert major influence over many biological aspects of organismal life. Transposable elements contribute significantly to the evolution of the genome by giving rise to genetic variations in both active and passive modes. Due to their intrinsic nature of mobility within the genome, TEs primarily cause gene disruption and large-scale genomic alterations including inversions, deletions, and duplications. Besides genomic instability, growing evidence also points to many physiologically important functions of TEs, such as gene regulation through cis-acting control elements and modulation of the transcriptome through epigenetic control. In this review, we discuss the latest evidence demonstrating the impact of TEs on genome stability and the underling mechanisms, including those developed to mitigate the deleterious impact of TEs on genomic stability and human health. We have also highlighted the potential therapeutic application of TEs.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

Paleozoic Protein Fossils Illuminate the Evolution of Vertebrate Genomes and Transposable Elements

Genomes hold a treasure trove of protein fossils: fragments of formerly protein-coding DNA, which mainly come from transposable elements (TEs) or host genes. These fossils reveal ancient evolution of TEs and genomes, and many fossils have been exapted to perform diverse functions important for the host's fitness. However, old and highly-degraded fossils are hard to identify, standard methods (e.g. BLAST) are not optimized for this task, and few Paleozoic protein fossils have been found. Here, a recently optimized method is used to find protein fossils in vertebrate genomes. It finds Paleozoic fossils predating the amphibian/amniote divergence from most major TE categories, including virus-related Polinton and Gypsy elements. It finds 10 fossils in the human genome (8 from TEs and 2 from host genes) that predate the last common ancestor of all jawed vertebrates, probably from the Ordovician period. It also finds types of transposon and retrotransposon not found in human before. These fossils have extreme sequence conservation, indicating exaptation: some have evidence of gene-regulatory function, and they tend to lienearest to developmental genes. Some ancient fossils suggest "genome tectonics", where two fragments of one TE have drifted apart by up to megabases, possibly explaining gene deserts and large introns. This paints a picture of great TE diversity in our aquatic ancestors, with patchy TE inheritance by later vertebrates, producing new genes and regulatory elements on the way. Host-gene fossils too have contributed anciently-conserved DNA segments. This paves the way to further studies of ancient protein fossils.

Epigenetic regulation of human non-coding RNA gene transcription

Recent investigations on the non-protein-coding transcriptome of human cells have revealed previously hidden layers of gene regulation relying on regulatory non-protein-coding (nc) RNAs, including the widespread ncRNA-dependent regulation of epigenetic chromatin states and of mRNA translation and stability. However, despite its centrality, the epigenetic regulation of ncRNA genes has received relatively little attention. In this mini-review, we attempt to provide a synthetic account of recent literature suggesting an unexpected complexity in chromatin-dependent regulation of ncRNA gene transcription by the three human nuclear RNA polymerases. Emerging common features, like the heterogeneity of chromatin states within ncRNA multigene families and their influence on 3D genome organization, point to unexplored issues whose investigation could lead to a better understanding of the whole human epigenomic network.

Mosaic cis-regulatory evolution drives transcriptional partitioning of HERVH endogenous retrovirus in the human embryo

The human endogenous retrovirus type-H (HERVH) family is expressed in the preimplantation embryo. A subset of these elements are specifically transcribed in pluripotent stem cells where they appear to exert regulatory activities promoting self-renewal and pluripotency. How HERVH elements achieve such transcriptional specificity remains poorly understood. To uncover the sequence features underlying HERVH transcriptional activity, we performed a phyloregulatory analysis of the long terminal repeats (LTR7) of the HERVH family, which harbor its promoter, using a wealth of regulatory genomics data. We found that the family includes at least eight previously unrecognized subfamilies that have been active at different timepoints in primate evolution and display distinct expression patterns during human embryonic development. Notably, nearly all HERVH elements transcribed in ESCs belong to one of the youngest subfamilies we dubbed LTR7up. LTR7 sequence evolution was driven by a mixture of mutational processes, including point mutations, duplications, and multiple recombination events between subfamilies, that led to transcription factor binding motif modules characteristic of each subfamily. Using a reporter assay, we show that one such motif, a predicted SOX2/3 binding site unique to LTR7up, is essential for robust promoter activity in induced pluripotent stem cells. Together these findings illuminate the mechanisms by which HERVH diversified its expression pattern during evolution to colonize distinct cellular niches within the human embryo.

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

Antisense overlapping long non-coding RNA regulates coding arachidonate 12-lipoxygenase gene by translational interference

The arachidonate 12-lipoxygenase (ALOX12) enzyme catalyzes polyunsaturated fatty acids and facilitates generation of bioactive lipid mediators associated with various biological processes and disease pathologies. The human genome assembly revealed that the ALOX12 gene overlaps an antisense non-coding gene designated as ALOX12-antisense 1 (ALOX12-AS1). This arrangement indicates that the uncharacterized ALOX12-AS1 long non-coding RNA (lncRNA) may bind to the sense coding ALOX12 mRNA to form an antisense-sense duplex providing the basis of a novel ALOX12 regulatory mechanism. Therefore, this study was designed to determine whether the interaction of ALOX12-AS1 with ALOX12 mRNA functions as an anti-sense/sense duplex-mediated regulatory mechanism controlling the cellular content of ALOX12. Our findings indicate that two major isoforms of ALOX12-AS1 lncRNA are ubiquitously expressed in a variety of primary adult human tissues and different transformed cell types. RNA-FISH revealed cell-type-specific cytosolic as well as nuclear and nucleolar localization of the lncRNA. Interestingly, phorbol ester-induced nucleo-cytoplasmic translocation of the lncRNA in monocytic THP-1 cells resulted in a reduction of ALOX12 protein without a concomitant change in its mRNA level. This indicated ALOX12-AS1 operates via an antisense-sense duplex-mediated translational downregulation mechanism. This deduction was validated by demonstrating sense/antisense duplex formation and an association of the duplex with ribosomal proteins in HEK293 cells. Overall, this study revealed a hitherto unknown mechanism of antisense lncRNA-mediated translational downregulation of ALOX12 that adds to the existing regulatory mechanisms for the modulation of potent bioactive lipid mediators that contribute to both health and disease.

Transcriptome Analyses Implicate Endogenous Retroviruses Involved in the Host Antiviral Immune System through the Interferon Pathway

Human endogenous retroviruses (HERVs) are the remains of ancient retroviruses that invaded our ancestors' germline cell and were integrated into the genome. The expression of HERVs has always been a cause for concern because of its association with various cancers and diseases. However, few previous studies have focused on specific activation of HERVs by viral infections. Our previous study has shown that dengue virus type 2 (DENV-2) infection induces the transcription of a large number of abnormal HERVs loci; therefore, the purpose of this study was to explore the relationship between exogenous viral infection and HERV activation further. In this study, we retrieved and reanalyzed published data on 21 transcriptomes of human cells infected with various viruses. We found that infection with different viruses could induce transcriptional activation of HERV loci. Through the comparative analysis of all viral datasets, we identified 43 key HERV loci that were up-regulated by DENV-2, influenza A virus, influenza B virus, Zika virus, measles virus, and West Nile virus infections. Furthermore, the neighboring genes of these HERVs were simultaneously up-regulated, and almost all such neighboring genes were interferon-stimulated genes (ISGs), which are enriched in the host's antiviral immune response pathways. Our data supported the hypothesis that activation of HERVs, probably via an interferon-mediated mechanism, plays an important role in innate immunity against viral infections.

Genome-wide bioinformatic analyses predict key host and viral factors in SARS-CoV-2 pathogenesis

The novel betacoronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused a worldwide pandemic (COVID-19) after emerging in Wuhan, China. Here we analyzed public host and viral RNA sequencing data to better understand how SARS-CoV-2 interacts with human respiratory cells. We identified genes, isoforms and transposable element families that are specifically altered in SARS-CoV-2-infected respiratory cells. Well-known immunoregulatory genes including CSF2, IL32, IL-6 and SERPINA3 were differentially expressed, while immunoregulatory transposable element families were upregulated. We predicted conserved interactions between the SARS-CoV-2 genome and human RNA-binding proteins such as the heterogeneous nuclear ribonucleoprotein A1 (hnRNPA1) and eukaryotic initiation factor 4 (eIF4b). We also identified a viral sequence variant with a statistically significant skew associated with age of infection, that may contribute to intracellular host-pathogen interactions. These findings can help identify host mechanisms that can be targeted by prophylactics and/or therapeutics to reduce the severity of COVID-19.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

Human Endogenous Retroviruses in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most aggressive and deadly brain tumor. It is primarily diagnosed in the elderly and has a 5-year survival rate of less than 6% even with the most aggressive therapies. The lack of biomarkers has made the development of immunotherapy for GBM challenging. Human endogenous retroviruses (HERVs) are a group of viruses with long terminal repeat (LTR) elements, which are believed to be relics from ancient viral infections. Recent studies have found that those repetitive elements play important roles in regulating various biological processes. The differentially expressed LTR elements from HERVs are potential biomarkers for immunotherapy to treat GBM. However, the understanding of the LTR element expression in GBM is greatly lacking.

METHODS

We obtained 1077.4 GB of sequencing data from public databases. These data were generated from 111 GBM tissue studies, 30 GBM cell lines studies, and 45 normal brain tissues studies. We analyzed repetitive elements that were differentially expressed in GBM and normal brain samples.

RESULTS

We found that 48 LTR elements were differentially expressed (p-value < 0.05) between GBM and normal brain tissues, of which 46 were HERV elements. Among these 46 elements, 34 significantly changed HERVs belong to the ERV1 superfamily. Furthermore, 43 out of the 46 differentially expressed HERV elements were upregulated.

CONCLUSION

Our results indicate significant differential expression of many HERV LTR elements in GBM and normal brain tissues. Expression levels of these elements could be developed as biomarkers for GBM treatments.

Hamster PIWI proteins bind to piRNAs with stage-specific size variations during oocyte maturation

In animal gonads, transposable elements are actively repressed to preserve genome integrity through the PIWI-interacting RNA (piRNA) pathway. In mice, piRNAs are abundantly expressed in male germ cells, and form effector complexes with three distinct PIWIs. The depletion of individual Piwi genes causes male-specific sterility with no discernible phenotype in female mice. Unlike mice, most other mammals have four PIWI genes, some of which are expressed in the ovary. Here, purification of PIWI complexes from oocytes of the golden hamster revealed that the size of the PIWIL1-associated piRNAs changed during oocyte maturation. In contrast, PIWIL3, an ovary-specific PIWI in most mammals, associates with short piRNAs only in metaphase II oocytes, which coincides with intense phosphorylation of the protein. An improved high-quality genome assembly and annotation revealed that PIWIL1- and PIWIL3-associated piRNAs appear to share the 5'-ends of common piRNA precursors and are mostly derived from unannotated sequences with a diminished contribution from TE-derived sequences, most of which correspond to endogenous retroviruses. Our findings show the complex and dynamic nature of biogenesis of piRNAs in hamster oocytes, and together with the new genome sequence generated, serve as the foundation for developing useful models to study the piRNA pathway in mammalian oocytes.