Transcriptome profiles of Trypanosoma brucei rhodesiense in Malawi reveal focus specific gene expression profiles associated with pathology.

BACKGROUND: Sleeping sickness caused by Trypanosoma brucei rhodesiense is a fatal disease and endemic in Southern and Eastern Africa. There is an urgent need to develop novel diagnostic and control tools to achieve elimination of rhodesiense sleeping sickness which might be achieved through a better understanding of trypanosome gene expression and genetics using endemic isolates. Here, we describe transcriptome profiles and population structure of endemic T. b. rhodesiense isolates in human blood in Malawi. METHODOLOGY: Blood samples of r-HAT cases from Nkhotakota and Rumphi foci were collected in PaxGene tubes for RNA extraction before initiation of r-HAT treatment. 100 million reads were obtained per sample, reads were initially mapped to the human genome reference GRCh38 using HiSat2 and then the unmapped reads were mapped against Trypanosoma brucei reference transcriptome (TriTrypDB54_TbruceiTREU927) using HiSat2. Differential gene expression analysis was done using the DeSeq2 package in R. SNP calling from reads that were mapped to the T. brucei genome was done using GATK in order to identify T.b. rhodesiense population structure. RESULTS: 24 samples were collected from r-HAT cases of which 8 were from Rumphi and 16 from Nkhotakota foci. The isolates from Nkhotakota were enriched with transcripts for cell cycle arrest and stumpy form markers, whereas isolates in Rumphi focus were enriched with transcripts for folate biosynthesis and antigenic variation pathways. These parasite focus-specific transcriptome profiles are consistent with the more virulent disease observed in Rumphi and a less symptomatic disease in Nkhotakota associated with the non-dividing stumpy form. Interestingly, the Malawi T.b. rhodesiense isolates expressed genes enriched for reduced cell proliferation compared to the Uganda T.b. rhodesiense isolates. PCA analysis using SNPs called from the RNAseq data showed that T. b. rhodesiense parasites from Nkhotakota are genetically distinct from those collected in Rumphi. CONCLUSION: Our results suggest that the differences in disease presentation in the two foci is mainly driven by genetic differences in the parasites in the two major endemic foci of Rumphi and Nkhotakota rather than differences in the environment or host response.

Impact of environmental factors on Biomphalaria pfeifferi vector capacity leading to human infection by Schistosoma mansoni in two regions of western Côte d'Ivoire.

ABSTRACT: BACKGROUND: Intestinal schistosomiasis remains a worrying health problem, particularly in western Côte d'Ivoire, despite control efforts. It is therefore necessary to understand all the factors involved in the development of the disease, including biotic and abiotic factors. The aim of this study was to examine the factors that could support the maintenance of the intermediate host and its vectorial capacity in western Côte d'Ivoire. METHODS: Data on river physicochemical, microbiological, and climatic parameters, the presence or absence of snails with Schistosoma mansoni, and human infections were collected between January 2020 and February 2021. Spearman rank correlation tests, Mann-Whitney, analysis of variance (ANOVA), and an appropriate model selection procedure were used to analyze the data. RESULTS: The overall prevalence of infected snails was 56.05%, with infection reaching 100% in some collection sites and localities. Of 26 sites examined, 25 contained thermophilic coliforms and 22 contained Escherichia coli. Biomphalaria pfeifferi was observed in environments with lower land surface temperature (LST) and higher relative air humidity (RAH), and B. pfeifferi infection predominated in more acidic environments. Thermal coliforms and E. coli preferred higher pH levels. Lower maximum LST (LST_Max) and higher RAH and minimum LST (LST_Min) were favorable to E. coli, and lower LST_Max favored coliforms. The presence of B. pfeifferi was positively influenced by water temperature (T °C), LST_Min, RAH, and precipitation (Pp) (P < 0.05) and negatively influenced by pH, total dissolved solids (TDS), electrical conductivity (EC), LST_Max, and mean land surface temperature (LST). The parameters pH, TDS, EC, LST_Min, LST, and Pp had a positive impact on snail infection, while LST_Max had a negative impact on infection. Only pH had a positive effect on coliform and E. coli abundance. Of the 701 people examined for human schistosomiasis, 73.13% were positive for the point-of-care circulating cathodic antigen (POC-CCA) test and 12.01% for the Kato-Katz (KK) test. A positive correlation was established between human infections and the abundance of Biomphalaria (r2 = 0.879, P = 0.04959). CONCLUSIONS: The results obtained reflect the environmental conditions that are conducive to the maintenance of S. mansoni infection in this part of the country. To combat this infection as effectively as possible, it will be necessary not only to redouble efforts but also to prioritize control according to the level of endemicity at the village level.

Phenome-wide analysis reveals epistatic associations between APOL1 variants and chronic kidney disease and multiple other disorders.

BACKGROUND: APOL1 variants G1 and G2 are common in populations with recent African ancestry. They are associated with protection from African sleeping sickness, however homozygosity or compound heterozygosity for these variants is associated with chronic kidney disease (CKD) and related conditions. What is not clear is the extent of associations with non-kidney-related disorders, and whether there are clusters of diseases associated with individual APOL1 genotypes. METHODS: Using a cohort of 7462 UK Biobank participants with recent African ancestry, we conducted a phenome-wide association study investigating associations between individual APOL1 genotypes and conditions identified by the International Classification of Disease phenotypes. FINDINGS: We identified 27 potential associations between individual APOL1 genotypes and a diverse range of conditions. G1/G2 compound heterozygotes were specifically associated with 26 of these conditions (all deleteriously), with an over-representation of infectious diseases (including hospitalisation and death resulting from COVID-19). The analysis also exposed complexities in the relationship between APOL1 and CKD that are not evident when risk variants are grouped together: G1 homozygosity, G2 homozygosity, and G1/G2 compound heterozygosity were each shown to be associated with distinct CKD phenotypes. The multi-locus nature of the G1/G2 genotype means that its associations would go undetected in a standard genome-wide association study. INTERPRETATION: Our findings have implications for understanding health risks and better-targeted detection, intervention, and therapeutic strategies, particularly in populations where APOL1 G1 and G2 are common such as in sub-Saharan Africa and its diaspora. FUNDING: This study was funded by the Wellcome Trust (209511/Z/17/Z) and H3Africa (H3A/18/004).

Differences in gene expression profiles in early and late stage rhodesiense HAT individuals in Malawi.

T. b. rhodesiense is the causative agent of Rhodesian human African trypanosomiasis (r-HAT) in Malawi. Clinical presentation of r-HAT in Malawi varies between foci and differs from East African HAT clinical phenotypes. The purpose of this study was to gain more insights into the transcriptomic profiles of patients with early stage 1 and late stage 2 HAT disease in Malawi. Whole blood from individuals infected with T. b. rhodesiense was used for RNA-Seq. Control samples were from healthy trypanosome negative individuals matched on sex, age range, and disease foci. Illumina sequence FASTQ reads were aligned to the GRCh38 release 84 human genome sequence using HiSat2 and differential analysis was done in R Studio using the DESeq2 package. XGR, ExpressAnalyst and InnateDB algorithms were used for functional annotation and gene enrichment analysis of significant differentially expressed genes. RNA-seq was done on 23 r-HAT case samples and 28 healthy controls with 7 controls excluded for downstream analysis as outliers. A total of 4519 genes were significant differentially expressed (p adjusted <0.05) in individuals with early stage 1 r-HAT disease (n = 12) and 1824 genes in individuals with late stage 2 r-HAT disease (n = 11) compared to controls. Enrichment of innate immune response genes through neutrophil activation was identified in individuals with both early and late stages of the disease. Additionally, lipid metabolism genes were enriched in late stage 2 disease. We further identified uniquely upregulated genes (log2 Fold Change 1.4-2.0) in stage 1 (ZNF354C) and stage 2 (TCN1 and MAGI3) blood. Our data add to the current understanding of the human transcriptome profiles during T. b. rhodesiense infection. We further identified biological pathways and transcripts enriched than were enriched during stage 1 and stage 2 r-HAT. Lastly, we have identified transcripts which should be explored in future research whether they have potential of being used in combination with other markers for staging or r-HAT.

Transcriptome analysis of peripheral blood of Schistosoma mansoni infected children from the Albert Nile region in Uganda reveals genes implicated in fibrosis pathology.

Over 290 million people are infected by schistosomes worldwide. Schistosomiasis control efforts focus on mass drug treatment with praziquantel (PZQ), a drug that kills the adult worm of all Schistosoma species. Nonetheless, re-infections have continued to be detected in endemic areas with individuals living in the same area presenting with varying infection intensities. Our objective was to characterize the transcriptome profiles in peripheral blood of children between 10-15 years with varying intensities of Schistosoma mansoni infection living along the Albert Nile in Uganda. RNA extracted from peripheral blood collected from 44 S. mansoni infected (34 high and 10 low by circulating anodic antigen [CAA] level) and 20 uninfected children was sequenced using Illumina NovaSeq S4 and the reads aligned to the GRCh38 human genome. Differential gene expression analysis was done using DESeq2. Principal component analysis revealed clustering of gene expression by gender when S. mansoni infected children were compared with uninfected children. In addition, we identified 14 DEGs between S. mansoni infected and uninfected individuals, 56 DEGs between children with high infection intensity and uninfected individuals, 33 DEGs between those with high infection intensity and low infection intensity and no DEGs between those with low infection and uninfected individuals. We also observed upregulation and downregulation of some DEGs that are associated with fibrosis and its regulation. These data suggest expression of fibrosis associated genes as well as genes that regulate fibrosis in S. mansoni infection. The relatively few significant DEGS observed in children with schistosomiasis suggests that chronic S. mansoni infection is a stealth infection that does not stimulate a strong immune response.

Variants of IL6, IL10, FCN2, RNASE3, IL12B and IL17B loci are associated with Schistosoma mansoni worm burden in the Albert Nile region of Uganda.

BACKGROUND: Individuals genetically susceptible to high schistosomiasis worm burden may contribute disproportionately to transmission and could be prioritized for control. Identifying genes involved may guide development of therapy. METHODOLOGY/PRINCIPAL FINDINGS: A cohort of 606 children aged 10-15 years were recruited in the Albert Nile region of Uganda and assessed for Schistosoma mansoni worm burden using the Up-Converting Particle Lateral Flow (UCP-LF) test detecting circulating anodic antigen (CAA), point-of-care Circulating Cathodic Antigen (POC-CCA) and Kato-Katz tests. Whole genome genotyping was conducted on 326 children comprising the top and bottom 25% of worm burden. Linear models were fitted to identify variants associated with worm burden in preselected candidate genes. Expression quantitative trait locus (eQTL) analysis was conducted for candidate genes with UCP-LF worm burden included as a covariate. Single Nucleotide Polymorphism loci associated with UCP-LF CAA included IL6 rs2066992 (OR = 0.43, p = 0.0006) and rs7793163 (OR = 2.0, p = 0.0007); IL21 SNP kgp513476 (OR 1.79, p = 0.0025) and IL17B SNP kgp708159 (OR = 0.35, p = 0.0028). A haplotype in the IL10 locus was associated with lower worm burden (OR = 0.53, p = 0.015) and overlapped SNPs rs1800896, rs1800871 and rs1800872. Significant haplotypes (p<0.05, overlapping significant SNP) associated with worm burden were observed in IL6 and the Th17 pathway IL12B and IL17B genes. There were significant eQTL in the IL6, IL5, IL21, IL25 and IFNG regions. CONCLUSIONS: Variants associated with S. mansoni worm burden were in IL6, FCN2, RNASE3, IL10, IL12B and IL17B gene loci. However only eQTL associations remained significant after Bonferroni correction. In summary, immune balance, pathogen recognition and Th17 pathways may play a role in modulating Schistosoma worm burden. Individuals carrying risk variants may be targeted first in allocation of control efforts to reduce the burden of schistosomiasis in the community.

Association between polymorphisms of IL4, IL13, IL10, STAT6 and IFNG genes, cytokines and immunoglobulin E levels with high burden of Schistosoma mansoni in children from schistosomiasis endemic areas of Cameroon.

Eliminating schistosomiasis as a public health problem by 2030 requires a better understanding of the disease transmission, especially the asymmetric distribution of worm burden in individuals living and sharing the same environment. It is in this light that this study was designed to identify human genetic determinants associated with high burden of S. mansoni and also with the plasma concentrations of IgE and four cytokines in children from two schistosomiasis endemic areas of Cameroon. In school-aged children of schistosomiasis endemic areas of Makenene and Nom-Kandi of Cameroon, S. mansoni infections and their infection intensities were evaluated in urine and stool samples using respectively the Point-of-care Circulating Cathodic Antigen test (POC-CCA) and the Kato Katz (KK) test. Thereafter, blood samples were collected in children harbouring high burden of schistosome infections as well as in their parents and siblings. DNA extracts and plasma were obtained from blood. Polymorphisms at 14 loci of five genes were assessed using PCR-restriction fragment length polymorphism and amplification-refractory mutation system. The ELISA test enabled to determine the plasma concentrations of IgE, IL-13, IL-10, IL-4 and IFN-γ. The prevalence of S. mansoni infections was significantly higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in Makenene (48.6% for POC-CCA and 7.9% for KK) compared to Nom-Kandi (31% for POC-CCA and 4.3% for KK). The infection intensities were also higher (P < 0.0001 for POC-CCA; P = 0.001 for KK) in children from Makenene than those from Nom-Kandi. The allele C of SNP rs3024974 of STAT6 was associated with an increased risk of bearing high burden of S. mansoni both in the additive (p = 0.009) and recessive model (p = 0.01) while the allele C of SNP rs1800871 of IL10 was protective (p = 0.0009) against high burden of S. mansoni. The alleles A of SNP rs2069739 of IL13 and G of SNP rs2243283 of IL4 were associated with an increased risk of having low plasma concentrations of IL-13 (P = 0.04) and IL-10 (P = 0.04), respectively. This study showed that host genetic polymorphisms may influence the outcome (high or low worm burden) of S. mansoni infections and also the plasma concentrations of some cytokines.

High prevalence of Schistosoma mansoni infection and stunting among school age children in communities along the Albert-Nile, Northern Uganda: A cross sectional study.

BACKGROUND: Knowing the prevalence of schistosomiasis is key to informing programmes to control and eliminate the disease as a public health problem. It is also important to understand the impact of infection on child growth and development in order to allocate appropriate resources and effort to the control of the disease. METHODS: We conducted a survey to estimate the prevalence of schistosomiasis among school aged children in villages along the Albert-Nile shore line in the district of Pakwach, North Western Uganda. A total of 914 children aged between 10-15 years were screened for Schistosoma mansoni using the POC-CCA and Kato Katz (KK) techniques. The infection intensities were assessed by POC-CCA and KK as well as CAA tests. The KK intensities were also correlated with POC-CCA and with CAA intensity. Anthropometric measurements were also taken and multivariate analysis was carried out to investigate their association with infection status. RESULTS: The prevalence of schistosomiasis using the POC-CCA diagnostic test was estimated at 85% (95% CI: 83-87), being highest amongst children living closer to the Albert-Nile shoreline. Visual scoring of the POC-CCA results was more sensitive than the Kato Katz test and was positively correlated with the quantified infection intensities by the CAA test. The majority of the children were underweight (BMI<18.5), and most notably, boys had significantly lower height for age (stunting) than girls in the same age range (p < 0.0001), but this was not directly associated with S. mansoni infection. CONCLUSION: High prevalence of S. mansoni infection in the region calls for more frequent mass drug administration with praziquantel. We observed high levels of stunting which was not associated with schistosomiasis. There is a need for improved nutrition among the children in the area.

SARS-CoV-2 seroprevalence at urban and rural sites in Kaduna State, Nigeria, during October/November 2021, immediately prior to detection of the Omicron variant.

BACKGROUND: Nigeria is Africa's most populated country. By November 2021 it had experienced three waves of SARS-CoV-2 infection. Peer-reviewed seroprevalence data assessing the proportion of the Nigerian population that have been infected were extremely limited. METHODS: We conducted a serosurvey in one urban site (n = 400) and one rural site (n = 402) in Kaduna State, Nigeria between 11 October 2021 and 8 November 2021. Z-tests were used to compare seroprevalence across age groups, locations and sexes. T tests were used to determine whether age or household size are associated with seropositivity. Associations between seropositivity and recent history of common Covid-19 symptoms were tested using logistic regression. RESULTS: SARS-CoV-2 antibodies were detected in 42.5% an 53.5% of participants at the urban and rural sites, respectively The overall age- and sex- stratified seroprevalence was 43.7% (42.2% for unvaccinated individuals). The data indicate an infection rate in Kaduna State ≥359-fold the rate derived from polymerase chain reaction-confirmed cases. In the urban site, seroprevalence among females and participants aged <20 was lower than other groups. Reporting loss of sense of taste and/or smell was strongly associated with seropositive status. Associations with seropositivity were also found for the reporting of dry cough, fever, headache, nausea and sore throat. CONCLUSIONS: This study provides baseline SARS-CoV-2 seroprevalence in Kaduna State, Nigeria, immediately prior to the spread of the Omicron variant. It indicates that in October/November 2021, approximately 56% of the population did not have detectable antibodies, and population subgroups with particularly low seroprevalence remain. It highlights limitations in using PCR-confirmed cases to estimate infection rates. The data will inform public health strategies in Nigeria and other sub-Saharan African countries with limited SARS-CoV-2 seroprevalence data.

Field assessment in Cameroon of a reader of POC-CCA lateral flow strips for the quantification of Schistosoma mansoni circulating cathodic antigen in urine.

BACKGROUND: Determining Schistosoma mansoni infection rate and intensity is challenging due to the low sensitivity of the Kato-Katz (KK) test that underestimates the true disease prevalence. Circulating cathodic antigen (CCA) excreted in urine is constantly produced by adult worms and has been used as the basis of a simple, non-invasive point of care test (POC-CCA) for Schistosoma mansoni infections. Although the abundance of CCA in urine is proportional to worm burden, the POC-CCA test is marketed as a qualitative test, making it difficult to investigate the wide range of infection intensities. This study was designed to compare the prevalence and intensity of S. mansoni by KK and POC-CCA and quantify, on fresh and frozen (<-20°C) urine samples, CCA using the visual scores and the ESEquant LR3 reader. METHODOLOGY: Stool and urine samples were collected from 759 school-aged children. The prevalence and intensity of S. mansoni were determined using KK and POC-CCA. The degree of the positivity of POC-CCA was estimated by quantifying CCA on fresh and frozen urine samples using visual scores and strip reader. The prevalence, the infection intensity as well the relative amounts of CCA were compared. RESULTS: The S. mansoni infection rates inferred from POC-CCA and KK were 40.7% and 9.4% respectively. Good correlations were observed between infection intensities recorded by; i) the reader and visual scoring system on fresh (Rho = 0.89) and frozen samples (Rho = 0.97), ii) the reader on fresh urine samples and KK (epg) (Rho = 0.44). Nevertheless, 238 POC-CCA positive children were negative for KK, and sixteen of them had high levels of CCA. The correlation between results from the reader on fresh and frozen samples was good (Rho = 0.85). On frozen samples, CCA was not detected in 55 samples that were positive in fresh urine samples. CONCLUSION: This study confirmed the low sensitivity of KK and the high capacity of POC-CCA to provide reliable data on the prevalence and intensity of S. mansoni infections. The lateral flow reader enabled accurate quantification of CCA under field conditions on fresh and frozen urine samples with less time and effort than KK.

Unmapped exome reads implicate a role for Anelloviridae in childhood HIV-1 long-term non-progression.

Human immunodeficiency virus (HIV) infection remains a significant public health burden globally. The role of viral co-infection in the rate of progression of HIV infection has been suggested but not empirically tested, particularly among children. We extracted and classified 42 viral species from whole-exome sequencing (WES) data of 813 HIV-infected children in Botswana and Uganda categorised as either long-term non-progressors (LTNPs) or rapid progressors (RPs). The Ugandan participants had a higher viral community diversity index compared to Batswana (p = 4.6 × 10-13), and viral sequences were more frequently detected among LTNPs than RPs (24% vs 16%; p = 0.008; OR, 1.9; 95% CI, 1.6-2.3), with Anelloviridae showing strong association with LTNP status (p = 3 × 10-4; q = 0.004, OR, 3.99; 95% CI, 1.74-10.25). This trend was still evident when stratified by country, sex, and sequencing platform, and after a logistic regression analysis adjusting for age, sex, country, and the sequencing platform (p = 0.02; q = 0.03; OR, 7.3; 95% CI, 1.6-40.5). Torque teno virus (TTV), which made up 95% of the Anelloviridae reads, has been associated with reduced immune activation. We identify an association between viral co-infection and prolonged AIDs-free survival status that may have utility as a biomarker of LTNP and could provide mechanistic insights to HIV progression in children, demonstrating the added value of interrogating off-target WES reads in cohort studies.

The Genetics of Human Schistosomiasis Infection Intensity and Liver Disease: A Review.

Schistosomiasis remains the fourth most prevalent parasitic disease affecting over 200 million people worldwide. Control efforts have focussed on the disruption of the life cycle targeting the parasite, vector and human host. Parasite burdens are highly skewed, and the majority of eggs are shed into the environment by a minority of the infected population. Most morbidity results from hepatic fibrosis leading to portal hypertension and is not well-correlated with worm burden. Genetics as well as environmental factors may play a role in these skewed distributions and understanding the genetic risk factors for intensity of infection and morbidity may help improve control measures. In this review, we focus on how genetic factors may influence parasite load, hepatic fibrosis and portal hypertension. We found 28 studies on the genetics of human infection and 20 studies on the genetics of pathology in humans. S. mansoni and S. haematobium infection intensity have been showed to be controlled by a major quantitative trait locus SM1, on chromosome 5q31-q33 containing several genes involved in the Th2 immune response, and three other loci of smaller effect on chromosomes 1, 6, and 7. The most common pathology associated with schistosomiasis is hepatic and portal vein fibroses and the SM2 quantitative trait locus on chromosome six has been linked to intensity of fibrosis. Although there has been an emphasis on Th2 cytokines in candidate gene studies, we found that four of the five QTL regions contain Th17 pathway genes that have been included in schistosomiasis studies: IL17B and IL12B in SM1, IL17A and IL17F in 6p21-q2, IL6R in 1p21-q23 and IL22RA2 in SM2. The Th17 pathway is known to be involved in response to schistosome infection and hepatic fibrosis but variants in this pathway have not been tested for any effect on the regulation of these phenotypes. These should be priorities for future studies.

Candidate gene family-based and case-control studies of susceptibility to high Schistosoma mansoni worm burden in African children: a protocol.

Background: Approximately 25% of the risk of Schistosoma mansoni is associated with host genetic variation. We will test 24 candidate genes, mainly in the T h2 and T h17 pathways, for association with S. mansoni infection intensity in four African countries, using family based and case-control approaches. Methods: Children aged 5-15 years will be recruited in S. mansoni endemic areas of Ivory Coast, Cameroon, Uganda and the Democratic Republic of Congo (DRC). We will use family based (study 1) and case-control (study 2) designs. Study 1 will take place in Ivory Coast, Cameroon, Uganda and the DRC. We aim to recruit 100 high worm burden families from each country except Uganda, where a previous study recruited at least 40 families. For phenotyping, cases will be defined as the 20% of children in each community with heaviest worm burdens as measured by the circulating cathodic antigen (CCA) assay. Study 2 will take place in Uganda. We will recruit 500 children in a highly endemic community. For phenotyping, cases will be defined as the 20% of children with heaviest worm burdens as measured by the CAA assay, while controls will be the 20% of infected children with the lightest worm burdens. Deoxyribonucleic acid (DNA) will be genotyped on the Illumina H3Africa SNP (single nucleotide polymorphisms) chip and genotypes will be converted to sets of haplotypes that span the gene region for analysis. We have selected 24 genes for genotyping that are mainly in the Th2 and Th17 pathways and that have variants that have been demonstrated to be or could be associated with Schistosoma infection intensity.   Analysis: In the family-based design, we will identify SNP haplotypes disproportionately transmitted to children with high worm burden. Case-control analysis will detect overrepresentation of haplotypes in extreme phenotypes with correction for relatedness by using whole genome principal components.

SNPs in IL4 and IFNG show no protective associations with human African trypanosomiasis in the Democratic Republic of the Congo: a case-control study.

Background: Human African trypanosomiasis (HAT) is a protozoal disease transmitted by tsetse flies. Infection with trypanosomes can lead directly to active HAT or latent infection with no detectable parasites, which may progress to active HAT or to spontaneous self-cure. Genetic variation could explain these differences in the outcome of infection. To test this hypothesis, polymorphisms in 17 candidate genes were tested ( APOL1 [ G1 and G2], CFH, HLA-A, HPR, HP, IL1B, IL12B, IL12RB1, IL10, IL4R, MIF, TNFA , IL6, IL4, IL8, IFNG, and HLA-G). Methods: Samples were collected in Democratic Republic of the Congo. 233 samples were genotyped: 100 active HAT cases, 33 from subjects with latent infections and 100 negative controls. Commercial service providers genotyped polymorphisms at 96 single nucleotide polymorphisms (SNPs) on 17 genes. Data were analyzed using Plink V1.9 software and R. Loci, with suggestive associations (uncorrected p < 0.05) validated using an additional 594 individuals, including 164 cases and 430 controls. Results: After quality control, 87 SNPs remained in the analysis. Two SNPs in IL4 and two in IFNG were suggestively associated (uncorrected p<0.05) with a differential risk of developing a Trypanosoma brucei gambiense infection in the Congolese population. The IFNG minor allele (rs2430561, rs2069718) SNPs were protective in comparison between latent infections and controls. Carriers of the rs2243258_T and rs2243279_A alleles of IL4 and the rs2069728_T allele of IFNG had a reduced risk of developing illness or latent infection, respectively. None of these associations were significant after Bonferroni correction for multiple testing. A validation study using more samples was run to determine if the absence of significant association was due to lack of power. Conclusions: This study showed no evidence of an association of HAT with IL4 and IFNG SNPs or with APOL1 G1 and G2 alleles, which have been found to be protective in other studies.

High Levels of Genetic Diversity within Nilo-Saharan Populations: Implications for Human Adaptation.

Africa contains more human genetic variation than any other continent, but the majority of the population-scale analyses of the African peoples have focused on just two of the four major linguistic groups, the Niger-Congo and Afro-Asiatic, leaving the Nilo-Saharan and Khoisan populations under-represented. In order to assess genetic variation and signatures of selection within a Nilo-Saharan population and between the Nilo-Saharan and Niger-Congo and Afro-Asiatic, we sequenced 50 genomes from the Nilo-Saharan Lugbara population of North-West Uganda and 250 genomes from 6 previously unsequenced Niger-Congo populations. We compared these data to data from a further 16 Eurasian and African populations including the Gumuz, another putative Nilo-Saharan population from Ethiopia. Of the 21 million variants identified in the Nilo-Saharan population, 3.57 million (17%) were not represented in dbSNP and included predicted non-synonymous mutations with possible phenotypic effects. We found greater genetic differentiation between the Nilo-Saharan Lugbara and Gumuz populations than between any two Afro-Asiatic or Niger-Congo populations. F3 tests showed that Gumuz contributed a genetic component to most Niger-Congo B populations whereas Lugabara did not. We scanned the genomes of the Lugbara for evidence of selective sweeps. We found selective sweeps at four loci (SLC24A5, SNX13, TYRP1, and UVRAG) associated with skin pigmentation, three of which already have been reported to be under selection. These selective sweeps point toward adaptations to the intense UV radiation of the Sahel.

Copy number variation in human genomes from three major ethno-linguistic groups in Africa.

BACKGROUND: Copy number variation is an important class of genomic variation that has been reported in 75% of the human genome. However, it is underreported in African populations. Copy number variants (CNVs) could have important impacts on disease susceptibility and environmental adaptation. To describe CNVs and their possible impacts in Africans, we sequenced genomes of 232 individuals from three major African ethno-linguistic groups: (1) Niger Congo A from Guinea and Côte d'Ivoire, (2) Niger Congo B from Uganda and the Democratic Republic of Congo and (3) Nilo-Saharans from Uganda. We used GenomeSTRiP and cn.MOPS to identify copy number variant regions (CNVRs). RESULTS: We detected 7608 CNVRs, of which 2172 were only deletions, 2384 were only insertions and 3052 had both. We detected 224 previously un-described CNVRs. The majority of novel CNVRs were present at low frequency and were not shared between populations. We tested for evidence of selection associated with CNVs and also for population structure. Signatures of selection identified previously, using SNPs from the same populations, were overrepresented in CNVRs. When CNVs were tagged with SNP haplotypes to identify SNPs that could predict the presence of CNVs, we identified haplotypes tagging 3096 CNVRs, 372 CNVRs had SNPs with evidence of selection (iHS > 3) and 222 CNVRs had both. This was more than expected (p < 0.0001) and included loci where CNVs have previously been associated with HIV, Rhesus D and preeclampsia. When integrated with 1000 Genomes CNV data, we replicated their observation of population stratification by continent but no clustering by populations within Africa, despite inclusion of Nilo-Saharans and Niger-Congo populations within our dataset. CONCLUSIONS: Novel CNVRs in the current study increase representation of African diversity in the database of genomic variants. Over-representation of CNVRs in SNP signatures of selection and an excess of SNPs that both tag CNVs and are subject to selection show that CNVs may be the actual targets of selection at some loci. However, unlike SNPs, CNVs alone do not resolve African ethno-linguistic groups. Tag haplotypes for CNVs identified may be useful in predicting African CNVs in future studies where only SNP data is available.

Variant antigen diversity in Trypanosoma vivax is not driven by recombination.

African trypanosomes (Trypanosoma) are vector-borne haemoparasites that survive in the vertebrate bloodstream through antigenic variation of their Variant Surface Glycoprotein (VSG). Recombination, or rather segmented gene conversion, is fundamental in Trypanosoma brucei for both VSG gene switching and for generating antigenic diversity during infections. Trypanosoma vivax is a related, livestock pathogen whose VSG lack structures that facilitate gene conversion in T. brucei and mechanisms underlying its antigenic diversity are poorly understood. Here we show that species-wide VSG repertoire is broadly conserved across diverse T. vivax clinical strains and has limited antigenic repertoire. We use variant antigen profiling, coalescent approaches and experimental infections to show that recombination plays little role in diversifying T. vivax VSG sequences. These results have immediate consequences for both the current mechanistic model of antigenic variation in African trypanosomes and species differences in virulence and transmission, requiring reconsideration of the wider epidemiology of animal African trypanosomiasis.

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.

Macrophage migrating inhibitory factor expression is associated with Trypanosoma brucei gambiense infection and is controlled by trans-acting expression quantitative trait loci in the Guinean population.

Infection by Trypanosoma brucei gambiense is characterized by a wide array of clinical outcomes, ranging from asymptomatic to acute disease and even spontaneous cure. In this study, we investigated the association between macrophage migrating inhibitory factor (MIF), an important pro-inflammatory cytokine that plays a central role in both innate and acquired immunity, and disease outcome during T. b. gambiense infection. A comparative expression analysis of patients, individuals with latent infection and controls found that MIF had significantly higher expression in patients (n = 141; 1.25 ±â€¯0.07; p < .0001) and latent infections (n = 25; 1.23 ±â€¯0.13; p = .0005) relative to controls (n = 46; 0.94 ±â€¯0.11). Furthermore, expression decreased significantly after treatment (patients before treatment n = 33; 1.40 ±â€¯0.18 versus patients after treatment n = 33; 0.99 ±â€¯0.10, p = .0001). We conducted a genome wide eQTL analysis on 29 controls, 128 cases and 15 latently infected individuals for whom expression and genotype data were both available. Four loci, including one containing the chemokine CXCL13, were found to associate with MIF expression. Genes at these loci are candidate regulators of increased expression of MIF after infection. Our study is the first data demonstrating that MIF expression is elevated in T. b. gambiense-infected human hosts but does not appear to contribute to pathology.

Association between IL1 gene polymorphism and human African trypanosomiasis in populations of sleeping sickness foci of southern Cameroon.

BACKGROUND: Human African Trypanosomiasis (HAT) is a neglected tropical disease caused by infections due to Trypanosoma brucei subspecies. In addition to the well-established environmental and behavioural risks of becoming infected, there is evidence for a genetic component to the response to trypanosome infection. We undertook a candidate gene case-control study to investigate genetic associations further. METHODOLOGY: We genotyped one polymorphism in each of seven genes (IL1A, IL1RN, IL4RN, IL6, HP, HPR, and HLA-G) in 73 cases and 250 controls collected from 19 ethno-linguistic subgroups stratified into three major ethno-linguistic groups, 2 pooled ethno-linguistic groups and 11 ethno-linguistic subgroups from three Cameroonian HAT foci. The seven polymorphic loci tested consisted of three SNPs, three variable numbers of tandem repeat (VNTR) and one INDEL. RESULTS: We found that the genotype (TT) and minor allele (T) of IL1A gene as well as the genotype 1A3A of IL1RN were associated with an increased risk of getting Trypanosoma brucei gambiense and develop HAT when all data were analysed together and also when stratified by the three major ethno-linguistic groups, 2 pooled ethno-linguistic subgroups and 11 ethno-linguistic subgroups. CONCLUSION: This study revealed that one SNP rs1800794 of IL1A and one VNTR rs2234663 of IL1RN were associated with the increased risk to be infected by Trypanosoma brucei gambiense and develop sleeping sickness in southern Cameroon. The minor allele T and the genotype TT of SNP rs1800794 in IL1A as well as the genotype 1A3A of IL1RN rs2234663 VNTR seem to increase the risk of getting Trypanosoma brucei gambiense infections and develop sleeping sickness in southern Cameroon.