Abstract
Free full text
![Logo of plospath](/http/europepmc.org/corehtml/pmc/pmcgifs/logo-plospath.png)
SIVagm Infection in Wild African Green Monkeys from South Africa: Epidemiology, Natural History, and Evolutionary Considerations
Abstract
Pathogenesis studies of SIV infection have not been performed to date in wild monkeys due to difficulty in collecting and storing samples on site and the lack of analytical reagents covering the extensive SIV diversity. We performed a large scale study of molecular epidemiology and natural history of SIVagm infection in 225 free-ranging AGMs from multiple locations in South Africa. SIV prevalence (established by sequencing pol, env, and gag) varied dramatically between infant/juvenile (7%) and adult animals (68%) (p<0.0001), and between adult females (78%) and males (57%). Phylogenetic analyses revealed an extensive genetic diversity, including frequent recombination events. Some AGMs harbored epidemiologically linked viruses. Viruses infecting AGMs in the Free State, which are separated from those on the coastal side by the Drakensberg Mountains, formed a separate cluster in the phylogenetic trees; this observation supports a long standing presence of SIV in AGMs, at least from the time of their speciation to their Plio-Pleistocene migration. Specific primers/probes were synthesized based on the pol sequence data and viral loads (VLs) were quantified. VLs were of 104–106 RNA copies/ml, in the range of those observed in experimentally-infected monkeys, validating the experimental approaches in natural hosts. VLs were significantly higher (107–108 RNA copies/ml) in 10 AGMs diagnosed as acutely infected based on SIV seronegativity (Fiebig II), which suggests a very active transmission of SIVagm in the wild. Neither cytokine levels (as biomarkers of immune activation) nor sCD14 levels (a biomarker of microbial translocation) were different between SIV-infected and SIV-uninfected monkeys. This complex algorithm combining sequencing and phylogeny, VL quantification, serology, and testing of surrogate markers of microbial translocation and immune activation permits a systematic investigation of the epidemiology, viral diversity and natural history of SIV infection in wild African natural hosts.
Author Summary
We simultaneously assessed, for the first time in a natural host, the epidemiology, diversity and natural history of SIVagmVer infection in wild vervet populations from South Africa. We report that African green monkeys (AGMs) have likely been infected with SIVagm for a long period, ranging from the time of their speciation to Plio-Pleistocene migrations, refuting previous molecular clock calculations suggesting SIVagm to be of recent occurrence. As a result of virus-host coadaptation, SIVagmVer infection is characterized by a lack of disease progression in spite of robust viral replication. We show that very active SIVagm transmission in adult AGMs contrasts with a very limited transmission to their offspring, in spite of massive exposure to SIVagm both in utero and through breastfeeding. The observation that some AGMs remain uninfected in spite of life-long exposure to SIVagm identifies wild vervets as an acceptable animal model for the exposed uninfected individuals, which can be used to identify correlates of resistance to HIV/SIV infection.
Introduction
Over 40 African nonhuman primate (NHP) species are naturally infected with simian immunodeficiency viruses (SIVs) [1]–[3]. Among these, African green monkeys (AGMs) (Chlorocebus genus) are the most numerous, most widely geographically spread and most commonly infected with SIV in the wild [1].
According to Groves [4], [5], AGMs are divided into species that are phenotypically and geographically distinct: vervets (C. pygerythrus) are the most widely spread in East and Southern Africa, ranging from the eastern Rift Valley in Ethiopia, Somalia and extreme southern Sudan to the Cape region in South Africa; grivets (C. aethiops) inhabit the area east of the White Nile in Ethiopia, Somalia from Khartoum to Mongalla, Eritrea and Ethiopia south to the Omo; tantalus monkeys (C. tantalus) are the prevalent species in north Central Africa from Ghana to Sudan and Kenya; green monkeys (C. sabaeus) reside in West Africa from Senegal to the Volta River, and have been introduced to Cape Verde Island as well as to the Caribbean [4]. A fifth AGM species, the malbrouck (C. cynosuros), is located in Central and South-Central Africa ranging from the Albertine Rift in the Democratic Republic of Congo to the Atlantic coast and Northern Namibia and Zambia. Finally, a sixth species, the Bale Mountain vervet (C. djamdjamensis) has a very limited distribution in the bamboo forests of highland Ethiopia. This taxonomic classification, however, is not universally accepted [6].
With the exception of the malbrouck and the Bale mountain vervet, each of the AGM species have been shown to be infected in the wild with species-specific SIVagm (defined as SIVagmVer, SIVagmGri, SIVagmTan and SIVagmSab) [7]–[11], with the reported prevalence rates ranging from30 to 50% [8], [12]–[19]. Phylogenetic analyses of SIVagm strains have shown that the spectrum of genetic diversity of SIVagm surpasses that of other primate lentiviruses [7], [9], [19]–[21], presumed to be due to codivergence/host-dependent evolution (i.e., the evolution of SIVagm from grivet, vervet, tantalus and green monkey species occurred in concert with host speciation) [8], [10], [11]. Together, these findings indicate that AGMs have coexisted with SIVagm for a long period of time, possibly even before their radiation and divergence from their common ancestor, approximately 1.5–3 million years ago [22].
Pathogenesis studies support an ancient coevolutionary relationship by showing that SIVagm are highly adapted to their hosts and that SIVagm-infected AGMs generally do not progress to AIDS [23], [24]. Key features of SIVagm infection in AGMs, as established by experimental studies performed on captive AGMs, include: (i) active viral replication, with set-point viral loads (VLs) similar to or higher than those found in HIV-infected patients [25]–[29]; (ii) significant depletion of CD4+ T cells during acute infection [25], [26], [30], [31], followed by rapid restoration to near preinfection levels in the peripheral blood [17], [25], [26], [30], [31] and delayed and incomplete restoration at mucosal sites [30]; (iii) maintenance of the balance between Th17 and T regulatory cells, due to preservation of the Th17 cell subset [32]; (iv) vigorous but transient inflammatory responses to the virus during acute infection, which are resolved with the transition from acute-to-chronic infection [30], [33], [34]; (v) productive infection of short-lived cells [31]; (vi) partial control of virus replication by the adaptive immune responses to SIV [35]–[41];(vii) no significant increase in CD4+ T cell apoptosis during either acute or chronic infection [30], [42], [43]. These normal levels of CD4+ T cell apoptosis together with preservation of Th17 cells probably allow AGMs to avoid enteropathy, breaches in the mucosal barrier and subsequent microbial translocation (MT) [30], as well as chronic immune activation and disease progression, while allowing CD4+ T cell recovery in the presence of high VLs [30], [44]. Furthermore, similar to other natural hosts of SIVs [45], AGMs have a number of adaptations that spare CD4+ T cells from virus-mediated killing in vivo. These features include a lower fraction of CD4+ T cells expressing the CCR5 coreceptor [46], [47] and down-regulation of the CD4 molecule by T helper cells as they enter the memory cell pool [48]–[50]. Altogether, these adaptations support the concept that the benign course of SIV infection in natural hosts is the result of coevolution over millennia [1], [51].
The ancestry of SIV infections and the codivergence of SIV types with their host species, however, could be accounted for by a model of preferential host-switching [52]. In the case of AGMs it has recently suggested that the phylogenetic profiles of SIVagm can be explained by a pattern of west-to-east transmission of the virus across existing AGM geographic ranges and that the viral divergence may have occurred more recently than previously suggested [53]. This study included only a small number of SIV strains covering a limited geographical distribution, which may have biased a full appreciation of the phylogenetic relationships and ancestry of SIVagm from different AGM species and likely altered the timing calculations for SIVagm emergence. A recent study in which calibration of molecular clocks was based on biogeographical features, reported that SIVs are actually older than estimated based on molecular clock calculations only [2].
The major focus of studies carried out thus far in wild African NHP species has been on identification and characterization of species-specific SIVs [54], [55]. A common denominator of these studies has been the difficulty of obtaining samples from wild animals and the lack of research infrastructure for pathogenesis studies in the field. Several approaches have been used to circumvent these limitations: (i) estimation of the SIV prevalence in captive monkeys, which significantly underestimates the prevalence in the wild, as most monkeys are captured at a young age and maternal-to-infant transmission rates are low in African NHP species [56]; (ii) use of bush meat samples. These are easy to obtain and thus this strategy has proven instrumental in identifying numerous SIVs and generating meaningful prevalence data [2], [57]–[62]. However, bush meat samples tend to be from adult populations and do not permit viral isolation or assessment of critical parameters of infection, such as VLs, immunophenotyping or assessment of the levels of immune activation; and (iii) use of noninvasive sampling (i.e., collection of fecal or urine samples), which has also been instrumental in assessing the prevalence of SIV infection in natural hosts, most notably in identifying and characterizing the animal reservoir of pandemic HIV strains [63]–[67]. However, similar to the use of bush meat samples, noninvasive samples cannot be used to assess key parameters of SIV infection. Thus, although SIVs have been identified and characterized for over two decades, studies of the natural history of SIV infection in the wild have only been performed at a cursory level [1], [54], [68].
Here, we performed for the first time in an African NHP host, a comprehensive large-scale characterization of the epidemiology and natural history of SIV infection in vervet monkeys living freely in different parks, nature reserves and farms from three provinces of the Republic of South Africa, to validate the virological and immunological features of experimental studies in natural hosts.
Results
Prevalence of retroviral infections in vervets from South Africa
Blood samples were collected from two hundred and twenty-five wild-trapped vervets over a 4.5 month period. Details on monkey capture and sample collection are presented as Supplementary information. The samples included represent 11 populations from 9 parks, game reserves and farms located in 3 provinces of South Africa–the Free State, KwaZulu-Natal, and Eastern Cape–covering majority significant proportion of vervet habitat in this region (Figure 1). Both sexes were equally represented. Samples from infants (<1 year old) and juvenile AGMs were included (Figure 1).
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g001.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g001.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g001.jpg)
Vervet samples from areas from three provinces (Free State, KwaZulu Natal and Eastern Cape) were included, as indicated. Sex and age-group coverage is illustrated.
Diagnosis of SIVagmVer infection was performed by PCR, which allows for the evaluation of SIVagm diversity in addition to prevalence estimates. Also, in infants, PCR permits the diagnosis of infected offspring (as opposed to serological assays which are not specific, i.e., may also detect passively transmitted maternal antibodies). RNA was extracted from all plasma samples and subjected to RT-PCR analysis using consensus primers designed to amplify a 600-bp fragment in pol integrase and a 900-bp fragment in the gp120 gene encompassing the V3–V5 regions and the 5′ end of the gp41 gene. In addition, gag PCR (846 bp fragment in the p24 gene) was performed on selected samples from each location. This analysis identified 103 AGMs that were virion RNA (vRNA) positive using one, two or three primer sets (Table 1), giving an overall prevalence of SIVagmVer infection of 59% (73/123) in females and 29% (30/102) in males, in the range of those found in previous reports in AGMs [8], [15], [17]. Prevalence levels were relatively similar between different locations (Table S1) and higher than in a cohort of semifree animals (Table S2). SIVagm prevalence was very uneven in different age groups: 7% (3/44) in infants (males: 4%, 1/26; females: 11%, 2/18), 16% (9/58) in juveniles (males: 15%, 5/34; females: 21%, 4/21) and 71% (91/128) in adults (males: 57%, 24/42; females: 78%, 67/86). Thus, we confirmed on a very large number of samples previous results reporting that a dramatic increase in the SIVagm prevalence and incidence occurs with the passage to sexual maturity in AGMs [16]. Interestingly, and different from previous reports, we identified cases of SIVagmVer infection in very young AGMs, which may be suggestive of vertical transmission of the virus in the wild.
Table 1
Females | Males | |||||
Total | SIV+ | Prevalence (%) | Total | SIV+ | Prevalence (%) | |
Infants | 18 | 2 | 11 | 26 | 1 | 4 |
Juveniles | 19 | 4 | 21 | 34 | 5 | 15 |
Adults | 86 | 67 | 78 | 42 | 24 | 57 |
Genetic diversity of SIVagm in vervets from South Africa
To determine the relationships of the newly identified SIVagm strains to each other and to previously characterized SIVagm strains, we constructed phylogenetic trees from amplified pol and env nucleotide sequences using a maximum likelihood approach (Figure 2 and Figure S1). While newly identified SIVagmVer strains naturally infecting vervets from South Africa exhibited a high genetic diversity, with average genetic distances in the pol gene of 16.2±4.8%, phylogenetic analyses also identified SIVagm strains that differed in less than 2% of their pol and env nucleotide sequences, indicating epidemiologically linked infections (Figure 2 and Figure S1). In general, strains originating from vervets from the same area clustered together (Figure 2) with a few exceptions that are probably due to male migration between groups when they reach sexual maturity.
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g002.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g002.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g002.jpg)
Maximum likelihood tree for pol gene. Sequences are colored depending on the region in which the sequences were sampled, with red indicating sequences from Free State, green indicating sequences from KwaZulu-Natal, and blue indicating sequences from Eastern Cape territories. Reference SIVagm strains from vervet monkeys in Kenya are shown in black. The tree is based on a 592 bp pol sequence alignment (after gap-containing sites were removed). Maximum likelihood estimates were performed using 1000 replicate bootstrap analysis. Internal nodes indicate level of support values for internal branching. Bar represents number of amino acid replacements (0.09) per site. The strain nomenclature includes the identification number, monkey status (M-male, F-female; I-infant; J-juvenile and A-adult),and the site and state of origin (FS-Free State; KZN-KwaZulu Natal; EC-Eastern Cape). Putative heterosexual transmission cases documented to be of recent occurrence based on VL and serology testing are underlined. Putative maternal-to-infant transmission cases are boxed.
As expected, phylogenetic analyses showed that SIVagmVer strains from South Africa clustered with the SIVagmVer reference strains (Figure 2 and Figure S1) and were only distantly related to the SIVagm strains from other AGM species (Figure 2). Furthermore, the newly characterized SIVagmVer strains from South Africa formed a subcluster within the SIVagmVer group, in agreement with the fact that all the reference strains available thus far have a distinct geographical origin, being collected from vervets in Kenya, Ethiopia and Tanzania [9], [19], [20]. Average genetic distances between the reference SIVagmVer strains (collected from Ethiopia, Kenya and Tanzania) and the newly characterized SIVagmVer strains were 25.4±2.1%, significantly higher than the overall genetic distances between South African strains (p<0.0001) (Table S3).
Similar to previous reports [66], comparison between the pol and env tree topologies showed numerous differences (Figure 2 and Figure S1), indicating extensive recombination within the SIVagm radiation (Figure S2). Mosaic branches occurred both near the tips of the trees and deeper in the trees, indicating that extensive recombination events have occurred throughout the evolutionary history of SIVagmVer [66]. This high frequency of recombinations between divergent SIVs indicates that SIVagm coinfection and/or superinfection occur frequently in wild vervets.
Remarkably, a geographically-associated clustering pattern can be defined for the South African SIVagmVer strains, with all but one of the SIVagmVer sequences from AGMs inhabiting different areas in the Free State forming a monophyletic cluster in the phylogenetic trees (Figure 2 and Figure S1).
Bootstrap analysis strongly supports a single origin of the Free State cluster, however parental nodes were associated with lower bootstrap support in the pol tree, suggesting ambiguity as to where the root of the Free State cluster should be placed. Lower bootstrap support may be indicative of recombination; a SplitsTree analysis shows that some KwaZulu-Natal sequences are likely to be recombinants of East Coast sequences, indicating a mixing between these two populations (Figure S2).
Interestingly, the SplitsTree analysis also identified one Free State sequence (VMT17083-IM-Gariep_FS) as potentially stemming from a distant recombination event with outgroup SIVagmVer sequences, this may suggest that this sequence is in fact a recombinant with an unknown parent (Figure S2).
The most likely explanation for these phylogenetic relationships is that the Drakensberg Mountains acted as a barrier to separate the vervet populations. We estimate that this separation could have occurred in a time frame ranging from 3 million years ago (during the AGM spread throughout sub-Saharan Africa) to 100,000 years ago (during the mass migrations that occurred in the Plio-Pleistocene glacial periods). We tested this hypothesis by calibrating a relaxed molecular clock to specify the time to the most recent ancestor (TMRCA) of the SIVagm samples surrounding the Drakensberg Mountains range to a time-frame spanning 100,000 to 3,000,000 years since the present for the env tree, and in the pol tree we additionally included SIV strains from several species and further calibrated the tree based on isolates from Bioko Island [2]. According to this analysis, divergence of the Free State lineage is likely to be closer to the 100,000 year mark, with mean most recent common ancestor (MRCA) for the South African SIVagm strains estimated at 329 thousands of years (kYa) (with 95% highest posterior density of 100–1,077 kYa) and 104 kYa (100–112 kYa), in the env tree (Table S4) and pol tree (Figure 3 and Table S4), respectively. This analysis places the root of the SIV tree older than previous studies [2], with the env tree giving a time estimate of SIV endemicity in simian species for 790 kYa (range: 151–2,595 kYa), and the pol tree indicating 248 kYa (range: 189–317 kYa) (Table S4).
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g003.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g003.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g003.jpg)
Relaxed molecular clock tree show the ultrametric divergence pattern between lineages. Calibration was based on internal nodes indicated with a red circle, by setting priors that span the presumed divergence MRCA date ranges, 100,000–3,000,000 years for SIVagm sequences surrounding the Drakensberg Mountains in South Africa, and 10,000±1,000 for SIVdrl sequences sampled from Bioko Island [2]. Internal nodes show divergence date estimates. Sequences are colored based on the territory from which they were sampled, with South African sequences colored red, green and blue, for Free State, KwaZulu Natal and Eastern Cape, respectively, and Bioko sequences colored in yellow.
Analyses of the env gene
SIVagmVer envelope sequences (909 bp) encompassing the V3–V5 regions of the surface (SU) envelope glycoprotein and the 5′ end of the transmembrane (TM) domain were amplified by nested PCR from 48 plasma samples. Alignment of their deduced amino acid sequences with those of published reference SIVagm strains from all four species of AGMs identified hypervariable (V3–V5) as well as conserved envelope domains (Figure 4), as reported previously [10].
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g004.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g004.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g004.jpg)
PCR-derived env nucleotide sequences were translated, aligned with previously reported SIVAGM env sequences [133] and compared with a consensus sequence generated by BioEdit. Dots denote sequence identity with the consensus sequence, while dashes represent gaps introduced to optimize the alignment. Triangles (∧) denote N-linked glycosylation sites. Bullets (•) mark CD4 binding sites. V3, V4 and V5 designate hypervariable SIVagm Env domains as previously described [10]. The envelope precursor cleavage site is indicated by an arrow above the consensus sequence. Strain nomenclature includes the identification number, monkey status (M-male, F-female; I-infant; J-juvenile and A-adult), and the site and the state of origin (FS-Free State; KZN-KwaZulu Natal; EC-Eastern Cape).
Regions corresponding to HIV-1 envelope domains of known function, including the envelope glycoprotein precursor cleavage site, the CD4 binding domain and the viral fusion peptide (N-terminus of gp41) were all highly conserved (Figure 4). Also conserved were the majority of cysteine residues (11 of 12) and, as reported previously [10], [18], there was very little sequence variability in the region that corresponds to the hypervariable V3 loop of HIV-1. Pairwise amino acid sequence comparisons were performed to examine the extent of SIVagmVer genetic diversity in the PCR amplified envelope fragment. As for the pol sequence analysis, env analysis revealed that SIVagmVer sequences collected from monkeys in the Free State form a distinct cluster from the monkeys living on coastal areas from both KwaZulu Natal and Eastern Cape provinces (Figure S1). This separate cluster was supported by both specific sequence signatures in the Free State-originating strains and by smaller genetic distances between these strains (Figure 4 and Figure S1).
Thus, sequence analysis identified: (i) a high diversity of SIVagm in South Africa, with the viruses collected from AGMs in the Free State forming an independent cluster, clearly separated from the strains collected from AGMs on the Indian Ocean coast; (ii) high frequency of recombination, suggestive of extensive transmission of the virus; and (iii) occurrence of highly related strains, suggestive of identifiable transmission clusters.
Assessment of the natural history of SIVagmVer infection in wild AGMs
It is generally assumed that SIV infection is nonpathogenic in natural hosts [23], [51], [56]. However, this paradigm was recently challenged by studies demonstrating that SIVcpz has a substantial negative impact on the health, reproduction and lifespan of chimpanzees in the wild [68]–[70]. We therefore took advantage of this unique sample set to assess for the first time the natural history of SIV infection in a natural host monkey species in the wild.
Clinical data
During sample collection, all included AGMs had a thorough clinical assessment. None of the vervets presented with any of the clinical signs associated with SIV infection, including fever, weight loss, lymphadenopathy or opportunistic infections. Although the cross-sectional nature of this study precludes assessment of changes in monkey weight and thus a direct assessment of weight loss, we performed comparative analyses of the body mass index (BMI) between SIV-positive and SIV-negative monkeys and showed that SIV status does not have an impact on normal weight of wild vervets (Figure S3).
Genotypic assessment of the coreceptor usage of SIVagmVer strains from South Africa
Coreceptor usage by the different SIV strains may have a significant impact on the natural history of SIV infection in the wild. Previous studies reported that SIVagmSab strains have the ability to use the CXCR4 coreceptor for viral entry, in addition to CCR5 [26]. Since the amounts of plasma available from the AGMs in South Africa precluded in vitro phenotypic testing of coreceptor usage by these new strains, we estimated coreceptor use in silico by analyzing the V3 region. V3 sequence impacts SIVmac and SIVagm tropism [71]–[73].The algorithms for genotypic assessment of coreceptor usage based on the HIV-1 V3 sequences cannot be implemented for SIV [74], [75]. We therefore assessed the overall net charges of the V3 loop for these sequences and found that all the viruses exhibited the same V3 loop net charge (+7), similar to that of other SIV strains that are documented to use CCR5 [1]. No significant differences in the potential N-glycosylation sites located in the V3–V5 region could be documented for the newly derived SIVagmVer strains from South Africa (Figure 4). We therefore inferred that SIVagmVer strains identified here likely use CCR5 as the major coreceptor for entry, similar to the majority of SIV isolates [26], [40], [76].
Staging of SIVagmVer infection in the wild: VL and serological testing
The levels of chronic viral replication are highly predictive of the outcome of HIV-1 infection [77], [78] and as we previously reported, may also predict the outcome of SIV infection in natural hosts, as the handful of African NHPs reported to progress to AIDS exhibited higher chronic VLs than nonprogressors [79]. To date, VLs in natural hosts of SIVs such as AGMs, sooty mangabeys and mandrills have only been assessed in either experimentally intravenously-infected monkeys [25], [80]–[84] or captive naturally infected monkeys [17], [28], [29], [79], [85], [86]. Only one study reported VLs in wild African monkeys on a very limited number of samples [85].
We used the generated SIVagmVer pol sequences to design specific primers and probes for viral quantification. The integrase region amplified by the pol primers used in our study is relatively well conserved; thus, the primers and probe had an excellent coverage of SIVagmVer diversity. We then quantified the VLs in 87 available samples from the 103 SIVagmVer-infected vervets. Results are shown in Figure 5. Plasma SIVagmVer VLs ranged from 104 to 107 copies/ml and were higher in the juvenile than in adult AGMs, albeit this difference was not statistically significant.
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g005.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g005.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g005.jpg)
No significant difference could be observed between infant, juvenile and adult AGM or between female and male monkeys. Samples with high VLs assessed as acutely infected are shown. VLs of lactating dams are depicted in lighter colors. Detection limit: 100 copies/ml.
At least 23 of the 49 SIVagm-infected dams for which the VLs were available (47%) were documented as lactating at the time of sampling (Figure 5). The VLs of lactating dams are detailed in Figure 5 and demonstrate a large exposure of infant AGMs to SIV through breastfeeding.
While VLs of the majority of SIVagm-infected animals fell within a relatively close range, consistent with the range previously reported in experimentally infected African NHPs that are natural hosts of SIV [27], VLs were higher than expected in approximately 10% of vervets (Figure 5). Since such high VLs were observed in the majority of juvenile AGMs, we reasoned that they may be recent infections [27].
High levels of viral replication observed around the peak of VLs during acute HIV/SIV infection generally occur prior to seroconversion [87] and correspond to the Fiebeg II stage of HIV-1 infection [87]. Hence, to stage the SIVagm infection in AGMs with high VLs, we further designed a gp41 peptide ELISA and assessed the levels of anti-SIVagm antibodies in the available samples. As shown in Figure 6, the majority of SIVagm-infected AGMs harbored detectable levels of anti-gp41 antibodies. Conversely, most of the AGMs presenting with high VLs were seronegative. Altogether, these results strongly suggest that AGMs with high VLs were in the acute stage of infection [87]. Note, however, that definitive proof in support of this conclusion would rely on serial sampling and demonstration of seroconversion and partial control of viral replication [87].
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g006.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g006.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g006.jpg)
Testing was done by peptide ELISA (mapping the immunodominant region of gp41 of SIVagmVer). Samples were grouped based on PCR and VL quantification results as originating from uninfected vervets, SIV-infected vervets with VLs>106 SIVagmVer RNA copies/ml, SIV-infected vervets with VLs ranging from 104 to 106 SIVagmVer RNA copies/ml and samples from infant vervets. Samples in the latter group are negative by PCR but show different levels of anti-SIVagm IgG, therefore being suggestive of transplacental crossing of maternal antibodies. Values on the Y-axis represent fold increases of the optical density of the samples (measured at 492 nm) over the cut-off value (established at an optical density of 0.2).
Interestingly, serological testing showed that the majority (72%) of the plasma samples collected from infant AGMs were seropositive for SIVagm (Figure 6), while SIVagm could be amplified from only 3 of the infant samples (Figure 2 and Table 1), suggesting passive transfer of maternal antibodies and documenting massive in utero exposure of AGM offspring to SIVagm.
In conclusion, monitoring the viral replication and the anti-gp41 antibodies in SIVagmVer-infected AGMs allowed us to demonstrate that: (i) naturally occurring chronic SIV infection in a natural host is characterized by high VLs, which are in the same range as in experimentally-infected animals [27]; (ii) SIVagm transmission is very active in the wild; and (iii) the relatively low prevalence of SIVagm infection in offspring strongly contrasts with their massive in utero and breastfeeding exposure.
Assessment of the integrity of mucosal barrier
Both SIV-infected AGMs and SMs were reported to maintain the integrity of mucosal barrier and thus to control MT [30], [84], which has been reported to be the main factor behind increased immune activation that drives disease progression in pathogenic HIV/SIV infections [88], [89]. To assess MT in wild vervets, we tested the levels of sCD14 as a surrogate biomarker. As shown in Figure 7, the levels of sCD14 were similar between SIV-infected and SIV-uninfected AGMs and were in the range of those we previously documented in captive uninfected and experimentally-infected AGMs. Even the samples collected from acutely infected monkeys did not show a significant increase in the levels of sCD14, in agreement with our previous results from captive experimentally infected AGMs [30], [90], and suggesting that wild SIVagm-infected AGMs maintain the integrity of the mucosal barrier.
Cytokine and chemokine testing
Finally, to assess the immune activation status between SIVagm-infected and uninfected AGMs in the wild, we measured and compared the levels of a large array of cytokines and chemokines in the available samples collected from the two groups in a Luminex assay. As shown in Figure 8, there was no significant increase in the cytokine/chemokine levels between SIV-infected and uninfected AGMs, with the exception of IL-6. In general, these cytokine/chemokine levels were in the same range of those previously reported by us and others [34] in captive uninfected and experimentally-infected AGMs during the chronic stage of SIVagm infection. This difference was mainly due to elevated levels of IL-6 in four of nine acutely-infected monkeys tested. When these samples were removed from the analysis, there were no significant differences in the levels of IL-6 between infected and uninfected AGMs. The increases in IL-6 probably reflect the transient increase in the levels of immune activation previously reported to occur during acute SIV infection in natural hosts [30], [33], [34], [91].
![An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g008.jpg An external file that holds a picture, illustration, etc.
Object name is ppat.1003011.g008.jpg](/http/europepmc.org/articles/PMC3547836/bin/ppat.1003011.g008.jpg)
P value was calculated by the Mann-Whitney U-test. Animals diagnosed as acutely infected are denoted by light red circles.
Thus, our analyses demonstrate the maintenance of the mucosal barrier and the lack of increased immune activation in acutely and chronically-infected AGMs in the wild. Immune activation being the major determinant of disease progression in pathogenic HIV/SIV infections, our results further support the clinical assessment of a benign outcome of SIV infection in wild AGMs.
Discussion
AGMs are the largest reservoir of SIV in the wild [23] and one of the main animal models for the study SIV pathogenesis in natural hosts [23], [92]. However, so far no study has systematically investigated the prevalence and incidence of SIVagmVer in vervets and whether the natural history of SIVagm infection in wild free living animals differs from that of captive populations. To address these two key questions, we performed a large-scale analysis of SIVagmVer in wild vervets from the Republic of South Africa. Using only 0.5 ml of plasma and a new complex algorithm combining sequencing and phylogeny with virus quantification, serology, and assessment of surrogate markers of microbial translocation and immune activation, we concomitantly investigated, for the first time in a wild natural host African NHP monkey species, SIV diversity and evolution simultaneously with epidemiology and natural history. The strategy outlined here enables a systematic investigation of the natural history of SIV infection in wild natural hosts.
Epidemiologic considerations
Vervets are abundant throughout East Africa along the Indian Ocean coast. However, most studies carried out thus far focused only on vervets from Northeast Africa (Kenya, Ethiopia and Tanzania) [14], and reported an SIVagmVer prevalence rate of 50% [17]. Conversely, our large-scale survey of SIVagm prevalence and diversity in vervets includes samples collected from all age groups and from multiple locations in South Africa. We confirm that the overall prevalence of SIVagm infection in the wild is high (46%), supporting the notion that SIVagm is uniformly distributed and widespread in AGMs [15], [16], [93]. These prevalence levels are similar to those previously reported for sooty mangabeys [66], [94], [95] and mandrills [96], which are among the most commonly infected NHP species known to date. For the first time, we document that one of the reasons for such extensive spread of SIV in wild vervets is the high steady-state SIVagm replication which probably facilitates virus transmission between AGMs through sexual contact [97].
SIVagm prevalence is higher in AGM females than in males. This difference is not surprising given that SIVagm transmission has lower efficacy through the foreskin compared to vaginal exposure [98] and that most males are expected to have significantly lower exposure, as only a few attain the dominance status necessary for sexual access to breeding females [99]. Note, however, that the robust levels of prevalence in males may be explained by the observation that despite the sexually restrictive nature of male dominance hierarchies in vervets, many subordinate males in multimale, multifemale groups can gain sexual access to females through a number of strategies [100]–[102], and thus become infected [103]. The relatively long breeding season, which lasts up to three months, may also reduce the ability of the alpha male to effectively monopolize fertile females, and increase the chances of transmission through sexual contact in subordinate males [104]. Furthermore, SIVagm transmission between males probably also occurs through aggressive, blood-exposing contacts for dominance, similar to other species [105], as many of the males sampled had open wounds or scars from such contact. Additionally, many females show marked aggression against dispersing males attempting to migrate into their social groups, often resulting in serious injury, providing another route of infection between males and infected females [106].
The very uneven distribution of SIVagm infection among different age groups strongly supports the sexual route as the main transmission route for SIVagm, together with a significantly less effective maternal-to-infant transmission. The low rates of mother-to-infant transmission are paradoxical, as our study documented large SIVagm exposure of the offspring both in utero and through breastfeeding. The rates of maternal-to-infant transmission in wild AGMs are significantly lower than the 35 to 40% mother-to-infant transmission rates reported in HIV infection (http://www.unaids.org) or 40–70% transmission rates in rhesus macaques [107] and are likely the result of virus-host coadaptation [97]. It is also possible that, similar to previous reports in captive sooty mangabeys and AGMs [108], [109], SIV infection of infant AGMs is associated with significantly lower levels of viremia than infection of adult monkeys, which may therefore go undetected by our assay.
Also conceivable is that, in addition to sexual and maternal-to-infant routes, a nonnegligible proportion of SIVagmVer transmissions occur through biting or fighting, as injuries are frequent (in our study sample, 10% of vervets showed signs of recent injuries that could result in exposure, such as deep lacerations). Injuries preponderantly occur during the mating season, due to both male transfer and to contests for dominance between males, and are usually cleaned by licking, which may involve many different individuals and may explain identification of near-identical viruses in vervets which are not direct mating partners. Studies have shown that monkeys are highly susceptible to SIV infection through oral exposure [110], [111], and the high levels of systemic viral replication documented in wild vervets may facilitate oral transmission during wound cleaning.
At least 10% of the cases of SIVagmVer infection in our group were recent infections (defined by high levels of viral replication close to the peak VL and negative serologies). Recent infections tend to occur in a narrow time frame (10–12 days postinoculation in experimentally infected macaques and African NHPs [27] and 3 weeks postinfection in HIV-1-infected patients [87]); therefore our results point to relatively high rates of SIVagm transmission in the wild. Note, however, that most transmissions are probably condensed in the time frame of the mating season, and thus the overall incidence of SIV infection in AGMs is not as high as suggested by our cross-sectional analysis.
Virus diversity and evolutionary considerations
Our study identified a significant degree of divergence between strains collected from relatively homogenous groups of vervets, similar to previous reports for other species of African NHPs hosts, such as mandrills and sooty mangabeys [61], [66], [95], [96]. Also, we have identified closely related strains, suggestive of identifiable chains of transmission. In some instances, such highly related strains were observed in mother-offspring pairs, and probably occurred as a result of direct mother-to-infant transmission of the virus. In other instances, clusters of closely related viruses included a strain collected from a monkey documented as being recently infected. No clear association could be established for several of the closely related strains, and, as previously mentioned, these infections might have occurred during other social interactions, such as grooming or wound licking.
SIVagmVer strains from different locations were generally separated on the phylogenetic trees. Intermixing of strains was also observed, suggestive of SIVagm spread between communities by dispersing males. Indeed, although in vervet populations males range with their social group, they also migrate between groups upon or shortly after reaching sexual maturity (sometimes multiple times during the males lifetime) [106], [112].
Interestingly, the SIVagmVer strains collected from Free State locations (west of the Drakensberg Mountains) clustered differently within pol and env trees (Figure 2 and S1) as well as in the gag tree (data not shown) from those collected from monkeys situated along the Indian Ocean coast (east of the Drakensberg Mountains). Given that the geographical spread between coastal monkeys and those from Free State is similar, the most likely explanation for these phylogenetic relationships is that the Drakensberg Mountains acted as a barrier (through low temperatures in the elevated areas) to separate the vervet populations on the two sides. The Drakensberg Mountains are very old, with an estimated age of 280 million years, which implies that a separation of the vervet populations could have occurred in a time frame ranging from 3 million years ago (during the AGM spread throughout sub-Saharan Africa) to 100,000 years go (during the mass migrations that occurred in the Plio-Pleistocene glacial periods). The Drakensberg Mountains may also have acted as an effective barrier to separate different populations of at least two other species, the common impala (Aepyceros melampus melampus) [113] and the chacma baboon (Papio ursinus) [114]. We therefore performed relaxed molecular clock analyses to specify TMRCA of the SIVagm samples surrounding the Drakensberg Mountains range to a time-frame spanning 100,000 to 3,000,000 years since the present for the env tree, and in the pol tree we additionally included SIV strains of several species and further calibrated the tree based on isolates from Bioko Island [2]. Molecular clock analysis using Bayesian techniques is highly dependent on priors that are specified [115]; however our analysis is in agreement with another study in which the use of biogeography showed that SIVs are very old [2], and refuted timing calculations based exclusively on molecular clocks without regards to biogeographical calibration points [53].
Taken together, these observations strongly suggest that SIVagm infection of vervets is more ancient than the most recent previous estimate [53] and suggests that the concept of host-dependent evolution (supported by the observation that the four SIVagm types cluster together in the SIV trees) cannot be discarded using current molecular clock calibrations, but may possibly be discernible through the integration of further biogeographical divergence estimates.
Natural history considerations
It is generally assumed that SIV infection is nonpathogenic in natural hosts [1], [23], [44], [51], [54]. Data supporting this assumption are derived from the study of captive monkeys and a very limited number of species (AGMs, sooty mangabeys and mandrills) [1], [23], [44], [51]. Sporadic cases of progression to AIDS were reported in each of these species upon infection with their species-specific viruses [116]–[118] or with cross-species transmitted SIVs [119]. Large scale studies in wild animals are needed to draw definitive conclusions about pathogenicity of SIVs in their natural hosts in the wild, as demonstrated by recent observation that SIVcpz infection is pathogenic in naturally-infected chimpanzees [68], [69], [120]. The fact that such a major effect went undetected for decades supports the need for systematic assessment of the natural history of SIV infection in African NHPs in their native habitat, and not in captive environments (where health status is controlled, nutrition is monitored and exposure to adventitious agents that may impact clinical status is limited), to provide significant information regarding the clinical outcome of SIV infection in natural hosts.
The main limitation of our assessment of natural history of SIVagmVer infection in vervets is that only plasma samples were available, and thus we could not perform immunophenotypic characterization or in vitro replication studies. We circumvented this limitation by employing a combination of biomarkers that could be monitored in plasma.
We report that SIVagmVer infection most likely has a benign outcome in vervets. BMI were not different between SIV-infected and SIV-uninfected vervets. Detected levels of viral replication were high, but in the range of those reported during experimental studies [27]. In the documented cases of acute SIVagmVer infection, the VL levels were also in the range of those observed in experimentally infected monkeys during acute infection [27].
We measured levels of sCD14 to assess the levels of MT, which is one of the proposed mechanisms of excessive immune activation that drives HIV/SIV infection to progress to AIDS [88], [89]. sCD14 is a surrogate marker for direct measurement of endotoxin or Gram negative bacteria, and numerous studies have shown a strong association of increasing sCD14 levels and increasing LPS levels in plasma in both pathogenic HIV and SIV infections (presumably due to loss of barrier function in the gut in turn due to HIV-mediated destruction of GALT) [88], [89], and lack of increase in both markers in the setting of nonpathogenic SIV infections in natural hosts [30], [84], [90], [121].
It is critical to assess the integrity of the mucosal barrier in a natural host in the wild, as we and others have previously reported that lack of MT during chronic infection in AGMs and sooty mangabeys may be one mechanism of controlling chronic immune activation and disease progression [30], [84]. These studies have been done using captive NHPs. Wild monkeys are exposed to many pathogens in a significantly more hostile environment than are captive animals, which may impact the ability of SIV-infected monkeys to maintain the barrier. However, we found no significant difference in the levels of sCD14 between infected and uninfected vervets in the wild.
In addition to being a surrogate marker of MT, sCD14 is a marker of macrophage activation, and thus our results also suggest that there are no significant differences in the levels of macrophage activation between SIV-infected and uninfected AGMs in the wild.
Lack of chronic immune activation was also supported by the similar levels of cytokines and chemokines tested in SIVagmVer-infected and uninfected wild vervets. The only exception was represented by an increased level of IL-6 in SIVagmVer-infected animals, which was driven by high levels of this cytokine in acutely infected monkeys. We and others reported that a transient increase in immune activation occurs during acute SIV infection in natural hosts [33], [35], [36], [91], [122], which is rapidly resolved with the transition to chronic infection [30]. Therefore, the similar levels of cytokine in plasma between SIV-infected and uninfected animals are somewhat surprising in the context of a significant number of acute infections in our study sample. This apparent discrepancy with the results of experimental studies may be due to the fact that most of the experimental studies involved intravenous administration of large amounts of virus that may trigger inflammation and immune activation. Alternatively, as we have previously shown, the acute increases in levels of plasma cytokines and chemokines are very transient in natural hosts, being thus likely that our cross-sectional samples missed such increases.
Altogether, sCD14 and cytokine/chemokine testing data support the conclusion that SIVagmVer infection of wild vervets is not associated with significant increases in levels of immune activation. Immune activation being the main driver of HIV disease progression, our results corroborate the clinical assessment and support a nonprogressive outcome of SIVagmVer infection in the wild.
Final considerations: Wild vervet populations from South Africa as an animal model for exposed-seronegative/exposed-uninfected (ESN/EU) patients
Both clinical experience and a growing medical literature indicate that some persons who have been exposed to HIV infection remain uninfected [123]. Although in some instances this may simply represent good fortune, cohorts of uninfected persons have been reported who are considered at high risk for infection. Clarifying the determinants of protection against HIV infection is a high priority that will require careful selection of high-risk uninfected cohorts, which should undergo targeted studies of plausible mediators and broad screening for unexpected determinants of protection. An animal model for ESN/EU may permit such an assessment while circumventing the major ethical boundaries to the study of ESN/EU humans. Such an animal model is not currently available. The SIVmac/RM model is not relevant for ESN/EU studies because this model was derived for pathogenesis and vaccine studies (thus employing strains selected for pathogenicity and not for their transmission characteristics) [124]. On the other hand, de novo development of such a model in RMs would be extremely expensive and face major macaque availability restrictions. In this context, the SIVagmVer-infected wild vervet cohorts that we identified in South Africa may serve as a useful model for ESN/EU studies. In these populations, there is an enormous exposure to SIVagm. However, while up to 80% of sexually-active adult females are infected with SIVagmVer, the remaining 20% are uninfected in spite of the presumably massive potential for exposure to SIV throughout their lives. Thus, field studies to assess the prevalence and incidence of SIV infection in the wild and prospectively confirm resistance to infection in a subset of ESN/EU monkeys, combined with genetic, immunologic and virologic assessment of the monkeys in the cohorts, may permit identification of correlates of ESN/EU, an instrumental leap for designing new therapeutic strategies aimed at preventing HIV infection.
In conclusion, we report that SIVagmVer is highly prevalent in wild vervet populations from South Africa, that the virus has long been present in this population and that coevolution between SIVagmVer and its natural host, the vervet monkey, resulted in an infection which has a benign outcome in the wild. More importantly, our study has identified a potential model for the study of ESN/EU cases, which has a strong potential to identify correlates of protection against HIV/SIV infection.
Materials and Methods
Ethics statement
All the animals sampled in this study were used according to regulations set forth by the Animal Welfare Act. The animal sampling protocols were approved by the University of Wisconsin-Milwaukee Institutional Animal Care and Use Committee (IACUC). Studies at the UCLA and University of Pittsburgh were carried out through an agreement with the UWM. At the University of the Free State UFS, ethical clearance was provided by the Interfaculty Animal Ethics Committee (project no. 13/2010).
Samples
The study was carried out in feral and semifree vervets (n=
225) living in different areas in South Africa, aged between 6 months and 10 years (Figure 1). Animals were individually trapped using established methods [125]. Details on animal capture are provided as Text S1. Blood was collected through venous puncture. Details are provided as Text S1. Each monkey included in the study had a microchip implanted for further identification and to prevent duplicate sampling. Detailed clinical assessment of the monkeys (including signs associated with immunodeficiency, i.e., lymphadenopathy, low weight, fever or rashes) was performed at the time of blood collection. Details are provided as Text S1. Upon collection, plasma was immediately stored at −80°C. A 500 µl aliquot of plasma from each monkey was available for this study.
Plasma processing and testing strategy
Upon arrival, plasma samples were thawed for RNA extraction and separated into three aliquots: the first aliquot was used for testing the prevalence levels of SIV by PCR. Additional SIV ELISA testing was performed using this aliquot to establish the stage of acute SIV infection and infant exposure to SIVagm. The second aliquot was used for quantification of SIVagm VLs by real-time PCR. Primers and probes were designed based on the SIVagm integrase sequences. Correlation between VLs and serology enabled identification of acutely SIV-infected AGMs [87]. Finally, the third aliquot was used to assess the natural history of SIVagm infection in the wild by testing surrogate markers for microbial translocation and immune activation (sCD14 and cytokine/chemokine levels).
Viral RNA extraction and cDNA synthesis
From each plasma specimen, viral RNA was extracted using the QIAamp Viral RNA Mini kit (QIAGEN, Germantown, MD). RNA was eluted and immediately subjected to cDNA synthesis. Reverse transcription of RNA to single-stranded cDNA was performed using SuperScript III reverse transcription according to the manufacturer's recommendations (Invitrogen, Carlsbad, CA). In brief, each cDNA reaction included 1×RT buffer, 0.5 mM of each deoxynucleoside triphosphate, 5 mM dithiothreitol, 2 U/ml RNaseOUT (RNase inhibitor), 10 U/ml of Super-Script III reverse transcriptase, and 0.25 µM antisense primers SIV-POL-OR (5′-ACB ACY GCN CCT TCH CCT TTC-3′), SIV-GAG-B (5′-CCT ACT CCC TGA CAG GCC GTC AGC ATT TCT TC-3′) and SIV-ENV-B (5′-AGA GCT GTG ACG CGG GCA TTG AGG TT-3′). The mixture was incubated at 50°C for 60 min, followed by an increase in temperature to 55°C for an additional 60 min. The reaction was then heat-inactivated at 70°C for 15 min and treated with 2 U of RNase H at 37°C for 20 min. The newly synthesized cDNA was used immediately or frozen at −80°C.
Nested PCR and DNA gel extraction
SIVagm RNA was amplified byPCR using Platinium Taq Polymerase High Fidelity (Invitrogen, Carlsbad, CA). A 600 bp pol integrase fragment was amplified by nested PCR using outer primers POL-IS4 (5′-CCA GCN CAC AAA GGN ATA GGA GG-3′) and POL-OR and inner primers SIV-POL-IS4 and SIV-POL-Unipol2 (5′-CCC CTA TTC CTC CCC TTC TTT TAA AA-3′) [59], [126]. An 846 bp gag fragment was amplified by nested PCR using outer primers GAG-A (5′-AGG TTA CGG CCC GGC GGA AAG AAA A-3′) and GAG-B, and inner primers GAG-C (5′-AGT ACA TGT TAA AAC ATG TAG TAT GGG C-3′) and GAG-F (5′-CCT TAA GCT TTT GTA GAA TCT ATC TAC ATA-3′) [127]. Finally, a 900-bp env fragment encompassing the V3-V5 gp120 region and the 5′ end of gp41 envelope regions was amplified by nested PCR using outer primers ENV-A (5′-GAA GCT TGT GAT AAA ACA TAT TGG AT-3′) and ENV-B, and inner primers ENV-C (5′-GTG CAT TGT ACA GGG TTA ATG AAT ACA ACA G-3′) and ENV-D (5′-TTC TTC TGC TGC AGTA TCC CAG CAA G-3′) [128]. Nested PCR products were subjected to 1% agarose (Invitrogen, Carlsbad, CA) gel electrophoresis and extracted and purified by QIAquick Gel Extraction kit (QIAGEN, Germantown, MD).
DNA sequencing
SIV gene amplicons were directly sequenced by the University of Pittsburgh Genomics and Proteomics Core laboratories (GPCL) using the nested PCR primers. Individual sequence for each amplicon was edited and inspected using BioEdit 7.1.3. Sequences were translated for analyses on the V3–V5 region of the SIVagmVer envelope with BioEdit.
Phylogenetic analysis
The gag, pol and env nucleotide sequence alignments were aligned with reference sequences obtained from the Los Alamos National Laboratory HIV Sequence Database (http://hiv-web.lanl.gov). The alignment for the pol region included nucleotide sequences sampled from the island of Bioko [2]. Newly derived SIV sequences were aligned using MUSCLE sequence alignment software [129] and alignments were edited manually where necessary.
Separate phylogenetic trees for each of the genes were inferred by maximum-likelihood using PhyML 3.0 [130], with a GTR substitution model, gamma-distributed rate heterogeneity at sites and an SPR tree search. Nonparametric bootstrap analysis with 1,000 replicates was performed to assess the reliability of branching order. In addition, calibrated molecular clock trees were calculated using BEAST v1.7.3 [115] to infer likely MRCA dates between different lineages. For molecular clock analysis, we used an HKY substitution model with gamma-distributed rate heterogeneity among sites and an uncorrelated lognormal relaxed molecular clock. A Yule tree prior was used for the pol region, and a constant size coalescent for the env region. Evolutionary rate was calibrated by setting a uniform prior bounded by 10,000 and 3,000,000 years, as TMRCA of SIVagm isolates sampled from areas surrounding the Drakensberg Mountains. For the pol region, a normal distribution with a mean of 10,000 years and a standard deviation of 1,000 years was additionally used to calibrate TMRCA of Bioko and mainland SIVdrl sequences. Chains were run for 107–108 states and ESS values were high (>300). Analysis of the molecular clock trees was performed using Tracer v1.5 (http://tree.bio.ed.ac.uk/software/tracer/), and trees were inferred by maximum clade credibility with mean node heights using TreeAnnotator v1.7.3, and visualized using FigTree v1.3.1 (http://tree.bio.ed.ac.uk/software/figtree/). To test for signals of recombination in the dataset, we ran the Neighbor-Net algorithm implemented in SplitsTree [131] to infer phylogenetic networks for the env and pol regions.
Viral load quantification
We used the generated pol alignment to design specific primers and probe for VL quantification in wild vervets: pol-F (5′-GTA GCC AGT GGG TTC ATA GAA GCA-3′), pol-R (5′-CAT TGC CTT TAC TTC CTG TGA AAT GA-3′) and pol-Probe (5′-/56-FAM/TAG AGA AAC/ZEN/AGG AAA AGA AAC AGC AA/3IABkF/-3′). Primers and probe were synthesized by IDTDNA (Coralville, IA) and were used in a one-step real-time PCR assay with TaqMan One Step PCR Master Mix (Applied Biosystems, Branchburg, NJ). Real-time PCR was performed in MicroAmp Optical 96-well plates (Applied Biosystems, Branchburg, NJ) by mixing the OneStep RT-PCR Enzyme Mix with 5 µl isolated RNA in a 50 µl final reaction volume. Real-time PCR conditions were as follows: 30 min at 50°C, 15 min at 95°C, followed by 45 cycles of 95°C for 15 s and 60°C for 1 min. Dilutions of all components were made using sterile RNase-free water. Data were collected and analyzed using the PE Applied Biosystems software provided. RNA copies/well were adjusted to copies per milliliter of original plasma. Samples were tested in duplicate and the number of RNA copies determined by comparison with a standard curve obtained using known amounts of SIV-pol RNA. The SIVagm RNA standard was produced by cloning one of the SIVagmVerpol gene fragments into Vector pBluescript SK+ and linearizing it by Xba I as a DNA template for RNA in vitro transcription. Standard RNA was produced by using MEGAscript kit (Applied Biosystems, Branchburgh, NJ), according to the manufacturer's instructions. Briefly, the transcription reaction system was assembled at room temperature and mixed thoroughly, followed by incubation at 37°C for 4 hours. One microliter of TURBO DNase was added to the reaction mixture followed by incubation at 37°C for 15 minutes to digest the DNA template. RNA was recovered by Lithium chloride precipitation. The amount of RNA was quantified at A260, aliquoted and immediately stored at −80°C.Detection limit of the SIVagmVer quantification assay was 100 copies/ml.
Serology
Based on the env sequence data, a peptide mapping the immunodominant region in the Gp41 transmembrane glycoprotein was synthesized and used in a peptide enzyme linked immunosorbent assay (ELISA), as described [132]. The inferred peptide sequence was TALEKYLEDQARLNVWGCAWKQVC, and this sequence was very well conserved between different SIVagmVer isolates, including previously reported reference strains from Kenya and Ethiopia [9], [10], [93]. The Gp36 peptide was synthesized to a purity of at least 90% (Fisher Scientific, USA) and the assay was performed as previously reported [132]. The cut-off of the reaction was arbitrarily set at 0.20, as per previous reports [132].
Microbial translocation (MT)
MT was assessed using the levels of soluble CD14 (sCD14) as a surrogate marker [89]. CD14 is a transmembranous protein, which also exists in a soluble form (sCD14; both as a shed membrane form and an alternatively spliced form), as a part of the complex that presents endotoxin (lipopolysaccharide or LPS) to TLR4 on monocytes. When monocytes are activated, ectodomain shedding results in increased sCD14 levels. sCD14 is therefore surrogate for direct measurement of endotoxin or Gram negative bacteria which translocate from the intestinal lumen to general circulation as a result of the immunologic injury inflicted at the mucosal level by the pathogenic HIV/SIV infection [88]–[90], [121]. sCD14 levels were measured using a quantitative sandwich enzyme immunoassay technique (Quantikine Human sCD14 Immunoassay, R&D Systems, Minneapolis, MN). The detection limit of this kit is 200 ng/mL and can range up to 5000 ng/mL, with an interassay coefficient of variability ranging between 7.19% and 10.9%.
Cytokine and chemokine testing
Cytokine testing in plasma was done using a sandwich immunoassay-based protein array system, Cytokine Monkey Magnetic 28-Plex Panel (Invitrogen, Camarillo, CA), as instructed by the manufacturer, and results were read by the Bio-Plex array reader (Bio-Rad Laboratories, Hercules, CA), which uses Luminex fluorescent-bead-based technology (Luminex Corporation, Austin, TX).
Statistical analyses
Statistical analyses were performed to analyze the results of viral load, sCD14 and cytokine testing using Mann-Whitney U-tests and Prism 4.0 software (Prism, Irvine, CA).
Nucleotide sequence accession numbers
The nucleotide sequences of the pol and env sequences from SIVagmVer-infected vervets from South Africa were deposited in the GeneBank (accession numbers: JX462308-JX462444).
Supporting Information
Figure S1
Maximum likelihood trees for env gene of the newly derived SIVagmVer sequences from wild vervets in South Africa. Maximum likelihood estimates were performed using 1000 replicate bootstrap analysis. Internal nodes indicate level of support values for internal branching. Sequences are colored depending on the region in which the sequences were sampled, with red indicating sequences from Free State, green indicating sequences from KwaZulu-Natal, and blue indicating sequences from Eastern Coast territories. In addition, sequences colored in gray are sampled from the Riverside Wildlife Rehabilitation and Environmental Education Centre (RWREC), Letsitele, Limpopo semifree colony that is housing released pets or animals with health problems, and their origin is unknown, and sequences colored in black indicate outgroup sequences. The strain nomenclature includes the identification number, monkey status (M-male, F-female; I-infant; J-juvenile and A-adult), and the site and the state of origin (FS-Free State; KZN-KwaZulu Natal; EC-Eastern Cape).
(TIFF)
Figure S2
SplitsTree showing the pol gene of South African SIVagmVer sequences. Nodes are colored by the region in which they were sampled, with red, blue and green denoting Free State, East Coast and KwaZulu-Natal, respectively. Potentially interesting recombinants have been circled, including VMT17083-IM-Gariep_FS and a cluster indicating possible mixing between East Coast and KwaZulu-Natal sequences.
(TIFF)
Figure S3
Comparative body mass index (BMI) assessment based on age and location in SIV-infected and SIV uninfected wild vervet monkeys ( Chlorocebus pygerithrus ) from South Africa.
(TIFF)
Table S1
SIV prevalence in vervet monkeys from South Africa based on location of collected samples.
(DOC)
Table S2
Age- and sex-related prevalence of SIVagmVer in semifree vervets from Riverside Wildlife Rehabilitation and Education Center, Letsitele, Limpopo.
(DOC)
Table S3
Estimates of evolutionary divergence over sequence pairs between groups. The number of base substitutions per site±standard error estimate(s) from averaging over all sequence pairs between groups are shown± Analyses were conducted using the Maximum Composite Likelihood model [1]. The rate variation among sites was modeled with a gamma distribution (shape parameter=
0.5). env analysis involved 81 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 822 positions in the final dataset for env. pol analysis involved 98 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 574 positions in the final pol dataset. Evolutionary analyses were conducted in MEGA5 [2].
(DOC)
Table S4
Most recent common ancestor (MRCA) estimates using relaxed molecular clocks to incorporate biogeographical assumptions. Calibration for the env gene is based on time to MRCA (TMRCA) of SIVagm sequences sampled around the Drakensberg Mountains, presumed to be 100,000–3,000,000 years of age. Calibration of the pol gene is based on TMRCA of SIVagm sequences as before, as well as SIVdrl sequences sampled from the Bioko Island of Equatorial Guinea, presumed to have been isolated for approximately 10,000 years. These estimates indicate that SIV is potentially considerably older than previous estimates, and refute timing estimates based exclusively on molecular clocks without calibration.
(DOC)
Text S1
Details on animal origin, animal capture, sampling methodology and clinical assessment.
(DOC)
Acknowledgments
We thank Drs. Beatrice Hahn, Vanessa Hirsch and Preston Marx for helpful discussion. Samples used in this study were collected as a part of the Systems Biology Sample Repository. We thank Department of Environmental Affairs, South Africa; Department of Economic Development and Environmental Affairs, Eastern Cape; Department of Tourism, Environmental and Economic Affairs, Free State Province; the Ezemvelo KZN Wildlife in KwaZulu Natal Province; and the Department of Economic Development, Environment and Tourism, Limpopo Province. Map of the Republic of South Africa used in Figure 1 was obtained from Wikipedia (http://en.wikipedia.org/wiki/File:South_Africa_sat.jpg).
Funding Statement
This work was supported by NIH grant R01OD010980 from the Office of Research Infrastructure Programs/OD (formerly R01RR016300 from NCRR) to NF and by NIH/NIAID/NCRR grants RO1 AI064066, R01 RR025781, P01 AI088564 to CA and IP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
Articles from PLOS Pathogens are provided here courtesy of PLOS
Full text links
Read article at publisher's site: https://doi.org/10.1371/journal.ppat.1003011
Read article for free, from open access legal sources, via Unpaywall:
https://journals.plos.org/plospathogens/article/file?id=10.1371/journal.ppat.1003011&type=printable
Citations & impact
Impact metrics
Citations of article over time
Alternative metrics
![Altmetric item for https://www.altmetric.com/details/1183239](https://api.altmetric.com/v1/donut/1183239_64.png)
Discover the attention surrounding your research
https://www.altmetric.com/details/1183239
Article citations
Adaptation of SIVmac to baboon primary cells results in complete absence of in vivo baboon infectivity.
Front Cell Infect Microbiol, 14:1408245, 28 Jun 2024
Cited by: 0 articles | PMID: 39006742
Making a Monkey out of Human Immunodeficiency Virus/Simian Immunodeficiency Virus Pathogenesis: Immune Cell Depletion Experiments as a Tool to Understand the Immune Correlates of Protection and Pathogenicity in HIV Infection.
Viruses, 16(6):972, 17 Jun 2024
Cited by: 0 articles | PMID: 38932264 | PMCID: PMC11209256
Review Free full text in Europe PMC
T cell activation is insufficient to drive SIV disease progression.
JCI Insight, 8(14):e161111, 24 Jul 2023
Cited by: 1 article | PMID: 37485874 | PMCID: PMC10443804
Walk on the wild side: SIV infection in African non-human primate hosts-from the field to the laboratory.
Front Immunol, 13:1060985, 12 Jan 2023
Cited by: 2 articles | PMID: 36713371 | PMCID: PMC9878298
Review Free full text in Europe PMC
Evidence of selection in the uncoupling protein 1 gene region suggests local adaptation to solar irradiance in savannah monkeys (Chlorocebus spp.).
Proc Biol Sci, 289(1982):20221254, 14 Sep 2022
Cited by: 1 article | PMID: 36100027 | PMCID: PMC9470266
Go to all (52) article citations
Data
Data behind the article
This data has been text mined from the article, or deposited into data resources.
BioStudies: supplemental material and supporting data
Nucleotide Sequences (3)
- (1 citation) ENA - VMT17083
- (1 citation) ENA - JX462444
- (1 citation) ENA - JX462308
Similar Articles
To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.
Factors associated with siman immunodeficiency virus transmission in a natural African nonhuman primate host in the wild.
J Virol, 88(10):5687-5705, 12 Mar 2014
Cited by: 38 articles | PMID: 24623416 | PMCID: PMC4019088
Simian immunodeficiency virus SIVagm.sab infection of Caribbean African green monkeys: a new model for the study of SIV pathogenesis in natural hosts.
J Virol, 80(10):4858-4867, 01 May 2006
Cited by: 99 articles | PMID: 16641277 | PMCID: PMC1472068
Mosaic genome structure of simian immunodeficiency virus from west African green monkeys.
EMBO J, 13(12):2935-2947, 01 Jun 1994
Cited by: 132 articles | PMID: 8026477 | PMCID: PMC395175
The history of SIVS and AIDS: epidemiology, phylogeny and biology of isolates from naturally SIV infected non-human primates (NHP) in Africa.
Front Biosci, 9:225-254, 01 Jan 2004
Cited by: 90 articles | PMID: 14766362
Review
SIVagm: genetic and biological features associated with replication.
Front Biosci, 8:d1170-85, 01 Sep 2003
Cited by: 27 articles | PMID: 12957815
Review
Funding
Funders who supported this work.
NCRR NIH HHS (3)
Grant ID: R01RR016300
Grant ID: R01 RR016300
Grant ID: R01 RR025781
NHLBI NIH HHS (1)
Grant ID: R01 HL117715
NIAID NIH HHS (2)
Grant ID: P01 AI088564
Grant ID: R01 AI064066
NIH HHS (2)
Grant ID: R01 OD010980
Grant ID: R01OD010980