Bonobos Exhibit Limited CD4 Diversity.

The CD4 receptor of chimpanzees contains amino acid substitutions at positions 25 (Q/R), 40 (Q/R), 52 (N/K), 55 (V/I), and 68 (P/T), all of which are the result of single nucleotide polymorphisms (SNPs) in exons 2 and 3 (55). These SNPs result in nine coding variants of the D1 domain, which appear to have evolved from an ancestral CD4 allele (QQNVP) by point mutation and recombination (55). Since bonobos (Pan paniscus) are the closest genetic relatives of chimpanzees (56), we asked whether they share CD4 polymorphisms. Bonobos live in the rain forests of the Democratic Republic of the Congo (DRC) south of the Congo River (SI Appendix, Fig. S1), and like chimpanzees, hunt and consume other primates (57); however, unlike chimpanzees, they do not appear to be SIV-infected (58). Extracting CD4 sequences from publicly available bonobo genome libraries (56), we obtained 14 full-length protein sequences, 13 of which were identical and differed from the ancestral chimpanzee QQNVP allele by a single amino acid replacement (Q429E) in the cytoplasmic domain (one sequence contained an additional R240W replacement in the D3 domain) (SI Appendix, Fig. S2). These results suggested that bonobos lack polymorphisms in the D1 domain.

To determine the full extent of the bonobo CD4 diversity, we used fecal samples from wild populations collected previously for molecular epidemiological studies of SIVcpz and ape malaria (58, 59). Specimens were available from 86 genotyped individuals sampled at 8 field sites (SI Appendix, Table S1), including from a genetically isolated population (TL) east of the Lomami River (SI Appendix, Fig. S1). Since the CD4 gene is comprised of multiple exons and introns spanning nearly 20 Kb on chromosome 12, we limited our analysis to exons 2 (165 bp) and 3 (159 bp), which together encode the entire D1 domain. Each exon was amplified separately because of a large intervening intron (>13 Kb) and sequenced without fragmentation to maintain linkage between polymorphic sites. In addition, each exon was amplified up to eight times to exclude allelic dropout (SI Appendix, Table S1). While exon 2 sequences were identical in all animals, exon 3 sequences differed by a single nonsynonymous substitution, which led to an amino acid replacement at position 83 (I83T) (Fig. 1A). This polymorphism was located at the periphery of the Env binding region at a site not previously reported to be under positive selection (Fig. 1B). Among the animals studied, the T83 variant was relatively uncommon (found in 14 bonobos at 4 field sites), while the I83 variant represented the predominant form (found in 72 bonobos at all field sites) (Fig. 1C and SI Appendix, Fig. S1). Thus, in contrast to chimpanzees, bonobos exhibit limited D1 domain diversity, with only two coding variants circulating in wild populations (Table 1), the most common of which represents the inferred ancestral CD4 allele of both chimpanzees and bonobos.

Fig. 1.

Allelic diversity of CD4 in bonobos. (A) CD4 coding variants identified in wild bonobo populations. D1 domain variants of bonobos (P. paniscus, Pp; derived from exon 2 and 3 sequences) are compared to human and chimpanzee CD4, with dots indicating identity to the human reference. A single amino acid replacement (I83T) in the bonobo CD4 is boxed. A potential N-linked glycosylation site is indicated by asterisks. The ancestral chimpanzee (P. troglodytes; Pt) allele is shown for reference, with positions that are polymorphic in chimpanzees underlined. (B) Crystal structure of the HIV-1 gp120 envelope domain (black) bound to human CD4 (gray) (PDB ID code 4R2G), with polymorphic positions indicated for chimpanzees (blue) and bonobos (green). (C) CD4 allele frequency in wild bonobos. The number of tested individuals is indicated. (D) Effects of bonobo CD4 polymorphisms on SIV Env mediated cell entry. The infectivity of pseudoviruses carrying different SIV Envs is shown for transiently transfected cells expressing human and bonobo CD4 variants and the cognate CC5 receptors. Values are scaled relative to the human CD4 (set to 100%). Bars represent the average of three independent transfections, each performed in triplicate, with SDs shown (fold-changes of Env infectivity for different CD4 alleles are shown in SI Appendix, Table S2).

Table 1.

D1 domain genotyping of African primates

Primate species	No. of individuals*	No. of D1 variants	SIV infection†	Non synonymous changes	Synonymous changes
Apes
Chimpanzee (Pan troglodytes)	544‡	9	+	6	1
Bonobo (Pan paniscus)	100§	2	−	1	0
Western gorilla (Gorilla gorilla)	97§	3	+	3	0
Eastern gorilla (Gorilla beringei)	25§	2	−	3	0
Cercopithecine monkeys
L’Hoest’s monkey (Allochrocebus lhoesti)	3	2	+	2	0
Sun-tailed monkey (Allochrocebus solatus)	2	3	+	2	1
Tantalus monkey (Chlorocebus tantalus	11§	3	+	3	0
Vervet monkey (Chlorocebus pygerythrus)	63§	4	+	4	0
Grivet (Chlorocebus aethiops)	26§	3	+	3	0
Green monkey (Chlorocebus sabaeus)	53§	4	+	3	0
Malbrouck (Chlorocebus cynosuros)	16§	4	+	3	0
Patas monkey (Erythrocebus patas)	5§	1	−	0	0
Red-tailed monkey (Cercopithecus ascanius)	4¶	2	+	2	1
Mustached monkey (Cercopithecus cephus)	15	11	+	8	0
Greater spot-nosed monkey (Cercopithecus nictitans)	17	9	+	8	1
Sykes’ monkey (Cercopithecus albogularis)	22	3	+	3	0
Blue monkey (Cercopithecus mitis)	7	3	+	3	0
De Brazza’s monkey (Cercopithecus neglectus)	3	1	+	0	0
Diana monkey (Cercopithecus diana)	5	3	?	2	0
Lesser spot-nosed monkey (Cercopithecus petaurista)	4	2	?	2	0
Crested mona monkey (Cercopithecus pogonias)	2	2	−	1	0
Red-capped mangabey (Cercocebus torquatus)	3	2	+	4#	0
Sooty mangabey (Cercocebus atys)	5‡	2	+	2#	0
Mandrill (Mandrillus sphinx)	7	2	+	1	1
Chacma baboon (Papio ursinus)	10	1	−	0	0
Olive baboon (Papio anubis)	13	1	−	0	0
Yellow baboon (Papio cynocephalus)	3	1	−	0	0
Colobine monkeys
Ugandan red colobus (Piliocolobus tephrosceles)	29§	6	+	7	1
Mantled guereza (Colobus guereza)	4	2	+	1	1

* In addition to the within-species diversity shown here, D1 domain sequences were also obtained for single individuals of the following species: Preuss’s monkey (Allochrocebus preussi), Hamlyn’s monkey Cercopithecus hamlyni), Lowe’s mona monkey (Cercopithecus lowei), Allen’s swamp monkey (Allenopithecus nigroviridis), Angolan talapoin (Miopithecus talapoin), Drill (Mandrillus leucophaeus), and Angola colobus (Colobus angolensis) (see SI Appendix, Table S5 for GenBank accession numbers).

^† Plus sign (+), naturally SIV infected; minus sign (−) not naturally SIV infected; question mark (?), insufficient sampling to determine SIV infection status.

^‡ D1 domain sequences were obtained from GenBank.

^§ All or some D1 domain sequences were extracted from whole genome or RNA-seq datasets (24, 48, 68). Note that the number of D1 domain variants represent minimum estimates, since alleles with low quality support were discarded.

^¶ Minimum number of individuals from 25 samples.

^# In addition to substitutions, a 3 bp deletion was observed in alleles of some red-capped and sooty mangabeys (Fig. 4A).

Although not naturally infected with SIV, bonobos are likely exposed to primate lentiviruses through their hunting behavior (57). Since the viruses infecting these prey species have not been characterized, we selected a panel of SIV Envs from diverse viral lineages and tested their ability to mediate entry into cells expressing the bonobo CD4 receptors (8, 11, 55). These included Envs from various SIVcpz and SIVgor strains, as well as SIVs infecting mustached monkeys, red-tailed monkeys (SIVasc; Cercopithecus ascanius), western red colobus (SIVwrc; Piliocolobus badius), tantalus monkeys (SIVtan; Chlorocebus tantalus), l’Hoest’s monkeys (SIVlho; Allochrocebus lhoesti), and sooty mangabeys. Briefly, 293T cells were transfected with plasmids expressing human and bonobo CD4 alleles as well as their cognate CCR5 coreceptors, and then infected with viruses pseudotyped with the various SIV Envs. Compared to human CD4, both bonobo alleles were less efficient in mediating cell entry (Fig. 1D and SI Appendix, Table S2), most likely because of a glycan at position 32 (N32), which sterically hinders SIV Env interaction in chimpanzees (55). However, all SIVcpz and SIVgor Envs, as well as most monkey SIV Envs, were able to utilize both bonobo CD4 variants. The I83 allele was as efficient in mediating cell entry as the chimpanzee QQNVP allele (SI Appendix, Fig. S3), which is not surprising given that the ectodomains of these two CD4 molecules are identical (SI Appendix, Fig. S2). The less frequent T83 allele appeared to enhance SIV Env-mediated cell entry (SI Appendix, Table S2), although this increase was only modest (on average 1.3-fold). Thus, bonobos lack the extensive CD4 diversity that appears to protect chimpanzees against primate lentiviruses. However, two Envs from SIVmus and SIVsmm utilized the bonobo CD4 variants only poorly (Fig. 1D). Although bonobos are not naturally exposed to these viruses, this finding suggests that bonobos may be resistant to some SIVs.

Western and Eastern Gorillas Encode Multiple D1 Domain Alleles.

The paucity of D1 domain polymorphisms in bonobos prompted us to examine western and eastern (Gorilla beringei) gorillas, which occupy nonoverlapping ranges in sub-Saharan Africa (SI Appendix, Fig. S1). Of these, only western lowland gorillas in Cameroon are SIVgor infected, having acquired this virus only once from sympatric chimpanzees (8, 11). Mining publicly available genome databases (56, 60), we extracted complete CD4 protein sequences from 14 western and eastern lowland gorillas. An alignment of these sequences revealed five polymorphic sites, four of which were located in the D1 domain at positions 18 (A/P), 27 (H/N/R), 31 (S/P), and 34 (R/M), while the fifth was located at position 142 (K/R) in the D2 domain (SI Appendix, Fig. S4). All D1 domain polymorphisms were caused by SNPs in exon 2, which yielded five variants, three of which were only found in western gorillas (AHSM, ARSM, ANSR), while the other two were only found in eastern gorillas (AHPM, PRSM) (Table 1).

To determine whether these alleles represent the entirety of D1 domain diversity, we used fecal DNA to sequence this region in 90 additional western and 18 eastern lowland (G. beringei graueri) gorillas sampled at 14 and 6 field sites throughout Cameroon and the DRC, respectively (SI Appendix, Fig. S1). Specimens were selected based on their geographic origin and individual information, with both SIVgor⁺ SIVgor⁻ western lowland gorillas included (SI Appendix, Table S3). Like for the bonobos, exons 2 and 3 were amplified separately because of the large intervening intron and sequenced without fragmentation. However, due to limited sample availability, some specimens could only be amplified once (SI Appendix, Table S3). Consistent with the genomic data, CD4 genotyping of wild gorillas confirmed the existence of AHSM, ARSM, and ANSR alleles in western and AHPM and PRSM alleles in eastern lowland gorillas, with no overlap between the two species (Fig. 2A). All D1 polymorphisms were located at the binding interface of Env and CD4 (Fig. 2B), with AHSM and AHPM representing the most frequent alleles in western and eastern gorillas, respectively (Fig. 2C). Although screening of wild gorillas did not uncover additional CD4 diversity, it is clear that additional intermediate alleles existed in these populations at some point in the past. For example, the ANSR allele of western gorillas is 2 amino acids different from its closest relatives, indicating three potential intermediates (AHSR, ANSM, ARSR), one of which must have existed in the past, if not now (SI Appendix, Fig. S5). Similarly, the AHPM and PRSM alleles in eastern gorillas differ by 3 amino acids, indicating six potential intermediates (AHSM, ARSM, ARPM, PHPM, PRPM, PHSM), two of which are present in western gorillas (SI Appendix, Fig. S5). Thus, there likely are additional CD4 alleles in both gorilla species that have not yet been identified.

Fig. 2.

Allelic diversity of CD4 in gorillas. (A) CD4 coding variants identified in western (G. gorilla gorilla; Ggg) and eastern (G. beringei graueri; Gbg) lowland gorillas. D1 domain variants in exon 2 are compared to the human CD4, with dots indicating sequence identity (exon 3 derived protein sequences were invariant). A potential N-linked glycosylation site is indicated by asterisks. Amino acid replacements are indicated in red, with allelic variants named based on the order of polymorphic amino acid residues. (B) Crystal structure of the HIV-1 gp120 envelope domain (black) bound to human CD4 (gray) (PDB ID code 4R2G) with polymorphic positions indicated for chimpanzees (blue) and gorillas (red). (C) CD4 allele frequencies in wild-living western (Upper) and eastern (Lower) gorillas. The number of tested individuals is indicated. (D) Effects of gorilla CD4 polymorphisms on SIV Env mediated cell entry. The infectivity of pseudoviruses carrying the SIV Envs indicated is shown for transiently transfected cells expressing human and gorilla CD4 variants and the cognate CCR5 receptors. Values are scaled relative to human CD4 (set to 100%). Bars represent the average of two or three independent transfections, each performed in triplicate, with SDs shown. MT145, MB897, EK505, LB715, and GAB2 represent SIVcpzPtt strains from central chimpanzees, while TAN2 and BF1167 represent SIVcpzPts strains from eastern chimpanzees. (E) Protective effect of the N15 glycan. The percentage of infected cells bearing the AHSM allele was compared to a mutant (N15T) lacking the N15 glycan (fold-changes of Env infectivity for different CD4 alleles are shown in SI Appendix, Table S4). Individual SIV Envs are color coded as in D.

To determine whether the D1 polymorphisms affected SIV Env-mediated cell entry, we synthesized one gorilla CD4 allele and generated the other D1 domain variants by site-directed mutagenesis. We then transfected 293T cells with plasmids expressing human and gorilla CD4, along with the corresponding CCR5 coreceptors, and infected these cells with pseudoviruses bearing the same SIV Envs used in the bonobo study. The results showed that all gorilla CD4 alleles were functional and facilitated pseudovirus entry, although infectivity levels were again lower compared to human CD4 (Fig. 2D). However, the five gorilla CD4 alleles differed in their ability to facilitate infection. For example, the ANSR allele mediated efficient entry of all SIVgor Envs, while the AHPM allele reduced infection of two of them (SI Appendix, Table S4). These allele-specific differences were even more pronounced when SIVcpz Envs were tested. For example, the MB897 Env utilized the AHSM and AHPM alleles, but was restricted by the ARSM, ANSR, and PRSM alleles, while the MT145 Env utilized the ARSM and PRSM alleles, but was restricted by the AHSM, ANSR, and AHPM alleles (Fig. 2D). Most, but not all, monkey SIV Envs utilized the five gorilla CD4 alleles with high efficiency (Fig. 2D and SI Appendix, Table S4). Thus, gorilla CD4 polymorphisms, like those in chimpanzees, inhibit SIV Env cell mediated entry in a strain and allele specific manner (55).

In addition to the five polymorphic sites, the gorilla CD4 encodes an invariant potential N-linked glycosylation site (PNGS) at position 15 (N15) that is absent in all other primate species (Fig. 2A). This PNGS, which has been experimentally confirmed to be glycosylated (61), appears to occupy a similar structural position as the chimpanzee glycan at position 66, which is known to interfere with SIV Env binding (55). To examine whether N15 has a protective effect, we changed the asparagine at position 15 in one representative gorilla CD4 allele (AHSM) to a threonine, which is present in all other African primates. Testing SIVgor, SIVcpz, and other SIV Envs, we found that removal of the N15 glycan increased Env mediated infectivity on average by 1.6-fold (SI Appendix, Table S4), with enhancement observed for nearly every Env (Fig. 2E). Thus, like the invariant N32 glycan in chimpanzees and bonobos, the invariant N15 glycan in gorillas provides some degree of protection against SIV infection.

Mustached Monkeys Exhibit an Extraordinary Degree of CD4 Diversity.

To examine whether CD4 diversity is unique to African apes, we next sequenced the D1 domain in mustached monkeys. We were particularly interested in this species because mustached monkeys harbor three different types of SIVmus, two of which represent recombinants with SIVgsn and other SIVs (14, 62, 63), but are infected at very low (0 to 6%) prevalence rates (2, 62, 6465–66). Using remnant DNA from previous blood collections (2, 64), we CD4-genotyped 15 members of this species (Table 1). Sequence analysis of exons 2 and 3 identified eight nonsynonymous SNPs, which resulted in amino acid replacements at positions 30, 34, 39, 40, 41, 42, 50, and 90, the combination of which indicated a minimum of 11 different CD4 coding variants (Fig. 3A). Of these, only nine D1 domain alleles could be unambiguously inferred, because amplicons from the remaining samples contained polymorphic sites in both exons 2 and 3. The combination of the latter resulted in at least two more alleles as well as a possible third (Fig. 3A), although their exact sequences could not be determined due to different possibilities of exon 2 and 3 linkage. One polymorphic site (I34T) affected a PNGS at position 32. Thus, analysis of a very small number of mustached monkeys identified an extraordinary degree of D1 domain diversity, with 11 of the 15 animals exhibiting a heterozygous CD4 genotype (SI Appendix, Table S5).

Fig. 3.

Allelic diversity of CD4 in mustached monkeys. (A) CD4 coding variants identified in mustached monkeys. Mustached monkey D1 domain variants derived from both exons 2 and 3 (indicated on the bottom) are compared to one representative allele, with dots indicating identity to this reference (sequences are trimmed to the polymorphic region). Polymorphic positions are highlighted in red and their position is indicated on the top. Alleles 1 to 9 could be unambiguously inferred; the remaining alleles remain ambiguous because of polymorphisms in both exons 2 and 3, which could not be linked. For individual 10, permutations of exons 2 and 3 combinations resulted in one new allele (either 10a or 10b), which was paired with an allele already identified in other individuals. For individual 11, permutations resulted in either one new allele (11a) combined with an already known allele, or a combination of two new alleles (11b and 12?). An N-linked glycosylation site is indicated by asterisks. An arrow marks allele 2, which is the inferred ancestral allele. (B) Effects of mustached monkey CD4 polymorphisms on SIV Env mediated cell entry. A heatmap displays the percentage of cells expressing the indicated CD4 variant that were infected by the corresponding SIV Env, averaged across two or three experiments each performed in triplicate (fold-changes of Env infectivity for different CD4 alleles are shown in SI Appendix, Table S6).

To examine whether the D1 domain polymorphisms affected SIV Env-mediated cell entry, we generated a subset of mustached monkey CD4 alleles (alleles 1 to 7) and tested these, together with the cognate CCR5 coreceptor (67), in the single-round infection assay (Fig. 3B). The results showed that all mustached monkey CD4 variants mediated cell entry of the two available SIVmus Envs. However, infectivity of the other SIV Envs varied, with many utilizing only certain alleles (Fig. 3B and SI Appendix, Fig. S6). A comparison of alleles that differed only at position 34 (alleles 1 vs. 3 and 2 vs. 4) showed that a threonine resulted in less-efficient entry than an isoleucine (SI Appendix, Table S6), thus confirming the protective effect of the glycan at position 32. Of all mustached monkey CD4 variants, allele 2 (the inferred ancestral state) was the most permissive, mediating infection even of an SIVdeb Env from De Brazza’s monkeys (Cercopithecus neglectus), the only SIV Env that failed to use human CD4. In contrast, allele 7 was the most restrictive, with only three Envs being able to utilize this variant, two of which were from SIVmus strains (SI Appendix, Table S6). Interestingly, this allele encodes a proline at position 40, which is an Env contact residue on human CD4 and the site of a protective Q40R polymorphism in chimpanzees (55).

Within-Species CD4 Diversity Is Common among African Primates.

To determine how many other primates exhibit CD4 polymorphisms, we used remnant blood or fecal DNA from 124 members of 23 additional African monkey species to amplify and sequence the D1 domain (Table 1). These included samples from 22 Cercopithecine and one Colobine monkey species, 13 of which are known to be naturally SIV-infected (SI Appendix, Table S5). In addition, we extracted D1 domain sequences from published RNA-sequencing (RNA-seq) data of 29 Ugandan red colobus (48) and 4 patas monkeys (Erythrocebus patas) (47), and mined whole-genome databases from 163 African green monkeys belonging to 5 different species (24). Finally, we obtained CD4 sequences from an Angolan colobus (Colobus angolensis), a drill (Mandrillus leucophaeus), and five sooty mangabeys from GenBank. Together with the previously published chimpanzee data (55), we were able to compare D1 domain sequences from 1,101 African primates representing 36 different species.

The various primate species exhibited a striking degree of CD4 diversity, both with respect to the number of species that had CD4 polymorphisms and the number of alleles that were detected. Of 29 species for which 2 or more individuals were genotyped (Table 1), 24 (83%) encoded more than 1 D1 domain variant, with many polymorphic sites located in putative Env–CD4 contact residues (Fig. 4A). In fact, the great majority of SNPs were nonsynonymous, with many species exhibiting no synonymous changes in the D1 domain (Table 1). Greater spot-nosed and mustached monkeys stood out as exhibiting the highest CD4 diversity, with the largest number of D1 variants identified (Table 1). Interestingly, identical polymorphisms were found in members of different species and genera. For example, the Q40R polymorphism, which was shown to affect Env–CD4 interaction in chimpanzees (55), was also found in two Cercopithecus, one Piliocolobus, and one Cercocebus species (Fig. 4A). In addition, two mustached monkey CD4 alleles exhibited a Q40P mutation (Fig. 3A), which was also highly restrictive in the infection assay (Fig. 3B). Although position 40 has not previously been reported to be under positive selection (52, 53), the frequent amino acid replacements at this site in distantly related primate species indicates parallel evolution, most likely driven by pathogenic SIV infections.

Fig. 4.

CD4 diversity in African primate species. (A) D1 domain positions exhibiting intraspecies polymorphisms. For each primate species, the positions of polymorphic residues are shown, with the amino acid residues indicated. Dots indicate identity to the CD4 consensus, while dashes indicate deletions. Contact residues between HIV-1 Env and human CD4 are indicated in red. (B and C) Maximum-likelihood trees of G6PD (B) and CD4 D1 domain (C) nucleotide sequences from different guenon species (color coded). Bootstrap values of >90% are shown; the scale bars indicate 0.001 and 0.004 nucleotide substitutions per site, respectively.

An alignment of D1 domain sequences also revealed that certain CD4 variants were present in multiple closely related primate species. L’Hoest’s and sun-tailed monkeys (genus Allochrocebus) shared certain D1 alleles, as did various African green monkey species, as well as different guenon species (SI Appendix, Fig. S7). To compare the genealogy of CD4 with that of a housekeeping gene, we used available genomic DNA to amplify a 1.6-kb fragment of the glucose-6-phosphate dehydrogenase (G6PD) gene from eight Cercopithecus species for which D1 domain sequences were also available (Fig. 4B). Phylogenetic analysis of these sequences revealed species-specific clustering and an overall topology that was very similar to phylogenies previously reported for these same primate species (68, 69). However, this was not the case for the corresponding CD4 D1 domain sequences (Fig. 4C). Unlike the G6PD sequences, most D1 domain sequences did not cluster according to their species, but were intermixed, with one D1 domain allele shared by six different species (Fig. 4C). These data suggest that the CD4 diversity found in these guenons predates their speciation, indicating trans-specific polymorphism. Since D1 domain sequences from different African green monkey species yielded similar results (SI Appendix, Fig. S8), it appears that CD4 variants have been maintained in various species by long-term balancing selection (70, 71).

Most Primate CD4 Receptors Are Glycan Shielded.

CD4 genotyping of the various primates revealed a total of six potential N-linked glycosylation sites within the D1 domain (SI Appendix, Fig. S7), all of which were predicted to project toward the Env trimer (Fig. 5A). A number of these seemed to be invariant within species, such as the N17 glycan that was present in the great majority of monkey species and the N15 glycan that was present in all eastern and western gorillas (Fig. 5B). Other glycans were polymorphic, such as the N21 glycan that was present in some, but not all, Chlorocebus species. In fact, African green monkeys exhibited three glycosylation sites in the D1 domain, but the maximal number in any one variant was only two (SI Appendix, Fig. S7). Certain glycans were unique to a particular species, including the N66 of chimpanzees and the N39 glycan of mandrills (Fig. 5B). The latter glycan is of interest since it is predicted to protrude directly into the Env binding site and may thus provide an explanation for the ability of SIVmnd1 to infect cells that do not express the CD4 receptor (72). Three guenon and two colobus species appeared to lack D1 domain glycans (Fig. 5B). However, these results are preliminary since in each case only very few individuals were genotyped. Thus, unlike humans, nearly all African primates encode at least one or more N-linked glycosylation sites in their D1 domain, which are positioned to impede SIV infection by creating steric clashes with glycans encoded by their Env glycoproteins.

Fig. 5.

D1 domain glycosylation in primates. (A) PNGS found in the D1 domain of different primate species. PNGS positions are modeled onto the structure of the HIV-1 envelope trimer, with different protomers highlighted in pink, green, and gray, respectively, and the bound human CD4 shown in black (PDB ID code 5U1F). (B) Phylogeny of African primate species adapted from Springer et al. (69). The tree highlights the presence of glycans within each species, which are color coded as in A. Solid and striped squares indicate invariant and polymorphic glycans, respectively. The scale bar indicates estimated primate divergence times as previously reported (69).

Nonhuman Primate Samples.

CD4 and G6PD genotyping was performed using remnant DNA from blood or fecal samples collected previously from both wild-caught and captive primates for molecular epidemiological studies of SIV and Plasmodium infections. All samples were obtained with the approval of the respective Institutional Animal Care and Use Committees. International samples were shipped in compliance with Convention on International Trade in Endangered Species of Wild Fauna and Flora regulations and country-specific import and export permits. Relevant information for all samples is summarized in SI Appendix, Tables S1, S3, and S5.

D1 Domain Genotyping.

CD4 genotyping was carried out as previously described (55). Briefly, CD4 exon 2 (165 bp) and exon 3 (159 bp) were PCR-amplified using primers in adjacent introns resulting in amplicons that were 247 bp and 222 bp in length, respectively. Amplicons were MiSeq sequenced without fragmentation to ensure linkage of variable sites for each exon.

Assembly of CD4 Genes from Bonobo and Gorilla Genomic Databases.

Whole-genome Illumina sequencing reads from bonobos and gorillas were downloaded from the National Center for Biotechnology Information (NCBI) database (56) and from the European Nucleotide Archive (60).

Extraction of D1 Domain Sequences from Monkey RNA-Seq and Genomic Databases.

RNA-seq reads from red colobus monkeys (48) were downloaded from the NCBI Sequence Read Archive (SRA) from BioProject PRJNA413051 and aligned as described in SI Appendix. Reads from African green monkey genome datasets (24) were gathered from NCBI SRA BioProjects PRJNA168472, PRJNA168520, PRJNA168527, PRJNA168522, PRJNA168521, PRJNA362602, and PRJNA407948. RNA-seq reads from patas monkeys, aligned to the Macaca mulatta chromosome 11 (accession NC_027903.1), were provided by the authors (47). Consensus sequence extraction proceeded as for the African green monkeys.

Expression Plasmids and Pseudovirus Stocks.

The construction of an SIVcpz∆env reporter backbone and several SIV Envs were previously reported (55). Codon optimized SIVgor, SIVgsn, SIVmon, and SIVdeb env genes were synthesized and cloned into pcDNA3.1. SIV Env pseudotypes were produced in 293T cells by cotransfection of SIVcpz∆env-GFP and SIV env expression plasmids and titered on TZM-bl cells. Full-length chimpanzee CD4 was mutagenized to obtain the bonobo CD4 clone. Full-length gorilla CD4 was synthesized based on genomic sequences. Chimpanzee CCR5 was mutagenized to match bonobo (XM_034956280) and gorilla (XM_004034016) GenBank sequences. The generation of mustached guenon CD4 and CCR5 expression plasmids has previously been described (67). CD4 polymorphisms were introduced by site-directed mutagenesis into one parental plasmid to obtain multiple D1 domain variants. All CD4 and CCR5 alleles were cloned into pcDNA3.1.

Transient CD4 Transfection Assay.

Env-mediated entry was measured as previously described (55). Briefly, 293T cells were transfected with the various CD4 alleles, along with the species-matched CCR5, infected with SIV Env bearing pseudoviruses, and then analyzed for GFP expression by flow cytometry 2 d post infection. CD4 and CCR5 expression levels were determined at the time of infection using receptor specific antibodies (e.g., SI Appendix, Fig. S9).

Phylogenetic Analyses.

Sequence alignments were trimmed to the same start and end position (297 bp or 99 aa for CD4 D1 domain, 1,584 bp for G6PD). For G6PD sequences, regions with indels causing ambiguity in the alignment were removed. Maximum-likelihood trees were constructed using RaxML v8.2.12 (88) with 100 rapid bootstrap replicates, with model GTRGAMMA for nucleotide sequences and PROTGAMMAJTTF for amino acid sequences.

Statistical Analyses.

A hierarchical Bayesian model was used to assess the effects of CD4 polymorphisms on SIV Env-mediated cell entry (SI Appendix). The effects of the various alleles on the infectivity of each Env, as well as an overall allele effect, was estimated by Markov chain Monte Carlo sampling using Stan (89).