The evolution of phylogeographic data sets.

Empirical phylogeographic studies have progressively sampled greater numbers of loci over time, in part motivated by theoretical papers showing that estimates of key demographic parameters improve as the number of loci increases. Recently, next-generation sequencing has been applied to questions about organismal history, with the promise of revolutionizing the field. However, no systematic assessment of how phylogeographic data sets have changed over time with respect to overall size and information content has been performed. Here, we quantify the changing nature of these genetic data sets over the past 20 years, focusing on papers published in Molecular Ecology. We found that the number of independent loci, the total number of alleles sampled and the total number of single nucleotide polymorphisms (SNPs) per data set has improved over time, with particularly dramatic increases within the past 5 years. Interestingly, uniparentally inherited organellar markers (e.g. animal mitochondrial and plant chloroplast DNA) continue to represent an important component of phylogeographic data. Single-species studies (cf. comparative studies) that focus on vertebrates (particularly fish and to some extent, birds) represent the gold standard of phylogeographic data collection. Based on the current trajectory seen in our survey data, forecast modelling indicates that the median number of SNPs per data set for studies published by the end of the year 2016 may approach ~20,000. This survey provides baseline information for understanding the evolution of phylogeographic data sets and underscores the fact that development of analytical methods for handling very large genetic data sets will be critical for facilitating growth of the field.

Wolbachia association with the tsetse fly, Glossina fuscipes fuscipes, reveals high levels of genetic diversity and complex evolutionary dynamics.

BACKGROUND: Wolbachia pipientis, a diverse group of α-proteobacteria, can alter arthropod host reproduction and confer a reproductive advantage to Wolbachia-infected females (cytoplasmic incompatibility (CI)). This advantage can alter host population genetics because Wolbachia-infected females produce more offspring with their own mitochondrial DNA (mtDNA) haplotypes than uninfected females. Thus, these host haplotypes become common or fixed (selective sweep). Although simulations suggest that for a CI-mediated sweep to occur, there must be a transient phase with repeated initial infections of multiple individual hosts by different Wolbachia strains, this has not been observed empirically. Wolbachia has been found in the tsetse fly, Glossina fuscipes fuscipes, but it is not limited to a single host haplotype, suggesting that CI did not impact its population structure. However, host population genetic differentiation could have been generated if multiple Wolbachia strains interacted in some populations. Here, we investigated Wolbachia genetic variation in G. f. fuscipes populations of known host genetic composition in Uganda. We tested for the presence of multiple Wolbachia strains using Multi-Locus Sequence Typing (MLST) and for an association between geographic region and host mtDNA haplotype using Wolbachia DNA sequence from a variable locus, groEL (heat shock protein 60). RESULTS: MLST demonstrated that some G. f. fuscipes carry Wolbachia strains from two lineages. GroEL revealed high levels of sequence diversity within and between individuals (Haplotype diversity = 0.945). We found Wolbachia associated with 26 host mtDNA haplotypes, an unprecedented result. We observed a geographical association of one Wolbachia lineage with southern host mtDNA haplotypes, but it was non-significant (p = 0.16). Though most Wolbachia-infected host haplotypes were those found in the contact region between host mtDNA groups, this association was non-significant (p = 0.17). CONCLUSIONS: High Wolbachia sequence diversity and the association of Wolbachia with multiple host haplotypes suggest that different Wolbachia strains infected G. f. fuscipes multiple times independently. We suggest that these observations reflect a transient phase in Wolbachia evolution that is influenced by the long gestation and low reproductive output of tsetse. Although G. f. fuscipes is superinfected with Wolbachia, our data does not support that bidirectional CI has influenced host genetic diversity in Uganda.

Implications of microfauna-host interactions for trypanosome transmission dynamics in Glossina fuscipes fuscipes in Uganda.

Tsetse flies (Diptera: Glossinidae) are vectors for African trypanosomes (Euglenozoa: kinetoplastida), protozoan parasites that cause African trypanosomiasis in humans (HAT) and nagana in livestock. In addition to trypanosomes, two symbiotic bacteria (Wigglesworthia glossinidia and Sodalis glossinidius) and two parasitic microbes, Wolbachia and a salivary gland hypertrophy virus (SGHV), have been described in tsetse. Here we determined the prevalence of and coinfection dynamics between Wolbachia, trypanosomes, and SGHV in Glossina fuscipes fuscipes in Uganda over a large geographical scale spanning the range of host genetic and spatial diversity. Using a multivariate analysis approach, we uncovered complex coinfection dynamics between the pathogens and statistically significant associations between host genetic groups and pathogen prevalence. It is important to note that these coinfection dynamics and associations with the host were not apparent by univariate analysis. These associations between host genotype and pathogen are particularly evident for Wolbachia and SGHV where host groups are inversely correlated for Wolbachia and SGHV prevalence. On the other hand, trypanosome infection prevalence is more complex and covaries with the presence of the other two pathogens, highlighting the importance of examining multiple pathogens simultaneously before making generalizations about infection and spatial patterns. It is imperative to note that these novel findings would have been missed if we had employed the standard univariate analysis used in previous studies. Our results are discussed in the context of disease epidemiology and vector control.

Influence of host phylogeographic patterns and incomplete lineage sorting on within-species genetic variability in Wigglesworthia species, obligate symbionts of tsetse flies.

Vertical transmission of obligate symbionts generates a predictable evolutionary history of symbionts that reflects that of their hosts. In insects, evolutionary associations between symbionts and their hosts have been investigated primarily among species, leaving population-level processes largely unknown. In this study, we investigated the tsetse (Diptera: Glossinidae) bacterial symbiont, Wigglesworthia glossinidia, to determine whether observed codiversification of symbiont and tsetse host species extends to a single host species (Glossina fuscipes fuscipes) in Uganda. To explore symbiont genetic variation in G. f. fuscipes populations, we screened two variable loci (lon and lepA) from the Wigglesworthia glossinidia bacterium in the host species Glossina fuscipes fuscipes (W. g. fuscipes) and examined phylogeographic and demographic characteristics in multiple host populations. Symbiont genetic variation was apparent within and among populations. We identified two distinct symbiont lineages, in northern and southern Uganda. Incongruence length difference (ILD) tests indicated that the two lineages corresponded exactly to northern and southern G. f. fuscipes mitochondrial DNA (mtDNA) haplogroups (P = 1.0). Analysis of molecular variance (AMOVA) confirmed that most variation was partitioned between the northern and southern lineages defined by host mtDNA (85.44%). However, ILD tests rejected finer-scale congruence within the northern and southern populations (P = 0.009). This incongruence was potentially due to incomplete lineage sorting that resulted in novel combinations of symbiont genetic variants and host background. Identifying these novel combinations may have public health significance, since tsetse is the sole vector of sleeping sickness and Wigglesworthia is known to influence host vector competence. Thus, understanding the adaptive value of these host-symbiont combinations may afford opportunities to develop vector control methods.

Trypanosoma brucei gambiense group 1 is distinguished by a unique amino acid substitution in the HpHb receptor implicated in human serum resistance.

Trypanosoma brucei rhodesiense (Tbr) and T. b. gambiense (Tbg), causative agents of Human African Trypanosomiasis (sleeping sickness) in Africa, have evolved alternative mechanisms of resisting the activity of trypanosome lytic factors (TLFs), components of innate immunity in human serum that protect against infection by other African trypanosomes. In Tbr, lytic activity is suppressed by the Tbr-specific serum-resistance associated (SRA) protein. The mechanism in Tbg is less well understood but has been hypothesized to involve altered activity and expression of haptoglobin haemoglobin receptor (HpHbR). HpHbR has been shown to facilitate internalization of TLF-1 in T.b. brucei (Tbb), a member of the T. brucei species complex that is susceptible to human serum. By evaluating the genetic variability of HpHbR in a comprehensive geographical and taxonomic context, we show that a single substitution that replaces leucine with serine at position 210 is conserved in the most widespread form of Tbg (Tbg group 1) and not found in related taxa, which are either human serum susceptible (Tbb) or known to resist lysis via an alternative mechanism (Tbr and Tbg group 2). We hypothesize that this single substitution contributes to reduced uptake of TLF and thus may play a key role in conferring serum resistance to Tbg group 1. In contrast, similarity in HpHbR sequence among isolates of Tbg group 2 and Tbb/Tbr provides further evidence that human serum resistance in Tbg group 2 is likely independent of HpHbR function.

Insight into the transmission biology and species-specific functional capabilities of tsetse (Diptera: glossinidae) obligate symbiont Wigglesworthia.

UNLABELLED: Ancient endosymbionts have been associated with extreme genome structural stability with little differentiation in gene inventory between sister species. Tsetse flies (Diptera: Glossinidae) harbor an obligate endosymbiont, Wigglesworthia, which has coevolved with the Glossina radiation. We report on the ~720-kb Wigglesworthia genome and its associated plasmid from Glossina morsitans morsitans and compare them to those of the symbiont from Glossina brevipalpis. While there was overall high synteny between the two genomes, a large inversion was noted. Furthermore, symbiont transcriptional analyses demonstrated host tissue and development-specific gene expression supporting robust transcriptional regulation in Wigglesworthia, an unprecedented observation in other obligate mutualist endosymbionts. Expression and immunohistochemistry confirmed the role of flagella during the vertical transmission process from mother to intrauterine progeny. The expression of nutrient provisioning genes (thiC and hemH) suggests that Wigglesworthia may function in dietary supplementation tailored toward host development. Furthermore, despite extensive conservation, unique genes were identified within both symbiont genomes that may result in distinct metabolomes impacting host physiology. One of these differences involves the chorismate, phenylalanine, and folate biosynthetic pathways, which are uniquely present in Wigglesworthia morsitans. Interestingly, African trypanosomes are auxotrophs for phenylalanine and folate and salvage both exogenously. It is possible that W. morsitans contributes to the higher parasite susceptibility of its host species. IMPORTANCE: Genomic stasis has historically been associated with obligate endosymbionts and their sister species. Here we characterize the Wigglesworthia genome of the tsetse fly species Glossina morsitans and compare it to its sister genome within G. brevipalpis. The similarity and variation between the genomes enabled specific hypotheses regarding functional biology. Expression analyses indicate significant levels of transcriptional regulation and support development- and tissue-specific functional roles for the symbiosis previously not observed in obligate mutualist symbionts. Retention of the genetically expensive flagella within these small genomes was demonstrated to be significant in symbiont transmission and tailored to the unique tsetse fly reproductive biology. Distinctions in metabolomes were also observed. We speculate an additional role for Wigglesworthia symbiosis where infections with pathogenic trypanosomes may depend upon symbiont species-specific metabolic products and thus influence the vector competence traits of different tsetse fly host species.