A panel of ancestry informative markers for the complex five-way admixed South African coloured population

Genetic admixture predictors of fetal alcohol spectrum disorders (FASD) in a South African population

Ancestrally admixed populations are underrepresented in genetic studies of complex diseases, which are still dominated by European-descent populations. This is relevant not only from a representation standpoint but also because of admixed populations' unique features, including being enriched for rare variants, for which effect sizes are disproportionately larger than common polymorphisms. Furthermore, results from these populations may be generalizable to other populations. The South African Cape Coloured (SACC) population is genetically admixed and has one of the highest prevalences of fetal alcohol spectrum disorders (FASD) worldwide. We profiled its admixture and examined associations between ancestry profiles and FASD outcomes using two longitudinal birth cohorts (N=308 mothers, 280 children) designed to examine effects of prenatal alcohol exposure on development. Participants were genotyped via MEGAex array to capture common and rare variants. Rare variants were overrepresented in our SACC cohorts, with numerous polymorphisms being monomorphic in other reference populations (e.g., ∼30,000 and ∼ 221,000 variants in gnomAD European and Asian populations, respectively). The cohorts showed global African (51 %; Bantu and San); European (26 %; Northern/Western); South Asian (18 %); and East Asian (5 %; largely Southern regions) ancestries. The cohorts exhibited high rates of homozygosity (6 %), with regions of homozygosity harboring more deleterious variants when lying within African local-ancestry genomic segments. Both maternal and child ancestry profiles were associated with higher FASD risk, and maternal and child ancestry-by-prenatal alcohol exposure interaction effects were seen on child cognition. Our findings indicate that the SACC population may be a valuable asset to identify novel disease-associated genetic loci for FASD and other diseases.

Blood donor return behavior in South Africa and the United States before and during the COVID-19 pandemic

BACKGROUND

Studies preceding the COVID-19 pandemic found that slower time-to-return was associated with first-time, deferred, and mobile drive blood donors. How donor return dynamics changed during the COVID-19 pandemic is not well understood.

METHODS

We analyzed visits by whole blood donors from 2017 to 2022 in South Africa (SA) and the United States (US) stratified by mobile and fixed environment, first-time and repeat donor status, and pre-COVID19 (before March 2020) and intra-COVID19. We used Kaplan-Meier curves to characterize time-to-return, cumulative incidence functions to analyze switching between donation environments, and Cox proportional hazards models to analyze factors influencing time-to-return.

RESULTS

Overall time-to-return was shorter in SA. Pre-COVID19, the proportion of donors returning within a year of becoming eligible was lower for deferred donors in both countries regardless of donation environment and deferral type. Intra-COVID19, the gap between deferred and non-deferred donors widened in the US but narrowed in SA, where efforts to schedule return visits from deferred donors were intensified, particularly for non-hemoglobin-related deferrals. Intra-COVID19, the proportion of donors returning within a year in SA was higher for deferred first-time donors (>81%) than for successful first-time donors (80% at fixed sites; 69% at mobile drives).

CONCLUSIONS

The pandemic complicated efforts to recruit new donors and schedule returning visits after completed donations. Concerted efforts to improve time-to-return for deferred donors helped mitigate donation loss in SA during the public health emergency.

Genetic admixture predictors of fetal alcohol spectrum disorders (FASD) in the South African Cape Coloured population

Ancestrally admixed populations are underrepresented in genetic studies of complex diseases, which are still dominated by European-descent populations. This is relevant not only from a representation standpoint but also because of admixed populations' unique features, including being enriched for rare variants, for which effect sizes are disproportionately larger than common polymorphisms. Furthermore, results from these populations may be generalizable to other populations. The South African Cape Coloured (SACC) population is genetically admixed, with one of the highest prevalences of fetal alcohol spectrum disorders (FASD) worldwide. We profiled its admixture and examined associations between ancestry profiles and FASD outcomes using two longitudinal birth cohorts ( N =308 mothers, 280 children) designed to examine effects of prenatal alcohol exposure on development. Participants were genotyped via MEGA-ex array to capture common and rare variants. Rare variants were overrepresented in our SACC cohorts, with numerous polymorphisms being monomorphic in other reference populations (e.g., ∼30,000 and ∼221,000 variants in gnomAD European and Asian populations, respectively). The cohorts showed global African (51%; Bantu and San); European (26%; Northern/Western); South Asian (18%); and East Asian (5%; largely Southern regions) ancestries. The cohorts exhibited high rates of homozygosity (6%), with regions of homozygosity harboring more deleterious variants when lying within African local-ancestry genomic segments. Both maternal and child ancestry profiles were associated with FASD risk and altered severity of prenatal alcohol exposure-related cognitive deficits in the child. Our findings indicate that the SACC population may be a valuable asset to identify novel disease-associated genetic loci for FASD and other diseases.

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.

Evolutionary Genetics and Admixture in African Populations

As the ancestral homeland of our species, Africa contains elevated levels of genetic diversity and substantial population structure. Importantly, African genomes are heterogeneous: They contain mixtures of multiple ancestries, each of which have experienced different evolutionary histories. In this review, we view population genetics through the lens of admixture, highlighting how multiple demographic events have shaped African genomes. Each of these historical vignettes paints a recurring picture of population divergence followed by secondary contact. First, we give a brief overview of genetic variation in Africa and examine deep population structure within Africa, including the evidence of ancient introgression from archaic "ghost" populations. Second, we describe the genetic legacies of admixture events that have occurred during the past 10,000 years. This includes gene flow between different click-speaking Khoe-San populations, the stepwise spread of pastoralism from eastern to southern Africa, multiple migrations of Bantu speakers across the continent, as well as admixture from the Middle East and Europe into the Sahel region and North Africa. Furthermore, the genomic signatures of more recent admixture can be found in the Cape Peninsula and throughout the African diaspora. Third, we highlight how natural selection has shaped patterns of genetic variation across the continent, noting that gene flow provides a potent source of adaptive variation and that selective pressures vary across Africa. Finally, we explore the biomedical implications of population structure in Africa on health and disease and call for more ethically conducted studies of genetic variation in Africa.

A custom hybridisation enrichment forensic intelligence panel to infer biogeographic ancestry, hair and eye colour, and Y chromosome lineage

Massively parallel sequencing can provide genetic data for hundreds to thousands of loci in a single assay for various types of forensic testing. However, available commercial kits require an initial PCR amplification of short-to-medium sized targets which limits their application for highly degraded DNA. Development and optimisation of large PCR multiplexes also prevents creation of custom panels that target different suites of markers for identity, biogeographic ancestry, phenotype, and lineage markers (Y-chromosome and mtDNA). Hybridisation enrichment, an alternative approach for target enrichment prior to sequencing, uses biotinylated probes to bind to target DNA and has proven successful on degraded and ancient DNA. We developed a customisable hybridisation capture method, that uses individually mixed baits to allow tailored and targeted enrichment to specific forensic questions of interest. To allow collection of forensic intelligence data, we assembled and tested a custom panel of hybridisation baits to infer biogeographic ancestry, hair and eye colour, and paternal lineage (and sex) on modern male and female samples with a range of self-declared ancestries and hair/eye colour combinations. The panel correctly estimated biogeographic ancestry in 9/12 samples (75%) but detected European admixture in three individuals from regions with admixed demographic history. Hair and eye colour were predicted correctly in 83% and 92% of samples respectively, where intermediate eye colour and blond hair were problematic to predict. Analysis of Y-chromosome SNPs correctly assigned sex and paternal haplogroups, the latter complementing and supporting biogeographic ancestry predictions. Overall, we demonstrate the utility of this hybridisation enrichment approach to forensic intelligence testing using a combined suite of biogeographic ancestry, phenotype, and Y-chromosome SNPs for comprehensive biological profiling.

A Pipeline for Phasing and Genotype Imputation on Mixed Human Data (Parents-Offspring Trios and Unrelated Subjects) by Reviewing Current Methods and Software

Genotype imputation has become an essential prerequisite when performing association analysis. It is a computational technique that allows us to infer genetic markers that have not been directly genotyped, thereby increasing statistical power in subsequent association studies, which consequently has a crucial impact on the identification of causal variants. Many features need to be considered when choosing the proper algorithm for imputation, including the target sample on which it is performed, i.e., related individuals, unrelated individuals, or both. Problems could arise when dealing with a target sample made up of mixed data, composed of both related and unrelated individuals, especially since the scientific literature on this topic is not sufficiently clear. To shed light on this issue, we examined existing algorithms and software for performing phasing and imputation on mixed human data from SNP arrays, specifically when related subjects belong to trios. By discussing the advantages and limitations of the current algorithms, we identified LD-based methods as being the most suitable for reconstruction of haplotypes in this specific context, and we proposed a feasible pipeline that can be used for imputing genotypes in both phased and unphased human data.

Unveiling forensically relevant biogeographic, phenotype and Y-chromosome SNP variation in Pakistani ethnic groups using a customized hybridisation enrichment forensic intelligence panel

Massively parallel sequencing following hybridisation enrichment provides new opportunities to obtain genetic data for various types of forensic testing and has proven successful on modern as well as degraded and ancient DNA. A customisable forensic intelligence panel that targeted 124 SNP markers (67 ancestry informative markers, 23 phenotype markers from the HIrisplex panel, and 35 Y-chromosome SNPs) was used to examine biogeographic ancestry, phenotype and sex and Y-lineage in samples from different ethnic populations of Pakistan including Pothwari, Gilgit, Baloach, Pathan, Kashmiri and Siraiki. Targeted sequencing and computational data analysis pipeline allowed filtering of variants across the targeted loci. Study samples showed an admixture between East Asian and European ancestry. Eye colour was predicted accurately based on the highest p-value giving overall prediction accuracy of 92.8%. Predictions were consistent with reported hair colour for all samples, using the combined highest p-value approach and step-wise model incorporating probability thresholds for light or dark shade. Y-SNPs were successfully recovered only from male samples which indicates the ability of this method to identify biological sex and allow inference of Y-haplogroup. Our results demonstrate practicality of using hybridisation enrichment and MPS to aid in human intelligence gathering and will open many insights into forensic research in South Asia.

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.

Perspectives on establishing a public cord blood inventory in South Africa

The South African population is highly diverse, both ethnically and genetically. This diversity is particularly true for the African ancestry and various mixed ancestry population groups. These groups are under-represented in national and international bone marrow and peripheral blood donor registries, making it challenging to identify HLA-matched and mismatched unrelated donors when patients from these groups require allogeneic hematopoietic stem and progenitor cell transplantation. In most high-income countries, banked cord blood (CB) units provide an attractive source of hematopoietic progenitor cells for genetically diverse populations. SA does not have a public CB inventory, leaving many patients without access to this important treatment modality. Haploidentical transplantation provides an alternative. In recent years, the use of post-transplant cyclophosphamide has significantly reduced the incidence of graft-versus-host disease after haploidentical transplantation and has improved transplantation outcomes. However, it is difficult to identify suitable haploidentical donors in SA because of family disruption and a high prevalence of HIV. Here the authors provide a brief historical overview of the ethnic and genetic diversity of the country and region. The authors provide a southern African perspective on HLA diversity, consider the allogeneic hematopoietic stem and progenitor cell transplantation landscape and explore the need to establish a public CB bank (CBB) in SA. The health policy and regulatory frameworks that will impact on a CBB in the country SA are also explored. Finally, the authors discuss several matters we believe require attention when considering the establishment of a sustainable public CBB in the South African context.

Prevalence of anti-SARS-CoV-2 antibodies among blood donors in Northern Cape, KwaZulu-Natal, Eastern Cape, and Free State provinces of South Africa in January 2021

Background

Population-level estimates of prevalence of anti-SARS-CoV-2 antibody positivity (seroprevalence) is a crucial epidemiological indicator for tracking the Covid-19 epidemic. Such data are in short supply, both internationally and in South Africa. The South African blood services (the South African National Blood Service, SANBS and the Western Cape Blood Service, WCBS) are coordinating a nationally representative survey of blood donors, which it is hoped can become a cost-effective surveillance method with validity for community-level seroprevalence estimation.

Methods

Leveraging existing arrangements, SANBS human research ethics committee permission was obtained to test blood donations collected on predefined days (7th, 10th, 12th, 15th, 20th, 23th and 25th January) for anti-SARS-CoV-2 antibodies, using the Roche Elecsys Anti-SARS-CoV-2 assay on the cobas e411 platform currently available in the blood services' donation testing laboratories. Using standard methods, prevalence analysis was done by province, age and race, allowing age to be regarded as either a continuous or categorical variable. Testing was performed in the Eastern Cape (EC), Free State (FS), KwaZulu Natal (ZN) and Northern Cape (NC) provinces.

Results

We report on data from 4858 donors - 1457 in EC; 463 in NC; 831 in FS and 2107 in ZN. Prevalence varied substantially across race groups and between provinces, with seroprevalence among Black donors consistently several times higher than among White donors, and the other main population groups (Coloured and Asian) not consistently represented in all provinces. There is no clear evidence that seroprevalence among donors varies by age. Weighted net estimates of prevalence (in the core age range 15-69) by province (compared with official clinically-confirmed COVID-19 case rates in mid-January 2021) are: EC-63%(2.8%), NC-32%(2.2%), FS-46%(2.4%), and ZN-52%(2.4%).

Conclusions

Our study demonstrates substantial differences in dissemination of SARS-CoV-2 infection between different race groups, most likely explained by historically based differences in socio-economic status and housing conditions. As has been seen in other areas, even such high seroprevalence does not guarantee population-level immunity against new outbreaks - probably due to viral evolution and waning of antibody neutralization. Despite its limitations, notably a 'healthy donor' effect, it seems plausible that these estimates are reasonably generalisable to actual population level anti-SARS-CoV-2 seroprevalence, but should be further verified.

Genetic variation in neuropeptide Y interacts with childhood trauma to influence anxiety sensitivity

BACKGROUND AND OBJECTIVES

Anxiety sensitivity (AS) refers to a fear of the negative implications of anxiety, and arises due to gene-environment interactions. We investigated whether genetic variation in two neuropeptides implicated in the stress response, neuropeptide Y (NPY) and pituitary adenylate cyclase-activating polypeptide receptor 1, interacted with childhood trauma (CT) to influence AS.

DESIGN AND METHODS

This cross-sectional study examined the CT x genetic variant effects on AS in 951 adolescents who self-identified as Xhosa or South African Colored (SAC) ethnicity.

RESULTS

In Xhosa females, the NPY rs5573 A allele and rs3037354 deletion variant were associated with increased (p = 0.035) and decreased (p = 0.034) AS, respectively. The interaction of CT and the NPY rs5574 A allele increased AS in SAC female participants (p = 0.043). The rs3037354 deletion variant protected against AS with increased CT in SAC male participants (p = 0.011).

CONCLUSIONS

The NPY rs5574 A allele and rs3037354 deletion variant interact with CT to act as risk and protective factors, respectively, for AS in an ethnicity- and sex- differentiated manner. Our results reaffirm the role of NPY and gene-environment interactions in anxiety-related behaviors and reinforce the need for psychiatric genetics studies in diverse populations.

Incidence of motor neuron disease/amyotrophic lateral sclerosis in South Africa: a 4-year prospective study

BACKGROUND AND PURPOSE

Little is known about the epidemiological features of amyotrophic lateral sclerosis (ALS) in sub-Saharan Africa, and data from the region are limited to clinical series or case reports. The aim of the study was to investigate the incidence rate and presentation of ALS in an ethnically diverse region of South Africa.

METHODS

We performed a 4-year prospective incidence study in the Western Cape Province of South Africa between 1 July 2014 and 30 June 2018, and used a two-source capture-recapture method for case ascertainment. Age- and sex-adjusted incidence rates (ASAIRs) were calculated using the 2010 US population as the reference.

RESULTS

A total of 203 incident cases were identified over the study period, resulting in a crude incidence rate (IR) of 1.09 [95% confidence interval (CI) 0.94-1.24] per 100 000 person-years in the at-risk population (aged >15 years). Capture-recapture analysis resulted in an estimated IR of 1.11 (95% CI 1.01-1.22) per 100 000 person-years. The ASAIR was 1.67 (95% CI 1.09-2.26) overall; 1.99 (95% CI 1.60-2.39) for men and 1.37 (95% CI 1.06-1.68) for women. When analysed separately, there was a substantial difference in ASAIRs between the different population groups, with the highest in the European ancestry group (2.62; 95% CI 2.49-2.75), the lowest in the African ancestry group (0.56, 95% CI 0.0-1.23), and an ASAIR in between these two in the mixed ancestry group (1.09, 95% CI 0.80-1.37).

CONCLUSION

The overall incidence of ALS in the Western Cape Province of South Africa appears to be lower than in North African and Western countries, but higher than in Asian countries. As suggested by previous epidemiological studies, ALS may be less frequent in people of African ancestry.

Prospective avenues for human population genomics and disease mapping in southern Africa

Population substructure within human populations is globally evident and a well-known confounding factor in many genetic studies. In contrast, admixture mapping exploits population stratification to detect genotype-phenotype correlations in admixed populations. Southern Africa has untapped potential for disease mapping of ancestry-specific disease risk alleles due to the distinct genetic diversity in its populations compared to other populations worldwide. This diversity contributes to a number of phenotypes, including ancestry-specific disease risk and response to pathogens. Although the 1000 Genomes Project significantly improved our understanding of genetic variation globally, southern African populations are still severely underrepresented in biomedical and human genetic studies due to insufficient large-scale publicly available data. In addition to a lack of genetic data in public repositories, existing software, algorithms and resources used for imputation and phasing of genotypic data (amongst others) are largely ineffective for populations with a complex genetic architecture such as that seen in southern Africa. This review article, therefore, aims to summarise the current limitations of conducting genetic studies on populations with a complex genetic architecture to identify potential areas for further research and development.

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.

Putting RFMix and ADMIXTURE to the test in a complex admixed population

Background

Global and local ancestry inference in admixed human populations can be performed using computational tools implementing distinct algorithms. The development and resulting accuracy of these tools has been tested largely on populations with relatively straightforward admixture histories but little is known about how well they perform in more complex admixture scenarios.

Results

Using simulations, we show that RFMix outperforms ADMIXTURE in determining global ancestry proportions even in a complex 5-way admixed population, in addition to assigning local ancestry with an accuracy of 89%. The ability of RFMix to determine global and local ancestry to a high degree of accuracy, particularly in admixed populations provides the opportunity for more accurate association analyses.

Conclusion

This study highlights the utility of the extension of computational tools to become more compatible to genetically structured populations, as well as the need to expand the sampling of diverse world-wide populations. This is particularly noteworthy as modern-day societies are becoming increasingly genetically complex and some genetic tools and commonly used ancestral populations are less appropriate. Based on these caveats and the results presented here, we suggest that RFMix be used for both global and local ancestry estimation in world-wide complex admixture scenarios particularly when including these estimates in association studies.

Genetic Affinities among Southern Africa Hunter-Gatherers and the Impact of Admixing Farmer and Herder Populations

Southern African indigenous groups, traditionally hunter-gatherers (San) and herders (Khoekhoe), are commonly referred to as “Khoe-San” populations and have a long history in southern Africa. Their ancestors were largely isolated up until ∼2,000 years ago before the arrival of pastoralists and farmers in southern Africa. Assessing relationships among regional Khoe-San groups has been challenging due to admixture with immigrant populations that obscure past population affinities and gene flow among these autochthonous communities. We re-evaluate a combined genome-wide data set of previously published southern Africa Khoe-San populations in conjunction with novel data from Khoe-San individuals collected in Xade (Central Kalahari Game Reserve, Botswana) prior to their resettlement outside the reserve. After excluding regions in the genome that trace their ancestry to recent migrant groups, the genetic diversity of 20 Khoe-San groups fitted an isolation-by-distance model. Even though isolation-by-distance explained most genetic affinities between the different autochthonous groups, additional signals of contact between Khoe-San groups could be detected. For instance, we found stronger genetic affinities, than what would be explained by isolation-by-distance gene flow, between the two geographically separated Khoe-San groups, who speak branches of the Kx’a-language family (ǂHoan and Ju). We also scanned the genome-wide data for signals of adaptive gene flow from farmers/herders into Khoe-San groups and identified a number of genomic regions potentially introduced by the arrival of the new groups. This study provides a comprehensive picture of affinities among Khoe-San groups, prior to the arrival of recent migrants, and found that these affinities are primarily determined by the geographic landscape.

An ancestry informative marker panel design for individual ancestry estimation of Hispanic population using whole exome sequencing data

Background

Europeans and American Indians were major genetic ancestry of Hispanics in the U.S. These ancestral groups have markedly different incidence rates and outcomes in many types of cancers. Therefore, the genetic admixture may cause biased genetic association study with cancer susceptibility variants specifically in Hispanics. For example, the incidence rate of liver cancer has been shown with substantial disparity between Hispanic, Asian and non-Hispanic white populations. Currently, ancestry informative marker (AIM) panels have been widely utilized with up to a few hundred ancestry-informative single nucleotide polymorphisms (SNPs) to infer ancestry admixture. Notably, current available AIMs are predominantly located in intron and intergenic regions, while the whole exome sequencing (WES) protocols commonly used in translational research and clinical practice do not cover these markers. Thus, it remains challenging to accurately determine a patient’s admixture proportion without additional DNA testing.

Results

In this study we designed an unique AIM panel that infers 3-way genetic admixture from three distinct and selective continental populations (African (AFR), European (EUR), and East Asian (EAS)) within evolutionarily conserved exonic regions. Initially, about 1 million exonic SNPs from selective three populations in the 1000 Genomes Project were trimmed by their linkage disequilibrium (LD), restricted to biallelic variants, and finally we optimized to an AIM panel with 250 SNP markers, or the UT-AIM250 panel, using their ancestral informativeness statistics. Comparing to published AIM panels, UT-AIM250 performed better accuracy when we tested with three ancestral populations (accuracy: 0.995 ± 0.012 for AFR, 0.997 ± 0.007 for EUR, and 0.994 ± 0.012 for EAS). We further demonstrated the performance of the UT-AIM250 panel to admixed American (AMR) samples of the 1000 Genomes Project and obtained similar results (AFR, 0.085 ± 0.098; EUR, 0.665 ± 0.182; and EAS, 0.250 ± 0.205) to previously published AIM panels (Phillips-AIM34: AFR, 0.096 ± 0.127, EUR, 0.575 ± 0.290, and EAS, 0.330 ± 0.315; Wei-AIM278: AFR, 0.070 ± 0.096, EUR, 0.537 ± 0.267, and EAS, 0.393 ± 0.300). Subsequently, we applied the UT-AIM250 panel to a clinical dataset of 26 self-reported Hispanic patients in South Texas with hepatocellular carcinoma (HCC). We estimated the admixture proportions using WES data of adjacent non-cancer liver tissues (AFR, 0.065 ± 0.043; EUR, 0.594 ± 0.150; and EAS, 0.341 ± 0.160). Similar admixture proportions were identified from corresponding tumor tissues. In addition, we estimated admixture proportions of The Cancer Genome Atlas (TCGA) collection of hepatocellular carcinoma (TCGA-LIHC) samples (376 patients) using the UT-AIM250 panel. The panel obtained consistent admixture proportions from tumor and matched normal tissues, identified 3 possible incorrectly reported race/ethnicity, and/or provided race/ethnicity determination if necessary.

Conclusions

Here we demonstrated the feasibility of using evolutionarily conserved exonic regions to infer admixture proportions and provided a robust and reliable control for sample collection or patient stratification for genetic analysis. R implementation of UT-AIM250 is available at https://github.com/chenlabgccri/UT-AIM250.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

GRAF-pop: A Fast Distance-Based Method To Infer Subject Ancestry from Multiple Genotype Datasets Without Principal Components Analysis

Inferring subject ancestry using genetic data is an important step in genetic association studies, required for dealing with population stratification. It has become more challenging to infer subject ancestry quickly and accurately since large amounts of genotype data, collected from millions of subjects by thousands of studies using different methods, are accessible to researchers from repositories such as the database of Genotypes and Phenotypes (dbGaP) at the National Center for Biotechnology Information (NCBI). Study-reported populations submitted to dbGaP are often not harmonized across studies or may be missing. Widely-used methods for ancestry prediction assume that most markers are genotyped in all subjects, but this assumption is unrealistic if one wants to combine studies that used different genotyping platforms. To provide ancestry inference and visualization across studies, we developed a new method, GRAF-pop, of ancestry prediction that is robust to missing genotypes and allows researchers to visualize predicted population structure in color and in three dimensions. When genotypes are dense, GRAF-pop is comparable in quality and running time to existing ancestry inference methods EIGENSTRAT, FastPCA, and FlashPCA2, all of which rely on principal components analysis (PCA). When genotypes are not dense, GRAF-pop gives much better ancestry predictions than the PCA-based methods. GRAF-pop employs basic geometric and probabilistic methods; the visualized ancestry predictions have a natural geometric interpretation, which is lacking in PCA-based methods. Since February 2018, GRAF-pop has been successfully incorporated into the dbGaP quality control process to identify inconsistencies between study-reported and computationally predicted populations and to provide harmonized population values in all new dbGaP submissions amenable to population prediction, based on marker genotypes. Plots, produced by GRAF-pop, of summary population predictions are available on dbGaP study pages, and the software, is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/Software.cgi.

The prevalence of human T-lymphotropic virus type 1 & 2 (HTLV-1/2) in South African blood donors

BACKGROUND AND OBJECTIVES

Donated blood is not currently screened for human T-cell lymphotropic virus (HTLV) in South Africa. Several small studies have detected HTLV-1 in South Africa, but prevalence by geographic region or population group is unavailable.

MATERIALS AND METHODS

We performed a large seroprevalence study of South African blood donors during 3 months in 2013. All geographic regions except the Western Cape were included, and Black and Coloured (local term for mixed race) donors were oversampled. Identity-unlinked plasma samples were screened with the Abbott Prism HTLV-1/2 assay, and repeatedly reactive samples were tested by the Inno-LIA HTLV-1/2 Score confirmatory assay. Odds ratios were calculated with multivariable logistic regression.

RESULTS

Of 46 752 donors tested, 133 (0·28%) were initially reactive, 111 (0·24%) repeatedly reactive and 57 (0·12%) confirmed positive for HTLV-1; none were HTLV-2 positive. Prevalence was 0·062% weighted to annual blood donations but highly concentrated in the Black population group (OR = 20·24 CI: 2·77-147·88); higher in females than males (OR = 1·81 CI: 1·06-3·08); and in donors aged >50 years compared to ages 16-19 (OR = 6·4 CI: 2·95-13·86). After controlling for age, sex and population group, there was no difference in prevalence between new and repeat blood donors or among geographic regions within South Africa.

CONCLUSIONS

We conclude that HTLV-1 infection is widespread among the Black population of South Africa and its epidemiology is similar to other endemic areas. Because South Africa is increasing its recruitment of Black blood donors, the implications for blood screening require further consideration.

The accuracy of the Transition Analysis of aging on a heterogenic South African population

Transition Analysis (TA) is an age estimation method for adult skeletons at death, using three skeletal traits together, namely, the cranial sutures (CS), the pubic symphyses (PS) and the auricular surfaces (AS). In forensic anthropology, estimating age at death is vital for positive identification of an individual; however, differences in skeletal elements have been observed between sex and population affinity. The aim of this study was to determine how accurately the TA could be applied to a population that differed from the original reference samples. Scoring methods, described by the authors of this analysis, were used to age individuals (n = 290) in the South African (SA) Mixed Ancestry population, which is a heterogenic population, consisting of many gene pools from different population groups. The results were compared to the actual ages of the skeletal remains, obtained through cadaver records, and the most likely age calculated by the TA, as well as age estimation by means of the three skeletal elements individually. Results showed that there was a significant correlation between the actual and estimated ages of the study population (Pearson correlation coefficient: 0.58). The AS was the most accurate trait, whereas the PS was the least accurate trait. In addition, females' estimated ages were more accurate than males. These findings correlated with other authors in terms of the significance of the TA when applied to different populations; however, the individual trait correlations were in disagreement with other authors. The results of this study are an important contribution to anthropological literature as they indicate that there are differences between populations and sex when determining age, but that the TA is still applicable for analysis on this heterogenic population.

Selection of Optimal Ancestry Informative Markers for Classification and Ancestry Proportion Estimation in Pigs

Using small sets of ancestry informative markers (AIMs) constitutes a cost-effective method to accurately estimate the ancestry proportions of individuals. This study aimed to generate a small and effective number of AIMs from ∼60 K single nucleotide polymorphism (SNP) data of porcine and estimate three ancestry proportions [East China pig (ECHP), South China pig (SCHP), and European commercial pig (EUCP)] from Asian breeds and European domestic breeds. A total of 186 samples of 10 pure breeds were divided into three groups: ECHP, SCHP, and EUCP. Using these samples and a one-vs.-rest SVM classifier, we found that using only seven AIMs could completely separate the three groups. Subsequently, we utilized supervised ADMIXTURE to calculate ancestry proportions and found that the 129 AIMs performed well on ancestry estimates when pseudo admixed individuals were used. Furthermore, another 969 samples of 61 populations were applied to evaluate the performance of the 129 AIMs. We also observed that the 129 AIMs were highly correlated with estimates using ∼60 K SNP data for three ancestry components: ECHP (Pearson correlation coefficient (r) = 0.94), SCHP (r = 0.94), and EUCP (r = 0.99). Our results provided an example of using a small number of pig AIMs for classifications and estimating ancestry proportions with high accuracy and in a cost-effective manner.

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.

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

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is a complex disease with a known human genetic component. Males seem to be more affected than females and in most countries the TB notification rate is twice as high in males than in females. While socio-economic status, behavior and sex hormones influence the male bias they do not fully account for it. Males have only one copy of the X chromosome, while diploid females are subject to X chromosome inactivation. In addition, the X chromosome codes for many immune-related genes, supporting the hypothesis that X-linked genes could contribute to TB susceptibility in a sex-biased manner. We report the first TB susceptibility genome-wide association study (GWAS) with a specific focus on sex-stratified autosomal analysis and the X chromosome. A total of 810 individuals (410 cases and 405 controls) from an admixed South African population were genotyped using the Illumina Multi Ethnic Genotyping Array, specifically designed as a suitable platform for diverse and admixed populations. Association testing was done on the autosome (8,27,386 variants) and X chromosome (20,939 variants) in a sex stratified and combined manner. SNP association testing was not statistically significant using a stringent cut-off for significance but revealed likely candidate genes that warrant further investigation. A genome wide interaction analysis detected 16 significant interactions. Finally, the results highlight the importance of sex-stratified analysis as strong sex-specific effects were identified on both the autosome and X chromosome.

Update on Inherited Retinal Disease in South Africa: Encouraging Diversity in Molecular Genetics

There is a glaring disparity in the populations included in genetic research; the majority of work involves European-derived cohorts, while other global populations - including Africans - are underrepresented. This is also true for the study of inherited retinal diseases. Being the most ancient of extant populations, African samples carry more variation than others, making them valuable for novel gene and variant discovery. The inclusion of diverse populations in research is essential to gain a more comprehensive understanding of genetic variation and molecular mechanisms of disease.

Toll-like receptor chaperone HSP90B1 and the immune response to Mycobacteria

Rationale

HSP90B1, also known as gp96, is a chaperone for multiple Toll-like receptors (TLRs) and is necessary for TLR-mediated inflammatory responses in murine myeloid cells. The molecule is also expressed in T-cells though its specific role is unknown. We hypothesized that human HSP90B1 regulates monocyte and T-cell responses to Mycobacterium tuberculosis (Mtb) and bacilli Calmette-Guerin (BCG) and that its variants are associated with susceptibility to TB disease.

Methods

We screened 17 haplotype-tagging SNPs in the HSP90B1 gene region for association with BCG-induced T-cell cytokine responses using both an ex-vivo whole blood assay (N = 295) and an intracellular cytokine staining assay (N = 180) on samples collected 10 weeks after birth. Using a case-control study design, we evaluated the same SNPs for association with TB disease in a South African pediatric cohort (N = 217 cases, 604 controls). A subset of these SNPs was evaluated for association with HSP90B1 expression in human monocytes, monocyte-derived dendritic cells, and T-cells using RT-PCR. Lastly, we used CRISPR/Cas9 gene editing to knock down HSP90B1 expression in a human monocyte cell line (U937). Knockdown and control cell lines were tested for TLR surface expression and control of Mtb replication.

Results

We identified three SNPs, rs10507172, rs10507173 and rs1920413, that were associated with BCG-induced IL-2 secretion (p = 0.017 for rs10507172 and p = 0.03 for rs10507173 and rs1920413, Mann-Whitney, dominant model). SNPs rs10507172 and rs10507173 were associated with TB disease in an unadjusted analysis (p = 0.036 and 0.025, respectively, dominant model) that strengthened with sensitivity analysis of the definite TB cases, which included only those patients with microbiologically confirmed Mtb (p = 0.007 and 0.012, respectively). Knockdowns of HSP90B1 in monocyte cell lines with CRISPR did not alter TLR2 surface expression nor influence Mtb replication relative to controls.

Conclusion

Among infants, an HSP90B1 gene-region variant is associated with BCG-induced IL-2 production and may be associated with protection from TB disease. HSP90B1 knockdown in human monocyte-like cell lines did not influence TLR2 surface localization nor Mtb replication. Together, these data suggest that HSP90B1 regulates T-cell, but not monocyte, responses to mycobacteria in humans.

The evolutionary history of Southern Africa

The genomic variability of Southern African groups is characterized by an exceptional degree of diversity, which is the result of long-term local evolutionary history, migrations and gene-flow. Over the last few years several investigations have identified and described signatures related to these processes, revealing how ancient and more recent events have shaped the structure and ancestry composition of local populations. Here we discuss recent insights into the genetic history of the Southernmost part of the African continent provided by the analysis of modern and ancient genomes. Future work is expected to clarify the population dynamics associated with the emergence of Homo sapiens across Africa and the details of the process of dispersion and admixture associated with the arrival of Bantu-speaking groups in the region.

The iron status of South African blood donors: balancing donor safety and blood demand

BACKGROUND

Several studies in developed countries have demonstrated high levels of iron deficiency (ID) among blood donors. There is a paucity of data for developing countries where blood shortages remain a major concern.

STUDY DESIGN AND METHODS

A total of 4412 donors were enrolled in the study. Specimens were collected for full blood count, iron, transferrin saturation, and ferritin assessment. Donor demographics were recorded. ID was indicated by a ferritin level of less than 20 ng/mL for men and less than 12 ng/mL for women. Anemia was defined as hemoglobin levels less than 12.5 g/dL. Regression models for predictors of ID were developed.

RESULTS

A total of 17.5% of all donors had ID, with 16.3% prevalence in women and 18.6% in men. Low hemoglobin had the highest association with ID (adjusted odds ratio [AOR], 11.078; 95% confidence interval [CI], 7.915-15.505); male donors had twice the odds of ID compared to female donors (AOR, 2.501; 95% CI, 1.964-3.185), while increasing age was associated with lower odds (AOD, 0.965; 95% CI, 0.956-0.975). Among male donors, an interdonation interval of less than 3 months (AOR, 2.679; 95% CI, 1.929-3.720) was associated with ID. Compared to other females combined, colored female donors (AOR, 2.335; 95% CI, 1.310-4.160) had higher odds and black female donors (AOR, 0.559; 95% CI, 0.369-0.845) lower odds of ID.

CONCLUSION

ID is common among South African donors; low hemoglobin, gender, ethnicity, and past donation history is independently associated with ID. Recommendations aimed at protecting donor health may increase blood shortages in South Africa.

Off-target phenotypes in forensic DNA phenotyping and biogeographic ancestry inference: A resource

With recent advances in DNA sequencing technologies it has become feasible and cost effective to genotype larger marker sets for forensic purposes. Two technologies that make use of the larger marker sets have come into focus in forensic research and applications; inference of biogeographic ancestry (BGA) and forensic DNA phenotyping (FDP). These methods hold the promise to reveal information about a yet unknown perpetrator from a DNA sample. In contrast, DNA-profiling, that is a standard practice in case work, relies on matching DNA-profiles between crime scene material and suspects on a database of DNA-profiles. Markers for DNA-profiling were developed under the premise to reveal as little additional information about the human source of the profile as possible, the rationale being that personal privacy rights have to be balanced against the public interest in solving a crime. The same argument holds for markers used in BGA and FDP; these markers might also reveal information on off-target phenotypes (OTPs), that go beyond BGA and the phenotypes targeted in FDP. In particular, health related OTPs might shift the balance between privacy protection and public interest. However, to our knowledge, there is currently no convenient resource available to incorporate knowledge on OTPs in BGA and FDP assay design and application. In order to provide such a resource, we performed a systematic search for OTPs associated with a comprehensive set of markers (1766 SNPs) used or suggested to be used for BGA inference and FDP. In this set, we identified a relatively small number of 27 SNPs (1.53%) that convey information on diverse health related OTPs such as cancer risk, induced asthma, or risk of alcoholism. Some of these SNPs are commonly used for FDP and BGA across different marker sets. We conclude that the effects of SNP markers used in FDP and BGA on OTPs are currently limited, with few exceptions that should be considered in a balanced decision on assay design and application.

Primary immunodeficiency diseases in a tuberculosis endemic region: challenges and opportunities

While individual primary immunodeficiency diseases (PIDs) are rare, collectively they represent a significant burden of disease. Recent estimates show that about one million people in Africa suffer from a PID. However, data from African PID registries reflect only a small percentage of the estimated prevalence. This disparity is partly due to the lack of PID awareness and the masking of PIDs by the endemic pathogens. Over three million tuberculosis (TB) cases were reported in Africa in 2016, with many of these from southern Africa. Despite concerted efforts to address this high burden of disease, the underlying genetic correlates of susceptibility to TB remain poorly understood. High penetrance mutations in immune system genes can cause PIDs that selectively predispose individuals to TB and other mycobacterial diseases. Additionally, the identification of individuals at a heightened risk of developing TB or of presenting with severe or disseminated TB due to their genetic ancestry is crucial to promote a positive treatment outcome. The screening for and identification of PID mutations in TB-endemic regions by next-generation sequencing (NGS) represents a promising approach to improve the understanding of what constitutes an effective immune response to TB, as well as the range of associated PIDs and phenotypes.

An Ancestry Informative Marker Set Which Recapitulates the Known Fine Structure of Populations in South Asia

The inference of genomic ancestry using ancestry informative markers (AIMs) can be useful for a range of studies in evolutionary genetics, biomedical research, and forensic analyses. However, the determination of AIMs for highly admixed populations with complex ancestries has remained a formidable challenge. Given the immense genetic heterogeneity and unique population structure of the Indian subcontinent, here we sought to derive AIMs that would yield a cohesive and faithful understanding of South Asian genetic origins. To discern the most optimal strategy for extracting AIMs for South Asians we compared three commonly used AIMs-determining methods namely, Infocalc, FST, and Smart Principal Component Analysis with ADMIXTURE, using previously published whole genome data from the Indian subcontinent. Our findings suggest that the Infocalc approach is likely most suitable for delineation of South Asian AIMs. In particular, Infocalc-2,000 (N = 2,000) appeared as the most informative South Asian AIMs panel that recapitulated the finer structure within South Asian genomes with high degree of sensitivity and precision, whereas a negative control with an equivalent number of randomly selected markers when used to interrogate the South Asian populations, failed to do so. We discuss the utility of all approaches under evaluation for AIMs derivation and interpreting South Asian genomic ancestries. Notably, this is the first report of an AIMs panel for South Asian ancestry inference. Overall these findings may aid in developing cost-effective resources for large-scale demographic analyses and foster expansion of our knowledge of human origins and disease, in the South Asian context.

Variation within voltage-gated calcium channel genes and antipsychotic treatment response in a South African first episode schizophrenia cohort

Voltage-gated calcium channels have been implicated in schizophrenia aetiology; however, little is known about their involvement in antipsychotic treatment response. This study investigated variants within the calcium channel subunit genes for association with antipsychotic treatment response in a first episode schizophrenia cohort. Twelve regulatory variants within seven genes were shown to be significantly associated with treatment outcome. Most notably, the CACNA1B rs2229949 CC genotype was associated with improved negative symptomology, where the C allele was predicted to abolish a miRNA-binding site (has-mir-5002-3p), suggesting a possible mechanism of action through which this variant may have an effect. These results implicate the calcium channel subunits in antipsychotic treatment response and suggest that increased activation of these channels may be explored to enhance or predict antipsychotic treatment outcome.

Hypothalamic-pituitary-adrenal axis variants and childhood trauma influence anxiety sensitivity in South African adolescents

Anxiety sensitivity (AS) is characterised by the fear of anxiety-related symptoms and is a risk factor for the development of anxiety-related disorders. We examined whether genetic variation in three stress response genes, CRHR1, NR3C1, and FKBP5, interact with childhood trauma (CT) to predict AS in South African adolescents. Xhosa (n = 634) and Coloured (n = 317) students completed self-report measures of AS and CT, and a total of eighteen polymorphisms within CRHR1, NR3C1, and FKBP5 were genotyped. Differences in AS based on genetic variation and CT were analysed within population and gender groups using multiple linear regression. Associations were found between AS and FKBP5 rs9296158 (p = 0.025) and rs737054 (p = 0.045) in Coloured males. Analysis of gene x CT interactions indicated that NR3C1 rs190488 CC-genotype, NR3C1 rs10482605 G-allele addition, and FKBP5 rs3800373 C-allele addition protect against AS with increasing CT in Xhosa females (p = 0.009), Xhosa males (p = 0.036) and Coloured males (p = 0.049), respectively. We identified two different protective single nucleotide polymorphism (SNP) combinations in a four-SNP CRHR1 haplotype in Coloured males. An analysis of the interaction between CT and a six-SNP FKBP5 haplotype in Coloured males revealed both protective and risk allelic combinations. Our results provide evidence for the influence of both genetic variation in CRHR1, NR3C1 and FKBP5, as well as CT x SNP interactions, on AS in South African adolescents. This study reinforces the importance of examining the influence of gene-environment (G X E) interactions within gender and population groups.

Sex-specific patterns in cortical and trabecular bone microstructure in the Kirsten Skeletal Collection, South Africa

OBJECTIVES

The purpose of this study was to provide bone histomorphometric reference data for South Africans of the Western Cape who likely dealt with health issues under the apartheid regime.

METHODS

The 206 adult individuals (_n female = 75, _n male = 131, mean = 47.9 ± 15.8 years) from the Kirsten Skeletal Collection, U. Stellenbosch, lived in the Cape Town metropole from the late 1960s to the mid-1990s. To study age-related changes in cortical and trabecular bone microstructure, photomontages of mid-thoracic rib cross-sections were quantitatively examined. Variables include relative cortical area (Rt.Ct.Ar), osteon population density (OPD), osteon area (On.Ar), bone volume fraction (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th), and trabecular spacing (Tb.Sp).

RESULTS

All cortical variables demonstrated significant relationships with age in both sexes, with women showing stronger overall age associations. Peak bone mass was compromised in some men, possibly reflecting poor nutritional quality and/or substance abuse issues throughout adolescence and early adulthood. In women, greater predicted decrements in On.Ar and Rt.Ct.Ar suggest a structural disadvantage with age, consistent with postmenopausal bone loss. Age-related patterns in trabecular bone microarchitecture are variable and difficult to explain. Except for Tb.Th, there are no statistically significant relationships with age in women. Men demonstrate significant negative correlations between BV/TV, Tb.N, and age, and a significant positive correlation between Tb.Sp and age.

CONCLUSIONS

This research highlights sex-specific differences in patterns of age-related bone loss, and provides context for discussion of contemporary South African bone health. While the study sample demonstrates indicators of poor bone quality, osteoporosis research continues to be under-prioritized in South Africa.

Whole-genome sequencing for an enhanced understanding of genetic variation among South Africans

The Southern African Human Genome Programme is a national initiative that aspires to unlock the unique genetic character of southern African populations for a better understanding of human genetic diversity. In this pilot study the Southern African Human Genome Programme characterizes the genomes of 24 individuals (8 Coloured and 16 black southeastern Bantu-speakers) using deep whole-genome sequencing. A total of ~16 million unique variants are identified. Despite the shallow time depth since divergence between the two main southeastern Bantu-speaking groups (Nguni and Sotho-Tswana), principal component analysis and structure analysis reveal significant (p < 10⁻⁶) differentiation, and F_ST analysis identifies regions with high divergence. The Coloured individuals show evidence of varying proportions of admixture with Khoesan, Bantu-speakers, Europeans, and populations from the Indian sub-continent. Whole-genome sequencing data reveal extensive genomic diversity, increasing our understanding of the complex and region-specific history of African populations and highlighting its potential impact on biomedical research and genetic susceptibility to disease.

African populations show a high level of genetic diversity and extensive regional admixture. Here, the authors sequence the whole genomes of 24 South African individuals of different ethnolinguistic origin and find substantive genomic divergence between two southeastern Bantu-speaking groups.

A Functional Toll-Interacting Protein Variant Is Associated with Bacillus Calmette-Guérin-Specific Immune Responses and Tuberculosis

RATIONALE

The molecular mechanisms that regulate tuberculosis susceptibility and bacillus Calmette-Guérin (BCG)-induced immunity are mostly unknown. However, induction of the adaptive immune response is a critical step in host control of Mycobacterium tuberculosis. Toll-interacting protein (TOLLIP) is a ubiquitin-binding protein that regulates innate immune responses, including Toll-like receptor signaling, which initiate adaptive immunity. TOLLIP variation is associated with susceptibility to tuberculosis, but the mechanism by which it regulates tuberculosis immunity is poorly understood.

OBJECTIVES

To identify functional TOLLIP variants and evaluate the role of TOLLIP variation on innate and adaptive immune responses to mycobacteria and susceptibility to tuberculosis.

METHODS

We used human cellular immunology approaches to characterize the role of a functional TOLLIP variant on monocyte mRNA expression and M. tuberculosis-induced monocyte immune functions. We also examined the association of TOLLIP variation with BCG-induced T-cell responses and susceptibility to latent tuberculosis infection.

MEASUREMENTS AND MAIN RESULTS

We identified a functional TOLLIP promoter region single-nucleotide polymorphism, rs5743854, which was associated with decreased TOLLIP mRNA expression in infant monocytes. After M. tuberculosis infection, TOLLIP-deficient monocytes demonstrated increased IL-6, increased nitrite, and decreased bacterial replication. The TOLLIP-deficiency G/G genotype was associated with decreased BCG-specific IL-2⁺ CD4⁺ T-cell frequency and proliferation. This genotype was also associated with increased susceptibility to latent tuberculosis infection.

CONCLUSIONS

TOLLIP deficiency is associated with decreased BCG-specific T-cell responses and increased susceptibility to tuberculosis. We hypothesize that the heightened antibacterial monocyte responses after vaccination of TOLLIP-deficient infants are responsible for decreased BCG-specific T-cell responses. Activating TOLLIP may provide a novel adjuvant strategy for BCG vaccination.

The role of human host genetics in tuberculosis resistance

INTRODUCTION

Tuberculosis (TB) remains a public health problem: the latest estimate of new incident cases per year is a staggering 10.4 million. Despite this overwhelming number, the majority of the immunocompetent population can control infection with Mycobacterium tuberculosis. The human genome underlies the immune response and contributes to the outcome of TB infection. Areas covered: Investigations of TB resistance in the general population have closely mirrored those of other infectious diseases and initially involved epidemiological observations. Linkage and association studies, including studies of VDR, SLC11A1 and HLA-DRB1 followed. Genome-wide association studies of common variants, not necessarily sufficient for disease, became possible after technological advancements. Other approaches involved the identification of those individuals with rare disease-causing mutations that strongly predispose to TB, epistasis and the role of ethnicity in disease. Despite these efforts, infection outcome, on an individual basis, cannot yet be predicted. Expert commentary: The early identification of future disease progressors is necessary to stem the TB epidemic. Human genetics may contribute to this endeavour and could in future suggest pathways to target for disease prevention. This will however require concerted efforts to establish large, well-phenotyped cohorts from different ethnicities, improved genomic resources and a better understanding of the human genome architecture.

The SIGLEC14 null allele is associated with Mycobacterium tuberculosis- and BCG-induced clinical and immunologic outcomes

Humans exposed to Mycobacterium tuberculosis (Mtb) have variable susceptibility to tuberculosis (TB) and its outcomes. Siglec-5 and Siglec-14 are members of the sialic-acid binding lectin family that regulate immune responses to pathogens through inhibitory (Siglec-5) and activating (Siglec-14) domains. The SIGLEC14 coding sequence is deleted in a high proportion of individuals, placing a SIGLEC5-like gene under the expression of the SIGLEC14 promoter (the SIGLEC14 null allele) and causing expression of a Siglec-5 like protein in monocytes and macrophages. We hypothesized that the SIGLEC14 null allele was associated with Mtb replication in monocytes, T-cell responses to the BCG vaccine, and clinical susceptibility to TB. The SIGLEC14 null allele was associated with protection from TB meningitis in Vietnamese adults but not with pediatric TB in South Africa. The null allele was associated with increased IL-2 and IL-17 production following ex-vivo BCG stimulation of blood from 10 week-old South African infants vaccinated with BCG at birth. Mtb replication was increased in THP-1 cells overexpressing either Siglec-5 or Siglec-14 relative to controls. To our knowledge, this is the first study to demonstrate an association between SIGLEC expression and clinical TB, Mtb replication, or BCG-specific T-cell cytokines.

Human leukocyte antigen class I (A, B and C) allele and haplotype variation in a South African Mixed ancestry population

South Africa has a large (∼53million), ethnically diverse population (black African, Caucasian, Indian/Asian and Mixed ancestry) and a high disease burden (particularly HIV-1 and Mycobacterium tuberculosis). The Mixed ancestry population constitutes ∼9% of the total population and was established ∼365years ago in the Western Cape region through interracial mixing of black Africans, Europeans and Asians. Admixed populations present unique opportunities to identify genetic factors involved in disease susceptibility. Since HLA genes are important mediators of host immunity, we investigated HLA-A, -B and -C allele and haplotype diversity in 50 healthy, unrelated individuals recruited from the Mixed ancestry population.

A post-GWAS analysis of predicted regulatory variants and tuberculosis susceptibility

Utilizing data from published tuberculosis (TB) genome-wide association studies (GWAS), we use a bioinformatics pipeline to detect all polymorphisms in linkage disequilibrium (LD) with variants previously implicated in TB disease susceptibility. The probability that these variants had a predicted regulatory function was estimated using RegulomeDB and Ensembl's Variant Effect Predictor. Subsequent genotyping of these 133 predicted regulatory polymorphisms was performed in 400 admixed South African TB cases and 366 healthy controls in a population-based case-control association study to fine-map the causal variant. We detected associations between tuberculosis susceptibility and six intronic polymorphisms located in MARCO, IFNGR2, ASHAS2, ACACA, NISCH and TLR10. Our post-GWAS approach demonstrates the feasibility of combining multiple TB GWAS datasets with linkage information to identify regulatory variants associated with this infectious disease.

Population structure and infectious disease risk in southern Africa

The KhoeSan populations are the earliest known indigenous inhabitants of southern Africa. The relatively recent expansion of Bantu-speaking agropastoralists, as well as European colonial settlement along the south-west coast, dramatically changed patterns of genetic diversity in a region which had been largely isolated for thousands of years. Owing to this unique history, population structure in southern Africa reflects both the underlying KhoeSan genetic diversity as well as differential recent admixture. This population structure has a wide range of biomedical and sociocultural implications; such as changes in disease risk profiles. Here, we consolidate information from various population genetic studies that characterize admixture patterns in southern Africa with an aim to better understand differences in adverse disease phenotypes observed among groups. Our review confirms that ancestry has a direct impact on an individual's immune response to infectious diseases. In addition, we emphasize the importance of collaborative research, especially for populations in southern Africa that have a high incidence of potentially fatal infectious diseases such as HIV and tuberculosis.

Fine-mapping of antipsychotic response genome-wide association studies reveals novel regulatory mechanisms

Aim: Noncoding variation has demonstrated regulatory effects on disease treatment outcomes. This study investigated the potential functionality of previously implicated noncoding variants on schizophrenia treatment response. Materials & methods: Predicted regulatory potential of variation identified from antipsychotic response genome-wide association studies was determined. Prioritized variants were assessed for association(s) with treatment outcomes in a South African first episode schizophrenia cohort (n = 103). Results: Bioinformatic and association results implicated a relationship between regulatory variants, expression of MANBA, COL9A2 and NFKB1, and treatment response. Three SNPs were associated with poor outcomes (rs230493: p = 1.88 × 10-6; rs3774959: p = 1.75 × 10-5; and rs230504: p = 1.48 × 10-4). Conclusion: This study has thoroughly investigated previous GWAS to pinpoint variants that may play a causal role in poor schizophrenia treatment outcomes, and provides potential candidate genes for further study in the field of antipsychotic response.

Genetic ancestry in relation to the metabolic response to a US versus traditional Mexican diet: a randomized crossover feeding trial among women of Mexican descent

BACKGROUND/OBJECTIVES

Certain populations with a large proportion of indigenous American (IA) genetic ancestry may be evolutionarily adapted to traditional diets high in legumes and complex carbohydrates, and may have a detrimental metabolic response to US diets high in refined carbohydrates and added sugars. We tested whether IA ancestry modified the metabolic response to a US versus traditional Mexican diet in a controlled dietary intervention.

SUBJECTS/METHODS

First and second generation Mexican immigrant women (n=53) completed a randomized crossover feeding trial testing the effects of a US versus traditional Mexican diet. The metabolic response to the diets was measured by fasting serum concentrations of glucose, insulin, insulin-like growth factor-1 (IGF-1), IGF-binding protein-3 (IGFBP-3), adiponectin, C-reactive protein, interleukin-6 and computed homeostasis model assessment for insulin resistance (HOMA_IR). Blood collected at baseline was used for genotyping, and estimation of African, European and IA ancestries with the use of 214 ancestry informative markers.

RESULTS

The genetic ancestral background was 56% IA, 38% European and 6% African. Women in the highest IA ancestry tertile (>62%) were shorter in height, less educated and less acculturated to the US lifestyle, and tended to have higher waist-to-hip ratio compared with women in the middle and lowest IA ancestry tertiles, respectively. Compared with the US diet, the traditional Mexican diet tended to reduce glucose, insulin, IGF-1, IGFBP-3 and HOMA_IR among women in the middle IA ancestry group (IA ancestry ⩽45-62%), whereas having no effect on biomarkers related to inflammation.

CONCLUSIONS

We observed modest interactions between IA ancestry and the metabolic response to a US versus traditional Mexican diet among Mexican immigrant women.