The Prokaryotic Roots of Eukaryotic Immune Systems.

Over the past two decades, studies have revealed profound evolutionary connections between prokaryotic and eukaryotic immune systems, challenging the notion of their unrelatedness. Immune systems across the tree of life share an operational framework, shaping their biochemical logic and evolutionary trajectories. The diversification of immune genes in the prokaryotic superkingdoms, followed by lateral transfer to eukaryotes, was central to the emergence of innate immunity in the latter. These include protein domains related to nucleotide second messenger-dependent systems, NAD+/nucleotide degradation, and P-loop NTPase domains of the STAND and GTPase clades playing pivotal roles in eukaryotic immunity and inflammation. Moreover, several domains orchestrating programmed cell death, ultimately of prokaryotic provenance, suggest an intimate link between immunity and the emergence of multicellularity in eukaryotes such as animals. While eukaryotes directly adopted some proteins from bacterial immune systems, they repurposed others for new immune functions from bacterial interorganismal conflict systems. These emerging immune components hold substantial biotechnological potential.

Identification of a proteolysis regulator for an essential enzyme in Mycobacterium tuberculosis.

In Mycobacterium tuberculosis proteins that are post-translationally modified with Pup, a prokaryotic ubiquitin-like protein, can be degraded by proteasomes. While pupylation is reversible, mechanisms regulating substrate specificity have not been identified. Here, we identify the first depupylation regulators: CoaX, a pseudokinase, and pantothenate, an essential, central metabolite. In a Δ coaX mutant, pantothenate synthesis enzymes were more abundant, including PanB, a substrate of the Pup-proteasome system. Media supplementation with pantothenate decreased PanB levels in a coaX and Pup-proteasome-dependent manner. In vitro , CoaX accelerated depupylation of Pup∼PanB, while addition of pantothenate inhibited this reaction. Collectively, we propose CoaX contributes to proteasomal degradation of PanB by modulating depupylation of Pup∼PanB in response to pantothenate levels. One Sentence Summary: A pseudo-pantothenate kinase regulates proteasomal degradation of a pantothenate synthesis enzyme in M. tuberculosis .

Functionally comparable but evolutionarily distinct nucleotide-targeting effectors help identify conserved paradigms across diverse immune systems.

While nucleic acid-targeting effectors are known to be central to biological conflicts and anti-selfish element immunity, recent findings have revealed immune effectors that target their building blocks and the cellular energy currency-free nucleotides. Through comparative genomics and sequence-structure analysis, we identified several distinct effector domains, which we named Calcineurin-CE, HD-CE, and PRTase-CE. These domains, along with specific versions of the ParB and MazG domains, are widely present in diverse prokaryotic immune systems and are predicted to degrade nucleotides by targeting phosphate or glycosidic linkages. Our findings unveil multiple potential immune systems associated with at least 17 different functional themes featuring these effectors. Some of these systems sense modified DNA/nucleotides from phages or operate downstream of novel enzymes generating signaling nucleotides. We also uncovered a class of systems utilizing HSP90- and HSP70-related modules as analogs of STAND and GTPase domains that are coupled to these nucleotide-targeting- or proteolysis-induced complex-forming effectors. While widespread in bacteria, only a limited subset of nucleotide-targeting effectors was integrated into eukaryotic immune systems, suggesting barriers to interoperability across subcellular contexts. This work establishes nucleotide-degrading effectors as an emerging immune paradigm and traces their origins back to homologous domains in housekeeping systems.

VirB, a key transcriptional regulator of Shigella virulence, requires a CTP ligand for its regulatory activities.

IMPORTANCE: Shigella species cause bacillary dysentery, the second leading cause of diarrheal deaths worldwide. There is a pressing need to identify novel molecular drug targets. Shigella virulence phenotypes are controlled by the transcriptional regulator, VirB. We show that VirB belongs to a fast-evolving, plasmid-borne clade of the ParB superfamily, which has diverged from versions with a distinct cellular role-DNA partitioning. We report that, like classic members of the ParB family, VirB binds a highly unusual ligand, CTP. Mutants predicted to be defective in CTP binding are compromised in a variety of virulence attributes controlled by VirB, likely because these mutants cannot engage DNA. This study (i) reveals that VirB binds CTP, (ii) provides a link between VirB-CTP interactions and Shigella virulence phenotypes, (iii) provides new insight into VirB-CTP-DNA interactions, and (iv) broadens our understanding of the ParB superfamily, a group of bacterial proteins that play critical roles in many bacteria.

VirB, a key transcriptional regulator of Shigella virulence, requires a CTP ligand for its regulatory activities.

The VirB protein, encoded by the large virulence plasmid of Shigella spp., is a key transcriptional regulator of virulence genes. Without a functional virB gene, Shigella cells are avirulent. On the virulence plasmid, VirB functions to offset transcriptional silencing mediated by the nucleoid structuring protein, H-NS, which binds and sequesters AT-rich DNA, making it inaccessible for gene expression. Thus, gaining a mechanistic understanding of how VirB counters H-NS-mediated silencing is of considerable interest. VirB is unusual in that it does not resemble classic transcription factors. Instead, its closest relatives are found in the ParB superfamily, where the best-characterized members function in faithful DNA segregation before cell division. Here, we show that VirB is a fast-evolving member of this superfamily and report for the first time that the VirB protein binds a highly unusual ligand, CTP. VirB binds this nucleoside triphosphate preferentially and with specificity. Based on alignments with the best-characterized members of the ParB family, we identify amino acids of VirB likely to bind CTP. Substitutions in these residues disrupt several well-documented activities of VirB, including its anti-silencing activity at a VirB-dependent promoter, its role in generating a Congo red positive phenotype in Shigella , and the ability of the VirB protein to form foci in the bacterial cytoplasm when fused to GFP. Thus, this work is the first to show that VirB is a bona fide CTP-binding protein and links Shigella virulence phenotypes to the nucleoside triphosphate, CTP. Importance: Shigella species cause bacillary dysentery (shigellosis), the second leading cause of diarrheal deaths worldwide. With growing antibiotic resistance, there is a pressing need to identify novel molecular drug targets. Shigella virulence phenotypes are controlled by the transcriptional regulator, VirB. We show that VirB belongs to a fast-evolving, primarily plasmid-borne clade of the ParB superfamily, which has diverged from versions that have a distinct cellular role - DNA partitioning. We are the first to report that, like classic members of the ParB family, VirB binds a highly unusual ligand, CTP. Mutants predicted to be defective in CTP binding are compromised in a variety of virulence attributes controlled by VirB. This study i) reveals that VirB binds CTP, ii) provides a link between VirB-CTP interactions and Shigella virulence phenotypes, and iii) broadens our understanding of the ParB superfamily, a group of bacterial proteins that play critical roles in many different bacteria.

Apprehending the NAD+-ADPr-Dependent Systems in the Virus World.

NAD+ and ADP-ribose (ADPr)-containing molecules are at the interface of virus-host conflicts across life encompassing RNA processing, restriction, lysogeny/dormancy and functional hijacking. We objectively defined the central components of the NAD+-ADPr networks involved in these conflicts and systematically surveyed 21,191 completely sequenced viral proteomes representative of all publicly available branches of the viral world to reconstruct a comprehensive picture of the viral NAD+-ADPr systems. These systems have been widely and repeatedly exploited by positive-strand RNA and DNA viruses, especially those with larger genomes and more intricate life-history strategies. We present evidence that ADP-ribosyltransferases (ARTs), ADPr-targeting Macro, NADAR and Nudix proteins are frequently packaged into virions, particularly in phages with contractile tails (Myoviruses), and deployed during infection to modify host macromolecules and counter NAD+-derived signals involved in viral restriction. Genes encoding NAD+-ADPr-utilizing domains were repeatedly exchanged between distantly related viruses, hosts and endo-parasites/symbionts, suggesting selection for them across the virus world. Contextual analysis indicates that the bacteriophage versions of ADPr-targeting domains are more likely to counter soluble ADPr derivatives, while the eukaryotic RNA viral versions might prefer macromolecular ADPr adducts. Finally, we also use comparative genomics to predict host systems involved in countering viral ADP ribosylation of host molecules.

Discovering Biological Conflict Systems Through Genome Analysis: Evolutionary Principles and Biochemical Novelty.

Biological replicators, from genes within a genome to whole organisms, are locked in conflicts. Comparative genomics has revealed a staggering diversity of molecular armaments and mechanisms regulating their deployment, collectively termed biological conflict systems. These encompass toxins used in inter- and intraspecific interactions, self/nonself discrimination, antiviral immune mechanisms, and counter-host effectors deployed by viruses and intragenomic selfish elements. These systems possess shared syntactical features in their organizational logic and a set of effectors targeting genetic information flow through the Central Dogma, certain membranes, and key molecules like NAD+. These principles can be exploited to discover new conflict systems through sensitive computational analyses. This has led to significant advances in our understanding of the biology of these systems and furnished new biotechnological reagents for genome editing, sequencing, and beyond. We discuss these advances using specific examples of toxins, restriction-modification, apoptosis, CRISPR/second messenger-regulated systems, and other enigmatic nucleic acid-targeting systems.

Bacterial death and TRADD-N domains help define novel apoptosis and immunity mechanisms shared by prokaryotes and metazoans.

Several homologous domains are shared by eukaryotic immunity and programmed cell-death systems and poorly understood bacterial proteins. Recent studies show these to be components of a network of highly regulated systems connecting apoptotic processes to counter-invader immunity, in prokaryotes with a multicellular habit. However, the provenance of key adaptor domains, namely those of the Death-like and TRADD-N superfamilies, a quintessential feature of metazoan apoptotic systems, remained murky. Here, we use sensitive sequence analysis and comparative genomics methods to identify unambiguous bacterial homologs of the Death-like and TRADD-N superfamilies. We show the former to have arisen as part of a radiation of effector-associated α-helical adaptor domains that likely mediate homotypic interactions bringing together diverse effector and signaling domains in predicted bacterial apoptosis- and counter-invader systems. Similarly, we show that the TRADD-N domain defines a key, widespread signaling bridge that links effector deployment to invader-sensing in multicellular bacterial and metazoan counter-invader systems. TRADD-N domains are expanded in aggregating marine invertebrates and point to distinctive diversifying immune strategies probably directed both at RNA and retroviruses and cellular pathogens that might infect such communities. These TRADD-N and Death-like domains helped identify several new bacterial and metazoan counter-invader systems featuring underappreciated, common functional principles: the use of intracellular invader-sensing lectin-like (NPCBM and FGS), transcription elongation GreA/B-C, glycosyltransferase-4 family, inactive NTPase (serving as nucleic acid receptors), and invader-sensing GTPase switch domains. Finally, these findings point to the possibility of multicellular bacteria-stem metazoan symbiosis in the emergence of the immune/apoptotic systems of the latter.

GREB1: An evolutionarily conserved protein with a glycosyltransferase domain links ERα glycosylation and stability to cancer.

What covalent modifications control the temporal ubiquitination of ERα and hence the duration of its transcriptional activity remain poorly understood. We show that GREB1, an ERα-inducible enzyme, catalyzes O-GlcNAcylation of ERα at residues T553/S554, which stabilizes ERα protein by inhibiting association with the ubiquitin ligase ZNF598. Loss of GREB1-mediated glycosylation of ERα results in reduced cellular ERα levels and insensitivity to estrogen. Higher GREB1 expression in ERα+ve breast cancer is associated with greater survival in response to tamoxifen, an ERα agonist. Mice lacking Greb1 exhibit growth and fertility defects reminiscent of phenotypes in ERα-null mice. In summary, this study identifies GREB1, a protein with an evolutionarily conserved domain related to DNA-modifying glycosyltransferases of bacteriophages and kinetoplastids, as the first inducible and the only other (apart from OGT) O-GlcNAc glycosyltransferase in mammalian cytoplasm and ERα as its first substrate.

Comprehensive classification of ABC ATPases and their functional radiation in nucleoprotein dynamics and biological conflict systems.

ABC ATPases form one of the largest clades of P-loop NTPase fold enzymes that catalyze ATP-hydrolysis and utilize its free energy for a staggering range of functions from transport to nucleoprotein dynamics. Using sensitive sequence and structure analysis with comparative genomics, for the first time we provide a comprehensive classification of the ABC ATPase superfamily. ABC ATPases developed structural hallmarks that unambiguously distinguish them from other P-loop NTPases such as an alternative to arginine-finger-based catalysis. At least five and up to eight distinct clades of ABC ATPases are reconstructed as being present in the last universal common ancestor. They underwent distinct phases of structural innovation with the emergence of inserts constituting conserved binding interfaces for proteins or nucleic acids and the adoption of a unique dimeric toroidal configuration for DNA-threading. Specifically, several clades have also extensively radiated in counter-invader conflict systems where they serve as nodal nucleotide-dependent sensory and energetic components regulating a diversity of effectors (including some previously unrecognized) acting independently or together with restriction-modification systems. We present a unified mechanism for ABC ATPase function across disparate systems like RNA editing, translation, metabolism, DNA repair, and biological conflicts, and some unexpected recruitments, such as MutS ATPases in secondary metabolism.

Mycobacterium tuberculosis Rv0991c Is a Redox-Regulated Molecular Chaperone.

The bacterial pathogen Mycobacterium tuberculosis is the leading cause of death by an infectious disease among humans. Here, we describe a previously uncharacterized M. tuberculosis protein, Rv0991c, as a molecular chaperone that is activated by oxidation. Rv0991c has homologs in most bacterial lineages and appears to function analogously to the well-characterized Escherichia coli redox-regulated chaperone Hsp33, despite a dissimilar protein sequence. Rv0991c is transcriptionally coregulated with hsp60 and hsp70 chaperone genes in M. tuberculosis, suggesting that Rv0991c functions with these chaperones in maintaining protein quality control. Supporting this hypothesis, we found that, like oxidized Hsp33, oxidized Rv0991c prevents the aggregation of a model unfolded protein in vitro and promotes its refolding by the M. tuberculosis Hsp70 chaperone system. Furthermore, Rv0991c interacts with DnaK and can associate with many other M. tuberculosis proteins. We therefore propose that Rv0991c, which we named "Ruc" (redox-regulated protein with unstructured C terminus), represents a founding member of a new chaperone family that protects M. tuberculosis and other species from proteotoxicity during oxidative stress.IMPORTANCEM. tuberculosis infections are responsible for more than 1 million deaths per year. Developing effective strategies to combat this disease requires a greater understanding of M. tuberculosis biology. As in all cells, protein quality control is essential for the viability of M. tuberculosis, which likely faces proteotoxic stress within a host. Here, we identify an M. tuberculosis protein, Ruc, that gains chaperone activity upon oxidation. Ruc represents a previously unrecognized family of redox-regulated chaperones found throughout the bacterial superkingdom. Additionally, we found that oxidized Ruc promotes the protein-folding activity of the essential M. tuberculosis Hsp70 chaperone system. This work contributes to a growing body of evidence that oxidative stress provides a particular strain on cellular protein stability.

Functional Innovation in the Evolution of the Calcium-Dependent System of the Eukaryotic Endoplasmic Reticulum.

The origin of eukaryotes was marked by the emergence of several novel subcellular systems. One such is the calcium (Ca2+)-stores system of the endoplasmic reticulum, which profoundly influences diverse aspects of cellular function including signal transduction, motility, division, and biomineralization. We use comparative genomics and sensitive sequence and structure analyses to investigate the evolution of this system. Our findings reconstruct the core form of the Ca2+-stores system in the last eukaryotic common ancestor as having at least 15 proteins that constituted a basic system for facilitating both Ca2+ flux across endomembranes and Ca2+-dependent signaling. We present evidence that the key EF-hand Ca2+-binding components had their origins in a likely bacterial symbiont other than the mitochondrial progenitor, whereas the protein phosphatase subunit of the ancestral calcineurin complex was likely inherited from the asgard archaeal progenitor of the stem eukaryote. This further points to the potential origin of the eukaryotes in a Ca2+-rich biomineralized environment such as stromatolites. We further show that throughout eukaryotic evolution there were several acquisitions from bacteria of key components of the Ca2+-stores system, even though no prokaryotic lineage possesses a comparable system. Further, using quantitative measures derived from comparative genomics we show that there were several rounds of lineage-specific gene expansions, innovations of novel gene families, and gene losses correlated with biological innovation such as the biomineralized molluscan shells, coccolithophores, and animal motility. The burst of innovation of new genes in animals included the wolframin protein associated with Wolfram syndrome in humans. We show for the first time that it contains previously unidentified Sel1, EF-hand, and OB-fold domains, which might have key roles in its biochemistry.

Corrigendum: The cyl Genes Reveal the Biosynthetic and Evolutionary Origins of the Group B Streptococcus Hemolytic Lipid, Granadaene.

[This corrects the article DOI: 10.3389/fmicb.2019.03123.].

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity.

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

Bacteria are the most numerous lifeforms on the planet. Most bacteria live as single cells that grow and multiply independently within larger communities of microbes. However, some bacterial cells assemble to form more complex structures where individual cells might perform distinct roles. Such bacteria are referred to as 'multicellular bacteria'. For example, cells of bacteria known as Streptomyces collectively form filaments that help the bacteria collect nutrients from their food sources, and aerial structures bearing reproductive spores. Bacteria in these filaments may come into contact with many other microbes in their surroundings including other bacteria within the same filament, other species of bacteria, and viruses. These contacts often lead to conflict, for example, if the microbes compete with each other for nutrients or if a virus tries to attack the bacteria. Bacteria have evolved immune systems that detect other microbes and use antibiotics, toxins and other defense mechanisms against them. Compared to single-celled bacteria, multicellular bacteria may be more vulnerable to threats from viruses because once a virus has overcome the defenses of one cell in the multicellular assembly, it may be easier for it to kill, or spread to the other cells. However, it is not clear how these systems evolved to deal with the unique problems of multicellular bacteria. Now, Kaur, Burroughs et al. have used computational approaches to search for new immune systems in diverse multicellular bacteria. The new classes of systems they found are each made of different molecular components, but all require a large input of energy to be activated. This activation barrier prevents the bacterial cells from deploying weapons unless the signal from the enemy microbe crosses a high enough threshold. Many tools used in molecular biology, and increasingly in medicine, have been derived from the immune systems of bacteria, such as the enzymes that cut or edit DNA. The findings of Kaur, Burroughs et al. may aid the development of new tools that specifically bind to viruses or other dangerous microbes, or inhibit their ability to interact with components in cells. The next step would be to perform experiments using some of the immune systems identified in this work.

The catalytic core of DEMETER guides active DNA demethylation in Arabidopsis.

The Arabidopsis DEMETER (DME) DNA glycosylase demethylates the maternal genome in the central cell prior to fertilization and is essential for seed viability. DME preferentially targets small transposons that flank coding genes, influencing their expression and initiating plant gene imprinting. DME also targets intergenic and heterochromatic regions, but how it is recruited to these differing chromatin landscapes is unknown. The C-terminal half of DME consists of 3 conserved regions required for catalysis in vitro. We show that this catalytic core guides active demethylation at endogenous targets, rescuing dme developmental and genomic hypermethylation phenotypes. However, without the N terminus, heterochromatin demethylation is significantly impeded, and abundant CG-methylated genic sequences are ectopically demethylated. Comparative analysis revealed that the conserved DME N-terminal domains are present only in flowering plants, whereas the domain architecture of DME-like proteins in nonvascular plants mainly resembles the catalytic core, suggesting that it might represent the ancestral form of the 5mC DNA glycosylase found in plant lineages. We propose a bipartite model for DME protein action and suggest that the DME N terminus was acquired late during land plant evolution to improve specificity and facilitate demethylation at heterochromatin targets.

The cyl Genes Reveal the Biosynthetic and Evolutionary Origins of the Group B Streptococcus Hemolytic Lipid, Granadaene.

Group B Streptococcus (GBS) is a ß-hemolytic, Gram-positive bacterium that commonly colonizes the female lower genital tract and is associated with fetal injury, preterm birth, spontaneous abortion, and neonatal infections. A major factor promoting GBS virulence is the ß-hemolysin/cytolysin, which is cytotoxic to several host cells. We recently showed that the ornithine rhamnolipid pigment, Granadaene, produced by the gene products of the cyl operon, is hemolytic. Here, we demonstrate that heterologous expression of the GBS cyl operon conferred hemolysis, pigmentation, and cytoxicity to Lactococcus lactis, a model non-hemolytic Gram-positive bacterium. Similarly, pigment purified from L. lactis is hemolytic, cytolytic, and identical in structure to Granadaene extracted from GBS, indicating the cyl operon is sufficient for Granadaene production in a heterologous host. Using a systematic survey of phyletic patterns and contextual associations of the cyl genes, we identify homologs of the cyl operon in physiologically diverse Gram-positive bacteria and propose undescribed functions of cyl gene products. Together, these findings bring greater understanding to the biosynthesis and evolutionary foundations of a key GBS virulence factor and suggest that such potentially toxic lipids may be encoded by other bacteria.

Unexpected Evolution of Lesion-Recognition Modules in Eukaryotic NER and Kinetoplast DNA Dynamics Proteins from Bacterial Mobile Elements.

The provenance of several components of major uniquely eukaryotic molecular machines are increasingly being traced back to prokaryotic biological conflict systems. Here, we demonstrate that the N-terminal single-stranded DNA-binding domain from the anti-restriction protein ArdC, deployed by bacterial mobile elements against their host, was independently acquired twice by eukaryotes, giving rise to the DNA-binding domains of XPC/Rad4 and the Tc-38-like proteins in the stem kinetoplastid. In both instances, the ArdC-N domain tandemly duplicated forming an extensive DNA-binding interface. In XPC/Rad4, the ArdC-N domains (BHDs) also fused to the inactive transglutaminase domain of a peptide-N-glycanase ultimately derived from an archaeal conflict system. Alongside, we delineate several parallel acquisitions from conjugative elements/bacteriophages that gave rise to key components of the kinetoplast DNA (kDNA) replication apparatus. These findings resolve two outstanding questions in eukaryote biology: (1) the origin of the unique DNA lesion-recognition component of NER and (2) origin of the unusual, plasmid-like features of kDNA.

Evolutionary convergence and divergence in archaeal chromosomal proteins and Chromo-like domains from bacteria and eukaryotes.

SH3-fold-ß-barrel domains of the chromo-like superfamily recognize epigenetic marks in eukaryotic proteins. Their provenance has been placed either in archaea, based on apparent structural similarity to chromatin-compacting Sul7d and Cren7 proteins, or in bacteria based on the presence of sequence homologs. Using sequence and structural evidence we establish that the archaeal Cren7/Sul7 proteins emerged from a zinc ribbon (ZnR) ancestor. Further, we show that the ancestral eukaryotic chromo-like domains evolved from bacterial versions, likely acquired from early endosymbioses, which already possessed an aromatic cage for recognition of modified amino-groups. These bacterial versions are part of a radiation of secreted SH3-fold domains, which spawned both chromo-like domains and classical SH3 domains in the context of peptide-recognition in the peptidoglycan or the extracellular matrix. This establishes that Cren7/Sul7 converged to a "SH3"-like state from a ZnR precursor via the loss of metal-chelation and acquisition of stronger hydrophobic interactions; it is unlikely to have participated in the evolution of the chromo-like domains. We show that archaea possess several Cren7/Sul7-related proteins with intact Zn-chelating ligands, which we predict to play previously unstudied roles in chromosome segregation during cell-division comparable to the PRC barrel and CdvA domain proteins.

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.

Expansions, diversification, and interindividual copy number variations of AID/APOBEC family cytidine deaminase genes in lampreys.

Cytidine deaminases of the AID/APOBEC family catalyze C-to-U nucleotide transitions in mRNA or DNA. Members of the APOBEC3 branch are involved in antiviral defense, whereas AID contributes to diversification of antibody repertoires in jawed vertebrates via somatic hypermutation, gene conversion, and class switch recombination. In the extant jawless vertebrate, the lamprey, two members of the AID/APOBEC family are implicated in the generation of somatic diversity of the variable lymphocyte receptors (VLRs). Expression studies linked CDA1 and CDA2 genes to the assembly of VLRA/C genes in T-like cells and the VLRB genes in B-like cells, respectively. Here, we identify and characterize several CDA1-like genes in the larvae of different lamprey species and demonstrate that these encode active cytidine deaminases. Structural comparisons of the CDA1 variants highlighted substantial differences in surface charge; this observation is supported by our finding that the enzymes require different conditions and substrates for optimal activity in vitro. Strikingly, we also found that the number of CDA-like genes present in individuals of the same species is variable. Nevertheless, irrespective of the number of different CDA1-like genes present, all lamprey larvae have at least one functional CDA1-related gene encoding an enzyme with predicted structural and chemical features generally comparable to jawed vertebrate AID. Our findings suggest that, similar to APOBEC3 branch expansion in jawed vertebrates, the AID/APOBEC family has undergone substantial diversification in lamprey, possibly indicative of multiple distinct biological roles.