History of tuberculosis disease is associated with genetic regulatory variation in Peruvians.

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

Single-cell eQTL models reveal dynamic T cell state dependence of disease loci.

Non-coding genetic variants may cause disease by modulating gene expression. However, identifying these expression quantitative trait loci (eQTLs) is complicated by differences in gene regulation across fluid functional cell states within cell types. These states-for example, neurotransmitter-driven programs in astrocytes or perivascular fibroblast differentiation-are obscured in eQTL studies that aggregate cells1,2. Here we modelled eQTLs at single-cell resolution in one complex cell type: memory T cells. Using more than 500,000 unstimulated memory T cells from 259 Peruvian individuals, we show that around one-third of 6,511 cis-eQTLs had effects that were mediated by continuous multimodally defined cell states, such as cytotoxicity and regulatory capacity. In some loci, independent eQTL variants had opposing cell-state relationships. Autoimmune variants were enriched in cell-state-dependent eQTLs, including risk variants for rheumatoid arthritis near ORMDL3 and CTLA4; this indicates that cell-state context is crucial to understanding potential eQTL pathogenicity. Moreover, continuous cell states explained more variation in eQTLs than did conventional discrete categories, such as CD4+ versus CD8+, suggesting that modelling eQTLs and cell states at single-cell resolution can expand insight into gene regulation in functionally heterogeneous cell types.

Multimodally profiling memory T cells from a tuberculosis cohort identifies cell state associations with demographics, environment and disease.

Multimodal T cell profiling can enable more precise characterization of elusive cell states underlying disease. Here, we integrated single-cell RNA and surface protein data from 500,089 memory T cells to define 31 cell states from 259 individuals in a Peruvian tuberculosis (TB) progression cohort. At immune steady state >4 years after infection and disease resolution, we found that, after accounting for significant effects of age, sex, season and genetic ancestry on T cell composition, a polyfunctional type 17 helper T (TH17) cell-like effector state was reduced in abundance and function in individuals who previously progressed from Mycobacterium tuberculosis (M.tb) infection to active TB disease. These cells are capable of responding to M.tb peptides. Deconvoluting this state-uniquely identifiable with multimodal analysis-from public data demonstrated that its depletion may precede and persist beyond active disease. Our study demonstrates the power of integrative multimodal single-cell profiling to define cell states relevant to disease and other traits.

CD1b Tetramers Broadly Detect T Cells That Correlate With Mycobacterial Exposure but Not Tuberculosis Disease State.

The non-polymorphic nature of CD1 proteins creates a situation in which T cells with invariant T cell receptors (TCRs), like CD1d-specific NKT cells, are present in all humans. CD1b is an abundant protein on human dendritic cells that presents M. tuberculosis (Mtb) lipid antigens to T cells. Analysis of T cell clones suggested that semi-invariant TCRs exist in the CD1b system, but their prevalence in humans is not known. Here we used CD1b tetramers loaded with mycolic acid or glucose monomycolate to study polyclonal T cells from 150 Peruvian subjects. We found that CD1b tetramers loaded with mycolic acid or glucose monomycolate antigens stained TRAV1-2+ GEM T cells or TRBV4-1+ LDN5-like T cells in the majority of subjects tested, at rates ~10-fold lower than NKT cells. Thus, GEM T cells and LDN5-like T cells are a normal part of the human immune system. Unlike prior studies measuring MHC- or CD1b-mediated activation, this large-scale tetramer study found no significant differences in rates of CD1b tetramer-mycobacterial lipid staining of T cells among subjects with Mtb exposure, latent Mtb infection or active tuberculosis (TB) disease. In all subjects, including "uninfected" subjects, CD1b tetramer+ T cells expressed memory markers at high levels. However, among controls with lower mycobacterial antigen exposure in Boston, we found significantly lower frequencies of T cells staining with CD1b tetramers loaded with mycobacterial lipids. These data link CD1b-specific T cell detection to mycobacterial exposure, but not TB disease status, which potentially explains differences in outcomes among CD1-based clinical studies, which used control subjects with low Mtb exposure.

Early progression to active tuberculosis is a highly heritable trait driven by 3q23 in Peruvians.

Of the 1.8 billion people worldwide infected with Mycobacterium tuberculosis, 5-15% will develop active tuberculosis (TB). Approximately half will progress to active TB within the first 18 months after infection, presumably because they fail to mount an effective initial immune response. Here, in a genome-wide genetic study of early TB progression, we genotype 4002 active TB cases and their household contacts in Peru. We quantify genetic heritability ([Formula: see text]) of early TB progression to be 21.2% (standard error 0.08). This suggests TB progression has a strong genetic basis, and is comparable to traits with well-established genetic bases. We identify a novel association between early TB progression and variants located in a putative enhancer region on chromosome 3q23 (rs73226617, OR = 1.18; P = 3.93 × 10-8). With in silico and in vitro analyses we identify rs73226617 or rs148722713 as the likely functional variant and ATP1B3 as a potential causal target gene with monocyte specific function.