Genetic Association and Transferability for Urinary Albumin-Creatinine Ratio as a Marker of Kidney Disease in four Sub-Saharan African Populations and non-continental Individuals of African Ancestry

BACKGROUND

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

METHODS

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

RESULTS

Two genome-wide significant (P<5x10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

CONCLUSION

This study contributes novel insights into the genetic architecture of kidney function in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations.

Expanding the human gut microbiome atlas of Africa

Population studies are crucial in understanding the complex interplay between the gut microbiome and geographical, lifestyle, genetic, and environmental factors. However, populations from low- and middle-income countries, which represent ~84% of the world population, have been excluded from large-scale gut microbiome research. Here, we present the AWI-Gen 2 Microbiome Project, a cross-sectional gut microbiome study sampling 1,803 women from Burkina Faso, Ghana, Kenya, and South Africa. By intensively engaging with communities that range from rural and horticultural to urban informal settlements and post-industrial, we capture population diversity that represents a far greater breadth of the world's population. Using shotgun metagenomic sequencing, we find that study site explains substantially more microbial variation than disease status. We identify taxa with strong geographic and lifestyle associations, including loss of Treponema and Cryptobacteroides species and gain of Bifidobacterium species in urban populations. We uncover a wealth of prokaryotic and viral novelty, including 1,005 new bacterial metagenome-assembled genomes, and identify phylogeography signatures in Treponema succinifaciens. Finally, we find a microbiome signature of HIV infection that is defined by several taxa not previously associated with HIV, including Dysosmobacter welbionis and Enterocloster sp. This study represents the largest population-representative survey of gut metagenomes of African individuals to date, and paired with extensive clinical biomarkers, demographic data, and lifestyle information, provides extensive opportunity for microbiome-related discovery and research.

Characterization of CYP2B6 and CYP2A6 Pharmacogenetic Variation in Sub-Saharan African Populations

Genetic variation in CYP2B6 and CYP2A6 is known to impact interindividual response to antiretrovirals, nicotine, and bupropion, among other drugs. However, the full catalogue of clinically relevant pharmacogenetic variants in these genes is yet to be established, especially across African populations. This study therefore aimed to characterize the star allele (haplotype) distribution in CYP2B6 and CYP2A6 across diverse and understudied sub-Saharan African (SSA) populations. We called star alleles from 961 high-depth full genomes using StellarPGx, Aldy, and PyPGx. In addition, we performed CYP2B6 and CYP2A6 star allele frequency comparisons between SSA and other global biogeographical groups represented in the new 1000 Genomes Project high-coverage dataset (n = 2,000). This study presents frequency information for star alleles in CYP2B6 (e.g., *6 and *18; frequency of 21-47% and 2-19%, respectively) and CYP2A6 (e.g., *4, *9, and *17; frequency of 0-6%, 3-10%, and 6-20%, respectively), and predicted phenotypes (for CYP2B6), across various African populations. In addition, 50 potentially novel African-ancestry star alleles were computationally predicted by StellarPGx in CYP2B6 and CYP2A6 combined. For each of these genes, over 4% of the study participants had predicted novel star alleles. Three novel star alleles in CYP2A6 (*54, *55, and *56) and CYP2B6 apiece, and several suballeles were further validated via targeted Single-Molecule Real-Time resequencing. Our findings are important for informing the design of comprehensive pharmacogenetic testing platforms, and are highly relevant for personalized medicine strategies, especially relating to antiretroviral medication and smoking cessation treatment in Africa and the African diaspora. More broadly, this study highlights the importance of sampling diverse African ethnolinguistic groups for accurate characterization of the pharmacogene variation landscape across the continent.

Genome-wide association study meta-analysis of blood pressure traits and hypertension in sub-Saharan African populations: an AWI-Gen study

Most hypertension-related genome-wide association studies (GWASs) focus on non-African populations, despite hypertension (a major risk factor for cardiovascular disease) being highly prevalent in Africa. The AWI-Gen study GWAS meta-analysis for blood pressure (BP)-related traits (systolic and diastolic BP, pulse pressure, mean-arterial pressure and hypertension) from three sub-Saharan African geographic regions (N = 10,775), identifies two novel genome-wide significant signals (p < 5E-08): systolic BP near P2RY1 (rs77846204; intergenic variant, p = 4.95E-08) and pulse pressure near LINC01256 (rs80141533; intergenic variant, p = 1.76E-08). No genome-wide signals are detected for the AWI-Gen GWAS meta-analysis with previous African-ancestry GWASs (UK Biobank (African), Uganda Genome Resource). Suggestive signals (p < 5E-06) are observed for all traits, with 29 SNPs associating with more than one trait and several replicating known associations. Polygenic risk scores (PRSs) developed from studies on different ancestries have limited transferability, with multi-ancestry PRS providing better prediction. This study provides insights into the genetics of BP variation in African populations.

Multimorbidity in African ancestry populations: a scoping review

OBJECTIVES

Multimorbidity (MM) is a growing concern linked to poor outcomes and higher healthcare costs. While most MM research targets European ancestry populations, the prevalence and patterns in African ancestry groups remain underexplored. This study aimed to identify and summarise the available literature on MM in populations with African ancestry, on the continent, and in the diaspora.

DESIGN

A scoping review was conducted in five databases (PubMed, Web of Science, Scopus, Science Direct and JSTOR) in July 2022. Studies were selected based on predefined criteria, with data extraction focusing on methodology and findings. Descriptive statistics summarised the data, and a narrative synthesis highlighted key themes.

RESULTS

Of the 232 publications on MM in African-ancestry groups from 2010 to June 2022-113 examined continental African populations, 100 the diaspora and 19 both. Findings revealed diverse MM patterns within and beyond continental Africa. Cardiovascular and metabolic diseases are predominant in both groups (80% continental and 70% diaspora). Infectious diseases featured more in continental studies (58% continental and 16% diaspora). Although many papers did not specifically address these features, as in previous studies, older age, being women and having a lower socioeconomic status were associated with a higher prevalence of MM, with important exceptions. Research gaps identified included limited data on African-ancestry individuals, inadequate representation, under-represented disease groups, non-standardised methodologies, the need for innovative data strategies, and insufficient translational research.

CONCLUSION

The growing global MM prevalence is mirrored in African-ancestry populations. Recognising the unique contexts of African-ancestry populations is essential when addressing the burden of MM. This review emphasises the need for additional research to guide and enhance healthcare approaches for African-ancestry populations, regardless of their geographic location.

Characterization of CYP2C19 pharmacogenetic variation in African populations and comparison with other global populations

Background: CYP2C19 is important in the metabolism of clopidogrel and several antidepressants. This study aimed to characterize the distribution of CYP2C19 star alleles (haplotypes) across diverse African populations compared with global populations. Methods: CYP2C19 star alleles and diplotypes were called from high coverage genomes using the StellarPGx pipeline. Results: CYP2C19*1 (51%), *2 (17%) and *17 (22%) were the most common star alleles across African populations in this study. It was observed that 3% of African participants had potentially novel CYP2C19 haplotypes. Conclusion: This study supports the necessity for CYP2C19 pharmacogenetic testing in African and global clinical settings, as well as the importance of comprehensive star allele characterization in the African context.

Performance and accuracy evaluation of reference panels for genotype imputation in sub-Saharan African populations

Based on evaluations of imputation performed on a genotype dataset consisting of about 11,000 sub-Saharan African (SSA) participants, we show Trans-Omics for Precision Medicine (TOPMed) and the African Genome Resource (AGR) to be currently the best panels for imputing SSA datasets. We report notable differences in the number of single-nucleotide polymorphisms (SNPs) that are imputed by different panels in datasets from East, West, and South Africa. Comparisons with a subset of 95 SSA high-coverage whole-genome sequences (WGSs) show that despite being about 20-fold smaller, the AGR imputed dataset has higher concordance with the WGSs. Moreover, the level of concordance between imputed and WGS datasets was strongly influenced by the extent of Khoe-San ancestry in a genome, highlighting the need for integration of not only geographically but also ancestrally diverse WGS data in reference panels for further improvement in imputation of SSA datasets. Approaches that integrate imputed data from different panels could also lead to better imputation.

Characterization of CYP2D6 Pharmacogenetic Variation in Sub-Saharan African Populations

Cytochrome P450 2D6 (CYP2D6) is a key enzyme in drug response owing to its involvement in the metabolism of ~ 25% of clinically prescribed medications. The encoding CYP2D6 gene is highly polymorphic, and many pharmacogenetics studies have been performed worldwide to investigate the distribution of CYP2D6 star alleles (haplotypes); however, African populations have been relatively understudied to date. In this study, the distributions of CYP2D6 star alleles and predicted drug metabolizer phenotypes-derived from activity scores-were examined across multiple sub-Saharan African populations based on bioinformatics analysis of 961 high-depth whole genome sequences. This was followed by characterization of novel star alleles and suballeles in a subset of the participants via targeted high-fidelity Single-Molecule Real-Time resequencing (Pacific Biosciences). This study revealed varying frequencies of known CYP2D6 alleles and predicted phenotypes across different African ethnolinguistic groups. Twenty-seven novel CYP2D6 star alleles were predicted computationally and two of them were further validated. This study highlights the importance of studying variation in key pharmacogenes such as CYP2D6 in the African context to better understand population-specific allele frequencies. This will aid in the development of better genotyping panels and star allele detection approaches with a view toward supporting effective implementation of precision medicine strategies in Africa and across the African diaspora.

Genome-wide Association Study Meta-analysis of Blood Pressure Traits and Hypertension in Sub-Saharan African Populations: An AWI-Gen Study

Most hypertension-related genome-wide association studies (GWAS) focus on non-African populations, despite hypertension (a major risk factor for cardiovascular disease) being highly prevalent in Africa. The AWI-Gen study GWAS meta-analysis for blood pressure-related traits (systolic and diastolic blood pressure, pulse pressure, mean-arterial pressure and hypertension) from three sub-Saharan African geographic regions (N=10,775), identified two genome-wide significant signals (p<5E-08): systolic blood pressure near P2RY1 (rs77846204; intergenic variant, p=4.25E-08) and pulse pressure near Linc01256 (rs80141533; intergenic variant, p=4.25E-08). No genome-wide signals were detected for the AWI-Gen GWAS meta-analysis with previous African-ancestry GWASs (UK Biobank (African), Uganda Genome Resource). Suggestive signals (p<5E-06) were observed for all traits, with 29 displaying pleiotropic effects and several replicating known associations. Polygenic risk scores developed from studies on different ancestries had limited transferability, with multi-ancestry models providing better prediction. This study provides insights into the genetics and physiology of blood pressure variation in African populations.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

H3AGWAS: a portable workflow for genome wide association studies

Background

Genome-wide association studies (GWAS) are a powerful method to detect associations between variants and phenotypes. A GWAS requires several complex computations with large data sets, and many steps may need to be repeated with varying parameters. Manual running of these analyses can be tedious, error-prone and hard to reproduce.

Results

The H3AGWAS workflow from the Pan-African Bioinformatics Network for H3Africa is a powerful, scalable and portable workflow implementing pre-association analysis, implementation of various association testing methods and post-association analysis of results.

Conclusions

The workflow is scalable—laptop to cluster to cloud (e.g., SLURM, AWS Batch, Azure). All required software is containerised and can run under Docker or Singularity.

Characterization of POR haplotype distribution in African populations and comparison with other global populations

Background & aim: POR is an enzyme that mediates electron transfer to enable the drug-metabolizing activity of CYP450 proteins. However, POR has been understudied in pharmacogenomics despite this vital role. This study aimed to characterize the genetic variation in POR across African populations and to compare the star allele (haplotype) distribution with that in other global populations. Materials & methods: POR star alleles were called from whole-genome sequencing data using the StellarPGx pipeline. Results: In addition to the common POR*1 and *28 (defined by rs1057868), five novel rare haplotypes were computationally inferred. No significant frequency differences were observed among the majority of African populations. However, POR*28 was observed at a higher frequency in individuals of non-African ancestry. Conclusion: This study highlights the distribution of POR alleles in Africa and across global populations with a view toward informing future precision medicine implementation.

Electronic Data Capture System (REDCap) for Health Care Research and Training in a Resource-Constrained Environment: Technology Adoption Case Study

Background

Electronic data capture (EDC) in academic health care organizations provides an opportunity for the management, aggregation, and secondary use of research and clinical data. It is especially important in resource-constrained environments such as the South African public health care sector, where paper records are still the main form of clinical record keeping.

Objective

The aim of this study was to describe the strategies followed by the University of the Witwatersrand Faculty of Health Sciences (Wits FHS) during the period from 2013 to 2021 to overcome resistance to, and encourage the adoption of, the REDCap (Research Electronic Data Capture; Vanderbilt University) system by academic and clinical staff. REDCap has found wide use in varying domains, including clinical studies and research projects as well as administrative, financial, and human resource applications. Given REDCap’s global footprint in >5000 institutions worldwide and potential for future growth, the strategies followed by the Wits FHS to support users and encourage adoption may be of importance to others using the system, particularly in resource-constrained settings.

Methods

The strategies to support users and encourage adoption included top-down organizational support; secure and reliable application, hosting infrastructure, and systems administration; an enabling and accessible REDCap support team; regular hands-on training workshops covering REDCap project setup and data collection instrument design techniques; annual local symposia to promote networking and awareness of all the latest software features and best practices for using them; participation in REDCap Consortium activities; and regular and ongoing mentorship from members of the Vanderbilt University Medical Center.

Results

During the period from 2013 to 2021, the use of the REDCap EDC system by individuals at the Wits FHS increased, respectively, from 129 active user accounts to 3447 active user accounts. The number of REDCap projects increased from 149 in 2013 to 12,865 in 2021. REDCap at Wits also supported various publications and research outputs, including journal articles and postgraduate monographs. As of 2020, a total of 233 journal articles and 87 postgraduate monographs acknowledged the use of the Wits REDCap system.

Conclusions

By providing reliable infrastructure and accessible support resources, we were able to successfully implement and grow the REDCap EDC system at the Wits FHS and its associated academic medical centers. We believe that the increase in the use of REDCap was driven by offering a dependable, secure service with a strong end-user training and support model. This model may be applied by other academic and health care organizations in resource-constrained environments planning to implement EDC technology.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Molecular dynamics of G6PD variants from sub-Saharan Africa

Precision medicine uses genomic guidance to improve drug treatment safety and efficacy. Prior knowledge of genetic variant impact can enable such strategies, but current knowledge of African variants remains scarce. G6PD variants are linked to haemolytic adverse effects for a number of drugs commonly used in African populations. We have investigated a set of G6PD variants with structural bioinformatics techniques to further characterise variants with known effect, and gain insights into variants with unknown impact. We observed wide variations in patterns of root-mean-square deviation between wild-type and variant structures. Variants with known, highly deleterious impact show structural effects which may likely result in the destabilisation of the G6PD homodimer. The V68M and N126D variants (which are both common across African populations, and together form the A- haplotype) induce large conformational shifts in the catalytic NADP+ binding domain. We observed a greater impact for the haplotype than for each of the individual variants in these cases. A novel African variant (M207T) shows the potential to disrupt interactions within the protein core, urging further investigation. We explore how characterising the molecular impact of African G6PD variants can enable advanced strategies for precision medicine, as well as impact the use of novel therapeutics aiming to treat G6PD deficiency. This knowledge can assist in bridging current knowledge gaps, and aid to facilitate precision medicine applications in African populations.

•

Assessment of African G6PD variation with structural bioinformatics.

•

Molecular dynamics of 500 ns to explore molecular motions.

•

Comparison of variants with known/unknown impact.

•

Exploring mechanisms of impact.

•

Knowledge building to enable G6PD precision medicine in Africa.

Meta-analysis of sub-Saharan African studies provides insights into genetic architecture of lipid traits

Genetic associations for lipid traits have identified hundreds of variants with clear differences across European, Asian and African studies. Based on a sub-Saharan-African GWAS for lipid traits in the population cross-sectional AWI-Gen cohort (N = 10,603) we report a novel LDL-C association in the GATB region (P-value=1.56 × 10⁻⁸). Meta-analysis with four other African cohorts (N = 23,718) provides supporting evidence for the LDL-C association with the GATB/FHIP1A region and identifies a novel triglyceride association signal close to the FHIT gene (P-value =2.66 × 10⁻⁸). Our data enable fine-mapping of several well-known lipid-trait loci including LDLR, PMFBP1 and LPA. The transferability of signals detected in two large global studies (GLGC and PAGE) consistently improves with an increase in the size of the African replication cohort. Polygenic risk score analysis shows increased predictive accuracy for LDL-C levels with the narrowing of genetic distance between the discovery dataset and our cohort. Novel discovery is enhanced with the inclusion of African data.

Genetic associations and polygenic scores for lipid traits have low transferability to African individuals. Here, the authors perform a large sub-Sarahan African lipid GWAS and find that larger datasets and better global representation in discovery GWAS help to bridge this gap.

Microbiomics: The Next Pillar of Precision Medicine and Its Role in African Healthcare

Limited access to technologies that support early monitoring of disease risk and a poor understanding of the geographically unique biological and environmental factors underlying disease, represent significant barriers to improved health outcomes and precision medicine efforts in low to middle income countries. These challenges are further compounded by the rich genetic diversity harboured within Southern Africa thus necessitating alternative strategies for the prediction of disease risk and clinical outcomes in regions where accessibility to personalized healthcare remains limited. The human microbiome refers to the community of microorganisms (bacteria, archaea, fungi and viruses) that co-inhabit the human body. Perturbation of the natural balance of the gut microbiome has been associated with a number of human pathologies, and the microbiome has recently emerged as a critical determinant of drug pharmacokinetics and immunomodulation. The human microbiome should therefore not be omitted from any comprehensive effort towards stratified healthcare and would provide an invaluable and orthogonal approach to existing precision medicine strategies. Recent studies have highlighted the overarching effect of geography on gut microbial diversity as it relates to human health. Health insights from international microbiome datasets are however not yet verified in context of the vast geographical diversity that exists throughout the African continent. In this commentary we discuss microbiome research in Africa and its role in future precision medicine initiatives across the African continent.

Short- and long-read metagenomics of urban and rural South African gut microbiomes reveal a transitional composition and undescribed taxa

Human gut microbiome research focuses on populations living in high-income countries and to a lesser extent, non-urban agriculturalist and hunter-gatherer societies. The scarcity of research between these extremes limits our understanding of how the gut microbiota relates to health and disease in the majority of the world's population. Here, we evaluate gut microbiome composition in transitioning South African populations using short- and long-read sequencing. We analyze stool from adult females living in rural Bushbuckridge (n = 118) or urban Soweto (n = 51) and find that these microbiomes are taxonomically intermediate between those of individuals living in high-income countries and traditional communities. We demonstrate that reference collections are incomplete for characterizing microbiomes of individuals living outside high-income countries, yielding artificially low beta diversity measurements, and generate complete genomes of undescribed taxa, including Treponema, Lentisphaerae, and Succinatimonas. Our results suggest that the gut microbiome of South Africans does not conform to a simple "western-nonwestern" axis and contains undescribed microbial diversity.

Genetic associations with carotid intima-media thickness link to atherosclerosis with sex-specific effects in sub-Saharan Africans

Atherosclerosis precedes the onset of clinical manifestations of cardiovascular diseases (CVDs). We used carotid intima-media thickness (cIMT) to investigate genetic susceptibility to atherosclerosis in 7894 unrelated adults (3963 women, 3931 men; 40 to 60 years) resident in four sub-Saharan African countries. cIMT was measured by ultrasound and genotyping was performed on the H3Africa SNP Array. Two new African-specific genome-wide significant loci for mean-max cIMT, SIRPA (p = 4.7E-08), and FBXL17 (p = 2.5E-08), were identified. Sex-stratified analysis revealed associations with one male-specific locus, SNX29 (p = 6.3E-09), and two female-specific loci, LARP6 (p = 2.4E-09) and PROK1 (p = 1.0E-08). We replicate previous cIMT associations with different lead SNPs in linkage disequilibrium with SNPs primarily identified in European populations. Our study find significant enrichment for genes involved in oestrogen response from female-specific signals. The genes identified show biological relevance to atherosclerosis and/or CVDs, sex-differences and transferability of signals from non-African studies.

nf-rnaSeqCount: A Nextflow pipeline for obtaining raw read counts from RNA-seq data

The rate of raw sequence production through Next-Generation Sequencing (NGS) has been growing exponentially due to improved technology and reduced costs. This has enabled researchers to answer many biological questions through “multi-omics” data analyses. Even though such data promises new insights into how biological systems function and understanding disease mechanisms, computational analyses performed on such large datasets comes with its challenges and potential pitfalls. The aim of this study was to develop a robust portable and reproducible bioinformatic pipeline for the automation of RNA sequencing (RNA-seq) data analyses. Using Nextflow as a workflow management system and Singularity for application containerisation, the nf-rnaSeqCount pipeline was developed for mapping raw RNA-seq reads to a reference genome and quantifying abundance of identified genomic features for differential gene expression analyses. The pipeline provides a quick and efficient way to obtain a matrix of read counts that can be used with tools such as DESeq2 and edgeR for differential expression analysis. Robust and flexible bioinformatic and computational pipelines for RNA-seq data analysis, from QC to sequence alignment and comparative analyses, will reduce analysis time, and increase accuracy and reproducibility of findings to promote transcriptome research.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Single Nucleotide Polymorphism Induces Divergent Dynamic Patterns in CYP3A5: A Microsecond Scale Biomolecular Simulation of Variants Identified in Sub-Saharan African Populations

Pharmacogenomics aims to reveal variants associated with drug response phenotypes. Genes whose roles involve the absorption, distribution, metabolism, and excretion of drugs, are highly polymorphic between populations. High coverage whole genome sequencing showed that a large proportion of the variants for these genes are rare in African populations. This study investigated the impact of such variants on protein structure to assess their functional importance. We used genetic data of CYP3A5 from 458 individuals from sub-Saharan Africa to conduct a structural bioinformatics analysis. Five missense variants were modeled and microsecond scale molecular dynamics simulations were conducted for each, as well as for the CYP3A5 wildtype and the Y53C variant, which has a known deleterious impact on enzyme activity. The binding of ritonavir and artemether to CYP3A5 variant structures was also evaluated. Our results showed different conformational characteristics between all the variants. No significant structural changes were noticed. However, the genetic variability seemed to act on the plasticity of the protein. The impact on drug binding might be drug dependant. We concluded that rare variants hold relevance in determining the pharmacogenomics properties of populations. This could have a significant impact on precision medicine applications in sub-Saharan Africa.

Artificial Intelligence-Based Screening for Mycobacteria in Whole-Slide Images of Tissue Samples

OBJECTIVES

This study aimed to develop and validate a deep learning algorithm to screen digitized acid fast-stained (AFS) slides for mycobacteria within tissue sections.

METHODS

A total of 441 whole-slide images (WSIs) of AFS tissue material were used to develop a deep learning algorithm. Regions of interest with possible acid-fast bacilli (AFBs) were displayed in a web-based gallery format alongside corresponding WSIs for pathologist review. Artificial intelligence (AI)-assisted analysis of another 138 AFS slides was compared to manual light microscopy and WSI evaluation without AI support.

RESULTS

Algorithm performance showed an area under the curve of 0.960 at the image patch level. More AI-assisted reviews identified AFBs than manual microscopy or WSI examination (P < .001). Sensitivity, negative predictive value, and accuracy were highest for AI-assisted reviews. AI-assisted reviews also had the highest rate of matching the original sign-out diagnosis, were less time-consuming, and were much easier for pathologists to perform (P < .001).

CONCLUSIONS

This study reports the successful development and clinical validation of an AI-based digital pathology system to screen for AFBs in anatomic pathology material. AI assistance proved to be more sensitive and accurate, took pathologists less time to screen cases, and was easier to use than either manual microscopy or viewing WSIs.

The Extent and Impact of Variation in ADME Genes in Sub-Saharan African Populations

Introduction: Investigating variation in genes involved in the absorption, distribution, metabolism, and excretion (ADME) of drugs are key to characterizing pharmacogenomic (PGx) relationships. ADME gene variation is relatively well characterized in European and Asian populations, but data from African populations are under-studied-which has implications for drug safety and effective use in Africa. Results: We identified significant ADME gene variation in African populations using data from 458 high-coverage whole genome sequences, 412 of which are novel, and from previously available African sequences from the 1,000 Genomes Project. ADME variation was not uniform across African populations, particularly within high impact coding variation. Copy number variation was detected in 116 ADME genes, with equal ratios of duplications/deletions. We identified 930 potential high impact coding variants, of which most are discrete to a single African population cluster. Large frequency differences (i.e., >10%) were seen in common high impact variants between clusters. Several novel variants are predicted to have a significant impact on protein structure, but additional functional work is needed to confirm the outcome of these for PGx use. Most variants of known clinical outcome are rare in Africa compared to European populations, potentially reflecting a clinical PGx research bias to European populations. Discussion: The genetic diversity of ADME genes across sub-Saharan African populations is large. The Southern African population cluster is most distinct from that of far West Africa. PGx strategies based on European variants will be of limited use in African populations. Although established variants are important, PGx must take into account the full range of African variation. This work urges further characterization of variants in African populations including in vitro and in silico studies, and to consider the unique African ADME landscape when developing precision medicine guidelines and tools for African populations.

Genetic substructure and complex demographic history of South African Bantu speakers

South Eastern Bantu-speaking (SEB) groups constitute more than 80% of the population in South Africa. Despite clear linguistic and geographic diversity, the genetic differences between these groups have not been systematically investigated. Based on genome-wide data of over 5000 individuals, representing eight major SEB groups, we provide strong evidence for fine-scale population structure that broadly aligns with geographic distribution and is also congruent with linguistic phylogeny (separation of Nguni, Sotho-Tswana and Tsonga speakers). Although differential Khoe-San admixture plays a key role, the structure persists after Khoe-San ancestry-masking. The timing of admixture, levels of sex-biased gene flow and population size dynamics also highlight differences in the demographic histories of individual groups. The comparisons with five Iron Age farmer genomes further support genetic continuity over ~400 years in certain regions of the country. Simulated trait genome-wide association studies further show that the observed population structure could have major implications for biomedical genomics research in South Africa.

A Digital Pathology Solution to Resolve the Tissue Floater Conundrum

CONTEXT.—

Pathologists may encounter extraneous pieces of tissue (tissue floaters) on glass slides because of specimen cross-contamination. Troubleshooting this problem, including performing molecular tests for tissue identification if available, is time consuming and often does not satisfactorily resolve the problem.

OBJECTIVE.—

To demonstrate the feasibility of using an image search tool to resolve the tissue floater conundrum.

DESIGN.—

A glass slide was produced containing 2 separate hematoxylin and eosin (H&E)-stained tissue floaters. This fabricated slide was digitized along with the 2 slides containing the original tumors used to create these floaters. These slides were then embedded into a dataset of 2325 whole slide images comprising a wide variety of H&E stained diagnostic entities. Digital slides were broken up into patches and the patch features converted into barcodes for indexing and easy retrieval. A deep learning-based image search tool was employed to extract features from patches via barcodes, hence enabling image matching to each tissue floater.

RESULTS.—

There was a very high likelihood of finding a correct tumor match for the queried tissue floater when searching the digital database. Search results repeatedly yielded a correct match within the top 3 retrieved images. The retrieval accuracy improved when greater proportions of the floater were selected. The time to run a search was completed within several milliseconds.

CONCLUSIONS.—

Using an image search tool offers pathologists an additional method to rapidly resolve the tissue floater conundrum, especially for those laboratories that have transitioned to going fully digital for primary diagnosis.

Author Correction: High-depth African genomes inform human migration and health

A Correction to this paper has been published: https://doi.org/10.1038/s41586-021-03286-9.

StellarPGx: A Nextflow Pipeline for Calling Star Alleles in Cytochrome P450 Genes

Bioinformatics pipelines for calling star alleles (haplotypes) in cytochrome P450 (CYP) genes are important for the implementation of precision medicine. Genotyping CYP genes using high throughput sequencing data is complicated, e.g., by being highly polymorphic, not to mention the structural variations especially in CYP2D6, CYP2A6, and CYP2B6. Genome graph-based variant detection approaches have been shown to be reliable for genotyping HLA alleles. However, their application to enhancing star allele calling in CYP genes has not been extensively explored. We present StellarPGx, a Nextflow pipeline for accurately genotyping CYP genes by combining genome graph-based variant detection, read coverage information from the original reference-based alignments, and combinatorial diplotype assignments. The implementation of StellarPGx using Nextflow facilitates its portability, reproducibility, and scalability on various user platforms. StellarPGx is currently able to genotype 12 important pharmacogenes belonging to the CYP1, 2, and 3 families. For purposes of validation, we use CYP2D6 as a model gene owing to its high degree of polymorphisms (over 130 star alleles defined to date, including complex structural variants) and clinical importance. We applied StellarPGx and three existing callers to 109 whole genome sequenced samples for which the Genetic Testing Reference Material Coordination Program (GeT-RM) has recently provided consensus truth CYP2D6 diplotypes. StellarPGx had the highest CYP2D6 diplotype concordance (99%) with GeT-RM compared with Cyrius (98%), Aldy (82%), and Stargazer (84%). This exemplifies the high accuracy of StellarPGx and highlights its importance for both research and clinical pharmacogenomics applications. The StellarPGx pipeline is open-source and available from https://github.com/SBIMB/StellarPGx.

An artificial intelligence algorithm for prostate cancer diagnosis in whole slide images of core needle biopsies: a blinded clinical validation and deployment study

BACKGROUND

There is high demand to develop computer-assisted diagnostic tools to evaluate prostate core needle biopsies (CNBs), but little clinical validation and a lack of clinical deployment of such tools. We report here on a blinded clinical validation study and deployment of an artificial intelligence (AI)-based algorithm in a pathology laboratory for routine clinical use to aid prostate diagnosis.

METHODS

An AI-based algorithm was developed using haematoxylin and eosin (H&E)-stained slides of prostate CNBs digitised with a Philips scanner, which were divided into training (1 357 480 image patches from 549 H&E-stained slides) and internal test (2501 H&E-stained slides) datasets. The algorithm provided slide-level scores for probability of cancer, Gleason score 7-10 (vs Gleason score 6 or atypical small acinar proliferation [ASAP]), Gleason pattern 5, and perineural invasion and calculation of cancer percentage present in CNB material. The algorithm was subsequently validated on an external dataset of 100 consecutive cases (1627 H&E-stained slides) digitised on an Aperio AT2 scanner. In addition, the AI tool was implemented in a pathology laboratory within routine clinical workflow as a second read system to review all prostate CNBs. Algorithm performance was assessed with area under the receiver operating characteristic curve (AUC), specificity, and sensitivity, as well as Pearson's correlation coefficient (Pearson's r) for cancer percentage.

FINDINGS

The algorithm achieved an AUC of 0·997 (95% CI 0·995 to 0·998) for cancer detection in the internal test set and 0·991 (0·979 to 1·00) in the external validation set. The AUC for distinguishing between a low-grade (Gleason score 6 or ASAP) and high-grade (Gleason score 7-10) cancer diagnosis was 0·941 (0·905 to 0·977) and the AUC for detecting Gleason pattern 5 was 0·971 (0·943 to 0·998) in the external validation set. Cancer percentage calculated by pathologists and the algorithm showed good agreement (r=0·882, 95% CI 0·834 to 0·915; p<0·0001) with a mean bias of -4·14% (-6·36 to -1·91). The algorithm achieved an AUC of 0·957 (0·930 to 0·985) for perineural invasion. In routine practice, the algorithm was used to assess 11 429 H&E-stained slides pertaining to 941 cases leading to 90 Gleason score 7-10 alerts and 560 cancer alerts. 51 (9%) cancer alerts led to additional cuts or stains being ordered, two (4%) of which led to a third opinion request. We report on the first case of missed cancer that was detected by the algorithm.

INTERPRETATION

This study reports the successful development, external clinical validation, and deployment in clinical practice of an AI-based algorithm to accurately detect, grade, and evaluate clinically relevant findings in digitised slides of prostate CNBs.

FUNDING

Ibex Medical Analytics.

Functional characterisation of the transcriptome from leaf tissue of the fluoroacetate-producing plant, Dichapetalum cymosum, in response to mechanical wounding

Dichapetalum cymosum produces the toxic fluorinated metabolite, fluoroacetate, presumably as a defence mechanism. Given the rarity of fluorinated metabolites in nature, the biosynthetic origin and function of fluoroacetate have been of particular interest. However, the mechanism for fluorination in D. cymosum was never elucidated. More importantly, there is a severe lack in knowledge on a genetic level for fluorometabolite-producing plants, impeding research on the subject. Here, we report on the first transcriptome for D. cymosum and investigate the wound response for insights into fluorometabolite production. Mechanical wounding studies were performed and libraries of the unwounded (control) and wounded (30 and 60 min post wounding) plant were sequenced using the Illumina HiSeq platform. A combined reference assembly generated 77,845 transcripts. Using the SwissProt, TrEMBL, GO, eggNOG, KEGG, Pfam, EC and PlantTFDB databases, a 69% annotation rate was achieved. Differential expression analysis revealed the regulation of 364 genes in response to wounding. The wound responses in D. cymosum included key mechanisms relating to signalling cascades, phytohormone regulation, transcription factors and defence-related secondary metabolites. However, the role of fluoroacetate in inducible wound responses remains unclear. Bacterial fluorinases were searched against the D. cymosum transcriptome but transcripts with homology were not detected suggesting the presence of a potentially different fluorinating enzyme in plants. Nevertheless, the transcriptome produced in this study significantly increases genetic resources available for D. cymosum and will assist with future research into fluorometabolite-producing plants.

Autozygosity influences cardiometabolic disease-associated traits in the AWI-Gen sub-Saharan African study

The analysis of the effects of autozygosity, measured as the change of the mean value of a trait among offspring of genetic relatives, reveals the existence of directional dominance or overdominance. In this study we detect evidence of the effect of autozygosity in 4 out of 13 cardiometabolic disease-associated traits using data from more than 10,000 sub-Saharan African individuals recruited from Ghana, Burkina Faso, Kenya and South Africa. The effect of autozygosity on these phenotypes is found to be sex-related, with inbreeding having a significant decreasing effect in men but a significant increasing effect in women for several traits (body mass index, subcutaneous adipose tissue, low-density lipoproteins and total cholesterol levels). Overall, the effect of inbreeding depression is more intense in men. Differential effects of inbreeding depression are also observed between study sites with different night-light intensity used as proxy for urban development. These results suggest a directional dominant genetic component mediated by environmental interactions and sex-specific differences in genetic architecture for these traits in the Africa Wits-INDEPTH partnership for Genomic Studies (AWI-Gen) cohort.

Gut microbiome profiling of a rural and urban South African cohort reveals biomarkers of a population in lifestyle transition

BACKGROUND

Comparisons of traditional hunter-gatherers and pre-agricultural communities in Africa with urban and suburban Western North American and European cohorts have clearly shown that diet, lifestyle and environment are associated with gut microbiome composition. Yet, little is known about the gut microbiome composition of most communities in the very diverse African continent. South Africa comprises a richly diverse ethnolinguistic population that is experiencing an ongoing epidemiological transition and concurrent spike in the prevalence of obesity, largely attributed to a shift towards more Westernized diets and increasingly inactive lifestyle practices. To characterize the microbiome of African adults living in more mainstream lifestyle settings and investigate associations between the microbiome and obesity, we conducted a pilot study, designed collaboratively with community leaders, in two South African cohorts representative of urban and transitioning rural populations. As the rate of overweight and obesity is particularly high in women, we collected single time-point stool samples from 170 HIV-negative women (51 at Soweto; 119 at Bushbuckridge), performed 16S rRNA gene sequencing on these samples and compared the data to concurrently collected anthropometric data.

RESULTS

We found the overall gut microbiome of our cohorts to be reflective of their ongoing epidemiological transition. Specifically, we find that geographical location was more important for sample clustering than lean/obese status and observed a relatively higher abundance of the Melainabacteria, Vampirovibrio, a predatory bacterium, in Bushbuckridge. Also, Prevotella, despite its generally high prevalence in the cohorts, showed an association with obesity. In comparisons with benchmarked datasets representative of non-Western populations, relatively higher abundance values were observed in our dataset for Barnesiella (log2fold change (FC) = 4.5), Alistipes (log2FC = 3.9), Bacteroides (log2FC = 4.2), Parabacteroides (log2FC = 3.1) and Treponema (log2FC = 1.6), with the exception of Prevotella (log2FC = - 4.7).

CONCLUSIONS

Altogether, this work identifies putative microbial features associated with host health in a historically understudied community undergoing an epidemiological transition. Furthermore, we note the crucial role of community engagement to the success of a study in an African setting, the importance of more population-specific studies to inform targeted interventions as well as present a basic foundation for future research.

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 surveyed¹. 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.

Accuracy and efficiency of an artificial intelligence tool when counting breast mitoses

BACKGROUND

The mitotic count in breast carcinoma is an important prognostic marker. Unfortunately substantial inter- and intra-laboratory variation exists when pathologists manually count mitotic figures. Artificial intelligence (AI) coupled with whole slide imaging offers a potential solution to this problem. The aim of this study was to accordingly critique an AI tool developed to quantify mitotic figures in whole slide images of invasive breast ductal carcinoma.

METHODS

A representative H&E slide from 320 breast invasive ductal carcinoma cases was scanned at 40x magnification. Ten expert pathologists from two academic medical centers labeled mitotic figures in whole slide images to train and validate an AI algorithm to detect and count mitoses. Thereafter, 24 readers of varying expertise were asked to count mitotic figures with and without AI support in 140 high-power fields derived from a separate dataset. Their accuracy and efficiency of performing these tasks were calculated and statistical comparisons performed.

RESULTS

For each experience level the accuracy, precision and sensitivity of counting mitoses by users improved with AI support. There were 21 readers (87.5%) that identified more mitoses using AI support and 13 reviewers (54.2%) that decreased the quantity of falsely flagged mitoses with AI. More time was spent on this task for most participants when not provided with AI support. AI assistance resulted in an overall time savings of 27.8%.

CONCLUSIONS

This study demonstrates that pathology end-users were more accurate and efficient at quantifying mitotic figures in digital images of invasive breast carcinoma with the aid of AI. Higher inter-pathologist agreement with AI assistance suggests that such algorithms can also help standardize practice. Not surprisingly, there is much enthusiasm in pathology regarding the prospect of using AI in routine practice to perform mundane tasks such as counting mitoses.

Interaction of the spike protein RBD from SARS-CoV-2 with ACE2: Similarity with SARS-CoV, hot-spot analysis and effect of the receptor polymorphism

The spread of COVID-19 caused by the SARS-CoV-2 outbreak has been growing since its first identification in December 2019. The publishing of the first SARS-CoV-2 genome made a valuable source of data to study the details about its phylogeny, evolution, and interaction with the host. Protein-protein binding assays have confirmed that Angiotensin-converting enzyme 2 (ACE2) is more likely to be the cell receptor through which the virus invades the host cell. In the present work, we provide an insight into the interaction of the viral spike Receptor Binding Domain (RBD) from different coronavirus isolates with host ACE2 protein. By calculating the binding energy score between RBD and ACE2, we highlighted the putative jump in the affinity from a progenitor form of SARS-CoV-2 to the current virus responsible for COVID-19 outbreak. Our result was consistent with previously reported phylogenetic analysis and corroborates the opinion that the interface segment of the spike protein RBD might be acquired by SARS-CoV-2 via a complex evolutionary process rather than a progressive accumulation of mutations. We also highlighted the relevance of Q493 and P499 amino acid residues of SARS-CoV-2 RBD for binding to human ACE2 and maintaining the stability of the interface. Moreover, we show from the structural analysis that it is unlikely for the interface residues to be the result of genetic engineering. Finally, we studied the impact of eight different variants located at the interaction surface of ACE2, on the complex formation with SARS-CoV-2 RBD. We found that none of them is likely to disrupt the interaction with the viral RBD of SARS-CoV-2.

G6PD variant distribution in sub-Saharan Africa and potential risks of using chloroquine/hydroxychloroquine based treatments for COVID-19

Chloroquine/hydroxychloroquine have been proposed as potential treatments for COVID-19. These drugs have warning labels for use in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency. Analysis of whole-genome sequence data of 458 individuals from sub-Saharan Africa showed significant G6PD variation across the continent. We identified nine variants, of which four are potentially deleterious to G6PD function, and one (rs1050828) that is known to cause G6PD deficiency. We supplemented data for the rs1050828 variant with genotype array data from over 11,000 Africans. Although this variant is common in Africans overall, large allele frequency differences exist between sub-populations. African sub-populations in the same country can show significant differences in allele frequency (e.g. 16.0% in Tsonga vs 0.8% in Xhosa, both in South Africa, ρ=2.4×10 -3 ). The high prevalence of variants in the G6PD gene found in this analysis suggests that it may be a significant interaction factor in clinical trials of chloroquine and hydrochloroquine for treatment of COVID-19 in Africans.

Novel and Known Gene-Smoking Interactions With cIMT Identified as Potential Drivers for Atherosclerosis Risk in West-African Populations of the AWI-Gen Study

Introduction

Atherosclerosis is a key contributor to the burden of cardiovascular diseases (CVDs) and many epidemiological studies have reported on the effect of smoking on carotid intima-media thickness (cIMT) and its subsequent effect on CVD risk. Gene-environment interaction studies have contributed towards understanding some of the missing heritability of genome-wide association studies. Gene-smoking interactions on cIMT have been studied in non-African populations (European, Latino-American, and African American) but no comparable African research has been reported. Our aim was to investigate smoking-SNP interactions on cIMT in two West African populations by genome-wide analysis.

Materials and methods

Only male participants from Burkina Faso (Nanoro = 993) and Ghana (Navrongo = 783) were included, as smoking was extremely rare among women. Phenotype and genotype data underwent stringent QC and genotype imputation was performed using the Sanger African Imputation Panel. Smoking prevalence among men was 13.3% in Nanoro and 42.5% in Navrongo. We analyzed gene-smoking interactions with PLINK after adjusting for covariates: age and 6 PCs (Model 1); age, BMI, blood pressure, fasting glucose, cholesterol levels, MVPA, and 6 PCs (Model 2). All analyses were performed at site level and for the combined data set.

Results

In Nanoro, we identified new gene-smoking interaction variants for cIMT within the previously described RCBTB1 region (rs112017404, rs144170770, and rs4941649) (Model 1: p = 1.35E-07; Model 2: p = 3.08E-08). In the combined sample, two novel intergenic interacting variants were identified, rs1192824 in the regulatory region of TBC1D8 (p = 5.90E-09) and rs77461169 (p = 4.48E-06) located in an upstream region of open chromatin. In silico functional analysis suggests the involvement of genes implicated in biological processes related to cell or biological adhesion and regulatory processes in gene-smoking interactions with cIMT (as evidenced by chromatin interactions and eQTLs).

Discussion

This is the first gene-smoking interaction study for cIMT, as a risk factor for atherosclerosis, in sub-Saharan African populations. In addition to replicating previously known signals for RCBTB1, we identified two novel genomic regions (TBC1D8, near BCHE) involved in this gene-environment interaction.

Organizing and running bioinformatics hackathons within Africa: The H3ABioNet cloud computing experience

The need for portable and reproducible genomics analysis pipelines is growing globally as well as in Africa, especially with the growth of collaborative projects like the Human Health and Heredity in Africa Consortium (H3Africa). The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for portable, reproducible pipelines adapted to heterogeneous computing environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. Building on the network's Standard Operating Procedures (SOPs) for common genomic analyses, H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon in 2016, with the purpose of translating those SOPs into analysis pipelines able to run on heterogeneous computing environments and meeting the needs of H3Africa research projects. This paper describes the preparations for this hackathon and reflects upon the lessons learned about its impact on building the technical and scientific expertise of African researchers. The workflows developed were made publicly available in GitHub repositories and deposited as container images on Quay.io.

Organizing and running bioinformatics hackathons within Africa: The H3ABioNet cloud computing experience

The need for portable and reproducible genomics analysis pipelines is growing globally as well as in Africa, especially with the growth of collaborative projects like the Human Health and Heredity in Africa Consortium (H3Africa). The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for portable, reproducible pipelines adapted to heterogeneous computing environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. Building on the network’s Standard Operating Procedures (SOPs) for common genomic analyses, H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon in 2016, with the purpose of translating those SOPs into analysis pipelines able to run on heterogeneous computing environments and meeting the needs of H3Africa research projects. This paper describes the preparations for this hackathon and reflects upon the lessons learned about its impact on building the technical and scientific expertise of African researchers. The workflows developed were made publicly available in GitHub repositories and deposited as container images on Quay.io.

Classical Cardiovascular Risk Factors and HIV are Associated With Carotid Intima-Media Thickness in Adults From Sub-Saharan Africa: Findings From H3Africa AWI-Gen Study

Background Studies on the determinants of carotid intima-media thickness ( CIMT ), a marker of sub-clinical atherosclerosis, mostly come from white, Asian, and diasporan black populations. We present CIMT data from sub-Saharan Africa, which is experiencing a rising burden of cardiovascular diseases and infectious diseases. Methods and Results The H3 (Human Hereditary and Health) in Africa's AWI-Gen (African-Wits-INDEPTH partnership for Genomic) study is a cross-sectional study conducted in adults aged 40 to 60 years from Burkina Faso, Kenya, Ghana, and South Africa. Cardiovascular disease risk and ultrasonography of the CIMT of right and left common carotids were measured. Multivariable linear and mixed-effect multilevel regression modeling was applied to determine factors related to CIMT. Data included 8872 adults (50.8% men), mean age of 50±6 years with age- and sex-adjusted mean (±SE) CIMT of 640±123μm. Participants from Ghana and Burkina Faso had higher CIMT compared with other sites. Age (β = 6.77, 95%CI [6.34-7.19]), body mass index (17.6[12.5-22.8]), systolic blood pressure (7.52[6.21-8.83]), low-density lipoprotein cholesterol (5.08[2.10-8.06]) and men (10.3[4.75- 15.9]) were associated with higher CIMT. Smoking was associated with higher CIMT in men. High-density lipoprotein cholesterol (-12.2 [-17.9- -6.41]), alcohol consumption (-13.5 [-19.1--7.91]) and HIV (-8.86 [-15.7--2.03]) were inversely associated with CIMT. Conclusions Given the rising prevalence of cardiovascular diseases risk factors in sub-Saharan Africa, atherosclerotic diseases may become a major pan-African epidemic unless preventive measures are taken particularly for prevention of hypertension, obesity, and smoking. HIV -specific studies are needed to fully understand the association between HIV and CIMT in sub-Saharan Africa.

Regional and sex-specific variation in BMI distribution in four sub-Saharan African countries: The H3Africa AWI-Gen study

Background: African populations are characterised by diversity at many levels including: demographic history, genetic ancestry, language, wealth, socio-political landscape, culture and behaviour. Several of these have a profound impact on body fat mass. Obesity, a key risk factor for cardiovascular and metabolic diseases, in the wake of the epidemiological and health transitions across the continent, requires detailed analysis together with other major risk factors.

Objective: To compare regional and sex-specific body mass index (BMI) distributions, using a cross-sectional study design, in adults aged 40–60 years across six study sites in four sub-Saharan African (SSA) countries and to compare the determinants of BMI at each.

Methods: Anthropometric measurements were standardised across sites and BMI calculated. Median BMI and prevalence of underweight, lean, overweight and obesity were compared between the sexes and across sites. Data from multivariable linear regression models for the principal determinants of BMI were summarised from the site-specific studies.

Results: BMI was calculated in 10,702 participants (55% female) and was significantly higher in women than men at nearly all sites. The highest prevalence of obesity was observed at the three South African sites (42.3–66.6% in women and 2.81–17.5% in men) and the lowest in West Africa (1.25–4.22% in women and 1.19–2.20% in men). Across sites, higher socio-economic status and educational level were associated with higher BMI. Being married and increased dietary intake were associated with higher BMI in some communities, whilst smoking and alcohol intake were associated with lower BMI, as was HIV infection in the regions where it was prevalent.

Conclusion: In SSA there is a marked variation in the prevalence of obesity both regionally and between men and women. Our data suggest that the drive for social upliftment within Africa will be associated with rising levels of obesity, which will require the initiation of targeted sex-specific intervention programmes across specific African communities.

Surveying Gut Microbiome Research in Africans: Toward Improved Diversity and Representation

Descriptive and translational investigations into the human gut microbiome (GM) are rapidly expanding; however, studies are largely restricted to industrialized populations in the USA and Europe. Little is known about microbial variability and its implications for health and disease in other parts of the world. Populations in Africa are particularly underrepresented. What limited research has been performed has focused on a few subject domains, including the impact of long-term lifestyle and dietary factors on GM ecology, its maturation during infancy, and the interrelationships between the microbiome, infectious disease, and undernutrition. Recently, international consortia have laid the groundwork for large-scale genomics and microbiome studies on the continent, with a particular interest in the epidemiologic transition to noncommunicable disease. Here, we survey the current landscape of GM scholarship in Africa and propose actionable recommendations to improve research capacity and output.

Developing reproducible bioinformatics analysis workflows for heterogeneous computing environments to support African genomics

Background

The Pan-African bioinformatics network, H3ABioNet, comprises 27 research institutions in 17 African countries. H3ABioNet is part of the Human Health and Heredity in Africa program (H3Africa), an African-led research consortium funded by the US National Institutes of Health and the UK Wellcome Trust, aimed at using genomics to study and improve the health of Africans. A key role of H3ABioNet is to support H3Africa projects by building bioinformatics infrastructure such as portable and reproducible bioinformatics workflows for use on heterogeneous African computing environments. Processing and analysis of genomic data is an example of a big data application requiring complex interdependent data analysis workflows. Such bioinformatics workflows take the primary and secondary input data through several computationally-intensive processing steps using different software packages, where some of the outputs form inputs for other steps. Implementing scalable, reproducible, portable and easy-to-use workflows is particularly challenging.

Results

H3ABioNet has built four workflows to support (1) the calling of variants from high-throughput sequencing data; (2) the analysis of microbial populations from 16S rDNA sequence data; (3) genotyping and genome-wide association studies; and (4) single nucleotide polymorphism imputation. A week-long hackathon was organized in August 2016 with participants from six African bioinformatics groups, and US and European collaborators. Two of the workflows are built using the Common Workflow Language framework (CWL) and two using Nextflow. All the workflows are containerized for improved portability and reproducibility using Docker, and are publicly available for use by members of the H3Africa consortium and the international research community.

Conclusion

The H3ABioNet workflows have been implemented in view of offering ease of use for the end user and high levels of reproducibility and portability, all while following modern state of the art bioinformatics data processing protocols. The H3ABioNet workflows will service the H3Africa consortium projects and are currently in use. All four workflows are also publicly available for research scientists worldwide to use and adapt for their respective needs. The H3ABioNet workflows will help develop bioinformatics capacity and assist genomics research within Africa and serve to increase the scientific output of H3Africa and its Pan-African Bioinformatics Network.

African genetic diversity provides novel insights into evolutionary history and local adaptations

Genetic variation and susceptibility to disease are shaped by human demographic history and adaptation. We can now study the genomes of extant Africans and uncover traces of population migration, admixture, assimilation and selection by applying sophisticated computational algorithms. There are four major ethnolinguistic divisions among present day Africans: Hunter-gatherer populations in southern and central Africa; Nilo-Saharan speakers from north and northeast Africa; Afro-Asiatic speakers from north and east Africa; and Niger-Congo speakers who are the predominant ethnolinguistic group spread across most of sub-Saharan Africa. The enormous ethnolinguistic diversity in sub-Saharan African populations is largely paralleled by extensive genetic diversity and until a decade ago, little was known about detailed origins and divergence of these groups. Results from large-scale population genetic studies, and more recently whole genome sequence data, are unravelling the critical role of events like migration and admixture and environmental factors including diet, infectious diseases and climatic conditions in shaping current population diversity. It is now possible to start providing quantitative estimates of divergence times, population size and dynamic processes that have affected populations and their genetic risk for disease. Finally, the availability of ancient genomes from Africa provides historical insights of unprecedented depth. In this review, we highlight some key interpretations that have emerged from recent African genome studies.

Genetic variants in SEC16B are associated with body composition in black South Africans

Objective

The latest genome-wide association studies of obesity-related traits have identified several genetic loci contributing to body composition (BC). These findings have not been robustly replicated in African populations, therefore, this study aimed to assess whether European BC-associated gene loci played a similar role in a South African black population.

Methods

A replication and fine-mapping study was performed in participants from the Birth to Twenty cohort (N = 1,926) using the Metabochip. Measurements included body mass index (BMI), waist and hip circumference, waist-to-hip ratio (WHR), total fat mass, total lean mass and percentage fat mass (PFM).

Results

SNPs in several gene loci, including SEC16B (P_adj <  9.48 × 10⁻⁷), NEGR1 (P_adj < 1.64 × 10⁻⁶), FTO (P_adj < 2.91 × 10⁻⁵), TMEM18 (P_adj < 2.27 × 10⁻⁵), and WARS2(P_adj < 3.25 × 10⁻⁵) were similarly associated (albeit not at array-wide signficance (P ≤ 6.7 × 10⁻⁷) with various phenotypes including fat mass, PFM, WHR linked to BC in this African cohort, however the associations were driven by different sentinel SNPs. More importantly, DXA-derived BC measures revealed stronger genetic associations than simple anthropometric measures. Association signals generated in this study were shared by European and African populations, as well as unique to this African cohort. Moreover, sophisticated estimates like DXA measures enabled an enhanced characterisation of genetic associations for BC traits.

Conclusion

Results from this study suggest that in-depth genomic studies in larger African cohorts may reveal novel SNPs for body composition and adiposity, which will provide greater insight into the aetiology of obesity.

Assessing runs of Homozygosity: a comparison of SNP Array and whole genome sequence low coverage data

BACKGROUND

Runs of Homozygosity (ROH) are genomic regions where identical haplotypes are inherited from each parent. Since their first detection due to technological advances in the late 1990s, ROHs have been shedding light on human population history and deciphering the genetic basis of monogenic and complex traits and diseases. ROH studies have predominantly exploited SNP array data, but are gradually moving to whole genome sequence (WGS) data as it becomes available. WGS data, covering more genetic variability, can add value to ROH studies, but require additional considerations during analysis.

RESULTS

Using SNP array and low coverage WGS data from 1885 individuals from 20 world populations, our aims were to compare ROH from the two datasets and to establish software conditions to get comparable results, thus providing guidelines for combining disparate datasets in joint ROH analyses. By allowing heterozygous SNPs per window, using the PLINK homozygosity function and non-parametric analysis, we were able to obtain non-significant differences in number ROH, mean ROH size and total sum of ROH between data sets using the different technologies for almost all populations.

CONCLUSIONS

By allowing 3 heterozygous SNPs per ROH when dealing with WGS low coverage data, it is possible to establish meaningful comparisons between data using SNP array and WGS low coverage technologies.

Genomic and environmental risk factors for cardiometabolic diseases in Africa: methods used for Phase 1 of the AWI-Gen population cross-sectional study

There is an alarming tide of cardiovascular and metabolic disease (CMD) sweeping across Africa. This may be a result of an increasingly urbanized lifestyle characterized by the growing consumption of processed and calorie-dense food, combined with physical inactivity and more sedentary behaviour. While the link between lifestyle and public health has been extensively studied in Caucasian and African American populations, few studies have been conducted in Africa. This paper describes the detailed methods for Phase 1 of the AWI-Gen study that were used to capture phenotype data and assess the associated risk factors and end points for CMD in persons over the age of 40 years in sub-Saharan Africa (SSA). We developed a population-based cross-sectional study of disease burden and phenotype in Africans, across six centres in SSA. These centres are in West Africa (Nanoro, Burkina Faso, and Navrongo, Ghana), in East Africa (Nairobi, Kenya) and in South Africa (Agincourt, Dikgale and Soweto). A total of 10,702 individuals between the ages of 40 and 60 years were recruited into the study across the six centres, plus an additional 1021 participants over the age of 60 years from the Agincourt centre. We collected socio-demographic, anthropometric, medical history, diet, physical activity, fat distribution and alcohol/tobacco consumption data from participants. Blood samples were collected for disease-related biomarker assays, and genomic DNA extraction for genome-wide association studies. Urine samples were collected to assess kidney function. The study provides base-line data for the development of a series of cohorts with a second wave of data collection in Phase 2 of the study. These data will provide valuable insights into the genetic and environmental influences on CMD on the African continent.