Expression of specific HLA class II alleles is associated with an increased risk for active tuberculosis and a distinct gene expression profile =============================================================================================================================================== * Leila Y. Chihab * Rebecca Kuan * Elizabeth J. Phillips * Simon A. Mallal * Virginie Rozot * Mark M. Davis * Thomas J. Scriba * Alessandro Sette * Bjoern Peters * Cecilia S. Lindestam Arlehamn * the SATVI Study Group ## Abstract Several HLA allelic variants have been associated with protection from or susceptibility to infectious and autoimmune diseases. Here, we examined whether specific HLA alleles would be associated with different *Mycobacterium tuberculosis (Mtb)* infection outcomes. The HLA alleles present at the -A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, and -DRB3/4/5 loci were determined in a cohort of 636 individuals with known *Mtb* infection outcomes from South Africa and the US. Among these individuals, 233 were QuantiFERON (QFT) negative, and 433 were QFT positive, indicating *Mtb* exposure. Of these, 99 QFT positive participants either had active tuberculosis (TB) upon enrollment or were diagnosed in the past. We found that DQA1*03:01, DPB1*04:02, and DRB4*01:01 were significantly more frequent in individuals with active TB (susceptibility alleles), as judged by Odds Ratios and associated p-values, while DPB1*105:01 was associated with protection from active TB. Peripheral blood mononuclear cells (PMBCs) from a subset of individuals were stimulated with *Mtb* antigens, revealing individuals who express any of the three susceptibility alleles were associated with lower magnitude of responses. Furthermore, we defined a gene signature associated with individuals expressing the susceptibility alleles that was characterized by lower expression of APC-related genes. In summary, we have identified specific HLA alleles associated with susceptibility to active TB and found that the expression of these alleles was associated with a decreased *Mtb*-specific T cell response and a specific gene expression signature. These results will help understand individual risk factors in progressing to active TB. ## Introduction *Mycobacterium tuberculosis* (*Mtb*), the cause of tuberculosis (TB) disease, infects about one-quarter of the world’s population and causes approximately 1.5 million deaths annually 1. Although *Mtb* infection is fairly prevalent worldwide, less than 10% of infections lead to active tuberculosis (ATB) disease1,2. In some cases individuals with *Mtb* infection can completely clear all living bacteria, while in most *Mtb* infection cases the bacteria persist and these individuals have an increased risk of bacterial reactivation leading to ATB 3. Latent infection can last weeks, months, or years with the likelihood of progression to ATB being highest within the first two years following infection. Certain conditions can increase the likelihood of progression to ATB disease (e.g. HIV co-infection, immune suppression therapy, malnutrition, and diabetes) 4. To reduce the incidence of TB the WHO has identified a strategy which includes identification and treatment of individuals with ATB, screening of individuals at high risk, as well as preventative therapy to those at high risk of progression to ATB 5. Importantly, the biological correlates that distinguish between different *Mtb* infection outcomes have yet to be fully elucidated and thus there is a need for improved prognostic and diagnostic tests to identify those at risk for ATB. Human Leukocyte Antigens (HLAs) are heterodimeric proteins that play a critical role in the immune recognition of antigens. HLA molecules present pathogen-derived peptides (epitopes) for recognition by T cells 6. HLA molecules are grouped into class I vs. class II, each encompassing several highly polymorphic loci with more than 15,000 different possible HLA allelic variants 7. In general, different HLA allelic variants present distinct epitope repertoires, and this high level of diversity results in a highly diverse set of peptide epitopes being presented in different individuals in a given population. In this context, genetic analyses and genome wide association studies (GWAS) have demonstrated that the presence of specific HLA allelic variants can exert either a protective or susceptibility effect for different autoimmune or infectious diseases 8-14. Here, we examined whether any specific HLA alleles are associated with ATB and thus could be used to identify individuals at risk for ATB, and whether expression of these HLA alleles can be linked to differences in the *Mtb*-specific immune response. ## Materials and Methods ### Ethics Statement This research was performed in accordance with approval from the Human Research Ethics Committees of the University of Cape Town (Cape Town, South Africa) and the La Jolla Institute for Immunology (La Jolla, USA). All participants provided written informed consent before participation. Adolescents provided written informed assent, with a parent or legal guardian providing written informed consent. ### Study Participants Study participants (n=620) from the Worcester region of South Africa were investigated. Participants were divided into a healthy cohort of QFT positive (n=318) and negative (203) individuals (QuantiFEROn-TB Gold In-Tube, Cellestis) and into individuals with ongoing (sputum XpertMTB/RIF-positive (Cepheid, Sunnyvale, CA)) TB disease (n=28) or a self-reported history of individuals with previously diagnosed active TB who underwent TB treatment (n=71). We also enrolled QFT+ participants (n=16) at the University of California, San Diego Anti-Viral Research Center (San Diego, USA). QFT status and whether individuals have or have had active TB is listed in **Supplementary table 1**. ### PBMC Isolation PBMC were purified from whole blood either using CPT tubes (BD), for the adolescents, or layered onto Ficoll (for adults) using 50ml Leukosep tubes (Greiner) by density-gradient centrifugation on Ficoll, according to the manufacturer’s instructions. Cells were suspended in FBS containing 10% (vol/vol) DMSO and cryopreserved in liquid nitrogen. The cryopreserved cells from South Africa were shipped to LJI for analysis. Cryopreserved PBMC were quickly thawed by incubating each cryovial at 37°C for 2 min, and cells transferred to cold medium (RPMI 1640 with L-glutamin and 25mM HEPES; Omega Scientific), supplemented with 5% human AB serum (GemCell), 1% penicillin streptomycin (Life Technologies), 1% glutamax (Life Technologies) and 20 U/ml benzonase nuclease (MilliporeSigma). Cells were centrifuged and resuspended in medium to determine cell concentration and viability using trypan blue and a hematocytometer. ### Peptides and other stimuli Peptides were synthesized as crude material on a small (1 mg) scale by A&A LLC (San Diego). Peptides were pooled into peptide pools, re-lyophilized and reconstituted at a concentration of 0.7mg/ml for MTB300 (PMID27409590). *Mtb* whole cell lysate (WCL) from strain H37Rv was obtained from BEI Resources, NIAID, NIH (NR-14822). ### HLA typing Participants were HLA-typed at the La Jolla Institute, by the American Society for Histocompatibility and Immunogenetics (ASHI)-accredited laboratory at Murdoch University (Western Australia), or Immucor, Sirona Genomics, CA. Typing at LJI was performed by next-generation sequencing 15. Specifically, amplicons were generated from the appropriate class II locus for exons 2 through 4 by PCR amplification. From these amplicons, sequencing libraries were generated (Illumina Nextera XT) and sequenced with MiSeq Reagent Kit v3 as per the manufacturer’s instructions (Illumina). Sequence reads were matched to HLA alleles and participant genotypes were assigned. Typing in Australia for class I (HLA A, B, C) and class II (DQA1, DQB1, DRB1, 3, 4, 5, DPB1) was performed using locus-specific PCR amplification of genomic DNA. Patient-specific, barcoded primers were used for amplification. Amplified products were quantitated and pooled by participant and up to 48 participants were pooled. An indexed (8 indexed MiSeq runs) library was then quantitated using Kappa universal QPCR library quantification kits. Sequencing was performed using an Illumina MiSeq using 2×300 paired-end chemistry. Reads were quality-filtered and passed through a proprietary allele calling algorithm and analysis pipeline using the latest IMGT HLA allele database as a reference. The algorithm was developed by E.J.P. and S.A.M. and relies on periodically updated versions of the freely available international immunogenetics information system ([http://www.imgt.org](http://www.imgt.org)) and an ASHI-accredited HLA allele caller software pipeline, IIID HLA analysis suite ([http://www.iiid.com.au/laboratory-testing/](http://www.iiid.com.au/laboratory-testing/)). Typing at Immucor, Sirona Genomics was performed using the MIA FORA NGS kit and analysis software (Immucor, Inc.) as previously described 16. The HLA type of each subject is listed in **Supplementary table 1**. Population coverage of HLA alleles were calculated using the population coverage tool from the Immune epitope database ([http://www.iedb.org](http://www.iedb.org)). ### ELISPOT assay PBMCs were thawed and antigen-specific cellular responses were measured by IFNγ ELISPOT assay with all antibodies and reagents from Mabtech (Nacka Strand, Sweden). Plates (96-well Immobilion-P; Millipore) were coated overnight at 4°C with 5µg/ml anti-IFNI (1-D1K; Mabtech). Briefly, 200,000 cells were plated in each well of the pre-coated plate and were stimulated in triplicate with 10µg/ml *Mtb* whole cell lysate, 1 µg/ml MTB300, 10µg/ml PHA as a positive control, or DMSO corresponding to the concentration present in the peptide pool. After 20h incubation at 37°C, wells were washed with PBS/0.05% Tween 20 and incubated with biotinylated anti-IFN□ (7-B6-1; Mabtech) for 2h. Spots were developed using Vectastain APC peroxidase (Vector Laboratories) and 3-amino-9-ethylcarbazole (Sigma-Aldrich). Spots were counted by computer-assisted image analysis (AID iSpot; Aid Diagnostica). Responses were considered positive if the net spot-forming cells (SFC) per 106 PBMC were ≥20, the stimulation index ≥2, and p≤0.05 (Student’s t-test, mean of triplicate values of the response against relevant stimuli vs. the DMSO control). ### Cell sorting and RNA purification PBMCs were thawed and added at a density of 1×106 cells per well to a round-bottom 96-well plate. They were stimulated for 24 hrs at 37°C with MTB300, *Mtb* whole cell lysate, anti-CD28 (1µg/ml, CD28.2; eBioscience) together with anti-CD3 (1µg/ml pre-coated overnight at 4°C, UCHT1; BioLegend), as well as left unstimulated as a negative control. After 24 hrs, cells were resuspended in PBS with 10% (v/v) FBS and incubated at 4°C for 10 min. Cells were then stained with fixable viability dye eFluor506 (eBioscience) and an antibody mixture containing anti-human CD3-AF700 (UCHT1; ThermoFisher), CD4-APCeFluor780 (RPA-T4; eBioscience), CD8-BV650 (RPA-T8; BioLegend), CCR7-PerCPCy5.5 (G043H7, Biolegend) CD14-V500 (M5E2; BD Bioscience), CD19-V500 (HIB19; BD Bioscience), CD25-FITC (M-A251; BD Bioscience), CD45RA-eFluor450 (HI100, eBiosciences), CD137-APC (4B4-1; BioLegend), OX40-PECy7 (Ber-ACT35; BioLegend), and PD-L1-PE (29E.2A3; BioLegend) for 20 min at room temperature. Cells were washed and resuspended in PBS before being transferred into a 5 ml polypropylene FACS tube (BD Bioscience). PBMCs were sorted, based on forward and side scatter excluding debris and doublets, on a FACSAria into TRIzol LS (Thermo Fisher). Acquisition files were analyzed using FlowJo. Total RNA was extracted from ∼100,000 cells in TRIzol LS using the miRNeasy Micro Kit (Qiagen) on a QIAcube (Qiagen). Total RNA was amplified according to Smart Seq protocol (2). cDNA was purified using AMPure XP beads. cDNA was used to prepare a standard barcoded sequencing library (Illumina). Samples were sequenced using an Illumina HiSeq2500 to obtain 50-bp single end reads. Samples that failed to be sequenced due to limited sample availability or failed the quality control were eliminated from further sequencing and analysis. The full protocol can be found at protocols.io ([http://dx.doi.org/10.17504/protocols.io.bxr6pm9e](http://dx.doi.org/10.17504/protocols.io.bxr6pm9e)). ### RNA-seq analysis Paired-end reads that passed Illumina filters were filtered for reads aligning to tRNA, rRNA, adapter sequences, and spike-in controls. The reads were aligned to the GRCh38 reference genome and Gencode v27 annotations using STAR (v2.6.1) 17. DUST scores were calculated with PRINSEQ Lite v0.20.3 18 and low-complexity reads (DUST > 4) were removed from BAM files. The alignment results were parsed via SAMtools 19 to generate SAM files. Read counts to each genomic feature were obtained with the featureCounts (v 1.6.5) 20 using the default option along with a minimum quality cut off (Phred > 10). After removing absent features (zero counts in all samples), the raw counts were imported into R v3.6.1 and genes with an average TPM < 1 were removed. R/Bioconductor package DESeq2 v.1.24.0 21 was used to normalize raw counts. Variance stabilizing transformation was applied to normalized counts to obtain log2 gene expression values. Quality control was performed using boxplots and Principal component analysis (PCA), using the ‘prcomp’ function in R, on log2 expression values. Differentially expressed genes were identified using the DESeq2 Wald test, and p-values were adjusted for multiple test correction using the Benjamini Hochberg algorithm 22. Genes with adjusted p values < 0.05 and log2 fold change >1 or < -1 were considered differentially expressed. Cell type enrichment was performed using DICE 23. The RNAseq data have been submitted to the Gene Expression Omnibus under accession number GSE207214 ([http://www.ncbi.nlm.nih.gov/geo/](http://www.ncbi.nlm.nih.gov/geo/)). ## Results ### HLA DQA1*03:01, DRB4*01:01, and DPB1*04:02 are associated with susceptibility for active TB and DPB1*105:01 with protection To determine whether specific HLA alleles are associated with active TB, samples were obtained from a cohort of 636 participants, of which 203 individuals were QFT negative, the remaining 433 individuals were QFT positive including 99 individuals with a history of active TB (ATB) or who were diagnosed with active TB upon enrollment. Participants were HLA typed by determining the specific alleles expressed at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1,-DRB1, and -DRB3/4/5 loci. The results of the HLA typing were utilized to investigate whether the expression of particular HLA alleles would be associated with *Mtb* infection outcomes. To identify allele associations with TB, odds ratios and p-values between active TB and QFT negative/QFT positive individuals were calculated utilizing the RATE algorithm 24. A type 1 error rate of 0.05 was accounted for via a Bonferroni correction on alleles expressed in at least 5% of the total study cohort (**Table 1**). View this table: [Table 1.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/T1) Table 1. HLA allele association with active TB Of the 324 HLA alleles expressed in one or more individuals of our cohort, we focused on the 87 alleles (37 class I and 50 class II alleles) that were present in at least 5% (i.e. 32 individuals) of the cohort to ensure we had sufficient power to detect differences in the HLA distribution. For these 87 alleles, 12 had a significant p-value when comparing their frequency in individuals with ATB compared to QFT negative and positive individuals (**Table 1**). Most of these variants are commonly expressed in the general population (Table 1). Of these variants, HLA DQA1*03:01 had the strongest association with ATB, with an odds ratio of 7.9 and p-value of 8.0E-10. A total of 45 participants carried this allele type, 24 of which (53%) had ATB. In addition, the DPB1*04:02 and DRB4*01:01 alleles were also associated with ATB, with odds ratios of 4.2 and 2.3 respectively. Finally, HLA DPB1*105:01 appeared to confer protection against ATB, since it was inversely correlated with ATB (OR= 0.4, p=0.0006). After Bonferroni correction, both the positive association between alleles DQA1*03:01, DPB1*04:02, DRB4*01:01 and ATB (p= 6.9E-8, 1.7E-4, 0.03), and the negative association between HLA DPB1*105:01 and ATB (p=0.05) were maintained. ### DQA1*04:01 and DRB1*03:02 are associated with protection from *Mtb* infection To protect from active disease, a factor like HLA could either protect from infection, or from progression to active TB. To investigate that we examined the alleles that were associated with susceptibility to active disease, and tested if these were enriched within the QFT positive. We found no alleles that were more frequent in QFT positive individuals compared to QFT negative (not shown). We next investigated the differences in HLA expression between *Mtb* infected (both QFT positive and ATB) participants and QFT negative participants, using the same approach as above (**Table 2**). We identified 19 alleles present in at least 5% of the participants with a significant p-value. In this comparison, DQA1*03:01 was again identified to have the closest association to *Mtb* infection with an odds ratio of 11.1 and a corrected p-value of 0.001. In addition, DQA1*04:01 and DRB1*03:02 were associated with protection from *Mtb* infection, both with odds ratios of 0.4 (p=0.033 and 0.050, respectively). In conclusion, we identified several alleles associated with ATB and *Mtb* infection susceptibility. We focused the rest of the study on the alleles associated with susceptibility for ATB. View this table: [Table 2.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/T2) Table 2. HLA allele association with *Mtb* infection ### DQB1*02:02, *03:03, *04:02, *06:01/02/09 co-expression with DQA1*03:01 increases susceptibility to ATB Mature HLA DQ molecules expressed on the cell surface are a dimer of DQ alpha and beta chains, and particular DQA alleles are preferentially co-expressed with particular DQB allelic variants. For this reason, we next sought to determine which DQ beta chain variants are most likely to be co-expressed with HLA DQA1*03:01 in ATB cases. To do so, we calculated the odds ratios and p-values of HLA DQ alpha and beta pairs to identify which pair is more closely associated to ATB (**Table 3**). Using the most prevalent DQ beta chains in our cohort, we paired various beta chains with DQA1*03:01 and determined the odds ratios with a confidence interval of 95%. When paired with beta chains DQB1*02:02, *04:02, *03:03, and *06:01/02, the odds ratios significantly increased from 7.9 to 54.2 (DQB1*02:02), 35.0 (DQB1*04:02), 8.4 (DQB1*03:03), 9.5 (DQB1*06:01), and infinity (DQB1*06:02; e.g. no expression in QFT positive or negative individuals), respectively. The odds ratio of association when DQA1*03:01 is paired to all possible DQB1 chains combined is 16, an increase from the original 7.9. Thus, when DQA1*03:01 is co-expressed with the beta chains highlighted in table 2, the association of DQA1*03:01 with ATB becomes even more pronounced. View this table: [Table 3.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/T3) Table 3. DQB1 alleles expressed with DQA1*03:01 ### Expression of susceptibility alleles is associated with lower magnitude of *Mtb*-specific response in an IFNγ ELISPOT assay We next investigated if the expression of the HLA alleles associated with ATB (**Table 1**) was associated with differences in the *Mtb*-specific immune response. PBMC samples from a subset of our cohort (all of which were QFT positive and 59 of which were ATB cases) were tested for IFNγ production in response to an *Mtb*-derived peptide pool (MTB300). For all three of the susceptibility alleles, (DQA1*03:01, DPB1*04:02, or DRB4*01:01) we found a trend towards a lower magnitude of IFNγ-specific responses, with the strongest trend being for DQA1*03:01 (**Figure 1a-c**). When the three alleles were considered together, we found a significantly lower magnitude of IFNγ-specific response in individuals with the susceptibility alleles compared to those without (**Figure 1d**). No significant difference was detected when comparing the responses from the ATB cohort with the QFT positive cohort (not shown). In conclusion, the expression of the three susceptibility alleles influences the *Mtb*-specific IFNγ response. ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/07/10/2022.07.06.22277293/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/F1) Figure 1. The expression of susceptibility alleles is associated with a lower magnitude of *Mtb*-specific IFNγ response. Magnitude of response against MTB300 as SFC per 106 PBMC measured in an IFNγ ELISPOT assay. Each point represents one individual (total n=144), median ± interquartile range is shown. Two-tailed Mann-Whitney U test. Individuals with susceptibility alleles are represented by closed circles and individuals without by open circles. a) DQA1*03:01 n=28 with and n=116 without, b) DPB1*04:02 n=36 with, n=108 without, c) DRB4*01:01 n=46 with, n=98 without, d) any susceptibility allele n=73 with, n=71 without. Red color indicates participants included in subsequent gene expression experiments. ### The presence of ATB susceptibility alleles is correlated with specific gene signatures To further investigate the biological mechanisms of HLA association with ATB susceptibility, PBMCs from a subset of QFT positive healthy individuals (**Fig 1d**, red color) were stimulated with MTB300 or a whole cell lysate (WCL) of the H37Rv TB strain. Live, singlet cells were sorted and their global gene expression was determined by RNA-Seq. The differential gene expression analysis of individuals who express the susceptibility alleles vs. individuals who lack the expression of these alleles in each stimulation condition from the ELISPOT assay (unstimulated, MTB300, and WC), revealed a specific gene signature. A total of 84 differentially expressed genes (DEGs) were identified across the three different stimulation conditions (**Supplementary table 2**; adjusted p values <0.05 and log2 fold change >1 or <-1). For the unstimulated samples, 5 genes were upregulated in individuals with the susceptibility alleles, and 29 genes were upregulated in those without. For the MTB300 stimulation condition, 8 genes were upregulated in with and 10 w/o, and for WCL it was 21 vs. 19 genes (**Figure 2a**). The 84 genes were filtered based on concordant expression across the different stimuli, e.g. if a gene was upregulated in individuals with the susceptibility allele in one condition but up in individuals without in another it was removed. This resulted in a gene signature of 63 genes (**Figure 2b and Supplementary table 2**). These genes could be grouped into 4 clusters based on their expression in individuals with vs. without the susceptibility alleles and across stimuli conditions (**Figure 2b**). Cluster 1 encompasses genes that were upregulated in the unstimulated and MTB300 condition, but with no difference between allele expression. The genes corresponding to Cluster 2 were upregulated in individuals with expression of the susceptibility alleles across the three different conditions compared to individuals lacking the alleles. Cluster 3 consisted of genes upregulated in individuals without the susceptibility alleles across different stimulations. Interestingly, cluster 3 included HLA-DQA1 and HLA-DQB1, as well as JUN and TIRAP. Cluster 4 was composed of genes with different patterns of upregulation in response to WCL, but no difference between allele expression. ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/07/10/2022.07.06.22277293/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/F2) Figure 2. Expression of the susceptibility alleles is associated with a specific gene signature. a) Volcano plots showing differentially expressed genes comparing individuals with susceptibility alleles (n=21) to those without (n=30) in three different stimuli conditions; unstimulated, MTB300, and WCL. Black indicates DEGs (adjusted p value <0.05 and log2 fold change >1 and <-1). b) Heatmap displaying the 63 DEGs from (a). Genes are ordered by hierarchical clustering, and samples are ordered according to allele expression and then according to stimulation condition. Black vertical lines indicate allele and stimuli separation and horizontal lines indicate cluster separation. c) Cell-type enrichment ([dice-database.org](http://dice-database.org)) for the 63 DEGs. Each color represents one gene. d) Gene expression at the mRNA level. Gene expression values in counts normalized by sequencing depth calculated by the DEseq2 package. Two-tailed Mann-Whitney U test. The 63 genes associated with the signature were enriched in genes expressed in B cells and monocytes (both classical and non-classical; **Figure 2c**). Finally, to visualize the gene expression changes between cohorts in each stimulus, 6 representative genes were chosen, 2 from cluster 2 (XPOT and CAB39) and 4 from cluster 3 (JUN, TIPARP, HLA-DQA1, HLA-DQB1; **Figure 2d**). The 4 genes from cluster 3 exhibited a trend for lower expression in individuals with the susceptibility alleles in each of the stimuli conditions. Conversely, the 2 genes from cluster 2 had higher expression in each stimuli condition in individuals with the susceptibility alleles. ### Gene signatures associated with ATB susceptibility alleles and low T cell responsiveness to *Mtb* antigens Finally, we investigated whether DEGs could be identified in PBMCs from participants who had different magnitudes of an *Mtb*-specific IFNγ response. Specifically, we compared the unstimulated samples from individuals expressing any of the three susceptibility alleles and who had a magnitude of response below the median for individuals included in the RNAseq (<568 SFCs) with individuals lacking those alleles and who had a magnitude of response above the median. Unstimulated samples were used for this analysis to avoid detecting genes differentially expressed solely due to the *Mtb*-specific magnitude i.e. potentially higher levels of T cell activation genes, etc. This analysis revealed 21 DEGs which were all upregulated in individuals without the susceptibility alleles (**Figure 3a**). This gene signature could distinguish between these two cohorts (with alleles; low expression and without alleles; high expression, **Figure 3b**,**c**). Furthermore, 5 of these 21 DEGs overlapped with the 63-gene signature described above (CXCL2, JUN, TNFAIP6, ZNF66, and AC068896.1). ![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/07/10/2022.07.06.22277293/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/F3) Figure 3. Differences in *Mtb*-specific magnitude of response is associated with a specific gene signature. a) Volcano plot showing differentially expressed genes comparing individuals with susceptibility alleles and low magnitude of expression (n=9) to those without and high magnitude of expression (n=13) in unstimulated samples. Black indicates DEGs (adjusted p value <0.05 and log2 fold change >1 and <-1). b) PCA plot illustrating differences between the two cohorts for the 21 DEGs. c) Heatmap displaying the 20 DEGs from (a). Genes are ordered by hierarchical clustering, and samples are ordered according to the two cohorts. Taken together, individuals expressing any of the susceptibility alleles have a specific gene signature that is enriched in genes related to APCs, possibly related to the lower magnitude of *Mtb*-specific IFNγ responses detected. ## Discussion In this study, we examined the potential association between the expression of specific HLA alleles and ATB disease. The identification of risk factors for ATB will help identify individuals at high risk of progression and could inform preventative treatment strategies. Using samples from 636 participants, we identified 3 HLA alleles associated with susceptibility for ATB. The expression of these alleles was associated with a lower magnitude of *Mtb*-specific T cell responses and also correlated with a lower expression of APC related genes. HLA alleles are highly polymorphic in general, and the genetic heterogeneity in HLA alleles also varies across different ethnic groups and geographical populations. However, the alleles identified as associated with differential disease outcomes in the present study are relatively common in most ethnic groups. While we identified an association of ATB with DQA1*03:01 in our study cohort from South Africa and the US, this allele was also found in association with pulmonary tuberculosis in a cohort from Iceland and Russia 25. We further identified the DQB1 alleles most likely expressed together with DQA1*03:01. These DQB1 alleles were DQB1*02:02, *04:02, *03:03, and *06:01/02. A meta-analysis of 12 studies found that HLA DQB1*06:01 was significantly associated with an increased risk of ATB as well 26. Sveinbjornsson et al. reported the association of the HLA class II allele DQA1*03 with *Mtb* infection and identified a noncoding variant, rs557011, located between the HLA-DQA1 and HLA-DRB1 alleles that also associated with *Mtb* infection 25. Many other groups have reported similar findings for different class I and class II HLA alleles 27-30. Due to the polymorphic nature and heterogeneity of HLA allele expression in diverse geographical locations, further studies in additional cohorts across different ethnicities may identify or confirm additional HLA alleles associated with susceptibility to ATB. In addition, HLA alleles present a strong linkage disequilibrium in which preferential combinations of alleles are inherited together in the genome more often than expected. The series of alleles at linked loci on a single chromosome is named “haplotype” 31. Specific haplotypes of association with ATB may be identified through more detailed HLA typing in which these linked loci are taken into consideration. In our study we did not determine these linkages. We further linked expression of the susceptibility alleles to changes in the *Mtb*-specific immune response. Individuals expressing DQA1*03:01, DPB1*04:02, and DRB4*01:01 had a lower *Mtb*-specific immune response compared to those lacking these alleles. Other groups have reported impaired immune responses due to the expression of specific HLA alleles in patients with human immunodeficiency virus (HIV) and hepatitis c virus (HCV) infection 32-35. Huang et al. identified an HLA-B*35 variant that was associated with accelerated HIV infection in a population of patients of European descent. They described that individuals who expressed the HLA-B*35 susceptibility allele had significantly impaired dendritic cell (DC) function than individuals lacking the susceptibility allele and hypothesized that this impaired DC function was the mechanism by which the susceptibility allele contributed to HIV disease 33. Additionally, HLA class I and class II variants have been found to confer protection against severe HCV infection and these protective alleles have been reported to better present immunodominant HCV epitopes 32,34,35. Finally, HLA-DQ variants have been reported to present streptococcal and staphylococcal superantigens at different strengths and influence the overall T cell-mediated response to bacterial superantigens 36. Overall, HLA allele variants have been linked previously to altered immune responses in the context of different diseases. In addition, using transcriptomics we defined gene signatures that could distinguish individuals with vs. without expression of the susceptibility alleles. The DEGs could be divided into 4 clusters that provided further insight into the role of the HLA susceptibility alleles and *Mtb* infection. Two clusters corresponded to genes that were differentially expressed in response to the specific stimuli unstim and MTB300 vs. WCL stimulation. One cluster was comprised of genes that were upregulated in individuals who express the susceptibility alleles. These genes included the transcription factors CEBPD and CEBPB, which are highly expressed on monocytes ([www.dice.org](http://www.dice.org) and [www.proteinatlas.org](http://www.proteinatlas.org)). Finally, the last gene cluster was comprised of genes that were upregulated in individuals who are not carriers of the susceptibility alleles. This cluster included both HLA-DQA1 and HLA-DQB1, as well as PAX5, JUN and CLEC17A which are all highly expressed on B cells ([www.dice.org](http://www.dice.org)). JUN is involved in toll-like receptor activity and antigen-cross presentation 37, and downregulation (as in individuals expressing the susceptibility alleles) leads to impairment of those functions. This cluster also contained TIRAP, which encodes for a molecule that is associated with toll-like receptors and also plays an important role in innate immunity 38. Thus, individuals expressing the susceptibility alleles have downregulation of important immune response related genes, which may contribute to an increased risk of ATB. Lastly, we identified a gene signature related to the magnitude of *Mtb*-specific T cell responses. This signature was upregulated in individuals lacking the susceptibility alleles and 5 genes overlapped with the signature detected when comparing DEGs irrespective of the *Mtb*-specific T cell response detected. This signature included JUN, as mentioned above, as well as, CCL19 which is involved in the homing of lymphocytes to the lung in *Mtb* infection 39. The signature also included CXCL2 which is a neutrophil-attracting chemokine 40. Taken together, these gene signature results suggest that the susceptibility alleles potentially increase the risk for ATB by impairing antigen loading and presentation to immune cells for proper clearance of *Mtb* infection. In conclusion, our study identified three HLA alleles (HLA DQA1*03:01, DPB1*04:02, and DRB4*01:01) associated with susceptibility to ATB. Further investigation revealed that individuals expressing these alleles have a decreased *Mtb*-specific immune response. The identified HLA alleles and associated gene signatures can help identify individuals who are at risk of ATB. ## Supporting information Supplemental Table 1 [[supplements/277293_file02.xlsx]](pending:yes) Supplemental table 2 [[supplements/277293_file03.xlsx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## Author contributions TJS, AS, BP, and CSLA participated in the design and direction of the study. LYC, RK, and CSLA performed literature search and participated in data collection, analysis and interpretation. VR and TJS and the SATVI-Ja Jolla Study Group recruited participants. EJP, SAM, and MMD coordinated and performed HLA typing. LYC and CSLA wrote the first draft. All authors reviewed, edited and approved the final version of the manuscript before submission. ## Declaration of interests The authors declare no competing interests. ## The SATVI study group View this table: [Table4](http://medrxiv.org/content/early/2022/07/10/2022.07.06.22277293/T4) ## Acknowledgements We thank the sequencing core and flow cytometry core at La Jolla Institute for Immunology for their help in these studies. This study was supported by the National Institutes of Health grants U19 AI118626, S10 RR027366, S10 OD016262. * Received July 6, 2022. * Revision received July 6, 2022. * Accepted July 10, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. WHO. Global Tuberculosis Report 2020. (2020). 2. Petruccioli, E., Scriba, T. J., Petrone, L., Hatherill, M., Cirillo, D. M., Joosten, S. A., Ottenhoff, T. H., Denkinger, C. M. & Goletti, D. Correlates of tuberculosis risk: predictive biomarkers for progression to active tuberculosis. Eur Respir J 48, 1751–1763, doi:10.1183/13993003.01012-2016 (2016). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjk6IjQ4LzYvMTc1MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA3LzEwLzIwMjIuMDcuMDYuMjIyNzcyOTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. Barry 3rd, C. E., Boshoff, H. I., Dartois, V., Dick, T., Ehrt, S., Flynn, J., Schnappinger, D., Wilkinson, R. J. & Young, D. The spectrum of latent tuberculosis: rethinking the biology and intervention strategies. Nat Rev Micro 7, 845–855 (2009). 4. Pai, M., Behr, M. A., Dowdy, D., Dheda, K., Divangahi, M., Boehme, C. C., Ginsberg, A., Swaminathan, S., Spigelman, M., Getahun, H., Menzies, D. & Raviglione, M. Tuberculosis. Nature reviews. Disease primers 2, 16076, doi:10.1038/nrdp.2016.76 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrdp.2016.76&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27784885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 5. WHO. The End TB Strategy. (2017). 6. 1. Anaya JM, 2. Shoenfeld Y, & 3. Rojas-Villarraga A Cruz-Tapias P, Castiblanco J & JM, A. in Autoimmunity: From Bench to Bedside (eds Anaya JM, Shoenfeld Y, & Rojas-Villarraga A) (El Rosario University Press, 2013). 7. Gao, J., Zhu, C., Zhu, Z., Tang, L., Liu, L., Wen, L. & Sun, L. The human leukocyte antigen and genetic susceptibility in human diseases. Journal of Bio-X Research 02, 112–120, doi:10.1097/JBR.0000000000000044 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/JBR.0000000000000044&link_type=DOI) 8. Blackwell, J. M., Jamieson, S. E. & Burgner, D. HLA and infectious diseases. Clinical microbiology reviews 22, 370–385, Table of Contents, doi:10.1128/CMR.00048-08 (2009). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjg6IjIyLzIvMzcwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDcvMTAvMjAyMi4wNy4wNi4yMjI3NzI5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. Dendrou, C. A., Petersen, J., Rossjohn, J. & Fugger, L. HLA variation and disease. Nat Rev Immunol 18, 325–339, doi:10.1038/nri.2017.143 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri.2017.143&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29292391&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 10. Sanchez-Mazas, A. A review of HLA allele and SNP associations with highly prevalent infectious diseases in human populations. Swiss Med Wkly 150, w20214, doi:10.4414/smw.2020.20214 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4414/smw.2020.20214&link_type=DOI) 11. Wang, L., Zou, Z. Q. & Wang, K. Clinical Relevance of HLA Gene Variants in HBV Infection. J Immunol Res 2016, 9069375, doi:10.1155/2016/9069375 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2016/9069375&link_type=DOI) 12. Gough, S. C. & Simmonds, M. J. The HLA Region and Autoimmune Disease: Associations and Mechanisms of Action. Curr Genomics 8, 453–465, doi:10.2174/138920207783591690 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2174/138920207783591690&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19412418&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000253956600005&link_type=ISI) 13. Eriksson, D., Royrvik, E. C., Aranda-Guillen, M., Berger, A. H., Landegren, N., Artaza, H., Hallgren, A., Grytaas, M. A., Strom, S., Bratland, E., Botusan, I. R., Oftedal, B. E., Breivik, L., Vaudel, M., Helgeland, O., Falorni, A., Jorgensen, A. P., Hulting, A. L., Svartberg, J., Ekwall, O., Fougner, K. J., Wahlberg, J., Nedrebo, B. G., Dahlqvist, P., Norwegian Addison Registry Study, G., Swedish Addison Registry Study, G., Knappskog, P. M., Wolff, A. S. B., Bensing, S., Johansson, S., Kampe, O. & Husebye, E. S. GWAS for autoimmune Addison’s disease identifies multiple risk loci and highlights AIRE in disease susceptibility. Nature communications 12, 959, doi:10.1038/s41467-021-21015-8 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-21015-8&link_type=DOI) 14. Tait, B. D., Colman, P. G., Morahan, G., Marchinovska, L., Dore, E., Gellert, S., Honeyman, M. C., Stephen, K. & Loth, A. HLA genes associated with autoimmunity and progression to disease in type 1 diabetes. Tissue Antigens 61, 146–153, doi:10.1034/j.1399-0039.2003.00013.x (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1399-0039.2003.00013.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12694582&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000182117300004&link_type=ISI) 15. McKinney, D. M., Fu, Z., Le, L., Greenbaum, J. A., Peters, B. & Sette, A. Development and validation of a sample sparing strategy for HLA typing utilizing next generation sequencing. Hum Immunol, doi:10.1016/j.humimm.2015.04.007 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.humimm.2015.04.007&link_type=DOI) 16. Thorstenson, Y. R., Creary, L. E., Huang, H., Rozot, V., Nguyen, T. T., Babrzadeh, F., Kancharla, S., Fukushima, M., Kuehn, R., Wang, C., Li, M., Krishnakumar, S., Mindrinos, M., Fernandez Vina, M. A., Scriba, T. J. & Davis, M. M. Allelic resolution NGS HLA typing of Class I and Class II loci and haplotypes in Cape Town, South Africa. Hum Immunol 79, 839–847, doi:10.1016/j.humimm.2018.09.004 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.humimm.2018.09.004&link_type=DOI) 17. Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., Batut, P., Chaisson, M. & Gingeras, T. R. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21, doi:10.1093/bioinformatics/bts635 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts635&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23104886&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312654600003&link_type=ISI) 18. Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864, doi:10.1093/bioinformatics/btr026 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr026&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21278185&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000288277300018&link_type=ISI) 19. Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., Marth, G., Abecasis, G., Durbin, R. & Genome Project Data Processing, S. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079, doi:10.1093/bioinformatics/btp352 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp352&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19505943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268808600014&link_type=ISI) 20. Liao Y S. G., Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930, doi:doi:10.1093/bioinformatics/btt656 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt656&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24227677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000334078300005&link_type=ISI) 21. Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15, 550, doi:10.1186/s13059-014-0550-8 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 22. Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. Series B (Methodological) 57, 289–300 (1995). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2346101&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:A1995QE4&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 23. Schmiedel, B. J., Singh, D., Madrigal, A., Valdovino-Gonzalez, A. G., White, B. M., Zapardiel-Gonzalo, J., Ha, B., Altay, G., Greenbaum, J. A., McVicker, G., Seumois, G., Rao, A., Kronenberg, M., Peters, B. & Vijayanand, P. Impact of Genetic Polymorphisms on Human Immune Cell Gene Expression. Cell 175, 1701–1715 e1716, doi:10.1016/j.cell.2018.10.022 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2018.10.022&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30449622&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 24. Paul, S., Dillon, M. B., Lindestam Arlehamn, C. S., Huang, H., Davis, M. M., McKinney, D. M., Scriba, T. J., Sidney, J., Peters, B. & Sette, A. A Population Response Analysis Approach To Assign Class II HLA-Epitope Restrictions. J Immunol 194, 6164–6176, doi:10.4049/jimmunol.1403074 (2015). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE5NC8xMi82MTY0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDcvMTAvMjAyMi4wNy4wNi4yMjI3NzI5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 25. Sveinbjornsson, G., Gudbjartsson, D. F., Halldorsson, B. V., Kristinsson, K. G., Gottfredsson, M., Barrett, J. C., Gudmundsson, L. J., Blondal, K., Gylfason, A., Gudjonsson, S. A., Helgadottir, H. T., Jonasdottir, A., Jonasdottir, A., Karason, A., Kardum, L. B., Knezevic, J., Kristjansson, H., Kristjansson, M., Love, A., Luo, Y., Magnusson, O. T., Sulem, P., Kong, A., Masson, G., Thorsteinsdottir, U., Dembic, Z., Nejentsev, S., Blondal, T., Jonsdottir, I. & Stefansson, K. HLA class II sequence variants influence tuberculosis risk in populations of European ancestry. Nat Genet 48, 318–322, doi:10.1038/ng.3498 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3498&link_type=DOI) 26. Li, C., Zhou, Y., Xiang, X., Zhou, Y. & He, M. The Relationship of HLA-DQ Alleles with Tuberculosis Risk: A Meta-analysis. Lung 193, 521–530, doi:10.1007/s00408-015-9747-1 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00408-015-9747-1&link_type=DOI) 27. Amirzargar, A. A., Yalda, A., Hajabolbaghi, M., Khosravi, F., Jabbari, H., Rezaei, N., Niknam, M. H., Ansari, B., Moradi, B. & Nikbin, B. The association of HLA-DRB, DQA1, DQB1 alleles and haplotype frequency in Iranian patients with pulmonary tuberculosis. Int J Tuberc Lung Dis 8, 1017–1021 (2004). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15305487&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223040300016&link_type=ISI) 28. Salie, M., van der Merwe, L., Moller, M., Daya, M., van der Spuy, G. D., van Helden, P. D., Martin, M. P., Gao, X. J., Warren, R. M., Carrington, M. & Hoal, E. G. Associations between human leukocyte antigen class I variants and the Mycobacterium tuberculosis subtypes causing disease. J Infect Dis 209, 216–223, doi:10.1093/infdis/jit443 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jit443&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23945374&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 29. Souza de Lima, D., Morishi Ogusku, M., Porto Dos Santos, M., de Melo Silva, C. M.,Alves de Almeida, V., Assumpcao Antunes, I., Boechat, A. L., Ramasawmy, R. & Sadahiro, A. Alleles of HLA-DRB1*04 Associated with Pulmonary Tuberculosis in Amazon Brazilian Population. PLoS One 11, e0147543, doi:10.1371/journal.pone.0147543 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0147543&link_type=DOI) 30. Wamala, D., Buteme, H. K., Kirimunda, S., Kallenius, G. & Joloba, M. Association between human leukocyte antigen class II and pulmonary tuberculosis due to mycobacterium tuberculosis in Uganda. BMC Infect Dis 16, 23, doi:10.1186/s12879-016-1346-0 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-016-1346-0&link_type=DOI) 31. Shiina, T., Hosomichi, K., Inoko, H. & Kulski, J. K. The HLA genomic loci map: expression, interaction, diversity and disease. J Hum Genet 54, 15–39, doi:10.1038/jhg.2008.5 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/jhg.2008.5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19158813&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265355700003&link_type=ISI) 32. Harcourt, G., Hellier, S., Bunce, M., Satsangi, J., Collier, J., Chapman, R., Phillips, R. & Klenerman, P. Effect of HLA class II genotype on T helper lymphocyte responses and viral control in hepatitis C virus infection. J Viral Hepat 8, 174–179, doi:10.1046/j.1365-2893.2001.00289.x (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2893.2001.00289.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11380794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000169351300003&link_type=ISI) 33. Huang, J., Goedert, J. J., Sundberg, E. J., Cung, T. D., Burke, P. S., Martin, M. P., Preiss, L., Lifson, J., Lichterfeld, M., Carrington, M. & Yu, X. G. HLA-B*35-Px-mediated acceleration of HIV-1 infection by increased inhibitory immunoregulatory impulses. J Exp Med 206, 2959–2966, doi:10.1084/jem.20091386 (2009). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjExOiIyMDYvMTMvMjk1OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA3LzEwLzIwMjIuMDcuMDYuMjIyNzcyOTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 34. Neumann-Haefelin, C., Spangenberg, H. C., Blum, H. E. & Thimme, R. Host and viral factors contributing to CD8+ T cell failure in hepatitis C virus infection. World J Gastroenterol 13, 4839–4847, doi:10.3748/wjg.v13.i36.4839 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3748/wjg.v13.i36.4839&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17828815&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 35. Wertheimer, A. M., Miner, C., Lewinsohn, D. M., Sasaki, A. W., Kaufman, E. & Rosen, H. R. Novel CD4+ and CD8+ T-cell determinants within the NS3 protein in subjects with spontaneously resolved HCV infection. Hepatology 37, 577–589, doi:10.1053/jhep.2003.50115 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/jhep.2003.50115&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12601356&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F07%2F10%2F2022.07.06.22277293.atom) 36. Llewelyn, M., Sriskandan, S., Peakman, M., Ambrozak, D. R., Douek, D. C., Kwok, W. W., Cohen, J. & Altmann, D. M. HLA class II polymorphisms determine responses to bacterial superantigens. J Immunol 172, 1719–1726, doi:10.4049/jimmunol.172.3.1719 (2004). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE3Mi8zLzE3MTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNy8xMC8yMDIyLjA3LjA2LjIyMjc3MjkzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 37. Novoszel, P., Drobits, B., Holcmann, M., Fernandes, C. S., Tschismarov, R., Derdak, S., Decker, T., Wagner, E. F. & Sibilia, M. The AP-1 transcription factors c-Jun and JunB are essential for CD8alpha conventional dendritic cell identity. Cell Death Differ 28, 2404–2420, doi:10.1038/s41418-021-00765-4 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41418-021-00765-4&link_type=DOI) 38. Belhaouane, I., Hoffmann, E., Chamaillard, M., Brodin, P. & Machelart, A. Paradoxical Roles of the MAL/Tirap Adaptor in Pathologies. Front Immunol 11, 569127, doi:10.3389/fimmu.2020.569127 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2020.569127&link_type=DOI) 39. Khader, S. A., Rangel-Moreno, J., Fountain, J. J., Martino, C. A., Reiley, W. W., Pearl, J. E., Winslow, G. M., Woodland, D. L., Randall, T. D. & Cooper, A. M. In a murine tuberculosis model, the absence of homeostatic chemokines delays granuloma formation and protective immunity. J Immunol 183, 8004–8014, doi:10.4049/jimmunol.0901937 (2009). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4My8xMi84MDA0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDcvMTAvMjAyMi4wNy4wNi4yMjI3NzI5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 40. Dorhoi, A., Iannaccone, M., Maertzdorf, J., Nouailles, G., Weiner, J., 3rd. & Kaufmann, S. H. Reverse translation in tuberculosis: neutrophils provide clues for understanding development of active disease. Front Immunol 5, 36, doi:10.3389/fimmu.2014.00036 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2014.00036&link_type=DOI)