Virus genomes are prone to extensive gene loss, gain, and exchange and share no universal genes. Therefore, in a broad-scale study of virus evolution, gene and genome network analyses can complement traditional phylogenetics. We performed an exhaustive comparative analysis of the genomes of double-stranded DNA (dsDNA) viruses by using the bipartite network approach and found a robust hierarchical modularity in the dsDNA virosphere. Bipartite networks consist of two classes of nodes, with nodes in one class, in this case genomes, being connected via nodes of the second class, in this case genes. Such a network can be partitioned into modules that combine nodes from both classes. The bipartite network of dsDNA viruses includes 19 modules that form 5 major and 3 minor supermodules. Of these modules, 11 include tailed bacteriophages, reflecting the diversity of this largest group of viruses. The module analysis quantitatively validates and refines previously proposed nontrivial evolutionary relationships. An expansive supermodule combines the large and giant viruses of the putative order "Megavirales" with diverse moderate-sized viruses and related mobile elements. All viruses in this supermodule share a distinct morphogenetic tool kit with a double jelly roll major capsid protein. Herpesviruses and tailed bacteriophages comprise another supermodule, held together by a distinct set of morphogenetic proteins centered on the HK97-like major capsid protein. Together, these two supermodules cover the great majority of currently known dsDNA viruses. We formally identify a set of 14 viral hallmark genes that comprise the hubs of the network and account for most of the intermodule connections.

Insertion sequences (IS) are the simplest and most abundant form of transposable DNA found in bacterial genomes. When present in multiple copies, it is thought that they can promote genomic plasticity and genetic exchange, thus being a major force of evolutionary change. The main processes that determine IS content in genomes are, though, a matter of debate. In this work, we take advantage of the large amount of genomic data currently available and study the abundance distributions of 33 IS families in 1811 bacterial chromosomes. This allows us to test simple models of IS dynamics and estimate their key parameters by means of a maximum likelihood approach. We evaluate the roles played by duplication, lateral gene transfer, deletion and purifying selection. We find that the observed IS abundances are compatible with a neutral scenario where IS proliferation is controlled by deletions instead of purifying selection. Even if there may be some cases driven by selection, neutral behavior dominates over large evolutionary scales. According to this view, IS and hosts tend to coexist in a dynamic equilibrium state for most of the time. Our approach also allows for a detection of recent IS expansions, and supports the hypothesis that rapid expansions constitute transient events-punctuations-during which the state of coexistence of IS and host becomes perturbated.

Ribonucleoside 5'-monophosphates (rNMPs) are the most common non-standard nucleotides found in DNA of eukaryotic cells, with over 100 million rNMPs transiently incorporated in the mammalian genome per cell cycle. Human ribonuclease (RNase) H2 is the principal enzyme able to cleave rNMPs in DNA. Whether RNase H2 may process abasic or oxidized rNMPs incorporated in DNA is unknown. The base excision repair (BER) pathway is mainly responsible for repairing oxidized and abasic sites into DNA. Here we show that human RNase H2 is unable to process an abasic rNMP (rAP site) or a ribose 8oxoG (r8oxoG) site embedded in DNA. On the contrary, we found that recombinant purified human apurinic/apyrimidinic endonuclease-1 (APE1) and APE1 from human cell extracts efficiently process an rAP site in DNA and have weak endoribonuclease and 3'-exonuclease activities on r8oxoG substrate. Using biochemical assays, our results provide evidence of a human enzyme able to recognize and process abasic and oxidized ribonucleotides embedded in DNA.

Two classes of RNase H hydrolyze RNA of RNA/DNA hybrids. In contrast to RNase H1 that requires four ribonucleotides for cleavage, RNase H2 can nick duplex DNAs containing a single ribonucleotide, suggesting different in vivo substrates. We report here the crystal structures of a type 2 RNase H in complex with substrates containing a (5')RNA-DNA(3') junction. They revealed a unique mechanism of recognition and substrate-assisted cleavage. A conserved tyrosine residue distorts the nucleic acid at the junction, allowing the substrate to function in catalysis by participating in coordination of the active site metal ion. The biochemical and structural properties of RNase H2 explain the preference of the enzyme for junction substrates and establish the structural and mechanistic differences with RNase H1. Junction recognition is important for the removal of RNA embedded in DNA and may play an important role in DNA replication and repair.

RNase H2 cleaves RNA sequences that are part of RNA/DNA hybrids or that are incorporated into DNA, thus, preventing genomic instability and the accumulation of aberrant nucleic acid, which in humans induces Aicardi-Goutières syndrome, a severe autoimmune disorder. The 3.1 Å crystal structure of human RNase H2 presented here allowed us to map the positions of all 29 mutations found in Aicardi-Goutières syndrome patients, several of which were not visible in the previously reported mouse RNase H2. We propose the possible effects of these mutations on the protein stability and function. Bacterial and eukaryotic RNases H2 differ in composition and substrate specificity. Bacterial RNases H2 are monomeric proteins and homologs of the eukaryotic RNases H2 catalytic subunit, which in addition possesses two accessory proteins. The eukaryotic RNase H2 heterotrimeric complex recognizes RNA/DNA hybrids and (5')RNA-DNA(3')/DNA junction hybrids as substrates with similar efficiency, whereas bacterial RNases H2 are highly specialized in the recognition of the (5')RNA-DNA(3') junction and very poorly cleave RNA/DNA hybrids in the presence of Mg(2+) ions. Using the crystal structure of the Thermotoga maritima RNase H2-substrate complex, we modeled the human RNase H2-substrate complex and verified the model by mutational analysis. Our model indicates that the difference in substrate preference stems from the different position of the crucial tyrosine residue involved in substrate binding and recognition.

Cellular levels of ribonucleoside triphosphates (rNTPs) are much higher than those of deoxyribonucleoside triphosphates (dNTPs), thereby influencing the frequency of incorporation of ribonucleoside monophosphates (rNMPs) by DNA polymerases (Pol) into DNA. RNase H2-initiated ribonucleotide excision repair (RER) efficiently removes single rNMPs in genomic DNA. However, processing of rNMPs by Topoisomerase 1 (Top1) in absence of RER induces mutations and genome instability. Here, we greatly increased the abundance of genomic rNMPs in Saccharomyces cerevisiae by depleting Rnr1, the major subunit of ribonucleotide reductase, which converts ribonucleotides to deoxyribonucleotides. We found that in strains that are depleted of Rnr1, RER-deficient, and harbor an rNTP-permissive replicative Pol mutant, excessive accumulation of single genomic rNMPs severely compromised growth, but this was reversed in absence of Top1. Thus, under Rnr1 depletion, limited dNTP pools slow DNA synthesis by replicative Pols and provoke the incorporation of high levels of rNMPs in genomic DNA. If a threshold of single genomic rNMPs is exceeded in absence of RER and presence of limited dNTP pools, Top1-mediated genome instability leads to severe growth defects. Finally, we provide evidence showing that accumulation of RNA/DNA hybrids in absence of RNase H1 and RNase H2 leads to cell lethality under Rnr1 depletion.

Viruses with RNA genomes dominate the eukaryotic virome, reaching enormous diversity in animals and plants. The recent advances of metaviromics prompted us to perform a detailed phylogenomic reconstruction of the evolution of the dramatically expanded global RNA virome. The only universal gene among RNA viruses is the gene encoding the RNA-dependent RNA polymerase (RdRp). We developed an iterative computational procedure that alternates the RdRp phylogenetic tree construction with refinement of the underlying multiple-sequence alignments. The resulting tree encompasses 4,617 RNA virus RdRps and consists of 5 major branches; 2 of the branches include positive-sense RNA viruses, 1 is a mix of positive-sense (+) RNA and double-stranded RNA (dsRNA) viruses, and 2 consist of dsRNA and negative-sense (-) RNA viruses, respectively. This tree topology implies that dsRNA viruses evolved from +RNA viruses on at least two independent occasions, whereas -RNA viruses evolved from dsRNA viruses. Reconstruction of RNA virus evolution using the RdRp tree as the scaffold suggests that the last common ancestors of the major branches of +RNA viruses encoded only the RdRp and a single jelly-roll capsid protein. Subsequent evolution involved independent capture of additional genes, in particular, those encoding distinct RNA helicases, enabling replication of larger RNA genomes and facilitating virus genome expression and virus-host interactions. Phylogenomic analysis reveals extensive gene module exchange among diverse viruses and horizontal virus transfer between distantly related hosts. Although the network of evolutionary relationships within the RNA virome is bound to further expand, the present results call for a thorough reevaluation of the RNA virus taxonomy.IMPORTANCE The majority of the diverse viruses infecting eukaryotes have RNA genomes, including numerous human, animal, and plant pathogens. Recent advances of metagenomics have led to the discovery of many new groups of RNA viruses in a wide range of hosts. These findings enable a far more complete reconstruction of the evolution of RNA viruses than was attainable previously. This reconstruction reveals the relationships between different Baltimore classes of viruses and indicates extensive transfer of viruses between distantly related hosts, such as plants and animals. These results call for a major revision of the existing taxonomy of RNA viruses.

The neuroinflammatory autoimmune disease Aicardi-Goutières syndrome (AGS) develops from mutations in genes encoding several nucleotide-processing proteins, including RNase H2. Defective RNase H2 may induce accumulation of self-nucleic acid species that trigger chronic type I interferon and inflammatory responses, leading to AGS pathology. We created a knock-in mouse model with an RNase H2 AGS mutation in a highly conserved residue of the catalytic subunit, Rnaseh2a(G37S/G37S) (G37S), to understand disease pathology. G37S homozygotes are perinatal lethal, in contrast to the early embryonic lethality previously reported for Rnaseh2b- or Rnaseh2c-null mice. Importantly, we found that the G37S mutation led to increased expression of interferon-stimulated genes dependent on the cGAS-STING signaling pathway. Ablation of STING in the G37S mice results in partial rescue of the perinatal lethality, with viable mice exhibiting white spotting on their ventral surface. We believe that the G37S knock-in mouse provides an excellent animal model for studying RNASEH2-associated autoimmune diseases.

Eukaryotic RNase H2 is a heterotrimeric enzyme. Here, we show that the biochemical composition and stoichiometry of the human RNase H2 complex is consistent with the properties previously deduced from genetic studies. The catalytic subunit of eukaryotic RNase H2, RNASEH2A, is well conserved and similar to the monomeric prokaryotic RNase HII. In contrast, the RNASEH2B and RNASEH2C subunits from human and Saccharomyces cerevisiae share very little homology, although they both form soluble B/C complexes that may serve as a nucleation site for the addition of RNASEH2A to form an active RNase H2, or for interactions with other proteins to support different functions. The RNASEH2B subunit has a PIP-box and confers PCNA binding to human RNase H2. Unlike Escherichia coli RNase HII, eukaryotic RNase H2 acts processively and hydrolyzes a variety of RNA/DNA hybrids with similar efficiencies, suggesting multiple cellular substrates. Moreover, of five analyzed mutations in human RNASEH2B and RNASEH2C linked to Aicardi-Goutières Syndrome (AGS), only one, R69W in the RNASEH2C protein, exhibits a significant reduction in specific activity, revealing a role for the C subunit in enzymatic activity. Near-normal activity of four AGS-related mutant enzymes was unexpected in light of their predicted impairment causing the AGS phenotype.

Understanding complex diseases will benefit the recognition of the properties of the gene networks that control biological functions. Here, we set out to model the gene network that controls T-cell activation in humans, which is critical for the development of autoimmune diseases such as Multiple Sclerosis (MS). The network was established on the basis of the quantitative expression from 104 individuals of 20 genes of the immune system, as well as on biological information from the Ingenuity database and Bayesian inference. Of the 31 links (gene interactions) identified in the network, 18 were identified in the Ingenuity database and 13 were new and we validated 7 of 8 interactions experimentally. In the MS patients network, we found an increase in the weight of gene interactions related to Th1 function and a decrease in those related to Treg and Th2 function. Indeed, we found that IFN-ss therapy induces changes in gene interactions related to T cell proliferation and adhesion, although these gene interactions were not restored to levels similar to controls. Finally, we identify JAG1 as a new therapeutic target whose differential behaviour in the MS network was not modified by immunomodulatory therapy. In vitro treatment with a Jagged1 agonist peptide modulated the T-cell activation network in PBMCs from patients with MS. Moreover, treatment of mice with experimental autoimmune encephalomyelitis with the Jagged1 agonist ameliorated the disease course, and modulated Th2, Th1 and Treg function. This study illustrates how network analysis can predict therapeutic targets for immune intervention and identified the immunomodulatory properties of Jagged1 making it a new therapeutic target for MS and other autoimmune diseases.

R-loops, three-strand structures consisting of mRNA hybridized to the complementary DNA and a single-stranded DNA loop, are formed in switch regions on the heavy-chain immunoglobulin locus. To determine if R-loops have a direct effect on any of the steps involved in isotype switching, we generated a transgenic mouse that over-expressed RNase H1, an enzyme that cleaves the RNA of RNA/DNA hybrids in B cells. R-loops in the switch μ region were depleted by 70% in ex vivo activated splenic B cells. Frequencies of isotype switching to IgG1, IgG2b, IgG2c, and IgG3 were the same as C57BL/6 control cells. However, somatic hypermutation was increased specifically on the transcribed strand from μ-γ joins, indicating that R-loops limit activation-induced (cytosine) deaminase access to the transcribed DNA strand. Our data suggest that, in the normal G+C-rich context of mammalian class switch recombination regions, R-loops are obligatory intermediates. Processing of the R-loops is needed to remove RNA allowing activation-induced (cytosine) deaminase to promote somatic hypermutation on both DNA strands to generate double-strand DNA breaks for efficient class switch recombination. One of the two cellular RNases H may assist in this process.

RNase H2 has two distinct functions: initiation of the ribonucleotide excision repair (RER) pathway by cleaving ribonucleotides (rNMPs) incorporated during DNA replication and processing the RNA portion of an R-loop formed during transcription. An RNase H2 mutant lacking RER activity but supporting R-loop removal revealed that rNMPs in DNA initiate p53-dependent DNA damage response and early embryonic arrest in mouse. However, an RNase H2 AGS-related mutant with residual RER activity develops to birth. Estimations of the number of rNMPs in DNA in these two mutants define a ribonucleotide threshold above which p53 induces apoptosis. Below the threshold, rNMPs in DNA trigger an innate immune response. Compound heterozygous cells, containing both defective enzymes, retain rNMPs above the threshold, indicative of competition for RER substrates between active and inactive enzymes, suggesting that patients with compound heterozygous mutations in RNASEH2 genes may not reflect the properties of recombinantly expressed proteins.

Cancer driver mutations often display mutual exclusion or co-occurrence, underscoring the key role of epistasis in carcinogenesis. However, estimating the magnitude of epistasis and quantifying its effect on tumor evolution remains a challenge. We develop a method (Coselens) to quantify conditional selection on the excess of nonsynonymous substitutions in cancer genes. Coselens infers the number of drivers per gene in different partitions of a cancer genomics dataset using covariance-based mutation models and determines whether coding mutations in a gene affect selection for drivers in any other gene. Using Coselens, we identify 296 conditionally selected gene pairs across 16 cancer types in the TCGA dataset. Conditional selection affects 25%-50% of driver substitutions in tumors with >2 drivers. Conditionally co-selected genes form modular networks, whose structures challenge the traditional interpretation of within-pathway mutual exclusivity and across-pathway synergy, suggesting a more complex scenario where gene-specific across-pathway epistasis shapes differentiated cancer subtypes.

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Ribonuclease H2 (RNase H2) protects genome integrity by its dual roles of resolving transcription-related R-loops and ribonucleotides incorporated in DNA during replication. To unlink these two functions, we generated a Saccharomyces cerevisiae RNase H2 mutant that can resolve R-loops but cannot cleave single ribonucleotides in DNA. This mutant definitively correlates the 2-5 bp deletions observed in rnh201Δ strains with single rNMPs in DNA. It also establishes a connection between R-loops and Sgs1-mediated replication reinitiation at stalled forks and identifies R-loops uniquely processed by RNase H2. In mouse, deletion of any of the genes coding for RNase H2 results in embryonic lethality, and in humans, RNase H2 hypomorphic mutations cause Aicardi-Goutières syndrome (AGS), a neuroinflammatory disorder. To determine the contribution of R-loops and rNMP in DNA to the defects observed in AGS, we characterized in yeast an AGS-related mutation, which is impaired in processing both substrates, but has sufficient R-loop degradation activity to complement the defects of rnh201Δ sgs1Δ strains. However, this AGS-related mutation accumulates 2-5 bp deletions at a very similar rate as the deletion strain.

"Candidatus Actinomarinales" was defined as a subclass of exclusively marine Actinobacteria with small cells and genomes. We have collected all the available genomes in databases to assess the diversity included in this group and analyzed it by comparative genomics. We have found the equivalent of five genera and 18 genomospecies. They have genome reduction parameters equal to those of freshwater actinobacterial "Candidatus Nanopelagicales" or marine alphaproteobacterial Pelagibacterales Genome recruitment shows that they are found only in the photic zone and mainly in surface waters, with only one genus that is found preferentially at or below the deep chlorophyll maximum. "Ca Actinomarinales" show a highly conserved core genome (80% of the gene families conserved for the whole order) with a saturation of genomic diversity of the flexible genome at the genomospecies level. We found only a flexible genomic island preserved throughout the order; it is related to the sugar decoration of the envelope and uses several tRNAs as hot spots to increase its genomic diversity. Populations had a discrete level of sequence diversity similar to other marine microbes but drastically different from the much higher levels found for Pelagibacterales Genomic analysis suggests that they are all aerobic photoheterotrophs with one type 1 rhodopsin and a heliorhodopsin. Like other actinobacteria, they possess the F420 coenzyme biosynthesis pathway, and its lower reduction potential could provide access to an increased range of redox chemical transformations. Last, sequence analysis revealed the first "Ca Actinomarinales" phages, including a prophage, with metaviromic islands related to sialic acid cleavage.IMPORTANCE Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge. Here, we study 182 genomes belonging to the only known exclusively marine pelagic group of the phylum Actinobacteria The aquatic branch of this phylum is largely known from environmental sequencing studies (single-amplified genomes [SAGs] and metagenome-assembled genomes [MAGs]), and we have collected and analyzed the available information present in databases about the "Ca. Actinomarinales." They are among the most streamlined microbes to live in the epipelagic zone of the ocean, and their study is critical to obtain a proper view of the diversity of Actinobacteria and their role in aquatic ecosystems.

Bacterial and archaeal evolution involve extensive gene gain and loss. Thus, phylogenetic trees of prokaryotes can be constructed both by traditional sequence-based methods (gene trees) and by comparison of gene compositions (genome trees). Comparing the branch lengths in gene and genome trees with identical topologies for 34 clusters of closely related bacterial and archaeal genomes, we show here that terminal branches of gene trees are systematically compressed compared to those of genome trees. Thus, sequence evolution is delayed compared to genome evolution by gene gain and loss. The extent of this delay differs widely among bacteria and archaea. Mathematical modeling shows that the divergence delay can result from sequence homogenization by homologous recombination. The model explains how homologous recombination maintains the cohesiveness of the core genome of a species while allowing extensive gene gain and loss within the accessory genome. Once evolving genomes become isolated by barriers impeding homologous recombination, gene and genome evolution processes settle into parallel trajectories, and genomes diverge, resulting in speciation.

The study of genes that evolve under conditional selection can shed light on the genomic underpinnings of adaptation, revealing epistasis and phenotypic plasticity. This protocol describes how to use the Coselens package to compare gene-level selection between two groups of samples. After installing Coselens and preparing the datasets, a typical run on a laptop takes less than 10 min. Coselens is best suited to analyze somatic mutations and data from experimental evolution, for which independently evolved samples are available. For complete details on the use and execution of this protocol, please refer to Iranzo et al. (2022).1.