Human SNP links differential outcomes in inflammatory and infectious disease to a FOXO3-regulated pathway.

The clinical course and eventual outcome, or prognosis, of complex diseases varies enormously between affected individuals. This variability critically determines the impact a disease has on a patient's life but is very poorly understood. Here, we exploit existing genome-wide association study data to gain insight into the role of genetics in prognosis. We identify a noncoding polymorphism in FOXO3A (rs12212067: T > G) at which the minor (G) allele, despite not being associated with disease susceptibility, is associated with a milder course of Crohn's disease and rheumatoid arthritis and with increased risk of severe malaria. Minor allele carriage is shown to limit inflammatory responses in monocytes via a FOXO3-driven pathway, which through TGFß1 reduces production of proinflammatory cytokines, including TNFα, and increases production of anti-inflammatory cytokines, including IL-10. Thus, we uncover a shared genetic contribution to prognosis in distinct diseases that operates via a FOXO3-driven pathway modulating inflammatory responses.

Host genotype-specific therapies can optimize the inflammatory response to mycobacterial infections.

Susceptibility to tuberculosis is historically ascribed to an inadequate immune response that fails to control infecting mycobacteria. In zebrafish, we find that susceptibility to Mycobacterium marinum can result from either inadequate or excessive acute inflammation. Modulation of the leukotriene A(4) hydrolase (LTA4H) locus, which controls the balance of pro- and anti-inflammatory eicosanoids, reveals two distinct molecular routes to mycobacterial susceptibility converging on dysregulated TNF levels: inadequate inflammation caused by excess lipoxins and hyperinflammation driven by excess leukotriene B(4). We identify therapies that specifically target each of these extremes. In humans, we identify a single nucleotide polymorphism in the LTA4H promoter that regulates its transcriptional activity. In tuberculous meningitis, the polymorphism is associated with inflammatory cell recruitment, patient survival and response to adjunctive anti-inflammatory therapy. Together, our findings suggest that host-directed therapies tailored to patient LTA4H genotypes may counter detrimental effects of either extreme of inflammation.

The lta4h locus modulates susceptibility to mycobacterial infection in zebrafish and humans.

Exposure to Mycobacterium tuberculosis produces varied early outcomes, ranging from resistance to infection to progressive disease. Here we report results from a forward genetic screen in zebrafish larvae that identify multiple mutant classes with distinct patterns of innate susceptibility to Mycobacterium marinum. A hypersusceptible mutant maps to the lta4h locus encoding leukotriene A(4) hydrolase, which catalyzes the final step in the synthesis of leukotriene B(4) (LTB(4)), a potent chemoattractant and proinflammatory eicosanoid. lta4h mutations confer hypersusceptibility independent of LTB(4) reduction, by redirecting eicosanoid substrates to anti-inflammatory lipoxins. The resultant anti-inflammatory state permits increased mycobacterial proliferation by limiting production of tumor necrosis factor. In humans, we find that protection from both tuberculosis and multibacillary leprosy is associated with heterozygosity for LTA4H polymorphisms that have previously been correlated with differential LTB(4) production. Our results suggest conserved roles for balanced eicosanoid production in vertebrate resistance to mycobacterial infection.

REL and BHLHE40 Variants Are Associated with IL-12 and IL-10 Responses and Tuberculosis Risk.

The major human genes regulating Mycobacterium tuberculosis-induced immune responses and tuberculosis (TB) susceptibility are poorly understood. Although IL-12 and IL-10 are critical for TB pathogenesis, the genetic factors that regulate their expression in humans are unknown. CNBP, REL, and BHLHE40 are master regulators of IL-12 and IL-10 signaling. We hypothesized that common variants in CNBP, REL, and BHLHE40 were associated with IL-12 and IL-10 production from dendritic cells, and that these variants also influence adaptive immune responses to bacillus Calmette-Guérin (BCG) vaccination and TB susceptibility. We characterized the association between common variants in CNBP, REL, and BHLHE40, innate immune responses in dendritic cells and monocyte-derived macrophages, BCG-specific T cell responses, and susceptibility to pediatric and adult TB in human populations. BHLHE40 single-nucleotide polymorphism (SNP) rs4496464 was associated with increased BHLHE40 expression in monocyte-derived macrophages and increased IL-10 from peripheral blood dendritic cells and monocyte-derived macrophages after LPS and TB whole-cell lysate stimulation. SNP BHLHE40 rs11130215, in linkage disequilibrium with rs4496464, was associated with increased BCG-specific IL-2+CD4+ T cell responses and decreased risk for pediatric TB in South Africa. SNPs REL rs842634 and rs842618 were associated with increased IL-12 production from dendritic cells, and SNP REL rs842618 was associated with increased risk for TB meningitis. In summary, we found that genetic variations in REL and BHLHE40 are associated with IL-12 and IL-10 cytokine responses and TB clinical outcomes. Common human genetic regulation of well-defined intermediate cellular traits provides insights into mechanisms of TB pathogenesis.

MUC5AC Genetic Variation Is Associated With Tuberculous Meningitis Cerebral Spinal Fluid Cytokine Responses and Mortality.

BACKGROUND: The purpose of this study was to assess if single nucleotide polymorphisms (SNPs) in lung mucins MUC5B and MUC5AC are associated with Mycobacterium tuberculosis outcomes. METHODS: Independent SNPs in MUC5B and MUC5AC (genotyped by Illumina HumanOmniExpress array) were assessed for associations with tumor necrosis factor (TNF) concentrations (measured by immunoassay) in cerebral spinal fluid (CSF) from tuberculous meningitis (TBM) patients. SNPs associated with CSF TNF concentrations were carried forward for analyses of pulmonary and meningeal tuberculosis susceptibility and TBM mortality. RESULTS: MUC5AC SNP rs28737416 T allele was associated with lower CSF concentrations of TNF (P = 1.8 × 10-8) and IFN-γ (P = 2.3 × 10-6). In an additive genetic model, rs28737416 T/T genotype was associated with higher susceptibility to TBM (odds ratio [OR], 1.24; 95% confidence interval [CI], 1.03-1.49; P = .02), but not pulmonary tuberculosis (OR, 1.11, 95% CI, .98-1.25; P = .10). TBM mortality was higher among participants with the rs28737416 T/T and T/C genotypes (35/119, 30.4%) versus the C/C genotype (11/89, 12.4%; log-rank P = .005) in a Vietnam discovery cohort (n = 210), an independent Vietnam validation cohort (n = 87; 9/87, 19.1% vs 1/20, 2.5%; log-rank P = .02), and an Indonesia validation cohort (n = 468, 127/287, 44.3% vs 65/181, 35.9%; log-rank P = .06). CONCLUSIONS: MUC5AC variants may contribute to immune changes that influence TBM outcomes.

A Bayesian approach for estimating typhoid fever incidence from large-scale facility-based passive surveillance data.

Decisions about typhoid fever prevention and control are based on estimates of typhoid incidence and their uncertainty. Lack of specific clinical diagnostic criteria, poorly sensitive diagnostic tests, and scarcity of accurate and complete datasets contribute to difficulties in calculating age-specific population-level typhoid incidence. Using data from the Strategic Typhoid Alliance across Africa and Asia program, we integrated demographic censuses, healthcare utilization surveys, facility-based surveillance, and serological surveillance from Malawi, Nepal, and Bangladesh to account for under-detection of cases. We developed a Bayesian approach that adjusts the count of reported blood-culture-positive cases for blood culture detection, blood culture collection, and healthcare seeking-and how these factors vary by age-while combining information from prior published studies. We validated the model using simulated data. The ratio of observed to adjusted incidence rates was 7.7 (95% credible interval [CrI]: 6.0-12.4) in Malawi, 14.4 (95% CrI: 9.3-24.9) in Nepal, and 7.0 (95% CrI: 5.6-9.2) in Bangladesh. The probability of blood culture collection led to the largest adjustment in Malawi, while the probability of seeking healthcare contributed the most in Nepal and Bangladesh; adjustment factors varied by age. Adjusted incidence rates were within or below the seroincidence rate limits of typhoid infection. Estimates of blood-culture-confirmed typhoid fever without these adjustments results in considerable underestimation of the true incidence of typhoid fever. Our approach allows each phase of the reporting process to be synthesized to estimate the adjusted incidence of typhoid fever while correctly characterizing uncertainty, which can inform decision-making for typhoid prevention and control.

Sources of Multidrug Resistance in Patients With Previous Isoniazid-Resistant Tuberculosis Identified Using Whole Genome Sequencing: A Longitudinal Cohort Study.

BACKGROUND: Meta-analysis of patients with isoniazid-resistant tuberculosis (TB) given standard first-line anti-TB treatment indicated an increased risk of multidrug-resistant TB (MDR-TB) emerging (8%), compared to drug-sensitive TB (0.3%). Here we use whole genome sequencing (WGS) to investigate whether treatment of patients with preexisting isoniazid-resistant disease with first-line anti-TB therapy risks selecting for rifampicin resistance, and hence MDR-TB. METHODS: Patients with isoniazid-resistant pulmonary TB were recruited and followed up for 24 months. Drug susceptibility testing was performed by microscopic observation drug susceptibility assay, mycobacterial growth indicator tube, and by WGS on isolates at first presentation and in the case of re-presentation. Where MDR-TB was diagnosed, WGS was used to determine the genomic relatedness between initial and subsequent isolates. De novo emergence of MDR-TB was assumed where the genomic distance was 5 or fewer single-nucleotide polymorphisms (SNPs), whereas reinfection with a different MDR-TB strain was assumed where the distance was 10 or more SNPs. RESULTS: Two hundred thirty-nine patients with isoniazid-resistant pulmonary TB were recruited. Fourteen (14/239 [5.9%]) patients were diagnosed with a second episode of TB that was multidrug resistant. Six (6/239 [2.5%]) were identified as having evolved MDR-TB de novo and 6 as having been reinfected with a different strain. In 2 cases, the genomic distance was between 5 and 10 SNPs and therefore indeterminate. CONCLUSIONS: In isoniazid-resistant TB, de novo emergence and reinfection of MDR-TB strains equally contributed to MDR development. Early diagnosis and optimal treatment of isoniazid-resistant TB are urgently needed to avert the de novo emergence of MDR-TB during treatment.

Human genetic variation in VAC14 regulates Salmonella invasion and typhoid fever through modulation of cholesterol.

Risk, severity, and outcome of infection depend on the interplay of pathogen virulence and host susceptibility. Systematic identification of genetic susceptibility to infection is being undertaken through genome-wide association studies, but how to expeditiously move from genetic differences to functional mechanisms is unclear. Here, we use genetic association of molecular, cellular, and human disease traits and experimental validation to demonstrate that genetic variation affects expression of VAC14, a phosphoinositide-regulating protein, to influence susceptibility to Salmonella enterica serovar Typhi (S Typhi) infection. Decreased VAC14 expression increased plasma membrane cholesterol, facilitating Salmonella docking and invasion. This increased susceptibility at the cellular level manifests as increased susceptibility to typhoid fever in a Vietnamese population. Furthermore, treating zebrafish with a cholesterol-lowering agent, ezetimibe, reduced susceptibility to S Typhi. Thus, coupling multiple genetic association studies with mechanistic dissection revealed how VAC14 regulates Salmonella invasion and typhoid fever susceptibility and may open doors to new prophylactic/therapeutic approaches.

Bacterial risk factors for treatment failure and relapse among patients with isoniazid resistant tuberculosis.

BACKGROUND: Drug resistant tuberculosis (TB) is increasing in prevalence worldwide. Treatment failure and relapse is known to be high for patients with isoniazid resistant TB treated with standard first line regimens. However, risk factors for unfavourable outcomes and the optimal treatment regimen for isoniazid resistant TB are unknown. This cohort study was conducted when Vietnam used the eight month first line treatment regimen and examined risk factors for failure/relapse among patients with isoniazid resistant TB. METHODS: Between December 2008 and June 2011 2090 consecutive HIV-negative adults (≥18 years of age) with new smear positive pulmonary TB presenting at participating district TB units in Ho Chi Minh City were recruited. Participants with isoniazid resistant TB identified by Microscopic Observation Drug Susceptibility (MODS) had extended follow-up for 2 years with mycobacterial culture to test for relapse. MGIT drug susceptibility testing confirmed 239 participants with isoniazid resistant, rifampicin susceptible TB. Bacterial and demographic factors were analysed for association with treatment failure and relapse. RESULTS: Using only routine programmatic sputum smear microscopy for assessment, (months 2, 5 and 8) 30/239 (12.6%) had an unfavourable outcome by WHO criteria. Thirty-nine patients were additionally detected with unfavourable outcomes during 2 year follow up, giving a total of 69/239 (28.9%) of isoniazid (INH) resistant cases with unfavourable outcome by 2 years of follow-up. Beijing lineage was the only factor significantly associated with unfavourable outcome among INH-resistant TB cases during 2 years of follow-up. (adjusted OR = 3.16 [1.54-6.47], P = 0.002). CONCLUSION: One third of isoniazid resistant TB cases suffered failure/relapse within 2 years under the old eight month regimen. Over half of these cases were not identified by standard WHO recommended treatment monitoring. Intensified research on early identification and optimal regimens for isoniazid resistant TB is needed. Infection with Beijing genotype of TB is a significant risk factor for bacterial persistence on treatment resulting in failure/relapse within 2 years. The underlying mechanism of increased tolerance for standard drug regimens in Beijing genotype strains remains unknown.

Leukotriene A4 Hydrolase Genotype and HIV Infection Influence Intracerebral Inflammation and Survival From Tuberculous Meningitis.

Background: Tuberculous meningitis (TBM) is the most devastating form of tuberculosis, yet very little is known about the pathophysiology. We hypothesized that the genotype of leukotriene A4 hydrolase (encoded by LTA4H), which determines inflammatory eicosanoid expression, influences intracerebral inflammation, and predicts survival from TBM. Methods: We characterized the pretreatment clinical and intracerebral inflammatory phenotype and 9-month survival of 764 adults with TBM. All were genotyped for single-nucleotide polymorphism rs17525495, and inflammatory phenotype was defined by cerebrospinal fluid (CSF) leukocyte and cytokine concentrations. Results: LTA4H genotype predicted survival of human immunodeficiency virus (HIV)-uninfected patients, with TT-genotype patients significantly more likely to survive TBM than CC-genotype patients, according to Cox regression analysis (univariate P = .040 and multivariable P = .037). HIV-uninfected, TT-genotype patients had high CSF proinflammatory cytokine concentrations, with intermediate and lower concentrations in those with CT and CC genotypes. Increased CSF cytokine concentrations correlated with more-severe disease, but patients with low CSF leukocytes and cytokine concentrations were more likely to die from TBM. HIV infection independently predicted death due to TBM (hazard ratio, 3.94; 95% confidence interval, 2.79-5.56) and was associated with globally increased CSF cytokine concentrations, independent of LTA4H genotype. Conclusions: LTA4H genotype and HIV infection influence pretreatment inflammatory phenotype and survival from TBM. LTA4H genotype may predict adjunctive corticosteroid responsiveness in HIV-uninfected individuals.

High prevalence of PI resistance in patients failing second-line ART in Vietnam.

BACKGROUND: There are limited data from resource-limited settings on antiretroviral resistance mutations that develop in patients failing second-line PI ART. METHODS: We performed a cross-sectional virological assessment of adults on second-line ART for ≥6 months between November 2006 and December 2011, followed by a prospective follow-up over 2 years of patients with virological failure (VF) at the Hospital for Tropical Diseases, Vietnam. VF was defined as HIV RNA concentrations ≥1000 copies/mL. Resistance mutations were identified by population sequencing of the pol gene and interpreted using the 2014 IAS-USA mutation list and the Stanford algorithm. Logistic regression modelling was performed to identify predictors of VF. RESULTS: Two hundred and thirty-one patients were enrolled in the study. The median age was 32 years; 81.0% were male, 95.7% were on a lopinavir/ritonavir-containing regimen and 22 (9.5%) patients had VF. Of the patients with VF, 14 (64%) carried at least one major protease mutation [median: 2 (IQR: 1-3)]; 13 (59%) had multiple protease mutations conferring intermediate- to high-level resistance to lopinavir/ritonavir. Mutations conferring cross-resistance to etravirine, rilpivirine, tipranavir and darunavir were identified in 55%, 55%, 45% and 27% of patients, respectively. Higher viral load, adherence <95% and previous indinavir use were independent predictors of VF. The 2 year outcomes of the patients maintained on lopinavir/ritonavir included: death, 7 (35%); worsening virological/immunological control, 6 (30%); and virological re-suppression, 5 (25%). Two patients were switched to raltegravir and darunavir/ritonavir with good HIV control. CONCLUSIONS: High-prevalence PI resistance was associated with previous indinavir exposure. Darunavir plus an integrase inhibitor and lamivudine might be a promising third-line regimen in Vietnam.

Common polymorphisms in the CD43 gene region are associated with tuberculosis disease and mortality.

CD43, a surface glycoprotein, regulates Mycobacterium tuberculosis macrophage binding, replication, and proinflammatory cytokine induction in a murine model. We hypothesized that single-nucleotide polymorphisms (SNPs) in the CD43 gene region are associated with human tuberculosis (TB) susceptibility. We performed a case-population study in discovery (352 TB cases and 382 control subjects) and validation cohorts (339 TB cases and 376 control subjects). We examined whether 11 haplotype-tagging SNPs in the CD43 gene region were associated with tuberculous meningitis (TBM) and pulmonary TB (PTB) in Vietnam. Three SNPs from the CD43 gene region were associated with TB susceptibility with a genotypic model. The association fit a recessive genetic model and was greater for TBM than for PTB (for TBM: rs4788172, odds ratio [OR], 1.64; 95% confidence interval [CI], 1.04-2.59, rs17842268 [OR, 2.20; 95% CI, 1.29-3.76, and rs12596308 [OR, 2.38; 95% CI, 1.47-3.89]). Among TBM cases, rs17842268 was associated with decreased survival (hazard ratio, 2.7; 95% CI, 1.1-6.5; P = 0.011). In addition, rs12596308 and rs17842268 were associated with focal neurologic deficit at TBM presentation. Our data suggest that CD43 polymorphisms are associated with TB susceptibility, disease manifestations, and worse outcomes. To our knowledge, this is the first report that links CD43 genetic variants with susceptibility and outcome from a disease.

Human TOLLIP regulates TLR2 and TLR4 signaling and its polymorphisms are associated with susceptibility to tuberculosis.

Tuberculosis, one of the leading causes of death worldwide, stimulates inflammatory responses with beneficial and pathologic consequences. The regulation and nature of an optimal inflammatory response to Mycobacterium tuberculosis remains poorly understood in humans. Insight into mechanisms of negative regulation of the TLR-mediated innate immune response to M. tuberculosis could provide significant breakthroughs in the design of new vaccines and drugs. We hypothesized that TOLLIP and its common variants negatively regulate TLR signaling in human monocytes and are associated with susceptibility to tuberculosis. Using short hairpin RNA knockdown of TOLLIP in peripheral blood human monocytes, we found that TOLLIP suppresses TNF and IL-6 production after stimulation with TLR2 and TLR4 ligands. In contrast, secretion of the anti-inflammatory cytokine IL-10 was induced by TOLLIP. We also discovered two common polymorphisms that are associated with either decreased levels of mRNA expression (rs3750920) or increased IL-6 production (rs5743899) in a sample of 56 healthy volunteers. Furthermore, in a case-population study in Vietnam with 760 cord blood samples and 671 TB case patients, we found that SNPs rs3750920 and rs5743899 were associated with susceptibility to tuberculosis (p = 7.03 × 10(-16) and 6.97 × 10(-7), respectively). These data demonstrate that TOLLIP has an anti-inflammatory effect on TLR signaling in humans and that TOLLIP deficiency is associated with an increased risk of tuberculosis. To our knowledge, these data also show the first associations of TOLLIP polymorphisms with any infectious disease. These data also implicate an unexpected mechanism of negative regulation of TLR signaling in human tuberculosis pathogenesis.

Identification of bacterial determinants of tuberculosis infection and treatment outcomes: a phenogenomic analysis of clinical strains.

BACKGROUND: Bacterial diversity could contribute to the diversity of tuberculosis infection and treatment outcomes observed clinically, but the biological basis of this association is poorly understood. The aim of this study was to identify associations between phenogenomic variation in Mycobacterium tuberculosis and tuberculosis clinical features. METHODS: We developed a high-throughput platform to define phenotype-genotype relationships in M tuberculosis clinical isolates, which we tested on a set of 158 drug-sensitive M tuberculosis strains sampled from a large tuberculosis clinical study in Ho Chi Minh City, Viet Nam. We tagged the strains with unique genetic barcodes in multiplicate, allowing us to pool the strains for in-vitro competitive fitness assays across 16 host-relevant antibiotic and metabolic conditions. Relative fitness was quantified by deep sequencing, enumerating output barcode read counts relative to input normalised values. We performed a genome-wide association study to identify phylogenetically linked and monogenic mutations associated with the in-vitro fitness phenotypes. These genetic determinants were further associated with relevant clinical outcomes (cavitary disease and treatment failure) by calculating odds ratios (ORs) with binomial logistic regressions. We also assessed the population-level transmission of strains associated with cavitary disease and treatment failure using terminal branch length analysis of the phylogenetic data. FINDINGS: M tuberculosis clinical strains had diverse growth characteristics in host-like metabolic and drug conditions. These fitness phenotypes were highly heritable, and we identified monogenic and phylogenetically linked variants associated with the fitness phenotypes. These data enabled us to define two genetic features that were associated with clinical outcomes. First, mutations in Rv1339, a phosphodiesterase, which were associated with slow growth in glycerol, were further associated with treatment failure (OR 5·34, 95% CI 1·21-23·58, p=0·027). Second, we identified a phenotypically distinct slow-growing subclade of lineage 1 strains (L1.1.1.1) that was associated with cavitary disease (OR 2·49, 1·11-5·59, p=0·027) and treatment failure (OR 4·76, 1·53-14·78, p=0·0069), and which had shorter terminal branch lengths on the phylogenetic tree, suggesting increased transmission. INTERPRETATION: Slow growth under various antibiotic and metabolic conditions served as in-vitro intermediate phenotypes underlying the association between M tuberculosis monogenic and phylogenetically linked mutations and outcomes such as cavitary disease, treatment failure, and transmission potential. These data suggest that M tuberculosis growth regulation is an adaptive advantage for bacterial success in human populations, at least in some circumstances. These data further suggest markers for the underlying bacterial processes that contribute to these clinical outcomes. FUNDING: National Health and Medical Research Council/A∗STAR, National Institutes of Allergy and Infectious Diseases, National Institute of Child Health and Human Development, and the Wellcome Trust Fellowship in Public Health and Tropical Medicine.

Common polymorphisms in the PKP3-SIGIRR-TMEM16J gene region are associated with susceptibility to tuberculosis.

BACKGROUND: Tuberculosis has been associated with genetic variation in host immunity. We hypothesized that single-nucleotide polymorphisms (SNPs) in SIGIRR, a negative regulator of Toll-like receptor/IL-1R signaling, are associated with susceptibility to tuberculosis. METHODS: We used a case-population study design in Vietnam with cases that had either tuberculous meningitis or pulmonary tuberculosis. We genotyped 6 SNPs in the SIGIRR gene region (including the adjacent genes PKP3 and TMEM16J) in a discovery cohort of 352 patients with tuberculosis and 382 controls. Significant associations were genotyped in a validation cohort (339 patients with tuberculosis, 376 controls). RESULTS: Three SNPs (rs10902158, rs7105848, rs7111432) were associated with tuberculosis in discovery and validation cohorts. The polymorphisms were associated with both tuberculous meningitis and pulmonary tuberculosis and were strongest with a recessive genetic model (odds ratios, 1.5-1.6; P = .0006-.001). Coinheritance of these polymorphisms with previously identified risk alleles in Toll-like receptor 2 and TIRAP was associated with an additive risk of tuberculosis susceptibility. CONCLUSIONS: These results demonstrate a strong association of SNPs in the PKP3-SIGIRR-TMEM16J gene region and tuberculosis in discovery and validation cohorts. To our knowledge, these are the first associations of polymorphisms in this region with any disease.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history.

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Direct inference and control of genetic population structure from RNA sequencing data.

RNAseq data can be used to infer genetic variants, yet its use for estimating genetic population structure remains underexplored. Here, we construct a freely available computational tool (RGStraP) to estimate RNAseq-based genetic principal components (RG-PCs) and assess whether RG-PCs can be used to control for population structure in gene expression analyses. Using whole blood samples from understudied Nepalese populations and the Geuvadis study, we show that RG-PCs had comparable results to paired array-based genotypes, with high genotype concordance and high correlations of genetic principal components, capturing subpopulations within the dataset. In differential gene expression analysis, we found that inclusion of RG-PCs as covariates reduced test statistic inflation. Our paper demonstrates that genetic population structure can be directly inferred and controlled for using RNAseq data, thus facilitating improved retrospective and future analyses of transcriptomic data.

Disruption of Aldehyde Dehydrogenase 2 protects against bacterial infection.

The ALDH2*2 (rs671) allele is one of the most common genetic mutations in humans, yet the positive evolutionary selective pressure to maintain this mutation is unknown, despite its association with adverse health outcomes. ALDH2 is responsible for the detoxification of metabolically produced aldehydes, including lipid-peroxidation end products derived from inflammation. Here, we demonstrate that host-derived aldehydes 4-hydroxynonenal (4HNE), malondialdehyde (MDA), and formaldehyde (FA), all of which are metabolized by ALDH2, are directly toxic to the bacterial pathogens Mycobacterium tuberculosis and Francisella tularensis at physiological levels. We find that Aldh2 expression in macrophages is decreased upon immune stimulation, and that bone marrow-derived macrophages from Aldh2 -/- mice contain elevated aldehydes relative to wild-type mice. Macrophages deficient for Aldh2 exhibited enhanced control of Francisella infection. Finally , mice lacking Aldh2 demonstrated increased resistance to pulmonary infection by M. tuberculosis , including in a hypersusceptible model of tuberculosis, and were also resistant to Francisella infection. We hypothesize that the absence of ALDH2 contributes to the host's ability to control infection by pathogens such as M. tuberculosis and F. tularensis , and that host-derived aldehydes act as antimicrobial factors during intracellular bacterial infections. One sentence summary: Aldehydes produced by host cells contribute to the control of bacterial infections.

High-throughput phenogenotyping of Mycobacteria tuberculosis clinical strains reveals bacterial determinants of treatment outcomes.

Background: Combatting the tuberculosis (TB) epidemic caused by Mycobacterium tuberculosis ( Mtb ) necessitates a better understanding of the factors contributing to patient clinical outcomes and transmission. While host and environmental factors have been evaluated, the impact of Mtb genetic background and phenotypic diversity is underexplored. Previous work has made associations between Mtb genetic lineages and some clinical and epidemiological features, but the bacterial traits underlying these connections are largely unknown. Methods: We developed a high-throughput functional genomics platform for defining genotype-phenotype relationships across a panel of Mtb clinical isolates. These phenotypic fitness profiles function as intermediate traits which can be linked to Mtb genetic variants and associated with clinical and epidemiological outcomes. We applied this approach to a collection of 158 Mtb strains from a study of Mtb transmission in Ho Chi Minh City, Vietnam. Mtb strains were genetically tagged in multiplicate, which allowed us to pool the strains and assess in vitro competitive fitness using deep sequencing across a set of 14 host-relevant antibiotic and metabolic conditions. Phylogenetic and monogenic associations with these intermediate traits were identified and then associated with clinical outcomes. Findings: Mtb clinical strains have a broad range of growth and drug response dynamics that can be clustered by their phylogenetic relationships. We identified novel monogenic associations with Mtb fitness in various metabolic and antibiotic conditions. Among these, we find that mutations in Rv1339 , a phosphodiesterase, which were identified through their association with slow growth in glycerol, are further associated with treatment failure. We also identify a previously uncharacterized subclade of Lineage 1 strains (L1.1.1.1) that is phenotypically distinguished by slow growth under most antibiotic and metabolic stress conditions in vitro . This clade is associated with cavitary disease, treatment failure, and demonstrates increased transmission potential. Interpretation: High-throughput phenogenotyping of Mtb clinical strains enabled bacterial intermediate trait identification that can provide a mechanistic link between Mtb genetic variation and patient clinical outcomes. Mtb strains associated with cavitary disease, treatment failure, and transmission potential display intermediate phenotypes distinguished by slow growth under various antibiotic and metabolic conditions. These data suggest that Mtb growth regulation is an adaptive advantage for host bacterial success in human populations, in at least some circumstances. These data further suggest markers for the underlying bacterial processes that govern these clinical outcomes. Funding: National Institutes of Allergy and Infectious Diseases: P01 AI132130 (SS, SMF); P01 AI143575 (XW, SMF); U19 AI142793 (QL, SMF); 5T32AI132120-03 (SS); 5T32AI132120-04 (SS); 5T32AI049928-17 (SS) Wellcome Trust Fellowship in Public Health and Tropical Medicine: 097124/Z/11/Z (NTTT) National Health and Medical Research Council (NHMRC)/A*STAR joint call: APP1056689 (SJD) The funding sources had no involvement in study methodology, data collection, analysis, and interpretation nor in the writing or submission of the manuscript. Research in context: Evidence before this study: We used different combinations of the words mycobacterium tuberculosis, tuberculosis, clinical strains, intermediate phenotypes, genetic barcoding, phenogenomics, cavitary disease, treatment failure, and transmission to search the PubMed database for all studies published up until January 20 th , 2022. We only considered English language publications, which biases our search. Previous work linking Mtb determinants to clinical or epidemiological data has made associations between bacterial lineage, or less frequently, genetic polymorphisms to in vitro or in vivo models of pathogenesis, transmission, and clinical outcomes such as cavitary disease, treatment failure, delayed culture conversion, and severity. Many of these studies focus on the global pandemic Lineage 2 and Lineage 4 Mtb strains due in part to a deletion in a polyketide synthase implicated in host-pathogen interactions. There are a number of Mtb GWAS studies that have led to novel genetic determinants of in vitro drug resistance and tolerance. Previous Mtb GWAS analyses with clinical outcomes did not experimentally test any predicted phenotypes of the clinical strains. Published laboratory-based studies of Mtb clinical strains involve relatively small numbers of strains, do not identify the genetic basis of relevant phenotypes, or link findings to the corresponding clinical outcomes. There are two recent studies of other pathogens that describe phenogenomic analyses. One study of 331 M. abscessus clinical strains performed one-by-one phenotyping to identify bacterial features associated with clearance of infection and another details a competition experiment utilizing three barcoded Plasmodium falciparum clinical isolates to assay antimalarial fitness and resistance. Added value of this study: We developed a functional genomics platform to perform high-throughput phenotyping of Mtb clinical strains. We then used these phenotypes as intermediate traits to identify novel bacterial genetic features associated with clinical outcomes. We leveraged this platform with a sample of 158 Mtb clinical strains from a cross sectional study of Mtb transmission in Ho Chi Minh City, Vietnam. To enable high-throughput phenotyping of large numbers of Mtb clinical isolates, we applied a DNA barcoding approach that has not been previously utilized for the high-throughput analysis of Mtb clinical strains. This approach allowed us to perform pooled competitive fitness assays, tracking strain fitness using deep sequencing. We measured the replicative fitness of the clinical strains in multiplicate under 14 metabolic and antibiotic stress condition. To our knowledge, this is the largest phenotypic screen of Mtb clinical isolates to date. We performed bacterial GWAS to delineate the Mtb genetic variants associated with each fitness phenotype, identifying monogenic associations with several conditions. We then defined Mtb phenotypic and genetic features associated with clinical outcomes. We find that a subclade of Mtb strains, defined by variants largely involved in fatty acid metabolic pathways, share a universal slow growth phenotype that is associated with cavitary disease, treatment failure and increased transmission potential in Vietnam. We also find that mutations in Rv1339 , a poorly characterized phosphodiesterase, also associate with slow growth in vitro and with treatment failure in patients. Implications of all the available evidence: Phenogenomic profiling demonstrates that Mtb strains exhibit distinct growth characteristics under metabolic and antibiotic stress conditions. These fitness profiles can serve as intermediate traits for GWAS and association with clinical outcomes. Intermediate phenotyping allows us to examine potential processes by which bacterial strain differences contribute to clinical outcomes. Our study identifies clinical strains with slow growth phenotypes under in vitro models of antibiotic and host-like metabolic conditions that are associated with adverse clinical outcomes. It is possible that the bacterial intermediate phenotypes we identified are directly related to the mechanisms of these outcomes, or they may serve as markers for the causal yet unidentified bacterial determinants. Via the intermediate phenotyping, we also discovered a surprising diversity in Mtb responses to the new anti-mycobacterial drugs that target central metabolic processes, which will be important in considering roll-out of these new agents. Our study and others that have identified Mtb determinants of TB clinical and epidemiological phenotypes should inform efforts to improve diagnostics and drug regimen design.