Pervasive mimicry in flight behavior among aposematic butterflies.

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Parallel evolution in mosquito vectors - a duplicated esterase locus is associated with resistance to pirimiphos-methyl in An. gambiae.

The primary control methods for the African malaria mosquito, Anopheles gambiae, are based on insecticidal interventions. Emerging resistance to these compounds is therefore of major concern to malaria control programmes. The organophosphate, pirimiphos-methyl, is a relatively new chemical in the vector control armoury but is now widely used in indoor residual spray campaigns. Whilst generally effective, phenotypic resistance has developed in some areas in malaria vectors. Here, we used a population genomic approach to identify novel mechanisms of resistance to pirimiphos-methyl in Anopheles gambiae s.l mosquitoes. In multiple populations, we found large and repeated signals of selection at a locus containing a cluster of detoxification enzymes, some of whose orthologs are known to confer resistance to organophosphates in Culex pipiens. Close examination revealed a pair of alpha-esterases, Coeae1f and Coeae2f, and a complex and diverse pattern of haplotypes under selection in An. gambiae, An. coluzzii and An. arabiensis. As in Cx. pipiens, copy number variants have arisen at this locus. We used diplotype clustering to examine whether these signals arise from parallel evolution or adaptive introgression. Using whole-genome sequenced phenotyped samples, we found that in West Africa, a copy number variant in Anopheles gambiae is associated with resistance to pirimiphos-methyl. Overall, we demonstrate a striking example of contemporary parallel evolution which has important implications for malaria control programmes.

The trichome pattern diversity of Cardamine shares genetic mechanisms with Arabidopsis but differs in environmental drivers.

Natural variation in trichome pattern (amount and distribution) is prominent among populations of many angiosperms. However, the degree of parallelism in the genetic mechanisms underlying this diversity and its environmental drivers in different species remain unclear. To address these questions, we analyzed the genomic and environmental bases of leaf trichome pattern diversity in Cardamine hirsuta, a relative of Arabidopsis (Arabidopsis thaliana). We characterized 123 wild accessions for their genomic diversity, leaf trichome patterns at different temperatures, and environmental adjustments. Nucleotide diversities and biogeographical distribution models identified two major genetic lineages with distinct demographic and adaptive histories. Additionally, C. hirsuta showed substantial variation in trichome pattern and plasticity to temperature. Trichome amount in C. hirsuta correlated positively with spring precipitation but negatively with temperature, which is opposite to climatic patterns in A. thaliana. Contrastingly, genetic analysis of C. hirsuta glabrous accessions indicated that, like for A. thaliana, glabrousness is caused by null mutations in ChGLABRA1 (ChGL1). Phenotypic genome-wide association studies (GWAS) further identified a ChGL1 haplogroup associated with low trichome density and ChGL1 expression. Therefore, a ChGL1 series of null and partial loss-of-function alleles accounts for the parallel evolution of leaf trichome pattern in C. hirsuta and A. thaliana. Finally, GWAS also detected other candidate genes (e.g. ChETC3, ChCLE17) that might affect trichome pattern. Accordingly, the evolution of this trait in C. hirsuta and A. thaliana shows partially conserved genetic mechanisms but is likely involved in adaptation to different environments.

Unpredictable repeatability in molecular evolution.

The extent of parallel evolution at the genotypic level is quantitatively linked to the distribution of beneficial fitness effects (DBFE) of mutations. The standard view, based on light-tailed distributions (i.e., distributions with finite moments), is that the probability of parallel evolution in duplicate populations is inversely proportional to the number of available mutations and, moreover, that the DBFE is sufficient to determine the probability when the number of available mutations is large. Here, we show that when the DBFE is heavy-tailed, as found in several recent experiments, these expectations are defied. The probability of parallel evolution decays anomalously slowly in the number of mutations or even becomes independent of it, implying higher repeatability of evolution. At the same time, the probability of parallel evolution is non-self-averaging-that is, it does not converge to its mean value, even when a large number of mutations are involved. This behavior arises because the evolutionary process is dominated by only a few mutations of high weight. Consequently, the probability varies widely across systems with the same DBFE. Contrary to the standard view, the DBFE is no longer sufficient to determine the extent of parallel evolution, making it much less predictable. We illustrate these ideas theoretically and through analysis of empirical data on antibiotic-resistance evolution.

Parallel adaptation to lower altitudes is associated with enhanced plasticity in Heliosperma pusillum (Caryophyllaceae).

High levels of phenotypic plasticity are thought to be inherently costly in stable or extreme environments, but enhanced plasticity may evolve as a response to new environments and foster novel phenotypes. Heliosperma pusillum forms glabrous alpine and pubescent montane ecotypes that diverged recurrently and polytopically (parallel evolution) and can serve as evolutionary replicates. The specific alpine and montane localities are characterized by distinct temperature conditions, available moisture, and light. Noteworthy, the ecotypes show a home-site fitness advantage in reciprocal transplantations. To disentangle the relative contribution of constitutive versus plastic gene expression to altitudinal divergence, we analyze the transcriptomic profiles of two parallely evolved ecotype pairs, grown in reciprocal transplantations at native altitudinal sites. In this incipient stage of divergence, only a minor proportion of genes appear constitutively differentially expressed between the ecotypes in both pairs, regardless of the growing environment. Both derived, montane populations bear comparatively higher plasticity of gene expression than the alpine populations. Genes that change expression plastically or constitutively underlie similar ecologically relevant pathways, related to response to drought and trichome formation. Other relevant processes, such as photosynthesis, rely mainly on plastic changes. The enhanced plasticity consistently observed in the montane ecotype likely evolved as a response to the newly colonized, drier, and warmer niche. We report a striking parallelism of directional changes in gene expression plasticity. Thus, plasticity appears to be a key mechanism shaping the initial stages of phenotypic evolution, likely fostering adaptation to novel environments.

Predictable and Divergent Change in the Multivariate P Matrix during Parallel Adaptation.

AbstractAdaptation to replicated environmental conditions can be remarkably predictable, suggesting that parallel evolution may be a common feature of adaptive radiation. An open question, however, is how phenotypic variation itself evolves during repeated adaptation. Here, we use a dataset of morphological measurements from 35 populations of threespine stickleback, consisting of 16 parapatric lake-stream pairs and three marine populations, to understand how phenotypic variation has evolved during transitions from marine to freshwater environments and during subsequent diversification across the lake-stream boundary. We find statistical support for divergent phenotypic covariance (P) across populations, with most diversification of P occurring among freshwater populations. Despite a close correspondence between within-population phenotypic variation and among-population divergence, we find that variation in P is unrelated to total variation in population means across the set of populations. For lake-stream pairs, we find that theoretical predictions for microevolutionary change can explain more than 30% of divergence in P matrices across the habitat boundary. Together, our results indicate that divergence in variance structure occurs primarily in dimensions of trait space with low phenotypic integration, correlated with disparate lake and stream environments. Our findings illustrate how conserved and divergent features of multivariate variation can underlie adaptive radiation.

Parallel evolution despite low genetic diversity in three-spined sticklebacks.

When populations repeatedly adapt to similar environments they can evolve similar phenotypes based on shared genetic mechanisms (parallel evolution). The likelihood of parallel evolution is affected by demographic history, as it depends on the standing genetic variation of the source population. The three-spined stickleback (Gasterosteus aculeatus) repeatedly colonized and adapted to brackish and freshwater. Most parallel evolution studies in G. aculeatus were conducted at high latitudes, where freshwater populations maintain connectivity to the source marine populations. Here, we analysed southern and northern European marine and freshwater populations to test two hypotheses. First, that southern European freshwater populations (which currently lack connection to marine populations) lost genetic diversity due to bottlenecks and inbreeding compared to their northern counterparts. Second, that the degree of genetic parallelism is higher among northern than southern European freshwater populations, as the latter have been subjected to strong drift due to isolation. The results show that southern populations exhibit lower genetic diversity but a higher degree of genetic parallelism than northern populations. Hence, they confirm the hypothesis that southern populations have lost genetic diversity, but this loss probably happened after they had already adapted to freshwater conditions, explaining the high degree of genetic parallelism in the south.

Comparative transcriptomics revealed parallel evolution and innovation of photosymbiosis molecular mechanisms in a marine bivalve.

Photosymbioses between heterotrophic hosts and autotrophic symbionts are evolutionarily prevalent and ecologically significant. However, the molecular mechanisms behind such symbioses remain less elucidated, which hinders our understanding of their origin and adaptive evolution. This study compared gene expression patterns in a photosymbiotic bivalve (Fragum sueziense) and a closely related non-symbiotic species (Trigoniocardia granifera) under different light conditions to detect potential molecular pathways involved in mollusc photosymbiosis. We discovered that the presence of algal symbionts greatly impacted host gene expression in symbiont-containing tissues. We found that the host immune functions were suppressed under normal light compared with those in the dark. In addition, we found that cilia in the symbiont-containing tissues play important roles in symbiont regulation or photoreception. Interestingly, many potential photosymbiosis genes could not be annotated or do not exhibit orthologues in T. granifera transcriptomes, indicating unique molecular functions in photosymbiotic bivalves. Overall, we found both novel and known molecular mechanisms involved in animal-algal photosymbiosis within bivalves. Given that many of the molecular pathways are shared among distantly related host lineages, such as molluscs and cnidarians, it indicates that parallel and/or convergent evolution is instrumental in shaping host-symbiont interactions and responses in these organisms.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs.

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization.

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

Comparative histology of abscission zones reveals the extent of convergence and divergence in seed shattering in weedy and cultivated rice.

The modification of seed shattering has been a recurring theme in rice evolution. The wild ancestor of cultivated rice disperses its seeds, but reduced shattering was selected during multiple domestication events to facilitate harvesting. Conversely, selection for increased shattering occurred during the evolution of weedy rice, a weed invading cultivated rice fields that has originated multiple times from domesticated ancestors. Shattering requires formation of a tissue known as the abscission zone (AZ), but how the AZ has been modified throughout rice evolution is unclear. We quantitatively characterized the AZ characteristics of relative length, discontinuity, and intensity in 86 cultivated and weedy rice accessions. We reconstructed AZ evolutionary trajectories and determined the degree of convergence among different cultivated varieties and among independent weedy rice populations. AZ relative length emerged as the best feature to distinguish high and low shattering rice. Cultivated varieties differed in average AZ morphology, revealing lack of convergence in how shattering reduction was achieved during domestication. In contrast, weedy rice populations typically converged on complete AZs, irrespective of origin. By examining AZ population-level morphology, our study reveals its evolutionary plasticity, and suggests that the genetic potential to modify the ecologically and agronomically important trait of shattering is plentiful in rice lineages.

Parallel shifts in flight-height associated with altitude across incipient Heliconius species.

Vertical gradients in microclimate, resource availability, and interspecific interactions are thought to underly stratification patterns in tropical insect communities. However, only a few studies have explored the adaptive significance of vertical space use during the early stages of reproductive isolation. We analysed flight-height variation across speciation events in Heliconius butterflies, representing parallel colonizations of high-altitude forest. We measured flight-height in wild H. erato venus and H. chestertonii, parapatric lowland and mountain specialists, respectively, and found that H. chestertonii consistently flies at a lower height. By comparing our data to previously published results for the ecologically equivalent H. e. cyrbia (lowland) and H. himera (high altitude), we found that the species flying closest to the ground are those that recently colonized high-altitude forests. We show that these repeated trends largely result from shared patterns of ecological selection producing parallel trait-shifts in H. himera and H. chestertonii. Although our results imply a signature of local adaptation, we did not find an association between resource distribution and flight-height in H. e. venus and H. chestertonii. We discuss how this pattern may be explained by variations in forest structure and microclimate. Overall, our findings underscore the importance of behavioural adjustments during early divergence mediated by altitude-shifts.

Introgression underlies phylogenetic uncertainty but not parallel plumage evolution in a recent songbird radiation.

Instances of parallel phenotypic evolution offer great opportunities to understand the evolutionary processes underlying phenotypic changes. However, confirming parallel phenotypic evolution and studying its causes requires a robust phylogenetic framework. One such example is the "black-and-white wagtails", a group of five species in the songbird genus Motacilla: one species, Motacilla alba, shows wide intra-specific plumage variation, while the four others form two pairs of very similar-looking species (M. aguimp + M. samveasnae and M. grandis + M. maderaspatensis, respectively). However, the two species in each of these pairs were not recovered as sisters in previous phylogenetic inferences. Their relationships varied depending on the markers used, suggesting that gene tree heterogeneity might have hampered accurate phylogenetic inference. Here, we use whole genome resequencing data to explore the phylogenetic relationships within this group, with a special emphasis on characterizing the extent of gene tree heterogeneity and its underlying causes. We first used multispecies coalescent methods to generate a "complete evidence" phylogenetic hypothesis based on genome-wide variants, while accounting for incomplete lineage sorting (ILS) and introgression. We then investigated the variation in phylogenetic signal across the genome, to quantify the extent of discordance across genomic regions, and test its underlying causes. We found that wagtail genomes are mosaics of regions supporting variable genealogies, because of ILS and inter-specific introgression. The most common topology across the genome, supporting M. alba and M. aguimp as sister species, appears to be influenced by ancient introgression. Additionally, we inferred another ancient introgression event, between M. alba and M. grandis. By combining results from multiple analyses, we propose a phylogenetic network for the black-and-white wagtails that confirms that similar phenotypes evolved in non-sister lineages, supporting parallel plumage evolution. Furthermore, the inferred reticulations do not connect species with similar plumage coloration, suggesting that introgression does not underlie parallel plumage evolution in this group. Our results demonstrate the importance of investigation of genome-wide patterns of gene tree heterogeneity to help understanding the mechanisms underlying phenotypic evolution.

Similar enzymatic functions in distinct bioluminescence systems: evolutionary recruitment of sulfotransferases in ostracod light organs.

Genes from ancient families are sometimes involved in the convergent evolutionary origins of similar traits, even across vast phylogenetic distances. Sulfotransferases are an ancient family of enzymes that transfer sulfate from a donor to a wide variety of substrates, including probable roles in some bioluminescence systems. Here, we demonstrate multiple sulfotransferases, highly expressed in light organs of the bioluminescent ostracod Vargula tsujii, transfer sulfate in vitro to the luciferin substrate, vargulin. We find luciferin sulfotransferases (LSTs) of ostracods are not orthologous to known LSTs of fireflies or sea pansies; animals with distinct and convergently evolved bioluminescence systems compared to ostracods. Therefore, distantly related sulfotransferases were independently recruited at least three times, leading to parallel evolution of luciferin metabolism in three highly diverged organisms. Reuse of homologous genes is surprising in these bioluminescence systems because the other components, including luciferins and luciferases, are completely distinct. Whether convergently evolved traits incorporate ancient genes with similar functions or instead use distinct, often newer, genes may be constrained by how many genetic solutions exist for a particular function. When fewer solutions exist, as in genetic sulfation of small molecules, evolution may be more constrained to use the same genes time and again.

A switch to feeding on cycads generates parallel accelerated evolution of toxin tolerance in two clades of Eumaeus caterpillars (Lepidoptera: Lycaenidae).

We assembled a complete reference genome of Eumaeus atala, an aposematic cycad-eating hairstreak butterfly that suffered near extinction in the United States in the last century. Based on an analysis of genomic sequences of Eumaeus and 19 representative genera, the closest relatives of Eumaeus are Theorema and Mithras We report natural history information for Eumaeus, Theorema, and Mithras Using genomic sequences for each species of Eumaeus, Theorema, and Mithras (and three outgroups), we trace the evolution of cycad feeding, coloration, gregarious behavior, and other traits. The switch to feeding on cycads and to conspicuous coloration was accompanied by little genomic change. Soon after its origin, Eumaeus split into two fast evolving lineages, instead of forming a clump of close relatives in the phylogenetic tree. Significant overlap of the fast evolving proteins in both clades indicates parallel evolution. The functions of the fast evolving proteins suggest that the caterpillars developed tolerance to cycad toxins with a range of mechanisms including autophagy of damaged cells, removal of cell debris by macrophages, and more active cell proliferation.

Genomic evolution of antibiotic resistance is contingent on genetic background following a long-term experiment with Escherichia coli.

Antibiotic resistance is a growing health concern. Efforts to control resistance would benefit from an improved ability to forecast when and how it will evolve. Epistatic interactions between mutations can promote divergent evolutionary trajectories, which complicates our ability to predict evolution. We recently showed that differences between genetic backgrounds can lead to idiosyncratic responses in the evolvability of phenotypic resistance, even among closely related Escherichia coli strains. In this study, we examined whether a strain's genetic background also influences the genotypic evolution of resistance. Do lineages founded by different genotypes take parallel or divergent mutational paths to achieve their evolved resistance states? We addressed this question by sequencing the complete genomes of antibiotic-resistant clones that evolved from several different genetic starting points during our earlier experiments. We first validated our statistical approach by quantifying the specificity of genomic evolution with respect to antibiotic treatment. As expected, mutations in particular genes were strongly associated with each drug. Then, we determined that replicate lines evolved from the same founding genotypes had more parallel mutations at the gene level than lines evolved from different founding genotypes, although these effects were more subtle than those showing antibiotic specificity. Taken together with our previous work, we conclude that historical contingency can alter both genotypic and phenotypic pathways to antibiotic resistance.

Sex-specific ornament evolution is a consistent feature of climatic adaptation across space and time in dragonflies.

Adaptation to different climates fuels the origins and maintenance of biodiversity. Detailing how organisms optimize fitness for their local climates is therefore an essential goal in biology. Although we increasingly understand how survival-related traits evolve as organisms adapt to climatic conditions, it is unclear whether organisms also optimize traits that coordinate mating between the sexes. Here, we show that dragonflies consistently adapt to warmer climates across space and time by evolving less male melanin ornamentation-a mating-related trait that also absorbs solar radiation and heats individuals above ambient temperatures. Continent-wide macroevolutionary analyses reveal that species inhabiting warmer climates evolve less male ornamentation. Community-science observations across 10 species indicate that populations adapt to warmer parts of species' ranges through microevolution of smaller male ornaments. Observations from 2005 to 2019 detail that contemporary selective pressures oppose male ornaments in warmer years; and our climate-warming projections predict further decreases by 2070. Conversely, our analyses show that female ornamentation responds idiosyncratically to temperature across space and time, indicating the sexes evolve in different ways to meet the demands of the local climate. Overall, these macro- and microevolutionary findings demonstrate that organisms predictably optimize their mating-related traits for the climate just as they do their survival-related traits.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds.

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Parallel evolution of UbiA superfamily proteins into aromatic O-prenyltransferases in plants.

Plants produce ∼300 aromatic compounds enzymatically linked to prenyl side chains via C-O bonds. These O-prenylated aromatic compounds have been found in taxonomically distant plant taxa, with some of them being beneficial or detrimental to human health. Although their O-prenyl moieties often play crucial roles in the biological activities of these compounds, no plant gene encoding an aromatic O-prenyltransferase (O-PT) has been isolated to date. This study describes the isolation of an aromatic O-PT gene, CpPT1, belonging to the UbiA superfamily, from grapefruit (Citrus × paradisi, Rutaceae). This gene was shown responsible for the biosynthesis of O-prenylated coumarin derivatives that alter drug pharmacokinetics in the human body. Another coumarin O-PT gene encoding a protein of the same family was identified in Angelica keiskei, an apiaceous medicinal plant containing pharmaceutically active O-prenylated coumarins. Phylogenetic analysis of these O-PTs suggested that aromatic O-prenylation activity evolved independently from the same ancestral gene in these distant plant taxa. These findings shed light on understanding the evolution of plant secondary (specialized) metabolites via the UbiA superfamily.

CD4 receptor diversity represents an ancient protection mechanism against primate lentiviruses.

Infection with human and simian immunodeficiency viruses (HIV/SIV) requires binding of the viral envelope glycoprotein (Env) to the host protein CD4 on the surface of immune cells. Although invariant in humans, the Env binding domain of the chimpanzee CD4 is highly polymorphic, with nine coding variants circulating in wild populations. Here, we show that within-species CD4 diversity is not unique to chimpanzees but found in many African primate species. Characterizing the outermost (D1) domain of the CD4 protein in over 500 monkeys and apes, we found polymorphic residues in 24 of 29 primate species, with as many as 11 different coding variants identified within a single species. D1 domain amino acid replacements affected SIV Env-mediated cell entry in a single-round infection assay, restricting infection in a strain- and allele-specific fashion. Several identical CD4 polymorphisms, including the addition of N-linked glycosylation sites, were found in primate species from different genera, providing striking examples of parallel evolution. Moreover, seven different guenons (Cercopithecus spp.) shared multiple distinct D1 domain variants, pointing to long-term trans-specific polymorphism. These data indicate that the HIV/SIV Env binding region of the primate CD4 protein is highly variable, both within and between species, and suggest that this diversity has been maintained by balancing selection for millions of years, at least in part to confer protection against primate lentiviruses. Although long-term SIV-infected species have evolved specific mechanisms to avoid disease progression, primate lentiviruses are intrinsically pathogenic and have left their mark on the host genome.