High-depth African genomes inform human migration and health

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed1. Here we performed whole-genome sequencing analyses of 426 individuals—comprising 50 ethnolinguistic groups, including previously unsampled populations—to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon—but in other genes, variants denoted as ‘likely pathogenic’ in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health

High-depth African genomes inform human migration and health

The African continent is regarded as the cradle of modern humans and African genomes contain more genetic variation than those from any other continent, yet only a fraction of the genetic diversity among African individuals has been surveyed1. Here we performed whole-genome sequencing analyses of 426 individuals—comprising 50 ethnolinguistic groups, including previously unsampled populations—to explore the breadth of genomic diversity across Africa. We uncovered more than 3 million previously undescribed variants, most of which were found among individuals from newly sampled ethnolinguistic groups, as well as 62 previously unreported loci that are under strong selection, which were predominantly found in genes that are involved in viral immunity, DNA repair and metabolism. We observed complex patterns of ancestral admixture and putative-damaging and novel variation, both within and between populations, alongside evidence that Zambia was a likely intermediate site along the routes of expansion of Bantu-speaking populations. Pathogenic variants in genes that are currently characterized as medically relevant were uncommon—but in other genes, variants denoted as ‘likely pathogenic’ in the ClinVar database were commonly observed. Collectively, these findings refine our current understanding of continental migration, identify gene flow and the response to human disease as strong drivers of genome-level population variation, and underscore the scientific imperative for a broader characterization of the genomic diversity of African individuals to understand human ancestry and improve health

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

INTRODUCTION Investment in Africa over the past year with regard to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing has led to a massive increase in the number of sequences, which, to date, exceeds 100,000 sequences generated to track the pandemic on the continent. These sequences have profoundly affected how public health officials in Africa have navigated the COVID-19 pandemic. RATIONALE We demonstrate how the first 100,000 SARS-CoV-2 sequences from Africa have helped monitor the epidemic on the continent, how genomic surveillance expanded over the course of the pandemic, and how we adapted our sequencing methods to deal with an evolving virus. Finally, we also examine how viral lineages have spread across the continent in a phylogeographic framework to gain insights into the underlying temporal and spatial transmission dynamics for several variants of concern (VOCs). RESULTS Our results indicate that the number of countries in Africa that can sequence the virus within their own borders is growing and that this is coupled with a shorter turnaround time from the time of sampling to sequence submission. Ongoing evolution necessitated the continual updating of primer sets, and, as a result, eight primer sets were designed in tandem with viral evolution and used to ensure effective sequencing of the virus. The pandemic unfolded through multiple waves of infection that were each driven by distinct genetic lineages, with B.1-like ancestral strains associated with the first pandemic wave of infections in 2020. Successive waves on the continent were fueled by different VOCs, with Alpha and Beta cocirculating in distinct spatial patterns during the second wave and Delta and Omicron affecting the whole continent during the third and fourth waves, respectively. Phylogeographic reconstruction points toward distinct differences in viral importation and exportation patterns associated with the Alpha, Beta, Delta, and Omicron variants and subvariants, when considering both Africa versus the rest of the world and viral dissemination within the continent. Our epidemiological and phylogenetic inferences therefore underscore the heterogeneous nature of the pandemic on the continent and highlight key insights and challenges, for instance, recognizing the limitations of low testing proportions. We also highlight the early warning capacity that genomic surveillance in Africa has had for the rest of the world with the detection of new lineages and variants, the most recent being the characterization of various Omicron subvariants. CONCLUSION Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve. This is important not only to help combat SARS-CoV-2 on the continent but also because it can be used as a platform to help address the many emerging and reemerging infectious disease threats in Africa. In particular, capacity building for local sequencing within countries or within the continent should be prioritized because this is generally associated with shorter turnaround times, providing the most benefit to local public health authorities tasked with pandemic response and mitigation and allowing for the fastest reaction to localized outbreaks. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century

Prevalence of arps10, fd, pfmdr-2, pfcrt and pfkelch13 gene mutations in Plasmodium falciparum parasite population in Uganda.

In Uganda, Artemether-Lumefantrine and Artesunate are recommended for uncomplicated and severe malaria respectively, but are currently threatened by parasite resistance. Genetic and epigenetic factors play a role in predisposing Plasmodium falciparum parasites to acquiring Pfkelch13 (K13) mutations associated with delayed artemisinin parasite clearance as reported in Southeast Asia. In this study, we report on the prevalence of mutations in the K13, pfmdr-2 (P. falciparum multidrug resistance protein 2), fd (ferredoxin), pfcrt (P. falciparum chloroquine resistance transporter), and arps10 (apicoplast ribosomal protein S10) genes in Plasmodium falciparum parasites prior to (2005) and after (2013) introduction of artemisinin combination therapies for malaria treatment in Uganda. A total of 200 P. falciparum parasite DNA samples were screened. Parasite DNA was extracted using QIAamp DNA mini kit (Qiagen, GmbH, Germany) procedure. The PCR products were sequenced using Sanger dideoxy sequencing method. Of the 200 P. falciparum DNA samples screened, sequencing for mutations in K13, pfmdr-2, fd, pfcrt, arps10 genes was successful in 142, 186, 141, 128 and 74 samples respectively. Overall, we detected six (4.2%, 6/142; 95%CI: 1.4-7.0) K13 single nucleotide polymorphisms (SNPs), of which 3.9% (2/51), 4.4% (4/91) occurred in 2005 and 2013 samples respectively. All four K13 SNPs in 2013 samples were non-synonymous (A578S, E596V, S600C and E643K) while of the two SNPs in 2005 samples, one (Y588N) is non-synonymous and the other (I587I) is synonymous. There was no statistically significant difference in the prevalence of K13 (p = 0.112) SNPs in the samples collected in 2005 and 2013. The overall prevalence of SNPs in pfmdr-2 gene was 39.8% (74/186, 95%CI: 25.1-50.4). Of this, 4.2% (4/95), 76.9% (70/91) occurred in 2005 and 2013 samples respectively. In 2005 samples only one SNP, Y423F (4.2%, 4/95) was found while in 2013, Y423F (38.5%, 35/91) and I492V (38.5%, 35/91) SNPs in the pfmdr-2 gene were found. There was a statistically significant difference in the prevalence of pfmdr-2 SNPs in the samples collected in 2005 and 2013 (p<0.001). The overall prevalence of arps10 mutations was 2.7% (2/72, 95%CI: 0.3-4.2). Two mutations, V127M (4.5%: 1/22) and D128H (4.5%: 1/22) in the arps10 gene were each found in P. falciparum parasite samples collected in 2013. There was no statistically significant difference in the prevalence of arps10 SNPs in the samples collected in 2005 and 2013 (p = 0.238). There were more pfmdr-2 SNPs in P. falciparum parasites collected after introduction of Artemisinin combination therapies in malaria treatment. This is an indicator of the need for continuous surveillance to monitor emergence of molecular markers of artemisinin resistance and its potential drivers in malaria affected regions globally

No association of a Vascular endothelial growth factor A (VEGFA) gene polymorphism with pre-eclampsia among pregnant women in Uganda

Background: Vascular endothelial growth factor A (VEGFA) is a major angiogenic factor that plays an important role in the formation of blood vessels during embryonic development. VEGFA has been implicated in the pathophysiology of pre-eclampsia (PE), since pre-eclamptic women present with reduced levels of free circulating VEGFA. The 3' untranslated region (3'-UTR) of the VEGFA gene consists of elements that regulate the transcription and hence expression of the VEGFA protein in circulation. Hence it is suggested that variations thereof could underlie the reduced VEGFA levels observed in pre-eclamptic women. The purpose of this study was to investigate presence of the + 936C/T polymorphism, a common single nucleotide polymorphism (SNP) in the 3'-UTR of the VEGFA gene, and determine its association with PE among pregnant women in Uganda. Results: There was no significant difference observed in the allele and genotype frequencies of the + 936C/T 3' UTR-VEGFA polymorphism between pre-eclamptic and normotensive pregnant women (P &gt; 0.05). Additionally, there was no significant difference in the median plasma levels of free VEGFA among women with the wild type, CT and TT genotypes of the + 936C/T VEGFA polymorphism (median = 0.84 pg/mL (IQR = 0.39-1.41) Vs 1.05 (0.61-1.18) Vs 1.05 (1.05-1.05) respectively, p-value = 0.7161). Conclusions: These study findings indicate that the + 936C/T 3' UTR-VEGFA polymorphism had no significant association with increased susceptibility to PE among women in Uganda. Further studies with a larger sample size are recommended

No association of a Vascular endothelial growth factor A (VEGFA) gene polymorphism with pre-eclampsia among pregnant women in Uganda.

BACKGROUND: Vascular endothelial growth factor A (VEGFA) is a major angiogenic factor that plays an important role in the formation of blood vessels during embryonic development. VEGFA has been implicated in the pathophysiology of pre-eclampsia (PE), since pre-eclamptic women present with reduced levels of free circulating VEGFA. The 3' untranslated region (3'-UTR) of the VEGFA gene consists of elements that regulate the transcription and hence expression of the VEGFA protein in circulation. Hence it is suggested that variations thereof could underlie the reduced VEGFA levels observed in pre-eclamptic women. The purpose of this study was to investigate presence of the + 936C/T polymorphism, a common single nucleotide polymorphism (SNP) in the 3'-UTR of the VEGFA gene, and determine its association with PE among pregnant women in Uganda. RESULTS: There was no significant difference observed in the allele and genotype frequencies of the + 936C/T 3' UTR-VEGFA polymorphism between pre-eclamptic and normotensive pregnant women (P > 0.05). Additionally, there was no significant difference in the median plasma levels of free VEGFA among women with the wild type, CT and TT genotypes of the + 936C/T VEGFA polymorphism (median = 0.84 pg/mL (IQR = 0.39-1.41) Vs 1.05 (0.61-1.18) Vs 1.05 (1.05-1.05) respectively, p-value = 0.7161). CONCLUSIONS: These study findings indicate that the + 936C/T 3' UTR-VEGFA polymorphism had no significant association with increased susceptibility to PE among women in Uganda. Further studies with a larger sample size are recommended

No association of a Vascular endothelial growth factor A (VEGFA) gene polymorphism with pre-eclampsia among pregnant women in Uganda.

Acknowledgements: We acknowledge the support given by Edgar Kigozi on behalf of the Medical and Molecular Laboratory, College of Health Sciences, Makerere University, Kampala, Uganda, towards administration of study laboratory events, optimisation and performing of laboratory assays, and shipping of generated PCR products to sequencing company. We recognise the support provided by the study team, including; Maxine Atuheirwe, Rosemary Byenkya, Margaret Sewagaba, Elizabeth Mutesi, Ruth Namubiru, Emily Nakirijja, Moreen Ahimbisibwe, Florence Walabyeki, Noela Kalyowa, Rita Kukundakwe, and Doreen Birungi.BACKGROUND: Vascular endothelial growth factor A (VEGFA) is a major angiogenic factor that plays an important role in the formation of blood vessels during embryonic development. VEGFA has been implicated in the pathophysiology of pre-eclampsia (PE), since pre-eclamptic women present with reduced levels of free circulating VEGFA. The 3' untranslated region (3'-UTR) of the VEGFA gene consists of elements that regulate the transcription and hence expression of the VEGFA protein in circulation. Hence it is suggested that variations thereof could underlie the reduced VEGFA levels observed in pre-eclamptic women. The purpose of this study was to investigate presence of the + 936C/T polymorphism, a common single nucleotide polymorphism (SNP) in the 3'-UTR of the VEGFA gene, and determine its association with PE among pregnant women in Uganda. RESULTS: There was no significant difference observed in the allele and genotype frequencies of the + 936C/T 3' UTR-VEGFA polymorphism between pre-eclamptic and normotensive pregnant women (P > 0.05). Additionally, there was no significant difference in the median plasma levels of free VEGFA among women with the wild type, CT and TT genotypes of the + 936C/T VEGFA polymorphism (median = 0.84 pg/mL (IQR = 0.39-1.41) Vs 1.05 (0.61-1.18) Vs 1.05 (1.05-1.05) respectively, p-value = 0.7161). CONCLUSIONS: These study findings indicate that the + 936C/T 3' UTR-VEGFA polymorphism had no significant association with increased susceptibility to PE among women in Uganda. Further studies with a larger sample size are recommended

Identification of a Clade-Specific HLA-C*03:02 CTL Epitope GY9 Derived from the HIV-1 p17 Matrix Protein.

Efforts towards an effective HIV-1 vaccine have remained mainly unsuccessful. There is increasing evidence for a potential role of HLA-C-restricted CD8+ T cell responses in HIV-1 control, including our recent report of HLA-C*03:02 among African children. However, there are no documented optimal HIV-1 CD8+ T cell epitopes restricted by HLA-C*03:02; additionally, the structural influence of HLA-C*03:02 on epitope binding is undetermined. Immunoinformatics approaches provide a fast and inexpensive method to discover HLA-restricted epitopes. Here, we employed immunopeptidomics to identify HLA-C*03:02 CD8+ T cell epitopes. We identified a clade-specific Gag-derived GY9 (GTEELRSLY) HIV-1 p17 matrix epitope potentially restricted to HLA-C*03:02. Residues E62, T142, and E151 in the HLA-C*03:02 binding groove and positions p3, p6, and p9 on the GY9 epitope are crucial in shaping and stabilizing the epitope binding. Our findings support the growing evidence of the contribution of HLA-C molecules to HIV-1 control and provide a prospect for vaccine strategies

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