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

Pulmonary tuberculosis epidemiology and genetics in Kazakhstan

Background

Tuberculosis (TB) is a major public health emergency in many countries, including Kazakhstan. Despite the decline in the incidence rate and having one of the highest treatment effectiveness in the world, the incidence rate of TB remains high in Kazakhstan. Social and environmental factors along with host genetics contribute to pulmonary tuberculosis (PTB) incidence. Due to the high incidence rate of TB in Kazakhstan, our research aimed to study the epidemiology and genetics of PTB in Kazakhstan.

Materials and methods

1,555 participants were recruited to the case-control study. The epidemiology data was taken during an interview. Polymorphisms of selected genes were determined by real-time PCR using pre-designed TaqMan probes.

Results

Epidemiological risk factors like diabetes (χ2 = 57.71, p < 0.001), unemployment (χ2 = 81.1, p < 0.001), and underweight-ranged BMI (<18.49, χ2 = 206.39, p < 0.001) were significantly associated with PTB. VDR FokI (rs2228570) and VDR BsmI (rs1544410) polymorphisms were associated with an increased risk of PTB. A/A genotype of the TLR8 gene (rs3764880) showed a significant association with an increased risk of PTB in Asians and Asian males. The G allele of the rs2278589 polymorphism of the MARCO gene increases PTB susceptibility in Asians and Asian females. VDR BsmI (rs1544410) polymorphism was significantly associated with PTB in Asian females. A significant association between VDR ApaI polymorphism and PTB susceptibility in the Caucasian population of Kazakhstan was found.

Conclusion

This is the first study that evaluated the epidemiology and genetics of PTB in Kazakhstan on a relatively large cohort. Social and environmental risk factors play a crucial role in TB incidence in Kazakhstan. Underweight BMI (<18.49 kg/m2), diabetes, and unemployment showed a statistically significant association with PTB in our study group. FokI (rs2228570) and BsmI (rs1544410) polymorphisms of the VDR gene can be used as possible biomarkers of PTB in Asian males. rs2278589 polymorphism of the MARCO gene may act as a potential biomarker of PTB in Kazakhs. BsmI polymorphism of the VDR gene and rs2278589 polymorphism of the MARCO gene can be used as possible biomarkers of PTB risk in Asian females as well as VDR ApaI polymorphism in Caucasians.

Human alveolar macrophage response to Mycobacterium tuberculosis: immune characteristics underlying large inter-individual variability

Background

Mycobacterium tuberculosis (M.tb), the causative bacterium of tuberculosis (TB), establishes residence and grows in human alveolar macrophages (AMs). Inter-individual variation in M.tb-human AM interactions can indicate TB risk and the efficacy of therapies and vaccines; however, we currently lack an understanding of the gene and protein expression programs that dictate this variation in the lungs.

Results

Herein, we systematically analyze interactions of a virulent M.tb strain H37Rv with freshly isolated human AMs from 28 healthy adult donors, measuring host RNA expression and secreted candidate proteins associated with TB pathogenesis over 72h. A large set of genes possessing highly variable inter-individual expression levels are differentially expressed in response to M.tb infection. Eigengene modules link M.tb growth rate with host transcriptional and protein profiles at 24 and 72h. Systems analysis of differential RNA and protein expression identifies a robust network with IL1B, STAT1, and IDO1 as hub genes associated with M.tb growth. RNA time profiles document stimulation towards an M1-type macrophage gene expression followed by emergence of an M2-type profile. Finally, we replicate these results in a cohort from a TB-endemic region, finding a substantial portion of significant differentially expressed genes overlapping between studies.

Conclusions

We observe large inter-individual differences in bacterial uptake and growth, with tenfold variation in M.tb load by 72h.The fine-scale resolution of this work enables the identification of genes and gene networks associated with early M.tb growth dynamics in defined donor clusters, an important step in developing potential biological indicators of individual susceptibility to M.tb infection and response to therapies.

Role of Nischarin in the pathology of diseases: a special emphasis on breast cancer

Nischarin has been demonstrated to have tumor suppressor functions. In this review, we comprehensively discuss up to date information about Nischarin. In addition, this paper aims to report the prognostic value, clinical relevance, and biological significance of the Nischarin gene (NISCH) in breast cancer (BCa) patients using the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) datasets. We evaluated NISCH gene expression and its correlation to patient survival, baseline expression, and expression variation based on age groups, tumor stage, tumor size, tumor grade, and lymph node status in different subtypes of BCa. Since NISCH has been extensively reported to inhibit EMT and cancer cell migration, we also checked for the correlation between NISCH and EMT genes in addition to the correlation between NISCH and cell migration genes. Our results indicate that NISCH is a tumor suppressor that plays a critical role in BCa initiation, progression, and tumor development. We find that there is a higher level of NISCH expression in normal breast tissues compared to breast cancer tissues. Also, aggressive subtypes of breast cancers, such as the triple negative/basal category, have decreased levels of NISCH as the disease progresses. Finally, we report that NISCH is inversely correlated with many EMT and cancer cell migration genes in BCa. Interestingly, we identified a significant negative correlation between NISCH expression and its methylation in breast cancer patients. Overall, the goal of this report is to establish a strong clinical basis for further investigation into the cellular, molecular, and physiological roles of NISCH in BCa. Ultimately, NISCH gene expression might be clinically harnessed as a biomarker or predictor of invasiveness and metastasis in BCa.

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.

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.

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

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

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.

Challenges of Immune Response Diversity in the Human Population Concerning New Tuberculosis Diagnostics, Therapies, and Vaccines

Universal approaches to the prevention and treatment of human diseases fail to take into account profound immune diversity resulting from genetic variations across populations. Personalized or precision medicine takes into account individual lifestyle, environment, and biology (genetics and immune status) and is being adopted in several disease intervention strategies such as cancer and heart disease. However, its application in infectious diseases, particularly global diseases such as tuberculosis (TB), is far more complex and in a state of infancy. Here, we discuss the impact of human genetic variations on immune responses and how they relate to failures seen in current TB diagnostic, therapy, and vaccine approaches across populations. We offer our perspective on the challenges and potential for more refined approaches going forward.

Single Nucleotide Polymorphisms in TLR4 Affect Susceptibility to Tuberculosis in Mexican Population from the State of Veracruz

Tuberculosis is still a global public health problem, with an estimated 10 million new cases and 1.6 million deaths in 2017. Of all humans infected with M. tuberculosis, only 10-15% will develop active tuberculosis disease during their lifetime, and data suggest that along with environmental factors, genetic factors influence susceptibility to develop active disease. Toll-like receptors (TLRs) are pattern recognition receptors that play a central role in the initiation and shaping of adaptive immune responses, and several TLRs have been shown to recognize mycobacterial components. In this work, we performed a case-control study to determine if common single nucleotide polymorphisms (SNPs) in genes encoding TLRs 1, 2, 4, 6, and 10 are associated with susceptibility to develop active tuberculosis in population from the state of Veracruz, Mexico. The study included 279 cases and 569 controls. The results show that the frequency of two SNPs in TLR4 was significantly higher in controls than in tuberculosis patients. The minor allele (G) of rs4986790 in TLR4 (D299G) decreased the risk of active tuberculosis in the allelic (A vs. G, OR = 0.31, 95%CI = 0.09‐0.81, p = 0.01) and in the dominant genetic model (AA vs. GG+AG, OR = 0.26, 95%CI = 0.09‐0.77, p = 0.02). Similarly, the minor allele (T) of rs4986791 in TLR4 (T399I) decreased the risk of active disease in the allelic model (C vs. T, OR = 0.29, 95%CI = 0.10‐0.90, p = 0.03). We did not find an association of SNPs in TLR1 (N248S), TLR2 (R753Q), TLR6 (S249P), and TLR10 (A153S and V298I) with tuberculosis disease. These results suggest that in this population, genetic variants of TLR4 affect the susceptibility for suffering active tuberculosis disease.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Simple Summary

Bovine tuberculosis is an infectious disease of cattle caused by Mycobacterium bovis characterized by the formation of tubercles in any organ or tissue. Bovine tuberculosis represents a significant veterinary and public health problem in many parts of the world. It is zoonotic, transmitted to humans through consumption of infected milk and other cattle products. Although many factors influence infection and progression of the disease, there must be an important host genetic component that explains why some animals get sick and others remain healty. We present evidence of genetic variants associated with resistance to tuberculosis in Mexican Holstein dairy cattle using a case-control approach with a selective DNA pooling. Here, we identified novel quantitative trait loci regions harboring genes involved in Mycobacterium spp. immune response. This is a first screening about resistance to tuberculosis infection on Mexican Holstein cattle based on a dense single nucleotide polymorphism chip. The identified genes belong to both, the already known, and the undisclosed quantitative trait loci regions.

Abstract

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Polymorphisms in Toll-Like Receptor 10 and Tuberculosis Susceptibility: Evidence from Three Independent Series

Background

The toll-like receptor 2 (TLR2)-mediated immune response is critical for host defense against Mycobacterium tuberculosis. There is evidence that TLR10, a TLR2 signaling modulator, may be involved in progression of tuberculosis (TB).

Methods

Using a self-validating case–control design, we tested for an association between seven TLR10 polymorphisms and susceptibility to TB in three independent series with two distinct populations. Single-nucleotide polymorphism (SNP) genotypes were determined by the SNPscan^TM method. Three genetic models (additive, dominant, and recessive) as well as multiple-SNP score analyses were used to evaluate the risk of TB associated with the TLR10 SNPs.

Results

By comparing TB patients with healthy controls, we observed two SNPs (rs11466617 and rs4129009) that were associated with decreased risk of TB in the Tibetan population, but did not in the Chinese Han population. Further analysis demonstrated that the rs11466617 Chengdu cohort genotype served as a protective factor against the progression of latent TB infection (LTBI) to active TB under the recessive model. None of the SNPs were significantly different in comparisons of TB-uninfected people with LTBI individuals. Additionally, when the underlying four TB-associated loci were considered together in a multiple-SNP score analysis, we observed an allele dose-dependent decrease in TB risk in Tibetans.

Conclusion

Variants of TLR10 may show an ethnic specificity on susceptibility to TB in Tibetan individuals. rs11466617 affected the susceptibility to progress from LTBI to active TB disease, but was not associated with the establishment of LTBI after M. tuberculosis exposure. More studies are needed to verify this genetic epidemiological result and unravel the role of TLR10 SNPs in the pathogenesis of TB.

Genomics of human pulmonary tuberculosis: from genes to pathways

Purpose of review

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), remains a major public health threat globally. Several lines of evidence support a role for host genetic factors in resistance/susceptibility to TB disease and MTB infection. However, results across candidate gene and genome-wide association studies (GWAS) are largely inconsistent, so a cohesive genetic model underlying TB risk has not emerged.

Recent Findings

Despite the difficulties in identifying consistent genetic associations, genetic studies of TB and MTB infection have revealed a few well-documented loci. These well validated genes are presented in this review, but there remains a large gap in how these genes translate into better understanding of TB. To address this, we present a pathway based extension of standard association analyses, seeding the results with the best validated genes from candidate gene and GWAS studies.

Summary

Several pathways were significantly enriched using pathway analyses that may help to explain population patterns of TB risk. In conclusion, we advocate for novel approaches to the study of host genetic analysis of TB that extend traditional association approaches.

Human Genomics of Mycobacterium tuberculosis Infection and Disease

Purpose of review

The study of the genetic basis of tuberculosis pathogenesis has benefited from powerful technological innovations, a more structured definition of latent and clinical manifestations of the disease, and the application of functional genomics approaches. This short review aims to summarize recent advances and to provide a link with results of previous human genetic studies of tuberculosis susceptibility.

Recent findings

Transcriptomics has been shown to be a useful tool to predict progression from latency to clinical disease while functional genomics has traced the molecular events that link pathogen-triggered gene expression and host genetics. Resistance to infection with Mycobacterium tuberculosis has been revealed to be strongly impacted by host genetics. Host genomics of clinical disease has been shown to be most powerful when focusing on carefully selected clinical entities and possibly by considering host pathogen combinations.

Summary

Future studies need to build on the latest molecular findings to define disease subtypes to successfully elucidate the human genetic component in tuberculosis pathogenesis.

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.