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.

Associations between Aminoquinoline Resistance Genotypes and Clinical Presentations of Plasmodium falciparum Infection in Uganda.

Mutations that mediate resistance of Plasmodium falciparum to aminoquinoline antimalarials are selected by prior drug use and may alter parasite fitness, but associations with clinical presentations are uncertain. We evaluated genotypes in samples from a case-control study of determinants of severe malaria in Ugandan children 4 months to 10 years of age. We studied 274 cases with severe malaria matched by age and geography to 275 uncomplicated malaria controls and 179 asymptomatic parasitemic controls. The overall prevalence of mutations of interest (considering mixed results as mutant) was 67.0% for PfCRT K76T, 8.5% for PfMDR1 N86Y, 71.5% for PfMDR1 Y184F, and 14.7% for PfMDR1 D1246Y. Compared to asymptomatic controls, the odds of mutant PfCRT 76T were lower for uncomplicated (odds ratio, 0.42 [95% confidence interval, 0.24 to 0.72]; P < 0.001) or severe (0.56 [0.32 to 0.97]; P = 0.031) malaria; the odds of mutant PfMDR1 86Y were lower for uncomplicated (0.33 [0.16 to 0.65]; P < 0.001) or severe (0.21 [0.09 to 0.45]; P < 0.001) malaria; and the odds of mutant PfMDR1 1246Y were higher for uncomplicated (1.83 [0.90 to 3.98]; P = 0.076) or severe (2.06 [1.01 to 4.55]; P = 0.033) malaria. The odds of mutant PfMDR1 184F were lower in severe than asymptomatic (0.59 [0.37 to 0.92]; P = 0.016) or uncomplicated (0.61 [0.41 to 0.90]; P = 0.009) malaria. Overall, the PfCRT 76T and PfMDR1 86Y mutations were associated with decreased risk of symptomatic malaria, PfMDR1 1246Y was associated with increased risk of symptomatic malaria, and PfMDR1 184F was associated with decreased risk of severe malaria. These results offer insights into parasite genotypes in children with different presentations, although the basis for the identified associations is likely complex.

Identification of Novel Potential Inhibitors of Pteridine Reductase 1 in Trypanosoma brucei via Computational Structure-Based Approaches and in Vitro Inhibition Assays.

Pteridine reductase 1 (PTR1) is a trypanosomatid multifunctional enzyme that provides a mechanism for escape of dihydrofolate reductase (DHFR) inhibition. This is because PTR1 can reduce pterins and folates. Trypanosomes require folates and pterins for survival and are unable to synthesize them de novo. Currently there are no anti-folate based Human African Trypanosomiasis (HAT) chemotherapeutics in use. Thus, successful dual inhibition of Trypanosoma brucei dihydrofolate reductase (TbDHFR) and Trypanosoma brucei pteridine reductase 1 (TbPTR1) has implications in the exploitation of anti-folates. We carried out molecular docking of a ligand library of 5742 compounds against TbPTR1 and identified 18 compounds showing promising binding modes. The protein-ligand complexes were subjected to molecular dynamics to characterize their molecular interactions and energetics, followed by in vitro testing. In this study, we identified five compounds which showed low micromolar Trypanosome growth inhibition in in vitro experiments that might be acting by inhibition of TbPTR1. Compounds RUBi004, RUBi007, RUBi014, and RUBi018 displayed moderate to strong antagonism (mutual reduction in potency) when used in combination with the known TbDHFR inhibitor, WR99210. This gave an indication that the compounds might inhibit both TbPTR1 and TbDHFR. RUBi016 showed an additive effect in the isobologram assay. Overall, our results provide a basis for scaffold optimization for further studies in the development of HAT anti-folates.

Assessing computational genomics skills: Our experience in the H3ABioNet African bioinformatics network.

The H3ABioNet pan-African bioinformatics network, which is funded to support the Human Heredity and Health in Africa (H3Africa) program, has developed node-assessment exercises to gauge the ability of its participating research and service groups to analyze typical genome-wide datasets being generated by H3Africa research groups. We describe a framework for the assessment of computational genomics analysis skills, which includes standard operating procedures, training and test datasets, and a process for administering the exercise. We present the experiences of 3 research groups that have taken the exercise and the impact on their ability to manage complex projects. Finally, we discuss the reasons why many H3ABioNet nodes have declined so far to participate and potential strategies to encourage them to do so.

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.

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.

Anti-Tick Vaccines: Current Advances and Future Prospects.

Ticks are increasingly a global public health and veterinary concern. They transmit numerous pathogens that are of veterinary and public health importance. Acaricides, livestock breeding for tick resistance, tick handpicking, pasture spelling, and anti-tick vaccines (ATVs) are in use for the control of ticks and tick-borne diseases (TTBDs); acaricides and ATVs being the most and least used TTBD control methods respectively. The overuse and misuse of acaricides has inadvertently selected for tick strains that are resistant to acaricides. Furthermore, vaccines are rare and not commercially available in sub-Saharan Africa (SSA). It doesn't help that many of the other methods are labor-intensive and found impractical especially for larger farm operations. The success of TTBD control is therefore dependent on integrating all the currently available methods. Vaccines have been shown to be cheap and effective. However, their large-scale deployment for TTBD control in SSA is hindered by commercial unavailability of efficacious anti-tick vaccines against sub-Saharan African tick strains. Thanks to advances in genomics, transcriptomics, and proteomics technologies, many promising anti-tick vaccine antigens (ATVA) have been identified. However, few of them have been investigated for their potential as ATV candidates. Reverse vaccinology (RV) can be leveraged to accelerate ATV discovery. It is cheap and shortens the lead time from ATVA discovery to vaccine production. This chapter provides a brief overview of recent advances in ATV development, ATVs, ATV effector mechanisms, and anti-tick RV. Additionally, it provides a detailed outline of vaccine antigen selection and analysis using computational methods.

Genome Sequences of Bacteriophages UPEC01, UPEC03, UPEC06, and UPEC07 Infecting Avian Pathogenic Escherichia coli.

Here, we present the genome sequences of four bacteriophages that infect avian pathogenic Escherichia coli. The phages were isolated from raw sewage in Kampala, Uganda. The genome sizes of the phages ranged between 143,140 bp and 178,307 bp, with an average G+C content of 41.25%.

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.

Virtual screening and in vitro validation identifies the first reported inhibitors of Salmonella enterica HPPK.

HPPK, which directly precedes DHPS in the folate biosynthetic pathway, is a promising but chronically under-exploited anti-microbial target. Here we report the identification of new S. enterica HPPK inhibitors, offering potential for new resistance circumventing S. enterica therapies as well as avenues for diversifying the current HPPK inhibitor space.

Identification of Ixodid Tick-Specific Aquaporin-1 Potential Anti-tick Vaccine Epitopes: An in-silico Analysis.

Ticks are arthropod vectors of pathogens of both Veterinary and Public health importance. Acaricide application, which is currently the mainstay of tick control, is hampered by high cost, the need for regular application and a selection of multi-acaricide resistant tick populations. In light of this, future tick control approaches are poised to rely on the integration of rational acaricide application and other methods, such as vaccination. To contribute to systematic research-guided efforts to produce anti-tick vaccines, we carried out an in-silico analysis of tick aquaporin-1 (AQP1) protein in order to identify tick-specific AQP1 peptide motifs that can be used in future peptide anti-tick vaccine development. We carried out multiple sequence alignment (MSA), motif analysis, homology modeling, and structural analysis to identify tick-specific AQP1 peptide motifs. BepiPred, Chou and Fasman-Turn, Karplus and Schulz Flexibility, and Parker-Hydrophilicity prediction models were used to predict these motifs' potential to induce B cell mediated immune responses. The tick AQP1 (GenBankID: QDO67142.1) protein was largely similar to the bovine AQP1 (PDB:1J4N) (23 % sequence similarity; Structural superimposition of the homology model and 14JN homotetramers gave an RMSD = 3.269 while superimposition of a single chain gave an RMSD = 1.475). Tick and bovine AQP1 transmembrane domains were largely similar while their extracellular and cytoplasmic domain loops showed variation. Two tick-specific AQP1 peptide motifs; M7 (residues 106-125, p = 5.4e-25), and M8 (residues 85-104, p = 3.3e-24) were identified. These two motifs are located on the extra-cellular AQP1 domain. Motifs; M7 and M8 showed the highest Parker-Hydrophilicity prediction immunogenicity scores of 1.784 and 1.536, respectively. These two motifs can be a good starting point for the development of potential tick AQP1 peptide-based anti-tick vaccines. Further analyses such as molecular dynamics, in vitro assays, and in vivo immunization assays are required to validate these findings.

No evidence for association between APOL1 kidney disease risk alleles and Human African Trypanosomiasis in two Ugandan populations.

BACKGROUND: Human African trypanosomiasis (HAT) manifests as an acute form caused by Trypanosoma brucei rhodesiense (Tbr) and a chronic form caused by Trypanosoma brucei gambiense (Tbg). Previous studies have suggested a host genetic role in infection outcomes, particularly for APOL1. We have undertaken candidate gene association studies (CGAS) in a Ugandan Tbr and a Tbg HAT endemic area, to determine whether polymorphisms in IL10, IL8, IL4, HLAG, TNFA, TNX4LB, IL6, IFNG, MIF, APOL1, HLAA, IL1B, IL4R, IL12B, IL12R, HP, HPR, and CFH have a role in HAT. METHODOLOGY AND RESULTS: We included 238 and 202 participants from the Busoga Tbr and Northwest Uganda Tbg endemic areas respectively. Single Nucleotide Polymorphism (SNP) genotype data were analysed in the CGAS. The study was powered to find odds ratios > 2 but association testing of the SNPs with HAT yielded no positive associations i.e. none significant after correction for multiple testing. However there was strong evidence for no association with Tbr HAT and APOL1 G2 of the size previously reported in the Kabermaido district of Uganda. CONCLUSIONS/SIGNIFICANCE: A recent study in the Soroti and Kaberamaido focus in Central Uganda found that the APOL1 G2 allele was strongly associated with protection against Tbr HAT (odds ratio = 0.2, 95% CI: 0.07 to 0.48, p = 0.0001). However, in our study no effect of G2 on Tbr HAT was found, despite being well powered to find a similar sized effect (OR = 0.9281, 95% CI: 0.482 to 1.788, p = 0.8035). It is possible that the G2 allele is protective from Tbr in the Soroti/Kabermaido focus but not in the Iganga district of Busoga, which differ in ethnicity and infection history. Mechanisms underlying HAT infection outcome and virulence are complex and might differ between populations, and likely involve several host, parasite or even environmental factors.