Neutrophil extracellular trap formation and gene programs distinguish TST/IGRA sensitization outcomes among Mycobacterium tuberculosis exposed persons living with HIV

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.


Introduction
Worldwide, tuberculosis (TB) remains the leading cause of death by a single bacterial agent [1].People who are living with HIV (PLWH) have an increased risk of TB [2,3].We previously identified a cohort of PLWH, living in a community with high TB burden in the Western Cape, South Africa, who are persistently TB, tuberculin skin test (TST) and interferon gamma release immune assay (IGRA) negative (HITTIN) [4][5][6].TST and IGRA measure persistent adaptive immune memory to Mycobacterium tuberculosis (Mtb) protein antigens and are used as a surrogate for prior Mtb infection.Despite these persons having a history of previously low CD4+ T-cell counts, prior to antiretroviral therapy (ART), they display no canonical evidence of prior Mtb infection or history of disease [4].However, we did detect circulating Mtb-specific antibodies, despite absence of canonical T-cell memory to dominant Mtb antigens, suggesting these individuals may possess a unique immune mechanism of TB protection, not reliant on conventional T-cell help.
We hypothesized that the innate immune system, and specifically neutrophils (PMN), play an inherent role in the protective control of Mtb infection in HITTIN.PMN are the most phenome Archive (EGA) as "Neutrophils as effector cells in resistance to infection by Mycobacterium tuberculosis in HIV-infected individuals" with the study ID EGAS00001007262 (https://ega-archive. org/studies/EGAS00001007262).Qualified researchers can request access to the data by contacting the listed contact person on https://egaarchive.org/datasets/EGAD00001010893and completing a Data Access Request form which will be evaluated by the Data Access Committee who will consult with the Health Research Ethics Committee of Stellenbosch and the Stellenbosch University Contracts Office.Data access will be for approved use and will be governed by the provisions laid out in the associated informed consent of individual data subjects and of the Health Research Ethics Committee of Stellenbosch University.All investigators will have to sign a Data Transfer Agreement as stipulated by the Stellenbosch University Contracts Office.A qualified researcher refers to a senior investigator who is employed or legitimately affiliated with an academic, non-profit or government institution and who has a track record in the field.Researchers granted access to Project Data must feedback the results of their research to the Data Access Committee, after publication for distribution to the community representative if required.Access is conditional upon availability of data and on signed agreement by the researcher(s) and the responsible employing institution to abide by the policies and conditions related to publication, data ownership, data return, intellectual property rights, data disposal, ethical approval, confidentiality and commercialization referred herein.
abundant leukocytes and among the first responders to Mtb infection in the lung in animal models as well as humans [7].They are armed with an arsenal of antimicrobial granules known to restrict Mtb growth and are key players in the inflammatory response against Mtb [8][9][10][11].PMN can control Mtb growth during acute infection [8,11,12].Household pulmonary TB contacts with higher initial peripheral neutrophil counts were less likely to become infected with Mtb [11].Despite lower RNA expression in PMN compared to other innate immune cells, pathogen-triggered gene expression changes underlie microbial responses by PMN [13,14].
The mechanism of cell death can influence the outcome of Mtb infection control.Neutrophils control inflammation, predominantly through apoptotic and cell clearance mechanisms with efferocytosis of Mtb-infected apoptotic neutrophils by macrophages favoring a beneficial host outcome [14][15][16][17][18]. Neutrophil extracellular traps (NETs) were recently highlighted as a potentially important mechanism of resistance to TST/IGRA conversion in a Ugandan cohort [19].As a primary cell death effector mechanism, neutrophils decondense their chromatin following histone modification, unraveling their DNA, which is then coated in cytoplasmic granular proteins before the DNA/histone/protein complex is expelled from the lysing neutrophil.This released web-like structured NET is able to capture extracellular microbes and kill them using the attached granule-derived antimicrobial peptides and histones [20,21].Mtb-induced NETs are phagocytosis and reactive oxidative species (ROS) dependent, but have been suggested to lack Mtb microbicidal activity [22,23].Recent studies point to potential crosstalk of NETs and forms of necrotic cell death pathways such as pyroptosis, despite the classical definition of NETosis as a unique mechanism of cell death [24].Necrotic neutrophil cell death, uncontrolled type 1 interferon (IFN) responses and abundant ROS drive a hyperinflammatory neutrophil phenotype and contribute to TB severity [18,[25][26][27][28].
The aim of this study was to gain insight into mechanisms underlying early neutrophil responses to Mtb in HITTIN persons with apparent protection from Mtb infection as inferred from the absence of canonical T-cell memory.Specifically, we tested if neutrophils from HIT-TIN study participants (PMN HITTIN ) displayed a transcriptional response to Mtb infection that was significantly different from that found in neutrophils from PLWH who have a robust Mtb-specific T-cell response, testing persistently TST and IGRA positive, but also with no history of TB (HIT).We demonstrated that the in vitro Mtb-induced gene expression changes in neutrophils from the HITTIN and HIT participants were significantly different.Genes differentially regulated by Mtb between the two groups were enriched for key regulators of ROS and NET formation.Fluorescence microscopy images corroborated the transcriptomic findings with PMN HIT demonstrating a greater response to Mtb with more NETs forming after infection.PMN HITTIN responded with less NETs and lower counts of key genes associated with nicotinamide adenine dinucleotide phosphate hydrogen (NADPH) activity and ROS formation possibly contributing to greater effector control of Mtb.

Study participants and samples
Individuals included in the HITTIN and HIT groups were a part of the large ResisTB cohort [4].Neutrophils obtained from 17 HITTIN and 11 HIT individuals, all of whom were PLWH and on ART, were used in the final analysis (Table 1).These individuals were part of stringently defined cohorts living in a high TB burden community who, despite low CD4+counts before ART initiation, never developed TB.
The average age of participants in the HITTIN group was 43.71 (±6.43) years and that of HIT was 44.09 (±6.89) (Table 1).There were no significant differences between the ratio of females to males in each group (p = 0.6, Fisher's exact test), age (p = 0.9, Wilcoxon rank sum test), or BMI (p = 0.6, Wilcoxon rank sum test) of the participants.All but one participant had CD4+ counts of < 200 cells/mm 3 prior to starting ART.Participants had controlled viral loads, with participants in the HITTIN group having been on ART for 7.35 (±3.15) and HIT for 8.97 (±3.2) years.Although HITTIN was on ART for a shorter period there was no significant difference to the time spent in ART in HIT (p = 0.2, Wilcoxon rank sum exact test).One participant was on isoniazid Mtb prophylactic therapy (IPT) at the time of the enrolment, this was corrected for during the analysis.All PLWH receive IPT according to national guidelines in South Africa.There was no significant difference between the groups in terms of having taken prior IPT and time spent on IPT (p = 0.2, Fisher's exact test).HITTIN and HIT have a similar distribution of chronic disease (p = 0.4, Fisher's exact test), and social habits such as cigarette smoking (p = 0.5, Fisher's exact test), alcohol use (p = 0.2, Fisher's exact test) and cannabis use (p = 0.4, Fisher's exact test).

Gene expression analysis of PMN in HITTIN and HIT in response to Mtb
Neutrophils were isolated from whole blood and infected with Mtb H37Rv for 1h and 6h.Purity of isolated neutrophils was confirmed by flow cytometry, with similar proportions of CD45+CD15+CD66b+CD16+ neutrophils isolated for HITTIN and HIT individuals (p = 0.28, Wilcoxon rank sum test) (S1 and S2 Tables).RNA was extracted from uninfected and infected neutrophils (PMN) at 1h and 6h timepoints and differential gene expression analyses performed by RNAseq.Significant differentially expressed genes (DEGs) were defined by an absolute cutoff of a log 2 FC of 0.2 and a false discovery rate (FDR) adjusted p-value of 5% (Tables 2 and S3).
Mtb infection triggered significant gene expression changes when compared to uninfected PMN at 1h: 151 up-and 40 downregulated genes for PMN HITTIN while 98 genes were up-and 11 were downregulated for PMN from the HIT group (PMN HIT ) (S1A and S1B Fig, S3   When comparing transcriptional responses of infected PMN HITTIN and PMN HIT directly, no significant DEG were identified at 1h (S1C Fig, S3 Table ).However, contrasting the effect differences for Mtb infection compared to no infection at 6h for each group (PMN HITTIN vs PMN HIT ), we identified 2285 genes with significant differential response between the two groups (Fig 1C, S3 Table ).Since the 1h time point showed limited differences in gene induction by Mtb, we focused on the 6h time point differences between PMN phenotypes in subsequent analyses.

Differential gene expression analysis showed a lower overall fold change difference during Mtb infection of PMN from HITTIN compared to HIT
We next compared the overall difference in transcriptomic response to Mtb between the two PMN groups.Consistent with different numbers of DEGs identified for each phenotypic group, when we evaluated the statistical significance of expression changes between the HIT-TIN and HIT groups, we identified an overall dampened 6h transcriptomic response in PMN HITTIN compared to PMN HIT (p < 2.2e 16

Pathway and gene ontology (GO) enrichment analysis for DEGs between PMN HITTIN and PMN HIT after 6h Mtb infection
We next investigated the 2285 genes with significant differential transcriptional response of PMN HITTIN and PMN HIT to Mtb after 6 hours infection by examining the different pathways and analyzing the biological processes that are characterized by these genes.We conducted separate term enrichment analyses of up-and downregulated genes and a third analysis that considered both up-and downregulated genes.GO terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome pathways were considered significant if a term was enriched for at least 5 DEGs and a FDR cutoff of 10%.An overview of the results is shown in Table 2.By focusing on all DEGs significantly different between PMN HITTIN and PMN HIT , we observed a total of 496 enriched KEGG, Reactome and GO terms (Table 2).When evaluating genes more strongly triggered in PMN HITTIN we detected 29 terms.Conversely, when focusing on genes less strongly triggered in PMN HITTIN and comparatively to PMN HIT , we detected 719 terms (Table 2).
Manhattan plots for the three term analyses at 6h are shown in Fig 3 with significantly different terms and pathways of interest indicated.Amongst the enriched terms in PMN HITTIN , were "Apoptosis", "Neutrophil extracellular trap formation", and "NADPH regeneration", which are terms that directly relate to microbicidal activity of PMN HITTIN (Fig 3B).Although  3), the overall enrichment of the terms is also influenced by genes which are less up-and downregulated (Fig 3A, Table 3).This was also the case for terms with genes triggered relatively less strongly in PMN HITTIN compared to PMN HIT .These were dominated by genes involved in neutrophil chemotaxis, neutrophil degranulation, necroptosis and necrotic death (Fig 3C).It is important to note that the overall response to Mtb was always lower in PMN HITTIN compared to PMN HIT .a The average log 2 FC Mtb infection compared to the uninfected response for neutrophils from HITTIN at 6h. b The average log 2 FC Mtb infection compared to the uninfected response for neutrophils from HIT at 6h. c The average log 2 FC response difference between neutrophils from HITTIN and HIT in response to Mtb infection at 6h (interaction test).
d The adjusted p-value after the Benjamini Hochberg correction for multiple testing for interaction test looking at the differential response of HITTIN to HIT at 6h infected ("PMN HITTIN x PMN HIT * ") Significant genes were defined as genes with an absolute log 2 FC � 0.2 and adjusted p value � 0.05 https://doi.org/10.1371/journal.pgen.1010888.t003

NET area change difference between HITTIN and HIT from 1 to 6h after Mtb infection
To directly evaluate the biological outcome in the total amount of NETs observed between HITTIN vs HIT, as is suggested by our transcriptional enrichments, we used fluorescence microscopy to quantitate NET formation.We stained fixed cells which were processed in parallel with the cells used for RNAseq.NETs were stained using anti-H2AH2B/DNA (PL2-3) which detects decondensed chromatin and nuclear DNA stained with Hoechst, and the area of both features was quantified.
As a measure of cellular viability over time we compared the total change in cell nuclei area between cells fixed at 1h and 6h of infection with Mtb.There was no significant two way interaction between PMN HITTIN and PMN HIT and the infection status, F(1,12) = 1.1870, p = 0.30, including a similar trend in direction of response for both PMN.However, there was a significant Mtb infection effect, F(1,12) = 9.7290, p = 0.009, with pairwise comparisons showing a significantly greater decrease in total cell nuclei area between 1h and 6h after Mtb infection for PMN HIT compared to PMN HITTIN (p = 0.04, pairwise t-test) (Figs 4 and S3).
When then comparing difference in NET area at 1h vs 6h, there was a statistically significant interaction between the PMN groups and infection status, F(1, 10.5924) = 5.3398, p < 0.0421).The simple main effect of phenotype group (considering the Bonferroni adjusted p-value) was significant for Mtb infection (p = 0.0007), but not for non-infection (p = 1).Consistent with the greater viability of PMN HITTIN at 6h of infection, pairwise comparisons show that PMN HITTIN also induce a significantly smaller change in NETs produced between 1h and 6 h of infection, compared to PMN HIT (p = 0.0003) (Figs 4 and S3).The transcriptional response for "Neutrophil extracellular trap formation" was driven by genes with significant positive FC (less downregulated), however the imaging data revealed overall lower NET formation in PMN HITTIN .These findings were in line with the overall lower transcriptional response observed in PMN HITTIN.

DEGs in HITTIN vs HIT after 6h Mtb infection in the "Neutrophil extracellular trap formation" pathway
We next investigated the DEGs associated with the "Neutrophil extracellular trap formation" enriched in the combined KEGG pathway (Fig 3A, Table 3).In particular, the lower transcriptional and NET response in PMN HITTIN prompted further evaluation of genes which showed a lower upregulated response to Mtb in PMN HITTIN compared to PMN HIT .Compared to PMN HIT , Mtb infection of PMN HITTIN , triggered a lower upregulation of genes involved in the multiple-protein NADPH oxidase complex including Rac family small GTPase 2 (RAC2), and the transmembrane catalytic [cytochrome b-245 -alpha (CYBA) and -beta (CYBB)].NADPH oxidase Nox2 (encoded by CYBB) and other cellular NADPH oxidases is involved with ROS production and NET formation [29][30][31].CYBB is a nuclear factor kappa B (NF-κB) transcriptional target and together with NFKB1 was also less upregulated in PMN HITTIN .Interestingly, NCF2, the gene encoding neutrophil cytosolic factor 2 (NCF-2 or p67-phox) was downregulated in PMN HITTIN and upregulated in PMN HIT .PMN HITTIN , also showed enrichment for DEGs related to "NADPH regeneration" (Fig 3A and 3B).NADPH can dually aid in ROS detoxification or production and is key for ROS mediated NET formation [30,32,33].
Histone deacetylases (HDACs) play a key role in NET formation and allow for peptidylarginine deiminase 4 (PAD4) mediated histone citrullination the initial step in chromatin decondensation [34,35].Compared to PMN HIT, PMN HITTIN had a less downregulated response in HDAC1, HDAC3, HDAC4 and PADI4 at 6h of Mtb infection.Gasdermin D (GSDMD) which plays a key role perforating the nuclear membrane to aid release of the decondensed chromatic during NET formation [36,37], also displayed the same pattern of expression regulation.Caspase 1 (CASP1) and 4 (CASP4) which activate GSDMD have a lower upregulation in PMN HITTIN compared to PMN HIT after 6h Mtb infection [38].
Other DEGs enriched in the KEGG NET pathway, are also involved in additional neutrophil functional responses.Cell membrane receptors TLR2 and TLR4 were less upregulated in PMN HITTIN in response to 6h Mtb infection.Downstream of TLR4, pathway activation of NF-κB, Protein Kinase B (AKT) and phosphoinositide 3-kinase (PI3-K) lead to pro-survival mechanisms [39].Integral to this TLR signaling system is mitogen-activated protein kinase (MAPK) and PI3-K.Dysregulation in especially the PI3-K/AKT signaling system contributes to an imbalance in neutrophil chemotaxis and can heighten inflammation and decrease pathogen clearance [40,41].MAP2K2 was less upregulated while MAPK1, MAPK3, AKT1, AKT2 and PIK3CD were less downregulated in PMN HITTIN .Azurocidin 1 (AZU1), the only antimicrobial peptide gene also included in the NET term, was downregulated after 6h Mtb infection in PMN HITTIN whilst upregulated in infected PMN HIT (Table 3).

Discussion
HITTIN remain persistently TST, IGRA and TB negative despite prolonged Mtb exposure and antibody evidence of prior infection.This suggests HITTIN have different protective mechanisms in the early response to Mtb infection to that which occurs in HIT who convert to a positive TST and IGRA following Mtb infection.Here we investigated the DEGs of PMN HITTIN to determine whether they are a distinct and previously undefined group that contribute to HIT-TINs unique ability to control Mtb infection, prevent progression to TB, and interact with the adaptive immune response while possibly limiting persisting Mtb-specific interferon gamma (IFN-γ) T-cell memory responses.Using RNAseq analysis of ex vivo Mtb infected PMN we found that PMN HITTIN had an overall lower transcriptional FC response to Mtb after 6h of infection relative to PMN HIT.Positively enriched terms and pathways included apoptosis, NETosis, NADPH regeneration and ROS formation, with pathways related to necrotic cell death, necroptosis, neutrophil chemotaxis, degranulation and immune exhaustion which were triggered less strongly in PMN HITTIN compared to PMN HIT after 6h infection with Mtb.Using fluorescence microscopy, we demonstrate that after 6h of infection PMN HITTIN have undergone less NETosis than PMN HIT, corroborating the overall lower transcriptional response to Mtb infection after 6h seen in PMN HITTIN .The question arises whether the less pronounced response by PMN HITTIN facilitates rapid Mtb control and a more contained neutrophil response, thereby preventing overt damage and an exacerbated inflammatory cascade, associated with TB.Their relatively lower induction of genes in response to Mtb suggests this is likely, requiring future functional studies for confirmation.
Neutrophil degranulation and NET formation processes were observed 1-6 months prior to persons developing TB and play an important role in Mtb infection progression to disease as well as potential lung destruction [42].This highlights the important regulatory mechanism of NET formation by neutrophils in the inflammatory response against Mtb infection.
Mtb-induced ROS triggered necrosis in neutrophils and decreased the ability of macrophages to control Mtb growth [18].Transcription of different groups of genes are regulated by NETosis-specific kinase cascades [43].NETs are triggered by a NADPH/ROS dependent or an independent mechanism through TLR2/TLR4/lipopolysaccharide (LPS) activation [44,45].The multi-protein NADPH oxidase complex consists of multiple subunits (including cytosolic p40-phox, p47phox and p67-phox, as well as the catalytic gp91-phox subunit) and other proteins (transmembrane p22-phox and nucleotide-binding Rac2) and is involved in the process to produce ROS [46].ROS is mostly known to modulate a pro-inflammatory effect, with increased levels positively correlated to higher levels of TNF and a greater control of Mtb infection [47].High levels of ROS can be damaging, but at lower levels ROS can mediate an anti-inflammatory effect in mouse neutrophils by inhibiting p-AKT and NF-κβ and inducing apoptosis [48][49][50][51][52][53].In our data, within the "Neutrophil extracellular trap formation" pathway in KEGG, several significantly under-stimulated genes in the pathway suggest potential relatively lower ROS levels in PMN HITTIN which could account for the relatively lower NET formation after 6h infection with Mtb.Another alternative mechanism PMN HITTIN displayed for inflammation control was enrichment of "NADPH regeneration".NADPH may play a role in neutralizing rather than producing ROS in PMN HITTIN.
The lower upregulation of CYBA, CYBB (aka NOX2) as well as RAC2 in PMN HITTIN could point to a mechanism to decrease ROS formation and prevent or decrease necrosis with consequently lower ROS mediated NET formation.Strikingly, NCF2 was downregulated in PMN HITTIN and upregulated in PMN HIT .The G allele of NCF2 rs10911362 is associated with protection against TB in the Western Chinese Han population [54].The rs3794624 polymorphism in CYBA was also linked with decreased TB susceptibility in Chinese Han persons [46,54].PMN HITTIN exhibit an interesting potential NCF2-mediated mechanism driving less ROS and NET formation resulting in lower necrotic cell death.
TLR2/4 signalling could also mediate NET formation independent of ROS.A lower upregulation of TLR2 and TLR4 in PMN HITTIN likely balance potential inflammatory effects of NET release with downstream upregulation of apoptosis.NETs can also be triggered by activated CASP1 and ELANE cleaving GSDMD [55].CASP1 was less up-and GSDMD less downregulated in PMN HITTIN .Although NETs have not been shown to control Mtb infection, the granular components associated with NETs could contribute [56].Azurocidin 1 (AZU1) can interact with pentraxin 3 (PTX3) to enhance microbial function.The mechanism of this interaction is not known.PTX3 is a soluble pattern recognition receptor with increased levels seen in sepsis and TB [57].PTX3 was significantly upregulated by Mtb infection in both PMN HITTIN and PMN HIT but there was not a statistically significant difference between these.The fact that DEGs for other granule proteins were not identified is likely because they are mostly transcribed in immature neutrophils and this degranulation is reflective of their potentially different action [58][59][60].
HDAC1 has been implicated as a key regulator of innate immunity in monocytes isolated from HIV negative persons who tested persistently TST negative after household TB exposure [61].HDAC inhibitors in macrophages such as phenylbutyrate improved Mtb control [61,62].This inhibition was synergistically improved with vitamin D [62].Less downregulation of HDAC1 in PMN HITTIN translated into a decreased inflammatory response, with a lower upregulation of genes in the NF-kappa β signaling pathway (Fig 3A and 3C).This corroborates a more repressive or anti-inflammatory role in PMN HITTIN .
Increased NETs can potentiate local damage due to release of granular components, but NET formation has also been shown to limit inflammation through degradation of chemokines and cytokines [63].This intricate balance is likely maintained by PMN HITTIN , which demonstrate lower NET formation in response to Mtb infection despite a positive enrichment of NET-related genes as compared to PMN HIT .The formation of NETs itself is likely not problematic as PMN HITTIN could potentially effectively localize and trap Mtb through NETs and activate other innate cells such as macrophages [22].The balance of NETosis with other cell death mechanisms may play an additional role to determine the fate of Mtb.

Study limitations
Use of ART can significantly affect the transcriptional responsiveness of alveolar macrophages to Mtb.It is important to note that although HITTIN were on ART for a shorter time, there was no significant difference to the time spent on ART in HIT and therefore ART is an unlikely factor driving the transcriptional differences observed between neutrophils from HIT-TIN and HIT [64].
The effect of contaminating T cells in PMN cultures cannot be completely excluded.T cells could possibly contribute to IFN-γ driven PMN transcriptional differences with IFN-γ known to upregulate key genes affecting NADPH activity [65,66].Per clinical classification, HITTIN lack IFN-γ T cell responses as measured by IGRA.In addition, good PMN purity (90%) vs 5% in T cells, as well as the 6h timepoint which is short for the IFN-γ driven responses, argue against the significant impact of contaminating cells.Further classification of non IFN-γ T cell subsets in HITTIN is needed.
To account for large inter-individual differences in response to Mtb, we employed subjectspecific fold changes.This approach precludes comparisons of group expression levels.Future work is needed to determine if there is a measured difference in ROS released by PMN HITTIN compared to PMN HIT .addition, neutrophil Mtb killing assays will determine if PMN HITTIN show improved Mtb infection clearance.Multiple studies have identified neutrophil subpopulations and it is possible that a specific subpopulation could be driving the response differences we observed [67][68][69].Flow cytometry for identification of potential subpopulations as well as single cell RNA sequencing would be highly informative.

Conclusion
In general, for TB, neutrophils are mostly linked to hyperinflammatory responses and more severe disease.Here we showed a distinctive gene expression profile for neutrophil responses from HITTIN individuals, who appeared protected from TB despite a lack of canonical Mtbspecific T-cell memory.These findings put neutrophils at the forefront of potential innate immune cell mechanisms of Mtb infection resistance.They highlight a distinct phenotypic response to Mtb in PMN HITTIN compared to neutrophils from persons otherwise defined as sensitized by Mtb.
We measured a reflection of the overall impaired transcriptional response and consequently no claim can be made about specific effect.Rather, we are reporting a less pronounced response to Mtb infection in HITTIN after 6h that also effects pathways known to represent anti-mycobacterial activity.Contrary to Mtb inducing necrotic cell death mechanisms as in most neutrophils, DEGs in PMN HITTIN showed decreased transcriptional responses for necrosis, with an enrichment of terms related to apoptotic cell death and NETosis.Fluorescence imaging corroborated significant reduction in NET formation between PMN HITTIN and PMN HIT, likely driven by lower ROS transcriptional pathways with a downregulation of NCF2, a key mediator in ROS formation, in PMN HITTIN .These molecular data implicated neutrophils as key effector cells in Mtb infection resistance.Further increased understanding of the crucial pathways of Mtb infection control highlighted in the study could be harnessed for the development of Mtb prevention and treatment strategies.

Ethics statement
The study was approved by the Health Research Ethics Committee of Stellenbosch University (S18/08/175(PhD)).Samples used in this study were leveraged from the ResisTB study.The Health Research Ethics Committee of Stellenbosch University (N16/03/033 and N16/03/033A) and the Faculty of Health Sciences Human Research Ethics Committee of the University of Cape Town (755/2016 and 702/2017) approved participant recruitment for ResisTB.Additional approval was obtained from the City of Cape Town and Western Cape government for access to the relevant clinics.Formal consent for research participation was given by each participant by signing written informed consent documents.

Participant recruitment
Samples were collected from 29 participants recruited to the ResisTB study in Cape Town, South Africa.The ResisTB study recruited a group of participants between the ages of 35-60 years old.The recruitment of this group for the ResisTB study has been fully described previously [4].Briefly, the participants had to be living with HIV in an area of high Mtb transmission and have no history of previous or current TB.They had to have a history of living with a low CD4+ count (either with two CD4+ <350 cells/mm 3 counts at least 6 months apart or a single CD4+ count <200 cells/mm 3 ) prior to initiating antiretroviral therapy (ART), and be immune reconstituted on ART for at least one year at time of enrolment with the last CD4 + count >200 cells/mm 3 [4].Exclusion criteria included pregnancy, previous TB, symptoms suggestive of active TB disease, participation in other interventional studies, and any AIDS defining illness in the year prior to enrolment.Participants were seen at three visits.During the enrolment visit, blood samples were taken for an IGRA using the QuantiFERON-TB Gold Plus (QFT-Plus) in tube test.For follow up, whole blood was collected in an EDTA tube for Ficoll gradient separation and neutrophil isolation as well as a second IGRA, followed by TST administration with PPD RT23 (Staten Serum Institute).After 3 days the TST reading was taken and participants who tested IGRA negative from the first visit had bloods taken for a third and final IGRA.This ensured that any T-cell response that would be boosted by TST administration would be identified by IGRA to ensure participants with low level T-cell memory were correctly identified.
For this study, samples were leveraged from 17 older HITTIN and 11 HIT participants.Age was used as a surrogate for increased exposure frequency to Mtb since most persons are infected with Mtb by the age of 30-35 in the Western Cape of South Africa, where the study was conducted [70,71].

Neutrophil isolation
Neutrophils were isolated by Ficoll gradient separation.Whole blood diluted 1:1 with 1x Phosphate buffered saline (PBS, Sigma-Aldrich, USA) was layered over the density gradient separation medium (Histopaque / Ficoll-Paque, Sigma-Aldrich, USA).Cells were centrifuged for 25 minutes at 400 x g.After this, cells were washed twice with 4˚C PBS and centrifuged each time at 400 x g for 10 minutes at 4˚C.Peripheral blood mononuclear cells were removed first and then the remaining plasma and Ficoll-Paque layer.A red blood cell (RBC) lysis buffer (component concentrations) was added in a ratio of 1:10 or topped up to 50ml, if a 50 ml centrifuge tube was used, to the remaining bottom layer and incubated for 10 minutes at 4˚C.RBC lysis buffer (8% Ammonium chloride [NH 4 Cl], 0.8% Sodium bicarbonate [NaHCO 3 , Sigma-Aldrich, USA] and 0.4% Ethylenediaminetetraacetic acid [EDTA, Sigma-Aldrich, USA]).
After RBC lysis, samples were centrifuged at 400 x g for 10 minutes at 4˚C and then washed twice with PBS (4˚C) and centrifuged at 400 x g for 10 minutes at 4˚C.After the cell count, 1 x 10 7 cells in RPMI-1640 with L-glutamine and sodium bicarbonate (Sigma-Aldrich, USA), were seeded evenly over each row of 3 wells in a 6-well plate cell culture plate (Nest Scientific USA Inc., USA) at 0.33 x 10 7 cells per well for the RNA sequencing experiment and 2 x 10 5 cells per well in a 96-well (M0562, Greiner, Sigma Aldrich, USA) for microscopy.For each participant two 6-well plates (one for each time-point i.e., 1h and 6h, with 3 wells for no infection and 3 wells for Mtb infection) were seeded.Cells were incubated for 1h-2h at 37˚C and 5% CO 2 .An aliquot of remaining neutrophils was fixed with 4% Paraformaldehyde (PFA, cat.no.43368, Alfa Aesar, USA) and then stored overnight at -80˚C before transferred to and stored in liquid nitrogen.

Staining for flow cytometry
The fixed neutrophil aliquots were thawed and washed in Dulbecco's phosphate buffered saline (DPBS, Cat.No. 14190144, Thermo Fisher Scientific, Australia).1.5 x 10 6 cells per participant were resuspended in surface staining buffer comprised of DPBS + 3% Foetal Bovine Serum (FBS, Cat.No. SFBS-AU, Bovogen Biologicals Pty Ltd, Australia).Fluorescently conjugated antibodies for cell surface staining (S1 Table ) were prepared in Brilliant Stain Buffer (Cat.no.556349, BD Biosciences, Australia) and incubated with cells in a 96 well plate (Cat.no.COR3894, Corning, In Vitro Technologies Pty Ltd, Australia) for 30 minutes at room temperature.Cells were washed using surface staining buffer, permeabilised using a 10-minute incubation in Perm/Wash Buffer (P/W, cat.no.554723, BD Biosciences, Australia) and washed once with P/W.The remaining antibodies listed in S1 Table were prepared in Brilliant Stain Buffer and incubated with cells in a 96 well plate for 45 minutes at 4˚C.Cells were washed twice with P/W and flow cytometry data was acquired using the Cytek Aurora Spectral Flow Cytometer (5 laser, 64 detector configuration).Data was unmixed using SpectroFlo (Version 3.0.3)and analysed using FlowJo (FlowJo 10.8.2) and GraphPad Prism 8.0.1 (see methods section Flow data analysis, S4 Fig) .Cell single-stained controls were prepared using the same protocol as experimental samples.

Mycobacterial cultures and neutrophil in vitro infection
Mtb single cells stocks for infection were prepared as previously described [62].Mtb H37Rv was grown in a liquid culture of Middlebrook 7H9 medium (Difco, Becton Dickinson, USA) with albumin-dextrose-catalase (ADC, Becton Dickinson, USA) and 0.05% Tween-80 (Sigma-Aldrich, USA) at 37˚C as a standing culture in a tray and mixed by swirling every few days to disperse clumps for 10 days.After 10 days, a liquid culture of 7H9 ADC, without Tween-80, was inoculated with 1/100 th volume of day 10 end-exponential growing phase Mtb and incubated in static standing culture, only swirled periodically, at 37˚C for 10 days.After this, cultures were centrifuged for 5 minutes at 2500RPM.Glass beads were used to break the pellet after which it was resuspended in PBS.The upper part of the bacterial suspension was harvested and spun for 10 minutes at 1400RPM.Then the upper part of the bacterial suspension was harvested again and mixed with glycerol (5% final volume) and aliquots were stored at -80˚C.Aliquots was serial diluted before and after freezing and plated on Middlebrook 7H10 agar (Becton Dickinson, USA) plates with oleic acid-albumin-dextrose-catalase (OADC, Becton Dickinson, USA) for colony forming unit (CFU) determination.Prior to infection, aliquots were thawed at room temperature.
Neutrophils were infected at a multiplicity of infection (MOI) of 1:1 for 1h, and 6h at 37˚C under 5% CO 2 .For the 3 wells in row 1 of each of the two plates, 1ml was removed from each well (3ml per plate giving a total of 6ml from the two plates).After this 2 x 10 7 Mtb H37Rv was added to the 6ml, carefully pipetted to mix and then 1ml (0.33 x 10 7 Mtb H37Rv) was returned to each of these wells for infection.For each infection experiment, Mtb infection MOI were confirmed by CFU counts of serially diluted inoculum plated on 7H10 OADC agar.

Staining for microscopy
Of the neutrophils isolated from 17 older HITTIN and 11 HIT participants, neutrophils from 12 HITTIN and 10 HIT were plated and stained for microscopy imaging.Cell staining was performed on 2 x 10 5 neutrophils per well for uninfected vs infected with Mtb H37Rv for each timepoint 1 and 6h.For each time point (1h, 6h) two wells (A, B) were plated for neutrophils to be infected with Mtb and two wells for no infection (C,D).At the designated timepoints the supernatant was discarded and replaced with 4% PFA (cat.no.43368, Alfa Aesar, USA).The plates were incubated at 4˚C for 24h, before removing the plates from the BSL-3.PFA was replaced with PBS and plates stored at 4˚C.
For immunofluorescence labelling, wells were stained with 1˚[1:500 mouse mAB PL2/3, kindly gifted by Arturo Zychlinsky (nucleosomal complex of Histone 2A, Histone 2B and chromatin)] and 2˚antibody cocktail [1:1000 α-mouse-Cy3 (cat.no.715-166-150, Amersham, UK) and 10μg/ml Hoechst 33342 (cat.no.14533, Sigma-Aldrich, USA)] [72][73][74].Prior to staining PBS was removed from all wells and replaced with Perm/Quench Buffer (50mM NH4Cl, 0.2% saponin in PBS) for a 15-minute incubation period.After removing Perm/Quench, PGAS Buffer (0.2% bovine serum albumin (BSA), 0.02% saponin and 0.02% Azide (NaN 3 ) was added for a 5-minute incubation.Once the PGAS was removed the 1˚antibody cocktail was added and left to incubate overnight at 4˚C.The following day the 1˚antibody cocktail was discarded, and two washes were completed with PGAS buffer before adding the 2˚antibody cocktail and incubating in the dark for 1h at room temperature.The staining was completed by removing the 2˚antibody cocktail and completing three washes with PGAS buffer.Each well was filled with PGAS buffer to the brim.The plates were covered with foil and stored at 4˚C for imaging.

Image acquisition
Image tile scans were acquired with the Zeiss AxioObserver Z1 microscope, equipped with a Colibri 7 light source for excitation of Hoechst 33342 (cat.no.14533, Sigma-Alrich) with LEDmodule 385nm, and Cy3 with LED-module 511nm.A quadruple band pass filter and triple band pass filter were used respectively for detection of Hoechst (wavelength range 412-438nm) and Cy3 (wavelength range 546-564 nm).Images were acquired with a LD A-Plan 40x/0.55 objective as a 6 x 8 tile scan to acquire a total area of 1.21mmx1.21mm(S1 Fig).

Imaging processing
Tiles were stitched together into single images using Zeiss ZenPro software (version 2.6), which were imported into FIJI/ImageJ (version 1.53t).Image tiles were split into separate channels, after which background subtraction (ranging from 50-100 pixels) and Otsu based thresholding was conducted on a per image basis to generate a binarized map of each image.Binary images were further processed by using the binary closing function in FIJI, as well as a top hat filter of 1-2 pixels to remove small non-specific pixels.For the nuclear quantification, an additional watershed function was applied to separate borders more accurately.Thereafter, morphometric data were obtained by specifying the area and circularity ranges of particles to be analysed through the Analyze Particles function in FIJI to determine the total area covered.

Preparation of RNASeq libraries
Uninfected and infected neutrophils (1x10 7 cells for each) were lysed with TRIzol (Invitrogen TRIzol Reagent, Fisher Scientific, USA) after 1h and 6h after infection, and stored at -80˚C.The miRNeasy kit (Qiagen, Germany) was used for total RNA extraction.One sample per participant and condition was used.RNA integrity (RIN) was assessed with the Agilent 2100 Bioanalyzer (Agilent Technologies, Germany).Samples with RIN >7 were selected for library preparation using TruSeq RNA Library Preparation Kit v2, Set A (Illumina, USA).Samples were sequenced in two batches (S4 Table ).Batch 1 was a preliminary test batch for exploratory analysis and to confirm data quality from PMN, consisting of samples from 3 HITTIN (18%) and 4 HIT (36%) participants.This was sequenced as unstranded, 100bp, single-end (SE) on an Illumina HiSeq4000 sequencer at Genome Quebec, Montreal, Canada.Batch 2 with samples from 14 HITTIN (82%) and 7 HIT (64%) participants was sequenced as unstranded, 150bp paired-end (PE) on an Illumina NovaSeq6000 sequencer at Genome Quebec, Montreal, Canada.There was no significant difference in the distribution of the samples used for each batch (p = 0.4, Fisher's exact test) (S4 Table ).

Quality control and raw data pre-processing
The quality of the raw sequence data was accessed by FastQC (v0.11.5) and MultiQC [75,76].The mean quality of reads was high with the mean sequence quality score (Phred Score) >35 for both batches.Duplicates were observed in both batches.For batch 1, sequenced as singleend reads (SE), we detected a fraction of 0.41-0.63duplicates which fall within and below the expected range (0.66-0.74) of duplicates according to the Universal Human Reference RNA (UHRR) [77].The fraction of duplicates (0.4-0.61) from batch 2 was much higher than expected (0.087-0.18) for paired-end (PE) reads [77].With single-end reads, fragmentation bias is usually the cause of these duplicates.For both SE and PE, de-duplication is not recommended [77].To minimize bias, duplicates were not removed from either of the batches.
Batch 1 and 2 were combined after read counts were generated.For read count generation the same method was applied to both batches.After initial raw sequence quality check by FastQC, data was filtered and trimmed using HTStream [75,78].The occurrences of rRNA read contamination was counted but not removed and contamination with PhiX, a control in Illumina runs, was removed by hts_SeqScreener [78].Adapters were trimmed with hts_Adap-terTrimmer and poly(A) tails were removed with hts_PolyATTrim [78].Any remaining N characters (unassigned bases) were removed with hts_NTrimmer and hts_QWindowTrim for quality trimming the end of the reads [78].Reads less than seventy-five base pairs long were removed by hts_LengthFilter [78].
The genome was indexed using GRCh38.p13v34 (ENSEMBLE v100) and reads were aligned to the genome with STAR (v2.5.3a) in a bam format [79][80][81].More than 80% of the reads were uniquely mapped.The reads per gene output files were combined before the final unstranded read counts matrix was extracted for all.The gene read count matrix was input in R (v4.0.3) for further analysis [82].The untransformed and raw count matrix was adjusted for batch effect using ComBat-seq while preserving the signal of the biological variables of interest namely the phenogroups (HITTIN vs HIT), the timepoints and infection [83] (S5 Fig).

Differential gene expression analysis
Raw counts were transformed to counts per million (CPM) and filtered using "filterByExpr" in the edgeR package [84].Briefly, the function kept genes at least 25 read counts or more in at least a minimum of samples (calculated as 70% of the samples in the smallest group).For this dataset genes with a CPM of 1.32 in at least 8 (70% of the smallest group of 11 samples) are retained, leaving 14602/60622 (24%) of the genes for further differential expression analysis (S6 Fig) .Normalization scaling factors were generated by "calcNormfactors" in edgeR using the method of trimmed mean of M-values (TMM) [84].
Outlier samples were observed in the multidimensional scaling (MDS) plots (S5 Fig).Instead of removing the outliers, limma's function "voomWithQualityWeights" allows for variations in samples by taking sample-specific variablity as well as the 'global intensity-dependent variability trends' as accounted for by 'voom', into consideration [85].The expression matrix was normalized and transformed to log 2 CPM and sample specific weights were incorporated with abundance dependent weights using limma's (v3.46.0) function "voom-WithQualityWeights" [85,86] (S7 Fig).

Data exploration
Data was visualized with multidimensional scaling (MDS) plots and principal components analysis (PCA) Scree plot (S8 and S9 Figs).Normalized counts were plotted without covariate correction to examine data for potential covariate effect.Initially clear separation was seen by batch effect, but this was corrected for (S5

Model design
For the analysis we used a paired design and created the model with subject IDs (for blocking in a paired design), group (HITTIN and HIT), time of infection (1 and 6h) and Mtb infection status (uninfected and infected) as factors in the model using model.matrixfrom the stats package (v4.0.3).
In the individual blocking, individuals were defined per timepoint, due to large variance introduced to each individual due to time effect (S8 Fig) .The three factors group, time of infection and Mtb infection status were grouped together as a single interaction term and was modelled with the subject effect as a means model: Where Gi is the log 2 (CPM) expression for each gene i (n = 14602) and E(Gi) is the expected gene expression.β 0 represents the intercept.βjkl represents the mean expression for each individual j (j = 1 to total of n) at time k (k = 1h or 6h) for Mtb infection state l (l = uninfected or infected).H represents HIT and N represents HITTIN.xj is a variable representing the n samples in the data.The β value is the mean expression for each specified group after blocking on the individual averages.� is the residual term and is assumed to be normally distributed with a constant variance across the range of data.
Contrasts were made using makeContrasts and defined as: iii) Differential response between groups at 1h infected ðb HITTIN1inf Þ À ðb HIT1inf Þ iv) Differential response between groups at 6h infected ðb HITTIN6inf Þ À ðb HIT6inf Þ Linear models were fitted with lmFit and eBayes to calculate the gene-wise test statistics (moderated t-statistic, p-values and B-statistic).Criteria for DEGs was defined by a Benjamini Hochberg FDR procedure as genes with an absolute logFC � 0.2 and adjusted p value � 0.05.
DEGs were used in a gene set enrichment analysis for pathways and Gene ontology (GO) enrichment using ReactomePA v1.34.0 and clusterProfiler v3.18.1 [87].enrichGO was used to test GO biological process, enrichKEGG for KEGG pathways and Reactome with enrichPathway.An FDR � 0.1 was used for Benjamini-Hochberg's multiple testing correction.Pathways and GO terms with less than five assigned genes were excluded.

Fluorescence image analysis
Image from total of 9 HITTIN and 5 HIT neutrophil were included in the final analysis.Results of the Hoechst cell nuclei area count as well as the total NETs area (chromatin channel) were analysed in R studio.The log ratio of total nuclei area (Hoechst channel) at 6h to 1h was used to compare the infection effect of Mtb between and within groups on the change in area of cells (as a proxy for cell counts, since masking often under or overestimated cell counts, especially with neutrophil nuclei with trinucleate structures and cell clumps) over 1h to 6h (Figs 4i and S1).The total NET area change was determined by the total NET area (chromatin channel) to the total cell area (Hoechst channel) as determined for each timepoint and noninfection or infection with Mtb.To investigate if there is a difference in NET area change from 1h to 6h between PMN HITTIN and PMN HIT in response to infection, we calculated the log transformed ratio of the NET area change at 6h to 1h in uninfected and Mtb infected for each subject (Figs 4ii and S1).
A two-way mixed ANOVA analysis approach was used with the phenotype group (HITTIN vs HIT) describing the between-subject factor and Mtb infection status (non-infected vs infected) the within subject factor.Assumptions of normality, homogeneity of variances and homogeneity of covariances were met.After the removal of the initial extreme outlier, there were still some 3 outliers in the remaining analysis.These outliers were not removed, and a robust ANOVA was performed using the WRS2 package in R [88].Further pairwise comparisons were done between PMN HITTIN and PMN HIT for non-infection as well as infection with Mtb, using a pairwise t-test and Bonferroni adjusted p-values."HITTIN 6h infected" is the same but for 6h."HIT 1h infected" represents the log 2 FC gene expression effect when comparing the 1h Mtb infection to the 1h uninfected response in PMN from HIV-1-positive, IGRA positive, tuberculin positive individuals (PMN HIT )(y-axis) ."HIT 6h infected" is the same but for 6h.Scatterplots display the correlation of the log 2 FC gene expression as defined for each group at 1h (A and B) or 6h (C and D).Each grey dot represents a differentially expressed gene (DEGs).The grey dots show combined DEGs from both groups represented on the x and y axis.Blue coloured dots indicate DEGs from the x-axis group and red dots indicate the y-axis group.Pearson correlation (R value) was calculated for the coloured DEGs.The log 2 FC for each group is plotted with the y-axis as the reference group.Differential gene expression is more pronounced in genes with the highest log 2 FC, with the majority of genes displaying good correlation between phenotype groups.

Fig 1 .
Fig 1. Volcano plots of differential gene expression at 6h Mtb infection by PMN from HITTIN and HIT.Volcano plot for transcriptional responses to Mtb challenge for neutrophils from HITTIN (PMN HITTIN ) and HIT (PMN HIT ) participants at 6h (A-C) post Mtb infection.The y-axis shows the negative log10 unadjusted P value and the x-axis the log 2 fold change (FC).The vertical dashed lines represent log 2 FC thresholds of -0.2 and 0.2.Each gene is represented by a dot.Genes which are downregulated or upregulated as determined by the FDR � 5% are shown in blue and red, respectively.Genes with non-significant expression changes and below the log 2FC threshold are shown in grey.Differentially expressed genes at 6h post-infection compared to 6h uninfected PMN from HITTIN (PMN HITTIN ) (A) and HIT participants (PMN HIT ) (B).Significant differentially triggered genes between PMN HITTIN and PMN HIT at 6h post infection (C).https://doi.org/10.1371/journal.pgen.1010888.g001

Fig 2 .
Fig 2. The absolute log fold change Mtb infection effect at 1 and 6 h for PMN HITTIN and PMN HIT .The density plot shows the absolute log fold change (logFC) Mtb infection effect of each group at 1h and 6h.All DEGs from Fig 1 had their logFC converted to absolute values and plotted using density function.Absolute values for neutrophil DEGs from HITTIN participants are shown in light blue at 1h (median = 0.051) and a darker blue at 6h (median = 0.383).Absolute values for neutrophil DEGs from HIT participants are shown in magenta at 1h (median = 0.065) and red at 6h (median = 0.550).Vertical coloured lines indicate the median of absolute values.PMN HITTIN show a significantly lower logFC response to Mtb infection at 6h (p < 2.2e-16, Wilcoxon rank sum test) compared to 1h of infection with Mtb (p = 0.8255, Wilcoxon rank sum test).https://doi.org/10.1371/journal.pgen.1010888.g002

Fig 3 .
Fig 3. Manhattan plot for enrichment tests of GO terms and Kegg and Reactome pathways.The Manhattan plot shows pathways and GO terms for the DEGs triggered significantly differentially by Mtb in neutrophils from HITTIN versus neutrophils from HIT participants after 6h of infection.The tested terms are distributed along the x-axis.The y-axis represents the negative log 10 false discovery rate (FDR) result with the horizontal dotted line indicating the 10% FDR significance cut-off.Grey dots represent terms not meeting the significance threshold.Coloured dots represent significant terms from GO (red), Kegg (green) and Reactome (blue).The number of DEGs in each term is represented by the scaled size of the dot.Panel (A) shows the terms detected by all significantly different DEGs in the HITTIN vs HIT contrast.Panel (B) shows significant terms which are driven by genes with significant positive FC/less downregulated by PMN from HITTIN participants while panel (C) displays terms for genes significantly less upregulated for PMN from HITTIN participants as compared to PMN HIT after 6h of Mtb infection.https://doi.org/10.1371/journal.pgen.1010888.g003

Fig 4 .
Fig 4. Representative images and quantification of the change in NET and nuclei area following Mtb infection of neutrophils from HITTIN and HIT.PMN HTTIN (A-H) and PMN HIT (I-P) following Mtb infection for 1h (A-D and I-L) and 6h (E-H and M-P).Mtb infected neutrophils at 1h and 6h were stained with Hoechst 33342 (A, E, I, M) and PL2-3 (B, F, J, N).Overlap between the two stains is shown in (C, G, K, O) with enlarged box panels in (D, H, L, P; blue, DNA; yellow, chromatin); (i) Boxplots for total log change in nuclei area (Hoescht stain, calculated as shown in S1 Fig) from 1h to 6h uninfected and infected with Mtb in HITTIN vs HIT.There was a significant Mtb infection effect, F(1,12) = 9.729, p = 0.009 with pairwise comparisons showing a significant difference between the total change in cell nuclei area from 1h to 6h after Mtb infection in HITTIN compared to HIT (p = 0.04*, pairwise t-test), (ii) Boxplots for total log change in NET area (corrected for by the change in cell nuclei area, calculated as shown in S1 Fig) (PL2-3 anti-chromatin stain) from 1h to 6h uninfected and infected with Mtb in HITTIN vs HIT.There was a significant Mtb infection effect, F(1, 10.5924) = 5.3398, p < 0.0421 with pairwise comparisons showing a significantly difference between the total changed in cell nuclei area from 1h to 6h after Mtb infection in HITTIN compared to HIT (p = 0.0003***, pairwise t-test).The scale bars represent 100 μm (A-C, E-G, I-K, M-O) and 50 μm (D, H, L, P). https://doi.org/10.1371/journal.pgen.1010888.g004 Fig).The post batch corrected data was reviewed (S10-S17 Figs).Sex and smoking showed some separation and would need further investigation (S10 and S11 Figs).No separation could be seen for age, BMI, chronic disease history, previous INH used, duration of INH and alcohol use (S12-S17 Figs).Using the paired design with blocking, accounted for the covariates effects by blocking on the effect of each individual and are discussed below.
i) Group specific Mtb infection compared to no infection effect at 1h (b HITTIN1inf ; b HIT1inf ) ii) Group specific Mtb infection compared to no infection effect at 6h (b HITTIN6inf ; b HIT6inf ) Mtb infection responses after 6h.(XLSX) S2 File.Table showing cell population distribution data as determined by flow cytometry.(XLS) S3 File.Table showing morphometric data from microscopy.(XLS) S1 Fig. Volcano plots of differential gene expression at 1h infection by PMN from HITTIN and HIT.Volcano plot for transcriptional responses to Mtb challenge for neutrophils from HITTIN (PMN HITTIN ) and HIT (PMN HIT ) participants at 1h (A-C) post Mtb infection.The yaxis shows the negative log10 unadjusted P value and the x-axis the log 2 fold change (FC).The vertical dashed lines represent log 2 FC thresholds of -0.2 and 0.2.Each gene is represented by a dot.Genes which are downregulated or upregulated as determined by the FDR � 5% are shown in blue and red, respectively.Genes with non-significant expression changes and below the log 2 FC threshold are shown in grey.Differentially expressed genes at 1h post-infection compared to 1h uninfected PMN from HITTIN (PMN HITTIN ) (A) and HIT participants (PMN HIT ) (B).Significant differentially triggered genes between PMN HITTIN and PMN HIT at 1h post infection (C).(PDF) S2 Fig. Correlation plots of log 2 FC gene expression after Mtb infection of PMN from HIT-TIN and HIT."HITTIN 1h infected" represents the log 2 FC gene expression effect when comparing the 1h Mtb infection to the 1h uninfected response in neutrophils (PMN) from HIV-1-positive persistently TB, tuberculin and IGRA negative individuals (PMN HITTIN )(x-axis) .
imaging process overview to calculate neutrophil nuclei and neutrophil extracellular trap (chromatin) area.Scale bars: 100 μm (A-C) and 50 μm (D-I).Isolated neutrophils were stimulated with Mtb and stained with Hoechst 33342 (A and boxed area D) and PL2-3 (B and boxed area E).C and the boxed area F show the overlap (blue, DNA; yellow, chromatin).The segmentation of fluorescent signals used to calculate the neutrophil nuclei area (G) and neutrophil extracellular trap area (H) with the overlap shown in (I).The scale bars represent 100 μm (A-C) and 50 μm (D-I).(PDF) S4 Fig. Flow Cytometry Analysis of Cell Populations.The first gate was applied to exclude debris (the low SSC-H, FSC-H values in the bottom left corner).Single cells were separated and gated by a CD45+ marker for leukocytes.After single cell gating, CD45+ cells were then grouped into CD15+CD66b+ (granulocytes) and CD15-CD66b-(non-granulocytes) cells.The CD15+CD66b+ cells were further classified as CD14-CD16+ (Neutrophils) and CD14-CD16-(Eosinophils). CD15-CD66b-were stratified as CD3+ (T-cells), CD3-CD14+ (Monocytes) and CD14-CD3-(Other) cells.(PDF) S5 Fig. Multidimensional scaling (MDS) plot before (A) and after (B) ComBat-seq batch correction.The multidimensional scaling (MDS) plots the Euclidian distances between samples with the x and y axis representing the sample distances between samples of read counts normalized by depth but not covariates.Each row of plots in (A) and (B) represent dimension 1 to 5 respectively (represented by the x-axis) and shown with the combination of the other dimensions on the y-axis.Samples are colored for batch effect with green showing samples sequenced on Illumina HiSeq2500 and blue for samples sequenced on Illumina NovaSeq6000.(A) shows the multiple combination of plots by dimension for the normalized counts prior to batch correction.Clear separation of samples by batch can be seen and is most notable in dimensions 4 and 5. (B) Appropriate batch correction after correction of raw counts with ComBat-seq.After batch correct outlier samples are most notably observed in dimension 3. (PDF) S6 Fig. Density plot showing successful filtering applied to counts.Density plot with the density represented on the y-axis for the log-cpm values before (A) and after (B) filtering.The dotted vertical line is equivalent to the counts per million (CPM) threshold of 1.32 which was used in the filtering step.(PDF) S7 Fig. Voom observational mean-variance trend and sample-specific weights.(A) shows rescaled (square-root of standard deviations) residual variances plotted against the mean expression (log 2 transformed with an offset of 2) of each gene.A decreasing trend is seen between the mean gene expression and the variance with higher expressed genes showing less variation.The sample specific weight for each sample used in the analysis is shown in (B).A total of 112 samples were analysed (28 participants and 4 conditions, uninfected and infected after 1 and 6h).(PDF) S8 Fig. Multidimensional scaling (MDS) plot of group separation.The multidimensional scaling (MDS) plots the Euclidian distances between samples with the x and y axis representing the sample distances between samples of read counts normalized by depth but not covariates.Each row of plots in represent dimension 1 to 5 respectively (represented by the x-axis) and shown with the combination of the other dimensions on the y-axis.Samples are colored for phenotype (HITTIN and HIT), timepoint (1h and 6h) and infection status (uninfected [NEG] and infected [INF]) and as depicted in the legend.Separation of the groups by time can be seen in dimension 1, by infection in dimension 2 and the phenotype in dimension 3. Additional separation is seen in dimension 4 and 5 due to sex and possibly smoking (see S8 and S9 Figs).(PDF) S9 Fig. Scree plot.The top 5 principal components are represented on the x-axis and the variance explained by each component on the y-axis.Principal component (PC) 1 contributes to

Table ) .
A higher number of Mtb activated gene expression changes were observed at 6h post infection compared to 6h uninfected with 3106 up-and 3548 downregulated DEGs for HITTIN and 3816 up-and 3794 downregulated DEGs for HIT (Fig 1A and 1B, S3 Table).

) (Fig 2). Despite overall lower log 2 FC among HITTIN, there was a strong correlation of log 2 FC values for each gene between PMN HITTIN
and PMN HIT , suggesting that Mtb responses globally are conserved across phenotypes (S2 Fig).Irrespective of up-or downregulation of specific genes, the absolute response to Mtb after 6h was always smaller in PMN HITTIN .