A blended genome and exome sequencing method captures genetic variation in an unbiased, high-quality, and cost-effective manner.

We deployed the Blended Genome Exome (BGE), a DNA library blending approach that generates low pass whole genome (1-4× mean depth) and deep whole exome (30-40× mean depth) data in a single sequencing run. This technology is cost-effective, empowers most genomic discoveries possible with deep whole genome sequencing, and provides an unbiased method to capture the diversity of common SNP variation across the globe. To evaluate this new technology at scale, we applied BGE to sequence >53,000 samples from the Populations Underrepresented in Mental Illness Associations Studies (PUMAS) Project, which included participants across African, African American, and Latin American populations. We evaluated the accuracy of BGE imputed genotypes against raw genotype calls from the Illumina Global Screening Array. All PUMAS cohorts had R 2 concordance ≥95% among SNPs with MAF≥1%, and never fell below ≥90% R 2 for SNPs with MAF<1%. Furthermore, concordance rates among local ancestries within two recently admixed cohorts were consistent among SNPs with MAF≥1%, with only minor deviations in SNPs with MAF<1%. We also benchmarked the discovery capacity of BGE to access protein-coding copy number variants (CNVs) against deep whole genome data, finding that deletions and duplications spanning at least 3 exons had a positive predicted value of ~90%. Our results demonstrate BGE scalability and efficacy in capturing SNPs, indels, and CNVs in the human genome at 28% of the cost of deep whole-genome sequencing. BGE is poised to enhance access to genomic testing and empower genomic discoveries, particularly in underrepresented populations.

Genetic structure correlates with ethnolinguistic diversity in eastern and southern Africa.

African populations are the most diverse in the world yet are sorely underrepresented in medical genetics research. Here, we examine the structure of African populations using genetic and comprehensive multi-generational ethnolinguistic data from the Neuropsychiatric Genetics of African Populations-Psychosis study (NeuroGAP-Psychosis) consisting of 900 individuals from Ethiopia, Kenya, South Africa, and Uganda. We find that self-reported language classifications meaningfully tag underlying genetic variation that would be missed with consideration of geography alone, highlighting the importance of culture in shaping genetic diversity. Leveraging our uniquely rich multi-generational ethnolinguistic metadata, we track language transmission through the pedigree, observing the disappearance of several languages in our cohort as well as notable shifts in frequency over three generations. We find suggestive evidence for the rate of language transmission in matrilineal groups having been higher than that for patrilineal ones. We highlight both the diversity of variation within Africa as well as how within-Africa variation can be informative for broader variant interpretation; many variants that are rare elsewhere are common in parts of Africa. The work presented here improves the understanding of the spectrum of genetic variation in African populations and highlights the enormous and complex genetic and ethnolinguistic diversity across Africa.

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.

Whole-genome sequencing of SARS-CoV-2 in Uganda: implementation of the low-cost ARTIC protocol in resource-limited settings.

Background: In January 2020, a previously unknown coronavirus strain was identified as the cause of a severe acute respiratory syndrome (SARS-CoV-2). The first viral whole-genome was sequenced using high-throughput sequencing from a sample collected in Wuhan, China. Whole-genome sequencing (WGS) is imperative in investigating disease outbreak transmission dynamics and guiding decision-making in public health. Methods: We retrieved archived SARS-CoV-2 samples at the Integrated Biorepository of H3Africa Uganda, Makerere University (IBRH3AU). These samples were collected previously from individuals diagnosed with coronavirus disease 2019 (COVID-19) using real-time reverse transcription quantitative polymerase chain reaction (RT-qPCR). 30 samples with cycle thresholds (Cts) values <25 were selected for WGS using SARS-CoV-2 ARTIC protocol at Makerere University Molecular Diagnostics Laboratory. Results: 28 out of 30 (93.3%) samples generated analyzable genomic sequence reads. We detected SARS-CoV-2 and lineages A (22/28) and B (6/28) from the samples. We further show phylogenetic relatedness of these isolates alongside other 328 Uganda (lineage A = 222, lineage B = 106) SARS-CoV-2 genomes available in GISAID by April 22, 2021 and submitted by the Uganda Virus Research Institute. Conclusions: Our study demonstrated adoption and optimization of the low-cost ARTIC SARS-CoV-2 WGS protocol in a resource limited laboratory setting. This work has set a foundation to enable rapid expansion of SARS-CoV-2 WGS in Uganda as part of the Presidential Scientific Initiative on Epidemics (PRESIDE) CoV-bank project and IBRH3AU.

Assessing a transmission network of Mycobacterium tuberculosis in an African city using single nucleotide polymorphism threshold analysis.

Tuberculosis (TB) is the leading cause of death in humans by a single infectious agent worldwide with approximately two billion humans latently infected with the bacterium Mycobacterium tuberculosis. Currently, the accepted method for controlling the disease is Tuberculosis Directly Observed Treatment Shortcourse (TB-DOTS). This program is not preventative and individuals may transmit disease before diagnosis, thus better understanding of disease transmission is essential. Using whole-genome sequencing and single nucleotide polymorphism analysis, we analyzed genomes of 145 M. tuberculosis clinical isolates from active TB cases from the Rubaga Division of Kampala, Uganda. We established that these isolates grouped into M. tuberculosis complex (MTBC) lineages 1, 2, 3, and 4, with the most isolates grouping into lineage 4. Possible transmission pairs containing ≤12 SNPs were identified in lineages 1, 3, and 4 with the prevailing transmission in lineages 3 and 4. Furthermore, investigating DNA codon changes as a result of specific SNPs in prominent virulence genes including plcA and plcB could indicate potentially important modifications in protein function. Incorporating this analysis with corresponding epidemiological data may provide a blueprint for the integration of public health interventions to decrease TB transmission in a region.

Low-coverage sequencing cost-effectively detects known and novel variation in underrepresented populations.

Genetic studies in underrepresented populations identify disproportionate numbers of novel associations. However, most genetic studies use genotyping arrays and sequenced reference panels that best capture variation most common in European ancestry populations. To compare data generation strategies best suited for underrepresented populations, we sequenced the whole genomes of 91 individuals to high coverage as part of the Neuropsychiatric Genetics of African Population-Psychosis (NeuroGAP-Psychosis) study with participants from Ethiopia, Kenya, South Africa, and Uganda. We used a downsampling approach to evaluate the quality of two cost-effective data generation strategies, GWAS arrays versus low-coverage sequencing, by calculating the concordance of imputed variants from these technologies with those from deep whole-genome sequencing data. We show that low-coverage sequencing at a depth of ≥4× captures variants of all frequencies more accurately than all commonly used GWAS arrays investigated and at a comparable cost. Lower depths of sequencing (0.5-1×) performed comparably to commonly used low-density GWAS arrays. Low-coverage sequencing is also sensitive to novel variation; 4× sequencing detects 45% of singletons and 95% of common variants identified in high-coverage African whole genomes. Low-coverage sequencing approaches surmount the problems induced by the ascertainment of common genotyping arrays, effectively identify novel variation particularly in underrepresented populations, and present opportunities to enhance variant discovery at a cost similar to traditional approaches.

Variations in Trim5α and Cyclophilin A genes among HIV-1 elite controllers and non controllers in Uganda: a laboratory-based cross-sectional study.

BACKGROUND: Tripartite Motif Containing 5 alpha (TRIM5α), a restriction factor produced ubiquitously in cells and tissues of the body plays an important role in the immune response against HIV. TRIM5α targets the HIV capsid for proteosomal destruction. Cyclophilin A, an intracellular protein has also been reported to influence HIV infectivity in a cell-specific manner. Accordingly, variations in TRIM5α and Cyclophilin A genes have been documented to influence HIV-1 disease progression. However, these variations have not been documented among Elite controllers in Uganda and whether they play a role in viral suppression remains largely undocumented. Our study focused on identifying the variations in TRIM5α and Cyclophilin A genes among HIV-1 Elite controllers and non-controllers in Uganda. RESULTS: From the sequence analysis, the rs10838525 G > A mutation in exon 2 of TRIM5α was only found among elite controllers (30%) while the rs3824949 in the 5'UTR was seen among 25% of the non-controllers. In the Cyclophilin A promoter, rs6850 was seen among 62.5% of the non-controllers and only among 10% elite controllers. Furthermore, rs17860048 in the Cyclophillin A promoter was predominantly seen among elite controllers (30%) and 12.5% non-controllers. From gene expression analysis, we noted that the respective genes were generally elevated among elite controllers, however, this difference was not statistically significant (TRIM5α p = 0.6095; Cyclophilin A p = 0.6389). CONCLUSION: Variations in TRIM5α and Cyclophillin A promoter may influence HIV viral suppression. The rs10838525 SNP in TRIM5α may contribute to viral suppression among HIV-1 elite controllers. The rs6850 in the cyclophillin A gene may be responsible for HIV-1 rapid progression among HIV-1 non-controllers. These SNPs should be investigated mechanistically to determine their precise role in HIV-1 viral suppression.

CA-MRSA and HA-MRSA coexist in community and hospital settings in Uganda.

Background: Methicillin resistant Staphylococcus aureus (MRSA) strains were once confined to hospitals however, in the last 20 years MRSA infections have emerged in the community in people with no prior exposure to hospitals. Strains causing such infections were novel and referred to as community-associated MRSA (CA-MRSA). The aim of this study was to determine the MRSA carriage rate in children in eastern Uganda, and to investigate coexistence between CA-MRSA and hospital-associated (HA-MRSA). Methods: Between February and October 2011, nasopharyngeal samples (one per child) from 742 healthy children under 5 years in rural eastern Uganda were processed for isolation of MRSA, which was identified based on inhibition zone diameter of ≤19 mm on 30 µg cefoxitin disk. SCCmec and spa typing were performed for MRSA isolates. Results: A total of 140 S. aureus isolates (18.9%, 140/742) were recovered from the children of which 5.7% (42/742) were MRSA. Almost all (95.2%, 40/42) MRSA isolates were multidrug resistant (MDR). The most prevalent SCCmec elements were types IV (40.5%, 17/42) and I (38.1%, 16/42). The overall frequency of SCCmec types IV and V combined, hence CA-MRSA, was 50% (21/42). Likewise, the overall frequency of SCCmec types I, II and III combined, hence HA-MRSA, was 50% (21/42). Spa types t002, t037, t064, t4353 and t12939 were detected and the most frequent were t064 (19%, 8/42) and t037 (12%, 5/42). Conclusion: The MRSA carriage rate in children in eastern Uganda is high (5.7%) and comparable to estimates for Mulago Hospital in Kampala city. Importantly, HA-MRSA (mainly of spa type t037) and CA-MRSA (mainly of spa type t064) coexist in children in the community in eastern Uganda, and due to high proportion of MDR detected, outpatient treatment of MRSA infection in eastern Uganda might be difficult.