<title>Background</title> <p>Telomeres are preferred genomic locations for gene families involved in virulence and pathogenicity in a wide range of microbial pathogens including many unicellular eukaryotic parasites [##REF##16132817##1##,##REF##11495803##2##]. High rates of recombination at telomeres presumably facilitate the generation of genetic diversity [##REF##17319749##3##], as has been postulated for the malaria parasite <italic>Plasmodium</italic>, where antigenically variable <italic>VAR </italic>genes are located predominantly at chromosome ends [##REF##11069183##4##,##REF##15979386##5##]. The prokaryote <italic>Borrelia hermsii</italic>, causative agent of relapsing fever, also has large numbers of variable antigen genes located on linear plasmids which can undergo gene conversion events allowing antigenic variation [##REF##16796672##6##]. Lastly, African trypanosomes have segregated large gene families involved in antigenic variation or host adaptation to the telomeres of a broad range of chromosomes [##REF##16132817##1##,##REF##12547344##7##].</p> <p>The African trypanosome <italic>Trypanosoma brucei </italic>causes African Sleeping Sickness in humans, which is transmitted by tsetse flies and is endemic to sub-Saharan Africa [##REF##16735165##8##]. West African trypanosomiasis, comprising more than 90% of human cases, is the chronic form of the disease and is caused by <italic>T. b. gambiense </italic>[##REF##16650633##9##]. In contrast, East African trypanosomiasis is an acute infection in humans caused by <italic>T. b. rhodesiense</italic>, with death occurring typically within 6–8 months in the absence of treatment [##REF##16735165##8##]. <italic>T. b. rhodesiense </italic>is considered to be a zoonotic pathogen with a very extensive animal reservoir including a wide range of large game animals [##REF##4400764##10##], while <italic>T. b. gambiense </italic>infects a more restricted range of reservoir hosts [##REF##16236560##11##]. African trypanosomes also cause debilitating disease in livestock. <italic>T. b. brucei </italic>infects a variety of livestock and wildlife species, but is not human infective. Genetically, <italic>T. b. brucei </italic>is extremely similar to <italic>T. b. rhodesiense</italic>, but lacks the serum resistance associated gene <italic>SRA </italic>conferring human serum resistance [##REF##15099811##12##,##REF##17643434##13##]. Another non-human trypanosomiasis is dourine in equines, which is caused by <italic>T. equiperdum</italic>. This trypanosome, is also closely related to <italic>T. brucei </italic>[##REF##18245376##14##], but has lost the ability to cycle through tsetse flies and instead is sexually transmitted [##REF##9806490##15##]. <italic>T. equiperdum </italic>therefore has a much more restricted host range than other vector born <italic>Trypanosoma </italic>species.</p> <p>Many subspecies of trypanosomes can multiply extracellularly in the mammalian bloodstream, where they are exposed to continuous attack both by antibodies and complement. Key to survival is a dense protective coat of Variant Surface Glycoprotein (VSG) [##REF##645##16##], which is essential for the bloodstream form trypanosome even <italic>in vitro </italic>[##REF##15937117##17##]. As trypanosomes can switch between expression of different VSGs, new trypanosome antigenic types can escape host antibodies raised against the old VSG, and can avoid eradication by antibody-mediated lysis allowing a chronic infection to be maintained [##REF##11461029##18##, ####REF##15288623##19##, ##REF##16908087##20####16908087##20##]. Although a single trypanosome has up to 1500 <italic>VSG </italic>genes and pseudogenes [##REF##16020726##21##,##REF##16246028##22##], only one <italic>VSG </italic>is expressed at a time from one of about twenty bloodstream form <italic>VSG </italic>expression sites (BESs) [##REF##11742402##23##,##REF##11955440##24##].</p> <p>BESs are polycistronic transcription units located at telomeres, which contain a variety of Expression Site Associated Gene (<italic>ESAG</italic>) families in addition to the telomeric <italic>VSG </italic>[##REF##2720787##25##,##REF##12106867##26##]. It has been hypothesised that some <italic>ESAG</italic>s encode surface proteins or receptors which play a role in adaptation of the trypanosome to life in different species of mammalian host [##REF##9461219##27##,##REF##11356509##28##]. For example, <italic>ESAG6 </italic>and <italic>ESAG7 </italic>encode polymorphic subunits of a trypanosome encoded receptor for host transferrin, the variable nature of which could affect the ability of the trypanosome to take up transferrin molecules from different species of mammalian host [##REF##9461219##27##,##REF##8013456##29##]. In addition, it has been shown that the <italic>SRA </italic>gene located within some BESs confers resistance to human serum [##REF##9865701##30##]. As the trypanosome switches between different BESs, it expresses different permutations and combinations of polymorphic <italic>ESAG</italic>s [##REF##9461219##27##,##REF##11356509##28##,##REF##10607631##31##]. The telomeric location of <italic>VSG </italic>expression sites could play a role in their control [##REF##16132817##1##,##REF##12547344##7##]. Additionally, high rates of recombination at telomere ends could play a role in generating genetic diversity, as has been proposed for variant antigen genes located at the telomeres of the malaria parasite <italic>Plasmodium falciparum </italic>[##REF##11495803##2##,##REF##11069183##4##].</p> <p>Although a number of BESs from <italic>T. b. brucei </italic>and <italic>T. b. rhodesiense </italic>have been sequenced [##REF##2720787##25##,##REF##12106867##26##,##REF##9865701##30##,##REF##11353069##32##], relatively little is known about BESs in <italic>T. brucei gambiense</italic>, although they are clearly very similar to those from <italic>T. b. brucei </italic>[##REF##1762631##33##]. We are interested in the genetic diversity of telomeric BESs, particularly with regards to whether BES repertoires derived from trypanosomes with a large host range show evidence for greater genetic diversity compared with BES repertoires isolated from trypanosomes which infect a more restricted number of mammalian species. As telomeric regions are notoriously difficult to isolate and characterise, and are typically very underrepresented in standard libraries, we previously developed a method to isolate entire repertoires of BES containing telomeres in yeast [##REF##15520294##34##]. Here, we present the isolation and analysis of the genetic diversity found within repertoires of BES containing telomeres from three different species or subspecies of trypanosomes.</p>
<title>Methods</title> <title>Trypanosome strains and DNA isolation and analysis</title> <p>The genomic DNA used in the <italic>T. b. gambiense </italic>TAR BES telomere library produced in this study was isolated from the insect form <italic>Trypanosoma brucei gambiense </italic>DAL972 genome strain (received from Bill Wickstead, Keith Gull and Wendy Gibson) according to [##REF##6101223##53##]. DNA from bloodstream form <italic>T. b. brucei </italic>EATRO 2340 was a gift of Keith Matthews (University of Edinburgh) and was isolated according to [##REF##8341329##54##]. Note that in the literature this strain is referred to as <italic>T. b. rhodesiense </italic>EATRO 2340 [##REF##17043361##35##]. However, as we did not find evidence for the presence of <italic>SRA </italic>which is considered diagnostic for <italic>T. b. rhodesiense </italic>using <italic>SRA </italic>specific primers (Additional file ##SUPPL##0##1##) [##REF##17643434##13##,##REF##12518862##55##,##REF##12798017##56##], this strain was provisionally redesignated <italic>T. brucei brucei </italic>EATRO 2340. <italic>T. equiperdum </italic>STIB 818 (gift of Wendy Gibson, University of Bristol) [##REF##12664158##40##,##REF##1463353##57##,##REF##1311051##58##] was amplified in mice, and after purification of trypanosomes from whole blood using DEAE columns, genomic DNA was isolated according to [##REF##8341329##54##]. All animal experiments were subject to ethical review by the University of Oxford, and were conducted according to the conditions of a Home Office project licence.</p> <title>TAR cloning of <italic>Trypanosoma sp</italic>. BES telomeres</title> <p>Transformation Associated Recombination (TAR) cloning was performed essentially according to [##REF##15520294##34##] using the <italic>S. cerevisiae </italic>strain TYC1 as recipient (<italic>MAT</italic>α, <italic>ura3-52</italic>, <italic>leu2Δ1</italic>, <italic>cyh</italic>2<sup>r</sup>). The <italic>T. brucei </italic>BES-specific TAR vector pEB4 contains the yeast positive selectable marker <italic>URA3</italic>, the negative selectable marker <italic>CYH2</italic>, as well as a yeast centromere, origin of replication, and one yeast telomere. The <italic>T. brucei </italic>specific TAR target within the pEB4 vector is a 560 bp fragment with a region of the <italic>T. brucei </italic>BES core promoter that is particularly conserved over a range of different <italic>T. brucei </italic>strains [##REF##15520294##34##]. Recombination in yeast between the pEB4 TAR vector linearised with <italic>Cla</italic>I and the cotransfected <italic>T. brucei </italic>genomic DNA, results in the production of stable "half-YAC" containing a <italic>T. brucei </italic>BES as well as one yeast telomere and one trypanosome telomere.</p> <p>Yeast colonies were screened for the presence of <italic>T. brucei </italic>BES TAR clones with primers specific for <italic>ESAG6</italic>/<italic>ESAG7 </italic>corresponding to sequences located within particularly conserved regions of these genes. <italic>ESAG6 </italic>and <italic>ESAG7 </italic>are conserved elements within trypanosome BESs [##REF##12106867##26##], and their presence was considered diagnostic for the presence of a <italic>T. brucei </italic>BES. For initial <italic>T. b. gambiense </italic>TAR clone screening <italic>ESAG6/7 </italic>screening primers ESAG7a (sense and anti-sense) were used (sequences in [##REF##15520294##34##]). For screening the rest of the <italic>T. b. gambiense </italic>TAR clones and all of the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>TAR clones, screening primers ESAG6/7-311s and ESAG6-809as were used. See the table in Additional file ##SUPPL##8##9## for all primer sequences used.</p> <title>TAR clone analysis and sequencing</title> <p>In most cases, positive yeast transformants were streaked out for single colonies, these were expanded on quarter petri dishes, and yeast genomic DNA was isolated [##REF##15520294##34##]. PCR products were generated by amplification for 35 cycles at 94°C for 30 seconds, 50°C for 30 seconds, 72°C for 1 minute for products smaller than 1 kb, and 1 minute 30 seconds for products longer than 1 kb. Fragments were precipitated and sequenced using BigDye terminator (Applied Biosystems). Sequence analysis was performed with modifications of the procedure described in [##REF##15520294##34##]. Sequences of all primers used in this study are presented in the table in Additional file ##SUPPL##8##9##.</p> <p>TAR clones were typed into BES sets by sequencing the promoter region and an approximately 700 bp central stretch of <italic>ESAG6</italic>. Promoter sequence was trimmed to delete sequence corresponding to the region of the BES promoter present in the pEB4 TAR vector target fragment. Identification of sequence polymorphisms within both the BES promoter and <italic>ESAG6 </italic>allowed the categorisation of the different TAR clones into different sets corresponding to different BESs.</p> <p>For typing of the different TAR clones using ES promoter sequence, primers ESP S1/S2 (sense and anti-sense) (sequence in [##REF##15520294##34##]) were used for both the PCR and the sequencing step. The sequence was trimmed to delete the 5' region corresponding to the TAR vector target fragment. For <italic>ESAG6 </italic>sequence typing primers ESAG6-287s and ESAG6-1045as were used for both the PCR and the sequencing step.</p> <p>After initial typing of the TAR clones into different BES sets, the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>genes were amplified and sequenced from two different TAR clones from each BES set. The strategy followed was:</p> <title>ESAG6</title> <p>The 5' end was amplified and sequenced using ESAG6-UPSs and ESAG6-652as. The 3' end was amplified and sequenced using ESAG6-796s and ESAG6-DNas.</p> <title>ESAG5</title> <p>Initial PCR to amplify the coding region was performed with ESAG5-UPSs and ESAG5-DNas. If these failed, two separate reactions using ESAG5-240s and ESAG5-DNas and ESAG5-UPSs and ESAG5-578as were tried. The above products were sequenced with ESAG5-UPSs, ESAG5-DNas, ESAG5-240s, ESAG5-827s, ESAG5-578as and ESAG5-1007as.</p> <title>ESAG2</title> <p>Initial PCR was performed with ESAG2-UPSs and ESAG2-DN1as. If this failed, primers ESAG2-UPS2s and ESAG2-DN3as were used. The PCR product was cleaned and sequenced with ESAG2-UPS1s or ESAG2-UPS2s, ESAG2-DN1as or ESAG2-DN3as, ESAG2-371s, ESAG2-871s, ESAG2-505as, ESAG2-900as, ESAG2-970as or ESAG2-1088as.</p> <p>Particular <italic>ESAG</italic>s were deemed absent from a given TAR clone if the amplification strategy detailed above was unsuccessful and if all additional primer pairs failed to give rise to an amplified product. For <italic>ESAG5 </italic>the additional primer pairs were ESAG5-1 (sense and anti-sense), ESAG5-2 (sense and anti-sense) and ESAG5-3 (sense and anti-sense) (sequences in [##REF##15520294##34##]). For <italic>ESAG2 </italic>the additional primer pairs were ESAG2-1 (sense and anti-sense), ESAG2-2 (sense and anti-sense) and ESAG2-3 (sense and anti-sense) (sequences in [##REF##15520294##34##]) and ESAG2-457s and ESAG2-900as (sequences in Additional file ##SUPPL##8##9##). The sequences of all primers used for PCR or sequencing are indicated in the table in Additional file ##SUPPL##8##9##.</p> <p>The GenBank database accession numbers for sequences from <italic>Trypanosoma brucei gambiense </italic>DAL972 are: BES promoter types 1–11 [<ext-link ext-link-type="gen" xlink:href="EU726336">EU726336</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726346">EU726346</ext-link>], <italic>ESAG6 </italic>type 1 [<ext-link ext-link-type="gen" xlink:href="EU726347">EU726347</ext-link>] and <italic>ESAG6 </italic>types 3–11 [<ext-link ext-link-type="gen" xlink:href="EU726348">EU726348</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726356">EU726356</ext-link>], <italic>ESAG5 </italic>types 1–7 [<ext-link ext-link-type="gen" xlink:href="EU726357">EU726357</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726363">EU726363</ext-link>] and <italic>ESAG2 </italic>types 1–10 [<ext-link ext-link-type="gen" xlink:href="EU726364">EU726364</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726373">EU726373</ext-link>]. The GenBank database accession numbers for sequences from <italic>Trypanosoma brucei brucei </italic>EATRO 2340 are: BES promoter types 1–18 [<ext-link ext-link-type="gen" xlink:href="EU726409">EU726409</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726426">EU726426</ext-link>], <italic>ESAG6 </italic>types 1–22 [<ext-link ext-link-type="gen" xlink:href="EU726427">EU726427</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726448">EU726448</ext-link>], <italic>ESAG5 </italic>types 1–18 [<ext-link ext-link-type="gen" xlink:href="EU726449">EU726449</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726466">EU726466</ext-link>] and <italic>ESAG2 </italic>types 1–17 [<ext-link ext-link-type="gen" xlink:href="EU726467">EU726467</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726483">EU726483</ext-link>]. The GenBank database accession numbers for sequences from <italic>Trypanosoma equiperdum </italic>STIB 818 are: BES promoter types 1–11 [<ext-link ext-link-type="gen" xlink:href="EU726374">EU726374</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726384">EU726384</ext-link>], <italic>ESAG6 </italic>type 2 [<ext-link ext-link-type="gen" xlink:href="EU726385">EU726385</ext-link>], <italic>ESAG6 </italic>type 5 [<ext-link ext-link-type="gen" xlink:href="EU726386">EU726386</ext-link>], <italic>ESAG6 </italic>type 6 [<ext-link ext-link-type="gen" xlink:href="EU726387">EU726387</ext-link>], <italic>ESAG6 </italic>types 8–13 [<ext-link ext-link-type="gen" xlink:href="EU726388">EU726388</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726393">EU726393</ext-link>], <italic>ESAG5 </italic>types 1–11 [<ext-link ext-link-type="gen" xlink:href="EU726394">EU726394</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726404">EU726404</ext-link>] and <italic>ESAG2 </italic>types 1–4 [<ext-link ext-link-type="gen" xlink:href="EU726405">EU726405</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU726408">EU726408</ext-link>].</p> <title>Sequence analysis</title> <p>For the phylogenetic analyses sequence alignments were generated using ClustalX [##REF##9396791##59##], and then manually edited and gap-stripped. Ambiguously-aligned regions were removed. Phylogenetic analysis of each alignment was then performed with PhyML [##REF##14530136##60##], assuming a general time reversible (GTR) nucleotide substitution model with gamma-distributed rate variation among sites and a class of invariant sites. Rate parameters and initial trees were estimated from the data. Bootstrap values were estimated from 100 bootstrap samples. To summarize the overall diversity of these repertoires, we calculated the proportion of segregating sites (%S) and the mean pairwise diversity (π) of each of the 12 nucleotide alignments and 9 amino acid alignments.</p> <p>Each alignment was tested for evidence of recombination using the likelihood-based method implemented in the program (GARD [##REF##16818476##61##]; <ext-link ext-link-type="uri" xlink:href="http://www.datamonkey.org/GARD"/>). At each locus, the nucleotide substitution model used in the subsequent GARD analysis was chosen using the model selection tool provided with the HyPhy software package [##REF##15509596##62##]; <ext-link ext-link-type="uri" xlink:href="http://www.hyphy.org"/>, and rate variation was modelled with a discrete distribution with up to three distinct rate classes. Analyses were run until incorporation of additional breakpoints did not lead to a decrease in the corrected Akaike information criterion (c-AIC). When the initial GARD analysis failed to converge in the time allocated, the alignment was divided into two subregions at the breakpoint with the highest c-AIC-support in the initial analysis, and each subregion was then separately analysed with GARD.</p> <p>In order to determine the adaptive evolution of ESAG sequences, we used the method described in [##REF##16895925##63##] to characterize the relative rates of non-synonymous and synonymous substitutions within the ESAG6, ESAG5 and ESAG2 gene families from each of the three trypanosome subspecies. This method extends existing likelihood-based approaches for detecting adaptive evolution of protein-coding sequences to recombinant sequences by allowing the genealogy to vary across the alignment while sharing the parameters of the codon substitution model between tracts [##REF##12871927##64##,##REF##12872913##65##]. Each ESAG family was codon-aligned (with pseudogenes excluded and ambiguously-aligned regions stripped), and the breakpoints identified by GARD analysis were adjusted by up to two bases to coincide with codon boundaries. The neighbour joining trees inferred during GARD analysis were used as estimates of the genealogies of the putative recombination tracts.</p> <p>To assess the evidence for adaptive evolution, HyPhy was used to obtain maximum likelihood estimates of the branch lengths, the transition-to-transversion rate ratio (κ), and the parameters of two models of the distribution of the non-synonymous/synonymous rate ratio (ω). Under the nearly neutral model (M1a), ω is drawn from a discrete distribution with two classes, ω<sub>0 </sub>&lt; 1 and ω<sub>1 </sub>= 1, with relative weights p<sub>0 </sub>and p<sub>1</sub>. The selection model (M2a) allows for a third category, ω<sub>2 </sub>&gt; 1, of adaptively-evolving sites with relative weight p<sub>2</sub>. Because the neutral model is nested within the selection model (setting p2 = 0 in M2a recovers M1a), a likelihood ratio test for positive selection can be performed by comparing twice the log-likelihood difference between M1a and M2a with a χ<sup>2</sup>-distribution on two degrees of freedom [##REF##15514074##66##]. The relative nonsynonymous substitution rates at individual codons were estimated using the empirical Bayes method described in [##REF##15514074##66##]: model M2a, with the maximum likelihood estimated parameters, was taken as a prior distribution for ω, and Bayes' formula was used to calculate the posterior distribution of ω at each codon, given the sequence data. All calculations were carried out using a HyPhy batch language script written by Konrad Scheffler and included in the current distribution of the HyPhy software package.</p>
<title>Results and Discussion</title> <title>Isolation of BES containing telomere libraries from three trypanosomatids</title> <p>We cloned repertoires of BES containing telomeres from three trypanosome subspecies: <italic>T. b. gambiense </italic>DAL 972 (genome strain currently being sequenced by the Sanger Research Institute), <italic>T. b. brucei </italic>EATRO 2340 and <italic>T. equiperdum </italic>STIB 818. <italic>T. b. brucei </italic>EATRO 2340 was originally thought to be a <italic>T. b</italic>. <italic>rhodesiense </italic>subspecies [##REF##17043361##35##]. However, as we did not find evidence for the presence of <italic>SRA </italic>which is considered diagnostic for <italic>T. b. rhodesiense </italic>[##REF##15099811##12##,##REF##17643434##13##], this strain was tentatively redesignated to be a <italic>T. b. brucei </italic>(see Additional file ##SUPPL##0##1## and Materials). The BES containing telomeres were isolated in yeast using a method relying on Transformation Associated Recombination (TAR) cloning [##REF##15520294##34##,##REF##18323808##36##]. Linearised yeast TAR vector pEB4, containing a BES promoter fragment and a yeast telomere, was cotransformed into yeast spheroplasts together with total trypanosome genomic DNA (Fig. ##FIG##0##1##). Recombination between the BES promoter within the vector and similar sequences within the trypanosome genomic DNA provides the yeast vector with a second (trypanosome derived) telomere, thereby stabilising the episome as a YAC (yeast artificial chromosome).</p> <p>As BES sequences are polymorphic, we typed the different BES TAR clones into different BES sets based on the promoter and <italic>ESAG6 </italic>sequence types as previously described for <italic>T. b. brucei </italic>427 [##REF##15520294##34##]. For the BES promoter sequence typing, an approximately 635 bp region immediately downstream of the BES promoter (not overlapping with the target fragment used for the TAR cloning) was amplified from each TAR clone and sequenced (See Additional file ##SUPPL##1##2## for BES promoter sequence alignment). For the <italic>ESAG6 </italic>sequence typing an approximately 760 bp region spanning the <italic>ESAG6 </italic>hypervariable region was amplified and sequenced from each TAR clone (see Additional file ##SUPPL##2##3## for <italic>ESAG6 </italic>sequence alignments). The presence of sequence polymorphisms within these different BES regions allowed the typing of the different TAR clones into different BES sets.</p> <p>A total of 204 <italic>T. b. gambiense </italic>TAR clones were sorted into 13 putative BES sets (Table ##TAB##0##1##), 208 <italic>T. b. brucei </italic>clones were sorted into 23 BES sets (Table ##TAB##1##2##) and 91 <italic>T. equiperdum </italic>clones were sorted into 16 BES sets (Table ##TAB##2##3##). Two TAR clones were chosen for each BES type and the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>open reading frames were isolated and sequenced from both of these TAR clones (one representative clone is listed in each Table). These <italic>ESAG</italic>s were chosen for further analysis as they are similar enough between different BESs to allow PCR amplification and sequencing directly from the yeast TAR clones.</p> <p>We have previously isolated a repertoire of 17 BES types from <italic>T. b. brucei </italic>427 from a total of 182 TAR clones [##REF##15520294##34##]. This resulted in us cloning all ten BESs known to be functionally active <italic>in vitro </italic>in <italic>T. b. brucei </italic>427 in our laboratory [##REF##16135228##37##]. Given the number of TAR clones analysed here (204 for <italic>T. b. gambiense </italic>and 208 for <italic>T. b. brucei </italic>EATRO 2340) it is likely that we have isolated complete or nearly complete BES sets from these two trypanosome subspecies.</p> <title>BES repertoire size and diversity</title> <p>The overall nucleotide diversity of each repertoire, summarized in Table ##TAB##3##4##, shows both the percentage of nucleotides in different loci that are polymorphic (% S) as well as the mean pairwise diversity (π). Similar measures of amino acid diversity are summarized in Table ##TAB##4##5##. Two patterns are apparent: First, in general the nucleotide and amino acid diversity of each BES locus repertoire is greatest in <italic>T. b. brucei</italic>, intermediate in <italic>T. equiperdum</italic>, and least in <italic>T. b. gambiense</italic>. Since this is true of both measures of diversity, these observations cannot be explained by the presence of a few exceptionally divergent sequences or a few exceptionally polymorphic sites in the <italic>T. b. brucei </italic>repertoires. Rather, there are broad-based, sequence-wide differences in the diversity of these repertoires. Second, despite having a smaller number of sequence types identified, the <italic>ESAG2 </italic>repertoire in each subspecies is more polymorphic than the repertoires of either the BES promoter region, <italic>ESAG6 </italic>or <italic>ESAG5</italic>.</p> <p>Sequence repertoires of loci from <italic>T. b. gambiense </italic>were consistently the least genetically diverse. In <italic>T. b. brucei </italic>the transferrin receptor subunit ESAG6 has been shown to be particularly polymorphic over a hypervariable domain, as well as over a domain implicated in binding polymorphic host transferrins [##REF##1697265##38##, ####REF##9405356##39##, ##REF##12664158##40####12664158##40##]. It has been proposed that ESAG6 sequence polymorphisms allow the trypanosome to express receptors with differing affinities for the polymorphic transferrin molecules from different species of mammalian host [##REF##9461219##27##,##REF##11356509##28##]. Sequence alignment of these polymorphic regions of ESAG6 shows that ESAG6 sequences from <italic>T. b. brucei </italic>EATRO 2340 have comparable levels of amino acid diversity to that found in other <italic>T. brucei </italic>species [##REF##12664158##40##] (Fig. ##FIG##1##2##). However, as shown previously, ESAG6 sequences from <italic>T. equiperdum</italic>, while diverse in the ESAG6 hypervariable region, have little amino acid diversity within the transferrin binding region (Fig. ##FIG##1##2##)[##REF##12664158##40##,##REF##15589799##41##]. Here, we show that there is very little amino acid diversity in ESAG6 in <italic>T. b. gambiense </italic>DAL 972, even over these normally quite polymorphic stretches of the protein (Fig. ##FIG##1##2##). Surprisingly, six out of the seven <italic>T. b. gambiense </italic>ESAG6 protein types identified were identical over this polymorphic stretch.</p> <title>Phylogenetic analyses of BES sequence families</title> <p>Phylogenetic analyses (Fig. ##FIG##2##3## and Fig. ##FIG##3##4##) show that the relationships between the BES repertoires in the different trypanosome subspecies differ between loci. Within the BES promoter region, as well as at <italic>ESAG6 </italic>and <italic>ESAG5</italic>, most of the sequences obtained from <italic>T. b. gambiense </italic>belong to nearly monophyletic groups (<italic>gambiense</italic>-like) with little nucleotide or amino acid diversity. Nonetheless, at each locus, reciprocal monophyly between the <italic>T. b. gambiense </italic>repertoire and homologous sequences in the other two subspecies is violated, both by the existence of <italic>T. b. gambiense </italic>sequences which fall outside the main <italic>gambiense</italic>-like groups (e.g. Tbg_P1 for the promoter region; Tbg_E6-8 for <italic>ESAG6</italic>) and by the existence of <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>sequences which lie within these groups (e.g. Tbb_P17 for the promoter region and Teq_E5-8 for <italic>ESAG5</italic>).</p> <p>In contrast, the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>sequences are both collectively more diverse and are mutually paraphyletic. Not only are there several groups containing sequences from both subspecies, but there are several <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>homologs which are nearly identical (e.g. promoter region: Tbb_P-5 and Teq_P-10, <italic>ESAG6</italic>: Tbb_E6-5 and Teq_E6-11, and <italic>ESAG5</italic>: Tbb_E5-10 and Teq_E5-3) (Fig. ##FIG##2##3## and Fig. ##FIG##3##4##). Furthermore, to the extent that there is local, strain-based structure outside of the <italic>gambiense</italic>-like groups, this mainly consists of groups of <italic>T. b. brucei </italic>sequences which have no closely related sequences in the <italic>T. equiperdum </italic>genome (e.g., the group of <italic>ESAG6 </italic>sequences including Tbb_E6-11, Tbb_E6-10 and 10 other <italic>T. b. brucei </italic>sequences).</p> <p>The <italic>ESAG2 </italic>phylogenetic tree reveals a rather different set of relationships between the sequences from the different subspecies (Fig. ##FIG##3##4##). Apart from Tbg_E2-1 and Teq_E2-3, all of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum ESAG2 </italic>sequences belong to a single group which is itself only partially strain-structured. Indeed, one pair of <italic>T. b. gambiense </italic>and <italic>T. equiperdum ESAG2 </italic>homologs is nearly identical (Tbg_E2-6 and Teq_E2-1). Furthermore, most <italic>ESAG2 </italic>sequences can be divided into two relatively divergent subfamilies, one consisting almost exclusively of <italic>T. b. brucei </italic>sequences apart from Teq_E2-3 (<italic>ESAG2 </italic>type A sequences), and the other containing almost all of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum </italic>sequences as well as six <italic>T. b. brucei </italic>sequences (<italic>ESAG2 </italic>type B sequences). Inspection of the ESAG2 amino acid sequence alignments (see Additional file ##SUPPL##3##4##) reveals that the ESAG2 consensus type A and type B sequences differ at 123 of 474 residues (25.9%). This alignment also suggests that the most basally-branching ESAG2 sequences may be recombinants between type A and type B sequences. For example, although the Tbb_E2-12 ESAG2 sequence has a type A N-terminal domain (residues 1–54), the remaining residues predominantly agree with the type B consensus sequence. Likewise, the Tbg_E2-1 ESAG2 sequence is most similar to the type A consensus, but contains a tract of type B amino acids between residues 302–359.</p> <p>The ESAG5 amino acid sequences can also be partitioned into two subfamilies; although the differentiation is less pronounced, with consensus sequences which differ at 72 of 480 residues (15%) (see Additional file ##SUPPL##4##5##). The group of six relatively homogeneous <italic>T. equiperdum </italic>and <italic>T. b. brucei </italic>sequences in the <italic>ESAG5 </italic>phylogenetic tree constitute the <italic>ESAG5 </italic>type A subfamily, while the remaining sequences constitute the type B subfamily. As with <italic>ESAG2</italic>, the most basally-branching <italic>ESAG5 </italic>sequence appears to be a type A-type B recombinant. Although this Teq_E5-11 <italic>ESAG5 </italic>sequence is more closely related to the <italic>ESAG5 </italic>type A subfamily, it includes several tracts of type B consensus residues (e.g., residues 59–106).</p> <title>Evidence of recombination within BES loci</title> <p>Having observed putatively recombinant sequences in both the <italic>ESAG2 </italic>and <italic>ESAG5 </italic>families, we used a phylogenetic method to rigorously test each of the sequences from the four BES loci in the three trypanosome subspecies for evidence of recombination. The number of recombination breakpoints was inferred using GARD analysis (Table ##TAB##3##4##). With the exception of the <italic>T. b. gambiense ESAG5 </italic>alignment, the GARD analysis detected at least one breakpoint in each alignment, with the density of breakpoints increasing towards the telomere. For example, in the <italic>T. b. brucei </italic>sequence alignments, adjacent breakpoints are separated on average by 207 bp in the promoter region, 200 bp in <italic>ESAG6</italic>, 179 bp in <italic>ESAG5</italic>, and 129 bp in <italic>ESAG2</italic>. While comparable numbers of breakpoints are detected in the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>alignments, fewer are present in the corresponding <italic>T. b. gambiense </italic>alignments (Table ##TAB##3##4##).</p> <p>Although these analyses suggest that allelic recombination rates within the trypanosome BES repertoire might vary between loci and between subspecies, it is important to point out that both trends, that of increasing recombination towards the telomeres, and of reduced recombination within <italic>T. b. gambiense</italic>, are confounded by subspecies and locus-specific differences in genetic variation. Indeed, the mean number of sites between adjacent breakpoints inferred by GARD is significantly negatively correlated with both measures of nucleotide diversity reported in Table ##TAB##3##4## (tract length ~%S: r<sup>2 </sup>= 0.466, P = 0.001; tract length ~π: r<sup>2 </sup>= 0.462, P = 0.001). This correlation reflects the fact that the power of GARD to infer recombination breakpoints is limited by the amount of polymorphism in the data. In particular, it is possible that the homogeneity of the <italic>T. b. gambiense ESAG5 </italic>repertoire is due to extensive gene conversion.</p> <title>Adaptive evolution and amino acid diversity in ESAG repertoires</title> <p>In order to characterize the role of selection in the molecular evolution of the <italic>ESAG </italic>repertoires of <italic>T. b. brucei</italic>, <italic>T. b. gambiense</italic>, and <italic>T. equiperdum</italic>, we estimated the distribution of the relative ratio of non-synonymous-to-synonymous substitutions (dN/dS ratio or ω) in ESAG6, ESAG5 and ESAG2 from the three trypanosome subspecies. Maximum likelihood (ML) estimates of the parameters of the nearly neutral model M1a and the selection model M2a, as well as the p-values of the likelihood ratio tests comparing the selection model with the nearly neutral model are shown in the table in Additional file ##SUPPL##5##6##. These analyses provide evidence for adaptive evolution within the ESAG6, ESAG5, and ESAG2 repertoires of <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>, and within the ESAG2 repertoire of <italic>T. b. gambiense</italic>. In contrast, there is no statistically significant evidence for diversifying selection within the ESAG6 and ESAG5 repertoires of <italic>T. b. gambiense</italic>.</p> <p>Because inference of subspecies-specific selection is complicated by the lack of reciprocal monophyly between the homologous ESAG repertoires of the different strains, the contrasting patterns seen in the ESAG6 and ESAG5 repertoires between <italic>T. b. gambiense </italic>and the other two subspecies are particularly striking (Fig. ##FIG##2##3## and ##FIG##3##4##). On the one hand, the paraphyly of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum </italic>ESAG2 repertoires could explain why this locus does not follow the same pattern (e.g. the amino acid replacements in ESAG2 could predate the origin of <italic>T. b. gambiense</italic>). Alternatively, diversifying selection on the ESAG2 sequences in the <italic>T. b. gambiense </italic>genome may have prevented the extensive homogenisation that has occurred within the ESAG5 and ESAG6 families. The methods employed here cannot distinguish between these possibilities.</p> <p>In order to look in more detail at the patterns of molecular evolution within individual ESAGs, relative nonsynonymous substitution rates at individual codons were estimated using the prior distributions determined by the ML estimates of the model M2a parameters. The estimated rates are listed in the table in Additional file ##SUPPL##6##7## and plotted against residue number in Additional file ##SUPPL##7##8## (details in the table in Additional file ##SUPPL##6##7##). Consistent with earlier surveys of ESAG6 diversity [##REF##12664158##40##,##REF##15589799##41##], several residues in the hypervariable region of ESAG6 appear to be under diversifying selection in all three trypanosome subspecies. Adaptive evolution within the two transferrin-binding regions (box 1 and box 2) is evident in the <italic>T. b. brucei </italic>repertoire but not in <italic>T. b. gambiense </italic>or in <italic>T. equiperdum</italic>. The <italic>T. equiperdum </italic>data is consistent with data published in [##REF##12664158##40##,##REF##15589799##41##]. These functional elements of ESAG6 are annotated in the table in Additional file ##SUPPL##6##7##.</p> <p>One unexpected finding is that several residues within the signal peptide of each ESAG appear to be under diversifying selection, despite the fact that these regions are cleaved from the nascent protein [##REF##10940569##42##]. There are three such residues in ESAG6, five in ESAG5, and four in ESAG2. However, there are many other residues under diversifying selection in these ESAGs which have no known structural or functional association. Several of these residues are inferred to be evolving adaptively in all three subspecies, including residues 186 in ESAG6 and 81 and 434 in ESAG2.</p> <p>Table ##TAB##4##5## shows the numbers of amino acid residues within each alignment at which the non-synonymous substitution rate is more than twice as large as the synonymous substitution rate (ω &gt; 2). Also shown are the average estimates of ω at such residues, which is a measure of the strength of diversifying selection. As can be seen in Additional file ##SUPPL##7##8##, there are notable differences in the patterns of adaptive evolution in the three ESAG families. While the absolute number and density of diversifying residues in the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>repertoires is greatest in ESAG5 and least in ESAG2, the apparent strength of diversifying selection at adaptively evolving residues is greatest in ESAG2.</p> <p>Curiously, almost the opposite pattern is evident in the <italic>T. b. gambiense </italic>ESAG repertoires, with comparable numbers of diversifying codons in ESAG6 and ESAG2, but almost no amino acid variation in the ESAG5 sequences. It is also surprising that the smaller and less diverse <italic>T. equiperdum </italic>ESAG repertoires contain at least as many apparently adaptively-evolving sites as the corresponding <italic>T. b. brucei </italic>repertoires. This difference is most pronounced in ESAG6, where there are about 50% more diversifying residues in <italic>T. equiperdum </italic>than in <italic>T. b. brucei</italic>, despite the apparently stronger selection at such residues in the latter subspecies due to the larger size of its host range. Of course, as remarked above, all inter-strain comparisons are confounded by the paraphyletic relationships between these repertoires.</p> <title>BES sequence diversity and trypanosome biology</title> <p>The three trypanosome subspecies analysed here are highly similar at the DNA sequence level, yet cause diseases with very different pathologies. Whereas <italic>T. b. gambiense </italic>infection produces a chronic form of human disease, <italic>T. b. brucei </italic>is not human infective. Similarly, these trypanosome subspecies have different susceptibilities to drugs, whereby <italic>T. b. gambiense </italic>is very susceptible to the drug DFMO unlike <italic>T. b. brucei </italic>[##REF##9303385##43##]. <italic>T. equiperdum</italic>, although closely related to <italic>T. b. brucei </italic>[##REF##18245376##14##], infects different environments within the mammalian host. Rather than multiplying extracellularly in the mammalian bloodstream like <italic>T. b. brucei </italic>and <italic>T. b. gambiense</italic>, <italic>T. equiperdum </italic>is a tissue parasite primarily localised in the mucous membranes of the urogenital tract, and is rarely observed in the bloodstream [##REF##9806490##15##]. The challenge will come in trying to understand how the limited genetic differences between these subspecies translate into their very different pathologies.</p> <p>If ESAGs mediate host-parasite interactions, one could expect the diversity of individual BES repertoires to be positively correlated with host range size. This is because trypanosomes infecting multiple host species may be under selection to evolve and then maintain distinct ESAG alleles that are adapted to the different environments encountered in different species of mammalian host. In fact, our data are not completely consistent with this prediction. Of the three trypanosomes analysed here, host range is greatest in <italic>T. b. brucei </italic>[##REF##4400764##10##], intermediate in <italic>T. b. gambiense </italic>[##REF##16236560##11##], and smallest in <italic>T. equiperdum </italic>[##REF##9806490##15##]. However, although ESAG diversity is consistently higher in <italic>T. b. brucei </italic>than in the other two subspecies, the <italic>T. b. gambiense </italic>BES repertoires are actually less diverse than those of <italic>T. equiperdum</italic>. Likewise, although we find evidence for adaptive evolution in all three of the ESAG repertoires sequenced in <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>, only the ESAG2 repertoire of <italic>T. b. gambiense </italic>appears to be under diversifying selection. In addition, we find no correlation between the number of BESs found in a particular trypanosome subspecies and the host range size. Lastly, we find evidence for adaptive evolution of the ESAG6 repertoire of <italic>T. equiperdum</italic>, where it would be expected to be minimal or absent due to the restricted size of the host range of this trypanosome.</p> <p>There are several reasons why ESAG diversity and host range size might not be positively correlated. The hypothesis that diverse ESAG repertoires could facilitate trypanosome infection of diverse mammalian hosts was first proposed for ESAG6 [##REF##9461219##27##], where different sequence polymorphisms affect binding affinity of the transferrin receptor to variable transferrin molecules [##REF##9405356##39##]. This proposal has subsequently proved controversial. It has more recently been argued that <italic>in vivo </italic>concentrations of transferrin are high enough that expression of even low-affinity transferrin receptors is more than adequate to provide a trypanosome with sufficient transferrin in a range of different mammalian hosts [##REF##12643995##44##,##REF##15713451##45##]. These results would imply that ESAG6 sequence diversity within a given trypanosome subspecies need not be correlated with size of host range.</p> <p>Alternatively, even if host range is an important determinant of the selection pressures on ESAG repertoires, observed levels of ESAG diversity within individual genomes are probably influenced by a variety of processes that can create or remove variation. Sequence variation can be both generated and removed through DNA recombination. One possibility is that <italic>T. b. gambiense </italic>has unusually high rates of telomeric DNA rearrangements, resulting in sequence homogenisation of <italic>ESAG </italic>repertoires through allelic gene conversion or gene duplication and deletion. Another confounding factor is that genetic exchange between ESAGs and their genome-internal paralogs could influence the polymorphism observed within the BES repertoires. However, we can not easily reconcile why high rates of gene conversion in <italic>T. b. gambiense </italic>would lead to the relative homogenisation of the <italic>ESAG5 </italic>gene family, but not of the ESAG2 repertoire. We will only be able to understand the selection pressures maintaining <italic>ESAG </italic>sequence diversity within a given trypanosome subspecies once we know more about the function of these various <italic>ESAG </italic>gene families.</p> <p>Another factor that might influence the generation and maintenance of ESAG sequence diversity is the degree of genetic exchange that any given trypanosome subspecies undergoes. Different <italic>T. b. brucei </italic>and <italic>T. b. rhodesiense </italic>strains have been shown to undergo sex in the laboratory after passage through tsetse flies [##REF##7935602##46##, ####REF##7723777##47##, ##REF##15941603##48####15941603##48##], although the frequency with which this occurs in the field is unclear [##REF##8905086##49##]. This genetic exchange could allow the continuous flow of new ESAG sequences into a strain, facilitating the accumulation of genetic diversity. However in the Type I group of <italic>T. b. gambiense </italic>trypanosomes including <italic>T. b. gambiense </italic>DAL972 analysed here, there is no evidence for significant genetic exchange occurring [##REF##9088421##50##]. Population structure analyses indicate that <italic>T. b. gambiense </italic>Type I populations form a highly homogeneous group of trypanosomes which appear to have expanded clonally [##REF##10602679##51##,##REF##8538703##52##], and could have lost the ability to undergo meiosis [##REF##15941603##48##]. Lack of or extremely infrequent genetic exchange together with significant rates of telomeric gene conversion would be expected to facilitate the homogenisation of repetitive gene families in the absence of a strong selection pressure maintaining diversity</p> <p>Why then are the ESAG sequences from <italic>T. b. gambiense </italic>more homogeneous than those from <italic>T. equiperdum</italic>, which could have lost the opportunity to undergo meiosis altogether when it ceased to be transmitted via tsetse flies? The diversity within these ESAG repertoires is presumably not at equilibrium. <italic>T. equiperdum </italic>is likely to have arisen from <italic>T. b. brucei </italic>relatively recently (as proposed in [##REF##18245376##14##]), and thus may not have had time to completely lose ancestral ESAG sequence variation that is not actively being selected for in the equine host. In contrast, <italic>T. b. gambiense </italic>(Type I) is thought to be much more distantly related to <italic>T. b. brucei </italic>[##REF##8538703##52##]. This is consistent with the genealogical relationships of the ESAG6, ESAG5 and BES promoter repertoires in these three subspecies, i.e., the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>repertoires are paraphyletic, while the <italic>T. b. gambiense </italic>repertoires are nearly reciprocally monophyletic, but does not explain the very different relationships seen at ESAG2.</p> <p>A point to remember is that BESs are dynamic groups of genes. Each of the BES repertoires cloned presents a snapshot of what is going on within a population of trypanosomes, and each repertoire will contain a subset of ESAG 'alleles' present within the given trypanosome population. Further population structure analysis will allow us to determine if what we see in these analysed subspecies is a consequence of sampling rather than due to selection or consequences of loss of parts of the life-cycle. While these preliminary analyses are not the end of the story, these BES libraries will be a useful resource for further analyses of the processes influencing genetic variation in the telomeres of African trypanosomes.</p>
<title>Results and Discussion</title> <title>Isolation of BES containing telomere libraries from three trypanosomatids</title> <p>We cloned repertoires of BES containing telomeres from three trypanosome subspecies: <italic>T. b. gambiense </italic>DAL 972 (genome strain currently being sequenced by the Sanger Research Institute), <italic>T. b. brucei </italic>EATRO 2340 and <italic>T. equiperdum </italic>STIB 818. <italic>T. b. brucei </italic>EATRO 2340 was originally thought to be a <italic>T. b</italic>. <italic>rhodesiense </italic>subspecies [##REF##17043361##35##]. However, as we did not find evidence for the presence of <italic>SRA </italic>which is considered diagnostic for <italic>T. b. rhodesiense </italic>[##REF##15099811##12##,##REF##17643434##13##], this strain was tentatively redesignated to be a <italic>T. b. brucei </italic>(see Additional file ##SUPPL##0##1## and Materials). The BES containing telomeres were isolated in yeast using a method relying on Transformation Associated Recombination (TAR) cloning [##REF##15520294##34##,##REF##18323808##36##]. Linearised yeast TAR vector pEB4, containing a BES promoter fragment and a yeast telomere, was cotransformed into yeast spheroplasts together with total trypanosome genomic DNA (Fig. ##FIG##0##1##). Recombination between the BES promoter within the vector and similar sequences within the trypanosome genomic DNA provides the yeast vector with a second (trypanosome derived) telomere, thereby stabilising the episome as a YAC (yeast artificial chromosome).</p> <p>As BES sequences are polymorphic, we typed the different BES TAR clones into different BES sets based on the promoter and <italic>ESAG6 </italic>sequence types as previously described for <italic>T. b. brucei </italic>427 [##REF##15520294##34##]. For the BES promoter sequence typing, an approximately 635 bp region immediately downstream of the BES promoter (not overlapping with the target fragment used for the TAR cloning) was amplified from each TAR clone and sequenced (See Additional file ##SUPPL##1##2## for BES promoter sequence alignment). For the <italic>ESAG6 </italic>sequence typing an approximately 760 bp region spanning the <italic>ESAG6 </italic>hypervariable region was amplified and sequenced from each TAR clone (see Additional file ##SUPPL##2##3## for <italic>ESAG6 </italic>sequence alignments). The presence of sequence polymorphisms within these different BES regions allowed the typing of the different TAR clones into different BES sets.</p> <p>A total of 204 <italic>T. b. gambiense </italic>TAR clones were sorted into 13 putative BES sets (Table ##TAB##0##1##), 208 <italic>T. b. brucei </italic>clones were sorted into 23 BES sets (Table ##TAB##1##2##) and 91 <italic>T. equiperdum </italic>clones were sorted into 16 BES sets (Table ##TAB##2##3##). Two TAR clones were chosen for each BES type and the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>open reading frames were isolated and sequenced from both of these TAR clones (one representative clone is listed in each Table). These <italic>ESAG</italic>s were chosen for further analysis as they are similar enough between different BESs to allow PCR amplification and sequencing directly from the yeast TAR clones.</p> <p>We have previously isolated a repertoire of 17 BES types from <italic>T. b. brucei </italic>427 from a total of 182 TAR clones [##REF##15520294##34##]. This resulted in us cloning all ten BESs known to be functionally active <italic>in vitro </italic>in <italic>T. b. brucei </italic>427 in our laboratory [##REF##16135228##37##]. Given the number of TAR clones analysed here (204 for <italic>T. b. gambiense </italic>and 208 for <italic>T. b. brucei </italic>EATRO 2340) it is likely that we have isolated complete or nearly complete BES sets from these two trypanosome subspecies.</p> <title>BES repertoire size and diversity</title> <p>The overall nucleotide diversity of each repertoire, summarized in Table ##TAB##3##4##, shows both the percentage of nucleotides in different loci that are polymorphic (% S) as well as the mean pairwise diversity (π). Similar measures of amino acid diversity are summarized in Table ##TAB##4##5##. Two patterns are apparent: First, in general the nucleotide and amino acid diversity of each BES locus repertoire is greatest in <italic>T. b. brucei</italic>, intermediate in <italic>T. equiperdum</italic>, and least in <italic>T. b. gambiense</italic>. Since this is true of both measures of diversity, these observations cannot be explained by the presence of a few exceptionally divergent sequences or a few exceptionally polymorphic sites in the <italic>T. b. brucei </italic>repertoires. Rather, there are broad-based, sequence-wide differences in the diversity of these repertoires. Second, despite having a smaller number of sequence types identified, the <italic>ESAG2 </italic>repertoire in each subspecies is more polymorphic than the repertoires of either the BES promoter region, <italic>ESAG6 </italic>or <italic>ESAG5</italic>.</p> <p>Sequence repertoires of loci from <italic>T. b. gambiense </italic>were consistently the least genetically diverse. In <italic>T. b. brucei </italic>the transferrin receptor subunit ESAG6 has been shown to be particularly polymorphic over a hypervariable domain, as well as over a domain implicated in binding polymorphic host transferrins [##REF##1697265##38##, ####REF##9405356##39##, ##REF##12664158##40####12664158##40##]. It has been proposed that ESAG6 sequence polymorphisms allow the trypanosome to express receptors with differing affinities for the polymorphic transferrin molecules from different species of mammalian host [##REF##9461219##27##,##REF##11356509##28##]. Sequence alignment of these polymorphic regions of ESAG6 shows that ESAG6 sequences from <italic>T. b. brucei </italic>EATRO 2340 have comparable levels of amino acid diversity to that found in other <italic>T. brucei </italic>species [##REF##12664158##40##] (Fig. ##FIG##1##2##). However, as shown previously, ESAG6 sequences from <italic>T. equiperdum</italic>, while diverse in the ESAG6 hypervariable region, have little amino acid diversity within the transferrin binding region (Fig. ##FIG##1##2##)[##REF##12664158##40##,##REF##15589799##41##]. Here, we show that there is very little amino acid diversity in ESAG6 in <italic>T. b. gambiense </italic>DAL 972, even over these normally quite polymorphic stretches of the protein (Fig. ##FIG##1##2##). Surprisingly, six out of the seven <italic>T. b. gambiense </italic>ESAG6 protein types identified were identical over this polymorphic stretch.</p> <title>Phylogenetic analyses of BES sequence families</title> <p>Phylogenetic analyses (Fig. ##FIG##2##3## and Fig. ##FIG##3##4##) show that the relationships between the BES repertoires in the different trypanosome subspecies differ between loci. Within the BES promoter region, as well as at <italic>ESAG6 </italic>and <italic>ESAG5</italic>, most of the sequences obtained from <italic>T. b. gambiense </italic>belong to nearly monophyletic groups (<italic>gambiense</italic>-like) with little nucleotide or amino acid diversity. Nonetheless, at each locus, reciprocal monophyly between the <italic>T. b. gambiense </italic>repertoire and homologous sequences in the other two subspecies is violated, both by the existence of <italic>T. b. gambiense </italic>sequences which fall outside the main <italic>gambiense</italic>-like groups (e.g. Tbg_P1 for the promoter region; Tbg_E6-8 for <italic>ESAG6</italic>) and by the existence of <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>sequences which lie within these groups (e.g. Tbb_P17 for the promoter region and Teq_E5-8 for <italic>ESAG5</italic>).</p> <p>In contrast, the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>sequences are both collectively more diverse and are mutually paraphyletic. Not only are there several groups containing sequences from both subspecies, but there are several <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>homologs which are nearly identical (e.g. promoter region: Tbb_P-5 and Teq_P-10, <italic>ESAG6</italic>: Tbb_E6-5 and Teq_E6-11, and <italic>ESAG5</italic>: Tbb_E5-10 and Teq_E5-3) (Fig. ##FIG##2##3## and Fig. ##FIG##3##4##). Furthermore, to the extent that there is local, strain-based structure outside of the <italic>gambiense</italic>-like groups, this mainly consists of groups of <italic>T. b. brucei </italic>sequences which have no closely related sequences in the <italic>T. equiperdum </italic>genome (e.g., the group of <italic>ESAG6 </italic>sequences including Tbb_E6-11, Tbb_E6-10 and 10 other <italic>T. b. brucei </italic>sequences).</p> <p>The <italic>ESAG2 </italic>phylogenetic tree reveals a rather different set of relationships between the sequences from the different subspecies (Fig. ##FIG##3##4##). Apart from Tbg_E2-1 and Teq_E2-3, all of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum ESAG2 </italic>sequences belong to a single group which is itself only partially strain-structured. Indeed, one pair of <italic>T. b. gambiense </italic>and <italic>T. equiperdum ESAG2 </italic>homologs is nearly identical (Tbg_E2-6 and Teq_E2-1). Furthermore, most <italic>ESAG2 </italic>sequences can be divided into two relatively divergent subfamilies, one consisting almost exclusively of <italic>T. b. brucei </italic>sequences apart from Teq_E2-3 (<italic>ESAG2 </italic>type A sequences), and the other containing almost all of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum </italic>sequences as well as six <italic>T. b. brucei </italic>sequences (<italic>ESAG2 </italic>type B sequences). Inspection of the ESAG2 amino acid sequence alignments (see Additional file ##SUPPL##3##4##) reveals that the ESAG2 consensus type A and type B sequences differ at 123 of 474 residues (25.9%). This alignment also suggests that the most basally-branching ESAG2 sequences may be recombinants between type A and type B sequences. For example, although the Tbb_E2-12 ESAG2 sequence has a type A N-terminal domain (residues 1–54), the remaining residues predominantly agree with the type B consensus sequence. Likewise, the Tbg_E2-1 ESAG2 sequence is most similar to the type A consensus, but contains a tract of type B amino acids between residues 302–359.</p> <p>The ESAG5 amino acid sequences can also be partitioned into two subfamilies; although the differentiation is less pronounced, with consensus sequences which differ at 72 of 480 residues (15%) (see Additional file ##SUPPL##4##5##). The group of six relatively homogeneous <italic>T. equiperdum </italic>and <italic>T. b. brucei </italic>sequences in the <italic>ESAG5 </italic>phylogenetic tree constitute the <italic>ESAG5 </italic>type A subfamily, while the remaining sequences constitute the type B subfamily. As with <italic>ESAG2</italic>, the most basally-branching <italic>ESAG5 </italic>sequence appears to be a type A-type B recombinant. Although this Teq_E5-11 <italic>ESAG5 </italic>sequence is more closely related to the <italic>ESAG5 </italic>type A subfamily, it includes several tracts of type B consensus residues (e.g., residues 59–106).</p> <title>Evidence of recombination within BES loci</title> <p>Having observed putatively recombinant sequences in both the <italic>ESAG2 </italic>and <italic>ESAG5 </italic>families, we used a phylogenetic method to rigorously test each of the sequences from the four BES loci in the three trypanosome subspecies for evidence of recombination. The number of recombination breakpoints was inferred using GARD analysis (Table ##TAB##3##4##). With the exception of the <italic>T. b. gambiense ESAG5 </italic>alignment, the GARD analysis detected at least one breakpoint in each alignment, with the density of breakpoints increasing towards the telomere. For example, in the <italic>T. b. brucei </italic>sequence alignments, adjacent breakpoints are separated on average by 207 bp in the promoter region, 200 bp in <italic>ESAG6</italic>, 179 bp in <italic>ESAG5</italic>, and 129 bp in <italic>ESAG2</italic>. While comparable numbers of breakpoints are detected in the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>alignments, fewer are present in the corresponding <italic>T. b. gambiense </italic>alignments (Table ##TAB##3##4##).</p> <p>Although these analyses suggest that allelic recombination rates within the trypanosome BES repertoire might vary between loci and between subspecies, it is important to point out that both trends, that of increasing recombination towards the telomeres, and of reduced recombination within <italic>T. b. gambiense</italic>, are confounded by subspecies and locus-specific differences in genetic variation. Indeed, the mean number of sites between adjacent breakpoints inferred by GARD is significantly negatively correlated with both measures of nucleotide diversity reported in Table ##TAB##3##4## (tract length ~%S: r<sup>2 </sup>= 0.466, P = 0.001; tract length ~π: r<sup>2 </sup>= 0.462, P = 0.001). This correlation reflects the fact that the power of GARD to infer recombination breakpoints is limited by the amount of polymorphism in the data. In particular, it is possible that the homogeneity of the <italic>T. b. gambiense ESAG5 </italic>repertoire is due to extensive gene conversion.</p> <title>Adaptive evolution and amino acid diversity in ESAG repertoires</title> <p>In order to characterize the role of selection in the molecular evolution of the <italic>ESAG </italic>repertoires of <italic>T. b. brucei</italic>, <italic>T. b. gambiense</italic>, and <italic>T. equiperdum</italic>, we estimated the distribution of the relative ratio of non-synonymous-to-synonymous substitutions (dN/dS ratio or ω) in ESAG6, ESAG5 and ESAG2 from the three trypanosome subspecies. Maximum likelihood (ML) estimates of the parameters of the nearly neutral model M1a and the selection model M2a, as well as the p-values of the likelihood ratio tests comparing the selection model with the nearly neutral model are shown in the table in Additional file ##SUPPL##5##6##. These analyses provide evidence for adaptive evolution within the ESAG6, ESAG5, and ESAG2 repertoires of <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>, and within the ESAG2 repertoire of <italic>T. b. gambiense</italic>. In contrast, there is no statistically significant evidence for diversifying selection within the ESAG6 and ESAG5 repertoires of <italic>T. b. gambiense</italic>.</p> <p>Because inference of subspecies-specific selection is complicated by the lack of reciprocal monophyly between the homologous ESAG repertoires of the different strains, the contrasting patterns seen in the ESAG6 and ESAG5 repertoires between <italic>T. b. gambiense </italic>and the other two subspecies are particularly striking (Fig. ##FIG##2##3## and ##FIG##3##4##). On the one hand, the paraphyly of the <italic>T. b. gambiense </italic>and <italic>T. equiperdum </italic>ESAG2 repertoires could explain why this locus does not follow the same pattern (e.g. the amino acid replacements in ESAG2 could predate the origin of <italic>T. b. gambiense</italic>). Alternatively, diversifying selection on the ESAG2 sequences in the <italic>T. b. gambiense </italic>genome may have prevented the extensive homogenisation that has occurred within the ESAG5 and ESAG6 families. The methods employed here cannot distinguish between these possibilities.</p> <p>In order to look in more detail at the patterns of molecular evolution within individual ESAGs, relative nonsynonymous substitution rates at individual codons were estimated using the prior distributions determined by the ML estimates of the model M2a parameters. The estimated rates are listed in the table in Additional file ##SUPPL##6##7## and plotted against residue number in Additional file ##SUPPL##7##8## (details in the table in Additional file ##SUPPL##6##7##). Consistent with earlier surveys of ESAG6 diversity [##REF##12664158##40##,##REF##15589799##41##], several residues in the hypervariable region of ESAG6 appear to be under diversifying selection in all three trypanosome subspecies. Adaptive evolution within the two transferrin-binding regions (box 1 and box 2) is evident in the <italic>T. b. brucei </italic>repertoire but not in <italic>T. b. gambiense </italic>or in <italic>T. equiperdum</italic>. The <italic>T. equiperdum </italic>data is consistent with data published in [##REF##12664158##40##,##REF##15589799##41##]. These functional elements of ESAG6 are annotated in the table in Additional file ##SUPPL##6##7##.</p> <p>One unexpected finding is that several residues within the signal peptide of each ESAG appear to be under diversifying selection, despite the fact that these regions are cleaved from the nascent protein [##REF##10940569##42##]. There are three such residues in ESAG6, five in ESAG5, and four in ESAG2. However, there are many other residues under diversifying selection in these ESAGs which have no known structural or functional association. Several of these residues are inferred to be evolving adaptively in all three subspecies, including residues 186 in ESAG6 and 81 and 434 in ESAG2.</p> <p>Table ##TAB##4##5## shows the numbers of amino acid residues within each alignment at which the non-synonymous substitution rate is more than twice as large as the synonymous substitution rate (ω &gt; 2). Also shown are the average estimates of ω at such residues, which is a measure of the strength of diversifying selection. As can be seen in Additional file ##SUPPL##7##8##, there are notable differences in the patterns of adaptive evolution in the three ESAG families. While the absolute number and density of diversifying residues in the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>repertoires is greatest in ESAG5 and least in ESAG2, the apparent strength of diversifying selection at adaptively evolving residues is greatest in ESAG2.</p> <p>Curiously, almost the opposite pattern is evident in the <italic>T. b. gambiense </italic>ESAG repertoires, with comparable numbers of diversifying codons in ESAG6 and ESAG2, but almost no amino acid variation in the ESAG5 sequences. It is also surprising that the smaller and less diverse <italic>T. equiperdum </italic>ESAG repertoires contain at least as many apparently adaptively-evolving sites as the corresponding <italic>T. b. brucei </italic>repertoires. This difference is most pronounced in ESAG6, where there are about 50% more diversifying residues in <italic>T. equiperdum </italic>than in <italic>T. b. brucei</italic>, despite the apparently stronger selection at such residues in the latter subspecies due to the larger size of its host range. Of course, as remarked above, all inter-strain comparisons are confounded by the paraphyletic relationships between these repertoires.</p> <title>BES sequence diversity and trypanosome biology</title> <p>The three trypanosome subspecies analysed here are highly similar at the DNA sequence level, yet cause diseases with very different pathologies. Whereas <italic>T. b. gambiense </italic>infection produces a chronic form of human disease, <italic>T. b. brucei </italic>is not human infective. Similarly, these trypanosome subspecies have different susceptibilities to drugs, whereby <italic>T. b. gambiense </italic>is very susceptible to the drug DFMO unlike <italic>T. b. brucei </italic>[##REF##9303385##43##]. <italic>T. equiperdum</italic>, although closely related to <italic>T. b. brucei </italic>[##REF##18245376##14##], infects different environments within the mammalian host. Rather than multiplying extracellularly in the mammalian bloodstream like <italic>T. b. brucei </italic>and <italic>T. b. gambiense</italic>, <italic>T. equiperdum </italic>is a tissue parasite primarily localised in the mucous membranes of the urogenital tract, and is rarely observed in the bloodstream [##REF##9806490##15##]. The challenge will come in trying to understand how the limited genetic differences between these subspecies translate into their very different pathologies.</p> <p>If ESAGs mediate host-parasite interactions, one could expect the diversity of individual BES repertoires to be positively correlated with host range size. This is because trypanosomes infecting multiple host species may be under selection to evolve and then maintain distinct ESAG alleles that are adapted to the different environments encountered in different species of mammalian host. In fact, our data are not completely consistent with this prediction. Of the three trypanosomes analysed here, host range is greatest in <italic>T. b. brucei </italic>[##REF##4400764##10##], intermediate in <italic>T. b. gambiense </italic>[##REF##16236560##11##], and smallest in <italic>T. equiperdum </italic>[##REF##9806490##15##]. However, although ESAG diversity is consistently higher in <italic>T. b. brucei </italic>than in the other two subspecies, the <italic>T. b. gambiense </italic>BES repertoires are actually less diverse than those of <italic>T. equiperdum</italic>. Likewise, although we find evidence for adaptive evolution in all three of the ESAG repertoires sequenced in <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>, only the ESAG2 repertoire of <italic>T. b. gambiense </italic>appears to be under diversifying selection. In addition, we find no correlation between the number of BESs found in a particular trypanosome subspecies and the host range size. Lastly, we find evidence for adaptive evolution of the ESAG6 repertoire of <italic>T. equiperdum</italic>, where it would be expected to be minimal or absent due to the restricted size of the host range of this trypanosome.</p> <p>There are several reasons why ESAG diversity and host range size might not be positively correlated. The hypothesis that diverse ESAG repertoires could facilitate trypanosome infection of diverse mammalian hosts was first proposed for ESAG6 [##REF##9461219##27##], where different sequence polymorphisms affect binding affinity of the transferrin receptor to variable transferrin molecules [##REF##9405356##39##]. This proposal has subsequently proved controversial. It has more recently been argued that <italic>in vivo </italic>concentrations of transferrin are high enough that expression of even low-affinity transferrin receptors is more than adequate to provide a trypanosome with sufficient transferrin in a range of different mammalian hosts [##REF##12643995##44##,##REF##15713451##45##]. These results would imply that ESAG6 sequence diversity within a given trypanosome subspecies need not be correlated with size of host range.</p> <p>Alternatively, even if host range is an important determinant of the selection pressures on ESAG repertoires, observed levels of ESAG diversity within individual genomes are probably influenced by a variety of processes that can create or remove variation. Sequence variation can be both generated and removed through DNA recombination. One possibility is that <italic>T. b. gambiense </italic>has unusually high rates of telomeric DNA rearrangements, resulting in sequence homogenisation of <italic>ESAG </italic>repertoires through allelic gene conversion or gene duplication and deletion. Another confounding factor is that genetic exchange between ESAGs and their genome-internal paralogs could influence the polymorphism observed within the BES repertoires. However, we can not easily reconcile why high rates of gene conversion in <italic>T. b. gambiense </italic>would lead to the relative homogenisation of the <italic>ESAG5 </italic>gene family, but not of the ESAG2 repertoire. We will only be able to understand the selection pressures maintaining <italic>ESAG </italic>sequence diversity within a given trypanosome subspecies once we know more about the function of these various <italic>ESAG </italic>gene families.</p> <p>Another factor that might influence the generation and maintenance of ESAG sequence diversity is the degree of genetic exchange that any given trypanosome subspecies undergoes. Different <italic>T. b. brucei </italic>and <italic>T. b. rhodesiense </italic>strains have been shown to undergo sex in the laboratory after passage through tsetse flies [##REF##7935602##46##, ####REF##7723777##47##, ##REF##15941603##48####15941603##48##], although the frequency with which this occurs in the field is unclear [##REF##8905086##49##]. This genetic exchange could allow the continuous flow of new ESAG sequences into a strain, facilitating the accumulation of genetic diversity. However in the Type I group of <italic>T. b. gambiense </italic>trypanosomes including <italic>T. b. gambiense </italic>DAL972 analysed here, there is no evidence for significant genetic exchange occurring [##REF##9088421##50##]. Population structure analyses indicate that <italic>T. b. gambiense </italic>Type I populations form a highly homogeneous group of trypanosomes which appear to have expanded clonally [##REF##10602679##51##,##REF##8538703##52##], and could have lost the ability to undergo meiosis [##REF##15941603##48##]. Lack of or extremely infrequent genetic exchange together with significant rates of telomeric gene conversion would be expected to facilitate the homogenisation of repetitive gene families in the absence of a strong selection pressure maintaining diversity</p> <p>Why then are the ESAG sequences from <italic>T. b. gambiense </italic>more homogeneous than those from <italic>T. equiperdum</italic>, which could have lost the opportunity to undergo meiosis altogether when it ceased to be transmitted via tsetse flies? The diversity within these ESAG repertoires is presumably not at equilibrium. <italic>T. equiperdum </italic>is likely to have arisen from <italic>T. b. brucei </italic>relatively recently (as proposed in [##REF##18245376##14##]), and thus may not have had time to completely lose ancestral ESAG sequence variation that is not actively being selected for in the equine host. In contrast, <italic>T. b. gambiense </italic>(Type I) is thought to be much more distantly related to <italic>T. b. brucei </italic>[##REF##8538703##52##]. This is consistent with the genealogical relationships of the ESAG6, ESAG5 and BES promoter repertoires in these three subspecies, i.e., the <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>repertoires are paraphyletic, while the <italic>T. b. gambiense </italic>repertoires are nearly reciprocally monophyletic, but does not explain the very different relationships seen at ESAG2.</p> <p>A point to remember is that BESs are dynamic groups of genes. Each of the BES repertoires cloned presents a snapshot of what is going on within a population of trypanosomes, and each repertoire will contain a subset of ESAG 'alleles' present within the given trypanosome population. Further population structure analysis will allow us to determine if what we see in these analysed subspecies is a consequence of sampling rather than due to selection or consequences of loss of parts of the life-cycle. While these preliminary analyses are not the end of the story, these BES libraries will be a useful resource for further analyses of the processes influencing genetic variation in the telomeres of African trypanosomes.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>African trypanosomes (including <italic>Trypanosoma brucei</italic>) are unicellular parasites which multiply in the mammalian bloodstream. <italic>T. brucei </italic>has about twenty telomeric bloodstream form Variant Surface Glycoprotein (<italic>VSG</italic>) 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 (<italic>ESAG</italic>s) in addition to the telomeric <italic>VSG</italic>. These polymorphic <italic>ESAG </italic>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.</p> <title>Results</title> <p>Here we describe the isolation of repertoires of BES containing telomeres from three trypanosome subspecies: <italic>Trypanosoma brucei gambiense </italic>DAL 972 (causative agent of West-African trypanosomiasis), <italic>T. b. brucei </italic>EATRO 2340 (a nonhuman infective strain) and <italic>T. equiperdum </italic>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, <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2</italic>) from these three trypanosome BES repertoires.</p> <title>Conclusion</title> <p>With the exception of ESAG2, the BES sequence repertoires derived from <italic>T. b. gambiense </italic>are both less diverse than and nearly reciprocally monophyletic relative to those from <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>. Furthermore, although we find evidence for adaptive evolution in all three ESAG repertoires in <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>, only ESAG2 appears to be under diversifying selection in <italic>T. b. gambiense</italic>. This low level of variation in the <italic>T. b. gambiense </italic>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 <italic>ESAG </italic>genetic diversity.</p>
<title>Authors' contributions</title> <p>RY carried out most of the experimental work including most of the TAR cloning of the <italic>T. brucei </italic>telomere sets, and most of the sequence analysis. In addition, RY performed some of the bioinformatic analysis. JET carried out bioinformatic analysis on the data and helped write the manuscript. A.K. helped with TAR cloning, sequencing and bioinformatic analysis. MB performed TAR cloning of some of the <italic>T. b. gambiense </italic>telomere set. EJL helped with design of the experimental approach and data analysis. GR helped with the design and coordination of the project, helped with data analysis and drafted the manuscript. All authors have read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We are very grateful to Wendy Gibson (University of Bristol) for provision of the <italic>T. equiperdum </italic>STIB 818 cell line and <italic>T. brucei rhodesiense </italic>DNA samples, Keith Matthews (University of Edinburgh) for generous amounts of <italic>T. b. brucei </italic>EATRO 2340 DNA, and to Bill Wickstead, Keith Gull (University of Oxford) and Wendy Gibson for the <italic>T. b. gambiense </italic>DAL 972 cell line. We thank Alexander Fullbrook and Manish Kushwaha for help with typing trypanosome strains. We thank members of the Maiden laboratory for advice with sequence protocols. We are grateful to Tara Stanne, Mani Narayanan, Manish Kushwaha, Nadina Vasileva and Andrew Voak for comments on the manuscript. G.R. is a Wellcome Senior Fellow in the Basic Biomedical Sciences. This research was funded by the Wellcome Trust.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Isolation of trypanosome BESs in yeast using Transformation Associated Recombination (TAR) cloning as described in</bold>[##REF##15520294## 34##]. A schematic of a typical BES is shown above with the promoter indicated with a flag, and different expression site associated genes (ESAGs) indicated with coloured boxes. Characteristic repeat arrays either upstream of the BES promoter (50 bp repeats) or adjacent to the telomeric <italic>VSG </italic>gene (70 bp repeats) are indicated with vertically hatched boxes. Trypanosome telomere repeats are indicated with white arrows. The linearised pEB4 TAR vector is shown below with yeast telomere repeats (black triangles) as well as a positive (+) and a negative (-) selectable marker flanking a BES promoter containing fragment. Recombination between the vector and genomic DNA on the BES promoter fragment results in a yeast artificial chromosome (YAC) which can be stably maintained.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Sequence alignment over a particularly polymorphic region of ESAG6 shows that <italic>T. b. gambiense </italic>ESAG6 types are significantly less diverse over the hypervariable region than ESAG6 sequences from <italic>T. b. brucei </italic>and <italic>T. equiperdum</italic>.</bold> The ESAG6 protein sequence types as listed in Tables 1-3 are indicated on the left with amino acid residue position above. The ESAG6 hypervariable region (HV) as described in [##REF##1697265##38##,##REF##12664158##40##] as well as boxes I and II of the ESAG6 transferrin binding site (Tf binding) as described in [##REF##9405356##39##] are indicated with red boxes. Residues that are the most dissimilar to the consensus within a given trypanosome subspecies are highlighted in white. Residues that are completely conserved within a given subspecies are indicated in dark blue, while residues with intermediate degrees of sequence conservation are highlighted with intermediate shades of blue.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Maximum likelihood phylogenetic trees show that BES promoter and <italic>ESAG6 </italic>sequences from <italic>T. b. gambiense </italic>972 (Tbg) cluster separately from those from <italic>T. b. brucei </italic>EATRO 2340 (Tbb) and <italic>T. equiperdum </italic>STIB 818 (Teq).</bold> The DNA sequence types are as listed in Tables 1-3. Sequence accession numbers are listed in the Materials and Methods. Alignments of the BES promoter sequences used are shown in the supplementary material (Additional file ##SUPPL##1##2##). Bootstrap support values are shown for key nodes only.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Maximum likelihood phylogenetic trees of <italic>ESAG5 </italic>and <italic>ESAG2 </italic>show that <italic>ESAG5 </italic>sequences from <italic>T. b. gambiense </italic>972 (Tbg) cluster separately from those from <italic>T. b. brucei </italic>EATRO 2340 (Tbb) and <italic>T. equiperdum </italic>STIB 818 (Teq)</bold>. Most of the <italic>T. b. brucei ESAG2 </italic>sequences also appear to cluster separately. The <italic>ESAG </italic>DNA sequence types are as listed in Tables 1-3. Sequence accession numbers are listed in the Materials and Methods. Bootstrap support values are shown for key nodes only. Type A and Type B sequences referred to in the text are indicated with brackets.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p><italic>T. brucei gambiense </italic>DAL 972 BES TAR clone library. </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>BES</bold></td><td align="center"><bold>Prom. type</bold></td><td align="center" colspan="2"><bold>ESAG6 type</bold></td><td align="center"><bold>No. TAR clones</bold></td><td align="center"><bold>Represent. TAR clone</bold></td><td align="center" colspan="2"><bold>ESAG5 type</bold></td><td align="center" colspan="2"><bold>ESAG2 type</bold></td></tr><tr><td/><td/><td colspan="2"><hr/></td><td/><td/><td colspan="4"><hr/></td></tr><tr><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td></tr></thead><tbody><tr><td align="center"><bold>1</bold></td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">18</td><td align="center">A41</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>2</bold></td><td align="center">2</td><td align="center">2<sup>a</sup></td><td align="center">-</td><td align="center">1</td><td align="center">2.2–37</td><td align="center">2</td><td align="center">2</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>3</bold></td><td align="center">2</td><td align="center">3</td><td align="center">1</td><td align="center">27</td><td align="center">A120</td><td align="center">2</td><td align="center">2</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>4</bold></td><td align="center">3</td><td align="center">4</td><td align="center">1</td><td align="center">15</td><td align="center">D18</td><td align="center">3</td><td align="center">3</td><td align="center">2<sup>c</sup></td><td align="center">1</td></tr><tr><td align="center"><bold>5</bold></td><td align="center">4</td><td align="center">5</td><td align="center">1</td><td align="center">17</td><td align="center">B49</td><td align="center">4</td><td align="center">3</td><td align="center">3</td><td align="center">2</td></tr><tr><td align="center"><bold>6</bold></td><td align="center">5</td><td align="center">5</td><td align="center">1</td><td align="center">35</td><td align="center">D30</td><td align="center">4</td><td align="center">3</td><td align="center">4</td><td align="center">3</td></tr><tr><td align="center"><bold>7</bold></td><td align="center">6</td><td align="center">6</td><td align="center">2</td><td align="center">22</td><td align="center">C10</td><td align="center">4</td><td align="center">3</td><td align="center">5</td><td align="center">4</td></tr><tr><td align="center"><bold>8</bold></td><td align="center">7</td><td align="center">7</td><td align="center">3</td><td align="center">8</td><td align="center">D102</td><td align="center">5</td><td align="center">4</td><td align="center">6</td><td align="center">5</td></tr><tr><td align="center"><bold>9</bold></td><td align="center">8</td><td align="center">8</td><td align="center">4</td><td align="center">3</td><td align="center">2.1–9</td><td align="center">4</td><td align="center">3</td><td align="center">7</td><td align="center">6</td></tr><tr><td align="center"><bold>10</bold></td><td align="center">8</td><td align="center">9</td><td align="center">5</td><td align="center">24</td><td align="center">B46</td><td align="center">6</td><td align="center">5</td><td align="center">8</td><td align="center">6</td></tr><tr><td align="center"><bold>11</bold></td><td align="center">9</td><td align="center">10</td><td align="center">6</td><td align="center">16</td><td align="center">C49</td><td align="center">4</td><td align="center">3</td><td align="center">9</td><td align="center">7</td></tr><tr><td align="center"><bold>12</bold></td><td align="center">10</td><td align="center">11</td><td align="center">7</td><td align="center">17</td><td align="center">98</td><td align="center">7</td><td align="center">6</td><td align="center">10</td><td align="center">8</td></tr><tr><td align="center"><bold>13</bold></td><td align="center">11</td><td align="center">11</td><td align="center">7</td><td align="center">1</td><td align="center">A121</td><td align="center">7</td><td align="center">6</td><td align="center">10</td><td align="center">8</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center"><bold>Total</bold></td><td align="center"><bold>11</bold></td><td align="center"><bold>11</bold></td><td align="center"><bold>7</bold></td><td align="center"><bold>204</bold></td><td/><td align="center"><bold>7</bold></td><td align="center"><bold>6</bold></td><td align="center"><bold>10</bold></td><td align="center"><bold>8</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p><italic>T. brucei brucei </italic>EATRO 2340 BES TAR clone library. </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>BES</bold></td><td align="center"><bold>Prom. type</bold></td><td align="center" colspan="2"><bold>ESAG6 type</bold></td><td align="center"><bold>No. TAR clones</bold></td><td align="center"><bold>Represent. TAR clone</bold></td><td align="center" colspan="2"><bold>ESAG5 type</bold></td><td align="center" colspan="2"><bold>ESAG2 type</bold></td></tr><tr><td/><td/><td colspan="2"><hr/></td><td/><td/><td colspan="4"><hr/></td></tr><tr><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td></tr></thead><tbody><tr><td align="center"><bold>1</bold></td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">11</td><td align="center">F44</td><td align="center">1</td><td align="center">1</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>2</bold></td><td align="center">2</td><td align="center">2</td><td align="center">2</td><td align="center">9</td><td align="center">E42</td><td align="center">2</td><td align="center">2</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>3</bold></td><td align="center">3</td><td align="center">3</td><td align="center">3</td><td align="center">9</td><td align="center">G37</td><td align="center">3</td><td align="center">3</td><td align="center">3</td><td align="center">3</td></tr><tr><td align="center"><bold>4</bold></td><td align="center">4</td><td align="center">4<sup>a</sup></td><td align="center">4<sup>a</sup></td><td align="center">3</td><td align="center">E16</td><td align="center">4</td><td align="center">4</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>5</bold></td><td align="center">5</td><td align="center">5</td><td align="center">5</td><td align="center">19</td><td align="center">G39</td><td align="center">5<sup>c</sup></td><td align="center">-</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>6</bold></td><td align="center">5</td><td align="center">6</td><td align="center">6</td><td align="center">5</td><td align="center">K13</td><td align="center">-<sup>d</sup></td><td align="center">-</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>7</bold></td><td align="center">6</td><td align="center">7</td><td align="center">7</td><td align="center">1</td><td align="center">V6</td><td align="center">6</td><td align="center">5</td><td align="center">4</td><td align="center">4</td></tr><tr><td align="center"><bold>8</bold></td><td align="center">7</td><td align="center">8</td><td align="center">8</td><td align="center">8</td><td align="center">E51</td><td align="center">7</td><td align="center">6</td><td align="center">5</td><td align="center">5</td></tr><tr><td align="center"><bold>9</bold></td><td align="center">8</td><td align="center">9</td><td align="center">9</td><td align="center">7</td><td align="center">S32</td><td align="center">8</td><td align="center">7</td><td align="center">6<sup>e</sup></td><td align="center">6<sup>e</sup></td></tr><tr><td align="center"><bold>10</bold></td><td align="center">9</td><td align="center">10</td><td align="center">10</td><td align="center">4</td><td align="center">U19</td><td align="center">9</td><td align="center">8</td><td align="center">7</td><td align="center">7</td></tr><tr><td align="center"><bold>11</bold></td><td align="center">9</td><td align="center">11 and 12<sup>f</sup></td><td align="center">11 and 12<sup>f</sup></td><td align="center">1</td><td align="center">S43</td><td align="center">10</td><td align="center">9</td><td align="center">8</td><td align="center">8</td></tr><tr><td align="center"><bold>12</bold></td><td align="center">10</td><td align="center">13</td><td align="center">13</td><td align="center">16</td><td align="center">E23</td><td align="center">11</td><td align="center">10</td><td align="center">9</td><td align="center">9</td></tr><tr><td align="center"><bold>13</bold></td><td align="center">11</td><td align="center">14</td><td align="center">14</td><td align="center">8</td><td align="center">J12</td><td align="center">12</td><td align="center">11</td><td align="center">10</td><td align="center">10</td></tr><tr><td align="center"><bold>14</bold></td><td align="center">11</td><td align="center">15</td><td align="center">15</td><td align="center">9</td><td align="center">F36</td><td align="center">13</td><td align="center">12</td><td align="center">11</td><td align="center">11</td></tr><tr><td align="center"><bold>15</bold></td><td align="center">12</td><td align="center">16</td><td align="center">16</td><td align="center">17</td><td align="center">E12</td><td align="center">14</td><td align="center">13</td><td align="center">12</td><td align="center">12</td></tr><tr><td align="center"><bold>16</bold></td><td align="center">13</td><td align="center">17</td><td align="center">17</td><td align="center">2</td><td align="center">R23</td><td align="center">10</td><td align="center">9</td><td align="center">8</td><td align="center">8</td></tr><tr><td align="center"><bold>17</bold></td><td align="center">13</td><td align="center">11</td><td align="center">11</td><td align="center">2</td><td align="center">E40</td><td align="center">10</td><td align="center">9</td><td align="center">8</td><td align="center">8</td></tr><tr><td align="center"><bold>18</bold></td><td align="center">13</td><td align="center">18</td><td align="center">18</td><td align="center">4</td><td align="center">J4</td><td align="center">10</td><td align="center">9</td><td align="center">13<sup>e</sup></td><td align="center">13<sup>e</sup></td></tr><tr><td align="center"><bold>19</bold></td><td align="center">14</td><td align="center">19</td><td align="center">19</td><td align="center">13</td><td align="center">E19</td><td align="center">15</td><td align="center">14</td><td align="center">14</td><td align="center">14</td></tr><tr><td align="center"><bold>20</bold></td><td align="center">15</td><td align="center">20</td><td align="center">20</td><td align="center">12</td><td align="center">F35</td><td align="center">16</td><td align="center">15</td><td align="center">15</td><td align="center">15</td></tr><tr><td align="center"><bold>21</bold></td><td align="center">16</td><td align="center">21</td><td align="center">21</td><td align="center">16</td><td align="center">E48</td><td align="center">17</td><td align="center">16</td><td align="center">16</td><td align="center">16</td></tr><tr><td align="center"><bold>22</bold></td><td align="center">17</td><td align="center">22</td><td align="center">22</td><td align="center">31</td><td align="center">F27</td><td align="center">18</td><td align="center">17</td><td align="center">17</td><td align="center">17</td></tr><tr><td align="center"><bold>23</bold></td><td align="center">18</td><td align="center">11</td><td align="center">11</td><td align="center">1</td><td align="center">E44</td><td align="center">-<sup>b</sup></td><td align="center">-</td><td align="center">-<sup>b</sup></td><td align="center">-</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center"><bold>Total</bold></td><td align="center"><bold>18</bold></td><td align="center"><bold>22</bold></td><td align="center"><bold>22</bold></td><td align="center"><bold>208</bold></td><td/><td align="center"><bold>18</bold></td><td align="center"><bold>17</bold></td><td align="center"><bold>17</bold></td><td align="center"><bold>17</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p><italic>T. equiperdum </italic>STIB 818 BES TAR clone library. </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>BES</bold></td><td align="center"><bold>Prom. type</bold></td><td align="center" colspan="2"><bold>ESAG6 type</bold></td><td align="center"><bold>No. TAR clones</bold></td><td align="center"><bold>Represent. TAR clone</bold></td><td align="center" colspan="2"><bold>ESAG5 type</bold></td><td align="center" colspan="2"><bold>ESAG2 type</bold></td></tr><tr><td/><td/><td colspan="2"><hr/></td><td/><td/><td colspan="4"><hr/></td></tr><tr><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td/><td/><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td><td align="center"><bold>DNA</bold></td><td align="center"><bold>Protein</bold></td></tr></thead><tbody><tr><td align="center"><bold>1</bold></td><td align="center">1</td><td align="center">1<sup>a</sup></td><td align="center">-</td><td align="center">9</td><td align="center">10</td><td align="center">1<sup>b</sup></td><td align="center">-</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>2</bold></td><td align="center">2</td><td align="center">2</td><td align="center">1</td><td align="center">1</td><td align="center">3</td><td align="center">-<sup>c</sup></td><td align="center">-</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>3</bold></td><td align="center">3</td><td align="center">3<sup>d</sup></td><td align="center">-</td><td align="center">11</td><td align="center">D27</td><td align="center">2</td><td align="center">1</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>4</bold></td><td align="center">4</td><td align="center">4<sup>d</sup></td><td align="center">-</td><td align="center">2</td><td align="center">25</td><td align="center">3</td><td align="center">2</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>5</bold></td><td align="center">5</td><td align="center">5<sup>b</sup></td><td align="center">-</td><td align="center">2</td><td align="center">17</td><td align="center">-<sup>c</sup></td><td align="center">-</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>6</bold></td><td align="center">6</td><td align="center">6</td><td align="center">2</td><td align="center">8</td><td align="center">14</td><td align="center">4</td><td align="center">3</td><td align="center">3</td><td align="center">3</td></tr><tr><td align="center"><bold>7</bold></td><td align="center">6</td><td align="center">7<sup>e</sup></td><td align="center">-</td><td align="center">7</td><td align="center">D44</td><td align="center">5</td><td align="center">4</td><td align="center">4</td><td align="center">4</td></tr><tr><td align="center"><bold>8</bold></td><td align="center">6</td><td align="center">8<sup>b</sup></td><td align="center">-</td><td align="center">3</td><td align="center">16</td><td align="center">6</td><td align="center">5</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>9</bold></td><td align="center">7</td><td align="center">9</td><td align="center">3</td><td align="center">7</td><td align="center">7</td><td align="center">7</td><td align="center">6</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>10</bold></td><td align="center">8</td><td align="center">9</td><td align="center">3</td><td align="center">9</td><td align="center">1</td><td align="center">8</td><td align="center">7</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="center"><bold>11</bold></td><td align="center">9</td><td align="center">10</td><td align="center">4</td><td align="center">10</td><td align="center">53</td><td align="center">9<sup>b</sup></td><td align="center">-</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>12</bold></td><td align="center">10</td><td align="center">11</td><td align="center">5</td><td align="center">8</td><td align="center">C23</td><td align="center">9<sup>b</sup></td><td align="center">-</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>13</bold></td><td align="center">10</td><td align="center">12</td><td align="center">6</td><td align="center">1</td><td align="center">K13</td><td align="center">9<sup>b</sup></td><td align="center">-</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>14</bold></td><td align="center">11</td><td align="center">13</td><td align="center">7</td><td align="center">4</td><td align="center">G47</td><td align="center">6</td><td align="center">5</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td align="center"><bold>15</bold></td><td align="center">11</td><td align="center">-<sup>f</sup></td><td align="center"><bold>-</bold></td><td align="center">8</td><td align="center">E6</td><td align="center">10</td><td align="center">8</td><td align="center">2</td><td align="center">2</td></tr><tr><td align="center"><bold>16</bold></td><td align="center">12<sup>g</sup></td><td align="center">-<sup>f</sup></td><td align="center">-</td><td align="center">1</td><td align="center">J46</td><td align="center">11</td><td align="center">9</td><td align="center">-<sup>c</sup></td><td align="center">-</td></tr><tr><td colspan="10"><hr/></td></tr><tr><td align="center"><bold>Total</bold></td><td align="center"><bold>12</bold></td><td align="center"><bold>13</bold></td><td align="center"><bold>7<sup>h</sup></bold></td><td align="center"><bold>91</bold></td><td/><td align="center"><bold>11</bold></td><td align="center"><bold>9</bold></td><td align="center"><bold>4</bold></td><td align="center"><bold>4</bold></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>BES nucleotide diversity. </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>Locus</bold></td><td align="left"><bold>Strain</bold></td><td align="center"><bold>N° seqs</bold></td><td align="center"><bold>Length</bold></td><td align="right"><bold>Polymorphic sites % S</bold></td><td align="center"><bold>Average nt. Diversity π</bold></td><td align="center"><bold>N° rec. break points</bold></td><td align="right"><bold>Tract length</bold></td></tr></thead><tbody><tr><td align="center"><bold>Promoter</bold></td><td align="left">T. b. g.</td><td align="center">13</td><td align="center">621</td><td align="right">5.6</td><td align="center">0.015</td><td align="center">1</td><td align="right">311</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">23</td><td/><td align="right">24.3</td><td align="center">0.058</td><td align="center">2</td><td align="right">207</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">15</td><td/><td align="right">13.8</td><td align="center">0.046</td><td align="center">2</td><td align="right">207</td></tr><tr><td align="center"><bold><italic>ESAG6</italic></bold></td><td align="left">T. b. g.</td><td align="center">12</td><td align="center">1197</td><td align="right">6.0</td><td align="center">0.012</td><td align="center">1</td><td align="right">599</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">23</td><td/><td align="right">17.3</td><td align="center">0.051</td><td align="center">5</td><td align="right">200</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">10</td><td/><td align="right">8.4</td><td align="center">0.030</td><td align="center">3</td><td align="right">299</td></tr><tr><td align="center"><bold><italic>ESAG5</italic></bold></td><td align="left">T. b. g.</td><td align="center">13</td><td align="center">1429</td><td align="right">1.5</td><td align="center">0.002</td><td align="center">0</td><td align="right">1429</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">21</td><td/><td align="right">20.0</td><td align="center">0.055</td><td align="center">7</td><td align="right">179</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">14</td><td/><td align="right">14.6</td><td align="center">0.061</td><td align="center">7</td><td align="right">179</td></tr><tr><td align="center"><bold><italic>ESAG2</italic></bold></td><td align="left">T. b. g.</td><td align="center">11</td><td align="center">1291</td><td align="right">15.3</td><td align="center">0.047</td><td align="center">7</td><td align="right">161</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">19</td><td/><td align="right">28.7</td><td align="center">0.101</td><td align="center">9</td><td align="right">129</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">10</td><td/><td align="right">16.9</td><td align="center">0.044</td><td align="center">9</td><td align="right">129</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>BES amino acid diversity. </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Locus</bold></td><td align="left"><bold>Strain</bold></td><td align="center"><bold>N° seqs</bold></td><td align="center"><bold>Codons</bold></td><td align="right"><bold>Polymorphic Sites % S</bold></td><td align="center"><bold>Average a.a. diversity π</bold></td><td align="center"><bold>dN/dS ratio ω &gt; 2</bold></td><td align="center"><bold>Average dN/dS ω (avg)</bold></td></tr></thead><tbody><tr><td align="left"><bold>ESAG6</bold></td><td align="left">T. b. g.</td><td align="center">12</td><td align="center">401</td><td align="right">9.2</td><td align="center">0.018</td><td align="center">15</td><td align="center">3.40</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">23</td><td/><td align="right">24.7</td><td align="center">0.079</td><td align="center">24</td><td align="center">3.20</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">8</td><td/><td align="right">9.2</td><td align="center">0.040</td><td align="center">38</td><td align="center">2.15</td></tr><tr><td align="left"><bold>ESAG5</bold></td><td align="left">T. b. g.</td><td align="center">13</td><td align="center">469</td><td align="right">2.8</td><td align="center">0.005</td><td align="center">0</td><td align="center">n.a.</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">19</td><td/><td align="right">27.5</td><td align="center">0.087</td><td align="center">48</td><td align="center">3.75</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">10</td><td/><td align="right">19.4</td><td align="center">0.087</td><td align="center">49</td><td align="center">2.29</td></tr><tr><td align="left"><bold>ESAG2</bold></td><td align="left">T. b. g.</td><td align="center">11</td><td align="center">424</td><td align="right">19.1</td><td align="center">0.059</td><td align="center">15</td><td align="center">3.38</td></tr><tr><td/><td align="left">T. b. b.</td><td align="center">19</td><td/><td align="right">37.0</td><td align="center">0.130</td><td align="center">11</td><td align="center">3.98</td></tr><tr><td/><td align="left">T. eq.</td><td align="center">10</td><td/><td align="right">19.6</td><td align="center">0.053</td><td align="center">11</td><td align="center">4.39</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Sup. Figure 1</bold>. Multiplex PCR typing of trypanosome DNA to establish the presence or absence of the Serum Resistance Associated gene <italic>SRA</italic>. Multiplex PCR was performed using the primer sets and conditions of Picozzi et al [##REF##17643434##13##]. Lanes indicate PCR reactions using no genomic DNA (lane 1) or genomic DNA from <italic>Trypanosoma brucei brucei </italic>427 (lane 2), <italic>T. b. brucei </italic>TREU 927/4 (lane 3), <italic>T. b. gambiense </italic>DAL 972 (lane 4), <italic>T. equiperdum </italic>STIB 818 (lane 5), <italic>T. b. brucei </italic>EATRO 2340 (lane 6), <italic>T. b. rhodesiense </italic>LVH 108 (lane 7) or <italic>T. b. rhodesiense </italic>WB58 (lane 8). PCR products amplifying the GPI-PLC gene (PLC), the SRA gene (SRA) or the SRA-like VSG (VSG SRA) are indicated on the right with arrows. A DNA ladder is on the left with sizes indicated in base pairs (bp).</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>Sup. Figure 2</bold>. Sequence alignment of BES promoter sequences analysed in this study. BES promoter sequences isolated from <italic>T. b. gambiense </italic>DAL 972, <italic>T. b. brucei </italic>EATRO 2340 or <italic>T. equiperdum </italic>were aligned using Vector NTI. Sequence types as indicated in Tables ##TAB##0##1##, ##TAB##1##2##, ##TAB##2##3## are indicated on the left.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>Sup. Figure 3</bold>. Sequence alignments of ESAG6 sequences analysed in this manuscript.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p><bold>Sup. Figure 4</bold>. Sequence alignments of ESAG2 sequences analysed in this manuscript.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p><bold>Sup. Figure 5</bold>. Sequence alignments of ESAG5 sequences analysed in this manuscript.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p><bold>Supplementary Table 2</bold>. Non-synonymous substitution rates for different genes located within trypanosome BESs isolated from <italic>Trypanosoma brucei gambiense </italic>DAL 972, <italic>T. b. brucei </italic>EATRO 2340, or <italic>T. equiperdum </italic>STIB 818.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p><bold>Supplementary Table 3</bold>. Table with the ratio of nonsynonymous-to-synonymous substitutions (ω) along ESAG6, ESAG5 and ESAG2.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p><bold>Sup. Figure 6</bold>. The dN/dS ratio (ω) of ESAG6, ESAG5 or ESAG2 calculated from sequence repertoires from <italic>T. b. gambiense</italic>, <italic>T. b. brucei </italic>and <italic>T. equiperdum </italic>plotted against codon number.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p><bold>Supplementary Table 1</bold>. Table of primers used.</p></caption></supplementary-material>
<table-wrap-foot><p>Overview of <italic>T. b. gambiense </italic>DAL 972 TAR clone library. A total of 204 <italic>T. b. gambiense </italic>TAR clones were isolated in yeast and typed into 13 different BES sets after sequencing over the promoter and <italic>ESAG6 </italic>regions. The different BES promoter and ESAG6 sequence types (either at the DNA or protein level) are indicated with numbers. Two TAR clones were chosen from each BES type (the name of one representative clone is indicated in the chart), and the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>ORFs were isolated and sequenced from both of these two clones. The sequence types of the ESAG5 and ESAG2 are indicated.</p><p><sup>a</sup>Partial sequence.</p><p><sup>b</sup>Presence not detected by PCR. See Materials and Methods for diagnostic primers that were used for each <italic>ESAG</italic>.</p><p><sup>c</sup>Approximately 300 bp is single stranded sequence.</p></table-wrap-foot> <table-wrap-foot><p>Overview of <italic>T. b. brucei </italic>EATRO 2340 TAR clone library. A total of 208 TAR clones were isolated and typed into 23 different BES sets. The sequence typing proceeded further as described for Table 1. As before, two TAR clones were chosen for each BES type (the name of one clone is indicated in the chart), and the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>ORFs were isolated and sequenced from both of these two clones.</p><p><sup>a</sup>First 1000 bp only.</p><p><sup>b</sup>Presence not detected by PCR. See Materials and Methods for diagnostic primers that were used for each <italic>ESAG</italic>.</p><p><sup>c</sup>Pseudogene.</p><p><sup>d</sup>Present but not sequenced.</p><p><sup>e</sup>One clone.</p><p><sup>f</sup>Two <italic>ESAG6 </italic>in one BES, both assembled.</p></table-wrap-foot> <table-wrap-foot><p>Overview of <italic>T. equiperdum </italic>STIB 818 TAR clone library. A total of 91 TAR clones were isolated and typed into 16 different BES sets. The sequence typing was performed as described for Table 1. As done previously, two TAR clones were chosen for each BES type (the name of one clone is indicated in the chart), and the full length <italic>ESAG6</italic>, <italic>ESAG5 </italic>and <italic>ESAG2 </italic>ORFs were isolated and sequenced from both of these two clones.</p><p><sup>a</sup>Double sequence in the centre indicating two sequence types.</p><p><sup>b</sup>Pseudogene.</p><p><sup>c</sup>Presence not detected by PCR. See Materials and Methods for diagnostic primers that were used for each <italic>ESAG</italic>.</p><p><sup>d</sup>Partial sequence.</p><p><sup>e</sup>Double sequence indicating two sequence types within the same BES.</p><p><sup>f </sup>Sequence failed.</p><p><sup>g</sup>Anti-sense sequence strand only.</p><p><sup>h</sup>Probably an underestimate.</p></table-wrap-foot> <table-wrap-foot><p>Nucleotide diversity and recombination in BES sequences from three trypanosomatid protozoa. The BES loci studied are indicated on the left. The strains are <italic>Trypanosoma brucei gambiense </italic>DAL 972 (T.b.g.), <italic>T. b. brucei </italic>EATRO 2340 (T.b.b.) and <italic>T. equiperdum </italic>STIB 818 (T. eq.). N° seqs. = number of nucleotide sequences in alignment; L = number of nucleotides in gap-stripped alignment; %S = percentage of polymorphic sites; π = average pairwise nucleotide diversity; N° rec. break points and tract length = number of recombination breakpoints and mean tract length inferred by GARD analysis.</p></table-wrap-foot> <table-wrap-foot><p>Amino acid diversity and selection in ESAG sequences from three trypanosomatid protozoa. The BES loci studied are indicated on the left. The strains are <italic>Trypanosoma brucei gambiense </italic>DAL 972 (T.b.g.), <italic>T. b. brucei </italic>EATRO 2340 (T.b.r.) and <italic>T. equiperdum </italic>STIB 818 (T. eq.). N° seqs. = number of protein sequences in alignment; codons = number of codons in gap-stripped alignment; % S = percentage of polymorphic sites; π = average pairwise amino acid diversity; dN/dS ratio ω &gt; 2 = number of amino acid residues with dN/dS ratio (ω) estimated to be greater than 2 indicating evidence for positive selection; Average dN/dS ω (avg.) = average ω estimate over residues with ω &gt; 2.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-385-1"/> <graphic xlink:href="1471-2164-9-385-2"/> <graphic xlink:href="1471-2164-9-385-3"/> <graphic xlink:href="1471-2164-9-385-4"/>
<media xlink:href="1471-2164-9-385-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S7.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S8.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-385-S9.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[]
BMC Genomics. 2008 Aug 12; 9:385