An Unexpectedly Complex Architecture for Skin Pigmentation in Africans.

Approximately 15 genes have been directly associated with skin pigmentation variation in humans, leading to its characterization as a relatively simple trait. However, by assembling a global survey of quantitative skin pigmentation phenotypes, we demonstrate that pigmentation is more complex than previously assumed, with genetic architecture varying by latitude. We investigate polygenicity in the KhoeSan populations indigenous to southern Africa who have considerably lighter skin than equatorial Africans. We demonstrate that skin pigmentation is highly heritable, but known pigmentation loci explain only a small fraction of the variance. Rather, baseline skin pigmentation is a complex, polygenic trait in the KhoeSan. Despite this, we identify canonical and non-canonical skin pigmentation loci, including near SLC24A5, TYRP1, SMARCA2/VLDLR, and SNX13, using a genome-wide association approach complemented by targeted resequencing. By considering diverse, under-studied African populations, we show how the architecture of skin pigmentation can vary across humans subject to different local evolutionary pressures.

A weakly structured stem for human origins in Africa.

Despite broad agreement that Homo sapiens originated in Africa, considerable uncertainty surrounds specific models of divergence and migration across the continent1. Progress is hampered by a shortage of fossil and genomic data, as well as variability in previous estimates of divergence times1. Here we seek to discriminate among such models by considering linkage disequilibrium and diversity-based statistics, optimized for rapid, complex demographic inference2. We infer detailed demographic models for populations across Africa, including eastern and western representatives, and newly sequenced whole genomes from 44 Nama (Khoe-San) individuals from southern Africa. We infer a reticulated African population history in which present-day population structure dates back to Marine Isotope Stage 5. The earliest population divergence among contemporary populations occurred 120,000 to 135,000 years ago and was preceded by links between two or more weakly differentiated ancestral Homo populations connected by gene flow over hundreds of thousands of years. Such weakly structured stem models explain patterns of polymorphism that had previously been attributed to contributions from archaic hominins in Africa2-7. In contrast to models with archaic introgression, we predict that fossil remains from coexisting ancestral populations should be genetically and morphologically similar, and that only an inferred 1-4% of genetic differentiation among contemporary human populations can be attributed to genetic drift between stem populations. We show that model misspecification explains the variation in previous estimates of divergence times, and argue that studying a range of models is key to making robust inferences about deep history.

Neutrophil extracellular trap formation and gene programs distinguish TST/IGRA sensitization outcomes among Mycobacterium tuberculosis exposed persons living with HIV.

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

Genome-wide association study of resistance to Mycobacterium tuberculosis infection identifies a locus at 10q26.2 in three distinct populations.

The natural history of tuberculosis (TB) is characterized by a large inter-individual outcome variability after exposure to Mycobacterium tuberculosis. Specifically, some highly exposed individuals remain resistant to M. tuberculosis infection, as inferred by tuberculin skin test (TST) or interferon-gamma release assays (IGRAs). We performed a genome-wide association study of resistance to M. tuberculosis infection in an endemic region of Southern Vietnam. We enrolled household contacts (HHC) of pulmonary TB cases and compared subjects who were negative for both TST and IGRA (n = 185) with infected individuals (n = 353) who were either positive for both TST and IGRA or had a diagnosis of TB. We found a genome-wide significant locus on chromosome 10q26.2 with a cluster of variants associated with strong protection against M. tuberculosis infection (OR = 0.42, 95%CI 0.35-0.49, P = 3.71×10-8, for the genotyped variant rs17155120). The locus was replicated in a French multi-ethnic HHC cohort and a familial admixed cohort from a hyper-endemic area of South Africa, with an overall OR for rs17155120 estimated at 0.50 (95%CI 0.45-0.55, P = 1.26×10-9). The variants are located in intronic regions and upstream of C10orf90, a tumor suppressor gene which encodes an ubiquitin ligase activating the transcription factor p53. In silico analysis showed that the protective alleles were associated with a decreased expression in monocytes of the nearby gene ADAM12 which could lead to an enhanced response of Th17 lymphocytes. Our results reveal a novel locus controlling resistance to M. tuberculosis infection across different populations.

Mycobacterium tuberculosis complex and human coadaptation: a two-way street complicating host susceptibility to TB.

For centuries, the Mycobacterium tuberculosis complex (MTBC) has infected numerous populations, both human and non-human, causing symptomatic tuberculosis (TB) in some hosts. Research investigating the MTBC and how it has evolved with its host over time is sparse and has not resulted in many significant findings. There are even fewer studies investigating adaptation of the human host susceptibility to TB and these have largely focused on genome-wide association and candidate gene association studies. However, results emanating from these association studies are rarely replicated and appear to be population specific. It is, therefore, necessary to relook at the approach taken to investigate the relationship between the MTBC and the human host. Understanding that the evolution of the pathogen is coupled to the evolution of the host might be the missing link needed to effectively investigate their relationship. We hypothesize that this knowledge will bolster future efforts in combating the disease.

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.

Rapid evolution of a skin-lightening allele in southern African KhoeSan.

Skin pigmentation is under strong directional selection in northern European and Asian populations. The indigenous KhoeSan populations of far southern Africa have lighter skin than other sub-Saharan African populations, potentially reflecting local adaptation to a region of Africa with reduced UV radiation. Here, we demonstrate that a canonical Eurasian skin pigmentation gene, SLC24A5, was introduced to southern Africa via recent migration and experienced strong adaptive evolution in the KhoeSan. To reconstruct the evolution of skin pigmentation, we collected phenotypes from over 400 ≠Khomani San and Nama individuals and high-throughput sequenced candidate pigmentation genes. The derived causal allele in SLC24A5, p.Ala111Thr, significantly lightens basal skin pigmentation in the KhoeSan and explains 8 to 15% of phenotypic variance in these populations. The frequency of this allele (33 to 53%) is far greater than expected from colonial period European gene flow; however, the most common derived haplotype is identical among European, eastern African, and KhoeSan individuals. Using four-population demographic simulations with selection, we show that the allele was introduced into the KhoeSan only 2,000 y ago via a back-to-Africa migration and then experienced a selective sweep (s = 0.04 to 0.05 in ≠Khomani and Nama). The SLC24A5 locus is both a rare example of intense, ongoing adaptation in very recent human history, as well as an adaptive gene flow at a pigmentation locus in humans.

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.

Identification of a novel WAS mutation in a South African patient presenting with atypical Wiskott-Aldrich syndrome: a case report.

BACKGROUND: The X-linked recessive primary immunodeficiency disease (PIDD) Wiskott-Aldrich syndrome (WAS) is identified by an extreme susceptibility to infections, eczema and thrombocytopenia with microplatelets. The syndrome, the result of mutations in the WAS gene which encodes the Wiskott-Aldrich protein (WASp), has wide clinical phenotype variation, ranging from classical WAS to X-linked thrombocytopaenia and X-linked neutropaenia. In many cases, the diagnosis of WAS in first affected males is delayed, because patients may not present with the classic signs and symptoms, which may intersect with other thrombocytopenia causes. CASE PRESENTATION: Here, we describe a three-year-old HIV negative boy presenting with recurrent infections, skin rashes, features of autoimmunity and atopy. However, platelets were initially reported as normal in numbers and morphology as were baseline immune investigations. An older male sibling had died in infancy from suspected immunodeficiency. Uncertainty of diagnosis and suspected severe PIDD prompted urgent further molecular investigation. Whole exome sequencing identified c. 397 G > A as a novel hemizygous missense mutation located in exon 4 of WAS. CONCLUSION: With definitive molecular diagnosis, we could target treatment and offer genetic counselling and prenatal diagnostic testing to the family. The identification of novel variants is important to confirm phenotype variations of a syndrome.

The X chromosome and sex-specific effects in infectious disease susceptibility.

The X chromosome and X-linked variants have largely been ignored in genome-wide and candidate association studies of infectious diseases due to the complexity of statistical analysis of the X chromosome. This exclusion is significant, since the X chromosome contains a high density of immune-related genes and regulatory elements that are extensively involved in both the innate and adaptive immune responses. Many diseases present with a clear sex bias, and apart from the influence of sex hormones and socioeconomic and behavioural factors, the X chromosome, X-linked genes and X chromosome inactivation mechanisms contribute to this difference. Females are functional mosaics for X-linked genes due to X chromosome inactivation and this, combined with other X chromosome inactivation mechanisms such as genes that escape silencing and skewed inactivation, could contribute to an immunological advantage for females in many infections. In this review, we discuss the involvement of the X chromosome and X inactivation in immunity and address its role in sexual dimorphism of infectious diseases using tuberculosis susceptibility as an example, in which male sex bias is clear, yet not fully explored.

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.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa.

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2+ HLA-C. The relatively high frequency of C2+ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

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.

Risk alleles for tuberculosis infection associate with reduced immune reactivity in a wild mammalian host.

Integrating biological processes across scales remains a central challenge in disease ecology. Genetic variation drives differences in host immune responses, which, along with environmental factors, generates temporal and spatial infection patterns in natural populations that epidemiologists seek to predict and control. However, genetics and immunology are typically studied in model systems, whereas population-level patterns of infection status and susceptibility are uniquely observable in nature. Despite obvious causal connections, organizational scales from genes to host outcomes to population patterns are rarely linked explicitly. Here we identify two loci near genes involved in macrophage (phagocyte) activation and pathogen degradation that additively increase risk of bovine tuberculosis infection by up to ninefold in wild African buffalo. Furthermore, we observe genotype-specific variation in IL-12 production indicative of variation in macrophage activation. Here, we provide measurable differences in infection resistance at multiple scales by characterizing the genetic and inflammatory variation driving patterns of infection in a wild mammal.

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.

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.

Exome sequencing identifies a novel TTC37 mutation in the first reported case of Trichohepatoenteric syndrome (THE-S) in South Africa.

BACKGROUND: Trichohepatoenteric syndrome (THE-S) or phenotypic diarrhoea of infancy is a rare autosomal recessive disorder characterised by severe infantile diarrhoea, facial dysmorphism, immunodeficiency and woolly hair. It was first described in 1982 in two infants with intractable diarrhoea, liver cirrhosis and abnormal hair structure on microscopy. We report on two siblings from a consanguineous family of Somali descent who, despite extensive clinical investigation, remained undiagnosed until their demise. The index patient died of fulminant cytomegalovirus pneumonitis at 3 months of age. METHODS: Whole exome sequencing (WES) was performed on a premortem DNA sample from the index case. Variants in a homozygous recessive state or compound heterozygous state were prioritized as potential candidate variants using TAPER™. Sanger sequencing was done to genotype the parents, unaffected sibling and a deceased sibling for the variant of interest. RESULTS: Exome sequencing identified a novel homozygous mutation (c.4507C > T, rs200067423) in TTC37 which was confirmed by Sanger sequencing in the index case. The identification of this mutation led to the diagnosis of THE-S in the proband and the same homozygous variant was confirmed in a male sibling who died 4 years earlier with severe chronic diarrhoea of infancy. The unaffected parents and sister were heterozygous for the identified variant. CONCLUSIONS: WES permitted definitive genetic diagnosis despite an atypical presentation in the index case and suggests that severe infection, likely secondary to immunodeficiency, may be a presenting feature. In addition definitive molecular diagnosis allows for genetic counseling and future prenatal diagnosis, and demonstrates the value of WES for post-mortem diagnosis of disorders with a non-specific clinical presentation in which a Mendelian cause is suspected.

Major Loci on Chromosomes 8q and 3q Control Interferon γ Production Triggered by Bacillus Calmette-Guerin and 6-kDa Early Secretory Antigen Target, Respectively, in Various Populations.

BACKGROUND: Interferon γ (IFN-γ) release assays (IGRAs) provide an in vitro measurement of antimycobacterial immunity that is widely used as a test for Mycobacterium tuberculosis infection. IGRA outcomes are highly heritable in various populations, but the nature of the involved genetic factors remains unknown. METHODS: We conducted a genome-wide linkage analysis of IGRA phenotypes in families from a tuberculosis household contact study in France and a replication study in families from South Africa to confirm the loci identified. RESULTS: We identified a major locus on chromosome 8q controlling IFN-γ production in response to stimulation with live bacillus Calmette-Guerin (BCG; LOD score, 3.81; P = 1.40 × 10(-5)). We also detected a second locus, on chromosome 3q, that controlled IFN-γ levels in response to stimulation with 6-kDa early secretory antigen target, when accounting for the IFN-γ production shared with that induced by BCG (LOD score, 3.72; P = 1.8 × 10(-5)). Both loci were replicated in South African families, where tuberculosis is hyperendemic. These loci differ from those previously identified as controlling the response to the tuberculin skin test (TST1 and TST2) and the production of TNF-α (TNF1). CONCLUSIONS: The identification of 2 new linkage signals in populations of various ethnic origins living in different M. tuberculosis exposure settings provides new clues about the genetic control of human antimycobacterial immunity.

The complete genome sequence of the African buffalo (Syncerus caffer).

BACKGROUND: The African buffalo (Syncerus caffer) is an important role player in the savannah ecosystem. It has become a species of relevance because of its role as a wildlife maintenance host for an array of infectious and zoonotic diseases some of which include corridor disease, foot-and-mouth disease and bovine tuberculosis. To date, no complete genome sequence for S. caffer had been available for study and the genomes of other species such as the domestic cow (Bos taurus) had been used as a proxy for any genetics analysis conducted on this species. Here, the high coverage genome sequence of the African buffalo (S. caffer) is presented. RESULTS: A total of 19,765 genes were predicted and 19,296 genes could be successfully annotated to S. caffer while 469 genes remained unannotated. Moreover, in order to extend a detailed annotation of S. caffer, gene clusters were constructed using twelve additional mammalian genomes. The S. caffer genome contains 10,988 gene clusters, of which 62 are shared exclusively between B. taurus and S. caffer. CONCLUSIONS: This study provides a unique genomic perspective for the S. caffer, allowing for the identification of novel variants that may play a role in the natural history and physiological adaptations.