Latent Tuberculosis Infection Is Associated with an Enrichment of Short-Chain Fatty Acid-Producing Bacteria in the Stool of Women Living with HIV.

Latent tuberculosis infection (LTBI) is common in people living with HIV (PLHIV) in high-TB-burden settings. Active TB is associated with specific stool taxa; however, little is known about the stool microbiota and LTBI in PLHIV. We characterised the stool microbiota of PLHIV with [interferon-γ release assay (IGRA)- and tuberculin skin test (TST)-positive] or without (IGRA- and TST-negative) LTBI (n = 25 per group). The 16S rRNA DNA sequences were analysed using QIIME2, Dirichlet-Multinomial Mixtures, DESeq2, and PICRUSt2. No α- or ß-diversity differences occurred by LTBI status; however, LTBI-positive people were Faecalibacterium-, Blautia-, Gemmiger-, and Bacteroides-enriched and Moryella-, Atopobium-, Corynebacterium-, and Streptococcus-depleted. Inferred metagenome data showed that LTBI-negative-enriched pathways included several metabolite degradation pathways. Stool from LTBI-positive people demonstrated differential taxa abundance based on a quantitative response to antigen stimulation. In LTBI-positive people, older people had different ß-diversities than younger people, whereas in LTBI-negative people, no differences occurred across age groups. Amongst female PLHIV, those with LTBI were, vs. those without LTBI, Faecalibacterium-, Blautia-, Gemmiger-, and Bacteriodes-enriched, which are producers of short-chain fatty acids. Taxonomic differences amongst people with LTBI occurred according to quantitative response to antigen stimulation and age. These data enhance our understanding of the microbiome's potential role in LTBI.

Latent tuberculosis infection is associated with an enrichment of short chain fatty acid producing bacteria in the stool of women living with HIV.

Background: Latent tuberculosis infection (LTBI) is common in people living with HIV (PLHIV) in high TB burden settings. Active TB is associated with specific stool taxa; however, little is known about the stool microbiota and LTBI, including in PLHIV. Method: Within a parent study that recruited adult females with HIV from Cape Town, South Africa into predefined age categories (18-25, 35-60 years), we characterised the stool microbiota of those with [interferon-γ release assay (IGRA)- and tuberculin skin test (TST)-positive] or without (IGRA- and TST- negative) LTBI (n=25 per group). 16S rRNA DNA sequences were analysed using QIIME2, Dirichlet Multinomial Mixtures, DESeq2 and PICRUSt2. Results: No α- or ß-diversity differences occurred by LTBI status; however, LTBI-positives were Faecalibacterium-, Blautia-, Gemmiger-, Bacteroides-enriched and Moryella-, Atopobium-, Corynebacterium-, Streptococcus-depleted. Inferred metagenome data showed LTBI-negative-enriched pathways included several involved in methylglyoxal degradation, L-arginine, putrescine, 4-aminobutanoate degradation and L-arginine and ornithine degradation. Stool from LTBI-positives demonstrated differential taxa abundance based on a quantitative response to antigen stimulation (Acidaminococcus-enrichment and Megamonas-, Alistipes-, and Paraprevotella-depletion associated with higher IGRA or TST responses, respectively). In LTBI-positives, older people had different ß-diversities than younger people whereas, in LTBI-negatives, no differences occurred across age groups. Conclusion: Amongst female PLHIV, those with LTBI had, vs. those without LTBI, Faecalibacterium, Blautia, Gemmiger, Bacteriodes-enriched, which are producers of short chain fatty acids. Taxonomic differences amongst people with LTBI occurred according to quantitative response to antigen stimulation and age. These data enhance our understanding of the microbiome's potential role in LTBI.

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.

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.

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.

Deciphering Genetic Susceptibility to Tuberculous Meningitis.

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

Local Ancestry Adjusted Allelic Association Analysis Robustly Captures Tuberculosis Susceptibility Loci.

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

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.

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.

An observational study identifying highly tuberculosis-exposed, HIV-1-positive but persistently TB, tuberculin and IGRA negative persons with M. tuberculosis specific antibodies in Cape Town, South Africa.

BACKGROUND: Mycobacterium tuberculosis (Mtb) infection is inferred from positive results of T-cell immune conversion assays measuring Mtb-specific interferon gamma production or tuberculin skin test (TST) reactivity. Certain exposed individuals do not display T-cell immune conversion in these assays and do not develop TB. Here we report a hitherto unknown form of this phenotype: HIV-1-positive persistently TB, tuberculin and IGRA negative (HITTIN). METHODS: A community-based case-control design was used to systematically screen and identify adults living with HIV (HIV+), aged 35-60 years, who met stringent study criteria, and then longitudinally followed up for repeat IGRA and TST testing. Participants had no history of TB despite living in TB hyper-endemic environments in Cape Town, South Africa with a provincial incidence of 681/100,000. Mtb-specific antibodies were measured using ELISA and Luminex. FINDINGS: We identified 48/286 (17%) individuals who tested persistently negative for Mtb-specific T-cell immunoreactivity (three negative Quantiferon results and one TST = 0mm) over 206±154 days on average. Of these, 97·2% had documented CD4 counts<200 prior to antiretroviral therapy (ART). They had received ART for 7·0±3·0 years with a latest CD4 count of 505·8±191·4 cells/mm3. All HITTIN sent for further antibody testing (n=38) displayed  Mtb-specific antibody titres. INTERPRETATION: Immune reconstituted HIV+ persons can be persistently non-immunoreactive to TST and interferon-γ T-cell responses to Mtb, yet develop species-specific antibody responses. Exposure is evidenced by Mtb-specific antibody titres. Our identification of HIV+ individuals displaying a persisting lack of response to TST and IGRA T-cell immune conversion paves the way for future studies to investigate this phenotype in the context of HIV-infection that so far have received only scant attention.

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.

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.

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.

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.

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.

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.