BACKGROUND
the mycobacterium avium complex comprises a group of environmental and opportunistic intracellular pathogens
<cit>  with the ability to persist within macrophages and escaping the host’s killing mechanisms
<cit> .
the complex is constituted by the subspecies of m. avium together with m. intracellulare <cit> .
two subspecies of clinical importance for humans and animals are m. avium subsp.
hominissuis  and m. avium subsp.
avium
<cit> .
mah is the subspecies mainly isolated from porcine and human cases of mac infection
<cit> , while maa is the causative agent of avian tuberculosis, but is occasionally responsible for m. avium infections in mammals
<cit> .
little is known about the routes of transmission and the virulence mechanisms of m. avium, and a zoonotic potential has not been ruled out
<cit> .
definite reasons for mah being the dominating subspecies in humans and pigs infected with m. avium have not been identified, although differences in exposure could be considered a contributing factor.
explanations to the higher prevalence of mah in the human and porcine population might also be found in differences between mah and maa at the molecular level of interaction between the host and mycobacterium, a process crucial to the outcome of infection, which therefore might affect the prevalence of the subspecies.
however, little has been done to investigate whether mah and maa interact differently with their target cells.
mycobacterial virulence depends on the ability of these bacteria to invade, persist and replicate within the hostile environment in macrophages of the host organism
<cit> .
this is achieved by the protective properties of the mycobacterial cell wall and by interference with the host’s immune response
<cit> .
to assess the ability of mah and maa to grow intracellularly, the use of primary cells will give a more representative result compared to immortalised macrophage cell-lines.
tissue resident macrophages represent a heterogeneous population and it is not clear whether m. avium infect a particular phenotype.
monocyte derived macrophages matured in absence of any stimuli give rise to a mixed culture of both the m <dig> and m <dig> phenotype
<cit>  and such cultures may thus be considered appropriate for infection studies with m. avium.
finally, using macrophages derived from the relevant target species is important since the global inflammatory response evoked in mice recently was demonstrated to be vastly different from humans
<cit> .
to the authors’ best knowledge, comparative studies of the ability of maa and mah to invade, replicate and trigger immune responses in human primary macrophages have not been published.
the aim of the present study was therefore to assess the growth of m. avium isolates in human primary macrophages and to characterise the response triggered in infected macrophages by selected isolates.
the isolates used in the present study represented the two different subspecies, were derived from different host species and differed in biofilm abilities, glycopeptidolipid  genes, and in the presence of the insertion sequence  element ismpa <dig> <cit> .
the present data showed that maa replicated to a larger extent in human macrophages suggesting that an inability to grow intracellularly is not likely the cause of the discrepancy of the prevalence of the two subspecies in humans.
furthermore, the response induced in the macrophages could not entirely explain the difference in growth rate.
methods
mycobacterial isolates
four clinical isolates of m. avium subsp.
hominissuis  and two clinical isolates of m. avium subsp.
avium  were used in the present study.
details on origin, biofilm producing abilities, presence of ismpa <dig>  nsgpl- and ssgpl genes were described by johansen et al.
<cit> , and are summarised in table
<dig>  the human isolates were from a strain collection at the the norwegian institute of public health.
all isolates were grown in middlebrook 7h <dig>  w/oadc , 1% tween  <dig> and 10% glycerol under agitation for two weeks before aliquots were stored at −70°c in an equivalent medium containing 10% glycerol.
the concentration of bacteria in the stock suspensions was determined by colony forming unit  counts.
bacterial inocula for the gene expression assay were prepared from seven days subcultures grown on plates of middlebrook 7h <dig>  with 10% oadc  at 37°c.
to adjust the number of bacteria in the inoculum, hsp <dig> real-time pcr
<cit>  was performed on serial dilutions in hbss  of a suspension adjusted to mcfarland standard  <dig> , homogenised through a 23 g needle to minimise clumping.
mpa1
isolates of mycobacterium avium used in the present study show different combinations of subspecies designation, origin, biofilm abilities and the genomic presence of ismpa <dig> and genes involved in the synthesis of non-specific  or serovar specific  gpls  <cit> .
cell culture
peripheral blood mononuclear cells  were isolated from buffy coat drawn from healthy human blood donors.
the material was commercially obtained from a blood donation centre, and the study was approved by the regional committee for medical research ethics, south-east norway.
fifty ml buffy coat was diluted 1: <dig> with rpmi  <dig>  with 1% l-glutamine .
pbmcs were harvested by density gradient centrifugation with lymphoprep  and resuspended in pbs with  <dig> % bsa and 2 mm edta at ph  <dig> .
isolation of cd14+ cells was performed using cd <dig> macs magnetic beads  according to the manufacturer’s instructions.
the cells were resuspended in rpmi  <dig> with 1% l-glutamine supplemented with 10% human serum to the desired concentration of the respective infection assay, added to the appropriate costar® cell culture plates  and incubated at 37°c with 5% co <dig>  the cd14+ cells were incubated overnight before used in the uptake and replication assays and allowed to mature into macrophages by incubation for five days before used in the gene expression study.
uptake and intracellular replication
uptake and intracellular replication of the six m. avium isolates in human blood derived cd14+ cells were assessed.
the assays were performed in 96-well cell culture plates , which in our hands gave a stable monolayer for most donors.
each isolate was tested in triplicates at an moi of 5: <dig>  the cells were washed twice with pre-warmed medium to remove extracellular bacteria before the quantification of intracellular bacteria was performed.
intracellular uptake of the bacterial isolates was assayed after three, six,  <dig> and  <dig> hours in cells from three donors.
intracellular replication of the same isolates was assayed in six other donors, following the protocol described above.
based on the results from the uptake studies, cd14+ cells were incubated with the bacterial isolates for six hours.
after removal of extracellular bacteria, fresh medium was added and the cells were incubated further.
cell lysates were harvested at six hours, and after one, four and seven days of infection.
quantitation of bacteria and cells
numbers of intracellular bacteria and cells was determined by enumeration of mycobacterial and human genomes in cell lysates.
replication of bacterial isolates in the supernatant had been tested and found insignificant.
lysis of cells was achieved through adding milliq water to the wells.
after  <dig> minutes, lysates were transferred into 2 ml o-ring vials  containing 200 μl silica beads  and stored at −20°c.
lysates from triplicate wells were pooled.
after completing all sampling, lysates were thawed and inactivated at 100°c for 30 min, prior to mechanical disruption using the minibeadbeater  for 45 sec.
enumeration of mycobacteria and cells was achieved by performing a duplex real time pcr of a 103 bp long segment of the mycobacterial hsp <dig> gene and a 150 bp long segment of the human β-globin gene on the lysates, as described by salte et al.
<cit> .
briefly, the segment of the mycobacterial gene hsp <dig>  was amplified using the primers mycofp <dig>  and tb  <dig> , and the vic tagged taqman-mgb probe mycopr <dig>
<cit> .
the segment of the human β-globin gene was amplified by using the primers bg-f  and bg-r , and the fam tagged taqman-mgb probe bg-taq
<cit> .
dilutions of known concentrations of commercially available mycobacterial dna  and human dna  were run together with the samples and used to create standard curves.
the conversion of the known amount of standard dna to number of bacterial and human genomes were described by salte et al.
<cit> , as well as the validation of the number of mycobacterial targets against cfu counts.
total volume of each reaction mixture was 20 μl, consisting of 8 μl lysate, 10 μl quantitect® mastermix  and 1 μl of a primer/probe mix for each target, with a stock concentration of 8 μm for each primer and 4 μm for the probe, resulting in a final primer- and probe concentration of  <dig>  μm and  <dig>  μm, respectively.
samples and standard dna dilutions were run on a stratagene mx3005p mxpro v  <dig>   and data analysed by the stratagene mxpro v  <dig>  software  with the following thermal cycling conditions: 95°c for 15 min, followed by  <dig> cycles of 95°c for 1 min and 60°c for 1 min.
gene expression profiling
isolates mah vi <dig>  mah  <dig> and maa  <dig>  were used to infect macrophages for gene expression analysis.
cd14+ cells were obtained as described above and allowed to mature into macrophages by incubation for five days in duplicate wells .
the cells were then infected with the three isolates at a moi of 10: <dig> for four hours.
total rna was isolated of the adherent macrophages using the rneasy mini kit® , following the manufacturer’s instructions.
rna from the different donors within each assay was pooled per bacterial isolate.
the first assay was performed in macrophages from three donors and subsequently repeated twice in two new donors per assay.
the quality and integrity of rna was confirmed by agilent  <dig> bioanalyzer  and nanodrop™ .
gene expression levels were profiled using illumina humanht- <dig> expression bead chip .
the gene expression data was exported from illumina’s genomestudio software package without background subtraction or normalization.
the raw data was preprocessed in the r language and environment for statistical computing
<cit>  employing the normexp by control  algorithm
<cit> .
probes that were not significantly expressed above background  were filtered out of the analysis.
differentially expressed genes were identified employing bayesian/moderated t-statistics with false discovery rate  control for multiple testing
<cit> .
differentially expressed genes were identified at fdr <  <dig> .
the raw microarray data have been made publicly available in the array express repository .
bioinformatics
biological functions enriched in the differentially expressed genes were identified using modular enrichment analysis with david
<cit> .
enrichment scores >  <dig>  were considered statistically significant.
biological pathways enriched in the data were identified using ipa software .
upstream regulator analysis was performed in ipa software to infer the putative pathways that give rise to the observed gene expression changes.
molecular interaction networks were reconstructed from the differentially expressed genes in ipa software using mechanistic data from prior studies
<cit> .
in this analysis differentially expressed genes were defined at fdr <  <dig> .
protein analysis
cd14+ cells were isolated and matured as described above for gene expression studies.
cell lawns were infected at an moi of 5: <dig> for  <dig> hours with maa  <dig>  mah vi <dig> and mah  <dig>  supernatants were removed for cytokine profiling, and subsequently the adherent cells were washed twice in pbs and collected with a cell scraper in 100 μl of 1x sds sample buffer for immunoblotting.
all samples were stored at −70°c until further processing.
sds-page and immunoblotting
proteins were separated by horizontal sds-page with precast gradient excelgel™ sds gradient  <dig> – 18%  and excelgel™ sds buffer strips  in a multiphor™ ii  <dig> , prior to diffusion blotting onto nitrocellulose membranes with pore size  <dig>  μm
<cit> .
proteins remaining within the gel were visualized by coomassie brilliant blue staining.
membranes were blocked with pbs with 2% bsa and incubated over night with primary antibodies  recognizing the following selected proteins involved in apoptosis: cleaved caspase- <dig> , cleaved caspase- <dig> , cleaved caspase- <dig> , cleaved caspase- <dig> , bcl- <dig>  phospho-bcl- <dig> .
detection of bound antibodies was performed by incubaction with ecl™ anti-rabbit igg, horseradish peroxidise labelled linked f <dig> fragment from donkey  and visualization by immun-star™ westernc™ kit .
cytokine profiling
the concentration of tnf-α, il- <dig>  il- <dig>  il- <dig> and il- <dig> was measured in supernatants by using bio-plex™ cytokine assay  within the bio-plex suspension array system , following the manufacturer’s instructions.
differences between isolates were tested for significance using the wilcoxon matched pairs signed-rank test, considering p-values ≤  <dig>  significant.
RESULTS
uptake and intracellular replication
all the isolates were taken up by the cd14+ cells, and intracellular bacteria were detected after  <dig> hours .
maa  <dig>  mah h <dig> and maa  <dig> appeared to be taken up to a greater extent, while the mah isolate vi <dig> showed the lowest degree of uptake in macrophages.
between six and  <dig> hours of infection some of the isolates showed an exponential increase in the number of intracellular bacteria, suggesting the event of onset of intracellular replication.
six hours was therefore chosen as the time between inoculation and removal of extracellular bacteria for the subsequent study of intracellular replication over time.
at this point in time vi <dig> had significantly  lower numbers of intracellular bacteria compared to the other isolates except for mah h <dig>  while no significant difference were observed between the other isolates using the paired t-test.
intracellular growth was only observed for the two maa isolates.
there was a significant increase in the ratio of mycobacteria to human genomes at day four for isolate maa  <dig>  and at day  <dig> for maa  <dig>  and maa  <dig>  and in absolute numbers of bacteria for isolate maa  <dig>  at day  <dig> using the paired t-test .
the remaining isolates, which were all mah, did not show significant increase in intracellular numbers throughout the duration of infection, however the count of intracellular bacteria was quite stable, indicating that the bacteria were able to persist within the macrophages.
enumeration of both cells and intracellular mycobacteria was performed by real-time pcr of single-copy genes.
this method does not distinguish between viable and non-viable bacteria and cells, but it has the advantage that the number of bacteria is not influenced by the degree of clumping which varies between maa and mah.
furthermore it enabled the measurement of the number of intracellular bacteria per cell in addition to absolute counts of bacteria.
some degree of cell death and loss of adherence were observed in both infected and non-infected cells, and varied considerably between donors.
similar results were observed whether absolute numbers of bacteria or the bacteria per cell ratio were used, however the difference in cell death between donors may explain why the absolute counts of bacteria only reached a significant level of p <  <dig>  for maa  <dig> at day  <dig>  while the ratio was significant for both maa isolates at day  <dig> and for maa  <dig> at day  <dig>
gene expression analysis and bioinformatics
based on the results from the initial replication assays, three isolates were selected for gene expression studies in macrophages.
maa  <dig>  which replicated intracellularly, and mah vi <dig>  which only persisted within cells, were selected, in addition to mah  <dig>  which, as opposed to the other two isolates, contains the is element ismpa <dig> and belong to an rflp type that has been demonstrated to be particularly virulent in pigs
<cit> .
the illumina® humanht- <dig> v <dig> expression bead allows whole-genome expression analysis of  <dig> samples on the same chip, each array targeting genomewide with more than  <dig>  probes.
of the  <dig>  probe sets with signal intensity above background,  <dig>  probe sets were differentially expressed  in the cells infected with m. avium , regardless of the individual isolate, compared to uninfected control cells .
the largest response was induced by maa  <dig> where  <dig>  probe sets  were differentially expressed.
of these probe sets,  <dig>  were up-regulated, and  <dig>  were down-regulated.
mah vi <dig> modulated expression of  <dig>  probe sets  representing  <dig>  unique genes.
in comparison, mah  <dig> altered the expression of only  <dig> probe sets  representing  <dig> unique genes.
a large overlap was observed in the responses elicited by the three isolates .
there were virtually no genes  that were uniquely modulated by isolate mah  <dig>  in contrast, hundreds of genes were uniquely differentially expressed in response to infection with maa  <dig> and mah vi <dig>
bioinformatics analyses demonstrated that the major biological functions up-regulated in response to the three isolates were apoptosis, inflammatory response, cytokine/chemokine activity, signal transduction, cell-proliferation and regulation of t-cell activation .
biological pathways up-regulated in the responses were related to apoptosis and innate immunity and cytokine signaling .
in addition, communication pathways between innate and adaptive immune responses and il- <dig> signaling were up-regulated.
upstream regulator analysis was performed to infer the putative molecular mechanisms that give rise to the observed gene expression changes, and is illustrated in table
<dig>  this analysis revealed major roles for tnf, nfkb, il-1b, ifng and others in the regulation of the responses.
altogether these analyses revealed a strong overlap between the impacts of all the three isolates.
the analysis was done using the database for annotation, visualization and integrated discovery .
an enrichment score >  <dig>  was considered statistically significant.
pathways induced by infection of human monocyte-derived macrophages with m. avium isolates.
the analysis was done in ingenuity systems software.
the ten top pathways are shown.
genes predicted to act as up-stream regulators for the observed responses after infection of human monocyte-derived macrophages with m. avium isolates.
the analysis was done in ingenuity systems software.
responses to chemicals/drugs are excluded.
the ten top genes are shown.
to identify the putative molecular mechanisms that are regulating the responses, we constructed molecular interaction networks for each isolate using mechanistic findings from previous studies.
this analysis showed that maa  <dig> induced the largest gene network in human macrophages, followed by mah vi <dig>  then by mah  <dig> .
the networks were organized around a series of hyper-connected hubs, and the most important hubs  are listed in the table
<dig>  notably, the network connectivity patterns were dominated by tnf, il-1b, il- <dig>  nfkb <dig>  nfkbia, ptgs <dig> and others, and many of these pathways were also highlighted by the upstream regulator analysis.
differential gene expression elicited by maa  <dig>  mah vi <dig> and mah 1655
at fdr <  <dig> , no significant differences were observed between cells infected with maa  <dig> and mah vi <dig>  while only two genes  were upregulated in maa  <dig> compared to mah  <dig> and one gene  was upregulated in mah vi <dig> compared to mah  <dig>
protein analysis
gene expression analysis pointed to apoptosis and immune response as two major biological functions upregulated in response to all three m. avium isolates, albeit at lower levels in mah  <dig> compared to maa  <dig> and mah vi <dig>  we thus attempted to look at proteins involved in apoptosis by western blot and cytokine by a multiplex assay.
however, differences between the isolates with regards to apoptotic proteins were not detected, probably due to the lower sensitivity in this assay compared to the microarray analyses .
macrophages infected with all three isolates consistently produced more il- <dig>  il- <dig> and tnf-α than uninfected controls .
the production of il- <dig> and il- <dig> in uninfected cells was often high, and in some donors infection led to a decrease in cytokine production.
the decrease in cytokine production was particularly evident for il- <dig> after infection with maa  <dig>  to a lesser extent by infection with mah  <dig> and rarely by mah vi <dig> .
there was a large donor to donor variation in the amount of cytokines produced, and this variation surpassed the difference induced by the isolates.
nevertheless, mah vi <dig> consistently induced more il- <dig>  il- <dig> and tnf-α than maa  <dig>  and more il- <dig> than mah  <dig> .
discussion
the present study describes the uptake and intracellular growth of maa and mah isolates in human monocytes, and characterises genome-wide transcriptional responses elicited by these isolates in macrophages.
mycobacterium avium isolates of different genotypic characteristics were all able to enter and persist within human primary monocytes.
intracellular replication was only seen by maa suggesting that the differences in the ability to survive and replicate in macrophages cannot explain why humans are more often infected with mah.
the gene expression program induced by infection with maa and mah isolates was broadly comparable in terms of the biological functions, pathways, and gene networks that were perturbed.
however, the degree of perturbation varied considerably.
the different level of expression led to an apparent discrepancy in the number of unique genes induced by the various isolates and the lack of differences when comparing the isolates with each other.
many of the genes significantly induced by maa <dig> were also induced at a lower level by the other isolates .
this intermediate expression was neither significantly different from uninfected cells nor from maa <dig>  a larger samples size might have given more significant differences between the isolates since large donor to donor variation is present when using primary cell lines.
in general, responses to maa <dig> were excessive, and were characterized by increased expression of genes downstream of type i and ii interferon signaling.
further studies will be required to determine if variations in these pathways account for the exaggerated responses and/or differences in prevalence.
similar pathways were induced in ruminant monocyte derived macrophages infected with the closely related mycobacterium avium subspecies paratuberculosis  <cit> , while a study using the murine cell line raw <dig>  showed that maa and mah induced higher levels of pro-inflammatory genes compared to map <cit> .
together with our data, these studies suggest that similar pathways were induced by all m. avium subspecies.
a direct comparison of the level of macrophage activation is however challenging considering the variation in cell cultures and isolates used in the various studies.
macrophage activation can be interpreted as an appropriate reaction to an invading pathogen or as a failed attempt of the mycobacteria to silence the host’s defences.
maa  <dig> and mah vi <dig> differed completely in their ability to replicate within host cells, although both isolates induced a strong activation of macrophages and innate immune responses on the transcriptional level.
analyses on the protein level revealed that the final synthesis of inflammatory cytokines were higher in response to mah vi <dig> compared to maa  <dig>  one could thus argue that the inflammatory response was sufficient to limit growth of mah vi <dig>  but not maa  <dig>  the complexity of the outcome of an infection, and in vitro culture, was illustrated by the fact that despite a low activation of the innate immune response mah  <dig> did not proliferate intracellularly.
the microarray data provided no clear explanation to the lower level of inflammatory cytokines induced by maa  <dig> on the protein level compared to mah vi <dig>  however, regulation of inflammation on the post-transcriptional level has been demonstrated to be important
<cit> , and a down-regulation of immune responses to avoid immunopathology can occur in chronic infections.
the presence of unidentified virulence factors in maa  <dig> more actively shutting down the innate responses cannot be excluded.
the genomes of the three isolates did for instance differ in the presence of gpl genes, and the involvement of gpls in m. avium colony morphology and of ssgpls in pathogenicity is increasingly recognised
<cit> .
the induction of pro-apoptotic genes, such as ripk <dig>  bid and tbid was seen after infection with all three isolates, however apoptotic pathways were affected to a lesser degree by mah  <dig> than by maa  <dig> and mah vi <dig>  the correlation between apoptosis and virulence of mycobacterial strains is debated.
inhibition of apoptosis has been described as a virulence factor in mycobacteria
<cit> , but it has also been suggested that m. avium can use apoptosis of macrophages as a tool for spreading
<cit> .
one of the other more pronounced differences between the isolates was the increased expression of genes involved in activation of t cells like cd <dig>  cd <dig>  cd <dig> and il <dig> in maa  <dig> and mah vi <dig>  cd <dig>  cd <dig> and cd <dig> are transmembrane surface proteins involved in activation of macrophages and t cells upon binding, while il- <dig> is essential for the persistence and function of th <dig> cells
<cit> .
these findings suggest that differential activation of apoptotic pathways and/or activation of t cells may occur in humans after infection with the different isolates.
the lower activation of macrophages by mah  <dig> may be a useful strategy to avoid the induction of an adaptive immune response in the host.
this isolate was included because it harbours ismpa <dig> and has a similar is <dig> rflp pattern to strains shown to be associated with virulence
<cit> .
further studies, including experimental infection models, could provide an answer to whether the “silencing approach” of mah  <dig> might be the most successful strategy.
other possible explanations for the low prevalence of maa observed in the human population, and also in pigs, include the lack of exposure or a reduced ability to cross the human intestinal barrier.
in the norwegian pig population only mah has been isolated
<cit> , but the virulence of maa in pigs has been demonstrated by others
<cit> , and previous research demonstrated that maa  <dig> and mah vi <dig> were equally able to infect piglets
<cit> .
it thus appears that a lack of exposure is the most likely explanation for the difference in prevalence in pigs.
the virulence of maa has also been demonstrated in mice
<cit> .
although the study of pedrosa et al.
was performed prior to the division of m. avium into subspecies, one could assume that the bird isolates used in the older studies are equivalent with maa, and the human isolates with mah.
the virulence of m. avium subspecies in humans can obviously not be determined in vivo.
however, altogether the studies in human primary cell lines and infection studies in other species, suggest that different exposure to mah and maa in humans is a plausible explanation for the difference in prevalence.
CONCLUSIONS
in conclusion, our findings demonstrate that the ability to replicate in macrophages cannot explain the difference in prevalence of the two subspecies in humans.
genome-wide expression profiling studies further showed that similar biological pathways and networks were elicited by the two subspecies, although responses to the maa subspecies were exaggerated.
availability of supporting data
the raw microarray data have been made publicly available in the array express repository http://www.ebi.ac.uk/arrayexpress/experiments/e-mtab-1101/.
the supporting data are included in the additional file.
competing interests
the authors state that there are no competing interests related to the present study.
authors’ contributions
aa was responsible for conception and design of the experiments, contributed to work on immunoblotting, performed work with cell-culture and bacterial quantification, data analysis and drafted the manuscript.
aj contributed to analysis of microarray data and critical revision of the manuscript.
af contributed to microarray analysis and critical revision of the manuscript.
tbj contributed to conception and design of the experiment and critical revision of the manuscript.
iah contributed to cell-culture, immunoblotting and critical revision of the manuscript.
bd contributed to conception and design and critical revision of the manuscript.
ab contributed to analysis of microarray data and drafting of the manuscript.
io was involved in conception and design of the experiment, performed cytokine profiling and contributed to drafting of the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1: table s1
list of genes differentially expressed in strain maa  <dig>  mah vi <dig> and mah  <dig> compared to uninfected cells.
click here for file
additional file 2: figure s1
volcano plots comparing gene expression patterns in human marophages cultured in the presence or absence of different isolates of mycobacterium avium.a: maa  <dig>  b: mah vi <dig>  c: mah  <dig>  in all three plots, the red data points show the location of genes that were upregulated in responses to isolate maa  <dig>
click here for file
acknowledgements
the project was funded by the research council of norway and the norwegian veterinary institute.
ab is supported by a brightspark foundation mccusker fellowship.
birgitta Åsjø and sharad pathak at gades insitute, haukeland hospital, bergen, are greatly acknowledged for sharing their experience with real-time pcr for determination of the number of mycobacteria per cell in infected macrophage cultures.
the microarray service was provided by the norwegian microarray consortium  at the national technology platform, and supported by the functional genomics program  in the research council of norway.

 the mycobacterium avium complex comprises a group of environmental and opportunistic intracellular pathogens with the ability to persist within macrophages and escaping the host’s killing mechanisms. the complex is constituted by the subspecies of m. avium together with m. intracellulare. two subspecies of clinical importance for humans and animals are m. avium subsp. hominissuis and m. avium subsp. avium. mah is the subspecies mainly isolated from porcine and human cases of mac infection, while maa is the causative agent of avian tuberculosis, but is occasionally responsible for m. avium infections in mammals. little is known about the routes of transmission and the virulence mechanisms of m. avium, and a zoonotic potential has not been ruled out. definite reasons for mah being the dominating subspecies in humans and pigs infected with m. avium have not been identified, although differences in exposure could be considered a contributing factor. explanations to the higher prevalence of mah in the human and porcine population might also be found in differences between mah and maa at the molecular level of interaction between the host and mycobacterium, a process crucial to the outcome of infection, which therefore might affect the prevalence of the subspecies. however, little has been done to investigate whether mah and maa interact differently with their target cells. mycobacterial virulence depends on the ability of these bacteria to invade, persist and replicate within the hostile environment in macrophages of the host organism. this is achieved by the protective properties of the mycobacterial cell wall and by interference with the host’s immune response. to assess the ability of mah and maa to grow intracellularly, the use of primary cells will give a more representative result compared to immortalised macrophage cell-lines. tissue resident macrophages represent a heterogeneous population and it is not clear whether m. avium infect a particular phenotype. monocyte derived macrophages matured in absence of any stimuli give rise to a mixed culture of both the m and m phenotype and such cultures may thus be considered appropriate for infection studies with m. avium. finally, using macrophages derived from the relevant target species is important since the global inflammatory response evoked in mice recently was demonstrated to be vastly different from humans. to the authors’ best knowledge, comparative studies of the ability of maa and mah to invade, replicate and trigger immune responses in human primary macrophages have not been published. the aim of the present study was therefore to assess the growth of m. avium isolates in human primary macrophages and to characterise the response triggered in infected macrophages by selected isolates. the isolates used in the present study represented the two different subspecies, were derived from different host species and differed in biofilm abilities, glycopeptidolipid genes, and in the presence of the insertion sequence element ismpa. the present data showed that maa replicated to a larger extent in human macrophages suggesting that an inability to grow intracellularly is not likely the cause of the discrepancy of the prevalence of the two subspecies in humans. furthermore, the response induced in the macrophages could not entirely explain the difference in growth rate. methods mycobacterial isolates four clinical isolates of m. avium subsp. hominissuis and two clinical isolates of m. avium subsp. avium were used in the present study. details on origin, biofilm producing abilities, presence of ismpa nsgpl- and ssgpl genes were described by johansen et al., and are summarised in table  the human isolates were from a strain collection at the the norwegian institute of public health. all isolates were grown in middlebrook 7h w/oadc, 1% tween and 10% glycerol under agitation for two weeks before aliquots were stored at −70°c in an equivalent medium containing 10% glycerol. the concentration of bacteria in the stock suspensions was determined by colony forming unit counts. bacterial inocula for the gene expression assay were prepared from seven days subcultures grown on plates of middlebrook 7h with 10% oadc at 37°c. to adjust the number of bacteria in the inoculum, hsp real-time pcr was performed on serial dilutions in hbss of a suspension adjusted to mcfarland standard, homogenised through a 23 g needle to minimise clumping. mpa1 isolates of mycobacterium avium used in the present study show different combinations of subspecies designation, origin, biofilm abilities and the genomic presence of ismpa and genes involved in the synthesis of non-specific or serovar specific gpls. cell culture peripheral blood mononuclear cells were isolated from buffy coat drawn from healthy human blood donors. the material was commercially obtained from a blood donation centre, and the study was approved by the regional committee for medical research ethics, south-east norway. fifty ml buffy coat was diluted 1: with rpmi with 1% l-glutamine. pbmcs were harvested by density gradient centrifugation with lymphoprep and resuspended in pbs with % bsa and 2 mm edta at ph . isolation of cd14+ cells was performed using cd macs magnetic beads according to the manufacturer’s instructions. the cells were resuspended in rpmi with 1% l-glutamine supplemented with 10% human serum to the desired concentration of the respective infection assay, added to the appropriate costar® cell culture plates and incubated at 37°c with 5% co the cd14+ cells were incubated overnight before used in the uptake and replication assays and allowed to mature into macrophages by incubation for five days before used in the gene expression study. uptake and intracellular replication uptake and intracellular replication of the six m. avium isolates in human blood derived cd14+ cells were assessed. the assays were performed in 96-well cell culture plates, which in our hands gave a stable monolayer for most donors. each isolate was tested in triplicates at an moi of 5: the cells were washed twice with pre-warmed medium to remove extracellular bacteria before the quantification of intracellular bacteria was performed. intracellular uptake of the bacterial isolates was assayed after three, six, and hours in cells from three donors. intracellular replication of the same isolates was assayed in six other donors, following the protocol described above. based on the results from the uptake studies, cd14+ cells were incubated with the bacterial isolates for six hours. after removal of extracellular bacteria, fresh medium was added and the cells were incubated further. cell lysates were harvested at six hours, and after one, four and seven days of infection. quantitation of bacteria and cells numbers of intracellular bacteria and cells was determined by enumeration of mycobacterial and human genomes in cell lysates. replication of bacterial isolates in the supernatant had been tested and found insignificant. lysis of cells was achieved through adding milliq water to the wells. after minutes, lysates were transferred into 2 ml o-ring vials containing 200 μl silica beads and stored at −20°c. lysates from triplicate wells were pooled. after completing all sampling, lysates were thawed and inactivated at 100°c for 30 min, prior to mechanical disruption using the minibeadbeater for 45 sec. enumeration of mycobacteria and cells was achieved by performing a duplex real time pcr of a 103 bp long segment of the mycobacterial hsp gene and a 150 bp long segment of the human β-globin gene on the lysates, as described by salte et al.. briefly, the segment of the mycobacterial gene hsp was amplified using the primers mycofp and tb, and the vic tagged taqman-mgb probe mycopr. the segment of the human β-globin gene was amplified by using the primers bg-f and bg-r, and the fam tagged taqman-mgb probe bg-taq. dilutions of known concentrations of commercially available mycobacterial dna and human dna were run together with the samples and used to create standard curves. the conversion of the known amount of standard dna to number of bacterial and human genomes were described by salte et al., as well as the validation of the number of mycobacterial targets against cfu counts. total volume of each reaction mixture was 20 μl, consisting of 8 μl lysate, 10 μl quantitect® mastermix and 1 μl of a primer/probe mix for each target, with a stock concentration of 8 μm for each primer and 4 μm for the probe, resulting in a final primer- and probe concentration of  μm and  μm, respectively. samples and standard dna dilutions were run on a stratagene mx3005p mxpro v and data analysed by the stratagene mxpro v software with the following thermal cycling conditions: 95°c for 15 min, followed by cycles of 95°c for 1 min and 60°c for 1 min. gene expression profiling isolates mah vi mah and maa were used to infect macrophages for gene expression analysis. cd14+ cells were obtained as described above and allowed to mature into macrophages by incubation for five days in duplicate wells. the cells were then infected with the three isolates at a moi of 10: for four hours. total rna was isolated of the adherent macrophages using the rneasy mini kit®, following the manufacturer’s instructions. rna from the different donors within each assay was pooled per bacterial isolate. the first assay was performed in macrophages from three donors and subsequently repeated twice in two new donors per assay. the quality and integrity of rna was confirmed by agilent bioanalyzer and nanodrop™. gene expression levels were profiled using illumina humanht- expression bead chip. the gene expression data was exported from illumina’s genomestudio software package without background subtraction or normalization. the raw data was preprocessed in the r language and environment for statistical computing employing the normexp by control algorithm. probes that were not significantly expressed above background were filtered out of the analysis. differentially expressed genes were identified employing bayesian/moderated t-statistics with false discovery rate control for multiple testing. differentially expressed genes were identified at fdr < . the raw microarray data have been made publicly available in the array express repository. bioinformatics biological functions enriched in the differentially expressed genes were identified using modular enrichment analysis with david. enrichment scores >  were considered statistically significant. biological pathways enriched in the data were identified using ipa software. upstream regulator analysis was performed in ipa software to infer the putative pathways that give rise to the observed gene expression changes. molecular interaction networks were reconstructed from the differentially expressed genes in ipa software using mechanistic data from prior studies. in this analysis differentially expressed genes were defined at fdr < . protein analysis cd14+ cells were isolated and matured as described above for gene expression studies. cell lawns were infected at an moi of 5: for hours with maa mah vi and mah supernatants were removed for cytokine profiling, and subsequently the adherent cells were washed twice in pbs and collected with a cell scraper in 100 μl of 1x sds sample buffer for immunoblotting. all samples were stored at −70°c until further processing. sds-page and immunoblotting proteins were separated by horizontal sds-page with precast gradient excelgel™ sds gradient – 18% and excelgel™ sds buffer strips in a multiphor™ ii, prior to diffusion blotting onto nitrocellulose membranes with pore size  μm. proteins remaining within the gel were visualized by coomassie brilliant blue staining. membranes were blocked with pbs with 2% bsa and incubated over night with primary antibodies recognizing the following selected proteins involved in apoptosis: cleaved caspase-, cleaved caspase-, cleaved caspase-, cleaved caspase-, bcl- phospho-bcl-. detection of bound antibodies was performed by incubaction with ecl™ anti-rabbit igg, horseradish peroxidise labelled linked f fragment from donkey and visualization by immun-star™ westernc™ kit. cytokine profiling the concentration of tnf-α, il- il- il- and il- was measured in supernatants by using bio-plex™ cytokine assay within the bio-plex suspension array system, following the manufacturer’s instructions. differences between isolates were tested for significance using the wilcoxon matched pairs signed-rank test, considering p-values ≤  significant.  uptake and intracellular replication all the isolates were taken up by the cd14+ cells, and intracellular bacteria were detected after hours. maa mah h and maa appeared to be taken up to a greater extent, while the mah isolate vi showed the lowest degree of uptake in macrophages. between six and hours of infection some of the isolates showed an exponential increase in the number of intracellular bacteria, suggesting the event of onset of intracellular replication. six hours was therefore chosen as the time between inoculation and removal of extracellular bacteria for the subsequent study of intracellular replication over time. at this point in time vi had significantly lower numbers of intracellular bacteria compared to the other isolates except for mah h while no significant difference were observed between the other isolates using the paired t-test. intracellular growth was only observed for the two maa isolates. there was a significant increase in the ratio of mycobacteria to human genomes at day four for isolate maa and at day for maa and maa and in absolute numbers of bacteria for isolate maa at day using the paired t-test. the remaining isolates, which were all mah, did not show significant increase in intracellular numbers throughout the duration of infection, however the count of intracellular bacteria was quite stable, indicating that the bacteria were able to persist within the macrophages. enumeration of both cells and intracellular mycobacteria was performed by real-time pcr of single-copy genes. this method does not distinguish between viable and non-viable bacteria and cells, but it has the advantage that the number of bacteria is not influenced by the degree of clumping which varies between maa and mah. furthermore it enabled the measurement of the number of intracellular bacteria per cell in addition to absolute counts of bacteria. some degree of cell death and loss of adherence were observed in both infected and non-infected cells, and varied considerably between donors. similar results were observed whether absolute numbers of bacteria or the bacteria per cell ratio were used, however the difference in cell death between donors may explain why the absolute counts of bacteria only reached a significant level of p <  for maa at day while the ratio was significant for both maa isolates at day and for maa at day gene expression analysis and bioinformatics based on the results from the initial replication assays, three isolates were selected for gene expression studies in macrophages. maa which replicated intracellularly, and mah vi which only persisted within cells, were selected, in addition to mah which, as opposed to the other two isolates, contains the is element ismpa and belong to an rflp type that has been demonstrated to be particularly virulent in pigs. the illumina® humanht- v expression bead allows whole-genome expression analysis of samples on the same chip, each array targeting genomewide with more than probes. of the probe sets with signal intensity above background, probe sets were differentially expressed in the cells infected with m. avium, regardless of the individual isolate, compared to uninfected control cells. the largest response was induced by maa where probe sets were differentially expressed. of these probe sets, were up-regulated, and were down-regulated. mah vi modulated expression of probe sets representing unique genes. in comparison, mah altered the expression of only probe sets representing unique genes. a large overlap was observed in the responses elicited by the three isolates. there were virtually no genes that were uniquely modulated by isolate mah in contrast, hundreds of genes were uniquely differentially expressed in response to infection with maa and mah vi bioinformatics analyses demonstrated that the major biological functions up-regulated in response to the three isolates were apoptosis, inflammatory response, cytokine/chemokine activity, signal transduction, cell-proliferation and regulation of t-cell activation. biological pathways up-regulated in the responses were related to apoptosis and innate immunity and cytokine signaling. in addition, communication pathways between innate and adaptive immune responses and il- signaling were up-regulated. upstream regulator analysis was performed to infer the putative molecular mechanisms that give rise to the observed gene expression changes, and is illustrated in table  this analysis revealed major roles for tnf, nfkb, il-1b, ifng and others in the regulation of the responses. altogether these analyses revealed a strong overlap between the impacts of all the three isolates. the analysis was done using the database for annotation, visualization and integrated discovery. an enrichment score >  was considered statistically significant. pathways induced by infection of human monocyte-derived macrophages with m. avium isolates. the analysis was done in ingenuity systems software. the ten top pathways are shown. genes predicted to act as up-stream regulators for the observed responses after infection of human monocyte-derived macrophages with m. avium isolates. the analysis was done in ingenuity systems software. responses to chemicals/drugs are excluded. the ten top genes are shown. to identify the putative molecular mechanisms that are regulating the responses, we constructed molecular interaction networks for each isolate using mechanistic findings from previous studies. this analysis showed that maa induced the largest gene network in human macrophages, followed by mah vi then by mah. the networks were organized around a series of hyper-connected hubs, and the most important hubs are listed in the table  notably, the network connectivity patterns were dominated by tnf, il-1b, il- nfkb nfkbia, ptgs and others, and many of these pathways were also highlighted by the upstream regulator analysis. differential gene expression elicited by maa mah vi and mah 1655 at fdr < , no significant differences were observed between cells infected with maa and mah vi while only two genes were upregulated in maa compared to mah and one gene was upregulated in mah vi compared to mah protein analysis gene expression analysis pointed to apoptosis and immune response as two major biological functions upregulated in response to all three m. avium isolates, albeit at lower levels in mah compared to maa and mah vi we thus attempted to look at proteins involved in apoptosis by western blot and cytokine by a multiplex assay. however, differences between the isolates with regards to apoptotic proteins were not detected, probably due to the lower sensitivity in this assay compared to the microarray analyses. macrophages infected with all three isolates consistently produced more il- il- and tnf-α than uninfected controls. the production of il- and il- in uninfected cells was often high, and in some donors infection led to a decrease in cytokine production. the decrease in cytokine production was particularly evident for il- after infection with maa to a lesser extent by infection with mah and rarely by mah vi. there was a large donor to donor variation in the amount of cytokines produced, and this variation surpassed the difference induced by the isolates. nevertheless, mah vi consistently induced more il- il- and tnf-α than maa and more il- than mah. discussion the present study describes the uptake and intracellular growth of maa and mah isolates in human monocytes, and characterises genome-wide transcriptional responses elicited by these isolates in macrophages. mycobacterium avium isolates of different genotypic characteristics were all able to enter and persist within human primary monocytes. intracellular replication was only seen by maa suggesting that the differences in the ability to survive and replicate in macrophages cannot explain why humans are more often infected with mah. the gene expression program induced by infection with maa and mah isolates was broadly comparable in terms of the biological functions, pathways, and gene networks that were perturbed. however, the degree of perturbation varied considerably. the different level of expression led to an apparent discrepancy in the number of unique genes induced by the various isolates and the lack of differences when comparing the isolates with each other. many of the genes significantly induced by maa were also induced at a lower level by the other isolates. this intermediate expression was neither significantly different from uninfected cells nor from maa a larger samples size might have given more significant differences between the isolates since large donor to donor variation is present when using primary cell lines. in general, responses to maa were excessive, and were characterized by increased expression of genes downstream of type i and ii interferon signaling. further studies will be required to determine if variations in these pathways account for the exaggerated responses and/or differences in prevalence. similar pathways were induced in ruminant monocyte derived macrophages infected with the closely related mycobacterium avium subspecies paratuberculosis, while a study using the murine cell line raw showed that maa and mah induced higher levels of pro-inflammatory genes compared to map. together with our data, these studies suggest that similar pathways were induced by all m. avium subspecies. a direct comparison of the level of macrophage activation is however challenging considering the variation in cell cultures and isolates used in the various studies. macrophage activation can be interpreted as an appropriate reaction to an invading pathogen or as a failed attempt of the mycobacteria to silence the host’s defences. maa and mah vi differed completely in their ability to replicate within host cells, although both isolates induced a strong activation of macrophages and innate immune responses on the transcriptional level. analyses on the protein level revealed that the final synthesis of inflammatory cytokines were higher in response to mah vi compared to maa one could thus argue that the inflammatory response was sufficient to limit growth of mah vi but not maa the complexity of the outcome of an infection, and in vitro culture, was illustrated by the fact that despite a low activation of the innate immune response mah did not proliferate intracellularly. the microarray data provided no clear explanation to the lower level of inflammatory cytokines induced by maa on the protein level compared to mah vi however, regulation of inflammation on the post-transcriptional level has been demonstrated to be important, and a down-regulation of immune responses to avoid immunopathology can occur in chronic infections. the presence of unidentified virulence factors in maa more actively shutting down the innate responses cannot be excluded. the genomes of the three isolates did for instance differ in the presence of gpl genes, and the involvement of gpls in m. avium colony morphology and of ssgpls in pathogenicity is increasingly recognised. the induction of pro-apoptotic genes, such as ripk bid and tbid was seen after infection with all three isolates, however apoptotic pathways were affected to a lesser degree by mah than by maa and mah vi the correlation between apoptosis and virulence of mycobacterial strains is debated. inhibition of apoptosis has been described as a virulence factor in mycobacteria, but it has also been suggested that m. avium can use apoptosis of macrophages as a tool for spreading. one of the other more pronounced differences between the isolates was the increased expression of genes involved in activation of t cells like cd cd cd and il in maa and mah vi cd cd and cd are transmembrane surface proteins involved in activation of macrophages and t cells upon binding, while il- is essential for the persistence and function of th cells. these findings suggest that differential activation of apoptotic pathways and/or activation of t cells may occur in humans after infection with the different isolates. the lower activation of macrophages by mah may be a useful strategy to avoid the induction of an adaptive immune response in the host. this isolate was included because it harbours ismpa and has a similar is rflp pattern to strains shown to be associated with virulence. further studies, including experimental infection models, could provide an answer to whether the “silencing approach” of mah might be the most successful strategy. other possible explanations for the low prevalence of maa observed in the human population, and also in pigs, include the lack of exposure or a reduced ability to cross the human intestinal barrier. in the norwegian pig population only mah has been isolated, but the virulence of maa in pigs has been demonstrated by others, and previous research demonstrated that maa and mah vi were equally able to infect piglets. it thus appears that a lack of exposure is the most likely explanation for the difference in prevalence in pigs. the virulence of maa has also been demonstrated in mice. although the study of pedrosa et al. was performed prior to the division of m. avium into subspecies, one could assume that the bird isolates used in the older studies are equivalent with maa, and the human isolates with mah. the virulence of m. avium subspecies in humans can obviously not be determined in vivo. however, altogether the studies in human primary cell lines and infection studies in other species, suggest that different exposure to mah and maa in humans is a plausible explanation for the difference in prevalence.  in conclusion, our findings demonstrate that the ability to replicate in macrophages cannot explain the difference in prevalence of the two subspecies in humans. genome-wide expression profiling studies further showed that similar biological pathways and networks were elicited by the two subspecies, although responses to the maa subspecies were exaggerated. availability of supporting data the raw microarray data have been made publicly available in the array express repository http://www.ebi.ac.uk/arrayexpress/experiments/e-mtab-1101/. the supporting data are included in the additional file. competing interests the authors state that there are no competing interests related to the present study. authors’ contributions aa was responsible for conception and design of the experiments, contributed to work on immunoblotting, performed work with cell-culture and bacterial quantification, data analysis and drafted the manuscript. aj contributed to analysis of microarray data and critical revision of the manuscript. af contributed to microarray analysis and critical revision of the manuscript. tbj contributed to conception and design of the experiment and critical revision of the manuscript. iah contributed to cell-culture, immunoblotting and critical revision of the manuscript. bd contributed to conception and design and critical revision of the manuscript. ab contributed to analysis of microarray data and drafting of the manuscript. io was involved in conception and design of the experiment, performed cytokine profiling and contributed to drafting of the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1: table s1 list of genes differentially expressed in strain maa mah vi and mah compared to uninfected cells. click here for file additional file 2: figure s1 volcano plots comparing gene expression patterns in human marophages cultured in the presence or absence of different isolates of mycobacterium avium.a: maa b: mah vi c: mah in all three plots, the red data points show the location of genes that were upregulated in responses to isolate maa click here for file acknowledgements the project was funded by the research council of norway and the norwegian veterinary institute. ab is supported by a brightspark foundation mccusker fellowship. birgitta Åsjø and sharad pathak at gades insitute, haukeland hospital, bergen, are greatly acknowledged for sharing their experience with real-time pcr for determination of the number of mycobacteria per cell in infected macrophage cultures. the microarray service was provided by the norwegian microarray consortium at the national technology platform, and supported by the functional genomics program in the research council of norway.
the aim of the present study was to investigate the ability of the m. avium subspecies to replicate intracellularly and to characterise the gene expression program triggered by infection of human primary macrophages. isolates of maa were able to replicate intracellularly suggesting that differences in exposure, uptake or induction of adaptive immunity are more likely explanations for the difference in prevalence between m. avium subspecies. mycobacterium avium isolates of different genetic characteristics invaded monocytes and induced different degree of macrophage activation. although similar pathways and networks were perturbed by the different isolates, the response to the maa subspecies was exaggerated, and there was evidence of increased activation of type i and ii interferon signaling pathways. transcriptional responses to the isolates were characterized by upregulation of genes involved in apoptosis, immune- and inflammatory response, signal transduction and nf-kb signaling, cell proliferation and t-cell activation. however, intracellular replication was only evident in cells infected with the two maa isolates. all isolates were able to invade and persist within human macrophages.
 mycobacterium avium subsp. avium and m. avium subsp. hominissuis are environmental mycobacteria and significant opportunistic pathogens. mycobacterium avium infections in humans and pigs are mainly due to mah. it is not known whether this is caused by a difference in virulence or difference in exposure to the two subspecies. the aim of the present study was to investigate the ability of the m. avium subspecies to replicate intracellularly and to characterise the gene expression program triggered by infection of human primary macrophages.  all isolates were able to invade and persist within human macrophages. however, intracellular replication was only evident in cells infected with the two maa isolates. transcriptional responses to the isolates were characterized by upregulation of genes involved in apoptosis, immune- and inflammatory response, signal transduction and nf-kb signaling, cell proliferation and t-cell activation. although similar pathways and networks were perturbed by the different isolates, the response to the maa subspecies was exaggerated, and there was evidence of increased activation of type i and ii interferon signaling pathways.  mycobacterium avium isolates of different genetic characteristics invaded monocytes and induced different degree of macrophage activation. isolates of maa were able to replicate intracellularly suggesting that differences in exposure, uptake or induction of adaptive immunity are more likely explanations for the difference in prevalence between m. avium subspecies. mycobacterium aviumhuman macrophagesgene expression 