Interspersed repeats and other mementos of transposable elements in mammalian genomes

DUX4 activates common and context-specific intergenic transcripts and isoforms

DUX4 regulates the expression of genic and nongenic elements and modulates chromatin accessibility during zygotic genome activation in cleavage stage embryos. Its misexpression in skeletal muscle causes facioscapulohumeral dystrophy (FSHD). By leveraging full-length RNA isoform sequencing with short-read RNA sequencing of DUX4-inducible myoblasts, we elucidate an isoform-resolved transcriptome featuring numerous unannotated isoforms from known loci and novel intergenic loci. While DUX4 activates similar programs in early embryos and FSHD muscle, the isoform usage of known DUX4 targets is notably distinct between the two contexts. DUX4 also activates hundreds of previously unannotated intergenic loci dominated by repetitive elements. The transcriptional and epigenetic profiles of these loci in myogenic and embryonic contexts indicate that the usage of DUX4-binding sites at these intergenic loci is influenced by the cellular environment. These findings demonstrate that DUX4 induces context-specific transcriptomic programs, enriching our understanding of DUX4-induced muscle pathology.

Posttranscriptional RNA stabilization of telomeric RNAs FRG2, DBE-T, D4Z4 at human 4q35 in response to genotoxic stress and D4Z4 macrosatellite repeat length

BACKGROUND

Reduced copy number of the D4Z4 macrosatellite at human chromosome 4q35 is associated with facioscapulohumeral muscular dystrophy (FSHD). A pervasive idea is that chromatin alterations at the 4q35 locus following D4Z4 repeat unit deletion lead to disease via inappropriate expression of nearby genes. Here, we sought to analyze transcription and chromatin characteristics at specific regions of 4q35 and how these are affected by D4Z4 deletions and exogenous stresses.

RESULTS

We found that the 4q subtelomere is subdivided into discrete domains, each with characteristic chromatin features associated with distinct gene expression profiles. Centromeric genes within 4q35 (SLC25A4, FAT1 and FRG1) display active histone marks at their promoters. In contrast, poised or repressed markings are present at telomeric loci including FRG2, DBE-T and D4Z4. We discovered that these discrete domains undergo region-specific chromatin changes upon treatment with chromatin enzyme inhibitors or genotoxic drugs. We demonstrated that the 4q35 telomeric FRG2, DBE-T and D4Z4-derived transcripts are induced upon DNA damage to levels inversely correlated with the D4Z4 repeat number, are stabilized through posttranscriptional mechanisms upon DNA damage and are bound to chromatin.

CONCLUSION

Our study reveals unforeseen biochemical features of RNAs from clustered transcription units within the 4q35 subtelomere. Specifically, the FRG2, DBE-T and D4Z4-derived transcripts are chromatin-associated and are stabilized posttranscriptionally after induction by genotoxic stress. Remarkably, the extent of this response is modulated by the copy number of the D4Z4 repeats, raising new hypotheses about their regulation and function in human biology and disease.

Chromatin accessibility landscape of mouse early embryos revealed by single-cell NanoATAC-seq2

In mammals, fertilized eggs undergo genome-wide epigenetic reprogramming to generate the organism. However, our understanding of epigenetic dynamics during preimplantation development at single-cell resolution remains incomplete. Here, we developed scNanoATAC-seq2, a single-cell assay for transposase-accessible chromatin using long-read sequencing for scarce samples. We present a detailed chromatin accessibility landscape of mouse preimplantation development, revealing distinct chromatin signatures in the epiblast, primitive endoderm, and trophectoderm during lineage segregation. Differences between zygotes and two-cell embryos highlight reprogramming in chromatin accessibility during the maternal-to-zygotic transition. Single-cell long-read sequencing enables in-depth analysis of chromatin accessibility in noncanonical imprinting, imprinted X chromosome inactivation, and low-mappability genomic regions, such as repetitive elements and paralogs. Our data provide insights into chromatin dynamics during mammalian preimplantation development and lineage differentiation.

The domesticated transposon protein L1TD1 associates with its ancestor L1 ORF1p to promote LINE-1 retrotransposition

Repression of retrotransposition is crucial for the successful fitness of a mammalian organism. The domesticated transposon protein L1TD1, derived from LINE-1 (L1) ORF1p, is an RNA-binding protein that is expressed only in some cancers and early embryogenesis. In human embryonic stem cells, it is found to be essential for maintaining pluripotency. In cancer, L1TD1 expression is highly correlative with malignancy progression and as such considered a potential prognostic factor for tumors. However, its molecular role in cancer remains largely unknown. Our findings reveal that DNA hypomethylation induces the expression of L1TD1 in HAP1 human tumor cells. L1TD1 depletion significantly modulates both the proteome and transcriptome and thereby reduces cell viability. Notably, L1TD1 associates with L1 transcripts and interacts with L1 ORF1p protein, thereby facilitating L1 retrotransposition. Our data suggest that L1TD1 collaborates with its ancestral L1 ORF1p as an RNA chaperone, ensuring the efficient retrotransposition of L1 retrotransposons, rather than directly impacting the abundance of L1TD1 targets. In this way, L1TD1 might have an important role not only during early development but also in tumorigenesis.

Syntenic lncRNA locus exhibits DNA regulatory functions with sequence evolution

Syntenic long non-coding RNAs (lncRNAs) often show limited sequence conservation across species, prompting concern in the field. This study delves into functional signatures of syntenic lncRNAs between humans and zebrafish. Syntenic lncRNAs are highly expressed in zebrafish, with ∼90 % located near protein-coding genes, either in sense or antisense orientation. During early zebrafish development and in human embryonic stem cells (H1-hESC), syntenic lncRNA loci are enriched with cis-regulatory repressor signatures, influencing the expression of development-associated genes. In later zebrafish developmental stages and specific human cell lines, these syntenic lncRNA loci function as enhancers or transcription start sites (TSS) for protein-coding genes. Analysis of transposable elements (TEs) in syntenic lncRNA sequences revealed intriguing patterns: human lncRNAs are enriched in simple repeat elements, while their zebrafish counterparts show enrichment in LTR elements. This sequence evolution likely arises from post-rearrangement mutations that enhance DNA elements or cis-regulatory functions. It may also contribute to vertebrate innovation by creating novel transcription factor binding sites within the locus. This study highlights the conserved functionality of syntenic lncRNA loci through DNA elements, emphasizing their conserved roles across species despite sequence divergence.

SINE retrotransposons import polyadenylation signals to 3'UTRs in dog (Canis familiaris)

BACKGROUND

Messenger RNA 3' untranslated regions (3'UTRs) control many aspects of gene expression and determine where the transcript will terminate. The polyadenylation signal (PAS) AAUAAA (AATAAA in DNA) is a key regulator of transcript termination and this hexamer, or a similar sequence, is very frequently found within 30 bp of 3'UTR ends. Short interspersed element (SINE) retrotransposons are found throughout genomes in high copy numbers. When inserted into genes they can disrupt expression, alter splicing, or cause nuclear retention of mRNAs. The genomes of the domestic dog and other carnivores carry hundreds of thousands of Can-SINEs, a tRNA-related SINE with transcription termination potential. Because of this we asked whether Can-SINEs may terminate transcript in some dog genes.

RESULTS

Each of the dog's nine Can-SINE consensus sequences carry an average of three AATAAA PASs on their sense strands but zero on their antisense strands. Consistent with the idea that Can-SINEs can terminate transcripts, we find that sense-oriented Can-SINEs are approximately ten times more frequent at 3' ends of 3'UTRs compared to further upstream within 3'UTRs. Furthermore, the count of AATAAA PASs on head-to-tail SINE sequences differs significantly between sense and antisense-oriented retrotransposons in transcripts. Can-SINEs near 3'UTR ends are likely to carry an AATAAA motif on the mRNA sense strand while those further upstream are not. We identified loci where Can-SINE insertion has truncated or altered a 3'UTR of the dog genome (dog 3'UTR) compared to the human ortholog. Dog 3'UTRs have peaks of AATAAA PAS frequency at 28, 32, and 36 bp from the end. The periodicity is partly explained by TAAA(n) repeats within Can-SINE AT-rich tails. We annotated all repeat-masked Can-SINE copies in the Boxer reference genome and found that the young SINEC_Cf type has a mode of 15 bp length for target site duplications (TSDs). All dog Can-SINE types favor integration at TSDs beginning with A(4).

CONCLUSION

Dog Can-SINE retrotransposition has imported AATAAA PASs into gene transcripts and led to alteration of 3'UTRs. AATAAA sequences are selectively removed from Can-SINEs in introns and upstream 3'UTR regions but are retained at the far downstream end of 3'UTRs, which we infer reflects their role as termination sequences for these transcripts.

The reconstruction of evolutionary dynamics of processed pseudogenes indicates deep silencing of "retrobiome" in naked mole rat

Approximately half of mammalian genomes are occupied by retrotransposons, highly repetitive interspersed genetic elements expanded through the mechanism of reverse transcription. The evolution of this "retrobiome" involved a series of explosive amplifications, presumably associated with high mutation rates, interspersed with periods of silencing. A by-product of retrotransposon activity is the formation of processed pseudogenes (PPGs)-intron-less, promoter-less DNA copies of messenger RNA (mRNA). We examined the proportion of PPGs with varying degrees of deviation from their ancestor mRNAs as an indicator of the intensity of retrotranspositions at different times in the past. Our analysis revealed a high proportion of "young'' (recently acquired) PPGs in the DNA of mice and rats, indicating significant retrobiome activity during the recent evolution of these species. The ongoing process of new PPG entries in mouse germ line DNA was confirmed by identifying diversity in PPG content within the single strain of mice, C57BL/6. In contrast, the highly abundant PPGs of the naked mole rat (NMR) exhibited substantial deviation from their mRNAs, with a near-complete lack of PPGs without mutations, indicative of the silencing of the retrobiome in the most recent evolutionary past, preceded by a period of high activity. This distinctive feature of the NMR genome was confirmed through the analysis of a broad range of mammalian species. The peculiar evolutionary dynamics of PPGs in the NMR, an organism with exceptional longevity and resistance to cancer, may reflect the role played by the retrobiome in aging and cancer.

Enhancing mRNA Interactions by Engineering the Arc Protein with Nucleocapsid Domain

Activity-regulated cytoskeleton-associated protein (Arc) forms virus-like capsids for mRNA transport between neurons. Unlike HIV-1 Group-specific Antigen (Gag), which uses its Nucleocapsid (NC) domain to bind HIV-1 genomic mRNA, mammalian Arc lacks the NC domain, and their direct mRNA binding interactions remain underexplored. This study examined rat Arc's binding to rat Arc 5' UTR (A5U), HIV-1 5' UTR (H5U), and GFP mRNAs, revealing weak binding with no significant preference. Adding the HIV-1 NC domain to rArc's C-terminus significantly improved binding to H5U, while also showing substantial binding to A5U at about 60% of its H5U level and exhibiting twice the affinity for A5U over GFP mRNA. Importantly, rArc-NC binds 3.4 times more A5U and H5U than GST-NC, indicating that rArc NTD-CA aids mRNA binding by HIV-1 NC. These findings suggest a conserved Gag protein-mRNA interaction mechanism, highlighting the potential for developing mRNA delivery systems that combine endogenous Gag NTD-CA with retroviral NC and UTRs.

TGIRT-seq of Inflammatory Breast Cancer Tumor and Blood Samples Reveals Widespread Enhanced Transcription Impacting RNA Splicing and Intronic RNAs in Plasma

Inflammatory breast cancer (IBC) is the most aggressive and lethal breast cancer subtype but lacks unequivocal genomic differences or robust biomarkers that differentiate it from non-IBC. Here, Thermostable Group II intron Reverse Transcriptase RNA-sequencing (TGIRT-seq) revealed myriad differences in tumor samples, Peripheral Blood Mononuclear Cells (PBMCs), and plasma that distinguished IBC from non-IBC patients and healthy donors across all tested receptor-based subtypes. These included numerous differentially expressed protein-coding gene and non-coding RNAs in all three sample types, a granulocytic immune response in IBC PBMCs, and over-expression of antisense RNAs, suggesting wide-spread enhanced transcription in both IBC tumors and PBMCs. By using TGIRT-seq to quantitate Intron-exon Depth Ratios (IDRs) and mapping reads to both genome and transcriptome reference sequences, we developed methods for parallel analysis of transcriptional and post-transcriptional gene regulation. This analysis identified numerous differentially and non-differentially expressed protein-coding genes in IBC tumors and PBMCs with high IDRs, the latter reflecting rate-limiting RNA splicing that negatively impacts mRNA production. Mirroring gene expression differences in tumors and PBMCs, over-represented protein-coding gene RNAs in IBC patient plasma were largely intronic RNAs, while those in non-IBC patients and healthy donor plasma were largely mRNA fragments. Potential IBC biomarkers in plasma included T-cell receptor pre-mRNAs and intronic, LINE-1, and antisense RNAs. Our findings provide new insights into IBC and set the stage for monitoring disease progression and response to treatment by liquid biopsy. The methods developed for parallel transcriptional and post-transcriptional gene regulation analysis have potentially broad RNA-seq and clinical applications.

Characterization of hAT DNA transposon superfamily in the genome of Neotropical fish Apareiodon sp

DNA transposons are diverse in fish genomes and have been described to generate genomic evolutionary novelties. hAT transposable element data are scarce in Teleostei genomes, making it challenging to conduct comparative genomic studies to understand their neutrality or function. This study aimed to perform a genomic and molecular characterization of hAT copies to assess the diversity of these elements and associate changes in these sequences to genomic and karyotypic novelties in Apareiodon sp. The data revealed that hAT TEs are highly abundant in the Apareiodon sp. genome, with few possibly autonomous copies. Highly conserved sequences with likely functional transposases were observed in nine hAT elements. A great diversity of hAT subgroups was observed, especially from Ac, Charlie, Blackjack, Tip100, hAT6, and hAT5, and a similar wave of hAT genomic invasion was identified in the genome for these six groups of hAT sequences. The data also revealed a distinct number of microsatellites within degenerated hAT copies. hAT sites were demonstrated to be dispersed in the Apareiodon sp. chromosomes and not involved in W chromosome-specific region differentiation. In conclusion, the genomic analysis revealed a great diversity of hAT elements, possible autonomous copies, and differentiation of degenerated transposable elements into tandem sequences.

Machine Learning Reveals the Diversity of Human 3D Chromatin Contact Patterns

Understanding variation in chromatin contact patterns across diverse humans is critical for interpreting noncoding variants and their effects on gene expression and phenotypes. However, experimental determination of chromatin contact patterns across large samples is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the variation in 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We apply this approach to thousands of human genomes from the 1000 Genomes Project and the inferred hominin ancestral genome. While patterns of 3D contact divergence genome wide are qualitatively similar to patterns of sequence divergence, we find substantial differences in 3D divergence and sequence divergence in local 1 megabase genomic windows. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, for 31% of genomic windows, a single individual has a rare divergent 3D contact map pattern. Using in silico mutagenesis, we find that most single nucleotide sequence changes do not result in changes to 3D chromatin contacts. However, in windows with substantial 3D divergence just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals variable contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

Age-associated changes in innate and adaptive immunity: role of the gut microbiota

Aging is generally regarded as an irreversible process, and its intricate relationship with the immune system has garnered significant attention due to its profound implications for the health and well-being of the aging population. As people age, a multitude of alterations occur within the immune system, affecting both innate and adaptive immunity. In the realm of innate immunity, aging brings about changes in the number and function of various immune cells, including neutrophils, monocytes, and macrophages. Additionally, certain immune pathways, like the cGAS-STING, become activated. These alterations can potentially result in telomere damage, the disruption of cytokine signaling, and impaired recognition of pathogens. The adaptive immune system, too, undergoes a myriad of changes as age advances. These include shifts in the number, frequency, subtype, and function of T cells and B cells. Furthermore, the human gut microbiota undergoes dynamic changes as a part of the aging process. Notably, the interplay between immune changes and gut microbiota highlights the gut's role in modulating immune responses and maintaining immune homeostasis. The gut microbiota of centenarians exhibits characteristics akin to those found in young individuals, setting it apart from the microbiota observed in typical elderly individuals. This review delves into the current understanding of how aging impacts the immune system and suggests potential strategies for reversing aging through interventions in immune factors.

DNA methylation governs the sensitivity of repeats to restriction by the HUSH-MORC2 corepressor

The human silencing hub (HUSH) complex binds to transcripts of LINE-1 retrotransposons (L1s) and other genomic repeats, recruiting MORC2 and other effectors to remodel chromatin. How HUSH and MORC2 operate alongside DNA methylation, a central epigenetic regulator of repeat transcription, remains largely unknown. Here we interrogate this relationship in human neural progenitor cells (hNPCs), a somatic model of brain development that tolerates removal of DNA methyltransferase DNMT1. Upon loss of MORC2 or HUSH subunit TASOR in hNPCs, L1s remain silenced by robust promoter methylation. However, genome demethylation and activation of evolutionarily-young L1s attracts MORC2 binding, and simultaneous depletion of DNMT1 and MORC2 causes massive accumulation of L1 transcripts. We identify the same mechanistic hierarchy at pericentromeric α-satellites and clustered protocadherin genes, repetitive elements important for chromosome structure and neurodevelopment respectively. Our data delineate the epigenetic control of repeats in somatic cells, with implications for understanding the vital functions of HUSH-MORC2 in hypomethylated contexts throughout human development.

Isogenic patient-derived organoids reveal early neurodevelopmental defects in spinal muscular atrophy initiation

Whether neurodevelopmental defects underlie postnatal neuronal death in neurodegeneration is an intriguing hypothesis only recently explored. Here, we focus on spinal muscular atrophy (SMA), a neuromuscular disorder caused by reduced survival of motor neuron (SMN) protein levels leading to spinal motor neuron (MN) loss and muscle wasting. Using the first isogenic patient-derived induced pluripotent stem cell (iPSC) model and a spinal cord organoid (SCO) system, we show that SMA SCOs exhibit abnormal morphological development, reduced expression of early neural progenitor markers, and accelerated expression of MN progenitor and MN markers. Longitudinal single-cell RNA sequencing reveals marked defects in neural stem cell specification and fewer MNs, favoring mesodermal progenitors and muscle cells, a bias also seen in early SMA mouse embryos. Surprisingly, SMN2-to-SMN1 conversion does not fully reverse these developmental abnormalities. These suggest that early neurodevelopmental defects may underlie later MN degeneration, indicating that postnatal SMN-increasing interventions might not completely amend SMA pathology in all patients.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

State-of-the-art in transposable element modulation affected by drugs in malignant prostatic cancer cells

Over recent years, the investigation of transposable elements (TEs) has granted researchers a deeper comprehension of their characteristics and functions, particularly regarding their significance in the mechanisms contributing to cancer development. This manuscript focuses on prostate carcinoma cell lines and offers a comprehensive review intended to scrutinize the associations and interactions between TEs and genes, as well as their response to treatment using various chemical drugs, emphasizing their involvement in cancer progression. We assembled a compendium of articles retrieved from the PubMed database to construct networks demonstrating correlations with genes and pharmaceuticals. In doing so, we linked the transposition of certain TE types to the expression of specific transcripts directly implicated in carcinogenesis. Additionally, we underline that treatment employing different drugs revealed unique patterns of TE reactivation. Our hypothesis gathers the current understanding and guides research toward evidence-based investigations, emphasizing the association between antiviral drugs, chemotherapy, and the reduced expression of TEs in patients affected by prostate cancer.

The Diverse Evolutionary Histories of Domesticated Metaviral Capsid Genes in Mammals

Selfish genetic elements comprise significant fractions of mammalian genomes. In rare instances, host genomes domesticate segments of these elements for function. Using a complete human genome assembly and 25 additional vertebrate genomes, we re-analyzed the evolutionary trajectories and functional potential of capsid (CA) genes domesticated from Metaviridae, a lineage of retrovirus-like retrotransposons. Our study expands on previous analyses to unearth several new insights about the evolutionary histories of these ancient genes. We find that at least five independent domestication events occurred from diverse Metaviridae, giving rise to three universally retained single-copy genes evolving under purifying selection and two gene families unique to placental mammals, with multiple members showing evidence of rapid evolution. In the SIRH/RTL family, we find diverse amino-terminal domains, widespread loss of protein-coding capacity in RTL10 despite its retention in several mammalian lineages, and differential utilization of an ancient programmed ribosomal frameshift in RTL3 between the domesticated CA and protease domains. Our analyses also reveal that most members of the PNMA family in mammalian genomes encode a conserved putative amino-terminal RNA-binding domain (RBD) both adjoining and independent from domesticated CA domains. Our analyses lead to a significant correction of previous annotations of the essential CCDC8 gene. We show that this putative RBD is also present in several extant Metaviridae, revealing a novel protein domain configuration in retrotransposons. Collectively, our study reveals the divergent outcomes of multiple domestication events from diverse Metaviridae in the common ancestor of placental mammals.

Population genetic analysis based on the polymorphisms mediated by transposons in the genomes of pig

Transposable elements (TEs) mobility is capable of generating a large number of structural variants (SVs), which can have considerable potential as molecular markers for genetic analysis and molecular breeding in livestock. Our results showed that the pig genome contains mainly TE-SVs generated by short interspersed nuclear elements (51,873/76.49%), followed by long interspersed nuclear elements (11,131/16.41%), and more than 84% of the common TE-SVs (Minor allele frequency, MAF > 0.10) were validated to be polymorphic. Subsequently, we utilized the identified TE-SVs to gain insights into the population structure, resulting in clear differentiation among the three pig groups and facilitating the identification of relationships within Chinese local pig breeds. In addition, we investigated the frequencies of TEs in the gene coding regions of different pig groups and annotated the respective TE types, related genes, and functional pathways. Through genome-wide comparisons of Large White pigs and Chinese local pigs utilizing the Beijing Black pigs, we identified TE-mediated SVs associated with quantitative trait loci and observed that they were mainly involved in carcass traits and meat quality traits. Lastly, we present the first documented evidence of TE transduction in the pig genome.

Retrotransposon-derived transcripts and their functions in immunity and disease

Retrotransposons, which account for approximately 42% of the human genome, have been increasingly recognized as "non-self" pathogen-associated molecular patterns (PAMPs) due to their virus-like sequences. In abnormal conditions such as cancer and viral infections, retrotransposons that are aberrantly expressed due to impaired epigenetic suppression display PAMPs, leading to their recognition by pattern recognition receptors (PRRs) of the innate immune system and triggering inflammation. This viral mimicry mechanism has been observed in various human diseases, including aging and autoimmune disorders. However, recent evidence suggests that retrotransposons possess highly regulated immune reactivity and play important roles in the development and function of the immune system. In this review, I discuss a wide range of retrotransposon-derived transcripts, their role as targets in immune recognition, and the diseases associated with retrotransposon activity. Furthermore, I explore the implications of chimeric transcripts formed between retrotransposons and known gene mRNAs, which have been previously underestimated, for the increase of immune-related gene isoforms and their influence on immune function. Retrotransposon-derived transcripts have profound and multifaceted effects on immune system function. The aim of this comprehensive review is to provide a better understanding of the complex relationship between retrotransposon transcripts and immune defense.

Mouse retrotransposons: sequence structure, evolutionary age, genomic distribution and function

Retrotransposons are transposable elements that are transposed via transcription and reverse transcription. Their copies have accumulated in the genome of mammals, occupying approximately 40% of mammalian genomic mass. These copies are often involved in numerous phenomena, such as chromatin spatial organization, gene expression, development and disease, and have been recognized as a driving force in evolution. Different organisms have gained specific retrotransposon subfamilies and retrotransposed copies, such as hundreds of Mus-specific subfamilies with diverse sequences and genomic locations. Despite this complexity, basic information is still necessary for present-day genomic and epigenomic studies. Herein, we describe the characteristics of each subfamily of Mus-specific retrotransposons in terms of sequence structure, phylogenetic relationships, evolutionary age, and preference for A or B compartments of chromatin.

Characterisation of a LINE-1 Insertion in the RP1 Gene by Targeted Adaptive Nanopore Sequencing in a Family with Retinitis Pigmentosa

Retinitis pigmentosa (RP) is a group of inherited degenerative retinal disorders affecting more than 1.5 million people worldwide. For 30-50% of individuals with RP, the genetic cause remains unresolved by current clinical diagnostic gene panels. It is likely explained by variants in novel RP-associated genes or noncoding regulatory regions, or by complex genetic alterations such as large structural variants. Recent developments in long-read sequencing techniques have opened an opportunity for efficient analysis of complex genetic variants. We analysed a Finnish family with dominantly inherited RP affecting six individuals in three generations. Two affected individuals underwent a comprehensive clinical examination in combination with a clinical diagnostic gene panel, followed by whole exome sequencing in our laboratory. They exhibited typical signs of RP, yet initial sequence analysis found no causative variants. Reanalysis of the sequencing data detected a LINE-1 (L1) retrotransposon insertion of unknown size in exon 4 of the RP1 axonemal microtubule-associated (RP1) gene. The large chimeric L1 insertion that segregated with the disease was further characterised using targeted adaptive nanopore sequencing of RP1, allowing us to identify a 5.6 kb L1 transposable element insertion in RP1 as the cause of RP in this family with dominantly inherited RP.

Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

Machine learning reveals the diversity of human 3D chromatin contact patterns

Understanding variation in chromatin contact patterns across human populations is critical for interpreting non-coding variants and their ultimate effects on gene expression and phenotypes. However, experimental determination of chromatin contacts at a population-scale is prohibitively expensive. To overcome this challenge, we develop and validate a machine learning method to quantify the diversity 3D chromatin contacts at 2 kilobase resolution from genome sequence alone. We then apply this approach to thousands of diverse modern humans and the inferred human-archaic hominin ancestral genome. While patterns of 3D contact divergence genome-wide are qualitatively similar to patterns of sequence divergence, we find that 3D divergence in local 1-megabase genomic windows does not follow sequence divergence. In particular, we identify 392 windows with significantly greater 3D divergence than expected from sequence. Moreover, 26% of genomic windows have rare 3D contact variation observed in a small number of individuals. Using in silico mutagenesis we find that most sequence changes to do not result in changes to 3D chromatin contacts. However in windows with substantial 3D divergence, just one or a few variants can lead to divergent 3D chromatin contacts without the individuals carrying those variants having high sequence divergence. In summary, inferring 3D chromatin contact maps across human populations reveals diverse contact patterns. We anticipate that these genetically diverse maps of 3D chromatin contact will provide a reference for future work on the function and evolution of 3D chromatin contact variation across human populations.

Mimicking the brain: Epstein-Barr virus and foreign agents as drivers of neuroimmune attack in multiple sclerosis

T cells have an essential role in adaptive immunity against pathogens and cancer, but failure of thymic tolerance mechanisms can instead lead to escape of T cells with the ability to attack host tissues. Multiple sclerosis (MS) occurs when structures such as myelin and neurons in the central nervous system (CNS) are the target of autoreactive immune responses, resulting in lesions in the brain and spinal cord which cause varied and episodic neurological deficits. A role for autoreactive T cell and antibody responses in MS is likely, and mounting evidence implicates Epstein-Barr virus (EBV) in disease mechanisms. In this review we discuss antigen specificity of T cells involved in development and progression of MS. We examine the current evidence that these T cells can target multiple antigens such as those from pathogens including EBV and briefly describe other mechanisms through which viruses could affect disease. Unravelling the complexity of the autoantigen T cell repertoire is essential for understanding key events in the development and progression of MS, with wider implications for development of future therapies.

The diverse evolutionary histories of domesticated metaviral capsid genes in mammals

Selfish genetic elements and their remnants comprise at least half of the human genome. Active transposons duplicate by inserting copies at new sites in a host genome. Following insertion, transposons can acquire mutations that render them inactive; the accrual of additional mutations can render them unrecognizable over time. However, in rare instances, segments of transposons become useful for the host, in a process called gene domestication. Using the first complete human genome assembly and 25 additional vertebrate genomes, we analyzed the evolutionary trajectories and functional potential of genes domesticated from the capsid genes of Metaviridae, a retroviral-like retrotransposon family. Our analysis reveals four families of domesticated capsid genes in placental mammals with varied evolutionary outcomes, ranging from universal retention to lineage-specific duplications or losses and from purifying selection to lineage-specific rapid evolution. The four families of domesticated capsid genes have divergent amino-terminal domains, inherited from four distinct ancestral metaviruses. Structural predictions reveal that many domesticated genes encode a previously unrecognized RNA-binding domain retained in multiple paralogs in mammalian genomes both adjacent to and independent from the capsid domain. Collectively, our study reveals diverse outcomes of domestication of diverse metaviruses, which led to structurally and evolutionarily diverse genes that encode important, but still largely-unknown functions in placental mammals. (207).

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Characterization of a strain-specific CD-1 reference genome reveals potential inter- and intra-strain functional variability

BACKGROUND

CD-1 is an outbred mouse stock that is frequently used in toxicology, pharmacology, and fundamental biomedical research. Although inbred strains are typically better suited for such studies due to minimal genetic variability, outbred stocks confer practical advantages over inbred strains, such as improved breeding performance and low cost. Knowledge of the full genetic variability of CD-1 would make it more useful in toxicology, pharmacology, and fundamental biomedical research.

RESULTS

We performed deep genomic DNA sequencing of CD-1 mice and used the data to identify genome-wide SNPs, indels, and germline transposable elements relative to the mm10 reference genome. We used multiple genome-wide sequencing data types and previously published CD-1 SNPs to validate our called variants. We used the called variants to construct a strain-specific CD-1 reference genome, which we show can improve mappability and reduce experimental biases from genome-wide sequencing data derived from CD-1 mice. Based on previously published ChIP-seq and ATAC-seq data, we find evidence that genetic variation between CD-1 mice can lead to alterations in transcription factor binding. We also identified a number of variants in the coding region of genes which could have effects on translation of genes.

CONCLUSIONS

We have identified millions of previously unidentified CD-1 variants with the potential to confound studies involving CD-1. We used the identified variants to construct a CD-1-specific reference genome, which can improve accuracy and reduce bias when aligning genomics data derived from CD-1 mice.

Arc protein, a remnant of ancient retrovirus, forms virus-like particles, which are abundantly generated by neurons during epileptic seizures, and affects epileptic susceptibility in rodent models

A product of the immediate early gene Arc (Activity-regulated cytoskeleton-associated protein or Arc protein) of retroviral ancestry resides in the genome of all tetrapods for millions of years and is expressed endogenously in neurons. It is a well-known protein, very important for synaptic plasticity and memory consolidation. Activity-dependent Arc expression concentrated in glutamatergic synapses affects the long-time synaptic strength of those excitatory synapses. Because it modulates excitatory-inhibitory balance in a neuronal network, the Arc gene itself was found to be related to the pathogenesis of epilepsy. General Arc knockout rodent models develop a susceptibility to epileptic seizures. Because of activity dependence, synaptic Arc protein synthesis also is affected by seizures. Interestingly, it was found that Arc protein in synapses of active neurons self-assemble in capsids of retrovirus-like particles, which can transfer genetic information between neurons, at least across neuronal synaptic boutons. Released Arc particles can be accumulated in astrocytes after seizures. It is still not known how capsid assembling and transmission timescale is affected by seizures. This scientific field is relatively novel and is experiencing swift transformation as it grapples with difficult concepts in light of evolving experimental findings. We summarize the emergent literature on the subject and also discuss the specific rodent models for studying Arc effects in epilepsy. We summarized both to clarify the possible role of Arc-related pseudo-viral particles in epileptic disorders, which may be helpful to researchers interested in this growing area of investigation.

Comparative genomics reveals insights into anuran genome size evolution

BACKGROUND

Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types.

RESULTS

Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats.

CONCLUSIONS

Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

Human DUX4 and porcine DUXC activate similar early embryonic programs in pig muscle cells: implications for preclinical models of FSHD

Human DUX4 and its mouse ortholog Dux are normally expressed in the early embryo-the 4-cell or 2-cell cleavage stage embryo, respectively-and activate a portion of the first wave of zygotic gene expression. DUX4 is epigenetically suppressed in nearly all somatic tissue, whereas facioscapulohumeral dystrophy (FSHD)-causing mutations result in its aberrant expression in skeletal muscle, transcriptional activation of the early embryonic program and subsequent muscle pathology. Although DUX4 and Dux both activate an early totipotent transcriptional program, divergence of their DNA binding domains limits the use of DUX4 expressed in mice as a preclinical model for FSHD. In this study, we identify the porcine DUXC messenger ribonucleic acid expressed in early development and show that both pig DUXC and human DUX4 robustly activate a highly similar early embryonic program in pig muscle cells. These results support further investigation of pig preclinical models for FSHD.

Insights into mammalian TE diversity through the curation of 248 genome assemblies

We examined transposable element (TE) content of 248 placental mammal genome assemblies, the largest de novo TE curation effort in eukaryotes to date. We found that although mammals resemble one another in total TE content and diversity, they show substantial differences with regard to recent TE accumulation. This includes multiple recent expansion and quiescence events across the mammalian tree. Young TEs, particularly long interspersed elements, drive increases in genome size, whereas DNA transposons are associated with smaller genomes. Mammals tend to accumulate only a few types of TEs at any given time, with one TE type dominating. We also found association between dietary habit and the presence of DNA transposon invasions. These detailed annotations will serve as a benchmark for future comparative TE analyses among placental mammals.

A Novel Alu Element Insertion in ATM Induces Exon Skipping in Suspected HBOC Patients

The vast majority of patients at risk of hereditary breast and/or ovarian cancer (HBOC) syndrome remain without a molecular diagnosis after routine genetic testing. One type of genomic alteration that is commonly missed by diagnostic pipelines is mobile element insertions (MEIs). Here, we reanalyzed multigene panel data from suspected HBOC patients using the MEI detection tool Mobster. A novel Alu element insertion in ATM intron 54 (ATM:c.8010+30_8010+31insAluYa5) was identified as a potential contributing factor in seven patients. Transcript analysis of patient-derived RNA from three heterozygous carriers revealed exon 54 skipping in 38% of total ATM transcripts. To manifest the direct association between the Alu element insertion and the aberrant splice pattern, HEK293T and MCF7 cells were transfected with wild-type or Alu element-carrying minigene constructs. On average, 77% of plasmid-derived transcripts lacked exon 54 in the presence of the Alu element insertion compared to only 4.7% of transcripts expressed by the wild-type minigene. These results strongly suggest ATM:c.8010+30_8010+31insAluYa5 as the main driver of ATM exon 54 skipping. Since this exon loss is predicted to cause a frameshift and a premature stop codon, mutant transcripts are unlikely to translate into functional proteins. Based on its estimated frequency of up to 0.05% in control populations, we propose to consider ATM:c.8010+30_8010+31insAluYa5 in suspected HBOC patients and to clarify its role in carcinogenesis through future epidemiological and functional analyses. Generally, the implementation of MEI detection tools in diagnostic sequencing pipelines could increase the diagnostic yield, as MEIs are likely underestimated contributors to genetic diseases.

Anti-HERV-K Drugs and Vaccines, Possible Therapies against Tumors

The footprint of human endogenous retroviruses (HERV), specifically HERV-K, has been found in malignancies, such as melanoma, teratocarcinoma, osteosarcoma, breast cancer, lymphoma, and ovary and prostate cancers. HERV-K is characterized as the most biologically active HERV due to possession of open reading frames (ORF) for all Gag, Pol, and Env genes, which enables it to be more infective and obstructive towards specific cell lines and other exogenous viruses, respectively. Some factors might contribute to carcinogenicity and at least one of them has been recognized in various tumors, including overexpression/methylation of long interspersed nuclear element 1 (LINE-1), HERV-K Gag, and Env genes themselves plus their transcripts and protein products, and HERV-K reverse transcriptase (RT). Therapies effective for HERV-K-associated tumors mostly target invasive autoimmune responses or growth of tumors through suppression of HERV-K Gag or Env protein and RT. To design new therapeutic options, more studies are needed to better understand whether HERV-K and its products (Gag/Env transcripts and HERV-K proteins/RT) are the initiators of tumor formation or just the disorder’s developers. Accordingly, this review aims to present evidence that highlights the association between HERV-K and tumorigenicity and introduces some of the available or potential therapies against HERV-K-induced tumors.

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

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

Strand asymmetries across genomic processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

LINE-associated cryptic splicing induces dsRNA-mediated interferon response and tumor immunity

RNA splicing plays a critical role in post-transcriptional gene regulation. Exponential expansion of intron length poses a challenge for accurate splicing. Little is known about how cells prevent inadvertent and often deleterious expression of intronic elements due to cryptic splicing. In this study, we identify hnRNPM as an essential RNA binding protein that suppresses cryptic splicing through binding to deep introns, preserving transcriptome integrity. Long interspersed nuclear elements (LINEs) harbor large amounts of pseudo splice sites in introns. hnRNPM preferentially binds at intronic LINEs and represses LINE-containing pseudo splice site usage for cryptic splicing. Remarkably, a subgroup of the cryptic exons can form long dsRNAs through base-pairing of inverted Alu transposable elements scattered in between LINEs and trigger interferon immune response, a well-known antiviral defense mechanism. Notably, these interferon-associated pathways are found to be upregulated in hnRNPM-deficient tumors, which also exhibit elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity. Targeting hnRNPM in tumors may be used to trigger an inflammatory immune response thereby boosting cancer surveillance.

Recurrent co-domestication of PIF/Harbinger transposable element proteins in insects

BACKGROUND

Transposable elements (TEs) are selfish DNA sequences capable of moving and amplifying at the expense of host cells. Despite this, an increasing number of studies have revealed that TE proteins are important contributors to the emergence of novel host proteins through molecular domestication. We previously described seven transposase-derived domesticated genes from the PIF/Harbinger DNA family of TEs in Drosophila and a co-domestication. All PIF TEs known in plants and animals distinguish themselves from other DNA transposons by the presence of two genes. We hypothesize that there should often be co-domestications of the two genes from the same TE because the transposase (gene 1) has been described to be translocated to the nucleus by the MADF protein (gene 2). To provide support for this model of new gene origination, we investigated available insect species genomes for additional evidence of PIF TE domestication events and explored the co-domestication of the MADF protein from the same TE insertion.

RESULTS

After the extensive insect species genomes exploration of hits to PIF transposases and analyses of their context and evolution, we present evidence of at least six independent PIF transposable elements proteins domestication events in insects: two co-domestications of both transposase and MADF proteins in Anopheles (Diptera), one transposase-only domestication event and one co-domestication in butterflies and moths (Lepidoptera), and two transposases-only domestication events in cockroaches (Blattodea). The predicted nuclear localization signals for many of those proteins and dicistronic transcription in some instances support the functional associations of co-domesticated transposase and MADF proteins.

CONCLUSIONS

Our results add to a co-domestication that we previously described in fruit fly genomes and support that new gene origination through domestication of a PIF transposase is frequently accompanied by the co-domestication of a cognate MADF protein in insects, potentially for regulatory functions. We propose a detailed model that predicts that PIF TE protein co-domestication should often occur from the same PIF TE insertion.

Comprehensive analysis of full-length transcripts reveals novel splicing abnormalities and oncogenic transcripts in liver cancer

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. Although long-reads sequencing technologies would enable the sequencing of full-length transcripts, the data analysis is difficult. In this study, we developed an analysis pipeline named SPLICE and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC) and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 46,663 transcripts from the protein-coding genes in the HCCs and the matched non-cancerous livers, of which 5,366 (11.5%) were novel. A comparison of expression levels identified 9,933 differentially expressed transcripts (DETs) in 4,744 genes. Interestingly, 746 genes with DETs, including the LINE1-MET transcript, were not found by a gene-level analysis. We also found that fusion transcripts of transposable elements and hepatitis B virus (HBV) were overexpressed in HCCs. In vitro experiments on DETs showed that LINE1-MET and HBV-human transposable elements promoted cell growth. Furthermore, fusion gene detection showed novel recurrent fusion events that were not detected in the short-reads. These results suggest the efficiency of full-length transcriptome studies and the importance of splicing variants in carcinogenesis.

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. In this study, we developed an analysis pipeline named SPLICE for long-read transcriptome sequencing and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC), and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 5,366 novel transcripts and 9,933 differentially expressed transcripts in 4,744 genes between HCCs and non-cancerous livers. An analysis of hepatitis B virus (HBV) transcripts showed that fusion transcripts of the HBV gene and human transposable elements were overexpressed in HBV-infected HCCs. We also identified fusion genes that were not found in the short-reads. These results suggest that long-reads sequencing technologies provide a fuller understanding of cancer transcripts and that our method contributes to the analysis of transcriptome sequences by such technologies.

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.

Somatic retrotransposition in the developing rhesus macaque brain

The retrotransposon LINE-1 (L1) is central to the recent evolutionary history of the human genome and continues to drive genetic diversity and germline pathogenesis. However, the spatiotemporal extent and biological significance of somatic L1 activity are poorly defined and are virtually unexplored in other primates. From a single L1 lineage active at the divergence of apes and Old World monkeys, successive L1 subfamilies have emerged in each descendant primate germline. As revealed by case studies, the presently active human L1 subfamily can also mobilize during embryonic and brain development in vivo. It is unknown whether nonhuman primate L1s can similarly generate somatic insertions in the brain. Here we applied approximately 40× single-cell whole-genome sequencing (scWGS), as well as retrotransposon capture sequencing (RC-seq), to 20 hippocampal neurons from two rhesus macaques (Macaca mulatta). In one animal, we detected and PCR-validated a somatic L1 insertion that generated target site duplications, carried a short 5' transduction, and was present in ∼7% of hippocampal neurons but absent from cerebellum and nonbrain tissues. The corresponding donor L1 allele was exceptionally mobile in vitro and was embedded in PRDM4, a gene expressed throughout development and in neural stem cells. Nanopore long-read methylome and RNA-seq transcriptome analyses indicated young retrotransposon subfamily activation in the early embryo, followed by repression in adult tissues. These data highlight endogenous macaque L1 retrotransposition potential, provide prototypical evidence of L1-mediated somatic mosaicism in a nonhuman primate, and allude to L1 mobility in the brain over the past 30 million years of human evolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Accuracy of multiple sequence alignment methods in the reconstruction of transposable element families

The construction of a high-quality multiple sequence alignment (MSA) from copies of a transposable element (TE) is a critical step in the characterization of a new TE family. Most studies of MSA accuracy have been conducted on protein or RNA sequence families, where structural features and strong signals of selection may assist with alignment. Less attention has been given to the quality of sequence alignments involving neutrally evolving DNA sequences such as those resulting from TE replication. Transposable element sequences are challenging to align due to their wide divergence ranges, fragmentation, and predominantly-neutral mutation patterns. To gain insight into the effects of these properties on MSA accuracy, we developed a simulator of TE sequence evolution, and used it to generate a benchmark with which we evaluated the MSA predictions produced by several popular aligners, along with Refiner, a method we developed in the context of our RepeatModeler software. We find that MAFFT and Refiner generally outperform other aligners for low to medium divergence simulated sequences, while Refiner is uniquely effective when tasked with aligning high-divergent and fragmented instances of a family.

Subfamily-specific differential contribution of individual monomers and the tether sequence to mouse L1 promoter activity

BACKGROUND

The internal promoter in L1 5'UTR is critical for autonomous L1 transcription and initiating retrotransposition. Unlike the human genome, which features one contemporarily active subfamily, four subfamilies (A_I, Gf_I and Tf_I/II) have been amplifying in the mouse genome in the last one million years. Moreover, mouse L1 5'UTRs are organized into tandem repeats called monomers, which are separated from ORF1 by a tether domain. In this study, we aim to compare promoter activities across young mouse L1 subfamilies and investigate the contribution of individual monomers and the tether sequence.

RESULTS

We observed an inverse relationship between subfamily age and the average number of monomers among evolutionarily young mouse L1 subfamilies. The youngest subgroup (A_I and Tf_I/II) on average carry 3-4 monomers in the 5'UTR. Using a single-vector dual-luciferase reporter assay, we compared promoter activities across six L1 subfamilies (A_I/II, Gf_I and Tf_I/II/III) and established their antisense promoter activities in a mouse embryonic fibroblast cell line and a mouse embryonal carcinoma cell line. Using consensus promoter sequences for three subfamilies (A_I, Gf_I and Tf_I), we dissected the differential roles of individual monomers and the tether domain in L1 promoter activity. We validated that, across multiple subfamilies, the second monomer consistently enhances the overall promoter activity. For individual promoter components, monomer 2 is consistently more active than the corresponding monomer 1 and/or the tether for each subfamily. Importantly, we revealed intricate interactions between monomer 2, monomer 1 and tether domains in a subfamily-specific manner. Furthermore, using three-monomer 5'UTRs, we established a complex nonlinear relationship between the length of the outmost monomer and the overall promoter activity.

CONCLUSIONS

The laboratory mouse is an important mammalian model system for human diseases as well as L1 biology. Our study extends previous findings and represents an important step toward a better understanding of the molecular mechanism controlling mouse L1 transcription as well as L1's impact on development and disease.

Methodologies for the De novo Discovery of Transposable Element Families

The discovery and characterization of transposable element (TE) families are crucial tasks in the process of genome annotation. Careful curation of TE libraries for each organism is necessary as each has been exposed to a unique and often complex set of TE families. De novo methods have been developed; however, a fully automated and accurate approach to the development of complete libraries remains elusive. In this review, we cover established methods and recent developments in de novo TE analysis. We also present various methodologies used to assess these tools and discuss opportunities for further advancement of the field.

Involvement of the DNA Phosphorothioation System in TorR Binding and Anaerobic TMAO Respiration in Salmonella enterica

Although the phosphorothioate (PT) modification, in which the nonbridging oxygen in the DNA sugar-phosphate backbone is replaced by sulfur, has been reported to play versatile roles in multiple cellular processes, very little data have been obtained to define the role of PT in epigenetic regulation. In this study, we report that the PT system in Salmonella enterica serovar Cerro 87 is involved in the transcriptional regulation of the torCAD operon encoding the trimethylamine N-oxide (TMAO) respiration machinery that enables the use of TMAO as a terminal electron acceptor for respiration when oxygen is not available. In vitro, PT enhanced the binding of the transcriptional activator of the torCAD operon, namely, TorR, to its DNA substrate (tor boxes). However, in vivo, the PT modification protein complex DndCDE downregulated torCAD transcription through competing with the binding of TorR to the tor boxes. The altered expression of torCAD caused by PT modification proteins affected cell growth that relied on TMAO respiration. To our knowledge, this is the first report supporting that PT proteins participate in transcriptional regulation, showing a new function of PT systems.

A hidden layer of structural variation in transposable elements reveals potential genetic modifiers in human disease-risk loci

Genome-wide association studies (GWAS) have been highly informative in discovering disease-associated loci but are not designed to capture all structural variations in the human genome. Using long-read sequencing data, we discovered widespread structural variation within SINE-VNTR-Alu (SVA) elements, a class of great ape-specific transposable elements with gene-regulatory roles, which represents a major source of structural variability in the human population. We highlight the presence of structurally variable SVAs (SV-SVAs) in neurological disease-associated loci, and we further associate SV-SVAs to disease-associated SNPs and differential gene expression using luciferase assays and expression quantitative trait loci data. Finally, we genetically deleted SV-SVAs in the BIN1 and CD2AP Alzheimer's disease-associated risk loci and in the BCKDK Parkinson's disease-associated risk locus and assessed multiple aspects of their gene-regulatory influence in a human neuronal context. Together, this study reveals a novel layer of genetic variation in transposable elements that may contribute to identification of the structural variants that are the actual drivers of disease associations of GWAS loci.

Characterizing mobile element insertions in 5675 genomes

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

LncRNA Biomarkers of Inflammation and Cancer

Long noncoding RNAs (lncRNAs) are promising candidates as biomarkers of inflammation and cancer. LncRNAs have several properties that make them well-suited as molecular markers of disease: (1) many lncRNAs are expressed in a tissue-specific manner, (2) distinct lncRNAs are upregulated based on different inflammatory or oncogenic stimuli, (3) lncRNAs released from cells are packaged and protected in extracellular vesicles, and (4) circulating lncRNAs in the blood are detectable using various RNA sequencing approaches. Here we focus on the potential for lncRNA biomarkers to detect inflammation and cancer, highlighting key biological, technological, and analytical considerations that will help advance the development of lncRNA-based liquid biopsies.

Transposable Elements in the Genome of the Lichen-Forming Fungus Umbilicaria pustulata and Their Distribution in Different Climate Zones along Elevation

Transposable elements (TEs) are an important source of genome plasticity across the tree of life. Drift and natural selection are important forces shaping TE distribution and accumulation. Fungi, with their multifaceted phenotypic diversity and relatively small genome size, are ideal models to study the role of TEs in genome evolution and their impact on the host's ecological and life history traits. Here we present an account of all TEs found in a high-quality reference genome of the lichen-forming fungus Umbilicaria pustulata, a macrolichen species comprising two climatic ecotypes: Mediterranean and cold temperate. We trace the occurrence of the newly identified TEs in populations along three elevation gradients using a Pool-Seq approach to identify TE insertions of potential adaptive significance. We found that TEs cover 21.26% of the 32.9 Mbp genome, with LTR Gypsy and Copia clades being the most common TEs. We identified 28 insertions displaying consistent insertion frequency differences between the two host ecotypes across the elevation gradients. Most of the highly differentiated insertions were located near genes, indicating a putative function. This pioneering study of the content and climate niche-specific distribution of TEs in a lichen-forming fungus contributes to understanding the roles of TEs in fungal evolution.