Footprint of APOBEC3 on the Genome of Human Retroelements

Almost half of the human genome is composed of transposable elements. The genomic structures and life cycles of some of these elements suggest they are a result of waves of retroviral infection and transposition over millions of years. The reduction of retrotransposition activity in primates compared to that in nonprimates, such as mice, has been attributed to the positive selection of several antiretroviral factors, such as apolipoprotein B mRNA editing enzymes. Among these, APOBEC3G is known to mutate G to A within the context of GG in the genome of endogenous as well as several exogenous retroelements (the underlining marks the G that is mutated). On the other hand, APOBEC3F and to a lesser extent other APOBEC3 members induce G-to-A changes within the nucleotide GA. It is known that these enzymes can induce deleterious mutations in the genome of retroviral sequences, but the evolution and/or inactivation of retroelements as a result of mutation by these proteins is not clear. Here, we analyze the mutation signatures of these proteins on large populations of long interspersed nuclear element (LINE), short interspersed nuclear element (SINE), and endogenous retrovirus (ERV) families in the human genome to infer possible evolutionary pressure and/or hypermutation events. Sequence context dependency of mutation by APOBEC3 allows investigation of the changes in the genome of retroelements by inspecting the depletion of G and enrichment of A within the APOBEC3 target and product motifs, respectively. Analysis of approximately 22,000 LINE-1 (L1), 24,000 SINE Alu, and 3,000 ERV sequences showed a footprint of GG→AG mutation by APOBEC3G and GA→AA mutation by other members of the APOBEC3 family (e.g., APOBEC3F) on the genome of ERV-K and ERV-1 elements but not on those of ERV-L, LINE, or SINE.     

APOBEC3 has not left an evolutionary footprint on the HIV-1 genome

It is known that the human immune proteins APOBEC3G and -F (hA3G/F) can inhibit Vif-deficient HIV by G-to-A mutation; however, the roles of these enzymes in the evolution of HIV are debated. We argue that if evolutionary pressure from hA3G/F exists there should be evidence of their imprint on the HIV genome in the form of (i) underrepresentation of hA3G/F target motifs (e.g., TGGG [targeted position is underlined]) and overrepresentation of product motifs (e.g., TAGG) and/or (ii) an increase in the ratio of nonsynonymous to synonymous (NS/S) G-to-A changes among hA3G/F target motifs and a decrease of NS/S A-to-G changes among hA3G/F product motifs. To test the first hypothesis, we studied the representation of hA3G/F target and product motifs in 1,932 complete HIV-1 genomes using Markov models. We found that the highly targeted motifs are not underrepresented and their product motifs are not overrepresented. To test the second hypothesis, we determined the NS/S G↔A changes among the hA3G/F target and product motifs in 1,540 complete sets of nine HIV-1 genes. The NS/S changes did not show an increasing/decreasing trend within the target/product motifs, but the NS/S changes within the motif AG was exceptionally low. We observed the same pattern by analyzing 740 human genes. Given that hA3G/F do not act on the human genome, this suggests a small NS/S change within AG has arisen by other mechanisms. We therefore find no evidence of an evolutionary footprint of hA3G/F. We postulate several mechanisms to explain why the HIV-1 genome does not contain the hA3G/F footprint.     

Twin gradients in APOBEC3 edited HIV-1 DNA reflect the dynamics of lentiviral replication

The human immunodeficiency virus (HIV) Vif protein blocks incorporation of two host cell cytidine deaminases, APOBEC3F and 3G, into the budding virion. Not surprisingly, on a vif background nascent minus strand DNA can be extensively edited leaving multiple uracil residues. Editing occurs preferentially in the context of TC (GA on the plus strand) and CC (GG) depending on the enzyme. To explore the distribution of APOBEC3F and –3G editing across the genome, a product/substrate ratio (AA + AG)/(GA + GG) was computed for a series of 30 edited genomes present in the data bases. Two highly polarized gradients were noted each with maxima just 5′ to the central polypurine tract (cPPT) and LTR proximal polypurine tract (3′PPT). The gradients are in remarkable agreement with the time the minus strand DNA remains single stranded. In vitro analyses of APOBEC3G deamination of nascent cDNA spanning the two PPTs showed no pronounced dependence on the PPT RNA:DNA heteroduplex ruling out the competing hypothesis of a PPT orientation effect. The degree of hypermutation varied smoothly among genomes indicating that the number of APOBEC3 molecules packaged varied considerably.     

Positioning of APOBEC3G/F Mutational Hotspots in the Human Immunodeficiency Virus Genome Favors Reduced Recognition by CD8+ T Cells

Due to constitutive expression in cells targeted by human immunodeficiency virus (HIV), and immediate mode of viral restriction upon HIV entry into the host cell, APOBEC3G (A3G) and APOBEC3F (A3F) have been considered primarily as agents of innate immunity. Recent bioinformatic and mouse model studies hint at the possibility that mutation of the HIV genome by these enzymes may also affect adaptive immunity but whether this occurs in HIV-infected individuals has not been examined. We evaluated whether APOBEC-mediated mutations within common HIV CD8⁺ T cell epitopes can potentially enhance or diminish activation of HIV-specific CD8⁺ T cells from infected individuals. We compared ex vivo activation of CD8⁺ T lymphocytes from HIV-infected individuals by wild type HIV peptide epitopes and synthetic variants bearing simulated A3G/F-induced mutations by measuring interferon-γ (IFN-γ) production. We found that A3G/F-induced mutations consistently diminished HIV-specific CD8⁺ T cell responses against the common epitopes we tested. If this reflects a significant trend in vivo, then adaptation by HIV to enrich sequences that are favored for mutation by A3G/F (A3G/F hotspots) in portions of its genome that encode immunogenic CD8⁺ T cell epitopes would favor CTL escape. Indeed, we found the most frequently mutated A3G motif (CCC) is enriched up to 6-fold within viral genomic sequences encoding immunodominant CD8⁺ T cell epitopes in Gag, Pol and Nef. Within each gene, A3G/F hotspots are more abundant in sequences encoding epitopes that are commonly recognized due to their HLA restriction. Thus, in our system, mutations of the HIV genome, mimicking A3G/F activity, appeared to abrogate or severely reduce CTL recognition. We suggest that the physiological significance of this potential effect in facilitating CTL escape is echoed in the adaptation of the HIV genome to enrich A3G/F hotspots in sequences encoding CTL epitopes that are more immunogenic at the population level.     

APOBEC3G-induced hypermutation of human immunodeficiency virus type-1 is typically a discrete "all or nothing" phenomenon

The rapid evolution of Human Immunodeficiency Virus (HIV-1) allows studies of ongoing host-pathogen interactions. One key selective host factor is APOBEC3G (hA3G) that can cause extensive and inactivating Guanosine-to-Adenosine (G-to-A) mutation on HIV plus-strand DNA (termed hypermutation). HIV can inhibit this innate anti-viral defense through binding of the viral protein Vif to hA3G, but binding efficiency varies and hypermutation frequencies fluctuate in patients. A pivotal question is whether hA3G-induced G-to-A mutation is always lethal to the virus or if it may occur at sub-lethal frequencies that could increase viral diversification. We show in vitro that limiting-levels of hA3G-activity (i.e. when only a single hA3G-unit is likely to act on HIV) produce hypermutation frequencies similar to those in patients and demonstrate in silico that potentially non-lethal G-to-A mutation rates are ～10-fold lower than the lowest observed hypermutation levels in vitro and in vivo. Our results suggest that even a single incorporated hA3G-unit is likely to cause extensive and inactivating levels of HIV hypermutation and that hypermutation therefore is typically a discrete "all or nothing" phenomenon. Thus, therapeutic measures that inhibit the interaction between Vif and hA3G will likely not increase virus diversification but expand the fraction of hypermutated proviruses within the infected host.     

The nucleotide targets of somatic mutation and the role of selection in immunoglobulin heavy chains of a teleost fish

Sequence analysis of H chain cDNA derived from the spleen of an individual catfish has shown that somatic mutation occurs within both the VH- and JH-encoded regions. Somatic mutation preferentially targets G and C nucleotides with approximately balanced frequencies, resulting in the predominant accumulation of G-to-A and C-to-T substitutions that parallel the activation-induced cytidine deaminase nucleotide exchanges known in mammals. The overall mutation rate of A nucleotides is not significantly different from that expected by sequence-insensitive mutations, and a significant bias exists against mutations occurring in T. Targeting of mutations is dependent upon the sequence of neighboring nucleotides, allowing statistically significant hotspot motifs to be identified. Dinucleotide, trinucleotide, and RGYW analyses showed that mutational targets in catfish are restricted when compared with the spectrum of targets known in mammals. The preferential targets for G and C mutation are the central GC positions in both AGCT and AGCA. The WA motif, recognized as a mammalian hotspot for A mutations, was not a significant target for catfish mutations. The only significant target for A mutations was the terminal position in AGCA. Lastly, comparisons of mutations located in framework region and CDR codons coupled with multinomial distribution studies found no substantial evidence in either independent or clonally related VDJ rearrangements to indicate that somatic mutation coevolved with mechanisms that select B cells based upon nonsynonymous mutations within CDR-encoded regions. These results suggest that the principal role of somatic mutation early in phylogeny was to diversify the repertoire by targeting hotspot motifs preferentially located within CDR-encoded regions.     

Role of APOBEC3 in genetic diversity among endogenous murine leukemia viruses

The ability of human and murine APOBECs (specifically, APOBEC3) to inhibit infecting retroviruses and retrotransposition of some mobile elements is becoming established. Less clear is the effect that they have had on the establishment of the endogenous proviruses resident in the human and mouse genomes. We used the mouse genome sequence to study diversity and genetic traits of nonecotropic murine leukemia viruses (polytropic [Pmv], modified polytropic [Mpmv], and xenotropic [Xmv] subgroups), the best-characterized large set of recently integrated proviruses. We identified 49 proviruses. In phylogenetic analyses, Pmvs and Mpmvs were monophyletic, whereas Xmvs were divided into several clades, implying a greater number of replication cycles between the integration events. Four distinct primer binding site types (Pro, Gln1, Gln2 and Thr) were dispersed within the phylogeny, indicating frequent mispriming. We analyzed the frequency and context of G-to-A mutations for the role of mA3 in formation of these proviruses. In the Pmv and Mpmv (but not Xmv) groups, mutations attributable to mA3 constituted a large fraction of the total. A significant number of nonsense mutations suggests the absence of purifying selection following mutation. A strong bias of G-to-A relative to C-to-T changes was seen, implying a strand specificity that can only have occurred prior to integration. The optimal sequence context of G-to-A mutations, TTC, was consistent with mA3. At least in the Pmv group, a significant 5′ to 3′ gradient of G-to-A mutations was consistent with mA3 editing. Altogether, our results for the first time suggest mA3 editing immediately preceding the integration event that led to retroviral endogenization, contributing to inactivation of infectivity.     

Conserved footprints of APOBEC3G on Hypermutated human immunodeficiency virus type 1 and human endogenous retrovirus HERV-K(HML2) sequences

The human polynucleotide cytidine deaminases APOBEC3G (hA3G) and APOBEC3F (hA3F) are antiviral restriction factors capable of inducing extensive plus-strand guanine-to-adenine (G-to-A) hypermutation in a variety of retroviruses and retroelements, including human immunodeficiency virus type 1 (HIV-1). They differ in target specificity, favoring plus-strand 5'GG and 5'GA dinucleotide motifs, respectively. To characterize their mutational preferences in detail, we analyzed single-copy, near-full-length HIV-1 proviruses which had been hypermutated in vitro by hA3G or hA3F. hA3-induced G-to-A mutation rates were significantly influenced by the wider sequence context of the target G. Moreover, hA3G, and to a lesser extent hA3F, displayed clear tetranucleotide preference hierarchies, irrespective of the genomic region examined and overall hypermutation rate. We similarly analyzed patient-derived hypermutated HIV-1 genomes using a new method for estimating reference sequences. The majority of these, regardless of subtype, carried signatures of hypermutation that strongly correlated with those induced in vitro by hA3G. Analysis of genome-wide hA3-induced mutational profiles confirmed that hypermutation levels were reduced downstream of the polypurine tracts. Additionally, while hA3G mutations were found throughout the genome, hA3F often intensely mutated shorter regions, the locations of which varied between proviruses. We extended our analysis to human endogenous retroviruses (HERVs) from the HERV-K(HML2) family, finding two elements that carried clear footprints of hA3G activity. This constitutes the most direct evidence to date for hA3G activity in the context of natural HERV infections, demonstrating the involvement of this restriction factor in defense against retroviral attacks over millions of years of human evolution.     

Human APOBEC3 Induced Mutation of Human Immunodeficiency Virus Type-1 Contributes to Adaptation and Evolution in Natural Infection

Human APOBEC3 proteins are cytidine deaminases that contribute broadly to innate immunity through the control of exogenous retrovirus replication and endogenous retroelement retrotransposition. As an intrinsic antiretroviral defense mechanism, APOBEC3 proteins induce extensive guanosine-to-adenosine (G-to-A) mutagenesis and inhibit synthesis of nascent human immunodeficiency virus-type 1 (HIV-1) cDNA. Human APOBEC3 proteins have additionally been proposed to induce infrequent, potentially non-lethal G-to-A mutations that make subtle contributions to sequence diversification of the viral genome and adaptation though acquisition of beneficial mutations. Using single-cycle HIV-1 infections in culture and highly parallel DNA sequencing, we defined trinucleotide contexts of the edited sites for APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H. We then compared these APOBEC3 editing contexts with the patterns of G-to-A mutations in HIV-1 DNA in cells obtained sequentially from ten patients with primary HIV-1 infection. Viral substitutions were highest in the preferred trinucleotide contexts of the edited sites for the APOBEC3 deaminases. Consistent with the effects of immune selection, amino acid changes accumulated at the APOBEC3 editing contexts located within human leukocyte antigen (HLA)-appropriate epitopes that are known or predicted to enable peptide binding. Thus, APOBEC3 activity may induce mutations that influence the genetic diversity and adaptation of the HIV-1 population in natural infection.     

APOBEC3D and APOBEC3F Potently Promote HIV-1 Diversification and Evolution in Humanized Mouse Model

Several APOBEC3 proteins, particularly APOBEC3D, APOBEC3F, and APOBEC3G, induce G-to-A hypermutations in HIV-1 genome, and abrogate viral replication in experimental systems, but their relative contributions to controlling viral replication and viral genetic variation in vivo have not been elucidated. On the other hand, an HIV-1-encoded protein, Vif, can degrade these APOBEC3 proteins via a ubiquitin/proteasome pathway. Although APOBEC3 proteins have been widely considered as potent restriction factors against HIV-1, it remains unclear which endogenous APOBEC3 protein(s) affect HIV-1 propagation in vivo. Here we use a humanized mouse model and HIV-1 with mutations in Vif motifs that are responsible for specific APOBEC3 interactions, DRMR/AAAA (4A) or YRHHY/AAAAA (5A), and demonstrate that endogenous APOBEC3D/F and APOBEC3G exert strong anti-HIV-1 activity in vivo. We also show that the growth kinetics of 4A HIV-1 negatively correlated with the expression level of APOBEC3F. Moreover, single genome sequencing analyses of viral RNA in plasma of infected mice reveal that 4A HIV-1 is specifically and significantly diversified. Furthermore, a mutated virus that is capable of using both CCR5 and CXCR4 as entry coreceptor is specifically detected in 4A HIV-1-infected mice. Taken together, our results demonstrate that APOBEC3D/F and APOBEC3G fundamentally work as restriction factors against HIV-1 in vivo, but at the same time, that APOBEC3D and APOBEC3F are capable of promoting viral diversification and evolution in vivo.     

Comparison of the differential context-dependence of DNA deamination by APOBEC enzymes: correlation with mutation spectra in vivo

To investigate the extent to which in vivo mutation spectra might reflect the intrinsic specificities of active mutators, genetic and biochemical assays were used to analyse the DNA target specificities of cytidine deaminases of the APOBEC family. The results reveal the critical importance of nucleotides immediately 5' of the targeted C for the specificity of all three enzymes studied (AID, APOBEC1 and APOBEC3G). At position -1, APOBEC1 showed a marked preference for dT, AID for dA/dG and APOBEC3G a strong preference for dC. Furthermore, AID and APOBEC3G showed distinct dependence on the nucleotide at position -2 with dA/dT being favoured by AID and dC by APOBEC3G. Most if not all activity of the recombinant deaminases on free dC could be attributed to low-level contamination by host enzymes. The target preference of APOBEC3G supports it being a major but possibly not sole contributor to HIV hypermutation without making it a dominant contribution to general HIV sequence variation. The specificity of AID as deduced from the genetic assay (which relies on inactivation of sacB of Bacillus subtilis) agrees well with that deduced by Pham et al. using an in vitro assay although we postulate that major intrinsic mutational hotspots in immunoglobulin V genes in vivo might reflect favoured sites of AID action being generated by proximal DNA targets located on opposite DNA strands. The target specificity of AID also accords with the spectrum of mutations observed in B lymphoma-associated oncogenes. The possibility of deaminase involvement in non-lymphoid human tumours is hinted at by tissue-specific differences in the spectra of dC transitions in tumour-suppressor genes. Thus, the patterns of hypermutation in antibodies and retroviruses owe much to the intrinsic sequence preferences of the AID/APOBEC family of DNA deaminases: analogous biases might also contribute to the spectra of cancer-associated mutation.     

APOBEC1 cytosine deaminase activity on single-stranded DNA is suppressed by replication protein A

Many APOBEC cytidine deaminase members are known to induce ‘off-target’ cytidine deaminations in 5′TC motifs in genomic DNA that contribute to cancer evolution. In this report, we characterized APOBEC1, which is a possible cancer related APOBEC since APOBEC1 mRNA is highly expressed in certain types of tumors, such as lung adenocarcinoma. We found a low level of APOBEC1-induced DNA damage, as measured by γH2AX foci, in genomic DNA of a lung cancer cell line that correlated to its inability to compete in vitro with replication protein A (RPA) for ssDNA. This suggests that RPA can act as a defense against off-target deamination for some APOBEC enzymes. Overall, the data support the model that the ability of an APOBEC to compete with RPA can better predict genomic damage than combined analysis of mRNA expression levels in tumors and analysis of mutation signatures.     

Transmissible Gastroenteritis Coronavirus Genome Packaging Signal Is Located at the 5′ End of the Genome and Promotes Viral RNA Incorporation into Virions in a Replication-Independent Process

Preferential RNA packaging in coronaviruses involves the recognition of viral genomic RNA, a crucial process for viral particle morphogenesis mediated by RNA-specific sequences, known as packaging signals. An essential packaging signal component of transmissible gastroenteritis coronavirus (TGEV) has been further delimited to the first 598 nucleotides (nt) from the 5′ end of its RNA genome, by using recombinant viruses transcribing subgenomic mRNA that included potential packaging signals. The integrity of the entire sequence domain was necessary because deletion of any of the five structural motifs defined within this region abrogated specific packaging of this viral RNA. One of these RNA motifs was the stem-loop SL5, a highly conserved motif in coronaviruses located at nucleotide positions 106 to 136. Partial deletion or point mutations within this motif also abrogated packaging. Using TGEV-derived defective minigenomes replicated in trans by a helper virus, we have shown that TGEV RNA packaging is a replication-independent process. Furthermore, the last 494 nt of the genomic 3′ end were not essential for packaging, although this region increased packaging efficiency. TGEV RNA sequences identified as necessary for viral genome packaging were not sufficient to direct packaging of a heterologous sequence derived from the green fluorescent protein gene. These results indicated that TGEV genome packaging is a complex process involving many factors in addition to the identified RNA packaging signal. The identification of well-defined RNA motifs within the TGEV RNA genome that are essential for packaging will be useful for designing packaging-deficient biosafe coronavirus-derived vectors and providing new targets for antiviral therapies.     

Evolution of HIV-1 isolates that use a novel Vif-independent mechanism to resist restriction by human APOBEC3G

The human APOBEC3G protein restricts the replication of Vif-deficient HIV-1 by deaminating nascent viral cDNA cytosines to uracils, leading to viral genomic strand G-to-A hypermutations. However, the HIV-1 Vif protein triggers APOBEC3G degradation, which helps to explain why this innate defense does not protect patients. The APOBEC3G-Vif interaction is a promising therapeutic target, but the benefit of the enabling of HIV-1 restriction in patients is unlikely to be known until Vif antagonists are developed. As a necessary prelude to such studies, cell-based HIV-1 evolution experiments were done to find out whether APOBEC3G can provide a long-term block to Vif-deficient virus replication and, if so, whether HIV-1 variants that resist restriction would emerge. APOBEC3G-expressing T cells were infected with Vif-deficient HIV-1. Virus infectivity was suppressed in 45/48 cultures for more than five weeks, but replication was eventually detected in three cultures. Virus-growth characteristics and sequencing demonstrated that these isolates were still Vif-deficient and that in fact, these viruses had acquired a promoter mutation and a Vpr null mutation. Resistance occurred by a novel tolerance mechanism in which the resistant viruses packaged less APOBEC3G and accumulated fewer hypermutations. These data support the development of antiretrovirals that antagonize Vif and thereby enable endogenous APOBEC3G to suppress HIV-1 replication.     

Natural variation in Vif: differential impact on APOBEC3G/3F and a potential role in HIV-1 diversification

The HIV-1 Vif protein counteracts the antiviral activity exhibited by the host cytidine deaminases APOBEC3G and APOBEC3F. Here, we show that defective vif alleles can readily be found in HIV-1 isolates and infected patients. Single residue changes in the Vif protein sequence are sufficient to cause the loss of Vif-induced APOBEC3 neutralization. Interestingly, not all the detected defects lead to a complete inactivation of Vif function since some mutants retained selective neutralizing activity against APOBEC3F but not APOBEC3G or vice versa. Concordantly, independently hypermutated proviruses with distinguishable patterns of G-to-A substitution attributable to cytidine deamination induced by APOBEC3G, APOBEC3F, or both enzymes were present in individuals carrying proviruses with completely or partly defective Vif variants. Natural variation in Vif function may result in selective and partial neutralization of cytidine deaminases and thereby promote viral sequence diversification within HIV-1 infected individuals.