RUNX1 is required in granulocyte-monocyte progenitors to attenuate inflammatory cytokine production by neutrophils.

The transcription factor RUNX1 is mutated in familial platelet disorder with associated myeloid malignancy (FPDMM) and in sporadic myelodysplastic syndrome and leukemia. RUNX1 was shown to regulate inflammation in multiple cell types. Here we show that RUNX1 is required in granulocyte-monocyte progenitors (GMPs) to epigenetically repress two inflammatory signaling pathways in neutrophils: Toll-like receptor 4 (TLR4) and type I interferon (IFN) signaling. RUNX1 loss in GMPs augments neutrophils' inflammatory response to the TLR4 ligand lipopolysaccharide through increased expression of the TLR4 coreceptor CD14. RUNX1 binds Cd14 and other genes encoding proteins in the TLR4 and type I IFN signaling pathways whose chromatin accessibility increases when RUNX1 is deleted. Transcription factor footprints for the effectors of type I IFN signaling-the signal transducer and activator of transcription (STAT1::STAT2) and interferon regulatory factors (IRFs)-were enriched in chromatin that gained accessibility in both GMPs and neutrophils when RUNX1 was lost. STAT1::STAT2 and IRF motifs were also enriched in the chromatin of retrotransposons that were derepressed in RUNX1-deficient GMPs and neutrophils. We conclude that a major direct effect of RUNX1 loss in GMPs is the derepression of type I IFN and TLR4 signaling, resulting in a state of fixed maladaptive innate immunity.

The Arms Race Between KRAB-Zinc Finger Proteins and Endogenous Retroelements and Its Impact on Mammals.

Nearly half of the human genome consists of endogenous retroelements (EREs) and their genetic remnants, a small fraction of which carry the potential to propagate in the host genome, posing a threat to genome integrity and cell/organismal survival. The largest family of transcription factors in tetrapods, the Krüppel-associated box domain zinc finger proteins (KRAB-ZFPs), binds to specific EREs and represses their transcription. Since their first appearance over 400 million years ago, KRAB-ZFPs have undergone dramatic expansion and diversification in mammals, correlating with the invasions of new EREs. In this article we review our current understanding of the structure, function, and evolution of KRAB-ZFPs and discuss growing evidence that the arms race between KRAB-ZFPs and the EREs they target is a major driving force for the evolution of new traits in mammals, often accompanied by domestication of EREs themselves.

Accelerated Evolution by Diversity-Generating Retroelements.

Diversity-generating retroelements (DGRs) create vast amounts of targeted, functional diversity by facilitating the rapid evolution of ligand-binding protein domains. Thousands of DGRs have been identified in bacteria, archaea, and their respective viruses. They are broadly distributed throughout the microbial world, with enrichment observed in certain taxa and environments. The diversification machinery works through a novel mechanism termed mutagenic retrohoming, whereby nucleotide sequence information is copied from an invariant DNA template repeat (TR) into an RNA intermediate, selectively mutagenized at TR adenines during cDNA synthesis by a DGR-encoded reverse transcriptase, and transferred to a variable repeat (VR) region within a variable-protein gene (54). This unidirectional flow of information leaves TR-DNA sequences unmodified, allowing for repeated rounds of mutagenic retrohoming to optimize variable-protein function. DGR target genes are often modular and can encode one or more of a wide variety of discrete functional domains appended to a diversifiable ligand-binding motif. Bacterial variable proteins often localize to cellsurfaces, although a subset appear to be cytoplasmic, while phage-encoded DGRs commonly diversify tail fiber-associated receptor-binding proteins. Here, we provide a comprehensive review of the mechanism and consequences of accelerated protein evolution by these unique and beneficial genetic elements.

Targeted hypermutation of putative antigen sensors in multicellular bacteria.

Diversity-generating retroelements (DGRs) are used by bacteria, archaea, and viruses as a targeted mutagenesis tool. Through error-prone reverse transcription, DGRs introduce random mutations at specific genomic loci, enabling rapid evolution of these targeted genes. However, the function and benefits of DGR-diversified proteins in cellular hosts remain elusive. We find that 82% of DGRs from one of the major monophyletic lineages of DGR reverse transcriptases are encoded by multicellular bacteria, which often have two or more DGR loci in their genomes. Using the multicellular purple sulfur bacterium Thiohalocapsa sp. PB-PSB1 as an example, we characterized nine distinct DGR loci capable of generating 10282 different combinations of target proteins. With environmental metagenomes from individual Thiohalocapsa aggregates, we show that most of PB-PSB1's DGR target genes are diversified across its biogeographic range, with spatial heterogeneity in the diversity of each locus. In Thiohalocapsa PB-PSB1 and other bacteria hosting this lineage of cellular DGRs, the diversified target genes are associated with NACHT-domain anti-phage defenses and putative ternary conflict systems previously shown to be enriched in multicellular bacteria. We propose that these DGR-diversified targets act as antigen sensors that confer a form of adaptive immunity to their multicellular consortia, though this remains to be experimentally tested. These findings could have implications for understanding the evolution of multicellularity, as the NACHT-domain anti-phage systems and ternary systems share both domain homology and conceptual similarities with the innate immune and programmed cell death pathways of plants and metazoans.

2-5A-Mediated decay (2-5AMD): from antiviral defense to control of host RNA.

Mammalian cells are exquisitely sensitive to the presence of double-stranded RNA (dsRNA), a molecule that they interpret as a signal of viral presence requiring immediate attention. Upon sensing dsRNA cells activate the innate immune response, which involves transcriptional mechanisms driving inflammation and secretion of interferons (IFNs) and interferon-stimulated genes (ISGs), as well as synthesis of RNA-like signaling molecules comprised of three or more 2'-5'-linked adenylates (2-5As). 2-5As were discovered some forty years ago and described as IFN-induced inhibitors of protein synthesis. The efforts of many laboratories, aimed at elucidating the molecular mechanism and function of these mysterious RNA-like signaling oligonucleotides, revealed that 2-5A is a specific ligand for the kinase-family endonuclease RNase L. RNase L decays single-stranded RNA (ssRNA) from viruses and mRNAs (as well as other RNAs) from hosts in a process we proposed to call 2-5A-mediated decay (2-5AMD). During recent years it has become increasingly recognized that 2-5AMD is more than a blunt tool of viral RNA destruction, but a pathway deeply integrated into sensing and regulation of endogenous RNAs. Here we present an overview of recently emerged roles of 2-5AMD in host RNA regulation.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents.

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900 bp full-length element of mys-related sequences (mysRS), an 8000 bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800 bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

The nature and genomic landscape of repetitive DNA classes in Chrysanthemum nankingense shows recent genomic changes.

BACKGROUND AND AIMS: Tandemly repeated DNA and transposable elements represent most of the DNA in higher plant genomes. High-throughput sequencing allows a survey of the DNA in a genome, but whole-genome assembly can miss a substantial fraction of highly repeated sequence motifs. Chrysanthemum nankingense (2n = 2x = 18; genome size = 3.07 Gb; Asteraceae), a diploid reference for the many auto- and allopolyploids in the genus, was considered as an ancestral species and serves as an ornamental plant and high-value food. We aimed to characterize the major repetitive DNA motifs, understand their structure and identify key features that are shaped by genome and sequence evolution. METHODS: Graph-based clustering with RepeatExplorer was used to identify and classify repetitive motifs in 2.14 millions of 250-bp paired-end Illumina reads from total genomic DNA of C. nankingense. Independently, the frequency of all canonical motifs k-bases long was counted in the raw read data and abundant k-mers (16, 21, 32, 64 and 128) were extracted and assembled to generate longer contigs for repetitive motif identification. For comparison, long terminal repeat retrotransposons were checked in the published C. nankingense reference genome. Fluorescent in situ hybridization was performed to show the chromosomal distribution of the main types of repetitive motifs. KEY RESULTS: Apart from rDNA (0.86 % of the total genome), a few microsatellites (0.16 %), and telomeric sequences, no highly abundant tandem repeats were identified. There were many transposable elements: 40 % of the genome had sequences with recognizable domains related to transposable elements. Long terminal repeat retrotransposons showed widespread distribution over chromosomes, although different sequence families had characteristic features such as abundance at or exclusion from centromeric or subtelomeric regions. Another group of very abundant repetitive motifs, including those most identified as low-complexity sequences (9.07 %) in the genome, showed no similarity to known sequence motifs or tandemly repeated elements. CONCLUSIONS: The Chrysanthemum genome has an unusual structure with a very low proportion of tandemly repeated sequences (~1.02 %) in the genome, and a high proportion of low-complexity sequences, most likely degenerated remains of transposable elements. Identifying the presence, nature and genomic organization of major genome fractions enables inference of the evolutionary history of sequences, including degeneration and loss, critical to understanding biodiversity and diversification processes in the genomes of diploid and polyploid Chrysanthemum, Asteraceae and plants more widely.

Diversity-Generating Retroelements in Prokaryotic Immunity.

Adaptive immunity systems found in different organisms fall into two major types. Prokaryotes possess CRISPR-Cas systems that recognize former invaders using memorized (captured) pieces of their DNA as pathogen signatures. Mammals possess a vast repertoire of antibodies and T-cell receptor variants generated in advance. In this second type of adaptive immunity, a pathogen presentation to the immune system specifically activates the cells that express matching antibodies or receptors. These cells proliferate to fight the infection and form the immune memory. The principle of preemptive production of diverse defense proteins for future use can hypothetically take place in microbes too. We propose a hypothesis that prokaryotes employ diversity-generating retroelements to prepare defense proteins against yet-unknown invaders. In this study, we test this hypothesis with the methods of bioinformatics and identify several candidate defense systems based on diversity-generating retroelements.

[The Involvement of the Gagr Gene, a Domesticated gag Retrovirus Gene, in the Stress-Response Pathway in Different Drosophila Species].

The Gagr gene is a domesticated gag retroelement gene in Drosophila melanogaster, whose function is associated with the stress response. The protein products of the Gagr gene and its homologues in different Drosophila species have a highly conserved structure; however, they demonstrate variability in the promoter region of the gene, which is apparently associated with the gradual acquisition of a new function and involvement in new signaling pathways. In this work, we studied the effect of oxidative stress induced by ammonium persulfate on the survival of various species of the genus Drosophila (D. melanogaster, D. mauritiana, D. simulans, D. yakuba, D. teissieri, and D. pseudoobscura), analyzed the correlation between the structure of promoter regions and stress-induced changes in the expression of the Gagr gene and its homologues in different Drosophila species, and compared the stress-induced changes in the expression of oxidative stress markers: the Jak-STAT signaling pathway activator gene upd3, Jak-STAT pathway effector vir-1, and IMD signaling pathway target Rel. It was found that in D. simulans and D. mauritiana sensitivity to ammonium persulfate is significantly increased, which correlates with a reduced level of transcription of vir-1 gene orthologues. The latter is due to a decrease in the number of binding sites for the transcription factor STAT92E, a component of the Jak-STAT signaling pathway, in the vir-1 promoter region. Consistent changes in the expression of the Gagr, upd3, and vir-1 genes are observed in all species of the melanogaster subgroup, except for D. pseudoobscura, which indicates an increase in the role of Gagr in the regulation of stress response pathways during the phylogenesis of the genus Drosophila.

Comparative repeatome analysis reveals new evidence on genome evolution in wild diploid Arachis (Fabaceae) species.

MAIN CONCLUSION: Opposing changes in the abundance of satellite DNA and long terminal repeat (LTR) retroelements are the main contributors to the variation in genome size and heterochromatin amount in Arachis diploids. The South American genus Arachis (Fabaceae) comprises 83 species organized in nine taxonomic sections. Among them, section Arachis is characterized by species with a wide genome and karyotype diversity. Such diversity is determined mainly by the amount and composition of repetitive DNA. Here we performed computational analysis on low coverage genome sequencing to infer the dynamics of changes in major repeat families that led to the differentiation of genomes in diploid species (x = 10) of genus Arachis, focusing on section Arachis. Estimated repeat content ranged from 62.50 to 71.68% of the genomes. Species with different genome composition tended to have different landscapes of repeated sequences. Athila family retrotransposons were the most abundant and variable lineage among Arachis repeatomes, with peaks of transpositional activity inferred at different times in the evolution of the species. Satellite DNAs (satDNAs) were less abundant, but differentially represented among species. High rates of evolution of an AT-rich superfamily of satDNAs led to the differential accumulation of heterochromatin in Arachis genomes. The relationship between genome size variation and the repetitive content is complex. However, largest genomes presented a higher accumulation of LTR elements and lower contents of satDNAs. In contrast, species with lowest genome sizes tended to accumulate satDNAs in detriment of LTR elements. Phylogenetic analysis based on repetitive DNA supported the genome arrangement of section Arachis. Altogether, our results provide the most comprehensive picture on the repeatome dynamics that led to the genome differentiation of Arachis species.

The Absence of Retroelement Activity Is Characteristic for Childhood Acute Leukemias and Adult Acute Lymphoblastic Leukemia.

Retroelements (RE) have been proposed as important players in cancerogenesis. Different cancer types are characterized by a different level of tumor-specific RE insertions. In previous studies, small cohorts of hematological malignancies, such as acute myeloid leukemia, multiple myeloma, and chronic lymphocytic leukemia have been characterized by a low level of RE insertional activity. Acute lymphoblastic leukemia (ALL) in adults and childhood acute leukemias have not been studied in this context. We performed a search for new RE insertions (Alu and L1) in 44 childhood ALL, 14 childhood acute myeloid leukemia, and 14 adult ALL samples using a highly sensitive NGS-based approach. First, we evaluated the method sensitivity revealing the 1% detection threshold for the proportion of cells with specific RE insertion. Following this result, we did not identify new tumor-specific RE insertions in the tested cohort of acute leukemia samples at the established level of sensitivity. Additionally, we analyzed the transcription levels of active L1 copies and found them increased. Thus, the increased transcription of active L1 copies is not sufficient for overt elevation of L1 retrotranspositional activity in leukemia.

In-Plant Persistence and Systemic Transport of Nicotiana benthamiana Retrozyme RNA.

Retrozymes are nonautonomous retrotransposons with hammerhead ribozymes in their long terminal repeats (LTRs). Retrozyme transcripts can be self-cleaved by the LTR ribozyme, circularized, and can undergo RNA-to-RNA replication. Here, we demonstrate that the Nicotiana benthamiana genome contains hundreds of retrozyme loci, of which nine represent full-length retrozymes. The LTR contains a promoter directing retrozyme transcription. Although retrozyme RNA is easily detected in plants, the LTR region is heavily methylated, pointing to its transcriptional silencing, which can be mediated by 24 nucleotide-long retrozyme-specific RNAs identified in N. benthamiana. A transcriptome analysis revealed that half of the retrozyme-specific RNAs in plant leaves have no exact matches to genomic retrozyme loci, containing up to 13% mismatches with the closest genomic sequences, and could arise as a result of many rounds of RNA-to-RNA replication leading to error accumulation. Using a cloned retrozyme copy, we show that retrozyme RNA is capable of replication and systemic transport in plants. The presented data suggest that retrozyme loci in the N. benthamiana genome are transcriptionally inactive, and that circular retrozyme RNA can persist in cells due to its RNA-to-RNA replication and be transported systemically, emphasizing functional and, possibly, evolutionary links of retrozymes to viroids-noncoding circular RNAs that infect plants.

Genetic Variants Associated with Long-Terminal Repeats Can Diagnostically Classify Cannabis Varieties.

Cannabis sativa (Cannabis) has recently been legalized in multiple countries globally for either its recreational or medicinal use. This, in turn, has led to a marked increase in the number of Cannabis varieties available for use in either market. However, little information currently exists on the genetic distinction between adopted varieties. Such fundamental knowledge is of considerable value and underpins the accelerated development of both a nascent pharmaceutical industry and the commercial recreational market. Therefore, in this study, we sought to assess genetic diversity across 10 Cannabis varieties by undertaking a reduced representation shotgun sequencing approach on 83 individual plants to identify variations which could be used to resolve the genetic structure of the assessed population. Such an approach also allowed for the identification of the genetic features putatively associated with the production of secondary metabolites in Cannabis. Initial analysis identified 3608 variants across the assessed population with phylogenetic analysis of this data subsequently enabling the confident grouping of each variety into distinct subpopulations. Within our dataset, the most diagnostically informative single nucleotide polymorphisms (SNPs) were determined to be associated with the long-terminal repeat (LTRs) class of retroelements, with 172 such SNPs used to fully resolve the genetic structure of the assessed population. These 172 SNPs could be used to design a targeted resequencing panel, which we propose could be used to rapidly screen different Cannabis plants to determine genetic relationships, as well as to provide a more robust, scientific classification of Cannabis varieties as the field moves into the pharmaceutical sphere.

Expression of APOBEC family members as regulators of endogenous retroelements and malignant transformation in systemic autoimmunity.

OBJECTIVE: To explore whether APOBEC family members are involved in the response to inappropriate expression of L1 retroelements in primary Sjögren's syndrome (SS) and systemic lupus erythematosus (SLE), as well as in SS related lymphomagenesis. METHODS: Minor salivary glands (MSG) and kidney biopsy (KB) specimens were obtained from 41 SS patients (10 with lymphoma) and 23 patients with SLE, respectively. PBMC and sera were also collected from 73 SLE patients. Full-length L1 transcripts, members of the APOBEC and IFN family were quantitated by real time PCR. Type I IFN activity was assessed in lupus plasma by a cell assay. RESULTS: APOBEC3A was increased in SS MSG, SLE KB and PBMC and correlated with L1. AID and APOBEC3G were particularly overexpressed in MSG tissues derived from SS lymphoma patients. CONCLUSION: These data reveal a previously unappreciated role of APOBEC family proteins in the pathogenesis of systemic autoimmunity and SS related lymphomagenesis.

Testing the retroelement invasion hypothesis for the emergence of the ancestral eukaryotic cell.

Phylogenetic evidence suggests that the invasion and proliferation of retroelements, selfish mobile genetic elements that copy and paste themselves within a host genome, was one of the early evolutionary events in the emergence of eukaryotes. Here we test the effects of this event by determining the pressures retroelements exert on simple genomes. We transferred two retroelements, human LINE-1 and the bacterial group II intron Ll.LtrB, into bacteria, and find that both are functional and detrimental to growth. We find, surprisingly, that retroelement lethality and proliferation are enhanced by the ability to perform eukaryotic-like nonhomologous end-joining (NHEJ) DNA repair. We show that the only stable evolutionary consequence in simple cells is maintenance of retroelements in low numbers, suggesting how retrotransposition rates and costs in early eukaryotes could have been constrained to allow proliferation. Our results suggest that the interplay between NHEJ and retroelements may have played a fundamental and previously unappreciated role in facilitating the proliferation of retroelements, elements of which became the ancestors of the spliceosome components in eukaryotes.

ERVmap analysis reveals genome-wide transcription of human endogenous retroviruses.

Endogenous retroviruses (ERVs) are integrated retroviral elements that make up 8% of the human genome. However, the impact of ERVs on human health and disease is not well understood. While select ERVs have been implicated in diseases, including autoimmune disease and cancer, the lack of tools to analyze genome-wide, locus-specific expression of proviral autonomous ERVs has hampered the progress in the field. Here we describe a method called ERVmap, consisting of an annotated database of 3,220 human proviral ERVs and a pipeline that allows for locus-specific genome-wide identification of proviral ERVs that are transcribed based on RNA-sequencing data, and provide examples of the utility of this tool. Using ERVmap, we revealed cell-type-specific ERV expression patterns in commonly used cell lines as well as in primary cells. We identified 124 unique ERV loci that are significantly elevated in the peripheral blood mononuclear cells of patients with systemic lupus erythematosus that represent an IFN-independent signature. Finally, we identified additional tumor-associated ERVs that correlate with cytolytic activity represented by granzyme and perforin expression in breast cancer tissue samples. The open-source code of ERVmap and the accompanied web tool are made publicly available to quantify proviral ERVs in RNA-sequencing data with ease. Use of ERVmap across a range of diseases and experimental conditions has the potential to uncover novel disease-associated antigens and effectors involved in human health that is currently missed by focusing on protein-coding sequences.

Diversification of AID/APOBEC-like deaminases in metazoa: multiplicity of clades and widespread roles in immunity.

AID/APOBEC deaminases (AADs) convert cytidine to uridine in single-stranded nucleic acids. They are involved in numerous mutagenic processes, including those underpinning vertebrate innate and adaptive immunity. Using a multipronged sequence analysis strategy, we uncover several AADs across metazoa, dictyosteliida, and algae, including multiple previously unreported vertebrate clades, and versions from urochordates, nematodes, echinoderms, arthropods, lophotrochozoans, cnidarians, and porifera. Evolutionary analysis suggests a fundamental division of AADs early in metazoan evolution into secreted deaminases (SNADs) and classical AADs, followed by diversification into several clades driven by rapid-sequence evolution, gene loss, lineage-specific expansions, and lateral transfer to various algae. Most vertebrate AADs, including AID and APOBECs1-3, diversified in the vertebrates, whereas the APOBEC4-like clade has a deeper origin in metazoa. Positional entropy analysis suggests that several AAD clades are diversifying rapidly, especially in the positions predicted to interact with the nucleic acid target motif, and with potential viral inhibitors. Further, several AADs have evolved neomorphic metal-binding inserts, especially within loops predicted to interact with the target nucleic acid. We also observe polymorphisms, driven by alternative splicing, gene loss, and possibly intergenic recombination between paralogs. We propose that biological conflicts of AADs with viruses and genomic retroelements are drivers of rapid AAD evolution, suggesting a widespread presence of mutagenesis-based immune-defense systems. Deaminases like AID represent versions "institutionalized" from the broader array of AADs pitted in such arms races for mutagenesis of self-DNA, and similar recruitment might have independently occurred elsewhere in metazoa.

A comparative study of long interspersed element-1 protein immunoreactivity in cutaneous malignancies.

BACKGROUND: Skin cancer is the most common cancer worldwide and commonly classified into malignant melanoma (MM) and Nonmelanoma skin cancers (NMSCs), which mainly include basal cell carcinoma (BCC) and squamous cell carcinoma (SCC). The extent to which Long Interspersed Element-1 (LINE-1, L1) ORF1p is expressed in cutaneous malignancies remains to be evaluated. This study aimed to assess LINE-1 ORF1p immunoreactivity in various skin cancer subtypes. METHOD: The expression level of LINE-1 ORF1p was evaluated in 95 skin cancer specimens comprising 36 (37.9%) BCC, 28 (29.5%) SCC, and 31 (32.6%) melanoma using the tissue microarray (TMA) technique. Then the association between expression of LINE-1 encoded protein and clinicopathological parameters was analyzed. RESULTS: We showed that LINE-1 ORF1p expression level was substantially higher in BCC and SCC patients compared with melanoma samples (p < 0.001). BCC cases had a higher LINE-1 histochemical score (H-score) compared with SCC cases (p = 0.004). In SCC samples, a lower level of LINE-1 ORF1p expression was associated with age younger than the mean (p = 0.041). At the same time, no significant correlation was found between LINE-1 ORF1p expression and other clinicopathological parameters (all p > 0.05). CONCLUSIONS: According to our observation, LINE-1 ORF1p immunoreactivity in various skin tumor subtypes extends previous studies of LINE-1 expression in different cancers. LINE-1ORF1p overexpression in NMSCs compared with MM can be considered with caution as a tumor-specific antigen for NMSCs.

ILS-Aware Analysis of Low-Homoplasy Retroelement Insertions: Inference of Species Trees and Introgression Using Quartets.

DNA sequence alignments have provided the majority of data for inferring phylogenetic relationships with both concatenation and coalescent methods. However, DNA sequences are susceptible to extensive homoplasy, especially for deep divergences in the Tree of Life. Retroelement insertions have emerged as a powerful alternative to sequences for deciphering evolutionary relationships because these data are nearly homoplasy-free. In addition, retroelement insertions satisfy the "no intralocus-recombination" assumption of summary coalescent methods because they are singular events and better approximate neutrality relative to DNA loci commonly sampled in phylogenomic studies. Retroelements have traditionally been analyzed with parsimony, distance, and network methods. Here, we analyze retroelement data sets for vertebrate clades (Placentalia, Laurasiatheria, Balaenopteroidea, Palaeognathae) with 2 ILS-aware methods that operate by extracting, weighting, and then assembling unrooted quartets into a species tree. The first approach constructs a species tree from retroelement bipartitions with ASTRAL, and the second method is based on split-decomposition with parsimony. We also develop a Quartet-Asymmetry test to detect hybridization using retroelements. Both ILS-aware methods recovered the same species-tree topology for each data set. The ASTRAL species trees for Laurasiatheria have consecutive short branch lengths in the anomaly zone whereas Palaeognathae is outside of this zone. For the Balaenopteroidea data set, which includes rorquals (Balaenopteridae) and gray whale (Eschrichtiidae), both ILS-aware methods resolved balaeonopterids as paraphyletic. Application of the Quartet-Asymmetry test to this data set detected 19 different quartets of species for which historical introgression may be inferred. Evidence for introgression was not detected in the other data sets.

DGR mutagenic transposition occurs via hypermutagenic reverse transcription primed by nicked template RNA.

Diversity-generating retroelements (DGRs) are molecular evolution machines that facilitate microbial adaptation to environmental changes. Hypervariation occurs via a mutagenic retrotransposition process from a template repeat (TR) to a variable repeat (VR) that results in adenine-to-random nucleotide conversions. Here we show that reverse transcription of the Bordetella phage DGR is primed by an adenine residue in TR RNA and is dependent on the DGR-encoded reverse transcriptase (bRT) and accessory variability determinant (Avd ), but is VR-independent. We also find that the catalytic center of bRT plays an essential role in site-specific cleavage of TR RNA for cDNA priming. Adenine-specific mutagenesis occurs during reverse transcription and does not involve dUTP incorporation, indicating it results from bRT-catalyzed misincorporation of standard deoxyribonucleotides. In vivo assays show that this hybrid RNA-cDNA molecule is required for mutagenic transposition, revealing a unique mechanism of DNA hypervariation for microbial adaptation.