Antimicrobial resistance in Bacillus-based biopesticide products.

The crisis of antimicrobial resistant bacterial infections is one of the most pressing public health issues. Common agricultural practices have been implicated in the generation of antimicrobial resistant bacteria. Biopesticides, live bacteria used for pest control, are non-pathogenic and considered safe for consumption. Application of bacteria-based pesticides to crops in high concentrations raises the possibility of unintentional contributions to the movement and generation of antimicrobial resistance genes in the environment. However, the presence of clinically relevant antimicrobial resistance genes and their resistance phenotypes are currently unknown. Here we use a combination of multiple bioinformatic and microbiological techniques to define resistomes of widely used biopesticides and determine how the presence of suspected antimicrobial resistance genes translates to observable resistance phenotypes in several biopesticide products. Our results demonstrate that biopesticide products are reservoirs of clinically relevant antimicrobial resistance genes and bear resistance to multiple drug classes.

Taxonomic revision of the Australian arid zone lizards Gehyra variegata and G. montium (Squamata, Gekkonidae) with description of three new species.

The taxonomy of central Australian populations of geckos of the genus Gehyra has been uncertain since chromosomal studies carried out in the 1970s and 1980s revealed considerable heterogeneity and apparently independent patterns of morphological and karyotypic diversity. Following detailed molecular genetic studies, species boundaries in this complex have become clearer and we here re-set the boundaries of the three named species involved, G. variegata (Duméril & Bibron, 1836), G. montium Storr, 1982, and G. nana King, 1982, and describe three new species. Two of the new species, G. moritzi and G. pulingka, include populations formerly assigned to either G. montium or G. nana Storr, 1982, while the third, G. versicolor, includes all of the eastern Australian populations formerly assigned to G. variegata.

Comparative genomics of Japanese encephalitis virus shows low rates of recombination and a small subset of codon positions under episodic diversifying selection.

Orthoflavivirus japonicum (JEV) is the dominant cause of viral encephalitis in the Asian region with 100,000 cases and 25,000 deaths reported annually. The genome is comprised of a single polyprotein that encodes three structural and seven non-structural proteins. We collated a dataset of 349 complete genomes from a number of public databases, and analysed the data for recombination, evolutionary selection and phylogenetic structure. There are low rates of recombination in JEV, subsequently recombination is not a major evolutionary force shaping JEV. We found a strong overall signal of purifying selection in the genome, which is the main force affecting the evolutionary dynamics in JEV. There are also a small number of genomic sites under episodic diversifying selection, especially in the envelope protein and non-structural proteins 3 and 5. Overall, these results support previous analyses of JEV evolutionary genomics and provide additional insight into the evolutionary processes shaping the distribution and adaptation of this important pathogenic arbovirus.

Anonymous nuclear loci in non-model organisms: making the most of high-throughput genome surveys.

MOTIVATION: When working with non-model organisms, few if any species-specific markers are available for phylogenetic, phylogeographic and population studies. Therefore, researchers often try to adapt markers developed in distantly related taxa, resulting in poor amplification and ascertainment bias in their target taxa. Markers can be developed de novo and anonymous nuclear loci (ANL) are proving to be a boon for researchers seeking large numbers of fast-evolving, independent loci. However, the development of ANL can be laboratory intensive and expensive. A workflow is described to identify suitable low-copy anonymous loci from high-throughput shotgun sequences, dramatically reducing the cost and time required to develop these markers and produce robust multilocus datasets. RESULTS: By successively removing repetitive and evolutionary conserved sequences from low coverage shotgun libraries, we were able to isolate thousands of potential ANL. Empirical testing of loci developed from two reptile taxa confirmed that our methodology yields markers with comparable amplification rates and nucleotide diversities to ANLs developed using other methodologies. Our approach capitalizes on next-generation sequencing technologies to enable the development of phylogenetic, phylogeographic and population markers for taxa lacking suitable genomic resources.

Delimiting species in recent radiations with low levels of morphological divergence: a case study in Australian Gehyra geckos.

Recent conceptual and methodological advances have increased the ability to apply multifaceted approaches to species delimitation, which is particularly useful in delimiting recently diversified species where single lines of evidence lead to incorrect species delimitation or assignment of individuals to species (e.g. cryptic, morphological species and paraphyletic, hybridizing species). Species of the Australian Gehyra gecko radiation have historically proven difficult to delimit due the group's uniform, almost continent-wide geographic distribution and conservative morphology, contrasting high chromosomal and genetic diversity. Using an integrated approach to species delimitation taking advantage of morphological, geographic distributional and multi-locus genetic data, we investigate the diversity within three Gehyra species from the Australian arid zone. Our results show that these species represent eight distinct phylogenetic lineages, which display different patterns of morphological distinction and reproductive isolation. Using a recently developed Bayesian species delimitation method, we also find different levels of support for putative species dependent on priors for population size and timing of diversification assumed. Our results show that the current taxonomy does not adequately account for the diversity of the group. Discrepancies between lines of evidence indicate that diversification of the group is recent and ongoing, thus posing challenges for both species concepts and delimitation.

Going through phages: a computational approach to revealing the role of prophage in Staphylococcus aureus.

Prophages have important roles in virulence, antibiotic resistance, and genome evolution in Staphylococcus aureus . Rapid growth in the number of sequenced S. aureus genomes allows for an investigation of prophage sequences at an unprecedented scale. We developed a novel computational pipeline for phage discovery and annotation. We combined PhiSpy, a phage discovery tool, with VGAS and PROKKA, genome annotation tools to detect and analyse prophage sequences in nearly 10 011 S . aureus genomes, discovering thousands of putative prophage sequences with genes encoding virulence factors and antibiotic resistance. To our knowledge, this is the first large-scale application of PhiSpy on a large-scale set of genomes (10 011 S . aureus ). Determining the presence of virulence and resistance encoding genes in prophage has implications for the potential transfer of these genes/functions to other bacteria via transduction and thus can provide insight into the evolution and spread of these genes/functions between bacterial strains. While the phage we have identified may be known, these phages were not necessarily known or characterized in S. aureus and the clustering and comparison we did for phage based on their gene content is novel. Moreover, the reporting of these genes with the S. aureus genomes is novel.

Morphological differentiation correlates with ecological but not with genetic divergence in a Gehyra gecko.

Body size affects life history, the ecological niche of an organism and its interactions with other organisms. Resultantly, marked differences in body size between related organisms are often an indication of a species boundary. This is particularly evident in the Gehyra variegata species complex of geckos, which displays differential body sizes between genetically divergent species, but high levels of intraspecific morphological conservatism. We report on a Gehyra population that displays extraordinary body size differentiation in comparison with other G. variegata species. We used morphological and environmental data to show this population is phenotypically and ecologically distinct from its parapatric congener Gehyra lazelli and that morphology and ecology are significantly correlated. Contrastingly, mtDNA analysis indicates paraphyly between the two groups, and allele frequencies at six microsatellite loci show no population structure concordant with morpho-/ecotype. These results suggest either ecological speciation or environmentally induced phenotypic polymorphism, in an otherwise morphologically conservative group.

Genomic Sequence Analysis of Methicillin- and Carbapenem-Resistant Bacteria Isolated from Raw Sewage.

Antibiotic resistance is one of the largest threats facing global health. Wastewater treatment plants are well-known hot spots for interaction between diverse bacteria, genetic exchange, and antibiotic resistance. Nonpathogenic bacteria theoretically act as reservoirs of antibiotic resistance subsequently transferring antibiotic resistance genes to pathogens, indicating that evolutionary processes occur outside clinical settings and may drive patterns of drug-resistant infections. We isolated and sequenced 100 bacterial strains from five wastewater treatment plants to analyze regional dynamics of antibiotic resistance in the California Central Valley. The results demonstrate the presence of a wide diversity of pathogenic and nonpathogenic bacteria, with an arithmetic mean of 5.1 resistance genes per isolate. Forty-three percent of resistance genes were located on plasmids, suggesting that large levels of gene transfer between bacteria that otherwise may not co-occur are facilitated by wastewater treatment. One of the strains detected was a Bacillus carrying pX01 and pX02 anthrax-like plasmids and multiple drug resistance genes. A correlation between resistance genes and taxonomy indicates that taxon-specific evolutionary studies may be useful in determining and predicting patterns of antibiotic resistance. Conversely, a lack of geographic correlation may indicate that landscape genetic studies to understand the spread of antibiotic resistance genes should be carried out at broader scales. This large data set provides insights into how pathogenic and nonpathogenic bacteria interact in wastewater environments and the resistance genes which may be horizontally transferred between them. This can help in determining the mechanisms leading to the increasing prevalence of drug-resistant infections observed in clinical settings. IMPORTANCE The reasons for the increasing prevalence of antibiotic-resistant infections are complex and associated with myriad clinical and environmental processes. Wastewater treatment plants operate as nexuses of bacterial interaction and are known hot spots for genetic exchange between bacteria, including antibiotic resistance genes. We isolated and sequenced 100 drug-resistant bacteria from five wastewater treatment plants in California's Central Valley, characterizing widespread gene sharing between pathogens and nonpathogens. We identified a novel, multiresistant Bacillus carrying anthrax-like plasmids. This empirical study supports the likelihood of evolutionary and population processes in the broader environment affecting the prevalence of clinical drug-resistant infections and identifies several taxa that may operate as reservoirs and vectors of antibiotic resistance genes.

A Curated, Comprehensive Database of Plasmid Sequences.

Plasmid sequences are central to a myriad of microbial functions and processes. Here, we have compiled a database of complete plasmid sequences and associated metadata curated from both NCBI's recent genome database update, which includes plasmids as organisms, and all available annotated bacterial genomes. The resultant database contains 10,892 complete plasmid sequences and associated metadata.

Where the plasmids roam: large-scale sequence analysis reveals plasmids with large host ranges.

Describing the role of plasmids and their contribution to the exchange of genetic material among bacteria is essential for understanding the fields of plasmid epidemiology, microbial ecology, and commercial and synthetic microbiology. Broad-host-range (BHR) plasmids are those that are found not only in a single bacterial species, but in members of different taxonomic groups and are of significant interest to researchers in many fields. We applied a novel approach to computationally identify new BHR plasmids, in which we searched for highly similar cognate plasmids within a comprehensive plasmid database. After identifying 125 plasmid groups with highly similar cognates found in multiple taxa, we closely examined BHR plasmids found in multiple families. The majority of our identified BHR plasmids are found in members of the Enterobacteriaceae and closely related taxa, while three BHR plasmids of potential commercial significance were found in two species of Cyanobacteria. One plasmid with an exceptionally broad host range was found in both Gram-positive and Gram-negative bacterial species. This analysis demonstrates the utility of this method in identifying new BHR plasmids while highlighting unknown ranges of previously documented plasmids.

Quantifying the Evolutionary Conservation of Genes Encoding Multidrug Efflux Pumps in the ESKAPE Pathogens To Identify Antimicrobial Drug Targets.

Increasing rates of antibiotic-resistant bacterial infection are one of the most pressing contemporary global health concerns. The ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) have been identified as the leading global cause of multidrug-resistant bacterial infections, and overexpression of multidrug efflux (MEX) transport systems has been identified as one of the most critical mechanisms facilitating the evolution of multidrug resistance in ESKAPE pathogens. Despite efforts to develop efflux pump inhibitors to combat antibiotic resistance, the need persists to identify additional targets for future investigations. We evaluated evolutionary pressures on 110 MEX-encoding genes from all annotated ESKAPE organism genomes. We identify several MEX genes under stabilizing selection-representing targets which can facilitate broad-spectrum treatments with evolutionary constraints limiting the potential emergence of escape mutants. We also examine MEX systems being evaluated as drug targets, demonstrating that divergent selection may underlie some of the problems encountered in the development of effective treatments-specifically in relation to the NorA system in S. aureus. This study provides a comprehensive evolutionary context to efflux in the ESKAPE pathogens, which will provide critical context to the evaluation of efflux systems as antibiotic targets. IMPORTANCE Increasing rates of antibiotic-resistant bacterial infection are one of the most pressing contemporary global health concerns. The ESKAPE pathogen group represents the leading cause of these infections, and upregulation of efflux pump expression is a significant mechanism of resistance in these pathogens. This has resulted in substantial interest in the development of efflux pump inhibitors to combat antibiotic-resistant infections; however, no widespread treatments have been developed to date. Our study evaluates an often-underappreciated aspect of resistance-the impact of evolutionary selection. We evaluate selection on all annotated efflux genes in all sequenced ESKAPE pathogens, providing critical context for and insight into current and future development of efflux-targeting treatments for resistant bacterial infections.

Genomic analyses of African Trypanozoon strains to assess evolutionary relationships and identify markers for strain identification.

African trypanosomes of the sub-genus Trypanozoon) are eukaryotic parasitesthat cause disease in either humans or livestock. The development of genomic resources can be of great use to those interested in studying and controlling the spread of these trypanosomes. Here we present a large comparative analysis of Trypanozoon whole genomes, 83 in total, including human and animal infective African trypanosomes: 21 T. brucei brucei, 22 T. b. gambiense, 35 T. b. rhodesiense and 4 T. evansi strains, of which 21 were from Uganda. We constructed a maximum likelihood phylogeny based on 162,210 single nucleotide polymorphisms (SNPs.) The three Trypanosoma brucei sub-species and Trypanosoma evansi are not monophyletic, confirming earlier studies that indicated high similarity among Trypanosoma "sub-species". We also used discriminant analysis of principal components (DAPC) on the same set of SNPs, identifying seven genetic clusters. These clusters do not correspond well with existing taxonomic classifications, in agreement with the phylogenetic analysis. Geographic origin is reflected in both the phylogeny and clustering analysis. Finally, we used sparse linear discriminant analysis to rank SNPs by their informativeness in differentiating the strains in our data set. As few as 84 SNPs can completely distinguish the strains used in our study, and discriminant analysis was still able to detect genetic structure using as few as 10 SNPs. Our results reinforce earlier results of high genetic similarity between the African Trypanozoon. Despite this, a small subset of SNPs can be used to identify genetic markers that can be used for strain identification or other epidemiological investigations.

Genetic diversity of Glossina fuscipes fuscipes along the shores of Lake Victoria in Tanzania and Kenya: implications for management.

BACKGROUND: Tsetse flies (Diptera: Glossinidae) are sole vectors for trypanosomiasis, which affect human health and livestock productivity in Africa. Little is known about the genetic diversity of Glossina fuscipes fuscipes, which is an important species in Tanzania and Kenya. The main objective of the study was to provide baseline data to determine the genetic variability and divergence of G. f. fuscipes in the Lake Victoria basin of Tanzania and Kenya in order to guide future vector control efforts in the region. FINDINGS: Two hundred and seventy five G. f. fuscipes from 8 sites along the shores of Lake Victoria were screened for genetic polymorphisms at 19 microsatellite loci. Samples were collected from two sites in Kenya and six sites in Tanzania. Four of the Tanzanian sites were located in the Rorya district, on the eastern shores of Lake Victoria, while the other two sites were from Ukerewe and Bukoba districts from the southern and western Lake Victoria shores, respectively. Four genetically distinct allopatric clusters were revealed by microsatellite analysis, which sorted the sampling sites according to geography, with sites separated by as little as ~65 km belonging to distinct genetic clusters, while samples located within ~35 km from each other group in the same cluster. CONCLUSION: Our results suggest that there is ongoing genetic admixture within sampling sites located ~35 km from each other, while sites located ~65 km apart are genetically isolated from each other. Similar patterns emerged from a parallel study on G. f. fuscipes analyzed from the Lake Victoria Uganda shores. From a control perspective these results suggest that for sites within the same genetic cluster, control efforts should be carried out in a coordinated fashion in order to avoid re-invasions. Future work should focus on better quantifying the extent and spatial patterns of the observed genetic discontinuities of the G. f. fuscipes populations along the Tanzanian shores. This will aid in their control by providing guidelines on the geographical extent of the area to be treated at the same time.

De Novo Genome Assembly Shows Genome Wide Similarity between Trypanosoma brucei brucei and Trypanosoma brucei rhodesiense.

BACKGROUND: Trypanosoma brucei is a eukaryotic pathogen which causes African trypanosomiasis. It is notable for its variant surface glycoprotein (VSG) coat, which undergoes antigenic variation enabled by a large suite of VSG pseudogenes, allowing for persistent evasion of host adaptive immunity. While Trypanosoma brucei rhodesiense (Tbr) and T. b gambiense (Tbg) are human infective, related T. b. brucei (Tbb) is cleared by human sera. A single gene, the Serum Resistance Associated (SRA) gene, confers Tbr its human infectivity phenotype. Potential genetic recombination of this gene between Tbr and non-human infective Tbb strains has significant epidemiological consequences for Human African Trypanosomiasis outbreaks. RESULTS: Using long and short read whole genome sequencing, we generated a hybrid de novo assembly of a Tbr strain, producing 4,210 scaffolds totaling approximately 38.8 megabases, which comprise a significant proportion of the Tbr genome, and thus represents a valuable tool for a comparative genomics analyses among human and non-human infective T. brucei and future complete genome assembly. We detected 5,970 putative genes, of which two, an alcohol oxidoreductase and a pentatricopeptide repeat-containing protein, were members of gene families common to all T. brucei subspecies, but variants specific to the Tbr strain sequenced in this study. Our findings confirmed the extremely high level of genomic similarity between the two parasite subspecies found in other studies. CONCLUSIONS: We confirm at the whole genome level high similarity between the two Tbb and Tbr strains studied. The discovery of extremely minor genomic differentiation between Tbb and Tbr suggests that the transference of the SRA gene via genetic recombination could potentially result in novel human infective strains, thus all genetic backgrounds of T. brucei should be considered potentially human infective in regions where Tbr is prevalent.

Phage selection restores antibiotic sensitivity in MDR Pseudomonas aeruginosa.

Increasing prevalence and severity of multi-drug-resistant (MDR) bacterial infections has necessitated novel antibacterial strategies. Ideally, new approaches would target bacterial pathogens while exerting selection for reduced pathogenesis when these bacteria inevitably evolve resistance to therapeutic intervention. As an example of such a management strategy, we isolated a lytic bacteriophage, OMKO1, (family Myoviridae) of Pseudomonas aeruginosa that utilizes the outer membrane porin M (OprM) of the multidrug efflux systems MexAB and MexXY as a receptor-binding site. Results show that phage selection produces an evolutionary trade-off in MDR P. aeruginosa, whereby the evolution of bacterial resistance to phage attack changes the efflux pump mechanism, causing increased sensitivity to drugs from several antibiotic classes. Although modern phage therapy is still in its infancy, we conclude that phages, such as OMKO1, represent a new approach to phage therapy where bacteriophages exert selection for MDR bacteria to become increasingly sensitive to traditional antibiotics. This approach, using phages as targeted antibacterials, could extend the lifetime of our current antibiotics and potentially reduce the incidence of antibiotic resistant infections.

Evolution in spatially mixed host environments increases divergence for evolved fitness and intrapopulation genetic diversity in RNA viruses.

Virus populations may be challenged to evolve in spatially heterogeneous environments, such as mixtures of host cells that pose differing selection pressures. Spatial heterogeneity may select for evolved polymorphisms, where multiple virus subpopulations coexist by specializing on a narrow subset of the available hosts. Alternatively, spatial heterogeneity may select for evolved generalism, where a single genotype dominates the virus population by occupying a relatively broader host niche. In addition, the extent of spatial heterogeneity should influence the degree of divergence among virus populations encountering identical environmental challenges. Spatial heterogeneity creates environmental complexity that should increase the probability of differing adaptive phenotypic solutions, thus producing greater divergence among replicate virus populations, relative to counterparts evolving in strictly homogeneous host environments. Here, we tested these ideas using experimental evolution of RNA virus populations grown in laboratory tissue culture. We allowed vesicular stomatitis virus (VSV) lineages to evolve in replicated environments containing BHK-21 (baby hamster kidney) cells, HeLa (human epithelial) cells, or spatially heterogeneous host cell mixtures. Results showed that generalist phenotypes dominated in evolved virus populations across all treatments. Also, we observed greater variance in host-use performance (fitness) among VSV lineages evolved under spatial heterogeneity, relative to lineages evolved in homogeneous environments. Despite measurable differences in fitness, consensus Sanger sequencing revealed no fixed genetic differences separating the evolved lineages from their common ancestor. In contrast, deep sequencing of evolved VSV populations confirmed that the degree of divergence among replicate lineages was correlated with a larger number of minority variants. This correlation between divergence and the number of minority variants was significant only when we considered variants with a frequency of at least 10 per cent in the population. The number of lower-frequency minority variants per population did not significantly correlate with divergence.

Whole genome sequencing shows sleeping sickness relapse is due to parasite regrowth and not reinfection.

The trypanosome Trypanosoma brucei gambiense (Tbg) is a cause of human African trypanosomiasis (HAT) endemic to many parts of sub-Saharan Africa. The disease is almost invariably fatal if untreated and there is no vaccine, which makes monitoring and managing drug resistance highly relevant. A recent study of HAT cases from the Democratic Republic of the Congo reported a high incidence of relapses in patients treated with melarsoprol. Of the 19 Tbg strains isolated from patients enrolled in this study, four pairs were obtained from the same patient before treatment and after relapse. We used whole genome sequencing to investigate whether these patients were infected with a new strain, or if the original strain had regrown to pathogenic levels. Clustering analysis of 5938 single nucleotide polymorphisms supports the hypothesis of regrowth of the original strain, as we found that strains isolated before and after treatment from the same patient were more similar to each other than to other isolates. We also identified 23 novel genes that could affect melarsoprol sensitivity, representing a promising new set of targets for future functional studies. This work exemplifies the utility of using evolutionary approaches to provide novel insights and tools for disease control.

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola.

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph--and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration--specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.

Genetic diversity and population structure of Trypanosoma brucei in Uganda: implications for the epidemiology of sleeping sickness and Nagana.

BACKGROUND: While Human African Trypanosomiasis (HAT) is in decline on the continent of Africa, the disease still remains a major health problem in Uganda. There are recurrent sporadic outbreaks in the traditionally endemic areas in south-east Uganda, and continued spread to new unaffected areas in central Uganda. We evaluated the evolutionary dynamics underpinning the origin of new foci and the impact of host species on parasite genetic diversity in Uganda. We genotyped 269 Trypanosoma brucei isolates collected from different regions in Uganda and southwestern Kenya at 17 microsatellite loci, and checked for the presence of the SRA gene that confers human infectivity to T. b. rhodesiense. RESULTS: Both Bayesian clustering methods and Discriminant Analysis of Principal Components partition Trypanosoma brucei isolates obtained from Uganda and southwestern Kenya into three distinct genetic clusters. Clusters 1 and 3 include isolates from central and southern Uganda, while cluster 2 contains mostly isolates from southwestern Kenya. These three clusters are not sorted by subspecies designation (T. b. brucei vs T. b. rhodesiense), host or date of collection. The analyses also show evidence of genetic admixture among the three genetic clusters and long-range dispersal, suggesting recent and possibly on-going gene flow between them. CONCLUSIONS: Our results show that the expansion of the disease to the new foci in central Uganda occurred from the northward spread of T. b. rhodesiense (Tbr). They also confirm the emergence of the human infective strains (Tbr) from non-infective T. b. brucei (Tbb) strains of different genetic backgrounds, and the importance of cattle as Tbr reservoir, as confounders that shape the epidemiology of sleeping sickness in the region.

Comparative genomics reveals multiple genetic backgrounds of human pathogenicity in the Trypanosoma brucei complex.

The Trypanosoma brucei complex contains a number of subspecies with exceptionally variable life histories, including zoonotic subspecies, which are causative agents of human African trypanosomiasis (HAT) in sub-Saharan Africa. Paradoxically, genomic variation between taxa is extremely low. We analyzed the whole-genome sequences of 39 isolates across the T. brucei complex from diverse hosts and regions, identifying 608,501 single nucleotide polymorphisms that represent 2.33% of the nuclear genome. We show that human pathogenicity occurs across a wide range of parasite genotypes, and taxonomic designation does not reflect genetic variation across the group, as previous studies have suggested based on a small number of genes. This genome-wide study allowed the identification of significant host and geographic location associations. Strong purifying selection was detected in genomic regions associated with cytoskeleton structure, and regulatory genes associated with antigenic variation, suggesting conservation of these regions in African trypanosomes. In agreement with expectations drawn from meiotic reciprocal recombination, differences in average linkage disequilibrium between chromosomes in T. brucei correlate positively with chromosome size. In addition to insights into the life history of a diverse group of eukaryotic parasites, the documentation of genomic variation across the T. brucei complex and its association with specific hosts and geographic localities will aid in the development of comprehensive monitoring tools crucial to the proposed elimination of HAT by 2020, and on a shorter term, for monitoring the feared merger between the two human infective parasites, T. brucei rhodesiense and T. b. gambiense, in northern Uganda.