Uganda Genome Resource: A rich research database for genomic studies of communicable and non-communicable diseases in Africa

Study Population – The General Population Cohort (GPC)

The GPC is a population-based study of approximately 22,000 individuals residing in 25 neighbouring villages in the Kyamulibwa sub-county, Kalungu district in rural south-western Uganda. The study was founded in 1989 by the Medical Research Council UK (MRC UK) in collaboration with the Uganda Virus Research Institute (UVRI) to study the epidemiology of HIV in a general population. The GPC population was initially recruited and assessed through annual house-to-house census and survey rounds until 2012, when biannual surveys commenced. Since its establishment, 26 rounds of survey and 29 rounds of census have been undertaken. Before any survey procedures are carried out, written informed consent is obtained from participants on the use of their clinical records for research purposes and sample storage for future use (Asiki et al., 2013). Data collected includes serological, demographic, and medical information from participants. Information regarding mortality, fertility, sexual behaviour migration, HIV infection perception are routinely collated.

The GPC Round 22 study of 2011 focused on the genetics and epidemiology of communicable and non-communicable disease. The survey round which was used to establish the Ugandan Genome Resource (UGR) consisted of five main stages, including mobilization (recruitment and consenting), mapping, census, survey, and results feedback and clinical follow-up. The specific objectives of this survey then were:

• To create a one-of-a-kind study for expanding on a large-scale prospective cohort research in an African population to evaluate a wide range of health indices—and to lay the platform for longer-term investigations.

• Using population, genetic and epidemiological techniques to provide aetiological insights into variance in cardio-metabolic and infectious risk factors.

• To help develop public health policies in other African countries by informing health policy and public health programs aimed at addressing the rise in NCDs in Uganda.

The cohort continues to offer as a valuable platform for investigating the relationship between communicable illnesses and NCDs in a regular annual survey of GPC.

Genotype generation, quality control and Imputation

The 2.5M Illumina chip array was used to genotype nearly 5,000 Ugandans at the Wellcome Trust Sanger Institute. Gurdasani and Fatumo et al 2019 have presented the quality control steps (Gurdasani et al., 2019). In summary, we used a strict quality control process to perform a series of steps in a logical order to eliminate a total of 39,368 autosomal markers that failed to meet the quality metrics for SNP call rate (>97 percent, 25,037 SNPs) and HWE (p<1×10⁻⁸, 14,331 SNPs). During sample QC, a total of 91 samples were eliminated because they failed the quality standards for sample call rate (>97%) or heterozygosity (HO=0.209333±0.007416 matching to the mean±3SD), or the sex extrapolated from the X-chromosome did not correspond to the reported sex. Three further samples were eliminated due to high relatedness (IBD>0.90). There were no samples that were classified as outliers in terms of population or ancestry. 2,230,258 autosomal markers and 4778 samples that met the stated threshold were subjected to further analysis. We carried out SNP Phasing with the aid of SHAPEIT2 (Delaneau, Coulonges and Zagury, 2008) using default settings, then imputation was done with IMPUTE2 (Howie, Donnelly and Marchini, 2009). All samples were imputed with a combined reference that was created by combining the UG2G sequence resource (n = 2,000, whole genome sequence data from the African Genome Variation Project (n = 320),), and the 1000 Genomes phase 3 project (n = 2,504). The principal components analysis plot for the GPC participants (n=4,778) was published (Gurdasani et al., 2019) and is shown here in supplementary material figure S1.

Uganda 2000 Genomes (UG2G)

The entire genomes of over 2,000 Ugandans from nine ethnolinguistic groups were sequenced using the Illumina HiSeq 2000 with 75bp paired end reads at low coverage, with an average coverage of 4x for each sample. 343 of these samples overlapped with people who had already been genotyped. An automated quality control process was used to bring down the data files that needs manual processing to ascertain the quality of BAM files produced. This method was based on the one developed for the UK10K project (Walter et al., 2015) which used a set of algorithmically derived standards to determine summary data computed from the input BAMs. Any line that fell below the “fail” standard for any of the metrics was deleted; while lines falling below the warn standard for any of the scores were manually investigated; and any line that passes any of these scores was given a status of “pass”. Overall, we deleted fourteen samples from the study. Full detailed on the quality control and how we computed the summary data has be described in Gurdasani et al., 2019.

Merging of Sequenced and Genotyped Data

We integrated sequenced and imputed genotyped data to produce an aggregated dataset to boost power for discovery in a genome-wide association studies. Because cryptic and family relatedness persisted across sequenced and genotyped data, we produced an aggregated dataset for analysis instead of separately meta-analysing the data, because data would be correlated rather than independent. As a result, conclusions from mixed model analysis that explicitly model this relationship are more likely to be true. We examined and deleted any consistent discrepancies between sequences and imputed genotype data after merging the two datasets. This was done by performing principal component analysis on the dataset to see if there was any distinction by data modality (imputed genotype data vs. sequenced data) among the 343 people who had their genotypes and sequences done in duplicate. On PCA, we noticed a strong separation of genotype imputed and sequence data points. For these 343 samples, we tested alternative concordance criteria between sequencing and imputed genotype data, screening out SNPs with a concordance of 0.80 and 0.90 the dataset. We discovered that to eliminate systematic effects detected between genotyping array and sequence data on PCA, a minimum concordance criterion of 0.90 was necessary.

There were no systematic changes between sequenced and genotyped data in PCAs after excluding 904,283 SNPs that exhibited 90 percent consonance in genotypes between the sequence and imputed genotype data. We examined the top ten PCs to confirm that systematic variations in the genomic data did not constitute an important axis of variation.

Phenotype and laboratory measurement

During survey round 22 which was conducted in 2011, several phenotypes based on clinical and physical examinations, laboratory tests, and self-reported questionnaires were collected from the respondents (Table 1) and these respondents who are still known to be alive and have not moved out of the GPC have been followed every year since then. A blood specimen was analysed for non-fasting blood lipids, blood cell traits (mean cell haemoglobin, red cell count, white cell count, mean cell haemoglobin concentration, haemoglobin, packed cell volume, mean cell volume and platelet), glycaemic characteristics, renal function, infectious biomarkers (HIV, hepatitis B and C). Basic demographics such as age, sex, marital status, and education level; anthropometrics such as BMI, weight, waist-to-hip ratio, height; blood pressure measurements; and lifestyle information such as smoking status, physical activities, and diet; sexual health & reproductive behaviour; sex education, condom use, pregnancy & outcome, and number of offspring were also collected (Table 1).

Genetic discoveries and polygenic prediction in UGR

A case in point in the use of our rich African genomics and phenotypic data, we undertook GWAS in 34 cardiometabolic traits including lipid, anthropometry traits, blood cell indices, HbA1c and reported novel loci associated with anthropometric, haematological, lipid, and glycaemic traits (Gurdasani et al., 2019). In another study (Fatumo et al 2020 - eGFR), we reported the first ever GWAS of kidney function in continental Africa. Leveraging clinical relatedness and correlations among phenotypes, we explored the power of multivariate GWAS to identify genetic risk factors implicating pleiotropic effects in blood cell traits (Fatumo et al., 2019; Soremekun et al., 2021), body shape (Nakabuye et al., 2022) and liver function. Recently, we showed that genetic risk score derived from data of African American individuals enhance polygenic prediction of lipid traits and T2DM in Sub-Sahara African, but prediction varied greatly between another dataset from South Africa and our East African genomic data (Chikowore et al., 2022, Kamiza et al., 2022). We have also demonstrated the Mendelian randomisation evidence of relation between lipid trait and T2DM (Soremekun et al., 2022), metabolic traits and stroke (Fatumo et al., 2021). Collectively, our studies show a need for improved representation of Africans in genomic studies and ensuring the generalisation of findings for genomic medicine. This is further supported by findings from another study as well (Martin et al., 2017). The UGR data has also been used to create a genotype imputation reference panel using UG2G available from the Sanger Imputation Service (imputation.sanger.ac.uk).

Contribution to collaborative studies

We contribute to global genetic studies through partnerships and consortia, such as the African Partnership for Chronic Disease Research (APCDR), an international network of research groups that collaborate to support and promote collaborative chronic disease research across Africa. An initiative created in response to the changing distribution of communicable diseases and the rising burden of noncommunicable diseases, as well as the recognition that low- and middle-income countries (LMICs), including those in Sub-Saharan Africa, will need to expand their health-care capacities to effectively respond to these epidemiological transitions.

We combine research expertise with three other MRC Units (MRC Integrative Epidemiology Unit, MRC Population Health Research Unit, and MRC Unit for Lifelong Health and Ageing) to investigate the potential to use Mendelian randomization (MR) to assess the generalisability of existing drugs (e.g., statins, anti-diabetics, and anti-hypertensives) and identify the potential to tailor drugs with pilot studies focusing on established pharmacological targets to specific subpopulations (e.g. CETP, HMGCR) and to see how changes in genetic architecture affect efficacy estimates in different groups.

We are part of the CARDINAL (CARDiometabolic Disorders IN African-ancestry PopuLations) which is study site of an NIH-funded Polygenic Risk Methods in Diverse Populations (PRIMED) Consortium https://primedconsortium.org/. CARDINAL (Adebamowo, C.A. et al., 2022) aims to integrate phenotype and genomic datasets from 50,000 African individuals from seven cohort studies and evaluate PRSs to develop a novel method that considers ancestry-specific genomic regions to improve PRS prediction in populations with genetic substructure.

Furthermore, we recently provided GWAS data to the Meta-Analyses of Glucose and Insulin-related Variables Consortium (MAGIC) in order to find additional loci that influence glycaemic and metabolic traits (Chen, J. et al 2021). We are aiming for opportunities to contribute key phenotypes such as lipids, blood cell traits, kidney, etc to other consortia. For GBMI we will contribute all phenotype in Table 1 when require, including opportunity to measure not previously collected phenotype using resources in our organisation. We believe that team science allows scientists to make the most progress toward breakthrough discoveries that benefit human health.

Future directions

The GPC is an active cohort of more than 22,000 participants. Genotype and sequence data is available for 6,657 respondents (N=5,000, 2,000 and 343 for genotype, sequence and overlapping samples respectively). We hope to genotype more samples to add on this resource. We also hope to sequence more samples at higher coverage in order to provide a reference panel with increased genome coverage. We also hope to extend our research into proteomics, metabolomics and single cell genomics in order to gain insights into the different mechanisms and pathways that could be implicated in different disease processes.

Data access and sharing of the UGR data

Request for resources and information should be directed to UGR’s Data Access Committee (DAC) via the Lead Contact, Dr. Segun Fatumo (segun.fatumo{at}mrcuganda.org; segun.fatumo{at}lshtm.ac.uk). UGR’s individual level data, genotype and sequence data are available under managed access to researchers. Requests for access will be granted for all research consistent with the consent provided by participants. This would include any research in the context of health and disease, that does not involve identifying the participants in any way.

The array and low and high depth sequence data have been deposited at the European Genome-phenome Archive (EGA, https://www.ebi.ac.uk/ega/, accession numbers EGAS00001001558/EGAD00010000965, EGAS00001000545/EGAD00001001639 and EGAS00001000545/EGAD00001005346 respectively. Requests for access to data may be directed to segun.fatumo{at}mrcuganda.org. Applications are reviewed by data access committee (DAC) and access is granted if the request is consistent with the consent provided by participants. The data producers may be consulted by the DAC to evaluate potential ethical conflicts. Requestors also sign an agreement which governs the terms on which access to data is granted.

However, full GWAS summary statistics of UGR is freely available on GWAS catalog https://www.ebi.ac.uk/gwas/ with study accssion numbers: GCST009041 (Eosinophil counts), GCST009042 (Total cholesterol levels), GCST009043 (LDL cholesterol levels), GCST0090414 (HDL cholesterol levels), GCST009045 (Triglyceride levels), GCST009046(Aspartate aminotransferase levels), GCST009047 (Alanine aminotransferase levels), GCST009048 (Serum albumin levels), GCST009049 (Serum alkaline phosphatase levels), GCST009050 (Gamma glutamyl transferase levels), GCST009051 (Bilirubin levels), GCST009052 (Diastolic blood pressure) GCST009053 (Systolic blood pressure), GCST009054 (Hemoglobin A1c levels), GCST009055 (Height), GCST009056 (Weight), GCST009057 (Body mass index), GCST009058 (Waist circumference), GCST009059 (Hip circumference), GCST009060 (Waist-hip ratio), GCST009061 (White blood cell count), GCST009062 (Red blood cell count), GCST009063 (mean corpuscular hemoglobin), GCST009064 (mean corpuscular hemoglobin concentration), GCST009065 (mean corpuscular volume), GCST009032 (red blood cell distribution width), GCST009033 (hematocrit), GCST009034 (hemoglobin measurement), GCST009035 (mean platelet volume), GCST009036 (platelet count), GCST009037 (lymphocyte count), GCST009038 (monocyte count), GCST009039 (basophil count), GCST009040 (neutrophil count)