# Evolution of the primate APOBEC3A cytidine deaminase gene and identification of related coding regions

Michel Henry, Christophe Terzian, Martine Peeters, Simon Wain-Hobson, Jean Pierre Vartanian

► To cite this version:

Michel Henry, Christophe Terzian, Martine Peeters, Simon Wain-Hobson, Jean Pierre Vartanian. Evolution of the primate APOBEC3A cytidine deaminase gene and identification of related coding regions. PLoS ONE, 2012, 7 (1), pp.e30036. 10.1371/journal.pone.0030036. [hal-02643187](https://hal.inrae.fr/hal-02643187)

HAL Id: hal-02643187

https://hal.inrae.fr/hal-02643187

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Evolution of the Primate APOBEC3A Cytidine Deaminase Gene and Identification of Related Coding Regions

Michel Henry¹, Christophe Terzian², Martine Peeters³, Simon Wain-Hobson¹, Jean-Pierre Vartanian¹*

1 Molecular Retrovirology Unit, Institut Pasteur, Paris, France, 2 UMR754, UCBL-INRA and Ecole Pratique des Hautes Etudes, Lyon, France, 3 Institut de Recherche pour le Developpement (IRD) and University of Montpellier 1, UMI 233, Montpellier, France

## Abstract

The APOBEC3 gene cluster encodes six cytidine deaminases (A3A-C, A3DE, A3F-H) with single stranded DNA (ssDNA) substrate specificity. For the moment A3A is the only enzyme that can initiate catabolism of both mitochondrial and nuclear DNA. Human A3A expression is initiated from two different methionine codons M1 or M13, both of which are in adequate but sub-optimal Kozak environments. In the present study, we have analyzed the genetic diversity among A3A genes across a wide range of 12 primates including New World monkeys, Old World monkeys and Hominids. Sequence variation was observed in exons 1-4 in all primates with up to 31% overall amino acid variation. Importantly for 3 hominids codon M1 was mutated to a threonine codon or valine codon, while for 5/12 primates strong Kozak M1 or M13 codons were found. Positive selection was apparent along a few branches which differed compared to positive selection in the carboxy-terminal of A3G that clusters with A3A among human cytidine deaminases. In the course of analyses, two novel non-functional A3A-related fragments were identified on chromosome 4 and 8 kb upstream of the A3 locus. This qualitative and quantitative variation among primate A3A genes suggest that subtle differences in function might ensue as more light is shed on this increasingly important enzyme.

**Citation:** Henry M, Terzian C, Peeters M, Wain-Hobson S, Vartanian J-P (2012) Evolution of the Primate APOBEC3A Cytidine Deaminase Gene and Identification of Related Coding Regions. PLoS ONE 7(1): e30036. doi:10.1371/journal.pone.0030036

**Editor:** Peter Sommer, Institut Pasteur Korea, Korea, Republic of Korea

**Received** October 14, 2011; **Accepted** December 8, 2011; **Published** January 17, 2012

**Copyright:** © 2012 Henry et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** Funding was received from the Institut Pasteur, Agence Nationale de Recherche sur le SIDA (ANRS), Agence Nationale de Recherche (ANR), l'Institut de Recherche pour le Développement (IRD), l'Institut National de la santé et de la Recherche (INSERM) and the Centre National de Recherche Scientifique (CNRS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: jean-pierre.vartanian@pasteur.fr

## Introduction

The APOBEC3 seven gene cluster (A3A-C, A3DE, A3F-H) encodes six cytidine deaminases with single stranded DNA (ssDNA) substrate specificity [1,2,3,4,5,6,7,8,9]. Several are clearly innate restriction factors for viruses, notably for retroviruses, hepadnaviruses or parvoviruses [3,5,6,10,11,12,13,14,15,16,17,18,19,20,21,22]. A3G and A3F constituted such a strong barrier for the lentiviral group of retroviruses that all but one encode a *vif* gene whose protein (Vif) is a powerful antagonist [3,5,6,7,19,23,24,25,26]. Hepatitis B virus (HBV) is restricted by at least two A3 enzymes while herpes simplex virus type 1 is restricted by A3C [20,22,27]. To date there are no reports of A3 antagonists encoded by these viral genomes. This antiviral role fits with the repeated observation that several A3 genes are up-regulated by type I and II interferons [10,28,29,30,31]. However, recent work has shown that this antiviral role is just part of a bigger picture [30,32,33]. For example, A3A can restrict Line transposition [11,34,35,36]. Several A3 enzymes can initiate catabolism of mitochondrial DNA, in which uracil N-glycosylase plays a major role downstream of editing [33]. For the moment A3A is the only enzyme that can initiate catabolism of both mitochondrial and nuclear DNA [33].

These A3 proteins mediate hydrolytic deamination at the C4 position that oxidises cytosine to uracil in ssDNA so generating C→U hyper-edited molecules [1,3,5,6,7,9,22,37,38]. The active sites of A3 enzymes are characterized by a conserved zinc-finger HAEX23-28PCX2-4C motif [4]. These A3 enzymes show a strong

preference for cytidine deamination occurring segment carboxy-terminal to the zinc finger impacts this dinucleotide specificity [39]. Human A3A expression is initiated from two different methionine codons (M1 and M13), both of which are in adequate but sub-optimal Kozak environments [40].

Even though a number of primate genomes are available, only the chimpanzee locus is colinear. For the orang-utan the *A3A* gene is incomplete while the entire locus contains 12 exon 3/exon 6 domains rather than 11. The *A3A* gene is missing in the Rhesus macaque assembly, while the marmoset locus doesn't exist *per se*, sequences being distributed over numerous contigs. As the *A3* locus shows signs of extensive gene conversion, the apparent gaps might reflect assembly problems.

We have analyzed the genetic diversity among *A3A* genes across a wide range of primates including New World monkeys, Old World monkeys and Hominids. There is variation among the Kozak motifs with the M1 initiator methionine being absent for chimpanzees, bonobos and gorillas. Some, but not all, *A3A* lineages show positive selection suggesting that A3A enzymes may not be truly orthologous.

## Results

### Primate A3A cytidine deaminases

Twelve primates A3A sequences spanning New and Old World monkeys were derived by amplification of genomic DNA and given aligned to the human sequence (Figure 1). The A3A protein

<figure><img src="image_2.png" /></figure>
is initiated at codons M1 or M13 giving rise to two different proteins both with ssDNA cytidine deaminase activity [40]. The Kozak context of both human A3A initiator codons is considered to be adequate. For 3 hominids, codon M1 was mutated to a threonine codon or valine codon which probably abrogates translation initiation (Figure 1, Table 1). For both New World monkey sequences, the M1 Kozak context was strong suggesting that translation initiation at M13 would be reduced. In addition, the Kozak context of the M13 codon was strong for 3/12 primates notably *C. guereza*, *C. aethiops* and *C. neglectus* (Table 1). For all the others, the context is considered to be adequate for translation initiation.

Sequence variation was observed in all exons apart from the very small exon 5. Some of the exon 5 sequences differ compared to some recently reported [35]. On a pairwise basis up to 31% amino acid divergence was observed overall, with 6%, 21% and

30% among hominids, Old World small monkeys and New World monkeys respectively. That the variation is as great overall as that between the New World monkeys, suggests that there has not been too much gene conversion in the New World lineage. Exon 3 encodes the hallmark HXEX<sub>23-28</sub>PCX<sub>2-4</sub>C motif for cytidine deaminases (Figure 1). Among all the human A3 enzymes only A3A encodes the PCX<sub>4</sub>C variant. Interestingly, the New World A3A sequences are singular in that they encode the PCX<sub>2</sub>C variant typical of all other A3 enzymes.

## A3A is under positive selection in Old World monkeys

In order to characterise whether this variation shows signs of selection, we estimated the relative numbers of non-synonymous (dN) and synonymous (dS) nucleotide substitutions per site and dN/dS ratios over the twelve primate species using the Hyphy

<figure><img src="image_3.png" /><figcaption>**Figure 1. Alignment of primate APOBEC3A proteins.** Twelve primate sequences were compared to *Homo sapiens* used as reference. Only differences are shown. Hyphens denote gaps introduced to maximize sequence identity. The numbering corresponds to that of the human sequence. The letters a, b, c are added to adjacent residue to accommodate insertions. Red denotes the first (M1, exon 1) and second initiation start codons (M13, exon 2). The crucial cytidine deaminase motif residues are highlighted in magenta. Positively and negatively selected codon sites are in blue and green respectively. The predicted secondary structure motifs for hA3A are underlined.
doi:10.1371/journal.pone.0030036.g001</figcaption></figure>

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1.** Kozak sequence contexts surrounding the A3A M1 and M13 initiation codons.</caption><thead><tr><th>Monkey</th><th>Met1</th><th>Kozak context</th><th>Met13</th><th>Kozak context</th></tr></thead><tbody><tr><td><i>H. sapiens</i></td><td>CACATGG</td><td>A (dequate)</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>P. troglodytes</i></td><td>CACACGG</td><td>N null)</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>P. paniscus</i></td><td>CACACGG</td><td>N</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>G. gorilla</i></td><td>CACGTGG</td><td>N</td><td>TTCATGG</td><td>A</td></tr><tr><td><i>P. pygmaeus</i></td><td>CACATGG</td><td>A</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>Hylobates lar</i></td><td>CACATGG</td><td>A</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>C. neglectus</i></td><td>CACATGG</td><td>A</td><td>GTGATGG</td><td>S (trong)</td></tr><tr><td><i>C. aethiops</i></td><td>CACATGG</td><td>A</td><td>GTGATGG</td><td>S</td></tr><tr><td><i>C. guereza</i></td><td>CACATGG</td><td>A</td><td>GTGATGG</td><td>S</td></tr><tr><td><i>M. sphinx</i></td><td>CACATGG</td><td>A</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>M. mulatta</i></td><td>CACATGG</td><td>A</td><td>TTGATGG</td><td>A</td></tr><tr><td><i>S. oedipus</i></td><td>GACATGG</td><td>S</td><td>CTGATGG</td><td>A</td></tr><tr><td><i>S. sciureus</i></td><td>GACATGG</td><td>S</td><td>CTGATGG</td><td>A</td></tr></tbody></table>

doi:10.1371/journal.pone.0030036.t001

package and FEL and REL methods [41]. We investigated models in which the dN/dS ratio is allowed to vary among the complete sequence using the GA-branch analysis. There was significant positive selection with estimated dN/dS ratios >1.0 (p>0.95), at five sites, notably D41, L62, C64, H160 & H168 (Figure 1, in blue). By contrast, several sites were under significant negative selection, notably S7, N24, V25c, A107, F125, E157 and W162 (Figure 1, in green).

A phylogenetic tree for the complete sequence of A3A was constructed using BioNJ (Figure 2A). The red internal branches denote those where dN/dS>1 (p>0.9) which are confined to a small fraction of the total number of branches. Among A3 enzymes, the A3A sequence is phylogenically closest to the carboxy-terminal domains of A3B and A3G. In view of a large collection of A3G sequences [42], a comparable analysis was made using the A3Gc sequences (Figure 2B). The branch-specific patterns of dN/dS variation for both A3A and A3Gc cytidine deaminases are different, a good example being the New World monkey lineage.

## An A3A exon 3-related sequence on human chromosome 4

When performing Blat searches for this study (UCSC Genome Bioinformatics), we identified a segment of 288 bp on human

A

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Phylogeny of primate A3A and A3Gc cytidine deaminases. A) Neighbor-joining tree based on A3A primate sequences presented in Figure 1. Only bootstrap values >70 are indicated. The nomenclature for the Kozak initiation sites (A: adequate, N: null, S: strong) are from Table 1. B) Neighbor-joining tree based on A3Gc primate sequences taken in the literature. For both, red branches correspond to class dN/dS>1 (p>0.9) while the remainder correspond to classes dN/dS<1. C) The accepted divergence of the Great Apes.
doi:10.1371/journal.pone.0030036.g002</figcaption></figure>

B

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>
chromosome 4 with strong homology to exon 3 of the A3A/A3Bc/A3Gc cluster (Figure 3A) which will be referred to as ΨA3chr4. Homology went out to a few hundred bases either side with the splice sites perfectly conserved. The sequence is present in human, chimpanzee, gorilla, orang-utan, macaque and marmoset genomes while absent in horse, dog, cat and rodent genomes. At the protein level, the exon revealed a HVEX<sub>n</sub>SCX<sub>2</sub>C motif similar to that for all A3 deaminases (HAEX<sub>n</sub>PCX<sub>2-4</sub>C) (Figure 3A). While the A→V substitution is found in AID and APOBEC1 sequences, the P→S substitution is without precedent. Phylogenetic analysis based on amino acid sequences showed that it emerged after the (A3A, A3Bc)A3Gc split (Figure 3B).

5' and 3' RACE failed to identify any transcripts while no EST was found in the databases. Nonetheless, to ascertain whether this exon encoded a functional domain, we synthesized a fusion gene with the exon surrounded by exons 1, 2, 4 and 5 of the human A3A gene. The construct was cloned in pcDNA3.1 TOPO resulting in addition of the V5 tag. When transfected into HeLa cells and stained with FITC-conjugated anti-V5 antibody the construct was viable and strongly nuclear, more so than hA3A indicating that residues impacting A3A localization lie in exon 3 (Figure 3C). In

order to demonstrate editing activity, HeLa cells were co-transfected by the reconstructed pΨA3chr4 clone and an infectious molecular clone of hepatitis B virus. Total DNA was analysed at 72 hrs by a nested PCR/3DPCR approach as previously described [22,43].

The minimal denaturation temperature (Td) for the HBV X gene segment analysed is 91.8°C (Figure 3D, [22]). When co-transfected with the reconstructed pΨA3chr4 clone, the lowest Td was equally 91.8°C indicating that the recombinant may not be packaged into assembling HBV virions. Accordingly, a non-viral region corresponding to *MT-COI* gene was been amplified by PCR/3DPCR [33]. The minimal Td for *MT-COI* DNA was 87°C with or without pΨA3chr4 (Figure 3E), suggesting that the chromosome 4 fragment is indeed devoid of ssDNA cytidine deaminase activity.

Finally, an additional ~1.1 kb A3A-related fragment was identified ~8 kb upstream of the human A3A gene in the same orientation as the entire A3 locus. For comparison, the A3A-A3B intergenic region is ~19 kb. The fragment comprises 104 bp (37%) of intron 4, exon 5 and downstream sequences. Overall it shows 96% nucleic acid homology to hA3A. As such it must

<figure><img src="image_8.png" /><figcaption>**Figure 3.** No cytidine deaminase activity associated with pΨA3chr4. A) Alignment of exon 3 of ΨA3chr4 proteins to human A3Bc, A3Gc and A3A. Differences with respect to ΨA3chr4 are highlighted in red. B) Neighbor-joining tree based on the exon 3 of human A3 genes. Only bootstrap values of >70 are given. C) Immunostaining of the V5-tagged pΨA3chr4 and A3A proteins with DAPI nuclear counterstain. D) Agarose gel of 3DPCR products across a denaturation temperature gradient from 93 to 85°C of X region of HBV. pCayw is an infectious molecular clone along with the empty expression vector, along with pΨA3chr4. M; molecular weight markers. The white asterisks denote the last amplification product obtained at 91.8°C. E) Agarose gel of 3DPCR products across a denaturation temperature gradient from 89 to 79°C on *MT-COI* gene. The white asterisks denote the last amplification product of *MT-COI* obtained at 87°C.
doi:10.1371/journal.pone.0030036.g003</figcaption></figure>

<figure><img src="image_9.png" /></figure>
represent a vestige of prior gene conversion. Indeed, the sequence is surrounded by repeat elements, some of which are found surrounding the *hA3A* and *hA3B* genes. This *A3A* remnant is found in the chimpanzee, orang-utan and rhesus macaque genomes.

## Discussion

The primate *A3A* gene shows considerable qualitative and quantitative genetic variation, with up to 31% amino acid variation. Translation initiation sites vary there being at least four different configurations (Table 1). Positive selection is apparent along a few but not all branches suggesting that differences may emerge when more attention is turned to this important enzyme.

From the outset, differences in the restriction patterns of primate A3G on HIV-1Δ*vir* were noted [44,45,46]. More recent reports show that several human and macaque A3 cytidine deaminases are not strictly equivalent when using HIV-1 as a readout [47]. Indeed, as several reports have shown subtle differences for A3B, A3DE and A3G [47,48,49], this should transpire for A3A. However, as this enzyme impacts the integrity of the human genome, it is possible that the variation in structure and evolution of the *A3A* gene could impact cell biology.

During data analyses, two *A3A* related fragments were identified. The ΨA3chr4 exon 3 fragment proved to be devoid of catalytic activity when spliced together with exons 1, 2, 4 and 5 from *A3A*. This solo A3 exon is reminiscent of the recent finding of an isolated APOBEC1 exon in the tetrapod lineage that was subsequently lost [50]. The second *A3A* fragment is particularly interesting in that it shows that the present organization of the primate A3 locus might well have come about via more gene conversion than previously thought [51]. In conclusion, there is subtle qualitative and quantitative variation among primate *A3A* genes. In turn, gene expression and perhaps interferon sensitivity might follow.

## Materials and Methods

### Animal samples

Faecal samples were collected from wild non-habituated western gorilla (*Gorilla gorilla gorilla*) and chimpanzees (*Pan troglodytes troglodytes*) in Cameroon with permission of the Cameroonian Ministries of Health, Research and Environment and Forestry and Wildlife, and from bonobo (*Pan paniscus*) in the Democratic Republic of Congo with the permission of the Ministries of Science and Technology and Forest Economy [52]. DNA was extracted as previously described [53]. For mantled guereza (*Colobus guereza*) and mandrills (*Mandrillus sphinx*), DNA was extracted from whole blood on samples that were collected on primate bushmeat with permission from Cameroonian Ministries of Health, Research and Environment and Forestry and Wildlife, as previously reported [54]. Primary cells and cells line were obtained for orang-utan (*Pongo pygmaeus*) that died of natural causes while housed at the Wanariset orang-utan Reintroduction Center in East Kalimantan, Indonesia [55] and white-handed gibbon (*Hylobates lar*, ATCC 57763) respectively, while samples of rhesus monkey (*Macaca mulatta*, ATCC CCL-7), vervet monkey (*Cercopithecus aethiops*, ATCC CCL-81), and necropsy tissue samples from a squirrel monkey (*Saimiri sciureus*) and cotton-top tamarin (*Saguinus Oedipus*) that died of natural causes while kept in a zoo have been already described [56]. Primary cells from De Brazza's monkey (*Cercopithecus neglectus*) came from an animal that died of natural causes while housed at the zoo de la Palmyre (France).

### PCR amplification, cloning and sequencing

Hot start PCR was performed with corresponding primers (Table 2). The first reaction involved standard amplification, the reaction parameters were 95°C for 5 min., followed by 35 cycles (95°C for 30 s., 50–55°C for 30 s. and 72°C for 1 min.) and finally for 10 min. at 72°C for the first round. Differential amplification occurred in the second round using the equivalent of 0.2 µL of the first round reaction as input. Conditions were identical to the first PCR. The buffer conditions for all amplification were 2.5 mM MgCl₂, 50 mM KCl, 10 mM Tris-HCl pH 8.3, 200 µM of each dNTP, 100 µM of each primer, and 2.5 units of BIOTaq polymerase (Bioline) in a final volume of 50 µL. PCR products were purified from agarose gels (Qiaex II kit, Qiagen, France) and ligated into the TOPO TA cloning vector (Invitrogen, France). After transformation of Top10 electrocompetent cells (Invitrogen), up to 15 clones were picked. Sequencing was outsourced to GATC biotech. All mutations were confirmed by inspection of the chromatogram. The pΨA3chr4 insert was synthetized by GeneCust and cloned into the pcDNA3.1 TOPO-V5 vector (Invitrogen).

### Cells and transfections

Briefly, 10⁵ QT6 cells (ATCC CRL 1708) were cotransfected with 1 µg of pΨA3chr4 plasmid DNA along with 1 µg pCayw, a plasmid encoding an infectious molecular genome of hepatitis B virus (HBV) using FuGENE 6 (Roche). Total DNA was extracted using the MasterPureTM complete DNA and RNA purification kit (Epicentre). QT6 cells were maintained in HAM's F40 medium, supplemented with 1% chicken serum, 10% FCS, 5% tryptose phosphate, 2 mM L-glutamine, 50 U/ml penicillin and 50 µg/ml streptomycin [15].

### Immunofluorescence

HeLa cells (ATCC CCL 2) were grown to a density of 5.10⁵ cells per dish [15] and transfected with 1 µg of pΨA3chr4 or

<table><caption>Table 2. Compendium of all PCR primers used in the study.</caption><thead><tr><td>Primer</td><td>Sequence</td></tr></thead><tbody><tr><td>E1out5</td><td>5' AGAACTTCCTTGTTCTGATGCTAATGTGGGTGG</td></tr><tr><td>E1out3</td><td>5'AAAGGGAGGCCCCGGGGTTTGAGGGTGA</td></tr><tr><td>E1int5</td><td>5' GGCCTGAGCTGGAGAAGGGGTGGGGCAG</td></tr><tr><td>E1int3</td><td>5' TCAGAAAAGCTCAAAGAGGAGGCTGAGG</td></tr><tr><td>E2out5</td><td>5' GGCCTGAGCACACTGAGCTGACCCTGGG</td></tr><tr><td>E2out3</td><td>5' GAATGTCCTGGAYTGGRAGGGCCCTG</td></tr><tr><td>E2int5</td><td>5' GAGACCCTGACAAGGCTTAGACA</td></tr><tr><td>E3out5</td><td>5' TGGTCCAGGCCCTCCCTCCCTGTTCAC</td></tr><tr><td>E3out3</td><td>5' CARGCARCAGMCCAAGGCCTTTCTCC</td></tr><tr><td>E3int5</td><td>5' ACCTCACACTCTGTTTCCTTTTCTA</td></tr><tr><td>E3int3</td><td>5' CCCACCCCACCCCACCTGAACCTTCC</td></tr><tr><td>E4out5</td><td>5' CGGGARYGKACTTATCTCYCCTGTC</td></tr><tr><td>E4out3</td><td>5' GGAGGWGCCTGGGGACAGKGACAATGA</td></tr><tr><td>E4int3</td><td>5' GGGAGAAGGAGGGAGGGATGCGGGAGG</td></tr><tr><td>E5out5</td><td>5' TGTGCCCTCTTTCCACTCTCTCACCTCC</td></tr><tr><td>E5out3</td><td>5' TCTGCTGCTCAACCCAGGTCTCT</td></tr><tr><td>E5int5</td><td>5' CCACTCTCTCACCTCCTGCTCCA</td></tr><tr><td>E5int3</td><td>5' TCTCTGCCTTCCTTAGAGACTGA</td></tr></tbody></table>

doi:10.1371/journal.pone.0030036.t002

<figure><img src="image_10.png" /></figure>
pA3A using FuGENE 6 (Roche). After 48 hours, the cells were washed twice with PBS, fixed for 45 minutes in a 50:50 methanol/ ethanol mix. As primary antibodies, a mouse monoclonal antibody specific for the V5 epitope tag (Invitrogen) was used at a 1:200 dilution for 1 hour at room temperature. Cells were washed twice with PBS, and FITC-conjugated anti-mouse antibody anti-mouse was used as second antibody (Sigma) at a dilution 1:200 for 30 minutes at room temperature. We used Vectashield, mounting medium for fluorescence with DAPI (Vector laboratories, Inc.). Immunofluorescence was observed by microscopy (Zeiss).

## Phylogenetic and computational analyses

Sequences were aligned using the MUSCLE program, and neighbor-joining trees were obtained using BioNJ as implemented in http://phylogeny.fr. The final output was edited using Treeview [57]. The relative numbers of non-synonymous (dN) and synonymous (dS) nucleotide substitutions per site were estimated using the random effects likelihood (REL) and the fixed effects likelihood (FEL) methods available via the Datamonkey web interface of the HyPhy package [58]. Estimates of dN/dS ratios were based on neighbor-joining trees obtained from phylogeny.fr.

We used the genetic algorithm (GA-Branch) method available in HyPhy [58] to detect lineage-specific variation in selection pressure. This assigns different classes of dN/dS ratios to each

lineage to determine the best-fit model of lineage-specific evolution, and it calculates the probability (≥90%) that along a specific lineage dN/dS>1 [41].

Accession numbers were deposited at GenBank: *Colobus guereza* (JN177339), *Cercopithecus aethiops* (JN177340), *Cercopithecus neglectus* (JN177341), *Mandrillus sphinx* (JN177342), *Macaca mulatta* (JN177343), *Hylobates lar* (JN177344), *Gorilla gorilla* (JN177345), *Pongo pygmaeus* (JN177346), *Pan paniscus* (JN177347), *Pan troglodytes troglodytes* (JN177348), *Saguinus oedipus* (JN177349) and *Saimiri sciureus* (JN177350).

## Acknowledgments

The Molecular Retrovirology Unit is “Equipe labelisée LIGUE 2010”. Faecal samples and blood from bushmeat were collected in Cameroon with the approval of Cameroonian Ministries of Health, Research, and Environment and Forestry and Wildlife. We thank the zoo de la Palmyre and Dr Pascal Pineau for providing primate samples.

## Author Contributions

Conceived and designed the experiments: MH SWH JPV. Performed the experiments: MH. Analyzed the data: CT SWH JPV. Contributed reagents/materials/analysis tools: MP. Wrote the paper: SWH JPV.

## References

1. Bishop KN, Holmes RK, Sheehy AM, Davidson NO, Cho SJ, et al. (2004) Cytidine deamination of retroviral DNA by diverse APOBEC proteins. Curr Biol 14: 1392-1396.

2. Chelico L, Pham P, Calabrese P, Goodman MF (2006) APOBEC3G DNA deaminase acts processively 3'→5' on single-stranded DNA. Nat Struct Mol Biol 13: 392-399.

3. Harris RS, Bishop KN, Sheehy AM, Craig HM, Petersen-Mahrt SK, et al. (2003) DNA deamination mediates innate immunity to retroviral infection. Cell 113: 803-809.

4. Jarmuz A, Chester A, Bayliss J, Gisbourne J, Dunham I, et al. (2002) An anthropoid-specific locus of orphan C to U RNA-editing enzymes on chromosome 22. Genomics 79: 285-296.

5. Lecossier D, Bouchonnet F, Clavel F, Hance AJ (2003) Hypermutation of HIV-1 DNA in the absence of the Vif protein. Science 300: 1112.

6. Mangeat B, Turelli P, Caron G, Friedli M, Perrin L, et al. (2003) Broad antiretroviral defence by human APOBEC3G through lethal editing of nascent reverse transcripts. Nature 424: 99-103.

7. Mariani R, Chen D, Schrofelbauer B, Navarro F, Konig R, et al. (2003) Species-specific exclusion of APOBEC3G from HIV-1 virions by Vif. Cell 114: 21-31.

8. Sheehy AM, Gaddis NC, Choi JD, Malim MH (2002) Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein. Nature 418: 646-650.

9. Suspène R, Sommer P, Henry M, Ferris S, Guétard D, et al. (2004) APOBEC3G is a single-stranded DNA cytidine deaminase and functions independently of HIV reverse transcriptase. Nucleic Acids Res 32: 2421-2429.

10. Bonvin M, Achermann F, Greeve I, Stroka D, Keogh A, et al. (2006) Interferon-inducible expression of APOBEC3 editing enzymes in human hepatocytes and inhibition of hepatitis B virus replication. Hepatology 43: 1364-1374.

11. Chen H, Lilley CE, Yu Q, Lee DV, Chou J, et al. (2006) APOBEC3A is a potent inhibitor of adeno-associated virus and retrotransposons. Curr Biol 16: 480-485.

12. Delebecque F, Suspène R, Calattini S, Casartelli N, Saib A, et al. (2006) Restriction of foamy viruses by APOBEC cytidine deaminases. J Virol 80: 605-614.

13. Derse D, Hill SA, Princler G, Lloyd P, Heidecker G (2007) Resistance of human T cell leukemia virus type 1 to APOBEC3G restriction is mediated by elements in nucleocapsid. Proc Natl Acad Sci USA 104: 2915-2920.

14. Fan J, Ma G, Nosaka K, Tanabe J, Satou Y, et al. (2010) APOBEC3G generates nonsense mutations in human T-cell leukemia virus type 1 proviral genomes in vivo. J Virol 84: 7278-7287.

15. Henry M, Guétard D, Suspène R, Rusniok C, Wain-Hobson S, et al. (2009) Genetic editing of HBV DNA by monodomain human APOBEC3 cytidine deaminases and the recombinant nature of APOBEC3G. PLoS One 4: e4277.

16. Lochelt M, Romen F, Bastone P, Muckenfuss H, Kirchner N, et al. (2005) The antiretroviral activity of APOBEC3 is inhibited by the foamy virus accessory Bet protein. Proc Natl Acad Sci USA 102: 7982-7987.

17. Mahieux R, Suspène R, Delebecque F, Henry M, Schwartz O, et al. (2005) Extensive editing of a small fraction of human T-cell leukemia virus type 1 genomes by four APOBEC3 cytidine deaminases. J Gen Virol 86: 2489-2494.

18. Russell RA, Wiegand HL, Moore MD, Schafer A, McClure MO, et al. (2005) Foamy virus Bet proteins function as novel inhibitors of the APOBEC3 family of innate antiretroviral defense factors. J Virol 79: 8724-8731.

19. Sheehy AM, Gaddis NC, Malim MH (2003) The antiretroviral enzyme APOBEC3G is degraded by the proteasome in response to HIV-1 Vif. Nat Med 9: 1404-1407.

20. Suspène R, Guétard D, Henry M, Sommer P, Wain-Hobson S, et al. (2005) Extensive editing of both hepatitis B virus DNA strands by APOBEC3 cytidine deaminases in vitro and in vivo. Proc Natl Acad Sci USA 102: 8321-8326.

21. Turelli P, Mangeat B, Jost S, Vianin S, Trono D (2004) Inhibition of hepatitis B virus replication by APOBEC3G. Science 303: 1829.

22. Vartanian JP, Henry M, Marchio A, Suspène R, Aynaud MM, et al. (2010) Massive APOBEC3 editing of hepatitis B viral DNA in cirrhosis. Plos Pathog 6: e1000928.

23. Bishop KN, Holmes RK, Sheehy AM, Malim MH (2004) APOBEC-mediated editing of viral RNA. Science 305: 645.

24. Liddament MT, Brown WL, Schumacher AJ, Harris RS (2004) APOBEC3F properties and hypermutation preferences indicate activity against HIV-1 in vivo. Curr Biol 14: 1385-1391.

25. Wiegand HL, Doehle BP, Bogerd HP, Cullen BR (2004) A second human antiretroviral factor, APOBEC3F, is suppressed by the HIV-1 and HIV-2 Vif proteins. EMBO J 23: 2451-2458.

26. Zheng YH, Irwin D, Kurosu T, Tokunaga K, Sata T, et al. (2004) Human APOBEC3F is another host factor that blocks human immunodeficiency virus type 1 replication. J Virol 78: 6073-6076.

27. Suspène R, Aynaud MM, Koch S, Pasdeloup D, Labetoulle M, et al. (2011) Genetic editing of herpes simplex virus 1 and Epstein-Barr herpesvirus genomes by human APOBEC3 cytidine deaminases in culture and in vivo. J Virol 85: 7594-7602.

28. Koning FA, Newman EN, Kim EY, Kunstman KJ, Wolinsky SM, et al. (2009) Defining APOBEC3 expression patterns in human tissues and hematopoietic cell subsets. J Virol 83: 9474-9485.

29. Refsland EW, Stenglein MD, Shindo K, Albin JS, Brown WL, et al. (2010) Quantitative profiling of the full APOBEC3 mRNA repertoire in lymphocytes and tissues: implications for HIV-1 restriction. Nucleic Acids Res 38: 4274-4284.

30. Stenglein MD, Burns MB, Li M, Lengyel J, Harris RS (2010) APOBEC3 proteins mediate the clearance of foreign DNA from human cells. Nat Struct Mol Biol 17: 222-229.

31. Wang FX, Huang J, Zhang H, Ma X (2008) APOBEC3G upregulation by alpha interferon restricts human immunodeficiency virus type 1 infection in human peripheral plasmacytoid dendritic cells. J Gen Virol 89: 722-730.

32. Landry S, Narvaiza I, Linfesty DC, Weitzman MD (2011) APOBEC3A can activate the DNA damage response and cause cell-cycle arrest. EMBO Rep 12: 444-450.

33. Suspène R, Aynaud M, Guétard D, Henry M, Eckhoff G, et al. (2011) Somatic hypermutation of human mitochondrial and nuclear DNA by APOBEC3 cytidine deaminases, a pathway for DNA catabolism. Proc Natl Acad Sci USA 108: 4858-4863.

<figure><img src="image_11.png" /></figure>
34. Bogerd HP, Wiegand HL, Hulme AE, Garcia-Perez JL, O'Shea KS, et al. (2006) Cellular inhibitors of long interspersed element 1 and Alu retrotransposition. Proc Natl Acad Sci USA 103: 8780-8785.

35. Bulliard Y, Narvaiza I, Bertero A, Peddi S, Rohrig UF, et al. (2011) Structure-function analyses point to a polynucleotide-accommodating groove essential for APOBEC3A restriction activities. J Virol 85: 1765-1776.

36. Muckenfuss H, Hamdorf M, Held U, Perkovic M, Lower J, et al. (2006) APOBEC3 proteins inhibit human LINE-1 retrotransposition. J Biol Chem 281: 22161-22172.

37. Beale RC, Petersen-Mahrt SK, Watt IN, Harris RS, Rada C, et al. (2004) Comparison of the differential context-dependence of DNA deamination by APOBEC enzymes: correlation with mutation spectra in vivo. J Mol Biol 337: 585-596.

38. Conticello SG, Thomas CJ, Petersen-Mahrt SK, Neuberger MS (2005) Evolution of the AID/APOBEC family of polynucleotide (deoxy)cytidine deaminases. Mol Biol Evol 22: 367-377.

39. Wang M, Rada C, Neuberger MS (2010) Altering the spectrum of immunoglobulin V gene somatic hypermutation by modifying the active site of AID. J Exp Med 207: 141-153.

40. Thielen BK, McNevin JP, McElrath MJ, Hunt BV, Klein KC, et al. (2010) Innate immune signaling induces high levels of TC-specific deaminase activity in primary monocyte-derived cells through expression of APOBEC3A isoforms. J Biol Chem 285: 27753-27766.

41. Pond SL, Frost SD (2005) A genetic algorithm approach to detecting lineage-specific variation in selection pressure. Mol Biol Evol 22: 478-485.

42. Sawyer SL, Emerman M, Malik HS (2004) Ancient adaptive evolution of the primate antiviral DNA-editing enzyme APOBEC3G. PLoS Biol 2: e275.

43. Suspène R, Henry M, Guillot S, Wain-Hobson S, Vartanian JP (2005) Recovery of APOBEC3-edited human immunodeficiency virus G→A hypermutants by differential DNA denaturation PCR. J Gen Virol 86: 125-129.

44. Bogerd HP, Doehle BP, Wiegand HL, Cullen BR (2004) A single amino acid difference in the host APOBEC3G protein controls the primate species specificity of HIV type 1 virion infectivity factor. Proc Natl Acad Sci USA 101: 3770-3774.

45. Schrofelbauer B, Chen D, Landau NR (2004) A single amino acid of APOBEC3G controls its species-specific interaction with virion infectivity factor (Vif). Proc Natl Acad Sci USA 101: 3927-3932.

46. Zennou V, Bieniasz PD (2006) Comparative analysis of the antiretroviral activity of APOBEC3G and APOBEC3F from primates. Virology 349: 31-40.

47. Hultquist JF, Lengyel JA, Refsland EW, Larue RS, Lackey L, et al. (2011) Human and Rhesus APOBEC3D, APOBEC3F, APOBEC3G, and APOBEC3H Demonstrate a Conserved Capacity to Restrict Vif-deficient HIV-1. J Virol 85: 11220-11234.

48. Dang Y, Abudu A, Son S, Harjes E, Spearman P, et al. (2011) Identification of a single amino acid required for APOBEC3 antiretroviral cytidine deaminase activity. J Virol 85: 5691-5695.

49. Wissing S, Montano M, Garcia-Perez JL, Moran JV, Greene WC (2011) Endogenous APOBEC3B restricts LINE-1 retrotransposition in transformed cells and human embryonic stem cells. J Biol Chem 286: 36427-36437.

50. Severi F, Chicca A, Conticello SG (2011) Analysis of reptilian APOBEC1 suggests that RNA editing may not be its ancestral function. Mol Biol Evol 28: 1125-1129.

51. Larue RS, Andresdottir V, Blanchard Y, Conticello SG, Derse D, et al. (2008) Guidelines for Naming Non-Primate APOBEC3 Genes and Proteins. J Virol 83: 494-497.

52. Liu W, Li Y, Learn GH, Rudicell RS, Robertson JD, et al. (2010) Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature 467: 420-425.

53. Keele BF, Van Heuverswyn F, Li Y, Bailes E, Takehisa J, et al. (2006) Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science 313: 523-526.

54. Aghokeng AF, Ayouba A, Mpoudi-Ngole E, Loul S, Liegeois F, et al. (2010) Extensive survey on the prevalence and genetic diversity of SIVs in primate bushmeat provides insights into risks for potential new cross-species transmissions. Infect Genet Evol 10: 386-396.

55. Warren KS, Heeney JL, Swan RA, Heriyanto, Verschoor EJ (1999) A new group of hepatnaviruses naturally infecting orangutans (Pongo pygmaeus). J Virol 73: 7860-7865.

56. Pineau P, Henry M, Suspene R, Marchio A, Dettai A, et al. (2005) A universal primer set for PCR amplification of nuclear histone H4 genes from all animal species. Mol Biol Evol 22: 582-588.

57. Page RDM (1996) Treeview: an application to display phylogenetic trees on personal computers. Computer Applications in the Biosciences 12: 357-358.

58. Pond SL, Frost SD, Muse SV (2005) HyPhy: hypothesis testing using phylogenies. Bioinformatics 21: 676-679.

<figure><img src="image_12.png" /></figure>