Comprehensive phylogeny of Myrmecocystus honey ants highlights cryptic diversity and infers evolution during aridification of the American Southwest.

The New World ant genus Myrmecocystus Wesmael, 1838 (Formicidae: Formicinae: Lasiini) is endemic to arid and semi-arid habitats of the western United States and Mexico. Several intriguing life history traits have been described for the genus, the best-known of which are replete workers, that store liquified food in their largely expanded crops and are colloquially referred to as "honeypots". Despite their interesting biology and ecological importance for arid ecosystems, the evolutionary history of Myrmecocystus ants is largely unknown and the current taxonomy presents an unsatisfactory systematic framework. We use ultraconserved elements to infer the evolutionary history of Myrmecocystus ants and provide a comprehensive, dated phylogenetic framework that clarifies the molecular systematics within the genus with high statistical support, reveals cryptic diversity, and reconstructs ancestral foraging activity. Using maximum likelihood, Bayesian and species tree approaches on a data set of 134 ingroup specimens (including samples from natural history collections and type material), we recover largely identical topologies that leave the position of only few clades uncertain and cover the intra- and interspecific variation of 28 of the 29 described and six undescribed species. In addition to traditional support values, such as bootstrap and posterior probability, we quantify genealogical concordance to estimate the effects of conflicting evolutionary histories on phylogenetic inference. Our analyses reveal that the current taxonomic classification of the genus is inconsistent with the molecular phylogenetic inference, and we identify cryptic diversity in seven species. Divergence dating suggests that the split between Myrmecocystus and its sister taxon Lasius occurred in the early Miocene. Crown group Myrmecocystus started diversifying about 14.08 Ma ago when the gradual aridification of the southwestern United States and northern Mexico led to formation of the American deserts and to adaptive radiations of many desert taxa.

The Effects of Sampling Location and Predictor Point Estimate Certainty on Posterior Support in Bayesian Phylogeographic Generalized Linear Models.

The use of generalized linear models in Bayesian phylogeography has enabled researchers to simultaneously reconstruct the spatiotemporal history of a virus and quantify the contribution of predictor variables to that process. However, little is known about the sensitivity of this method to the choice of the discrete state partition. Here we investigate this question by analyzing a data set containing 299 sequences of the West Nile virus envelope gene sampled in the United States and fifteen predictors aggregated at four spatial levels. We demonstrate that although the topology of the viral phylogenies was consistent across analyses, support for the predictors depended on the level of aggregation. In particular, we found that the variance of the predictor support metrics was minimized at the most precise level for several predictors and maximized at more sparse levels of aggregation. These results suggest that caution should be taken when partitioning a region into discrete locations to ensure that interpretable, reproducible posterior estimates are obtained. These results also demonstrate why researchers should use the most precise discrete states possible to minimize the posterior variance in such estimates and reveal what truly drives the diffusion of viruses.

Longitudinal analysis of Plasmodium falciparum genetic variation in Turbo, Colombia: implications for malaria control and elimination.

BACKGROUND: Malaria programmes estimate changes in prevalence to evaluate their efficacy. In this study, parasite genetic data was used to explore how the demography of the parasite population can inform about the processes driving variation in prevalence. In particular, how changes in treatment and population movement have affected malaria prevalence in an area with seasonal malaria. METHODS: Samples of Plasmodium falciparum collected over 8 years from a population in Turbo, Colombia were genotyped at nine microsatellite loci and three drug-resistance loci. These data were analysed using several population genetic methods to detect changes in parasite genetic diversity and population structure. In addition, a coalescent-based method was used to estimate substitution rates at the microsatellite loci. RESULTS: The estimated mean microsatellite substitution rates varied between 5.35 × 10(-3) and 3.77 × 10(-2) substitutions/locus/month. Cluster analysis identified six distinct parasite clusters, five of which persisted for the full duration of the study. However, the frequencies of the clusters varied significantly between years, consistent with a small effective population size. CONCLUSIONS: Malaria control programmes can detect re-introductions and changes in transmission using rapidly evolving microsatellite loci. In this population, the steadily decreasing diversity and the relatively constant effective population size suggest that an increase in malaria prevalence from 2004 to 2007 was primarily driven by local rather than imported cases.

Accelerated diversification of nonhuman primate malarias in Southeast Asia: adaptive radiation or geographic speciation?

Although parasitic organisms are found worldwide, the relative importance of host specificity and geographic isolation for parasite speciation has been explored in only a few systems. Here, we study Plasmodium parasites known to infect Asian nonhuman primates, a monophyletic group that includes the lineage leading to the human parasite Plasmodium vivax and several species used as laboratory models in malaria research. We analyze the available data together with new samples from three sympatric primate species from Borneo: The Bornean orangutan and the long-tailed and the pig-tailed macaques. We find several species of malaria parasites, including three putatively new species in this biodiversity hotspot. Among those newly discovered lineages, we report two sympatric parasites in orangutans. We find no differences in the sets of malaria species infecting each macaque species indicating that these species show no host specificity. Finally, phylogenetic analysis of these data suggests that the malaria parasites infecting Southeast Asian macaques and their relatives are speciating three to four times more rapidly than those with other mammalian hosts such as lemurs and African apes. We estimate that these events took place in approximately a 3-4-Ma period. Based on the genetic and phenotypic diversity of the macaque malarias, we hypothesize that the diversification of this group of parasites has been facilitated by the diversity, geographic distributions, and demographic histories of their primate hosts.

The evolutionary history of Plasmodium vivax as inferred from mitochondrial genomes: parasite genetic diversity in the Americas.

Plasmodium vivax is the most prevalent human malaria parasite in the Americas. Previous studies have contrasted the genetic diversity of parasite populations in the Americas with those in Asia and Oceania, concluding that New World populations exhibit low genetic diversity consistent with a recent introduction. Here we used an expanded sample of complete mitochondrial genome sequences to investigate the diversity of P. vivax in the Americas as well as in other continental populations. We show that the diversity of P. vivax in the Americas is comparable to that in Asia and Oceania, and we identify several divergent clades circulating in South America that may have resulted from independent introductions. In particular, we show that several haplotypes sampled in Venezuela and northeastern Brazil belong to a clade that diverged from the other P. vivax lineages at least 30,000 years ago, albeit not necessarily in the Americas. We propose that, unlike in Asia where human migration increases local genetic diversity, the combined effects of the geographical structure and the low incidence of vivax malaria in the Americas has resulted in patterns of low local but high regional genetic diversity. This could explain previous views that P. vivax in the Americas has low genetic diversity because these were based on studies carried out in limited areas. Further elucidation of the complex geographical pattern of P. vivax variation will be important both for diversity assessments of genes encoding candidate vaccine antigens and in the formulation of control and surveillance measures aimed at malaria elimination.

The effect of fluctuating selection on the genealogy at a linked site.

The genealogical consequences of temporally fluctuating selection at linked neutrally-evolving sites are studied using coalescent processes structured by genetic backgrounds. Surprisingly, although between-generation fluctuating selection and within-generation fecundity variance polymorphism lead to indistinguishable allele frequency dynamics at the selected site, I show that these two scenarios affect the genealogical structure of the population in distinctive ways. Whereas coalescence times are elevated at the selected site when heterozygotes have lower within-generation fecundity variance than either homozygote, fluctuating selection typically depresses the depth of the genealogy at tightly-linked sites. More importantly, these results indicate that fluctuating selection will have a characteristic signature combining an excess of selected variation at the affected site and a slight reduction in neutral variation at tightly linked sites. This observation suggests that it may be possible to distinguish balanced polymorphisms maintained by environmental variation from those maintained by heterozygote advantage.

TbISWI regulates multiple polymerase I (Pol I)-transcribed loci and is present at Pol II transcription boundaries in Trypanosoma brucei.

The unicellular eukaryote Trypanosoma brucei is unusual in having very little transcriptional control. The bulk of the T. brucei genome is constitutively transcribed by RNA polymerase II (Pol II) as extensive polycistronic transcription units. Exceptions to this rule include several RNA Pol I transcription units such as the VSG expression sites (ESs), which are mono-allelically expressed. TbISWI, a member of the SWI2/SNF2 related chromatin remodeling ATPases, plays a role in repression of Pol I-transcribed ESs in both bloodstream- and procyclic-form T. brucei. We show that TbISWI binds both active and silent ESs but is depleted from the ES promoters themselves. TbISWI knockdown results in an increase in VSG transcripts from the silent VSG ESs. In addition to its role in the repression of the silent ESs, TbISWI also contributes to the downregulation of the Pol I-transcribed procyclin loci, as well as nontranscribed VSG basic copy arrays and minichromosomes. We also show that TbISWI is enriched at a number of strand switch regions which form the boundaries between Pol II transcription units. These strand switch regions are the presumed sites of Pol II transcription initiation and termination and are enriched in modified histones and histone variants. Our results indicate that TbISWI is a versatile chromatin remodeler that regulates transcription at multiple Pol I loci and is particularly abundant at many Pol II transcription boundaries in T. brucei.

A coalescent dual process in a Moran model with genic selection, and the lambda coalescent limit.

The genealogical structure of neutral populations in which reproductive success is highly-skewed has been the subject of many recent studies. Here we derive a coalescent dual process for a related class of continuous-time Moran models with viability selection. In these models, individuals can give birth to multiple offspring whose survival depends on both the parental genotype and the brood size. This extends the dual process construction for a multi-type Moran model with genic selection described in Etheridge and Griffiths (2009). We show that in the limit of infinite population size the non-neutral Moran models converge to a Markov jump process which we call the lamda-Fleming-Viot process with viability selection and we derive a coalescent dual for this process directly from the generator and as a limit from the Moran models. The dual is a branching-coalescing process similar to the Ancestral Selection Graph which follows the typed ancestry of genes backwards in time with real and virtual lineages. As an application, the transition functions of the non-neutral Moran and lamda-coalescent models are expressed as mixtures of the transition functions of the dual process.

The genealogical consequences of fecundity variance polymorphism.

The genealogical consequences of within-generation fecundity variance polymorphism are studied using coalescent processes structured by genetic backgrounds. I show that these processes have three distinctive features. The first is that the coalescent rates within backgrounds are not jointly proportional to the infinitesimal variance, but instead depend only on the frequencies and traits of genotypes containing each allele. Second, the coalescent processes at unlinked loci are correlated with the genealogy at the selected locus; i.e., fecundity variance polymorphism has a genomewide impact on genealogies. Third, in diploid models, there are infinitely many combinations of fecundity distributions that have the same diffusion approximation but distinct coalescent processes; i.e., in this class of models, ancestral processes and allele frequency dynamics are not in one-to-one correspondence. Similar properties are expected to hold in models that allow for heritable variation in other traits that affect the coalescent effective population size, such as sex ratio or fecundity and survival schedules.

Telomeric expression sites are highly conserved in Trypanosoma brucei.

Subtelomeric regions are often under-represented in genome sequences of eukaryotes. One of the best known examples of the use of telomere proximity for adaptive purposes are the bloodstream expression sites (BESs) of the African trypanosome Trypanosoma brucei. To enhance our understanding of BES structure and function in host adaptation and immune evasion, the BES repertoire from the Lister 427 strain of T. brucei were independently tagged and sequenced. BESs are polymorphic in size and structure but reveal a surprisingly conserved architecture in the context of extensive recombination. Very small BESs do exist and many functioning BESs do not contain the full complement of expression site associated genes (ESAGs). The consequences of duplicated or missing ESAGs, including ESAG9, a newly named ESAG12, and additional variant surface glycoprotein genes (VSGs) were evaluated by functional assays after BESs were tagged with a drug-resistance gene. Phylogenetic analysis of constituent ESAG families suggests that BESs are sequence mosaics and that extensive recombination has shaped the evolution of the BES repertoire. This work opens important perspectives in understanding the molecular mechanisms of antigenic variation, a widely used strategy for immune evasion in pathogens, and telomere biology.

Isolation and analysis of the genetic diversity of repertoires of VSG expression site containing telomeres from Trypanosoma brucei gambiense, T. b. brucei and T. equiperdum.

BACKGROUND: African trypanosomes (including Trypanosoma brucei) are unicellular parasites which multiply in the mammalian bloodstream. T. brucei has about twenty telomeric bloodstream form Variant Surface Glycoprotein (VSG) expression sites (BESs), of which one is expressed at a time in a mutually exclusive fashion. BESs are polycistronic transcription units, containing a variety of families of expression site associated genes (ESAGs) in addition to the telomeric VSG. These polymorphic ESAG families are thought to play a role in parasite-host adaptation, and it has been proposed that ESAG diversity might be related to host range. Analysis of the genetic diversity of these telomeric gene families has been confounded by the underrepresentation of telomeric sequences in standard libraries. We have previously developed a method to selectively isolate sets of trypanosome BES containing telomeres using Transformation associated recombination (TAR) cloning in yeast. RESULTS: Here we describe the isolation of repertoires of BES containing telomeres from three trypanosome subspecies: Trypanosoma brucei gambiense DAL 972 (causative agent of West-African trypanosomiasis), T. b. brucei EATRO 2340 (a nonhuman infective strain) and T. equiperdum STIB 818 (which causes a sexually transmitted disease in equines). We have sequenced and analysed the genetic diversity at four BES loci (BES promoter region, ESAG6, ESAG5 and ESAG2) from these three trypanosome BES repertoires. CONCLUSION: With the exception of ESAG2, the BES sequence repertoires derived from T. b. gambiense are both less diverse than and nearly reciprocally monophyletic relative to those from T. b. brucei and T. equiperdum. Furthermore, although we find evidence for adaptive evolution in all three ESAG repertoires in T. b. brucei and T. equiperdum, only ESAG2 appears to be under diversifying selection in T. b. gambiense. This low level of variation in the T. b. gambiense BES sequence repertoires is consistent both with the relatively narrow host range of this subspecies and its apparent long-term clonality. However, our data does not show a clear correlation between size of trypanosome host range and either number of BESs or extent of ESAG genetic diversity.

Environmental variation, fluctuating selection and genetic drift in subdivided populations.

Although there have many studies of the population genetical consequences of environmental variation, little is known about the combined effects of genetic drift and fluctuating selection in structured populations. Here we use diffusion theory to investigate the effects of temporally and spatially varying selection on a population of haploid individuals subdivided into a large number of demes. Using a perturbation method for processes with multiple time scales, we show that as the number of demes tends to infinity, the overall frequency converges to a diffusion process that is also the diffusion approximation for a finite, panmictic population subject to temporally fluctuating selection. We find that the coefficients of this process have a complicated dependence on deme size and migration rate, and that changes in these demographic parameters can determine both the balance between the dispersive and stabilizing effects of environmental variation and whether selection favors alleles with lower or higher fitness variance.

Gene flow and the genealogical history of Heliconius heurippa.

BACKGROUND: The neotropical butterfly Heliconius heurippa has a hybrid colour pattern, which also contributes to reproductive isolation, making it a likely example of hybrid speciation. Here we used phylogenetic and coalescent-based analyses of multilocus sequence data to investigate the origin of H. heurippa. RESULTS: We sequenced a mitochondrial region (CoI and CoII), a sex-linked locus (Tpi) and two autosomal loci (w and sd) from H. heurippa and the putative parental species, H. cydno and H. melpomene. These were analysed in combination with data from two previously sequenced autosomal loci, Dll and Inv. H. heurippa was monophyletic at mtDNA and Tpi, but showed a shared distribution of alleles derived from both parental lineages at all four autosomal loci. Estimates of genetic differentiation showed that H. heurippa is closer to H. cydno at mtDNA and three autosomal loci, intermediate at Tpi, and closer to H. melpomene at Dll. Using coalescent simulations with the Isolation-Migration model (IM), we attempted to establish the incidence of gene flow in the origin of H. heurippa. This analysis suggested that ongoing introgression is frequent between all three species and variable in extent between loci. CONCLUSION: Introgression, which is a necessary precursor of hybrid speciation, seems to have also blurred the coalescent history of these species. The origin of Heliconius heurippa may have been restricted to introgression of few colour pattern genes from H. melpomene into the H. cydno genome, with little evidence of genomic mosaicism.

Switching trypanosome coats: what's in the wardrobe?

The African trypanosome Trypanosoma brucei is best known for its extraordinarily sophisticated antigenic variation of a protective variant surface glycoprotein (VSG) coat. T. brucei has >1000 VSG genes and pseudogenes, of which one is transcribed at a time from one of multiple telomeric VSG expression sites. Switching the active VSG gene can involve DNA rearrangements replacing the old VSG with a new one, or alternatively transcriptional control. The astonishing revelation from the T. brucei genome sequence is that <7% of the sequenced VSGs seem to have fully functional coding regions. This preponderance of pseudogenes in the VSG gene repertoire will necessitate a rethink of how antigenic variation in African trypanosomes operates.

HIV-1 intra-subtype superinfection rates: estimates using a structured coalescent with recombination.

Because recombination between a pair of viral genomes can occur only when both viruses are present in the same host, genealogical evidence of recombination is influenced by the rate of viral migration between infected hosts. If superinfection is rare, then recombining viral genomes will usually be more closely related to each other than to viral genomes present in other hosts and the impact of recombination on the genealogy of a sample of viruses from different hosts may be weak. We used this relationship to estimate the relative rate of intra-subtype HIV-1 superinfection in six urban populations. Comparisons of the population recombination rates estimated from the HIV-1 sequence data with population recombination rates estimated for sets of sequences simulated using a structured coalescent process suggest that intra-subtype superinfection rates in all but one of these populations may be as high as 15% of the corresponding infection rate. However, we caution that these estimates might be upwardly biased if variation in contact and mixing rates between infected hosts causes viral lineages to be concentrated in groups with higher than average superinfection rates.

Sperm competition and the dynamics of X chromosome drive: stability and extinction.

Several empirical studies of sperm competition in populations polymorphic for a driving X chromosome have revealed that Sex-ratio males (those carrying a driving X) are at a disadvantage relative to Standard males. Because the frequency of the driving X chromosome determines the population-level sex ratio and thus alters male and female mating rates, the evolutionary consequences of sperm competition for sex chromosome meiotic drive are subtle. As the SR allele increases in frequency, the ratio of females to males also increases, causing an increase in the male mating rate and a decrease in the female mating rate. While the former change may exacerbate the disadvantage of Sex-ratio males during sperm competition, the latter change decreases the incidence of sperm competition within the population. We analyze a model of the effects of sperm competition on a driving X chromosome and show that these opposing trends in male and female mating rates can result in two coexisting locally stable equilibria, one corresponding to a balanced polymorphism of the SR and ST alleles and the second to fixation of the ST allele. Stochastic fluctuations of either the population sex ratio or the SR frequency can then drive the population away from the balanced polymorphism and into the basin of attraction for the second equilibrium, resulting in fixation of the SR allele and extinction of the population.