Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease

Humoral pathways of innate immune regulation in granuloma formation

The humoral arm of mammalian innate immunity regulates several molecular mechanisms involved in resistance to pathogens, inflammation, and tissue repair. Recent studies highlight the crucial role played by humoral mediators in granulomatous inflammation. However the molecular mechanisms linking the function of these soluble molecules to the initiation and maintenance of granulomas remain elusive. We propose that humoral innate immunity coordinates fundamental physiological processes in macrophages which, in turn, initiate activation and transformation events that enable granuloma formation. We discuss the involvement of humoral mediators in processes such as immune activation, phagocytosis, metabolism, and tissue remodeling, and how these can dictate macrophage functionality during granuloma formation. These advances present opportunities for discovering novel disease factors and developing targeted, more effective treatments for granulomatous diseases.

Genetic Diversity of Mycobacterium tuberculosis Strains Isolated from HIV-Infected Patients in Mexico

There has been very limited investigation regarding the genetic diversity of Mycobacterium tuberculosis (MTb) strains isolated from human immunodeficiency virus (HIV)-infected patients in Mexico. In this study, we isolated 93 MTb strains from pulmonary and extrapulmonary samples of HIV-infected patients treated in a public hospital in Mexico City to evaluate the genetic diversity using spoligotyping and mycobacterial interspersed repetitive unit-variable-number tandem-repeat (MIRU-VNTR) typing (based on 24 loci). The cohort comprised 80 male and 13 female individuals. There was a positive correlation between a high HIV viral load (>100,000 copies) and extrapulmonary tuberculosis (TB) (r = 0.306, p = 0.008). Lineage 4 was the most frequent lineage (79 strains). In this lineage, we found the H clade (n = 24), including the Haarlem, H3, and H1 families; the T clade (n = 22), including T1 and T2; the X clade (n = 15), including X1 and X3; the LAM clade (n = 14), including LAM1, LAM2, LAM3, LAM6, and LAM9; the S clade (n = 2); Uganda (n = 1); and Ghana (n = 1). We also found 12 strains in the EAI clade belonging to lineage 1, including the EAI2-Manila and EAI5 families. Interestingly, we identified one strain belonging to the Beijing family, which is part of lineage 2. One strain could not be identified. This study reports high genetic diversity among MTb strains, highlighting the need for a molecular epidemiological surveillance system that can help to monitor the spread of these strains, leading to more appropriate measures for TB control in HIV-infected patients.

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.

Genome-wide association and functional genomic analyses for various hoof health traits in North American Holstein cattle

Hoof diseases are a major welfare and economic issue in the global dairy cattle production industry, which can be minimized through improved management and breeding practices. Optimal genetic improvement of hoof health could benefit from a deep understanding of the genetic background and biological underpinning of indicators of hoof health. Therefore, the primary objectives of this study were to perform genome-wide association studies, using imputed high-density genetic markers data from North American Holstein cattle, for 8 hoof-related traits: digital dermatitis, sole ulcer, sole hemorrhage, white line lesion, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, and toe ulcer, and a hoof health index. De-regressed estimated breeding values from 25,580 Holstein animals were used as pseudo-phenotypes for the association analyses. The genomic quality control, genotype phasing, and genotype imputation were performed using the PLINK (version 1.9), Eagle (version 2.4.1), and Minimac4 software, respectively. The functional genomic analyses were performed using the GALLO R package and the DAVID platform. We identified 22, 34, 14, 22, 28, 33, 24, 43, and 15 significant markers for digital dermatitis, heel horn erosion, interdigital dermatitis, interdigital hyperplasia, sole hemorrhage, sole ulcer, toe ulcer, white line lesion disease, and the hoof health index, respectively. The significant markers were located across all autosomes, except BTA10, BTA12, BTA20, BTA26, BTA27, and BTA28. Moreover, the genomic regions identified overlap with various previously reported quantitative trait loci for exterior, health, meat and carcass, milk, production, and reproduction traits. The enrichment analyses identified 44 significant gene ontology terms. These enriched genomic regions harbor various candidate genes previously associated with bone development, metabolism, and infectious and immunological diseases. These findings indicate that hoof health traits are highly polygenic and influenced by a wide range of biological processes.

An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

Multi-platform whole genome sequencing for tuberculosis clinical and surveillance applications

Whole genome sequencing (WGS) of Mycobacterium tuberculosis offers valuable insights for tuberculosis (TB) control. High throughput platforms like Illumina and Oxford Nanopore Technology (ONT) are increasingly used globally, although ONT is known for higher error rates and is less established for genomic studies. Here we present a study comparing the sequencing outputs of both Illumina and ONT platforms, analysing DNA from 59 clinical isolates in highly endemic TB regions of Thailand. The resulting sequence data were used to profile the M. tuberculosis pairs for their lineage, drug resistance and presence in transmission chains, and were compared to publicly available WGS data from Thailand (n = 1456). Our results revealed isolates that are predominantly from lineages 1 and 2, with consistent drug resistance profiles, including six multidrug-resistant strains; however, analysis of ONT data showed longer phylogenetic branches, emphasising the technologies higher error rate. An analysis incorporating the larger dataset identified fifteen of our samples within six potential transmission clusters, including a significant clade of 41 multi-drug resistant isolates. ONT's extended sequences also revealed strain-specific structural variants in pe/ppe genes (e.g. ppe50), which are candidate loci for vaccine development. Despite some limitations, our results show that ONT sequencing is a promising approach for TB genomic research, supporting precision medicine and decision-making in areas with less developed infrastructure, which is crucial for tackling the disease's significant regional burden.

Did the human genome project affect research on Schizophrenia?

The Human Genome Project was undertaken primarily to discover genetic causes and better treatments for human diseases. Schizophrenia was targeted since three of the project`s principal architects had a personal interest and also because, based on family, adoption, and twin studies, schizophrenia was widely believed to be a genetic disorder. Extensive studies using linkage analysis, candidate genes, genome wide association studies [GWAS], copy number variants, exome sequencing and other approaches have failed to identify causal genes. Instead, they identified almost 300 single nucleotide polymorphisms [SNPs] associated with altered risks of developing schizophrenia as well as some rare variants associated with increased risk in a small number of individuals. Risk genes play a role in the clinical expression of most diseases but do not cause the disease in the absence of other factors. Increasingly, observers question whether schizophrenia is strictly a genetic disorder. Beginning in 1996 NIMH began shifting its research resources from clinical studies to basic research based on the promise of the Human Genome Project. Consequently, three decades later NIMH's genetics investment has yielded almost nothing clinically useful for individuals currently affected. It is time to review NIMH`s schizophrenia research portfolio.

Assessing the association between the circulating levels of inflammatory cytokines and the risk of tuberculosis: A bidirectional two-sample mendelian randomization study

BACKGROUND

Numerous observational studies have previously reported an association between inflammatory cytokines and tuberculosis (TB). However, the causal relationship between these factors remains unclear. Consequently, we conducted two-sample Mendelian randomization (MR) analyses to ascertain the causal link between levels of inflammatory cytokines and the risk of TB.

METHODS

Single nucleotide polymorphisms (SNPs) robustly associated with the cytokines, located in or close to their coding gene. SNP was obtained from genome-wide association studies (GWAS) of 8293 individuals of Finnish. TB data was obtained from the UK Biobank, which included 46,293 individuals of European ancestry (comprising 2277 TB cases and 46,056 controls). Two-sample, bi-directional MR analyses using inverse-variance weighted (IVW) method as the primary analysis. Followed by comprehensive sensitivity analyses to validate the robustness of results.

RESULT

The study showed that the causal relationship between circulating levels of interleukin (IL)-7 and risk of TB (odds ratio [OR] = 1.001, 95% confidence intervals [CIs]: 1.000, 1.003. p = 0.047). No causal associations were observed between other influencing factors and the occurrence of TB. Furthermore, the analysis revealed that TB infection exhibited negative causal associations with macrophage inflammatory protein 1 alpha ([MIP-1α], OR = 0.007, 95% CI: 0.000, 0.192. p = 0.004), IL-2 (OR = 0.014, 95% CI: 0.010, 0.427. p = 0.014), interleukin-2 receptor alpha chain([IL-2rα], OR = 0.019, 95% CI: 0.001, 0.525. p = 0.019) and basic fibroblast growth factor ([bFGF], OR = 0.066, 95% CI: 0.006, 0.700. p = 0.024).

CONCLUSION

The study has illuminated the causal link between inflammatory cytokines and TB, thereby enhancing our comprehension of the potential mechanisms underlying TB pathogenesis. This discovery offers promising avenues for the identification of novel therapeutic targets in TB treatment. These insights may ultimately pave the way for more effective treatment approaches, thereby improving patient outcomes.

BurkHostGEN: a study protocol for evaluating variations in the Burkholderia pseudomallei and host genomes associated with melioidosis infection

Background

Melioidosis is a frequently fatal disease caused by an environmental bacterium Burkholderia pseudomallei. The disease is prevalent in northeast Thailand, particularly among rice field farmers who are at risk of bacterial exposure through contact with contaminated soil and water. However, not all exposure results in disease, and infection can manifest diverse outcomes. We postulate that genetic factors, whether from the bacterium, the host or the combination of both, may influence disease outcomes. To address this hypothesis, we aim to collect, sequence, and analyse genetic data from melioidosis patients and controls, along with isolates of B. pseudomallei obtained from patients. Additionally, we will study the metagenomics of the household water supply for both patients and controls, including the presence of B. pseudomallei.

Methods

BurkHostGEN is an ongoing observational study being conducted at Sunpasitthiprasong Hospital, Ubon Ratchathani, Thailand. We are obtaining consent from 600 melioidosis patients and 700 controls, spanning both sexes, to collect 1 mL of blood for host DNA analysis, 3 mL of blood for RNA analysis, as well as 5 L of household water supply for metagenomic analysis. Additionally, we are isolating B. pseudomallei from the melioidosis patients to obtain bacterial DNA. This comprehensive approach will allow us to identify B. pseudomallei and their paired host genetic factors associated with disease acquisition and severity. Ethical approvals have been obtained for BurkHostGEN. Host and bacterial genetic data will be uploaded to European Genome-Phenome Archive (EGA) and European Nucleotide Archive (ENA), respectively.

Conclusions

BurkHostGEN holds the potential to discover bacterial and host genetic factors associated with melioidosis infection and severity of illness. It can also support various study designs, including biomarker validation, disease pathogenesis, and epidemiological analysis not only for melioidosis but also for other infectious diseases.

G2GSnake: a Snakemake workflow for host-pathogen genomic association studies

Summary

Joint analyses of paired host and pathogen genome sequences have the potential to enhance our understanding of host-pathogen interactions. A systematic approach to conduct such a joint analysis is through a "genome-to-genome" (G2G) association study, which involves testing for associations between all host and pathogen genetic variants. Significant associations reveal host genetic factors that might drive pathogen variation, highlighting biological mechanisms likely to be involved in host control and pathogen escape. Here, we present a Snakemake workflow that allows researchers to conduct G2G studies in a reproducible and scalable manner. In addition, we have developed an intuitive R Shiny application that generates custom summaries of the results, enabling users to derive relevant insights.

Availability and implementation

G2GSnake is freely available at: https://github.com/zmx21/G2GSnake under the MIT license.

Relationship between ambient temperature at sampling and the interferon gamma test result for bovine tuberculosis in cattle

Bovine tuberculosis (bTB) is a disease of significant economic and zoonotic importance, therefore, optimising tests for the identification of Mycobacterium bovis infected cattle is essential. The Interferon Gamma (IFN-γ) Release Assay (IGRA) can diagnose M. bovis infected cattle at an early stage, is easy to perform and can be used alongside skin tests for confirmatory purposes or to increase diagnostic sensitivity. It is known that IGRA performance is sensitive to environmental conditions under which samples are taken and transported. In this study, the association between the ambient temperature on the day of bleeding and the subsequent IGRA result for bTB was quantified using field samples from Northern Ireland (NI). Results of 106,434 IGRA results (2013-2018) were associated with temperature data extracted from weather stations near tested cattle herds. Model dependent variables were the levels of IFN-γ triggered by avian purified protein derivative (PPDa), M. bovis PPD (PPDb), their difference (PPD(b-a)) as well as the final binary outcome (positive or negative for M. bovis infection). IFN-γ levels after both PPDa and PPDb stimulation were lowest at the extremes of the temperature distribution for NI. The highest IGRA positive probability (above 6%) was found on days with moderate maximum temperatures (6-16 °C) or moderate minimum temperatures (4-7 °C). Adjustment for covariates did not lead to major changes in the model estimates. These data suggest that IGRA performance can be affected when samples are taken at high or low temperatures. Whilst it is difficult to exclude physiological factors, the data nonetheless supports the temperature control of samples from bleeding through to laboratory to help mitigate post-collection confounders.