A Sex-Stratified Genome-Wide Association Study of Tuberculosis Using a Multi-Ethnic Genotyping Array

History of tuberculosis disease is associated with genetic regulatory variation in Peruvians

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Multi-ancestry meta-analysis of host genetic susceptibility to tuberculosis identifies shared genetic architecture

The heritability of susceptibility to tuberculosis (TB) disease has been well recognized. Over 100 genes have been studied as candidates for TB susceptibility, and several variants were identified by genome-wide association studies (GWAS), but few replicate. We established the International Tuberculosis Host Genetics Consortium to perform a multi-ancestry meta-analysis of GWAS, including 14,153 cases and 19,536 controls of African, Asian, and European ancestry. Our analyses demonstrate a substantial degree of heritability (pooled polygenic h2 = 26.3%, 95% CI 23.7-29.0%) for susceptibility to TB that is shared across ancestries, highlighting an important host genetic influence on disease. We identified one global host genetic correlate for TB at genome-wide significance (p<5 × 10-8) in the human leukocyte antigen (HLA)-II region (rs28383206, p-value=5.2 × 10-9) but failed to replicate variants previously associated with TB susceptibility. These data demonstrate the complex shared genetic architecture of susceptibility to TB and the importance of large-scale GWAS analysis across multiple ancestries experiencing different levels of infection pressure.

Sex and Gender Differences in Tuberculosis Pathogenesis and Treatment Outcomes

Tuberculosis remains a daunting public health concern in many countries of the world. A consistent observation in the global epidemiology of tuberculosis is an excess of cases of active pulmonary tuberculosis among males compared with females. Data from both humans and animals also suggest that males are more susceptible than females to develop active pulmonary disease. Similarly, male sex has been associated with poor treatment outcomes. Despite this growing body of evidence, little is known about the mechanisms driving sex bias in tuberculosis disease. Two dominant hypotheses have been proposed to explain the predominance of active pulmonary tuberculosis among males. The first is based on the contribution of biological factors, such as sex hormones and genetic factors, on host immunity during tuberculosis. The second is focused on non-biological factors such as smoking, professional exposure, and health-seeking behaviors, known to be influenced by gender. In this chapter, we review the literature regarding these two prevailing hypotheses by presenting human but also experimental animal studies. In addition, we presented studies aiming at examining the impact of sex and gender on other clinical forms of tuberculosis such as latent tuberculosis infection and extrapulmonary tuberculosis, which both appear to have their own specificities in relation to sex. We also highlighted potential intersections between sex and gender in the context of tuberculosis and shared future directions that could guide in elucidating mechanisms of sex-based differences in tuberculosis pathogenesis and treatment outcomes.

A systematic review of analytical methods used in genetic association analysis of the X-chromosome

Genetic association studies have been very successful at elucidating the genetic background of many complex diseases/traits. However, the X-chromosome is often neglected in these studies because of technical difficulties and the fact that most tools only utilize genetic data from autosomes. In this review, we aim to provide an overview of different practical approaches that are followed to incorporate the X-chromosome in association analysis, such as Genome-Wide Association Studies and Expression Quantitative Trait Loci Analysis. In general, the choice of which test statistics is most appropriate will depend on three main criteria: (1) the underlying X-inactivation model, (2) if Hardy-Weinberg equilibrium holds and sex-specific allele frequencies are expected and (3) whether adjustment for confounding variables is required. All in all, it is recommended that a combination of different association tests should be used for the analysis of X-chromosome.

Monocyte Transcriptional Responses to Mycobacterium tuberculosis Associate with Resistance to Tuberculin Skin Test and Interferon Gamma Release Assay Conversion

Heavy exposure to Mycobacterium tuberculosis, the etiologic agent of tuberculosis (TB) and among the top infectious killers worldwide, results in infection that is cleared, contained, or progresses to disease. Some heavily exposed tuberculosis contacts show no evidence of infection using the tuberculin skin test (TST) and interferon gamma release assay (IGRA); yet the mechanisms underlying this "resister" (RSTR) phenotype are unclear. To identify transcriptional responses that distinguish RSTR monocytes, we performed transcriptome sequencing (RNA-seq) on monocytes isolated from heavily exposed household contacts in Uganda and gold miners in South Africa after ex vivo M. tuberculosis infection. Gene set enrichment analysis (GSEA) revealed several gene pathways that were consistently enriched in response to M. tuberculosis among RSTR subjects compared to controls with positive TST/IGRA testing (latent TB infection [LTBI]) across Uganda and South Africa. The most significantly enriched gene set in which expression was increased in RSTR relative to LTBI M. tuberculosis-infected monocytes was the tumor necrosis factor alpha (TNF-α) signaling pathway whose core enrichment (leading edge) substantially overlapped across RSTR populations. These leading-edge genes included candidate resistance genes (ABCA1 and DUSP2) with significantly increased expression among Uganda RSTRs (false-discovery rate [FDR], <0.1). The distinct monocyte transcriptional response to M. tuberculosis among RSTR subjects, including increased expression of the TNF signaling pathway, highlights genes and inflammatory pathways that may mediate resistance to TST/IGRA conversion and provides therapeutic targets to enhance host restriction of M. tuberculosis intracellular infection. IMPORTANCE After heavy M. tuberculosis exposure, the events that determine why some individuals resist TST/IGRA conversion are poorly defined. Enrichment of the TNF signaling gene set among RSTR monocytes from multiple distinct cohorts suggests an important role for the monocyte TNF response in determining this alternative immune outcome. These TNF responses to M. tuberculosis among RSTRs may contribute to antimicrobial programs that result in early clearance or the priming of alternative (gamma interferon-independent) cellular responses.

Major sex differences in allele frequencies for X chromosomal variants in both the 1000 Genomes Project and gnomAD

An unexpectedly high proportion of SNPs on the X chromosome in the 1000 Genomes Project phase 3 data were identified with significant sex differences in minor allele frequencies (sdMAF). sdMAF persisted for many of these SNPs in the recently released high coverage whole genome sequence of the 1000 Genomes Project that was aligned to GRCh38, and it was consistent between the five super-populations. Among the 245,825 common (MAF>5%) biallelic X-chromosomal SNPs in the phase 3 data presumed to be of high quality, 2,039 have genome-wide significant sdMAF (p-value <5e-8). sdMAF varied by location: non-pseudo-autosomal region (NPR) = 0.83%, pseudo-autosomal regions (PAR1) = 0.29%, PAR2 = 13.1%, and X-transposed region (XTR)/PAR3 = 0.85% of SNPs had sdMAF, and they were clustered at the NPR-PAR boundaries, among others. sdMAF at the NPR-PAR boundaries are biologically expected due to sex-linkage, but have generally been ignored in association studies. For comparison, similar analyses found only 6, 1 and 0 SNPs with significant sdMAF on chromosomes 1, 7 and 22, respectively. Similar sdMAF results for the X chromosome were obtained from the high coverage whole genome sequence data from gnomAD V 3.1.2 for both the non-Finnish European and African/African American samples. Future X chromosome analyses need to take sdMAF into account.

The human X chromosome contains over 800 genes and is the 8^th largest human chromosome. Genome-wide associations studies have generally failed to examine variants on the X chromosome for association with diseases and traits, partly due to complexities of the data analysis, and challenges with genotype imputation. We examined X chromosomal variants from the 1000 Genomes Project for sex differences in allele frequency and found that many variants showed significant differences. These variants cluster at the centromeric parts of the pseudoautosomal regions 1 and 2, as well as the putative pseudo-autosomal region 3 (also termed X-transposed region). This pattern was observed in high coverage whole genome sequence data from the same subjects that was aligned to GRCh38, suggesting that is not an artefact of low coverage sequencing or problems specific to GRCh37. In addition, we replicated this phenomenon in high coverage whole genome sequence aligned to GRCh38 from the gnomAD database in both the non-Finnish European and African/African American samples. These findings have implications for the analysis of X chromosomal variants for disease and trait associations.

Deciphering Genetic Susceptibility to Tuberculous Meningitis

Tuberculous meningitis (TBM) is the most severe form of extrapulmonary tuberculosis (TB) that arises when a caseating meningeal granuloma discharges its contents into the subarachnoid space. It accounts for ~1% of all disease caused by Mycobacterium tuberculosis and the age of peak incidence is from 2-4 years. The exact pathogenesis of TBM is still not fully understood and the mechanism(s) by which the bacilli initially invade the blood-brain-barrier are still to be elucidated. This study investigated the involvement of the host genome in TBM susceptibility, by considering common variants (minor allele frequency (MAF) >5%) using microarray genotyping and rare variants (MAF <1%) via exome sequencing. A total of 123 TBM cases, 400 pulmonary TB (pTB) cases and 477 healthy controls were genotyped on the MEGA array. A genome-wide association study (GWAS) comparing 114 TBM cases to 395 healthy controls showed no association with TBM susceptibility. A second analysis comparing 114 TBM cases to 382 pTB cases was conducted to investigate variants associated with different TB phenotypes. No significant associations were found with progression from pTB to TBM. Ten TBM cases and 10 healthy controls were exome sequenced. Gene set association tests SKAT-O and SKAT Common Rare were used to assess the association of rare SNPs and the cumulative effect of both common and rare SNPs with susceptibility to TBM, respectively. Ingenuity Pathway Analysis (IPA) of the top-hits of the SKAT-O analysis showed that NOD2 and CYP4F2 are both important in TBM pathogenesis and highlighted these as targets for future study. For the SKAT Common Rare analysis Centriolar Coiled-Coil Protein 110 (CCP110) was nominally associated (p = 5.89x10^-6) with TBM susceptibility. In addition, several top-hit genes ascribed to the development of the central nervous system (CNS) and innate immune system regulation were identified. Exome sequencing and GWAS of our TBM cohort has identified a single previously undescribed association of CCP110 with TBM susceptibility. These results advance our understanding of TBM in terms of both variants and genes that influence susceptibility. In addition, several candidate genes involved in innate immunity have been identified for further genotypic and functional investigation.

Local Ancestry Adjusted Allelic Association Analysis Robustly Captures Tuberculosis Susceptibility Loci

Pulmonary tuberculosis (TB), caused by Mycobacterium tuberculosis, is a complex disease. The risk of developing active TB is in part determined by host genetic factors. Most genetic studies investigating TB susceptibility fail to replicate association signals particularly across diverse populations. South African populations arose because of multi-wave genetic admixture from the indigenous KhoeSan, Bantu-speaking Africans, Europeans, Southeast Asian-and East Asian populations. This has led to complex genetic admixture with heterogenous patterns of linkage disequilibrium and associated traits. As a result, precise estimation of both global and local ancestry is required to prevent both false positive and false-negative associations. Here, 820 individuals from South Africa were genotyped on the SNP-dense Illumina Multi-Ethnic Genotyping Array (∼1.7M SNPs) followed by local and global ancestry inference using RFMix. Local ancestry adjusted allelic association (LAAA) models were utilized owing to the extensive genetic heterogeneity present in this population. Hence, an interaction term, comprising the identification of the minor allele that corresponds to the ancestry present at the specific locus under investigation, was included as a covariate. One SNP (rs28647531) located on chromosome 4q22 was significantly associated with TB susceptibility and displayed a SNP minor allelic effect (G allele, frequency = 0.204) whilst correcting for local ancestry for Bantu-speaking African ancestry (p-value = 5.518 × 10^-7; OR = 3.065; SE = 0.224). Although no other variants passed the significant threshold, clear differences were observed between the lead variants identified for each ancestry. Furthermore, the LAAA model robustly captured the source of association signals in multi-way admixed individuals from South Africa and allowed the identification of ancestry-specific disease risk alleles associated with TB susceptibility that have previously been missed.

The role of the X chromosome in infectious diseases

Many infectious diseases in humans present with a sex bias. This bias arises from a combination of environmental factors, hormones and genetics. In this study, we review the contribution of the X chromosome to the genetic factor associated with infectious diseases. First, we give an overview of the X-linked genes that have been described in the context of infectious diseases and group them in four main pathways that seem to be dysregulated in infectious diseases: nuclear factor kappa-B, interleukin 2 and interferon γ cascade, toll-like receptors and programmed death ligand 1. Then, we review the infectious disease associations in existing genome-wide association studies (GWAS) from the GWAS Catalog and the Pan-UK Biobank, describing the main associations and their possible implications for the disease. Finally, we highlight the importance of including the X chromosome in GWAS analysis and the importance of sex-specific analysis.

A case-control study on correlation between the single nucleotide polymorphism of CLEC4E and the susceptibility to tuberculosis among Han people in Western China

Background

Tuberculosis (TB) is one of the leading causes of morbidity and mortality in Western China. Preclinical studies have suggested the protective effect of the C-type lectin receptor of family 4 member E (CLEC4E) from TB. Herein, we investigated the association between CLEC4E gene variants and TB susceptibility in a western Chinese Han population.

Methods

We genotyped four single nucleotide polymorphisms (SNPs) rs10841856, rs10770847, rs10770855 and rs4480590 in the CLEC4E gene using the improved multiplex ligation detection reaction (iMLDR) assay in 900 TB cases and 1534 healthy controls.

Results

After stratifying the whole data by sex, it was found that males exhibited mutant allele G of rs10841856 was more strongly associated with increased TB risk after Bonferroni correction (OR = 1.334, 95% CI: 1.142–1.560; P < 0.001 after adjusting for age; p = 0.001 after Bonferroni correction). The genetic model analysis found that rs10841856 was associated with the increased risk of TB among males under the dominant model (OR = 1.557, 95% CI = 1.228–1.984, P < 0.001 after adjusting for age, P < 0.001 after Bonferroni correction). Bioinformatics analysis suggested that rs10841856 might fall in putative functional regions and might be the expression quantitative trait loci (eQTL) for CLEC4E and long noncoding RNA RP11-561P12.5.

Conclusions

Our study revealed that rs10841856 in the CLEC4E gene might be related to increased TB risk, especially the dominant genetic model among male Han individuals from Western China

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06448-2.

A multi-phenotype genome-wide association study of clades causing tuberculosis in a Ghanaian- and South African cohort

Despite decades of research and advancements in diagnostics and treatment, tuberculosis remains a major public health concern. New computational methods are needed to interrogate the intersection of host- and bacterial genomes. Paired host genotype datum and infecting bacterial isolate information were analysed for associations using a multinomial logistic regression framework implemented in SNPTest. A cohort of 853 admixed South African participants and a Ghanaian cohort of 1359 participants were included. Two directly genotyped variants, namely rs529920 and rs41472447, were identified in the Ghanaian cohort as being statistically significantly associated with risk for infection with strains of different members of the MTBC. Thus, a multinomial logistic regression using paired host-pathogen data may prove valuable for investigating the complex relationships driving infectious disease.

CLEC4E (Mincle) genetic variation associates with pulmonary tuberculosis in Guinea-Bissau (West Africa)

Tuberculosis (TB) is the leading cause of death from a single infectious agent. According to the WHO, 85% of cases in 2018 were pulmonary tuberculosis (PTB), making it the most prevalent form of the disease. Although the bacillus responsible for disease, Mycobacterium tuberculosis (MTB), is estimated to infect 1.7 billion people worldwide, only a small portion of those infected (5-10%) will transition into active TB. Because such a small fraction of infected people develop active disease, we hypothesized that underlying host genetic variation associates with developing active pulmonary disease. Variation in CLEC4E has been of interest in previous association studies showing either no effect or protection from PTB. For our study we assessed 60 SNPs in 11 immune genes, including CLEC4E, using a case-control study from Guinea-Bissau. The 289 cases and 322 controls differed in age, sex, and ethnicity all of which were included in adjusted models. Initial association analysis with unadjusted logistic regression revealed putative association with seven SNPs (p < 0.05). All SNPs were then assessed in an adjusted model. Of the six SNPs that remained significant, three of them were assigned to the CLEC4E gene (rs12302046, rs10841847, and rs11046143). Of these, only rs10841847 passed FDR adjustment for multiple testing. Adjusted regression analyses showed that the minor allele at rs10841847 associated with higher risk of developing PTB (OR = 1.55, CI = 1.22-1.96, p-value = 0.00036). Based on these initial association tests, CLEC4E seemed to be the predictor of interest for PTB risk in this population. Haplotype analysis (2-SNP and 3-SNP windows) showed that minor alleles in segments including rs10841847 were the only ones to pass the threshold of global significance, compared to other haplotypes (p-value < 0.05). Linkage disequilibrium patterns showed that rs12302046 is in high LD with rs10841847 (r2 = 0.67), and all other SNPs lost significance when adjusted for rs10841847 effects. These findings indicate that rs10841847 in CLEC4E is the single best predictor of pulmonary tuberculosis risk in our study population. These results provide evidence for the hypothesis that genetic variation of CLEC4E influences risk to TB in Guinea-Bissau.

Integrative genomics analysis identifies promising SNPs and genes implicated in tuberculosis risk based on multiple omics datasets

More than 10 GWASs have reported numerous genetic loci associated with tuberculosis (TB). However, the functional effects of genetic variants on TB remains largely unknown. In the present study, by combining a reported GWAS summary dataset (N = 452,264) with 3 independent eQTL datasets (N = 2,242) and other omics datasets downloaded from public databases, we conducted an integrative genomics analysis to highlight SNPs and genes implicated in TB risk. Based on independent biological and technical validations, we prioritized 26 candidate genes with eSNPs significantly associated with gene expression and TB susceptibility simultaneously; such as, CDC16 (rs7987202, rs9590408, and rs948182) and RCN3 (rs2946863, rs2878342, and rs3810194). Based on the network-based enrichment analysis, we found these 26 highlighted genes were jointly connected to exert effects on TB susceptibility. The co-expression patterns among these 26 genes were remarkably changed according to Mycobacterium tuberculosis (MTB) infection status. Based on 4 independent gene expression datasets, 21 of 26 genes (80.77%) showed significantly differential expressions between TB group and control group in mesenchymal stem cells, mice blood and lung tissues, as well as human alveolar macrophages. Together, we provide robust evidence to support 26 highlighted genes as important candidates for TB.

A regulatory variant in the C1Q gene cluster is associated with tuberculosis susceptibility and C1qA plasma levels in a South African population

Several genetic studies have implicated genes that encode for components of the innate immune response in tuberculosis (TB) susceptibility. The complement system is an early player in the innate immune response and provides the host with initial protection by promoting phagocytosis of apoptotic or necrotic cells. The C1q molecule is the first component of the classical pathway that leads to the activation of complement by binding to immune complexes and is encoded by the C1Q gene cluster. We investigated variants in this region to determine its association with TB susceptibility. Five single nucleotide polymorphisms (SNPs) (rs12033074, rs631090, rs172378, rs587585, and rs665691) were genotyped using TaqMan® SNP assays in 456 TB cases and 448 healthy controls and analysed by logistic regression models. The rs587585 variant showed a significant additive allelic association where the minor G allele was found more frequently in TB cases than in controls in both the discovery (p = 0.023; OR = 1.30; 95% CI, 1.04-1.64) and validation cohort (p = 0.038; OR = 1.31; 95% CI, 1.22-1.40). In addition, we detected increased C1qA expression when comparing cases and controls (p = 0.037) and linked this to a dosage effect of the G allele, which increased C1qA expression in TB cases. This is the first study to report the association of C1Q gene polymorphisms with progression to tuberculosis.

Human global and population-specific genetic susceptibility to Mycobacterium tuberculosis infection and disease

PURPOSE OF REVIEW

Multiple lines of evidence support a role of the host genetic component in Mycobacterium tuberculosis infection and disease progression. However, genomic studies of tuberculosis susceptibility have been disappointing compared with that of other complex disorders. Recently the field has explored alternative strategies to facilitate locus discovery. Results emanating from these efforts during the last 18 months are addressed in this review.

RECENT FINDINGS

There has been a renewed focus on the refinement of phenotypic definitions of infection and disease as well as on age-related, sex-specific and population-specific effects. Genome-wide association studies have yielded candidate genes but the findings have not always been transferable to all population groups. Candidate gene association studies remain popular as it is used for GWAS replication and is affordable, particularly in lower and middle-income countries. Pharmacogenetic studies involving tuberculosis drugs may locate variants that can be cost-effectively genotyped to identify individuals at risk of developing adverse events during treatment.

SUMMARY

Additional GWAS and candidate gene association studies of crudely defined study participants are unlikely to make further important contributions to the TB susceptibility field. Instead refined phenotyping will allow the elucidation of genetic mechanisms contributing to infection and disease in distinct populations and the calculation of polygenic risk scores.

Sex Disparity in Severity of Lung Lesions in Newly Identified Tuberculosis Is Age-Associated

Background: The age-associated characteristic of computed tomography (CT) images of tuberculosis (TB) and the reason for male bias in TB are still not clear. Methods: We compared the CT images, clinical inflammatory indices and sputum bacterial counts between 594 non-smoking men and women with newly diagnosed TB with matched large span of ages from 15 to 92 years old. Logistic regression analyses were used to identify the cavity-associated factors of men and women, separately and in combination. Results: Sputum bacterial counts, ratio of cavities, lung injury scores, and level of C reactive protein were significantly higher in men than in women with ages from 15 to 74, but not in cases older than 75. In CT images, thick walled cavity, cicatricial emphysema and parenchymal bands were present in men at ages of 15-74 more than matched women. Ratios of cases with lobular emphysema and pleural effusion were higher in men after age of 56. While ratios of cases with parenchymal bands, calcification, pleural effusion, pleural thickening, lobular emphysema and bronchovascular distortion increased with aging, those of centrilobular nodules, micronodules and tree in bud decreased with aging in men. Erythrocyte sedimentation rate (ESR) increased with aging, but no differences were found between men and women in ESR or T-SPOT TB tests. Higher complement C4 and lower body mass index in men and positive result in anti-TB antibody test in women were strongly associated with the presence of cavity. Conclusions: The sex bias in TB is age-associated. TB prevention, treatment and research should take differences of sex and age into account.

Homozygosity for TYK2 P1104A underlies tuberculosis in about 1% of patients in a cohort of European ancestry

The human genetic basis of tuberculosis (TB) has long remained elusive. We recently reported a high level of enrichment in homozygosity for the common TYK2 P1104A variant in a heterogeneous cohort of patients with TB from non-European countries in which TB is endemic. This variant is homozygous in ∼1/600 Europeans and ∼1/5,000 people from other countries outside East Asia and sub-Saharan Africa. We report a study of this variant in the UK Biobank cohort. The frequency of P1104A homozygotes was much higher in patients with TB (6/620, 1%) than in controls (228/114,473, 0.2%), with an odds ratio (OR) adjusted for ancestry of 5.0 [95% confidence interval (CI): 1.96-10.31, P = 2 × 10-3]. Conversely, we did not observe enrichment for P1104A heterozygosity, or for TYK2 I684S or V362F homozygosity or heterozygosity. Moreover, it is unlikely that more than 10% of controls were infected with Mycobacterium tuberculosis, as 97% were of European genetic ancestry, born between 1939 and 1970, and resided in the United Kingdom. Had all of them been infected, the OR for developing TB upon infection would be higher. These findings suggest that homozygosity for TYK2 P1104A may account for ∼1% of TB cases in Europeans.

Evaluating the Accuracy of Imputation Methods in a Five-Way Admixed Population

Genotype imputation is a powerful tool for increasing statistical power in an association analysis. Meta-analysis of multiple study datasets also requires a substantial overlap of SNPs for a successful association analysis, which can be achieved by imputation. Quality of imputed datasets is largely dependent on the software used, as well as the reference populations chosen. The accuracy of imputation of available reference populations has not been tested for the five-way admixed South African Colored (SAC) population. In this study, imputation results obtained using three freely-accessible methods were evaluated for accuracy and quality. We show that the African Genome Resource is the best reference panel for imputation of missing genotypes in samples from the SAC population, implemented via the freely accessible Sanger Imputation Server.