The preferred nucleotide contexts of the AID/APOBEC cytidine deaminases have differential effects when mutating retrotransposon and virus sequences compared to host genes

Role and molecular mechanism of APOBEC3B in the development and progression of gastric cancer

Gastric cancer is a common malignant tumor with a high mortality rate. Abnormal APOBEC3B (apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3B) expression increases tumor susceptibility. However, the exact molecular mechanism of APOBEC3B expression in the development of gastric cancer is still unknown. We investigated the effect of APOBEC3B on the malignant biological behavior of gastric cancer cells and discussed the role of APOBEC3B in the development and progression of gastric cancer. APOBEC3B protein levels were measured in 161 gastric cancer samples using western blotting and immunohistochemistry. Both in vitro and in vivo assays were performed, and molecules were analyzed using bioinformatics analysis and western blotting. APOBEC3B was overexpressed in gastric cancer. Moreover, APOBEC3B significantly enhanced cell proliferation in vitro and tumorigenicity in vivo. Regarding the underlying mechanism, APOBEC3B promoted the proliferation of gastric cancer cells by upregulating P53, MCM2 (minichromosome maintenance protein 2), and cyclin D1. Our results suggest that APOBEC3B is involved in cancer progression, providing a new theoretical basis for the prevention and treatment of gastric cancer.

Engineered deaminases as a key component of DNA and RNA editing tools

Over recent years, zinc-dependent deaminases have attracted increasing interest as key components of nucleic acid editing tools that can generate point mutations at specific sites in either DNA or RNA by combining a targeting module (such as a catalytically impaired CRISPR-Cas component) and an effector module (most often a deaminase). Deaminase-based molecular tools are already being utilized in a wide spectrum of therapeutic and research applications; however, their medical and biotechnological potential seems to be much greater. Recent reports indicate that the further development of nucleic acid editing systems depends largely on our ability to engineer the substrate specificity and catalytic activity of the editors themselves. In this review, we summarize the current trends and achievements in deaminase engineering. The presented data indicate that the potential of these enzymes has not yet been fully revealed or understood. Several examples show that even relatively minor changes in the structure of deaminases can give them completely new and unique properties.

Apobec-Mediated Retroviral Hypermutation In Vivo is Dependent on Mouse Strain

Replication of the complex retrovirus mouse mammary tumor virus (MMTV) is antagonized by murine Apobec3 (mA3), a member of the Apobec family of cytidine deaminases. We have shown that MMTV-encoded Rem protein inhibits proviral mutagenesis by the Apobec enzyme, activation-induced cytidine deaminase (AID) during viral replication in BALB/c mice. To further study the role of Rem in vivo , we have infected C57BL/6 (B6) mice with a superantigen-independent lymphomagenic strain of MMTV (TBLV-WT) or a mutant strain (TBLV-SD) that is defective in Rem and its cleavage product Rem-CT. Unlike MMTV, TBLV induced T-cell tumors in µMT mice, indicating that mature B cells, which express the highest AID levels, are not required for TBLV replication. Compared to BALB/c, B6 mice were more susceptible to TBLV infection and tumorigenesis. The lack of Rem expression accelerated B6 tumorigenesis at limiting doses compared to TBLV-WT in either wild-type B6 or AID-deficient mice. However, unlike proviruses from BALB/c mice, high-throughput sequencing indicated that proviral G-to-A or C-to-T changes did not significantly differ in the presence and absence of Rem expression. Ex vivo stimulation showed higher levels of mA3 relative to AID in B6 compared to BALB/c splenocytes, but effects of agonists differed in the two strains. RNA-Seq revealed increased transcripts related to growth factor and cytokine signaling in TBLV-SD-induced tumors relative to those from TBLV-WT, consistent with a third Rem function. Thus, Rem-mediated effects on tumorigenesis in B6 mice are independent of Apobec-mediated proviral hypermutation.

Accelerated evolution of SARS-CoV-2 in free-ranging white-tailed deer

The zoonotic origin of the COVID-19 pandemic virus highlights the need to fill the vast gaps in our knowledge of SARS-CoV-2 ecology and evolution in non-human hosts. Here, we detected that SARS-CoV-2 was introduced from humans into white-tailed deer more than 30 times in Ohio, USA during November 2021-March 2022. Subsequently, deer-to-deer transmission persisted for 2-8 months, disseminating across hundreds of kilometers. Newly developed Bayesian phylogenetic methods quantified how SARS-CoV-2 evolution is not only three-times faster in white-tailed deer compared to the rate observed in humans but also driven by different mutational biases and selection pressures. The long-term effect of this accelerated evolutionary rate remains to be seen as no critical phenotypic changes were observed in our animal models using white-tailed deer origin viruses. Still, SARS-CoV-2 has transmitted in white-tailed deer populations for a relatively short duration, and the risk of future changes may have serious consequences for humans and livestock.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational 'hotspots', can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at TC motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of TC hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more TC hotspots than expected by chance (over-representation) and fewer GC hotspots than expected (under-representation). These results suggest the MPXV genome: (1) may have evolved in a host with a particular APOBEC GC hotspot preference, (2) has inverted terminal repeat (ITR) regions-which may be exposed to APOBEC3 for longer during viral replication-and longer genes likely to evolve faster, and therefore (3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational "hotspots," can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at T C motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of T C hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more T C hotspots than expected by chance (over-representation) and fewer G C hotspots than expected (under-representation). These results suggest the MPXV genome: 1) may have evolved in a host with a particular APOBEC G C hotspot preference, 2) has inverted terminal repeat (ITR) regions -which may be exposed to APOBEC3 for longer during viral replication- and longer genes likely to evolve faster, and therefore 3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

Accelerated evolution of SARS-CoV-2 in free-ranging white-tailed deer

While SARS-CoV-2 has sporadically infected a wide range of animal species worldwide1, the virus has been repeatedly and frequently detected in white-tailed deer in North America2â€"7. The zoonotic origins of this pandemic virus highlight the need to fill the vast gaps in our knowledge of SARS-CoV-2 ecology and evolution in non-human hosts. Here, we detected SARS-CoV-2 was introduced from humans into white-tailed deer more than 30 times in Ohio, USA during November 2021-March 2022. Subsequently, deer-to-deer transmission persisted for 2-8 months, which disseminated across hundreds of kilometers. We discovered that alpha and delta variants evolved in white-tailed deer at three-times the rate observed in humans. Newly developed Bayesian phylogenetic methods quantified how SARS-CoV-2 evolution is not only faster in white-tailed deer but driven by different mutational biases and selection pressures. White-tailed deer are not just short-term recipients of human viral diversity but serve as reservoirs for alpha and other variants to evolve in new directions after going extinct in humans. The long-term effect of this accelerated evolutionary rate remains to be seen as no critical phenotypic changes were observed in our animal model experiments using viruses isolated from white-tailed deer. Still, SARS-CoV-2 viruses have transmitted in white-tailed deer populations for a relatively short duration, and the risk of future changes may have serious consequences for humans and livestock.

Characterizing microRNA editing and mutation sites in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder whose pathogenesis is still unclear. MicroRNAs (miRNAs) are a kind of endogenous small non-coding RNAs that play important roles in the post-transcriptional regulation of genes. Recent researches show that miRNAs are edited in multiple ways especially in central nervous systems. A-to-I editing of RNA catalyzed by Adenosine deaminases acting on RNA (ADARs) happens intensively in brain and is also noticed in other organs and tissues. Although miRNAs are widely edited in human brain, miRNA editing in ASD is still largely unexplored. In order to reveal the editing events of miRNAs in ASD, we analyzed 131 miRNA-seq samples from 8 different brain regions of ASD patients and normal controls. We identified 834 editing sites with significant editing levels, of which 70 sites showed significantly different editing levels in the superior frontal gyrus samples of ASD patients (ASD-SFG) when compared with those of control samples. The editing level of an A-to-I editing site in hsa-mir-376a-1 (hsa-mir-376a-1_9_A_g) in ASD-SFG is higher than that of normal controls, and the difference is exaggerated in individuals under 10 years. The increased expression of ADAR1 is consistent with the increased editing level of hsa-mir-376a-1_9_A_g in ASD-SFG samples compared to normal SFG samples. Furthermore, we verify that A-to-I edited hsa-mir-376a-5p directly represses GPR85 and NAPB, which may contribute to the abnormal neuronal development of ASD patients. These results provide new insights into the mechanism of ASD.

Diversification of immunoglobulin genes by gene conversion in the domestic chicken (Gallus gallus domesticus)

Sustainable modern poultry production depends on effective protection against infectious diseases and a diverse range of antibodies is key for an effective immune response. In the domestic chicken, somatic gene conversion is the dominant process in which the antibody immunoglobulin genes are diversified. Affinity maturation by somatic hypermutation (SHM) also occurs, but the relative contribution of gene conversion versus somatic hypermutation to immunoglobulin (Ig) gene diversity is poorly understood. In this study, we use high throughput long-read sequencing to study immunoglobulin diversity in multiple immune-associated tissues in Rhode Island Red chickens. To better understand the impact of genetic diversification in the chicken, a novel gene conversion identification software was developed (BrepConvert). In this study, BrepConvert enabled the identification of over 1 million gene conversion events. Mapping the occurrence of putative somatic gene conversion (SGC) events throughout the variable gene region revealed repetitive and highly restricted patterns of genetic insertions in both the antibody heavy and light chains. These patterns coincided with the locations of genetic variability in available pseudogenes and align with antigen binding sites, predominately the complementary determining regions (CDRs). We found biased usage of pseudogenes during gene conversion, as well as immunoglobulin heavy chain diversity gene (IGHD) preferences during V(D)J gene rearrangement, suggesting that antibody diversification in chickens is more focused than the genetic potential for diversity would suggest.

Characterizing Relevant MicroRNA Editing Sites in Parkinson's Disease

MicroRNAs (miRNAs) are extensively edited in human brains. However, the functional relevance of the miRNA editome is largely unknown in Parkinson's disease (PD). By analyzing small RNA sequencing profiles of brain tissues of 43 PD patients and 88 normal controls, we found that the editing levels of five A-to-I and two C-to-U editing sites are significantly correlated with the ages of normal controls, which is disrupted in PD patients. We totally identified 362 miRNA editing sites with significantly different editing levels in prefrontal cortices of PD patients (PD-PC) compared to results of normal controls. We experimentally validated that A-to-I edited miR-497-5p, with significantly higher expression levels in PD-PC compared to normal controls, directly represses OPA1 and VAPB. Furthermore, overexpression of A-to-I edited miR-497-5p downregulates OPA1 and VAPB in two cell lines, and inhibits proliferation of glioma cells. These results suggest that the hyperediting of miR-497-5p in PD contributes to enhanced progressive neurodegeneration of PD patients. Our results provide new insights into the mechanistic understanding, novel diagnostics, and therapeutic clues of PD.

Control of Human Anelloviruses by Cytosine to Uracil Genome Editing

Anelloviruses are the most common viruses infecting humans. Every human carries a nonpathogenic personal anellovirus virome (anellome), yet it is unknown which mechanisms contribute to its stability. Here, we assessed the dynamics and impact of a host antiviral defense mechanism-cytidine deaminase activity leading to C to U editing in anelloviruses-on the stability of the anellome. We investigated anellome sequence data obtained from serum samples collected every 6 months from two healthy subjects followed for more than 30 years. The subjects were infected by a total of 64 anellovirus lineages. Minus-stranded C to U editing was observed in lineages belonging to the Alpha-, Beta-, and Gammatorquevirus genera. The edited genomes were present within virus particles, therefore editing must have occurred at the late stages of the virus life cycle. Editing was favored by 5'-TC contexts in the virus genome, indicating that apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like, catalytic subunit 3 or A3 (APOBEC3) proteins are involved. Within a lineage, mutational dynamics varied over time and few fixations of mutations were detected, indicating that C to U editing is a dead end for a virus genome. We detected an editing coldspot in the GC-rich regions, suggesting that the GC-rich region is crucial for genome packaging, since only packaged virus particles were included in the analysis. Finally, we noticed a lineage-specific reduced concentration after an editing event, yet no clearance. In conclusion, cytidine deaminase activity does not clear anelloviruses, nor does it play a major role in virus evolution, but it does contribute to the stability of the anellome. IMPORTANCE Despite significant attention on anellovirus research, the interaction between the anellovirus virome and the human host remains unknown. We show the dynamics of APOBEC3-mediated cytidine deaminase activity on anelloviruses during a 30-year period of chronic infection and postulate that this antiviral mechanism controls anelloviruses. These results expand our knowledge of anellovirus-host interactions, which may be important for the design of gene therapies.

Chickens, more than humans, focus the diversity of their immunoglobulin genes on the complementarity-determining region but utilise amino acids, indicative of a more cross-reactive antibody repertoire

The mechanisms of B-cell diversification differ greatly between aves and mammals, but both produce B cells and antibodies capable of supporting an effective immune response. To see how differences in the generation of diversity might affect overall repertoire diversity, we have compared the diversity characteristics of immunoglobulin genes from domestic chickens to those from humans. Both use V(D)J gene rearrangement and somatic hypermutation, but only chickens use somatic gene conversion. A range of diversity analysis tools were used to investigate multiple aspects of amino acid diversity at both the germline and repertoire levels. The effect of differing amino acid usages on antibody characteristics was assessed. At both the germline and repertoire levels, chickens exhibited lower amino acid diversity in comparison to the human immunoglobulin genes, especially outside of the complementarity-determining region (CDR). Chickens were also found to possess much larger and more hydrophilic CDR3s with a higher predicted protein binding potential, suggesting that the antigen-binding site in chicken antibodies is more flexible and more polyreactive than that seen in human antibodies.

Identification of natural product 3, 5-diiodotyrosine as APOBEC3B inhibitor to prevent somatic mutation accumulation and cancer progression

BACKGROUND

The development of cancer is largely dependent on the accumulation of somatic mutations, indicating the potential to develop cancer chemoprevention agents targeting mutation drivers. However, ideal cancer chemoprevention agents that can effectively inhibit the mutation drivers have not been identified yet.

METHODS

The somatic mutation signatures and expression analyses of APOBEC3B were performed in patient with pan-cancer. The computer-aided screening and skeleton-based searching were performed to identify natural products that can inhibit the activity of APOBEC3B. 4-nitroquinoline-1-oxide (4-NQO)-induced spontaneous esophageal squamous cell carcinoma (ESCC) and azoxymethane/dextran sulfate sodium (AOM/DSS)-induced spontaneous colon cancer mouse models were conducted to investigate the influences of APOBEC3B inhibitor on the prevention of somatic mutation accumulation and cancer progression.

RESULTS

Here, we discovered that the cytidine deaminase APOBEC3B correlated somatic mutations were widely observed in a variety of cancers, and its overexpression indicated poor survival. SMC247 (3, 5-diiodotyrosine), as a source of kelp iodine without side effects, could strongly bind APOBEC3B (K_D=65 nM) and effectively inhibit its deaminase activity (IC₅₀=1.69 µM). Interestingly, 3, 5-diiodotyrosine could significantly reduce the clusters of mutations, prevent the precancerous lesion progression, and prolong the survival in 4-NQO-induced spontaneous ESCC and AOM/DSS-induced spontaneous colon cancer mouse models. Furthermore, 3, 5-diiodotyrosine could reduce colitis, increase the proportion and function of T lymphocytes via IL-15 in tumor microenvironment. The synergistic cancer prevention effects were observed when 3, 5-diiodotyrosine combined with PD-1/PD-L1 blockade.

CONCLUSIONS

This is the first prove-of-concept study to elucidate that the natural product 3, 5-diiodotyrosine could prevent somatic mutation accumulation and cancer progression through inhibiting the enzymatic activity of APOBEC3B. In addition, 3, 5-diiodotyrosine could reduce the colitis and increase the infiltration and function of T lymphocytes via IL-15 in tumor microenvironment. 3, 5-diiodotyrosine combined with PD-1/PD-L1 blockade could elicit synergistic cancer prevention effects, indicating a novel strategy for both prevent the somatic mutation accumulation and the immune-suppressive microenvironment exacerbation.

DYW deaminase domain has a distinct preference for neighboring nucleotides of the target RNA editing sites

C-to-U RNA editing sites in plant organelles show a strong bias for neighboring nucleotides. The nucleotide upstream of the target cytidine is typically C or U, whereas A and G are less common and rare, respectively. In pentatricopeptide repeat (PPR)-type RNA editing factors, the PPR domain specifically binds to the 5' sequence of target cytidines, whereas the DYW domain catalyzes the C-to-U deamination. We comprehensively analyzed the effects of neighboring nucleotides of the target cytidines using an Escherichia coli orthogonal system. Physcomitrium PPR56 efficiently edited target cytidines when the nucleotide upstream was U or C, whereas it barely edited when the position was G or the nucleotide downstream was C. This preference pattern, which corresponds well with the observed nucleotide bias for neighboring nucleotides in plant organelles, was altered when the DYW domain of OTP86 or DYW1 was adopted. The PPR56 chimeric proteins edited the target sites even when the -1 position was G. Our results suggest that the DYW domain possesses a distinct preference for the neighboring nucleotides of the target sites, thus contributing to target selection in addition to the existing selection determined by the PPR domain.

Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus

The largest monkeypox virus (MPXV) outbreak described so far in non-endemic countries was identified in May 2022 (refs. 1-6). In this study, shotgun metagenomics allowed the rapid reconstruction and phylogenomic characterization of the first MPXV outbreak genome sequences, showing that this MPXV belongs to clade 3 and that the outbreak most likely has a single origin. Although 2022 MPXV (lineage B.1) clustered with 2018-2019 cases linked to an endemic country, it segregates in a divergent phylogenetic branch, likely reflecting continuous accelerated evolution. An in-depth mutational analysis suggests the action of host APOBEC3 in viral evolution as well as signs of potential MPXV human adaptation in ongoing microevolution. Our findings also indicate that genome sequencing may provide resolution to track the spread and transmission of this presumably slow-evolving double-stranded DNA virus.

Unusual mammalian usage of TGA stop codons reveals that sequence conservation need not imply purifying selection

The assumption that conservation of sequence implies the action of purifying selection is central to diverse methodologies to infer functional importance. GC-biased gene conversion (gBGC), a meiotic mismatch repair bias strongly favouring GC over AT, can in principle mimic the action of selection, this being thought to be especially important in mammals. As mutation is GC→AT biased, to demonstrate that gBGC does indeed cause false signals requires evidence that an AT-rich residue is selectively optimal compared to its more GC-rich allele, while showing also that the GC-rich alternative is conserved. We propose that mammalian stop codon evolution provides a robust test case. Although in most taxa TAA is the optimal stop codon, TGA is both abundant and conserved in mammalian genomes. We show that this mammalian exceptionalism is well explained by gBGC mimicking purifying selection and that TAA is the selectively optimal codon. Supportive of gBGC, we observe (i) TGA usage trends are consistent at the focal stop codon and elsewhere (in UTR sequences); (ii) that higher TGA usage and higher TAA→TGA substitution rates are predicted by a high recombination rate; and (iii) across species the difference in TAA <-> TGA substitution rates between GC-rich and GC-poor genes is largest in genomes that possess higher between-gene GC variation. TAA optimality is supported both by enrichment in highly expressed genes and trends associated with effective population size. High TGA usage and high TAA→TGA rates in mammals are thus consistent with gBGC’s predicted ability to “drive” deleterious mutations and supports the hypothesis that sequence conservation need not be indicative of purifying selection. A general trend for GC-rich trinucleotides to reside at frequencies far above their mutational equilibrium in high recombining domains supports the generality of these results.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

MicroRNA editing patterns in Huntington's disease

Huntington’s disease (HD) is a neurodegenerative disease. MicroRNAs (miRNAs) are small non-coding RNAs that mediate post-transcriptional regulation of target genes. Although miRNAs are extensively edited in human brains, the editome of miRNAs in brains of HD patients is largely unknown. By analyzing the small RNA sequencing profiles of brain tissues of 28 HD patients and 83 normal controls, 1182 miRNA editing sites with significant editing levels were identified. In addition to 27 A-to-I editing sites, we identified 3 conserved C-to-U editing sites in miRNAs of HD patients. 30 SNPs in the miRNAs of HD patients were also identified. Furthermore, 129 miRNA editing events demonstrated significantly different editing levels in prefrontal cortex samples of HD patients (HD-PC) when compared to those of healthy controls. We found that hsa-mir-10b-5p was edited to have an additional cytosine at 5’-end in HD-PC, and the edited hsa-mir-10b repressed GTPBP10 that was often downregulated in HD. The down-regulation of GTPBP10 might contribute to the progression of HD by causing gradual loss of function of mitochondrial. These results provide the first endeavor to characterize the miRNA editing events in HD and their potential functions.

Understanding memory B cell selection

The mammalian adaptive immune system has evolved over millions of years to become an incredibly effective defense against foreign antigens. The adaptive immune system's humoral response creates plasma B cells and memory B cells, each with their own immunological objectives. The affinity maturation process is widely viewed as a heuristic to solve the global optimization problem of finding B cells with high affinity to the antigen. However, memory B cells appear to be purposely selected earlier in the affinity maturation process and have lower affinity. We propose that this memory B cell selection process may be an approximate solution to two optimization problems: optimizing for affinity to similar antigens in the future despite mutations or other minor differences, and optimizing to warm start the generation of plasma B cells in the future. We use simulations to provide evidence for our hypotheses, taking into account data showing that certain B cell mutations are more likely than others. Our findings are consistent with memory B cells having high-affinity to mutated antigens, but do not provide strong evidence that memory B cells will be more useful than selected naive B cells for seeding the secondary germinal centers.

Stepwise Evolution and Exceptional Conservation of ORF1a/b Overlap in Coronaviruses

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for the propagation of coronaviruses. The combination of genomic features that make up PFE-the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements-places severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess the evolutionary dynamics of PFE in great detail. Here, we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and ORF1b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well-defined phylogenetic specificity. We further examined sequencing data from over 1,500,000 complete genomes and 55,000 raw read data sets to demonstrate exceptional conservation and detect signatures of selection within the PFE region.

Global quantification exposes abundant low-level off-target activity by base editors

Base editors are dedicated engineered deaminases that enable directed conversion of specific bases in the genome or transcriptome in a precise and efficient manner, and hold promise for correcting pathogenic mutations. A major concern limiting application of this powerful approach is the issue of off-target edits. Several recent studies have shown substantial off-target RNA activity induced by base editors and demonstrated that off-target mutations may be suppressed by improved deaminases versions or optimized guide RNAs. Here, we describe a new class of off-target events that are invisible to the established methods for detection of genomic variations and were thus far overlooked. We show that nonspecific, seemingly stochastic, off-target events affect a large number of sites throughout the genome or the transcriptome, and account for the majority of off-target activity. We develop and employ a different, complementary approach that is sensitive to the stochastic off-target activity and use it to quantify the abundant off-target RNA mutations due to current, optimized deaminase editors. We provide a computational tool to quantify global off-target activity, which can be used to optimize future base editors. Engineered base editors enable directed manipulation of the genome or transcriptome at single-base resolution. We believe that implementation of this computational approach would facilitate design of more specific base editors.

APOBEC3C, a nucleolar protein induced by genotoxins, is excluded from DNA damage sites

The human genome contains 11 APOBEC (apolipoprotein B mRNA editing catalytic polypeptide‐like) cytidine deaminases classified into four families. These proteins function mainly in innate antiviral immunity and can also restrict endogenous retrotransposable element multiplication. The present study focuses on APOBEC3C (A3C), a member of the APOBEC3 subfamily. Some APOBEC3 proteins use their enzymatic activity on genomic DNA, inducing mutations and DNA damage, while other members facilitate DNA repair. Our results show that A3C is highly expressed in cells treated with DNA‐damaging agents. Its expression is regulated by p53. Depletion of A3C slightly decreases proliferation and does not affect DNA repair via homologous recombination or nonhomologous end joining. The A3C interactomes obtained from control cells and cells exposed to the genotoxin etoposide indicated that A3C is a nucleolar protein. This was confirmed by the detection of either endogenous or ectopic A3C in nucleoli. Interestingly, we show that A3C is excluded from areas of DNA breaks in live cells. Our data also indicate that the C‐terminal part of A3C is responsible for its nucleolar localization and exclusion from DNA damage sites.

AID and APOBECs as Multifaceted Intrinsic Virus-Restricting Factors: Emerging Concepts in the Light of COVID-19

The AID (activation-induced cytidine deaminase)/APOBEC (apolipoprotein B mRNA editing enzyme catalytic subunit) family with its multifaceted mode of action emerges as potent intrinsic host antiviral system that acts against a variety of DNA and RNA viruses including coronaviruses. All family members are cytosine-to-uracil deaminases that either have a profound role in driving a strong and specific humoral immune response (AID) or restricting the virus itself by a plethora of mechanisms (APOBECs). In this article, we highlight some of the key aspects apparently linking the AID/APOBECs and SARS-CoV-2. Among those is our discovery that APOBEC4 shows high expression in cell types and anatomical parts targeted by SARS-CoV-2. Additional focus is given by us to the lymphoid structures and AID as the master regulator of germinal center reactions, which result in antibody production by plasma and memory B cells. We propose the dissection of the AID/APOBECs gene signature towards decisive determinants of the patient-specific and/or the patient group-specific antiviral response. Finally, the patient-specific mapping of the AID/APOBEC polymorphisms should be considered in the light of COVID-19.

Stepwise evolution and exceptional conservation of ORF1a/b overlap in coronaviruses

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for propagation of coronaviruses. A combination of genomic features that make up PFE-the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements-puts severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess evolutionary dynamics of PFE in great detail. Here we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well defined phylogenetic specificity. We further examined sequencing data from over 350,000 complete genomes and 55,000 raw read datasets to demonstrate exceptional conservation of the PFE region.

Infection of Bronchial Epithelial Cells by the Human Adenoviruses A12, B3, and C2 Differently Regulates the Innate Antiviral Effector APOBEC3B

Human adenoviruses (HAdVs) are a large family of DNA viruses that include more than 100 genotypes divided into seven species (A to G) and induce respiratory tract infections, gastroenteritis, and conjunctivitis. Genetically modified adenoviruses are also used as vaccines, gene therapies, and anticancer treatments. The APOBEC3s are a family of cytidine deaminases that restrict viruses by introducing mutations in their genomes. Viruses developed different strategies to cope with the APOBEC3 selection pressure, but nothing is known on the interplay between the APOBEC3s and the HAdVs. In this study, we focused on three HAdV strains: the B3 and C2 strains, as they are very frequent, and the A12 strain, which is less common but is oncogenic in animal models. We demonstrated that the three HAdV strains induce a similar APOBEC3B upregulation at the transcriptional level. At the protein level, however, APOBEC3B is abundantly expressed during HAdV-A12 and -C2 infection and shows a nuclear distribution. On the contrary, APOBEC3B is barely detectable in HAdV-B3-infected cells. APOBEC3B deaminase activity is detected in total protein extracts upon HAdV-A12 and -C2 infection. Bioinformatic analysis demonstrates that the HAdV-A12 genome bears a stronger APOBEC3 evolutionary footprint than that of the HAdV-C2 and HAdV-B3 genomes. Our results show that HAdV infection triggers the transcriptional upregulation of the antiviral innate effector APOBEC3B. The discrepancies between the APOBEC3B mRNA and protein levels might reflect the ability of some HAdV strains to antagonize the APOBEC3B protein. These findings point toward an involvement of APOBEC3B in HAdV restriction and evolution. IMPORTANCE The APOBEC3 family of cytosine deaminases has important roles in antiviral innate immunity and cancer. Notably, APOBEC3A and APOBEC3B are actively upregulated by several DNA tumor viruses and contribute to transformation by introducing mutations in the cellular genome. Human adenoviruses (HAdVs) are a large family of DNA viruses that cause generally asymptomatic infections in immunocompetent adults. HAdVs encode several oncogenes, and some HAdV strains, like HAdV-A12, induce tumors in hamsters and mice. Here, we show that HAdV infection specifically promotes the expression of the APOBEC3B gene. We report that infection with the A12 strain induces a strong expression of an enzymatically active APOBEC3B protein in bronchial epithelial cells. We provide bioinformatic evidence that HAdVs' genomes and notably the A12 genome are under APOBEC3 selection pressure. Thus, APOBEC3B might contribute to adenoviral restriction, diversification, and oncogenic potential of particular strains.

Read between the lines: Diversity of non-translational selection pressures on local codon usage

Protein coding genes can contain specific motifs within their nucleotide sequence that function as a signal for various biological pathways. The presence of such sequence motifs within a gene can have beneficial or detrimental effects on the phenotype and fitness of an organism, and this can lead to the enrichment or avoidance of this sequence motif. The degeneracy of the genetic code allows for the existence of alternative synonymous sequences that exclude or include these motifs, while keeping the encoded amino acid sequence intact. This implies that locally, there can be a selective pressure for preferentially using a codon over its synonymous alternative in order to avoid or enrich a specific sequence motif. This selective pressure could -in addition to mutation, drift and selection for translation efficiency and accuracy- contribute to shape the codon usage bias. In this review, we discuss patterns of avoidance of (or enrichment for) the various biological signals contained in specific nucleotide sequence motifs: transcription and translation initiation and termination signals, mRNA maturation signals, and antiviral immune system targets. Experimental data on the phenotypic or fitness effects of synonymous mutations in these sequence motifs confirm that they can be targets of local selection pressures on codon usage. We also formulate the hypothesis that transposable elements could have a similar impact on codon usage through their preferred integration sequences. Overall, selection on codon usage appears to be a combination of a global selection pressure imposed by the translation machinery, and a patchwork of local selection pressures related to biological signals contained in specific sequence motifs.

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

The significance of bioengineered nanoplatforms against SARS-CoV-2: From detection to genome editing

COVID-19 outbreak can impose serious negative impacts on the infrastructures of societies including the healthcare systems. Despite the increasing research efforts, false positive or negative results that may be associated with serologic or even RT-PCR tests, inappropriate or variable immune response, and high rates of mutations in coronavirus may negatively affect virus detection process and effectiveness of the vaccines or drugs in development. Nanotechnology-based research attempts via developing state-of-the-art techniques such as nanomechatronics ones and advanced materials including the sensors for detecting the pathogen loads at very low concentrations or site-specific delivery of therapeutics, and real-time protections against the pandemic outbreaks by nanorobots can provide outstanding biomedical breakthroughs. Considering the unique characteristics of pathogens particularly the newly-emerged ones and avoiding the exaggerated optimism or simplistic views on the prophylactic and therapeutic approaches including the one-size-fits-all ones or presenting multiple medications that may be associated with synergistic toxicities rather than enhanced efficiencies might pave the way towards the development of more appropriate treatment strategies with reduced safety concerns. This paper highlights the significance of nanoplatforms against the viral disorders and their capabilities of genome editing that may facilitate taking more appropriate measures against SARS-CoV-2.

Evidence for Strong Mutation Bias toward, and Selection against, U Content in SARS-CoV-2: Implications for Vaccine Design

Large-scale re-engineering of synonymous sites is a promising strategy to generate vaccines either through synthesis of attenuated viruses or via codon-optimized genes in DNA vaccines. Attenuation typically relies on deoptimization of codon pairs and maximization of CpG dinucleotide frequencies. So as to formulate evolutionarily informed attenuation strategies that aim to force nucleotide usage against the direction favored by selection, here, we examine available whole-genome sequences of SARS-CoV-2 to infer patterns of mutation and selection on synonymous sites. Analysis of mutational profiles indicates a strong mutation bias toward U. In turn, analysis of observed synonymous site composition implicates selection against U. Accounting for dinucleotide effects reinforces this conclusion, observed UU content being a quarter of that expected under neutrality. Possible mechanisms of selection against U mutations include selection for higher expression, for high mRNA stability or lower immunogenicity of viral genes. Consistent with gene-specific selection against CpG dinucleotides, we observe systematic differences of CpG content between SARS-CoV-2 genes. We propose an evolutionarily informed approach to attenuation that, unusually, seeks to increase usage of the already most common synonymous codons. Comparable analysis of H1N1 and Ebola finds that GC3 deviated from neutral equilibrium is not a universal feature, cautioning against generalization of results.

Coronavirus genomes carry the signatures of their habitats

Coronaviruses such as SARS-CoV-2 regularly infect host tissues that express antiviral proteins (AVPs) in abundance. Understanding how they evolve to adapt or evade host immune responses is important in the effort to control the spread of infection. Two AVPs that may shape viral genomes are the zinc finger antiviral protein (ZAP) and the apolipoprotein B mRNA editing enzyme-catalytic polypeptide-like 3 (APOBEC3). The former binds to CpG dinucleotides to facilitate the degradation of viral transcripts while the latter frequently deaminates C into U residues which could generate notable viral sequence variations. We tested the hypothesis that both APOBEC3 and ZAP impose selective pressures that shape the genome of an infecting coronavirus. Our investigation considered a comprehensive number of publicly available genomes for seven coronaviruses (SARS-CoV-2, SARS-CoV, and MERS infecting Homo sapiens, Bovine CoV infecting Bos taurus, MHV infecting Mus musculus, HEV infecting Sus scrofa, and CRCoV infecting Canis lupus familiaris). We show that coronaviruses that regularly infect tissues with abundant AVPs have CpG-deficient and U-rich genomes; whereas those that do not infect tissues with abundant AVPs do not share these sequence hallmarks. Among the coronaviruses surveyed herein, CpG is most deficient in SARS-CoV-2 and a temporal analysis showed a marked increase in C to U mutations over four months of SARS-CoV-2 genome evolution. Furthermore, the preferred motifs in which these C to U mutations occur are the same as those subjected to APOBEC3 editing in HIV-1. These results suggest that both ZAP and APOBEC3 shape the SARS-CoV-2 genome: ZAP imposes a strong CpG avoidance, and APOBEC3 constantly edits C to U. Evolutionary pressures exerted by host immune systems onto viral genomes may motivate novel strategies for SARS-CoV-2 vaccine development.

The Role of APOBECs in Viral Replication

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) proteins are a diverse and evolutionarily conserved family of cytidine deaminases that provide a variety of functions from tissue-specific gene expression and immunoglobulin diversity to control of viruses and retrotransposons. APOBEC family expansion has been documented among mammalian species, suggesting a powerful selection for their activity. Enzymes with a duplicated zinc-binding domain often have catalytically active and inactive domains, yet both have antiviral function. Although APOBEC antiviral function was discovered through hypermutation of HIV-1 genomes lacking an active Vif protein, much evidence indicates that APOBECs also inhibit virus replication through mechanisms other than mutagenesis. Multiple steps of the viral replication cycle may be affected, although nucleic acid replication is a primary target. Packaging of APOBECs into virions was first noted with HIV-1, yet is not a prerequisite for viral inhibition. APOBEC antagonism may occur in viral producer and recipient cells. Signatures of APOBEC activity include G-to-A and C-to-T mutations in a particular sequence context. The importance of APOBEC activity for viral inhibition is reflected in the identification of numerous viral factors, including HIV-1 Vif, which are dedicated to antagonism of these deaminases. Such viral antagonists often are only partially successful, leading to APOBEC selection for viral variants that enhance replication or avoid immune elimination.

Footprint of the host restriction factors APOBEC3 on the genome of human viruses

APOBEC3 enzymes are innate immune effectors that introduce mutations into viral genomes. These enzymes are cytidine deaminases which transform cytosine into uracil. They preferentially mutate cytidine preceded by thymidine making the 5'TC motif their favored target. Viruses have evolved different strategies to evade APOBEC3 restriction. Certain viruses actively encode viral proteins antagonizing the APOBEC3s, others passively face the APOBEC3 selection pressure thanks to a depleted genome for APOBEC3-targeted motifs. Hence, the APOBEC3s left on the genome of certain viruses an evolutionary footprint. The aim of our study is the identification of these viruses having a genome shaped by the APOBEC3s. We analyzed the genome of 33,400 human viruses for the depletion of APOBEC3-favored motifs. We demonstrate that the APOBEC3 selection pressure impacts at least 22% of all currently annotated human viral species. The papillomaviridae and polyomaviridae are the most intensively footprinted families; evidencing a selection pressure acting genome-wide and on both strands. Members of the parvoviridae family are differentially targeted in term of both magnitude and localization of the footprint. Interestingly, a massive APOBEC3 footprint is present on both strands of the B19 erythroparvovirus; making this viral genome one of the most cleaned sequences for APOBEC3-favored motifs. We also identified the endemic coronaviridae as significantly footprinted. Interestingly, no such footprint has been detected on the zoonotic MERS-CoV, SARS-CoV-1 and SARS-CoV-2 coronaviruses. In addition to viruses that are footprinted genome-wide, certain viruses are footprinted only on very short sections of their genome. That is the case for the gamma-herpesviridae and adenoviridae where the footprint is localized on the lytic origins of replication. A mild footprint can also be detected on the negative strand of the reverse transcribing HIV-1, HIV-2, HTLV-1 and HBV viruses. Together, our data illustrate the extent of the APOBEC3 selection pressure on the human viruses and identify new putatively APOBEC3-targeted viruses.

Association between APOBEC3H-Mediated Demethylation and Immune Landscape in Head and Neck Squamous Carcinoma

Immunotherapy has been demonstrated as a promising strategy in controlling head and neck squamous cell carcinoma (HNSC). The AID/APOBEC family is well characterized as DNA mutator and considered to play critical roles in immune responses in HNSC. However, the expression pattern and deamination-dependent demethylation roles of AID/APOBECs in HNSC are unclear. In this study, the RNA-seq and DNA methylation profiles of HNSC from TCGA database and cell-based experiments were applied to analyze the relationships between AID/APOBEC expression levels, patients' clinical outcomes, methylation alterations, and immune responses. Here, we found that APOBEC3H was abnormally upregulated in HNSC patients. HPV+ patients tended to have higher APOBEC3H levels than HPV- patients. Remarkably, patients with high APOBEC3H levels showed a favorable overall survival. Furthermore, tumors with high APOBEC3H levels exhibited a genome-wide DNA hypomethylation pattern. APOBEC3H was identified to demethylate and upregulate CXCL10 and improve CD8+ T cell tumor infiltration in the tumor microenvironment. Collectively, APOBEC3H plays critical roles in CD8+ T cell immune infiltration and activation in HNSC, which may be a potential biomarker for oncoimmunotherapy in HNSC.

Is the tryptophan codon of gene vif the Achilles' heel of HIV-1?

To evaluate the impact of hypermutation on the HIV-1 dissemination at the population level we studied 7072 sequences HIV-1 gene vif retrieved from the public databank. From this dataset 854 sequences were selected because they had associated values of CD4+ T lymphocytes counts and viral loads and they were used to assess the correlation between clinical parameters and hypermutation. We found that the frequency of stop codons at sites 5, 11 and 79 ranged from 2.8x10^-4 to 4.2x10^-4. On the other hand, at codons 21, 38, 70, 89 and 174 the frequency of stop codons ranged from 1.4x10^-3 to 2.5x10^-3. We also found a correlation between clinical parameters and hypermutation where patients harboring proviruses with one or more stop codons at the tryptophan sites of the gene vif had higher CD4+ T lymphocytes counts and lower viral loads compared to the population. Our findings indicate that A3 activity potentially restrains HIV-1 replication because individuals with hypermutated proviruses tend to have lower numbers of RNA copies. However, owing to the low frequency of hypermutated sequences observed in the databank (44 out of 7072), it is unlikely that A3 has a significant impact to curb HIV-1 dissemination at the population level.

A-to-I editing of Malacoherpesviridae RNAs supports the antiviral role of ADAR1 in mollusks

Background

Adenosine deaminase enzymes of the ADAR family are conserved in metazoans. They convert adenine into inosine in dsRNAs and thus alter both structural properties and the coding potential of their substrates. Acting on exogenous dsRNAs, ADAR1 exerts a pro- or anti-viral role in vertebrates and Drosophila.

Results

We traced 4 ADAR homologs in 14 lophotrochozoan genomes and we classified them into ADAD, ADAR1 or ADAR2, based on phylogenetic and structural analyses of the enzymatic domain. Using RNA-seq and quantitative real time PCR we demonstrated the upregulation of one ADAR1 homolog in the bivalve Crassostrea gigas and in the gastropod Haliotis diversicolor supertexta during Ostreid herpesvirus-1 or Haliotid herpesvirus-1 infection. Accordingly, we demonstrated an extensive ADAR-mediated editing of viral RNAs. Single nucleotide variation (SNV) profiles obtained by pairing RNA- and DNA-seq data from the viral infected individuals resulted to be mostly compatible with ADAR-mediated A-to-I editing (up to 97%). SNVs occurred at low frequency in genomic hotspots, denoted by the overlapping of viral genes encoded on opposite DNA strands. The SNV sites and their upstream neighbor nucleotide indicated the targeting of selected adenosines. The analysis of viral sequences suggested that, under the pressure of the ADAR editing, the two Malacoherpesviridae genomes have evolved to reduce the number of deamination targets.

Conclusions

We report, for the first time, evidence of an extensive editing of Malacoherpesviridae RNAs attributable to host ADAR1 enzymes. The analysis of base neighbor preferences, structural features and expression profiles of molluscan ADAR1 supports the conservation of the enzyme function among metazoans and further suggested that ADAR1 exerts an antiviral role in mollusks.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1472-6) contains supplementary material, which is available to authorized users.

Suboptimal T-cell Therapy Drives a Tumor Cell Mutator Phenotype That Promotes Escape from First-Line Treatment

Antitumor T-cell responses raised by first-line therapies such as chemotherapy, radiation, tumor cell vaccines, and viroimmunotherapy tend to be weak, both quantitatively (low frequency) and qualitatively (low affinity). We show here that T cells that recognize tumor-associated antigens can directly kill tumor cells if used at high effector-to-target ratios. However, when these tumor-reactive T cells were present at suboptimal ratios, direct T-cell-mediated tumor cell killing was reduced and the ability of tumor cells to evolve away from a coapplied therapy (oncolytic or suicide gene therapy) was promoted. This T-cell-mediated increase in therapeutic resistance was associated with C to T transition mutations that are characteristic of APOBEC3 cytosine deaminase activity and was induced through a TNFα and protein kinase C-dependent pathway. Short hairpin RNA inhibition of endogenous APOBEC3 reduced rates of tumor escape from oncolytic virus or suicide gene therapy to those seen in the absence of antitumor T-cell coculture. Conversely, overexpression of human APOBEC3B in tumor cells enhanced escape from suicide gene therapy and oncolytic virus therapy both in vitro and in vivo Our data suggest that weak affinity or low frequency T-cell responses against tumor antigens may contribute to the ability of tumor cells to evolve away from first-line therapies. We conclude that immunotherapies need to be optimized as early as possible so that, if they do not kill the tumor completely, they do not promote treatment resistance.

A Needle in A Haystack: Tracing Bivalve-Associated Viruses in High-Throughput Transcriptomic Data

Bivalve mollusks thrive in environments rich in microorganisms, such as estuarine and coastal waters, and they tend to accumulate various particles, including viruses. However, the current knowledge on mollusk viruses is mainly centered on few pathogenic viruses, whereas a general view of bivalve-associated viromes is lacking. This study was designed to explore the viral abundance and diversity in bivalve mollusks using transcriptomic datasets. From analyzing RNA-seq data of 58 bivalve species, we have reconstructed 26 nearly complete and over 413 partial RNA virus genomes. Although 96.4% of the predicted viral proteins refer to new viruses, some sequences belong to viruses associated with bivalve species or other marine invertebrates. We considered short non-coding RNAs (sncRNA) and post-transcriptional modifications occurring specifically on viral RNAs as tools for virus host-assignment. We could not identify virus-derived small RNAs in sncRNA reads obtained from the oyster sample richest in viral reads. Single Nucleotide Polymorphism (SNP) analysis revealed 938 A-to-G substitutions occurring on the 26 identified RNA viruses, preferentially impacting the AA di-nucleotide motif. Under-representation analysis revealed that the AA motif is under-represented in these bivalve-associated viruses. These findings improve our understanding of bivalve viromes, and set the stage for targeted investigations on the specificity and dynamics of identified viruses.

Evolutionary effects of the AID/APOBEC family of mutagenic enzymes on human gamma-herpesviruses

The human gamma-herpesviruses, Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, establish lifelong latency in B cells and are associated with multiple malignancies. Virus-host coevolution often drive changes in both host immunity and in the viral genome. We consider one host immune mechanism, the activation-induced deaminase (AID)/APOBEC family of cytidine deaminases, that induces mutations in viral DNA. AID, the ancestral gene in the family has a conserved role in somatic hypermutation, a key step in antibody affinity maturation. The APOBEC3 subfamily, of which there are seven genes in human, have evolved antiviral functions and have diversified in terms of their expression pattern, subcellular localization, and DNA mutation motifs (hotspots). In this study, we investigated how the human gamma-herpesviruses have evolved to avoid the action of the AID/APOBEC enzymes and determine if these enzymes are contributing to the ongoing evolution of the viruses. We used computational methods to evaluate observed versus expected frequency of AID/APOBEC hotspots in viral genomes and found that the viruses have evolved to limit the representation of AID and certain APOBEC3 motifs. At the same time, the remaining hotspots were highly likely to cause amino acid changes, suggesting prolonged evolutionary pressure of the enzymes on the viruses. To study current hypermutation, as opposed to historical mutation processes, we also analyzed putative mutations derived from alignments of published viral genomes and found again that AID and APOBEC3 appear to target the genome most frequently. New protein variants resulting from AID/APOBEC activity may have important consequences in health, including vaccine development (epitope evolution) and host immune evasion.

The cytidine deaminase under-representation reporter (CDUR) as a tool to study evolution of sequences under deaminase mutational pressure

Background

Activation induced deaminase (AID) and apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3 (APOBEC3) are deaminases that mutate C to U on single-stranded DNA (ssDNA). AID is expressed primarily in germinal center B-cells, where it facilitates affinity maturation and class-switch recombination. APOBEC3 are a family of anti-viral proteins that act as part of the intrinsic immune response. In both cases, there are particular sequence motifs, also known as “mutation motifs”, to which these deaminases prefer to bind and mutate.

Results

We present a program, the cytidine deaminase under-representation reporter (CDUR) designed to statistically determine whether a given sequence has an under/over-representation of these mutation motifs. CDUR shows consitency with other studies of mutation motifs, as we show by analyzing sequences from the adeno-associated virus 2 (AAV2) and human papillomavirus (HPV).

Conclusion

Using various shuffling mechanisms to generate different null model distributions, we can tailor CDUR to correct for metrics such as GC-content, dinucleotide frequency, and codon bias.

Molecular cloning and characterization of APOBEC3 family in tree shrew

The APOBEC3 family is a series antiviral factors that inhibit the replication of many viruses, such as HIV-1 and HBV. Tree shrews (Tupaia belangeri) possess great potential as an animal model for human diseases and therapeutic responses. However, the APOBEC3 family is unknown in tree shrews. Recent work has showed the presence of the APOBEC3 family in tree shrews. In this work, the cDNA sequences of five APOBEC3 members were identified in tree shrews, namely, tsAPOBEC3A, -3C, -3F, -3G and -3H. The results showed that their sequences encoded a zinc (Z)-coordinating-domain as a characteristic of APOBEC3 proteins. Phylogenetic analysis revealed that the tree shrew APOBEC3 (tsAPOBEC3) genes have occurred independently and that they are clustered with other mammalian APOBEC3 members. Transcript expression analysis indicated that tsAPOBEC3 genes are constitutively expressed, and high in immune-related tissues. tsAPOBEC3 gene expression was up-regulated in hepatocytes and PBMCs by IFN-α stimulation. Finally, tsAPOBEC3 proteins could edit both sides of DNA by inserting G→A and C→T hypermutations. Overall, the results suggest that the tsAPOBEC3 family could play a key role in defense immunity through distinct editing mechanisms. Our results provided insights into the genetic basis for the development of a tree shrew model for studying viral infection. Future studies will focus on deepening our understanding on the antiviral functions of these editing enzymes in tree shrew.

The CpG dinucleotide content of the HIV-1 envelope gene may predict disease progression

The clinical course of HIV-1 varies greatly among infected individuals. Despite extensive research, virus factors associated with slow-progression remain poorly understood. Identification of unique HIV-1 genomic signatures linked to slow-progression remains elusive. We investigated CpG dinucleotide content in HIV-1 envelope gene as a potential virus factor in disease progression. We analysed 1808 HIV-1 envelope gene sequences from three independent longitudinal studies; this included 1280 sequences from twelve typical-progressors and 528 sequences from six slow-progressors. Relative abundance of CpG dinucleotides and relative synonymous codon usage (RSCU) for CpG-containing codons among HIV-1 envelope gene sequences from typical-progressors and slow-progressors were analysed. HIV-1 envelope gene sequences from slow-progressors have high-CpG dinucleotide content and increased number of CpG-containing codons as compared to typical-progressors. Our findings suggest that observed differences in CpG-content between typical-progressors and slow-progressors is not explained by differences in the mononucleotide content. Our results also highlight that the high-CpG content in HIV-1 envelope gene from slow-progressors is observed immediately after seroconversion. Thus CpG dinucleotide content of HIV-1 envelope gene is a potential virus-related factor that is linked to disease progression. The CpG dinucleotide content of HIV-1 envelope gene may help predict HIV-1 disease progression at early stages after seroconversion.