The Immunological Conundrum of Endogenous Retroelements.

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

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon.

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

Cryo-EM Structures of a Group II Intron Reverse Splicing into DNA.

Group II introns are a class of retroelements that invade DNA through a copy-and-paste mechanism known as retrotransposition. Their coordinated activities occur within a complex that includes a maturase protein, which promotes splicing through an unknown mechanism. The mechanism of splice site exchange within the RNA active site during catalysis also remains unclear. We determined two cryo-EM structures at 3.6-Å resolution of a group II intron reverse splicing into DNA. These structures reveal that the branch-site domain VI helix swings 90°, enabling substrate exchange during DNA integration. The maturase assists catalysis through a transient RNA-protein contact with domain VI that positions the branch-site adenosine for lariat formation during forward splicing. These findings provide the first direct evidence of the role the maturase plays during group II intron catalysis. The domain VI dynamics closely parallel spliceosomal branch-site helix movement and provide strong evidence for a retroelement origin of the spliceosome.

Breaching Self-Tolerance to Alu Duplex RNA Underlies MDA5-Mediated Inflammation.

Aberrant activation of innate immune receptors can cause a spectrum of immune disorders, such as Aicardi-Goutières syndrome (AGS). One such receptor is MDA5, a viral dsRNA sensor that induces antiviral immune response. Using a newly developed RNase-protection/RNA-seq approach, we demonstrate here that constitutive activation of MDA5 in AGS results from the loss of tolerance to cellular dsRNAs formed by Alu retroelements. While wild-type MDA5 cannot efficiently recognize Alu-dsRNAs because of its limited filament formation on imperfect duplexes, AGS variants of MDA5 display reduced sensitivity to duplex structural irregularities, assembling signaling-competent filaments on Alu-dsRNAs. Moreover, we identified an unexpected role of an RNA-rich cellular environment in suppressing aberrant MDA5 oligomerization, highlighting context dependence of self versus non-self discrimination. Overall, our work demonstrates that the increased efficiency of MDA5 in recognizing dsRNA comes at a cost of self-recognition and implicates a unique role of Alu-dsRNAs as virus-like elements that shape the primate immune system.

Competition between two HUSH complexes orchestrates the immune response to retroelement invasion.

The human silencing hub (HUSH) preserves genome integrity through the epigenetic repression of invasive genetic elements. However, despite our understanding of HUSH as an obligate complex of three subunits, only loss of MPP8 or Periphilin, but not TASOR, triggers interferon signaling following derepression of endogenous retroelements. Here, we resolve this paradox by characterizing a second HUSH complex that shares MPP8 and Periphilin but assembles around TASOR2, an uncharacterized paralog of TASOR. Whereas HUSH represses LINE-1 retroelements marked by the repressive histone modification H3K9me3, HUSH2 is recruited by the transcription factor IRF2 to repress interferon-stimulated genes. Mechanistically, HUSH-mediated retroelement silencing sequesters the limited pool of the shared subunits MPP8 and Periphilin, preventing TASOR2 from forming HUSH2 complexes and hence relieving the HUSH2-mediated repression of interferon-stimulated genes. Thus, competition between two HUSH complexes intertwines retroelement silencing with the induction of an immune response, coupling epigenetic and immune aspects of genome defense.

Contribution and therapeutic implications of retroelement insertions in ataxia telangiectasia.

Certain classes of genetic variation still escape detection in clinical sequencing analysis. One such class is retroelement insertion, which has been reported as a cause of Mendelian diseases and may offer unique therapeutic implications. Here, we conducted retroelement profiling on whole-genome sequencing data from a cohort of 237 individuals with ataxia telangiectasia (A-T). We found 15 individuals carrying retroelement insertions in ATM, all but one of which integrated in noncoding regions. Systematic functional characterization via RNA sequencing, RT-PCR, and/or minigene splicing assays showed that 12 out of 14 intronic insertions led or contributed to ATM loss of function by exon skipping or activating cryptic splice sites. We also present proof-of-concept antisense oligonucleotides that suppress cryptic exonization caused by a deep intronic retroelement insertion. These results provide an initial systematic estimate of the contribution of retroelements to the genetic architecture of recessive Mendelian disorders as ∼2.1%-5.5%. Our study highlights the importance of retroelement insertions as causal variants and therapeutic targets in genetic diseases.

Z-RNA biology: a central role in the innate immune response?

Z-RNA is a higher-energy, left-handed conformation of RNA, whose function has remained elusive. A growing body of work alludes to regulatory roles for Z-RNA in the immune response. Here, we review how Z-RNA features present in cellular RNAs-especially containing retroelements-could be recognized by a family of winged helix proteins, with an impact on host defense. We also discuss how mutations to specific Z-contacting amino acids disrupt their ability to stabilize Z-RNA, resulting in functional losses. We end by highlighting knowledge gaps in the field, which, if addressed, would significantly advance this active area of research.

Marsupial genome analysis suggests that satellite DNA formation from walb endogenous retrovirus is an event specific to the red-necked wallaby.

We recently identified walbRep, a satellite DNA residing in the genome of the red-necked wallaby Notamacropus rufogriseus. It originates from the walb endogenous retrovirus and is organized in a manner in which the provirus structure is retained. The walbRep repeat units feature an average pairwise nucleotide identity as high as 99.5%, raising the possibility of a recent origin. The tammar wallaby N. eugenii is a species estimated to have diverged from the red-necked wallaby 2-3 million years ago. In PCR analyses of these two and other related species, walbRep-specific fragment amplification was observed only in the red-necked wallaby. Sequence database searches for the tammar wallaby resulted in sequence alignment lists that were sufficiently powerful to exclude the possibility of walbRep existence. These results suggested that the walbRep formation occurred in the red-necked wallaby lineage after its divergence from the tammar wallaby lineage, thus in a time span of maximum 3 million years.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3' ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

Separable structural requirements for cDNA synthesis, nontemplated extension, and template jumping by a non-LTR retroelement reverse transcriptase.

Broad evolutionary expansion of polymerase families has enabled specialization of their activities for distinct cellular roles. In addition to template-complementary synthesis, many polymerases extend their duplex products by nontemplated nucleotide addition (NTA). This activity is exploited for laboratory strategies of cloning and sequencing nucleic acids and could have important biological function, although the latter has been challenging to test without separation-of-function mutations. Several retroelement and retroviral reverse transcriptases (RTs) support NTA and also template jumping, by which the RT performs continuous complementary DNA (cDNA) synthesis using physically separate templates. Previous studies that aimed to dissect the relationship between NTA and template jumping leave open questions about structural requirements for each activity and their interdependence. Here, we characterize the structural requirements for cDNA synthesis, NTA, template jumping, and the unique terminal transferase activity of Bombyx mori R2 non-long terminal repeat retroelement RT. With sequence alignments and structure modeling to guide mutagenesis, we generated enzyme variants across motifs generally conserved or specific to RT subgroups. Enzyme variants had diverse NTA profiles not correlated with other changes in cDNA synthesis activity or template jumping. Using these enzyme variants and panels of activity assay conditions, we show that template jumping requires NTA. However, template jumping by NTA-deficient enzymes can be rescued using primer duplex with a specific length of 3' overhang. Our findings clarify the relationship between NTA and template jumping as well as additional activities of non-long terminal repeat RTs, with implications for the specialization of RT biological functions and laboratory applications.

Diversity at single nucleotide to pangenome scales among sulfur cycling bacteria in salt marshes.

IMPORTANCE: Salt marshes are known for their significant carbon storage capacity, and sulfur cycling is closely linked with the ecosystem-scale carbon cycling in these ecosystems. Sulfate reducers are key for the decomposition of organic matter, and sulfur oxidizers remove toxic sulfide, supporting the productivity of marsh plants. To date, the complexity of coastal environments, heterogeneity of the rhizosphere, high microbial diversity, and uncultured majority hindered our understanding of the genomic diversity of sulfur-cycling microbes in salt marshes. Here, we use comparative genomics to overcome these challenges and provide an in-depth characterization of sulfur-cycling microbial diversity in salt marshes. We characterize communities across distinct sites and plant species and uncover extensive genomic diversity at the taxon level and specific genomic features present in MAGs affiliated with uncultivated sulfur-cycling lineages. Our work provides insights into the partnerships in salt marshes and a roadmap for multiscale analyses of diversity in complex biological systems.

Tentaclins-A Novel Family of Phage Receptor-Binding Proteins That Can Be Hypermutated by DGR Systems.

Diversity-generating retroelements (DGRs) are prokaryotic systems providing rapid modification and adaptation of target proteins. In phages, the main targets of DGRs are receptor-binding proteins that are usually parts of tail structures and the variability of such host-recognizing structures enables phage adaptation to changes on the bacterial host surface. Sometimes, more than one target gene containing a hypermutated variable repeat (VR) can be found in phage DGRs. The role of mutagenesis of two functionally different genes is unclear. In this study, several phage genomes that contain DGRs with two target genes were found in the gut virome of healthy volunteers. Bioinformatics analysis of these genes indicated that they encode proteins with different topology; however, both proteins contain the C-type lectin (C-lec) domain with a hypermutated beta-hairpin on its surface. One of the target proteins belongs to a new family of proteins with a specific topology: N-terminal C-lec domain followed by one or more immunoglobulin domains. Proteins from the new family were named tentaclins after TENTACLe + proteIN. The genes encoding such proteins were found in the genomes of prophages and phages from the gut metagenomes. We hypothesized that tentaclins are involved in binding either to bacterial receptors or intestinal/immune cells.

Conserved microRNAs and Flipons Shape Gene Expression during Development by Altering Promoter Conformations.

The classical view of gene regulation draws from prokaryotic models, where responses to environmental changes involve operons regulated by sequence-specific protein interactions with DNA, although it is now known that operons are also modulated by small RNAs. In eukaryotes, pathways based on microRNAs (miR) regulate the readout of genomic information from transcripts, while alternative nucleic acid structures encoded by flipons influence the readout of genetic programs from DNA. Here, we provide evidence that miR- and flipon-based mechanisms are deeply connected. We analyze the connection between flipon conformation and the 211 highly conserved human miR that are shared with other placental and other bilateral species. The direct interaction between conserved miR (c-miR) and flipons is supported by sequence alignments and the engagement of argonaute proteins by experimentally validated flipons as well as their enrichment in promoters of coding transcripts important in multicellular development, cell surface glycosylation and glutamatergic synapse specification with significant enrichments at false discovery rates as low as 10-116. We also identify a second subset of c-miR that targets flipons essential for retrotransposon replication, exploiting that vulnerability to limit their spread. We propose that miR can act in a combinatorial manner to regulate the readout of genetic information by specifying when and where flipons form non-B DNA (NoB) conformations, providing the interactions of the conserved hsa-miR-324-3p with RELA and the conserved hsa-miR-744 with ARHGAP5 genes as examples.

Gene-tree misrooting drives conflicts in phylogenomic coalescent analyses of palaeognath birds.

Phylogenomic analyses of ancient rapid radiations can produce conflicting results that are driven by differential sampling of taxa and characters as well as the limitations of alternative analytical methods. We re-examine basal relationships of palaeognath birds (ratites and tinamous) using recently published datasets of nucleotide characters from 20,850 loci as well as 4301 retroelement insertions. The original studies attributed conflicting resolutions of rheas in their inferred coalescent and concatenation trees to concatenation failing in the anomaly zone. By contrast, we find that the coalescent-based resolution of rheas is premised upon extensive gene-tree estimation errors. Furthermore, retroelement insertions contain much more conflict than originally reported and multiple insertion loci support the basal position of rheas found in concatenation trees, while none were reported in the original publication. We demonstrate how even remarkable congruence in phylogenomic studies may be driven by long-branch misplacement of a divergent outgroup, highly incongruent gene trees, differential taxon sampling that can result in gene-tree misrooting errors that bias species-tree inference, and gross homology errors. What was previously interpreted as broad, robustly supported corroboration for a single resolution in coalescent analyses may instead indicate a common bias that taints phylogenomic results across multiple genome-scale datasets. The updated retroelement dataset now supports a species tree with branch lengths that suggest an ancient anomaly zone, and both concatenation and coalescent analyses of the huge nucleotide datasets fail to yield coherent, reliable results in this challenging phylogenetic context.

Evolution and dosage compensation of nucleolar organizing regions (NORs) mediated by mobile elements in turtles with female (ZZ/ZW) but not with male (XX/XY) heterogamety.

Understanding the evolution and regulation of nucleolar organizing regions (NORs) is important to elucidate genome structure and function. This is because ribosomal gene (rDNA) copy number and activity mediate protein biosynthesis, stress response, ageing, disease, dosage compensation and genome stability. Here, we found contrasting dosage compensation of sex-linked NORs in turtles with male and female heterogamety. Most taxa examined exhibit homomorphic rRNA gene clusters in a single autosome pair (determined by 28S rDNA fluorescence in situ hybridization), whereas NORs are sex-linked in Apalone spinifera, Pelodiscus sinensis and Staurotypus triporcatus. Full-dosage compensation upregulates the male X-NOR (determined via silver staining-AgNOR) in Staurotypus (who lacks Y-NOR) compared with female X-AgNORs. In softshell Apalone and Pelodiscus, who share homologous ZZ/ZW micro-chromosomes, their enlarged W-NOR is partially active (due to 28S rDNA invasion by R2 retroelements), whereas their smaller Z-NOR is silent in females but active in both male-Zs (presumably because the W-NOR meets cellular demands and excessive NOR activity is costly). We hypothesize that R2 disruption favoured W enlargement to add intact 28S-units, perhaps facilitated by reduced recombination during sex chromosome evolution. The molecular basis of the potentially adaptive female Z-silencing is likely intricate and perhaps epigenetic, as non-ribosomal Z genes are active in Apalone females. Yet, Emydura maquarii exhibit identical heteromorphism in their autosomal NOR (R2 invaded 28S-units and the small-autosome NOR is silent), suggesting that the softshell turtle pattern can evolve independent of sex chromosome evolution. Our study illuminates the complex sex chromosome evolution and dosage compensation of non-model systems that challenges classic paradigms.

Full-length LTR retroelements in Capsicum annuum revealed a few species-specific family bursts with insertional preferences.

Capsicum annuum is a species that has undergone an expansion of the size of its genome caused mainly by the amplification of repetitive DNA sequences, including mobile genetic elements. Based on information obtained from sequencing the genome of pepper, the estimated fraction of retroelements is approximately 81%, and previous results revealed an important contribution of lineages derived from Gypsy superfamily. However, the dynamics of the retroelements in the C. annuum genome is poorly understood. In this way, the present work seeks to investigate the phylogenetic diversity and genomic abundance of the families of autonomous (complete and intact) LTR retroelements from C. annuum and inspect their distribution along its chromosomes. In total, we identified 1151 structurally full-length retroelements (340 Copia; 811 Gypsy) grouped in 124 phylogenetic families in the base of their retrotranscriptase. All the evolutive lineages of LTR retroelements identified in plants were present in pepper; however, three of them comprise 83% of the entire LTR retroelements population, the lineages Athila, Del/Tekay, and Ale/Retrofit. From them, only three families represent 70.8% of the total number of the identified retroelements. A massive family-specific wave of amplification of two of them occurred in the last 0.5 Mya (GypsyCa_16; CopiaCa_01), whereas the third is more ancient and occurred 3.0 Mya (GypsyCa_13). Fluorescent in situ hybridization performed with family and lineage-specific probes revealed contrasting patterns of chromosomal affinity. Our results provide a database of the populations LTR retroelements specific to C. annuum genome. The most abundant families were analyzed according to chromosome insertional preferences, suppling useful tools to the design of retroelement-based markers specific to the species.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi.

BACKGROUND: The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction. RESULTS: At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss. CONCLUSIONS: Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Whole-Genome Analyses Resolve the Phylogeny of Flightless Birds (Palaeognathae) in the Presence of an Empirical Anomaly Zone.

Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totaling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4274 CR1 retroelements identified from multispecies whole-genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other nonostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations and topology tests indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.

Surveillance of Retroelement Expression and Nucleic-Acid Immunity by Histone Methyltransferase SETDB1.

In human cancers, histone methyltransferase SETDB1 (SET domain, bifurcated 1) is frequently overexpressed but its significance in carcinogenesis remains elusive. A recent study shows that SETDB1 downregulation induces de-repression of retroelements and innate immunity in cancer cells. The possibility of SETDB1 functioning as a surveillant of retroelement expression is discussed in this study: the cytoplasmic presence of retroelement-derived nucleic acids (RdNAs) drives SETDB1 into the nucleus by the RNA-interference route, rendering the corresponding retroelement transcriptionally inert. These RdNAs could, therefore, be signals of genome instability sent out for SETDB1 present in the cytoplasm to maintain genome integrity.

Characterization and Evolution of Germ1, an Element that Undergoes Diminution in Lampreys (Cyclostomata: Petromyzontidae).

The sea lamprey (Petromyzon marinus) undergoes substantial genomic alterations during embryogenesis in which specific sequences are deleted from the genome of somatic cells yet retained in cells of the germ line. One element that undergoes diminution in P. marinus is Germ1, which consists of a somatically rare (SR) region and a fragment of 28S rDNA. Although the SR-region has been used as a marker for genomic alterations in lampreys, the evolutionary significance of its diminution is unknown. We examined the Germ1 element in five additional species of lamprey to better understand its evolutionary significance. Each representative species contained sequences similar enough to the Germ1 element of P. marinus to be detected via PCR and Southern hybridizations, although the SR-regions of Lampetra aepyptera and Lethenteron appendix are quite divergent from the homologous sequences of Petromyzon and three species of Ichthyomyzon. Lamprey Germ1 sequences have a number of features characteristic of the R2 retrotransposon, a mobile element that specifically targets 28S rDNA. Phylogenetic analyses of the SR-regions revealed patterns generally consistent with relationships among the species included in our study, although the 28S-fragments of each species/genus were most closely related to its own functional rDNA, suggesting that the two components of Germ1 were assembled independently in each lineage. Southern hybridizations showed evidence of genomic alterations involving Germ1 in each species. Our results suggest that Germ1 is a R2 retroelement that occurs in the genome of P. marinus and other petromyzontid lampreys, and that its diminution is incidental to the reduction in rDNA copies during embryogenesis.