A comprehensive genetic analysis of candidate genes regulating response to Trypanosoma congolense infection in mice

Elimination of Reference Mapping Bias Reveals Robust Immune Related Allele-Specific Expression in Crossbred Sheep

Pervasive allelic variation at both gene and single nucleotide level (SNV) between individuals is commonly associated with complex traits in humans and animals. Allele-specific expression (ASE) analysis, using RNA-Seq, can provide a detailed annotation of allelic imbalance and infer the existence of cis-acting transcriptional regulation. However, variant detection in RNA-Seq data is compromised by biased mapping of reads to the reference DNA sequence. In this manuscript, we describe an unbiased standardized computational pipeline for allele-specific expression analysis using RNA-Seq data, which we have adapted and developed using tools available under open license. The analysis pipeline we present is designed to minimize reference bias while providing accurate profiling of allele-specific expression across tissues and cell types. Using this methodology, we were able to profile pervasive allelic imbalance across tissues and cell types, at both the gene and SNV level, in Texel×Scottish Blackface sheep, using the sheep gene expression atlas data set. ASE profiles were pervasive in each sheep and across all tissue types investigated. However, ASE profiles shared across tissues were limited, and instead, they tended to be highly tissue-specific. These tissue-specific ASE profiles may underlie the expression of economically important traits and could be utilized as weighted SNVs, for example, to improve the accuracy of genomic selection in breeding programs for sheep. An additional benefit of the pipeline is that it does not require parental genotypes and can therefore be applied to other RNA-Seq data sets for livestock, including those available on the Functional Annotation of Animal Genomes (FAANG) data portal. This study is the first global characterization of moderate to extreme ASE in tissues and cell types from sheep. We have applied a robust methodology for ASE profiling to provide both a novel analysis of the multi-dimensional sheep gene expression atlas data set and a foundation for identifying the regulatory and expressed elements of the genome that are driving complex traits in livestock.

Association of APOL1 renal disease risk alleles with Trypanosoma brucei rhodesiense infection outcomes in the northern part of Malawi

Trypanosoma brucei (T.b.) rhodesiense is the cause of the acute form of human African trypanosomiasis (HAT) in eastern and southern African countries. There is some evidence that there is diversity in the disease progression of T.b. rhodesiense in different countries. HAT in Malawi is associated with a chronic haemo-lymphatic stage infection compared to other countries, such as Uganda, where the disease is acute with more marked neurological impairment. This has raised the question of the role of host genetic factors in infection outcomes. A candidate gene association study was conducted in the northern region of Malawi. This was a case-control study involving 202 subjects, 70 cases and 132 controls. All individuals were from one area; born in the area and had been exposed to the risk of infection since birth. Ninety-six markers were genotyped from 17 genes: IL10, IL8, IL4, HLA-G, TNFA, IL6, IFNG, MIF, APOL, HLA-A, IL1B, IL4R, IL12B, IL12R, HP, HPR, and CFH. There was a strong significant association with APOL1 G2 allele (p = 0.0000105, OR = 0.14, CI95 = [0.05-0.41], BONF = 0.00068) indicating that carriers of the G2 allele were protected against T.b. rhodesiense HAT. SNP rs2069845 in IL6 had raw p < 0.05, but did not remain significant after Bonferroni correction. There were no associations found with the other 15 candidate genes. Our finding confirms results from other studies that the G2 variant of APOL1 is associated with protection against T.b. rhodesiense HAT.

Animal African Trypanosomiasis: Time to Increase Focus on Clinically Relevant Parasite and Host Species

Animal African trypanosomiasis (AAT), caused by Trypanosoma congolense and Trypanosoma vivax, remains one of the most important livestock diseases in sub-Saharan Africa, particularly affecting cattle. Despite this, our detailed knowledge largely stems from the human pathogen Trypanosoma brucei and mouse experimental models. In the postgenomic era, the genotypic and phenotypic differences between the AAT-relevant species of parasite or host and their model organism counterparts are increasingly apparent. Here, we outline the timeliness and advantages of increasing the research focus on both the clinically relevant parasite and host species, given that improved tools and resources for both have been developed in recent years. We propose that this shift of emphasis will improve our ability to efficiently develop tools to combat AAT.

The within-host dynamics of African trypanosome infections

African trypanosomes are single-celled protozoan parasites that are capable of long-term survival while living extracellularly in the bloodstream and tissues of mammalian hosts. Prolonged infections are possible because trypanosomes undergo antigenic variation-the expression of a large repertoire of antigenically distinct surface coats, which allows the parasite population to evade antibody-mediated elimination. The mechanisms by which antigen genes become activated influence their order of expression, most likely by influencing the frequency of productive antigen switching, which in turn is likely to contribute to infection chronicity. Superimposed upon antigen switching as a contributor to trypanosome infection dynamics is the density-dependent production of cell-cycle arrested parasite transmission stages, which limit the infection while ensuring parasite spread to new hosts via the bite of blood-feeding tsetse flies. Neither antigen switching nor developmental progression to transmission stages is driven by the host. However, the host can contribute to the infection dynamic through the selection of distinct antigen types, the influence of genetic susceptibility or trypanotolerance and the potential influence of host-dependent effects on parasite virulence, development of transmission stages and pathogenicity. In a zoonotic infection cycle where trypanosomes circulate within a range of host animal populations, and in some cases humans, there is considerable scope for a complex interplay between parasite immune evasion, transmission potential and host factors to govern the profile and outcome of infection.

The QTL within the H2 Complex Involved in the Control of Tuberculosis Infection in Mice Is the Classical Class II H2-Ab1 Gene

The level of susceptibility to tuberculosis (TB) infection depends upon allelic variations in numerous interacting genes. In our mouse model system, the whole-genome quantitative trait loci (QTLs) scan revealed three QTLs involved in TB control on chromosomes 3, 9, and in the vicinity of the H2 complex on chromosome 17. For the present study, we have established a panel of new congenic, MHC-recombinant mouse strains bearing differential small segments of chromosome 17 transferred from the TB-susceptible I/St (H2^j) strain onto the genetic background of TB-resistant C57BL/6 (B6) mice (H2^b). This allowed narrowing the QTL interval to 17Ch: 33, 77–34, 34 Mb, containing 36 protein-encoding genes. Cloning and sequencing of the H2^j allelic variants of these genes demonstrated profound polymorphic variations compare to the H2^b haplotype. In two recombinant strains, B6.I-249.1.15.100 and B6.I-249.1.15.139, recombination breakpoints occurred in different sites of the H2-Aβ 1 gene (beta-chain of the Class II heterodimer H2-A), providing polymorphic variations in the domain β1 of the Aβ-chain. These variations were sufficient to produce different TB-relevant phenotypes: the more susceptible B6.I-249.1.15.100 strain demonstrated shorter survival time, more rapid body weight loss, higher mycobacterial loads in the lungs and more severe lung histopathology compared to the more resistant B6.I-249.1.15.139 strain. CD4⁺ T cells recognized mycobacterial antigens exclusively in the context of the H2-A Class II molecule, and the level of IFN-γ-producing CD4⁺ T cells in the lungs was significantly higher in the resistant strain. Thus, we directly demonstrated for the first time that the classical H2- Ab1 Class II gene is involved in TB control. Molecular modeling of the H2-A^j product predicts that amino acid (AA) substitutions in the Aβ-chain modify the motif of the peptide–MHC binding groove. Moreover, unique AA substitutions in both α- and β-chains of the H2-A^j molecule might affect its interactions with the T-cell receptor (TCR).

Many genes of the host regulate interactions with Mycobacterium tuberculosis and determine the level of susceptibility to, and severity of, tuberculosis (TB). Identification of these genes and their alleles is continuing and contributes new knowledge about the host-pathogen interactions. So far, forward genetic approaches (from phenotype to gene) have identified several chromosomal segments involved in genetic control of TB in mice (quantitative trait loci—QTL), but only one particular gene, Ipr1, has been identified. Here, we report the identification of a second TB-controlling gene. On the basis of a pair of mouse inbred strains with polar susceptibility to TB infection (susceptible I/St and more resistant C57BL/6) we established a panel of recombinant strains carrying small segments of Chromosome 17 from I/St on the genetic background of C57BL/6. A combination of genetic mapping, gene sequencing, TB phenotypes assessment and immunological approaches demonstrates that the H2-Ab1 gene encoding the beta-chain of the Class II heterodimer H2-A determines susceptibility to TB infection. The importance of allelic polymorphisms in Class II genes encoding antigen-presenting molecules in susceptibility to infection has been suspected. This is the first prove of this role obtained by the methods of classical forward genetics.

Tolerance to Trypanosomatids: A Threat, or a Key for Disease Elimination?

So far, research on trypanosomatid infections has been driven by 'disease by disease' approaches, leading to different concepts and control strategies. It is, however, increasingly clear that they share common features such as the ability to generate long-lasting asymptomatic infections in their mammalian hosts. Trypanotolerance, long integrated in animal African trypanosomiasis control, historically refers to the ability of cattle breeds to limit Trypanosoma infection and pathology, but has only recently been recognized in humans. Whilst trypanotolerance is absent from the vocabulary on leishmaniasis and Chagas disease, asymptomatic infections also occur. We review the concept of trypanotolerance across the trypanosomatids and discuss the importance of asymptomatic carriage in the current context of elimination.

The quality of methods reporting in parasitology experiments

There is a growing concern both inside and outside the scientific community over the lack of reproducibility of experiments. The depth and detail of reported methods are critical to the reproducibility of findings, but also for making it possible to compare and integrate data from different studies. In this study, we evaluated in detail the methods reporting in a comprehensive set of trypanosomiasis experiments that should enable valid reproduction, integration and comparison of research findings. We evaluated a subset of other parasitic (Leishmania, Toxoplasma, Plasmodium, Trichuris and Schistosoma) and non-parasitic (Mycobacterium) experimental infections in order to compare the quality of method reporting more generally. A systematic review using PubMed (2000-2012) of all publications describing gene expression in cells and animals infected with Trypanosoma spp was undertaken based on PRISMA guidelines; 23 papers were identified and included. We defined a checklist of essential parameters that should be reported and have scored the number of those parameters that are reported for each publication. Bibliometric parameters (impact factor, citations and h-index) were used to look for association between Journal and Author status and the quality of method reporting. Trichuriasis experiments achieved the highest scores and included the only paper to score 100% in all criteria. The mean of scores achieved by Trypanosoma articles through the checklist was 65.5% (range 32-90%). Bibliometric parameters were not correlated with the quality of method reporting (Spearman's rank correlation coefficient <-0.5; p>0.05). Our results indicate that the quality of methods reporting in experimental parasitology is a cause for concern and it has not improved over time, despite there being evidence that most of the assessed parameters do influence the results. We propose that our set of parameters be used as guidelines to improve the quality of the reporting of experimental infection models as a pre-requisite for integrating and comparing sets of data.

A parent-of-origin effect determines the susceptibility of a non-informative F1 population to Trypanosoma cruzi infection in vivo

The development of Chagas disease is determined by a complex interaction between the genetic traits of both the protozoan parasite, T. cruzi, and the infected host. This process is regulated by multiple genes that control different aspects of the host-parasite interaction. While determination of the relevant genes in humans is extremely difficult, it is feasible to use inbred mouse strains to determine the genes and loci responsible for host resistance to infection. In this study, we investigated the susceptibility of several inbred mouse strains to infection with the highly virulent Y strain of T. cruzi and found a considerable difference in susceptibility between A/J and C57BL/6 mice. We explored the differences between these two mouse strains and found that the A/J strain presented higher mortality, exacerbated and uncontrolled parasitemia and distinct histopathology in the target organs, which were associated with a higher parasite burden and more extensive tissue lesions. We then employed a genetic approach to assess the pattern of inheritance of the resistance phenotype in an F1 population and detected a strong parent-of-origin effect determining the susceptibility of the F1 male mice. This effect is unlikely to result from imprinted genes because the inheritance of this susceptibility was affected by the direction of the parental crossing. Collectively, our genetic approach of using the F1 population suggests that genes contained in the murine chromosome X contribute to the natural resistance against T. cruzi infection. Future linkage studies may reveal the locus and genes participating on the host resistance process reported herein.

Genetic control of resistance to Trypanosoma brucei brucei infection in mice

Background

Trypanosoma brucei brucei infects livestock, with severe effects in horses and dogs. Mouse strains differ greatly in susceptibility to this parasite. However, no genes controlling these differences were mapped.

Methods

We studied the genetic control of survival after T. b. brucei infection using recombinant congenic (RC) strains, which have a high mapping power. Each RC strain of BALB/c-c-STS/A (CcS/Dem) series contains a different random subset of 12.5% genes from the parental “donor” strain STS/A and 87.5% genes from the “background” strain BALB/c. Although BALB/c and STS/A mice are similarly susceptible to T. b. brucei, the RC strain CcS-11 is more susceptible than either of them. We analyzed genetics of survival in T. b. brucei-infected F₂ hybrids between BALB/c and CcS-11. CcS-11 strain carries STS-derived segments on eight chromosomes. They were genotyped in the F₂ hybrid mice and their linkage with survival was tested by analysis of variance.

Results

We mapped four Tbbr (Trypanosoma brucei brucei response) loci that influence survival after T. b. brucei infection. Tbbr1 (chromosome 3) and Tbbr2 (chromosome 12) have effects on survival independent of inter-genic interactions (main effects). Tbbr3 (chromosome 7) influences survival in interaction with Tbbr4 (chromosome 19). Tbbr2 is located on a segment 2.15 Mb short that contains only 26 genes.

Conclusion

This study presents the first identification of chromosomal loci controlling susceptibility to T. b. brucei infection. While mapping in F₂ hybrids of inbred strains usually has a precision of 40–80 Mb, in RC strains we mapped Tbbr2 to a 2.15 Mb segment containing only 26 genes, which will enable an effective search for the candidate gene. Definition of susceptibility genes will improve the understanding of pathways and genetic diversity underlying the disease and may result in new strategies to overcome the active subversion of the immune system by T. b. brucei.

Trypanosoma brucei are extracellular protozoa transmitted to mammalian host by the tsetse fly. They developed several mechanisms that subvert host's immune defenses. Therefore analysis of genes affecting host's resistance to infection can reveal critical aspects of host-parasite interactions. Trypanosoma brucei brucei infects many animal species including livestock, with particularly severe effects in horses and dogs. Mouse strains differ greatly in susceptibility to T. b. brucei. However, genes controlling susceptibility to this parasite have not been mapped. We analyzed the genetic control of survival after T. b. brucei infection using CcS/Dem recombinant congenic (RC) strains, each of which contains a different random set of 12.5% genes of their donor parental strain STS/A on the BALB/c genetic background. The RC strain CcS-11 is even more susceptible to parasites than BALB/c or STS/A. In F₂ hybrids between BALB/c and CcS-11 we detected and mapped four loci, Tbbr1-4 (Trypanosoma brucei brucei response 1–4), that control survival after T. b. brucei infection. Tbbr1 (chromosome 3) and Tbbr2 (chromosome 12) have independent effects, Tbbr3 (chromosome 7) and Tbbr4 (chromosome 19) were detected by their mutual inter-genic interaction. Tbbr2 was precision mapped to a segment of 2.15 Mb that contains 26 genes.

Multigenic control and sex bias in host susceptibility to spore-induced pulmonary anthrax in mice

Mechanisms underlying susceptibility to anthrax infection are unknown. Using a phylogenetically diverse panel of inbred mice and spores of Bacillus anthracis Ames, we investigated host susceptibility to pulmonary anthrax. Susceptibility profiles for survival time and organ pathogen load differed across strains, indicating distinct genetic controls. Tissue infection kinetics analysis showed greater systemic dissemination in susceptible DBA/2J (D) mice but a higher terminal bacterial load in resistant BALB/cJ (C) mice. Interestingly, the most resistant strains, C and C57BL/6J (B), demonstrated a sex bias for susceptibility. For example, BALB/cJ females had a significantly higher survival time and required 4-fold more spores for 100% mortality compared to BALB/cJ males. To identify genetic regions associated with differential susceptibility, survival time and extent of organ infection were assessed using mice derived from two susceptibility models: (i) BXD advanced recombinant inbred strains and (ii) F2 offspring generated from polar responding C and D strains. Genome-wide analysis of BXD strain survival identified linkage on chromosomes 5, 6, 9, 11, and 14. Quantitative trait locus (QTL) analysis of the C×DF2 population revealed a significant QTL (designated Rpai1 for resistance to pulmonary anthrax infection, locus 1) for survival time on chromosome 17 and also identified a chromosome 11 locus for lung pathogen burden. The striking difference between genome-wide linkage profiles for these two mouse models of anthrax susceptibility supports our hypothesis that these are multigenic traits. Our data provide the first evidence for a differential sex response to anthrax resistance and further highlight the unlikelihood of a single common genetic contribution for this response across strains.

Genetic and expression analysis of cattle identifies candidate genes in pathways responding to Trypanosoma congolense infection

African bovine trypanosomiasis caused by Trypanosoma sp., is a major constraint on cattle productivity in sub-Saharan Africa. Some African Bos taurus breeds are highly tolerant of infection, but the potentially more productive Bos indicus zebu breeds are much more susceptible. Zebu cattle are well adapted for plowing and haulage, and increasing their tolerance of trypanosomiasis could have a major impact on crop cultivation as well as dairy and beef production. We used three strategies to obtain short lists of candidate genes within QTL that were previously shown to regulate response to infection. We analyzed the transcriptomes of trypanotolerant N'Dama and susceptible Boran cattle after infection with Trypanosoma congolense. We sequenced EST libraries from these two breeds to identify polymorphisms that might underlie previously identified quantitative trait loci (QTL), and we assessed QTL regions and candidate loci for evidence of selective sweeps. The scan of the EST sequences identified a previously undescribed polymorphism in ARHGAP15 in the Bta2 trypanotolerance QTL. The polymorphism affects gene function in vitro and could contribute to the observed differences in expression of the MAPK pathway in vivo. The expression data showed that TLR and MAPK pathways responded to infection, and the former contained TICAM1, which is within a QTL on Bta7. Genetic analyses showed that selective sweeps had occurred at TICAM1 and ARHGAP15 loci in African taurine cattle, making them strong candidates for the genes underlying the QTL. Candidate QTL genes were identified in other QTL by their expression profile and the pathways in which they participate.