A Major Genetic Locus in Trypanosoma brucei Is a Determinant of Host Pathology

The progression and variation of pathology during infections can be due to components from both host or pathogen, and/or the interaction between them. The influence of host genetic variation on disease pathology during infections with trypanosomes has been well studied in recent years, but the role of parasite genetic variation has not been extensively studied. We have shown that there is parasite strain-specific variation in the level of splenomegaly and hepatomegaly in infected mice and used a forward genetic approach to identify the parasite loci that determine this variation. This approach allowed us to dissect and identify the parasite loci that determine the complex phenotypes induced by infection. Using the available trypanosome genetic map, a major quantitative trait locus (QTL) was identified on T. brucei chromosome 3 (LOD = 7.2) that accounted for approximately two thirds of the variance observed in each of two correlated phenotypes, splenomegaly and hepatomegaly, in the infected mice (named TbOrg1). In addition, a second locus was identified that contributed to splenomegaly, hepatomegaly and reticulocytosis (TbOrg2). This is the first use of quantitative trait locus mapping in a diploid protozoan and shows that there are trypanosome genes that directly contribute to the progression of pathology during infections and, therefore, that parasite genetic variation can be a critical factor in disease outcome. The identification of parasite loci is a first step towards identifying the genes that are responsible for these important traits and shows the power of genetic analysis as a tool for dissecting complex quantitative phenotypic traits.


Introduction
Forward genetics provides a powerful tool for analysing phenotypes and identifying genes that are responsible for a number of important traits. The importance of the linkage mapping approach is particularly appropriate for the analysis of phenotypes for which there is no obvious candidate gene, or when more than one gene is predicted to be involved in a phenotype, so called complex traits or quantitative trait loci. Such quantitative traits include disease severity, where there is a range of pathogenesis phenotypes caused by a particular pathogen.
African trypanosomes cause a disease syndrome of high morbidity and mortality across large areas of sub-Saharan Africa in both humans (sleeping sickness) and domesticated animals (Nagana). The parasites are spread by blood-feeding tsetse fly vectors (Glossina ssp.), which inject the organisms into the mammalian host's bloodstream where they replicate extracellularly resulting in a chronic wasting condition. Two subspecies of Trypanosoma brucei, T. b. rhodesiense and T. b. gambiense, are able to infect humans causing a disease that is normally fatal and affects approximately 70,000 people per year [1] (although this number is undoubtedly a gross underestimate [2]). In livestock, tens of millions of animals are affected each year with T. brucei and the related pathogens T. congolense and T. vivax [3].
In any host-pathogen relationship, variation in disease outcome can arise from differences between either hosts, pathogens, or both. In trypanosome biology, variation in parasite virulence has been well documented but the genetic basis for this has been largely unexplored. Classically, the two T. brucei subspecies have been described as causing different pathology; T. b. rhodesiense causes a short, acute disease, while that seen with T. b. gambiense is more chronic and less severe [4]. However, the clinical distinctions are unquestionably less clear than textbooks suggest. Laboratory experiments have demonstrated differences in pathogenicity during mouse infections with different strains of T. b. gambiense [5,6], and recently it has been shown that there is markedly different pathology elicited in terms of disease severity and progression between T. b. rhodesiense patients in geographically distinct areas in Uganda and Malawi [7,8]. In the major cattle pathogens, T. congolense and T. vivax, variation between isolates in clinical pathology are also well documented [9,10,11].
Variation in parasite virulence phenotypes that can be directly attributed to trypanosome genetic variation has been experimentally demonstrated in mice in several studies using T. brucei and T. congolense [12,13]. Thus, variation is well known and clinically important. The genetic basis for this variation however, has not been investigated. This is in marked contrast to the very considerable effort that has been directed towards dissecting the genetic basis behind variation in pathology ('trypanotolerance') attributable to the mammalian host in response to infections with the major veterinary pathogen, T. congolense in both mice [14,15] and cattle [16]. These studies have resulted in the identification of loci in both host systems that contribute to the control of infection.
As all three major African trypanosome species are responsible for a wide range of virulence phenotypes, understanding the genetic determinants of this variation is an essential factor to be integrated into any model of pathogenesis of trypanosomiasis [7,12,17]. Therefore, the contribution of the host to the control of disease and the parasite to the pathogenesis must both be examined in order to produce a holistic picture of host-parasite interactions and the survival of the host. Observed differences in virulence between trypanosome strains point to a genetic basis for the spectrum of disease caused in the host and therefore provide a route for identifying the parasite factors that cause disease in the mammalian host. The use of a classical genetic linkage mapping approach, in which the inheritance of phenotypic traits are analysed in progeny of genetic crosses and examined for co-segregation with genetic markers, is a valuable route to identifying genes and loci that has been used for a number traits in pathogenic parasites. This approach is dependent upon the development of a genetic map for each species. Recently, genetic maps have been generated for a number of haploid parasites, for example, Plasmodium falciparum [18], Plasmodium chabaudi [19], Toxoplasma gondii [20], and Eimeria tenella [21], opening up the possibility of using classical genetic analysis to identify genes involved in important phenotypes. For example, the identification of protective antigens that enable survival upon exposure to the host immune response in Eimeria species [22] and P. falciparum [23], a gene that influences the ability to invade red blood cells in P. falciparum [24], and genes that confer resistance to chloroquine and quinine in P. falciparum [25,26]. The haploid nature of these organisms means that the contribution of a single allele to a phenotype can be measured in isolation without the complication of an effect of a second allele, allowing for the identification of genes involved in various phenotypes via a linkage mapping approach with remarkably few progeny.
Perhaps the most compelling case for the power of genetic analysis is that of parasite virulence in T. gondii, whereby Quantitative Trait Loci (QTL) associated with virulence were identified [20], which directly led to the identification of secreted kinases as the key genes involved in strain specific pathogenesis of toxoplasmosis [27,28,29]. Indeed, the relevance of these genes has been confirmed in T. gondii field strains [30] and these latter studies elegantly underline the application of the classical genetic approach to identifying genes of large effect in parasites [31,32]. The development of genetic maps for two sub-species of T. brucei, T. b. brucei and T. b. gambiense [33,34] provides the potential for using a classical genetic approach in this organism, although the diploid nature of this parasite means that only phenotypes that are heterozygous in the parental strain for which the map was generated can be mapped using F1 progeny.
Pathology in trypanosome-infected hosts is multi-faceted and we have chosen to focus on four indices of pathology; anaemia, reticulocytosis, splenomegaly and hepatomegaly. These traits are all important in both humans and cattle in infections with all three species of trypanosome [11,35,36]. Here, we report a genetic analysis of these inherited host pathology traits and show that two, splenomegaly and hepatomegaly, resulted in a highly significant QTL on chromosome 3 of the T. brucei genome.

Ethics statement
The maintenance and care of experimental animals complied with the appropriate legislation; the UK Animals (Scientific Procedures) Act, 1986, and with the national and University of Glasgow maintenance and care guidelines.

Trypanosomes
A panel of 39 independent F1 T. brucei progeny clones had been previously generated from a genetic cross between a strain originally isolated from a tsetse fly (TREU927) and a second isolated from a hartebeeste (STIB247) [37,38]. Parental and progeny clones were previously genotyped for 182 microsatellite markers [33] and thirty one progeny were used in the current study (Fig. 1). The eight further clones isolated from the TREU927/STIB247 cross are only available at present as the insect life cycle stage procyclic forms and are thus not suitable for in vivo mammalian studies. The two parental lines were each transmitted through tsetse flies and recloned (using methods described in [37]) so that their transmission/passage history was very similar to that of the hybrid progeny. All parasites used readily infect rodents and passages in our laboratory are kept to a maximum of six times in mice from tsetse transmission, in order to maintain tsetse transmissibility and minimise selection by growth in rodents [39]. Five mini-and microsatellite markers, all on different chromosomes of T. brucei, were used for routine genotyping of all infections (Ch1/MS42, Ch2/PLC, Ch3/IJ15/ 1, Ch5/JS2 and Ch9/49; for primer sequences and PCR conditions see [33]).

Animal infections
Parasites were grown from stabilates (cryopreserved in liquid nitrogen) in a donor ICR mouse (Harlan, UK). Parasites in

Author Summary
Trypanosomes are single-celled organisms that are transmitted between animal hosts by the tsetse fly. These parasites infect a wide range of mammals and in sub-Saharan Africa are extensively debilitating to livestock, and some species are also able to infect humans causing a disease, sleeping sickness, that is usually fatal unless treated. Some trypanosome strains cause more severe disease than others, and studying these differences may allow the identification of how serious disease is caused. We approached this problem by looking at how differences in disease symptoms (enlarged spleen and liver, and reduced blood cell numbers) that are caused in infections in mice with two strains of Trypanosoma brucei, TREU927 and STIB247. These disease manifestations are clinically relevant in human and livestock trypanosome infections. Examining how the symptoms are inherited in infections with offspring of a cross between the two strains allowed the identification of a region of the T. brucei genome that contains a gene (or several genes) that contributes significantly towards the enlarged spleen and liver observed in infected mice. This is a first step towards identifying the parasite genes that cause disease in the host (virulence factors), which may provide routes for developing novel therapies against the disease.
logarithmic growth were harvested from the donor mouse by terminal exsanguination. The trypanosomes were counted in triplicate in an improved Neubauer haemocytometer, and diluted in Carter's Balanced Salt Solution (CBSS) to 1610 4 trypanosomes per 0.2 ml inoculum. Parasites were then inoculated via the intraperitoneal route into experimental BALB/c mice (Harlan, UK). Five mice were infected per strain, and an equivalent number of uninfected mice were included as controls. Infections with all of the progeny were completed in nine different batches and, in each case, infections with both parental strains were used to ensure consistency of phenotypes across all of the mouse batches, and to provide a means of minimising batch-to-batch variation (see below). At day 10 post-infection, mice were euthanased. Day 10 was chosen as it was the time point at which the greatest number of pathogenesis phenotypes were significantly different between the parental strains [Morrison et al., submitted].

Phenotype measurements
We chose to measure two features that can be readily related to known clinical pathology in multiple host speciesorganomegaly and anaemia. At day 10 post-infection, mice were euthanased and the spleen and liver dissected and weighed. After weighing, the spleens were snap frozen in liquid nitrogen, and stored at 280uC. A FACs based assay was developed to measure anaemia in the mouse model [Morrison et al., submitted]. Briefly, parameters for FACs analysis were determined for each cell population (trypanosomes, red blood cells, reticulocyctes and white blood cells), by assaying each population in isolation with a known number of fluorescent beads, allowing the gating and calculation of cell populations from whole infected blood. For each strain 5 ml blood was collected and added to 2 ml of in CBSS containing 100 U/ml heparin, followed by the addition of 198 ml 1% paraformaldehyde in PBS pH 7.4, mixed thoroughly to ensure no aggregates formed and fixed for 30 minutes at room temperature. To stain reticulocytes, 745 ml thiazole orange (100 ng/ml in PBS pH 7.4) was added to the blood/fixative mixture and allowed to stain for 1 hour in the dark at room temperature [40]. Fluorescent beads (Countbright TM absolute counting beads for flow cytometry, Invitrogen, UK) were prepared by diluting to 1000 beads/ml and ensuring even mixing by sonication. Fifty ml of the bead mixture was then added to the blood/thiazole orange preparation (total volume 1 ml; 50,000 beads/ml). FACS was carried out using a Becton Dickinson FACSCalibur, using detector FL1-H for thiazole orange. Data were analyzed using CellQuest version 3.3. For each sample 10,000 events were counted five times. This assay differs from that recently described by Antoine-Moussiaux et al [41] and is described in full in Morrison et al (submitted).

Data and QTL analysis
The data for all four phenotypes were normally distributed (D'Agostino-Pearson normality test, p.0.05 in all cases), and this allowed calculation of Pairwise correlations of phenotype data using Pearson correlation coefficient and a two-tailed P-value, and from this a value of r 2 was generated. For QTL analysis, the mean phenotype value for each progeny clone was calculated as a percentage of the phenotype for the TREU927 parental strain in that particular batch to take batch-to-batch variation into account. QTL analysis was performed using MapManager QTX software [42] for the T. brucei genetic map [33]. The variance was calculated from the five biological replicates of each phenotype for each progeny. The Likelihood Ratio Statistic (LRS) significance values were calculated by 1000 random permutations of the phenotypes relative to the genotypes that were intrinsic to the experimental data set (to obtain the equivalent LOD score, the LRS is divided by approximately 4.6). Thresholds of statistical significance (x 2 statistic) were calculated from each permutation test for each phenotype. The calculated threshold values were in line with those suggested for genome-wide scan studies [43]. For a simplified explanation of the inheritance model and relevant features of the TREU927 6 STIB247 cross see Fig S1.

Marker generation and update of the genetic map
An additional genetic marker, (TB3/19), was developed for T. brucei chromosome 3 based on the allele specific amplification of a 1 kb fragment of DNA of the TREU927 parental strain and the failure to amplify the other TREU927 allele or STIB247 parental alleles, due to a single nucleotide polymorphism in the primer site(s). As a positive control for DNA integrity, primers directed against a second non-polymorphic 1 kb fragment were used (TB3/ 20), which amplified both TREU927 alleles. This allele-specific PCR (TB3/19) segregated in the progeny and was incorporated into the genetic map. Primer sequences were Ch3/19A TGA-GTTCCTCTTGCACTCCC, Ch3/19B TCTTTACGTGT-GCGCGCTAGG, Ch3/20A GATGTAATGTCCCTTCGGA-TTGCG and Ch3/20B GTGTCTTGTAGTTTATGACGGC; PCR conditions were described previously [33]; Genotyping gaps in the original map were also filled (for version 2 of the genetic map, see http://tinyurl.com/trypmap).

Phenotype measurements
The trypanosomes strains TREU927 and STIB247 consistently differ in the pathogenesis they induce, including the degree to which they cause anaemia, reticulocytosis, splenomegaly and hepatomegaly (Fig 1). We measured these parameters in nine batches of mice and consistently observed the same pattern with greater splenomegaly and reticulocytosis in 247-infected mice and marked anaemia in 927-infected mice (Fig 1 and Fig S2). Whilst the patterns were consistent, the absolute values showed some batch-to-batch variation and to account for this, subsequent data are presented as percentages of the 927 parental levels. A detailed analysis of hepatomegaly development over multiple time points permitted robust statistical analysis by 2-way ANOVA and showed hepatomegaly to be greater in STIB247-infected mice [Morrison et al, submitted], although the data for day 10 show a smaller difference than was observed for the other three phenotypes ( Fig  S2).
In order to investigate if there is a genetic basis for this variation in these phenotypes, experimental mouse infections were carried out with 31 independent F1 progeny from a genetic cross between TREU927 and STIB247. The mean data for each phenotype for each progeny clone are illustrated in Figure 2, revealing a continuous distribution for each phenotype. All four phenotypes show marked transgressive segregation with between 10 and 20 hybrids displaying average values outside the parental range. Interestingly, few values fall below that observed for TREU927 and a large number of the progeny display values far greater than either parent. This heterosis was particularly marked for splenomegaly, hepatomegaly and reticulocytosis but less so for red cell numbers. All four phenotypes show classical quantitative inheritance with no clear segregation into different classes in the progeny, but the progeny clones at the top and bottom ends of each distribution show clear and reproducible differences, suggesting that parasite derived genetic determinants contribute to each phenotype.
Pairwise correlations were calculated for all phenotypes, revealing three combinations that show significant correlation (Fig 3). The most significant correlation was between red blood cell numbers and percentage reticulocytes, which suggested that 57% of the variance is shared between the two phenotypes. There were also significant correlations between liver and spleen weights (44% of shared variance) and reticulocyte percentage and liver weight (21%).

Linkage analysis
The relatively continuous distribution of phenotypic data from the F1 progeny in all cases suggests that multiple loci are involved in determining the phenotype. In order to investigate the genetic basis determining these traits a quantitative linkage mapping approach was therefore applied to the data using the available genetic map. The T. brucei genetic map created from a TREU927 6 STIB247 cross [33] comprising 196 markers spread across the 11 megabase chromosomes (but not covering subtelomeric regions, or mini -and intermediate chromosomes), was used to analyse whether the phenotypes (spleen weight, liver weight, reticulocyte percentage and red blood cell numbers) co-segregated with genetic markers, using QTL analysis.
Splenomegaly and hepatomegaly. Two phenotypes, splenomegaly and hepatomegaly, showed evidence for highly significant QTLs, mapping to the same region on chromosome 3 with equivalent LOD scores of 7.2 and 7.5, respectively (Fig 4; Table 1). The hepatomegaly QTL had a significance level of p#0.001 for a genome-wide analysis, and was estimated to account for approximately 64% of the phenotype variance. The splenomegaly QTL was also highly significant (p#0.001), and contributed a similar amount (66%) to the phenotype variance. Interval mapping confirmed in both cases that the QTL lies between markers TB3/3 and TB3/10, a region that covers 56.3 cM (physical region between 462529 bp and 1460975 bp, a physical distance of 998446 bp) and which contains 383 genes (Fig 5). For hepatomegaly there is no suggestion of more than one locus involved within the interval of the QTL. However, for splenomegaly although the QTL remains significant between TB3/3 and TB3/10, there are two distinct regions within the QTL showing highly significant linkage (between TB3/3 and TB3/5, and TB3/6 and TB3/10, respectively, Fig 5). This may indicate that there is more than one locus involved in splenomegaly within the context of the QTL (and indeed this could also be the case for hepatomegaly), but this will require further work to fully elucidate. As well as the linkage described above, the hepatomegaly phenotype also shows significant Measurements taken were A) spleen weight, B) liver weight, C) red blood cell numbers, D) reticulocyte percentage (error bars are 95% confidence interval). Phenotypes are expressed as percentage of the value of the mice infected with TREU927 in the respective batch of infections, and are arranged in ascending order. The mean value for the STIB247 infections is shown (grey bar) and the TREU927 value (100%) is indicated by a solid line. Progeny names are as in the TREU927 6STIB247 genetic map (see [38] Table 1).

Anaemia (red blood cells and reticulocytosis).
Linkage analysis of percentage reticulocytes revealed significant linkage to three loci on chromosomes 2, 8 and 11, but linkage hinged upon a  (Fig 4; Table 1). Interestingly these overlapping minor QTLs on chromosome 2 appear to be significantly linked to three of the four phenotypes measured in this study: hepatomegaly, splenomegaly and reticulocytosis (Fig 4;  Table 1).
For the red blood cell phenotype, two loci were identified as being significantly linked on chromosomes 10 and 11. However,  Table 1).
In summary, we have identified one highly significant locus on T. brucei chromosome 3 that accounts for approximately two thirds of the variance observed in splenomegaly and hepatomegaly, which we name TbOrg1 (for T. brucei organomegaly QTL 1). The significance of the QTL on chromosome 3 is reinforced by the fact that both phenotypes show a strong correlation with each other (Fig 3) and link to this region. There is also a smaller, but still significant, locus on chromosome 2 that contributes to strainspecific splenomegaly, hepatomegaly and reticulocytosis, named TbOrg2. The importance of TbOrg2 is enhanced by the fact that three phenotypes are all significantly linked to the locus (Fig 4), contributing approximately 30-35% to each of the phenotypes, and this therefore represents a parasite locus that has a minor but multifactorial effect on pathogenesis in the mouse host. We conducted a secondary screen after fixing TbOrg1 on chromosome 3, for both splenomegaly and hepatomegaly, but this did not reveal any further significant loci (data not shown).

Discussion
We have used genetic linkage analysis to investigate a variety of pathogenesis phenotypes in infections with T. brucei. Our approach, using genetically distinct parasite strains in one inbred host background effectively means that observed differences in infection profile must originate from parasite genotype differences. Differences caused by host factors also exist, but have been minimised by the use of inbred mice and controlling for any batchto-batch variation. Our results clearly show that variation in pathology has a parasite determined genetic basis and that using a linkage mapping approach it was possible to identify a region of the genome on chromosome 3 that has a major effect on the gross pathology in mice. This study is an important step, therefore, towards identifying the parasite virulence factors that underlie differential pathogenesis.
The phenotype data in this study do not segregate in a manner that would indicate they are determined by a single gene, unlike that of previous studies in haploid parasites, where traits in progeny could definitively be scored as (for example) virulent, intermediate or non-virulent [20], or drug resistant or sensitive [25]. The phenotypes we have examined are, or undoubtedly can be, the end product of multiple processes, and therefore will have multiple genes contributing to the variance. They are therefore true 'quantitative traits', and for this reason QTL analysis was the most appropriate tool to search for linkage between trypanosome genotype and host infection phenotype. The identification of a highly significant QTL, despite the complex nature of the phenotypes, suggests that we have identified a parasite locus that determines a process responsible for a large proportion of the variance observed.
The correlation of liver and spleen weights, and linkage of each to the same region of chromosome 3, means that it is reasonable to assume that the same locus is responsible for the differential pathogenesis in both phenotypes. Whether this is due to a parasite gene product directly influencing both processes individually, or whether the gene product defines the polarisation or sequence in which multiple downstream host and parasite pathways proceed, is one of the questions that remain to be answered. Certainly the spleen and liver are significant organs with respect to the response to infection, and splenomegaly is one of the classical clinical signs observed in experimental cases of animal trypanosomiasis and also in natural cases of animal and human disease [35,44,45]. The spleen and liver obviously both have roles with respect to the innate immune response [46] and, with respect to trypanosome infections, the direction in which polarisation of the immune response proceeds seems to determine the progression of disease [47,48,49,50]. We have previously shown that in infections with the two parental strains, pathways of the innate immune response, macrophage polarisation (IL10, RXR/LXR pathways), and alternative macrophage activation are the most significantly differentially regulated [Morrison et al., submitted]. The strain specific pathogenesis is therefore not only due to a difference in magnitude of the same pathways, but distinct pathways. It seems reasonable, therefore, to postulate that there is a parasite gene product that determines this polarisation, and the gene is located within the locus on chromosome 3. The confidence interval of the major chromosome 3 QTL for splenomegaly and hepatomegaly (TbOrg1) includes 383 genes. Of the genes within the QTL, 206 are conserved hypothetical proteins, 109 have function assigned by sequence homology, 62 have a predicted signal peptide, 10 are sequence orphans, 3 are pseudogenes and only 16 have been experimentally characterised (Table S1). Given the number of candidate genes, a fine-mapping approach is required to narrow down the region of linkage and reduce the number of genes to a level where candidates can be experimentally examined. There is scope for this approach, as the QTL interval for both phenotypes contains 12 crossovers in the 31 progeny, and therefore informative markers will further define the linkage boundaries. Fine mapping, in which more informative genetic markers within the QTL region are identified, can be undertaken in a number of ways, the most definitive of which will be to sequence the genomes of the parental strains to identify all single nucleotide polymorphisms within the locus, which can then be used to define the crossover boundaries in the progeny more precisely. There are additional and complementary techniques that will also aid in defining candidate genes, for example; the generation of further progeny in order to incorporate a greater number of informative crossovers in the genetic map and the identification of genes that are heterozygous for 927 (by identifying heterozygous single nucleotide polymorphisms (SNPs) in the relevant chromosome sequence) within the QTL (given the necessity for the gene to be heterozygous to obtain segregation in this cross). Additionally, gene expression data could be used to rule out, for example, those genes only expressed in the procyclic life cycle stage, and inferences of gene function could also be made from the genome sequence data. Similar post-genomic approaches were used successfully to identify T. gondii virulence factors, when starting from QTLs spanning regions of a similar magnitude [27]. Additionally, while the approach using F1 progeny clearly allows the identification of dominant heterozygous loci (Fig S1), a limitation is that it will not detect co-dominant or recessive loci. Therefore, using the same approach, but with an F2 cross, would allow the dissection of further loci that are not detectable with the current approach, and provide a more comprehensive picture of the parasite genes contributing to pathogenesis. However, no successful F2 cross has been reported in trypanosomes thus far.
The QTLs identified on chromosome 2 for splenomegaly, hepatomegaly and reticulocytosis (TbOrg2) are potentially interesting. In all three cases there are multiple markers linked and the loci of the three phenotypes partially overlap, suggesting that the correlations between the traits (two of which are highly significant and the third approaches significance) may have some causal basis worthy of further investigation. In contrast, there are no QTLs in common for red blood cell numbers and reticulocytosis despite the very strong correlation between these two traits. One explanation for this is that the spleen, and to a lesser extent the liver, are haematopoietic organs, and therefore linked to reticulocyte production. In contrast, the change in red blood cell numbers will largely be a function of the erythrophagocytosis that occurs in early infection [36]. The implication is that, whilst the traits may be linked mechanistically, the parasite determinants that drive organomegaly/reticulocytosis and red blood cell numbers may be different.
In summary, we have demonstrated for the first time that forward genetics provides a powerful tool to map genes in the trypanosome genome that are responsible for causing complex phenotypes in the mouse host. We have defined two loci, one a major locus on Chromosome 3 (TbOrg1) that contributes to a significant amount of the variance observed in splenomegaly and hepatomagaly during mouse infections with T. brucei, and a second QTL on Chromosome 2 that contributes to splenomegaly, hepatomegaly and reticulocytosis in these infections (TbOrg2). Identification of the trypanosome genes that mediate the differential pathogenesis will provide a fundamental insight into the mechanisms by which the parasite causes disease in the mammalian host, and these findings may provide potential for developing therapeutic interventions to alleviate the disease associated with this pathogen. Figure S1 Inheritance model for TREU927 6 STIB247 F1 cross. In this simplified model the phenotype is treated as a binary trait determined by a single gene. (Spleen +++ = significantly enlarged spleen, spleen + = slightly enlarged spleen; A = dominant TREU927 allele, a = recessive TREU927 allele, x = STIB247 allele). The cross will allow the detection of traits encoded by a dominant heterozygous allele in TREU927. The genome of STIB247 is predominantly homozygous, with 94% of markers shown to be homozygous [33]. Therefore, the map is constructed based on segregation of alleles in the F1 progeny for loci heterozygous in the TREU 927 parent only. In order to detect co-dominant or recessive alleles, it will be necessary to generate a panel of F2 progeny. For segregation data see http://tinyurl.com/ trypmap.