BACKGROUND
bacteria from the genus brucella are the etiological agents of brucellosis, a worldwide zoonotic infectious disease that has a negative economic impact on animal production and human public health  <cit> .
based on its 16s rrna sequence, brucella is included in the α <dig> subclass of the proteobacteria, along with plant  and other mammalian  symbionts  <cit> .
the genus brucella consists of six recognized species, grouped according to their primary host preferences, i.e.
b. abortus : cattle, b. melitensis : sheep and goats, b. suis : hogs, b. ovis : sheep, b. canis : dogs and b. neotomae : wood desert rats  <cit> .
due to their high virulence to humans, b. abortus, b. melitensis and b. suis are considered potential bioterrorist agents, having been classified as major biodefense/biothreat pathogens, and their possession and use is strictly regulated in the united states  <cit> .
natural brucella infections occur primarily through adhesion to and penetration of mucosal epithelia.
the mucosal surface of the alimentary tract is a major route for b. melitensis and b. abortus invasion, while the mucosa of the genital tract is the principal route of entry for b. ovis, b. suis and b. canis  <cit> .
in vitro studies have shown that within a few minutes after binding non-professional phagocytic cells, brucella are actively internalized via receptor-mediated phagocytosis without inducing obvious damage to the cells  <cit> .
brucella bind sialic acid residues present on eukaryotic cell membranes  <cit>  and are internalized by epitheloid-like cells in an active mechanism in which the organism induces its own internalization via activation of small gtpases of the rho subfamily and rearrangements of the host cell actin cytoskeleton and microtubules  <cit> .
bacteria have the ability to express surface molecules able to recognize unique or common receptor components present on many eukaryotic cell surface.
three brucella gene products have been characterized as important for invasion in non-phagocytic cells: a two-component regulatory system  that modulates the expression of outer membrane proteins necessary for recruiting small gtpase proteins required for actin polymerization and penetration  <cit> , a brucella surface protein, called sp <dig>  which enables brucella to adhere to non-phagocytic cells  <cit> , and a hypothetical protein encoded by the bmei <dig> gene, which is critical for brucella melitensis internalization in hela cells after  <dig> h post-infection  <cit> .
these few examples are all that is currently known about the molecular mechanisms underlying brucella adhesion and internalization in eukaryotic cells.
hela cells have extensively been used as a model to investigate the internalization of brucellae of epithelial cells during the colonization of the susceptible host  <cit> .
here, we employed this cell line to evaluate the rate of invasion of b. melitensis at different growth phases.
our results indicate that cultures of b. melitensis in the late-log phase of growth were more invasive in non-professional phagocytic cells than cultures at mid-log and stationary growth phases.
using cdna microarrays, we characterized the transcriptome of the most  and the least  invasive growth phases of b. melitensis cultures as a preliminary approach for identifying pathogen candidate genes involved in epithelial cell invasion process.
microarray analysis revealed a greater number of genes up-regulated in these cultures than in stationary phase cultures.
consistent with the expected differences due to growth, there was a more active metabolism and invasiveness of cultures in late-log phase than cultures in stationary phase.
given the role that some of these genes have in pathogenesis in other bacterial species, we believe that these data may offer insight into potential growth-phase regulated brucella virulence genes involved in the initial host:pathogen interactions.
RESULTS
b. melitensis  <dig> m at late-log phase of growth were more invasive to epithelial cells than were bacteria at mid-log and stationary growth phases
as described in the methods section, b. melitensis was grown to mid-log growth phase, late-log growth phase, or stationary growth phase.
at each of these growth phases, bacteria were enumerated, used to infect a representative epithelial cell line , and rna was extracted and microarrays were performed to identify altered gene expression.
under our experimental conditions, there were  <dig>  ×  <dig> cfu/ml  at the mid-log growth phase,  <dig> ×  <dig> cfu/ml  at late-log phase, and  <dig> ×  <dig> cfu/ml  at stationary phase .
for invasion experiments, a consistent multiplicity of infection  factor of  <dig>  b. melitensis cells per hela cell was used to normalize the number of bacteria used.
the average number of intracellular bacteria recovered was  <dig> cfu at mid-log phase,  <dig> cfu at late-log phase of growth and  <dig> cfu at stationary growth phase per  <dig> cells inoculated .
these values represent the average of three independent experiments.
b. melitensis  <dig> m cultures grown to late-log phase and then co-incubated with hela cells for  <dig> min were therefore  <dig>   and  <dig>   times more invasive than were cultures at mid-log and stationary growth phases.
whole-genome expression analysis of the most and the least b. melitensis  <dig> m invasive growth phases: reliability of array data
to analyze the molecular differences between the most and the least invasive phenotype, four biological replicates of cultures at late-log and stationary growth phases were analyzed using cdna microarrays.
genomic dna was used as an internal control for each experiment in order to allow experiment-to-experiment comparisons  <cit> .
as expected, there was little variability between gdna signals from array to array, even under the two different conditions examined .
the r <dig> value for any two arrays  was between  <dig>  and  <dig> , even before normalization.
when the values for each conditional replicate were averaged , the resulting r <dig> value was  <dig>  .
comparisons of rna cy <dig> fluorescent signals  yielded similar r <dig> values .
in order to further minimize the incidence of false positives and increase the consistency and reliability of the microarray analysis results, the data were analyzed separately using four different techniques: genespring combinatorial analysis, spotfire decisionsite  <dig>  pairwise comparisons, sam two-class unpaired comparisons, and anova.
a change in gene expression was considered significant if the p value was less than  <dig> , the fold-change was at least  <dig> , and the gene expression alteration occurred for all replicate experiments.
we further expected each gene to be significantly differentially expressed for at least two of the three replicate spots for each experimental array set .
based on these criteria, genes that were deemed significant by all four analytical methods  were organized by cogs functional categories  <cit>  and compiled into a list that included  <dig> genes  that were up- or down-regulated when b. melitensis was grown to late-log phase, compared to stationary phase .
a direct comparison of the signal intensity values of these genes indicated that the difference between log and stationary phases was specifically due to differential gene expression and not array spatial bias, as indicated in figure  <dig>  when the average gdna intensity values for these  <dig> genes were plotted , the r <dig> value was  <dig>  .
however, the r <dig> value for the same genes comparing the cy <dig> fluorescence values instead  was extremely low .
to confirm the microarray results, we randomly chose  <dig> differentially expressed genes  and performed qrt-pcr.
based on qrt-pcr results, transcript levels of  <dig> of these genes  were altered greater than  <dig> -fold and in the same direction as was determined by microarray analysis.
two other genes  were determined to be differentially expressed and in the same direction of microarray analysis, but the fold change was lower than  <dig>  no significant difference in the expression level of bmei <dig> was observed by qrt-pcr .
b. melitensis  <dig> m express different sets of genes in late-log and stationary phases of growth in f12k tissue culture medium
of the  <dig> genes significantly altered in b. melitensis during late-log phase ,  <dig>  were up- and  <dig>  were down-regulated, compared to when the bacteria were allowed to reach stationary phase .
the relative changes in gene expression ranged from a  <dig> -fold induction of the glycerol-3-phosphate regulon repressor gene  to a  <dig> -fold down-regulation of the locus bmeii <dig> .
as expected, the majority of gene expression changes were associated with growth and metabolism.
among the up-regulated genes were those associated with dna replication, transcription and translation , nucleotide, amino acid, lipid and carbohydrate metabolism , energy production and conversion , membrane transport  and cell envelope, biogenesis and outer membrane , while the  <dig> down-regulated genes were distributed among several cogs .
hierarchical cluster analysis using spotfire decisionsite  <dig>  was performed on normalized cy <dig>  signal intensity values of the microarray data from the four log phase and four stationary phase samples.
all four samples from the log phase of growth clustered together, apart from those collected at stationary phase .
moreover, genes that clustered together were indeed differentially expressed between the two growth conditions.
the higher number of genes up-regulated in late-log growth phase coincides with a more active metabolism of late-log phase cultures compared to those at stationary phase.
in the following sections, we will focus our comments on those genes differentially expressed by microarray analysis that encode or are predicted to encode virulence factors, some of which may be involved in brucella:host interaction.
protein-encoded genes which play a role in brucella invasiveness in non-phagocytic cells did not have differential expression between the most and the least invasive cultures
currently, only three brucella gene products have been characterized as important for invasion in non-phagocytic cells.
the b. abortus two-component regulatory system bvrr/bvrs encoded by bvrr/bvrs genes, regulates the structure of outer membrane components and plays a critical role in cell penetration and intracellular survival  <cit> .
this two-component system is highly conserved in the genus brucella  <cit> , with chvi/chvg  representing the b. melitensis homolog.
in this study, neither of the two genes that encode this two-component system were differentially expressed between the most and the least invasive b. melitensis cultures.
another brucella invasive-characterized gene product is sp <dig>  a surface protein that enables b. suis to attach and penetrate non-phagocytic cells  <cit> .
the role of this gene has not been evaluated in b. melitensis, although a homolog is encoded by the ugpb gene present on the chromosome ii of the b. melitensis  <dig> m genome .
in this study, ugpb was not differentially expressed when global gene expression of b. melitensis cultures at late-log phase was compared to cultures at stationary growth phase.
recently, a third gene product was reported to be involved in brucella internalization in non-phagocytic cells  <cit> .
in that study, a b. melitensis mutant with interruption in the bmei <dig> gene exhibited a marked decrease in its ability to invade hela cells at  <dig> and  <dig> h post-infection, suggesting the relevance of this gene in the brucella invasion process after  <dig> h p.i.
in this study, bmei <dig> was not found altered due to growth-phase.
collectively, these results indicate that the higher invasiveness observed in b. melitensis cultures at late-log phase of growth under our experimental conditions was not due to the differential expression of these three characterized gene products.
also, these results suggest that genes encoding these three products were not growth-phase regulated in b. melitensis under our experimental conditions.
however, they might be transcribed at a time that we did not measure, they could be constitutively expressed and act in concert with other factors, or they could be expressed following epithelial cell contact.
it is perhaps worth noting that only one of these three gene products  has been effectively demonstrated to contribute to b. melitensis virulence, although after one hour post infection rather than the  <dig> minutes used in this study.
well-known b. melitensis virulence genes had different expression profiles in late-log phase of growth compared to stationary growth phase
several genes whose products are known to be associated with brucella melitensis virulence , were differentially expressed between the most and the least invasive cultures.
these included genes that encode t4ss proteins and the flagellar apparatus.
the virb locus, for instances, encodes the type iv secretion system  and plays a critical role in brucella virulence and intracellular multiplication  <cit> .
three genes encoding components for the virb operon, such as virb <dig> , virb <dig>  and virb <dig>  were up-regulated in b. melitensis cultures at late-log phase compared to stationary growth phase.
pathogenic bacteria produce flagella to promote colonization and invasion of mucosa.
brucellae are traditionally characterized as non-motile bacteria, yet the sequence of the b. melitensis genome contains three clusters of flagellar genes  <cit>  and their participation in establishing chronic brucellosis has been established  <cit> .
in our study, five genes such as flic , flif , flin , flha  and flgd  which encode parts of the flagellar apparatus or regulate its expression, were differentially expressed in late-log phase cultures compared to stationary phase cultures.
previous studies reported scant influence of t4ss and flagella in the invasion process  <cit> .
thus, the highest penetration observed in late-log phase cultures was probably not due to the expression of these genes.
several transcriptional regulator genes were differentially expressed in late-log phase compared to stationary growth phase
transcriptional regulators control bacterial gene expression in response to specific signals.
twenty-two genes encoding transcriptional regulators belonging to the arac , asnc , beti , deor , gntr , iclr , lysr , luxr , marr , merr , and rpir  families were differentially expressed in late-log phase b. melitensis cultures compared to stationary phase cultures.
some of these transcription factors are known to be involved in positive regulation of gene expression .
others are involved in repression , while members of iclr and lysr families could be activators or repressors of gene expression  <cit> .
nevertheless, the contribution of these regulators and their targets to b. melitensis internalization epithelial cells has not been fully examined.
the locus encoding the alternative sigma  <dig> factor  that allows brucella to survive under general stress situations was up-regulated in stationary phase cultures.
the bmei <dig> locus that encodes a subunit of the other alternative sigma  <dig> factor , which allows transcription of those genes involved in utilization of nitrogen and carbon sources and energy metabolism, was up-regulated in late-log phase cultures compared to stationary phase cultures.
two-component transcriptional regulators are comprised of a cytoplasmic membrane-located sensor protein and a cytoplasmic response regulator protein  <cit> .
eight orfs encoding for two-component response regulators have been identified in the b. melitensis  <dig> m genome  <cit> .
one of the signal transduction-encoded genes up-regulated in late-log phase cultures , was previously identified in b. melitensis attenuated mutants  <cit> .
the other  is a central regulator of carbohydrate metabolism genes and also plays a role in virulence development of certain pathogens  <cit> .
although the molecular regulation of these response regulators in b. melitensis is currently unknown, understanding how vsr, hprk and others are regulated, could offer insight into b. melitensis virulence.
identifying the target genes of these transcriptional regulators would significantly clarify the role of growth-phase in brucella physiology, metabolism and virulence regulation.
almost all differentially expressed genes encoding cell envelope and outer membrane components were up-regulated in late-log phase cultures
the ability of brucella to invade cells has been linked to its outer membrane  properties, as well as to structures built within the cell envelope  <cit> .
twenty-six genes directly involved in cell envelope and outer membrane biogenesis were differentially expressed at late-log compared to stationary phase of growth.
these included genes that encode outer membrane proteins , lipoproteins , lps , and peptidoglycan biosynthesis .
the main cogs functional category of genes that were up-regulated in b. melitensis cultures at late-log compare to stationary phase of growth were orfs encoding membrane transport proteins.
these included genes encoding transporters specific for amino acids , carbohydrates  and uncharacterized transporters .
collectively, these data indicate an active conversion of metabolites to components of the cell envelope structure at late-log phase, which given the importance of the om in virulence, could influence the initial brucella:host cell interactions, facilitating attachment and entry into host cells.
discussion
the molecular mechanisms involved in the initial interactions between brucella and epithelial cells have not been well characterized.
previous studies have used hela cells as a model for studying adhesion and internalization of brucella spp.
in non-professional phagocytic cells  <cit> .
these studies found that brucellae bind to cellular receptors containing sialic acid residues and induce their own uptake by a local rearrangement of the host cell cytoskeleton around the invading organisms.
the ability of the bacteria to adhere to and penetrate eukaryotic cells is a well orchestrated process that requires several factors/gene products in order to be successful  <cit> .
to date, only a few brucella gene products involved in non-phagocytic cell invasion have been identified  <cit> .
this study was performed with the goal of better understanding initial molecular interactions between brucella and its host through the molecular analysis of growth phase-specific gene regulation.
our initial experiment indicated that cultures of b. melitensis at late-log growth phase in cell culture medium were more invasive to non-phagocytic cells than cultures at mid-log and stationary growth phases.
similar results have been observed for other invasive pathogens, such as salmonella spp.
or yersinia enterocolitica  <cit> .
even with the high moi used , b. melitensis were internalized in lower numbers by epithelioid-like hela cells at  <dig> min p.i.
than reported in another study  <cit> .
the difference in invasion may have been influenced by the f12k cell culture medium used to growth the agent.
b. melitensis reach stationary phase at a lower od  in f12k cell culture medium than in rich bacterial culture medium  or another cell culture medium  .
these results suggest that f12k medium apparently contains suboptimal nutrients for brucella development.
even though, we grew b. melitensis in f12k medium and immediately added the bacteria to hela cells with the goal of reducing bacterial pre-infection manipulations , which had probably modified the original transcriptome of the cultures, since bacterial gene expression changes quickly in response to environmental modification  <cit> .
the relationship between growth phase and invasiveness is dependent upon the expression of bacterial virulence factors at different growth-phase.
for instance, shigella flexneri is more invasive during the early phase of their exponential growth, because invasion proteins  are secreted in higher amounts during this growth phase  <cit> .
salmonella enterica and legionella pneumophila have their secretion systems assembled and effector proteins properly stored in the cytoplasm only at the late exponential and stationary growth phases, respectively  <cit> .
in order to understand why our system evoked greater invasiveness in b. melitensis cultures at late-log phase in the first  <dig> min p.i., we conducted a global gene expression detection study using cdna microarray technology.
microarray analysis revealed that  <dig> genes were significantly differentially expressed between the most  and the least  invasive cultures .
as expected, the majority of the observed changes in gene expression were related to the bacterial response under the increased growth conditions in tissue culture media.
for example, the up-regulation of genes associated with transcription and translation, nutrient metabolism, transport, and energy production and conversion all correspond to a more active metabolism of late-log phase cultures, compared to cultures at stationary phase.
as was expected, several cell division- and dna synthesis-related genes were also up-regulated at late-log phase, when the bacterial population was still actively growing.
alternatively, genes down-regulated in late-log phase were more heterogeneous in nature, demonstrating no predominant functional category.
as expected, an increased expression of the locus bmei <dig>  encoding the alternative sigma  <dig> factor was observed in stationary phase cultures  <cit> .
sigma  <dig> factor regulates the transcription of heat shock genes, which allow the bacteria to survive not only an abrupt increase in temperature, but also general stress situations, such as nutrient limitations during stationary growth phase  <cit> .
previous work identified a role in b. melitensis invasion of hela cells for the hypothetical protein encoded by bmei <dig> orf, which increases invasiveness only after  <dig> h p.i.
<cit> .
that study clearly showed that the presence or absence of the gene transcript did not modify the ability of b. melitensis to invade hela cells during the first  <dig> min p.i., i.e.
the brucella-hela co-incubation time used in our study.
under our experimental conditions, bmei <dig> was not found phase growth regulated.
these data suggest that bmei <dig> may be transcribed after prolonged host cell contact, thereby facilitating the invasion process at later time points.
further characterization of the regulation of this gene and its product is clearly warranted.
in seeking to identify possible contributors to the increased invasiveness of b. melitensis at late-log phase, the conversion of metabolites to components that alter cell envelope structure were evaluated.
altered outer membrane/cell wall topology would be expected to influence the initial bacteria:host cell interaction that may facilitate attachment and entry into host cells.
in our study, a significant number of genes directly involved in cell envelope/outer membrane biogenesis were differentially expressed .
one of the genes up-regulated at late-log growth phase was the locus bmei <dig>  the product of this gene has not yet been characterized in b. melitensis; however, it has high homology  to an immunogenic outer membrane protein, omp <dig>   <cit> .
omp <dig> is a haemin-binding protein  <cit> , which binds to, and extracts iron from, the host.
iron has been identified as a required element for epithelial invasion in microbial pathogens  <cit> , and the expression of this locus, along with other iron-related genes in late-log phase cultures , may influence the internalization ability of brucellae.
sp <dig> is another surface-exposed outer membrane protein with a critical role in brucella suis adherence to, and invasion of, non-phagocytic cells  <cit> .
the role of this protein, which is encoded by the ugpb gene  present in the chromosome ii of b. melitensis  <dig> m genome, was not previously described for b. melitensis adhesion to and/or penetration of epithelial cells.
the transcript from the ugpb gene was not identified as differentially expressed in our cdna microarray analysis between the most and the least invasive cultures.
therefore, under our experimental conditions, this omp seems not to be involved in the higher invasiveness of the late-log phase cultures.
it is possible that the composition of the cell culture medium does not induce the expression of ugpb, or it is also possible that ugpb is constitutively expressed and/or act in concert with other factors.
although genetic analysis reveals that ugpb may belong to an operon  that encodes for a sn-glycerol-3-phosphate abc transporter  <cit> , the experimental evidence does not support this hypothesis.
a previous study showed that the product of ugpb in b. suis is indeed a surface-exposed protein with adhesion and invasion activity  <cit> .
in fact, in this study, three of the transcripts predicted to encode the transport system  were highly up-regulated  in late-log phase cultures, when compared to stationary phase cultures.
in concordance with previous experimental evidence, our microarray data would support the finding of others that ugpb does not belong to an operon that encodes for a sn-glycerol-3-phosphate abc transporter.
in addition, our results support growth-phase regulation of the sn-glycerol-3-phosphate abc transport system, which has been implicated in brucella pathogenesis  <cit> .
the ability of brucella to invade host cells is linked to its om properties.
b. melitensis omp profile changes during culture growth  <cit> , as gene expression is transcriptional regulated by environmental conditions  <cit> .
although the different expression pattern observed in gene-encoding cell envelope products between late-log and stationary phase cultures here may have been due to medium composition, it is more likely that they are due to growth phase regulation.
in this study, we did not evaluate the role of the omp in internalization in epithelial cells and therefore their individual participation in increased invasiveness of late-log phase cultures could not be determined.
only two differentially expressed genes encoding for o-chain and peptidoglycan layer biosynthesis from this study , were previously evaluated in brucella pathogenesis , although not in epithelial cells internalization  <cit> .
due to the importance that the cell envelope in initial host:pathogen interaction, the regulation and role of gene-encoding om products differentially expressed in this study should be addressed in future studies.
rapid adaptive physiological response to multiple environmental and cellular signals in bacteria is mainly mediated by transcriptional regulators and two-component regulatory systems.
prokaryotic genes putatively coding for transcriptional regulators are grouped in families based on sequence similarity and functional criteria.
twenty-two transcripts, belonging to  <dig> families of transcriptional regulators, were differentially expressed in our study .
it was recently reported that b. melitensis mutants for  <dig> of these  <dig> transcriptional regulators were not attenuated after one-week of infection in mice  <cit> .
however, effects of these transcriptional regulators on internalization of b. melitensis by non-phagocytic cells have not been examined.
their contribution to invasion therefore remains unknown.
luxr is a well-known family of transcriptional activators that regulates various functions in microbes  <cit> .
there are two loci  that encode transcripts belonging to this family of transcriptional regulators in the b. melitensis genome, and their expression is required for transcription of virulence factors such as virb operon and flagella  <cit> .
the transcriptional regulator vjbr was not differentially expressed in our study, but the other luxr homolog , was 221-fold up-regulated in the late-log phase of growth, compared to stationary phase cultures.
the targets of blxr are currently unidentified, but regulatory effects on other transcriptional-regulatory proteins and proteins predicted to be involved in cell envelope biogenesis was observed  <cit> .
it may be possible that some of these gene products regulated by blxr positively influence b. melitensis invasion of hela cells.
analysis of the invasive phenotype of a b. melitensis blxr deletion mutant in hela cells would be the first step in determining the importance of this transcriptional regulator during the initial host:brucella interaction, followed by the identification of the effector gene it regulates.
bvrr/bvrs is a well characterized two-component regulatory system that controls the expression of genes essential for brucella abortus invasion to non-phagocytic cells  <cit> .
high level of identity is present between b. melitensis chvi/chvg  and the b. abortus bvrr/bvrs proteins  <cit> .
in our study, no transcriptional change was observed in bmei2036/i <dig> orfs between the most  and the least  invasive cultures.
likely, brucella maintain a basal expression level of the regulatory locus, as a change in the phosphorylation of the protein required for activity rather than transcription.
twenty orfs dedicated to signal transduction were identified in b. melitensis genome  <cit> .
the importance of some of them in brucella virulence had been characterized lately, including blxr, vjbr, ftcr and bvrr/bvrs  <cit> .
however, their contribution to internalization in non-phagocytic cells is less known.
recently, mutants with defective expression in two transcriptional regulators  had an altered pattern in initial host:pathogen interaction due to surface modifications  <cit> .
future identification of the target genes of these regulators would clarify brucella physiology, metabolism and virulence regulation.
several motility-related genes were more highly expressed at late-log phase compared to stationary phase, including kinesin-like protein, chemotaxis motd protein and genes related to flagellum apparatus synthesis and functions, e.g.
flagellin itself .
flagellin has been well-characterized as a contributor to bacterial virulence through chemotaxis and adhesion to and invasion of host cells  <cit> .
the extent to which flagellar machinery participates in the process of invasion seems to depend at least partly on the species of bacteria and/or the host cell type.
for instance, flagellar-associated motility in salmonella is not required but accelerates invasion of caco- <dig> colonic epithelial cells  <cit> , whereas the invasion of acanthamoeba astronyxis by burkholderia pseudomallei absolutely requires an intact flagellum apparatus  <cit> .
in the case of b. melitensis, a previous study demonstrated that expression of flagella is growth curve-dependent and required for persistent disease in a mouse model but not for invasion in cellular models  <cit> .
that study reported that a functional flagellum was assembled in early-log growth phase cultures but not at later time points.
in our study, we did not analyze gene expression at early time points of the growth curve, but the results indicated that some flagellar genes were expressed more in late-log phase cultures as compared to stationary phase cultures.
the differences in flagellum gene expression between the study of fretin et al.
and ours could be attributed to evaluation of different steps of the process , different culture media used, or post-transcriptional regulatory mechanisms.
the flagellar apparatus is built hierarchically under complex regulation.
thirty-one flagellar genes distributed in three clusters on chromosome ii and along with three transcriptional regulators of flagellar system expression have been identified in b. melitensis  <cit> .
however, the order of flagella gene expression and the whole system regulation in brucellae has not been established.
here, only five genes from two loci encoding different parts of the flagellar apparatus were differentially expressed in late-log phase cultures compared to stationary phase cultures.
detection of expression of some but not all genes from an operon is not uncommon with microarray data, due to the inherent nature of microarrays  that produces variation that leads to some false negatives  <cit> .
in a previous study, rambow-larsen et al.
using a cdna microarray, also identified only  <dig> of the  <dig> flagellar genes, belonging to different flagellar loci and encoding for distinct parts of the flagellar apparatus, expressed under a putative quorum-sensing regulator blxr  <cit> .
similarly, microarray detected changes in expression of only some of the genes of the flagellar operon in salmonella enterica serovar typhimurium, which is transcribed with a polycistronic message, despite a 10-fold difference in some genes of each operon  <cit> .
two different functions, motility and protein secretion have been ascribed to flagella, but these roles have yet to be demonstrated in brucellae.
we were not able to evaluate the role of b. melitensis flagellar gene expression in invasion under our experimental conditions, but undoubtedly, the presence of flagellar machinery and other adhesion/motility factors at late-log phase, and their exact contribution to the brucella invasion process warrant further studies.
the virb operon has been reported to be essential for intracellular survival and multiplication of brucella  <cit> , but its role in adherence and internalization is contradictory  <cit> .
in our study, three genes from the operon  were up-regulated in late-log growth phase cultures compared to the stationary phase of growth.
virb is transcribed as an operon, with no secondary promoters.
it is maximally expressed in b. melitensis at the early exponential phase of the growth curve, and its expression decays as the bacteria reach the stationary phase  <cit> .
however, the half-lives of the individual segments of the virb transcript are not known.
under our experimental conditions, it is possible that virb was expressed earlier in the growth curve, and the different rate of transcript degradation allowed the detection of expression of some genes of the operon in late-log phase but not in stationary phase cultures.
distinct half-lives of the individual segments of the polycistronic mrna, results in the differential expression of each gene in the operon.
this phenomenon has been well characterized in other bacteria  <cit> , and is worthy to additional evaluation of b. melitensis virb operon.
in addition, and similar to mention for flagellar genes, microarray could detect expression of some but not all genes from an operon, due to the inherent nature of the technique.
further, our analysis method was particularly stringent in order to greatly reduce false positives at the risk of additional false negatives.
thus, other genes in the virb operon were increased in expression such as virb <dig>  virb <dig>  virb <dig>  virb <dig> and virb <dig>  although not statistically significant because of the stringency of our statistical analysis.
finally, genes with uncharacterized function that were differentially expressed at late-log phase compared with the stationary phase also deserve some special consideration.
this group of "hidden genes" represents 22% of the differentially expressed genes identified in this study, and it may contain some of the heretofore unknown virulence factors utilized for b. melitensis to invade and infect the host, as was previously suggested  <cit> .
conversely, brucella internalization should not be disregarded as a product of synergistic action among several gene products in non-phagocytic cells.
CONCLUSIONS
our study reveals that b. melitensis grown in cell culture medium at late-log phase are more invasive in non-phagocytic cells than cultures grown at mid-log or stationary growth phases.
cdna microarrays provide informative differential transcriptional profiles of the most  and the least  invasive b. melitensis cultures.
we consider these data a platform for conducting further studies on the brucella:host initial interaction.
since the roles of the majority of differentially expressed genes in this study are not well defined in brucella pathogenesis, future studies on brucella virulence can now be specifically focused to more precisely delineate the roles of candidate genes identified in this study.

 bacteria from the genus brucella are the etiological agents of brucellosis, a worldwide zoonotic infectious disease that has a negative economic impact on animal production and human public health. based on its 16s rrna sequence, brucella is included in the α subclass of the proteobacteria, along with plant and other mammalian symbionts. the genus brucella consists of six recognized species, grouped according to their primary host preferences, i.e. b. abortus: cattle, b. melitensis: sheep and goats, b. suis: hogs, b. ovis: sheep, b. canis: dogs and b. neotomae: wood desert rats. due to their high virulence to humans, b. abortus, b. melitensis and b. suis are considered potential bioterrorist agents, having been classified as major biodefense/biothreat pathogens, and their possession and use is strictly regulated in the united states. natural brucella infections occur primarily through adhesion to and penetration of mucosal epithelia. the mucosal surface of the alimentary tract is a major route for b. melitensis and b. abortus invasion, while the mucosa of the genital tract is the principal route of entry for b. ovis, b. suis and b. canis. in vitro studies have shown that within a few minutes after binding non-professional phagocytic cells, brucella are actively internalized via receptor-mediated phagocytosis without inducing obvious damage to the cells. brucella bind sialic acid residues present on eukaryotic cell membranes and are internalized by epitheloid-like cells in an active mechanism in which the organism induces its own internalization via activation of small gtpases of the rho subfamily and rearrangements of the host cell actin cytoskeleton and microtubules. bacteria have the ability to express surface molecules able to recognize unique or common receptor components present on many eukaryotic cell surface. three brucella gene products have been characterized as important for invasion in non-phagocytic cells: a two-component regulatory system that modulates the expression of outer membrane proteins necessary for recruiting small gtpase proteins required for actin polymerization and penetration, a brucella surface protein, called sp which enables brucella to adhere to non-phagocytic cells, and a hypothetical protein encoded by the bmei gene, which is critical for brucella melitensis internalization in hela cells after h post-infection. these few examples are all that is currently known about the molecular mechanisms underlying brucella adhesion and internalization in eukaryotic cells. hela cells have extensively been used as a model to investigate the internalization of brucellae of epithelial cells during the colonization of the susceptible host. here, we employed this cell line to evaluate the rate of invasion of b. melitensis at different growth phases. our results indicate that cultures of b. melitensis in the late-log phase of growth were more invasive in non-professional phagocytic cells than cultures at mid-log and stationary growth phases. using cdna microarrays, we characterized the transcriptome of the most and the least invasive growth phases of b. melitensis cultures as a preliminary approach for identifying pathogen candidate genes involved in epithelial cell invasion process. microarray analysis revealed a greater number of genes up-regulated in these cultures than in stationary phase cultures. consistent with the expected differences due to growth, there was a more active metabolism and invasiveness of cultures in late-log phase than cultures in stationary phase. given the role that some of these genes have in pathogenesis in other bacterial species, we believe that these data may offer insight into potential growth-phase regulated brucella virulence genes involved in the initial host:pathogen interactions.  b. melitensis m at late-log phase of growth were more invasive to epithelial cells than were bacteria at mid-log and stationary growth phases as described in the methods section, b. melitensis was grown to mid-log growth phase, late-log growth phase, or stationary growth phase. at each of these growth phases, bacteria were enumerated, used to infect a representative epithelial cell line, and rna was extracted and microarrays were performed to identify altered gene expression. under our experimental conditions, there were × cfu/ml at the mid-log growth phase, × cfu/ml at late-log phase, and × cfu/ml at stationary phase. for invasion experiments, a consistent multiplicity of infection factor of b. melitensis cells per hela cell was used to normalize the number of bacteria used. the average number of intracellular bacteria recovered was cfu at mid-log phase, cfu at late-log phase of growth and cfu at stationary growth phase per cells inoculated. these values represent the average of three independent experiments. b. melitensis m cultures grown to late-log phase and then co-incubated with hela cells for min were therefore and times more invasive than were cultures at mid-log and stationary growth phases. whole-genome expression analysis of the most and the least b. melitensis m invasive growth phases: reliability of array data to analyze the molecular differences between the most and the least invasive phenotype, four biological replicates of cultures at late-log and stationary growth phases were analyzed using cdna microarrays. genomic dna was used as an internal control for each experiment in order to allow experiment-to-experiment comparisons. as expected, there was little variability between gdna signals from array to array, even under the two different conditions examined. the r value for any two arrays was between and, even before normalization. when the values for each conditional replicate were averaged, the resulting r value was. comparisons of rna cy fluorescent signals yielded similar r values. in order to further minimize the incidence of false positives and increase the consistency and reliability of the microarray analysis results, the data were analyzed separately using four different techniques: genespring combinatorial analysis, spotfire decisionsite pairwise comparisons, sam two-class unpaired comparisons, and anova. a change in gene expression was considered significant if the p value was less than, the fold-change was at least, and the gene expression alteration occurred for all replicate experiments. we further expected each gene to be significantly differentially expressed for at least two of the three replicate spots for each experimental array set. based on these criteria, genes that were deemed significant by all four analytical methods were organized by cogs functional categories and compiled into a list that included genes that were up- or down-regulated when b. melitensis was grown to late-log phase, compared to stationary phase. a direct comparison of the signal intensity values of these genes indicated that the difference between log and stationary phases was specifically due to differential gene expression and not array spatial bias, as indicated in figure when the average gdna intensity values for these genes were plotted, the r value was. however, the r value for the same genes comparing the cy fluorescence values instead was extremely low. to confirm the microarray results, we randomly chose differentially expressed genes and performed qrt-pcr. based on qrt-pcr results, transcript levels of of these genes were altered greater than -fold and in the same direction as was determined by microarray analysis. two other genes were determined to be differentially expressed and in the same direction of microarray analysis, but the fold change was lower than no significant difference in the expression level of bmei was observed by qrt-pcr. b. melitensis m express different sets of genes in late-log and stationary phases of growth in f12k tissue culture medium of the genes significantly altered in b. melitensis during late-log phase, were up- and were down-regulated, compared to when the bacteria were allowed to reach stationary phase. the relative changes in gene expression ranged from a -fold induction of the glycerol-3-phosphate regulon repressor gene to a -fold down-regulation of the locus bmeii. as expected, the majority of gene expression changes were associated with growth and metabolism. among the up-regulated genes were those associated with dna replication, transcription and translation, nucleotide, amino acid, lipid and carbohydrate metabolism, energy production and conversion, membrane transport and cell envelope, biogenesis and outer membrane, while the down-regulated genes were distributed among several cogs. hierarchical cluster analysis using spotfire decisionsite was performed on normalized cy signal intensity values of the microarray data from the four log phase and four stationary phase samples. all four samples from the log phase of growth clustered together, apart from those collected at stationary phase. moreover, genes that clustered together were indeed differentially expressed between the two growth conditions. the higher number of genes up-regulated in late-log growth phase coincides with a more active metabolism of late-log phase cultures compared to those at stationary phase. in the following sections, we will focus our comments on those genes differentially expressed by microarray analysis that encode or are predicted to encode virulence factors, some of which may be involved in brucella:host interaction. protein-encoded genes which play a role in brucella invasiveness in non-phagocytic cells did not have differential expression between the most and the least invasive cultures currently, only three brucella gene products have been characterized as important for invasion in non-phagocytic cells. the b. abortus two-component regulatory system bvrr/bvrs encoded by bvrr/bvrs genes, regulates the structure of outer membrane components and plays a critical role in cell penetration and intracellular survival. this two-component system is highly conserved in the genus brucella, with chvi/chvg representing the b. melitensis homolog. in this study, neither of the two genes that encode this two-component system were differentially expressed between the most and the least invasive b. melitensis cultures. another brucella invasive-characterized gene product is sp a surface protein that enables b. suis to attach and penetrate non-phagocytic cells. the role of this gene has not been evaluated in b. melitensis, although a homolog is encoded by the ugpb gene present on the chromosome ii of the b. melitensis m genome. in this study, ugpb was not differentially expressed when global gene expression of b. melitensis cultures at late-log phase was compared to cultures at stationary growth phase. recently, a third gene product was reported to be involved in brucella internalization in non-phagocytic cells. in that study, a b. melitensis mutant with interruption in the bmei gene exhibited a marked decrease in its ability to invade hela cells at and h post-infection, suggesting the relevance of this gene in the brucella invasion process after h p.i. in this study, bmei was not found altered due to growth-phase. collectively, these results indicate that the higher invasiveness observed in b. melitensis cultures at late-log phase of growth under our experimental conditions was not due to the differential expression of these three characterized gene products. also, these results suggest that genes encoding these three products were not growth-phase regulated in b. melitensis under our experimental conditions. however, they might be transcribed at a time that we did not measure, they could be constitutively expressed and act in concert with other factors, or they could be expressed following epithelial cell contact. it is perhaps worth noting that only one of these three gene products has been effectively demonstrated to contribute to b. melitensis virulence, although after one hour post infection rather than the minutes used in this study. well-known b. melitensis virulence genes had different expression profiles in late-log phase of growth compared to stationary growth phase several genes whose products are known to be associated with brucella melitensis virulence, were differentially expressed between the most and the least invasive cultures. these included genes that encode t4ss proteins and the flagellar apparatus. the virb locus, for instances, encodes the type iv secretion system and plays a critical role in brucella virulence and intracellular multiplication. three genes encoding components for the virb operon, such as virb, virb and virb were up-regulated in b. melitensis cultures at late-log phase compared to stationary growth phase. pathogenic bacteria produce flagella to promote colonization and invasion of mucosa. brucellae are traditionally characterized as non-motile bacteria, yet the sequence of the b. melitensis genome contains three clusters of flagellar genes and their participation in establishing chronic brucellosis has been established. in our study, five genes such as flic, flif, flin, flha and flgd which encode parts of the flagellar apparatus or regulate its expression, were differentially expressed in late-log phase cultures compared to stationary phase cultures. previous studies reported scant influence of t4ss and flagella in the invasion process. thus, the highest penetration observed in late-log phase cultures was probably not due to the expression of these genes. several transcriptional regulator genes were differentially expressed in late-log phase compared to stationary growth phase transcriptional regulators control bacterial gene expression in response to specific signals. twenty-two genes encoding transcriptional regulators belonging to the arac, asnc, beti, deor, gntr, iclr, lysr, luxr, marr, merr, and rpir families were differentially expressed in late-log phase b. melitensis cultures compared to stationary phase cultures. some of these transcription factors are known to be involved in positive regulation of gene expression. others are involved in repression, while members of iclr and lysr families could be activators or repressors of gene expression. nevertheless, the contribution of these regulators and their targets to b. melitensis internalization epithelial cells has not been fully examined. the locus encoding the alternative sigma factor that allows brucella to survive under general stress situations was up-regulated in stationary phase cultures. the bmei locus that encodes a subunit of the other alternative sigma factor, which allows transcription of those genes involved in utilization of nitrogen and carbon sources and energy metabolism, was up-regulated in late-log phase cultures compared to stationary phase cultures. two-component transcriptional regulators are comprised of a cytoplasmic membrane-located sensor protein and a cytoplasmic response regulator protein. eight orfs encoding for two-component response regulators have been identified in the b. melitensis m genome. one of the signal transduction-encoded genes up-regulated in late-log phase cultures, was previously identified in b. melitensis attenuated mutants. the other is a central regulator of carbohydrate metabolism genes and also plays a role in virulence development of certain pathogens. although the molecular regulation of these response regulators in b. melitensis is currently unknown, understanding how vsr, hprk and others are regulated, could offer insight into b. melitensis virulence. identifying the target genes of these transcriptional regulators would significantly clarify the role of growth-phase in brucella physiology, metabolism and virulence regulation. almost all differentially expressed genes encoding cell envelope and outer membrane components were up-regulated in late-log phase cultures the ability of brucella to invade cells has been linked to its outer membrane properties, as well as to structures built within the cell envelope. twenty-six genes directly involved in cell envelope and outer membrane biogenesis were differentially expressed at late-log compared to stationary phase of growth. these included genes that encode outer membrane proteins, lipoproteins, lps, and peptidoglycan biosynthesis. the main cogs functional category of genes that were up-regulated in b. melitensis cultures at late-log compare to stationary phase of growth were orfs encoding membrane transport proteins. these included genes encoding transporters specific for amino acids, carbohydrates and uncharacterized transporters. collectively, these data indicate an active conversion of metabolites to components of the cell envelope structure at late-log phase, which given the importance of the om in virulence, could influence the initial brucella:host cell interactions, facilitating attachment and entry into host cells. discussion the molecular mechanisms involved in the initial interactions between brucella and epithelial cells have not been well characterized. previous studies have used hela cells as a model for studying adhesion and internalization of brucella spp. in non-professional phagocytic cells. these studies found that brucellae bind to cellular receptors containing sialic acid residues and induce their own uptake by a local rearrangement of the host cell cytoskeleton around the invading organisms. the ability of the bacteria to adhere to and penetrate eukaryotic cells is a well orchestrated process that requires several factors/gene products in order to be successful. to date, only a few brucella gene products involved in non-phagocytic cell invasion have been identified. this study was performed with the goal of better understanding initial molecular interactions between brucella and its host through the molecular analysis of growth phase-specific gene regulation. our initial experiment indicated that cultures of b. melitensis at late-log growth phase in cell culture medium were more invasive to non-phagocytic cells than cultures at mid-log and stationary growth phases. similar results have been observed for other invasive pathogens, such as salmonella spp. or yersinia enterocolitica. even with the high moi used, b. melitensis were internalized in lower numbers by epithelioid-like hela cells at min p.i. than reported in another study. the difference in invasion may have been influenced by the f12k cell culture medium used to growth the agent. b. melitensis reach stationary phase at a lower od in f12k cell culture medium than in rich bacterial culture medium or another cell culture medium. these results suggest that f12k medium apparently contains suboptimal nutrients for brucella development. even though, we grew b. melitensis in f12k medium and immediately added the bacteria to hela cells with the goal of reducing bacterial pre-infection manipulations, which had probably modified the original transcriptome of the cultures, since bacterial gene expression changes quickly in response to environmental modification. the relationship between growth phase and invasiveness is dependent upon the expression of bacterial virulence factors at different growth-phase. for instance, shigella flexneri is more invasive during the early phase of their exponential growth, because invasion proteins are secreted in higher amounts during this growth phase. salmonella enterica and legionella pneumophila have their secretion systems assembled and effector proteins properly stored in the cytoplasm only at the late exponential and stationary growth phases, respectively. in order to understand why our system evoked greater invasiveness in b. melitensis cultures at late-log phase in the first min p.i., we conducted a global gene expression detection study using cdna microarray technology. microarray analysis revealed that genes were significantly differentially expressed between the most and the least invasive cultures. as expected, the majority of the observed changes in gene expression were related to the bacterial response under the increased growth conditions in tissue culture media. for example, the up-regulation of genes associated with transcription and translation, nutrient metabolism, transport, and energy production and conversion all correspond to a more active metabolism of late-log phase cultures, compared to cultures at stationary phase. as was expected, several cell division- and dna synthesis-related genes were also up-regulated at late-log phase, when the bacterial population was still actively growing. alternatively, genes down-regulated in late-log phase were more heterogeneous in nature, demonstrating no predominant functional category. as expected, an increased expression of the locus bmei encoding the alternative sigma factor was observed in stationary phase cultures. sigma factor regulates the transcription of heat shock genes, which allow the bacteria to survive not only an abrupt increase in temperature, but also general stress situations, such as nutrient limitations during stationary growth phase. previous work identified a role in b. melitensis invasion of hela cells for the hypothetical protein encoded by bmei orf, which increases invasiveness only after h p.i.. that study clearly showed that the presence or absence of the gene transcript did not modify the ability of b. melitensis to invade hela cells during the first min p.i., i.e. the brucella-hela co-incubation time used in our study. under our experimental conditions, bmei was not found phase growth regulated. these data suggest that bmei may be transcribed after prolonged host cell contact, thereby facilitating the invasion process at later time points. further characterization of the regulation of this gene and its product is clearly warranted. in seeking to identify possible contributors to the increased invasiveness of b. melitensis at late-log phase, the conversion of metabolites to components that alter cell envelope structure were evaluated. altered outer membrane/cell wall topology would be expected to influence the initial bacteria:host cell interaction that may facilitate attachment and entry into host cells. in our study, a significant number of genes directly involved in cell envelope/outer membrane biogenesis were differentially expressed. one of the genes up-regulated at late-log growth phase was the locus bmei the product of this gene has not yet been characterized in b. melitensis; however, it has high homology to an immunogenic outer membrane protein, omp. omp is a haemin-binding protein, which binds to, and extracts iron from, the host. iron has been identified as a required element for epithelial invasion in microbial pathogens, and the expression of this locus, along with other iron-related genes in late-log phase cultures, may influence the internalization ability of brucellae. sp is another surface-exposed outer membrane protein with a critical role in brucella suis adherence to, and invasion of, non-phagocytic cells. the role of this protein, which is encoded by the ugpb gene present in the chromosome ii of b. melitensis m genome, was not previously described for b. melitensis adhesion to and/or penetration of epithelial cells. the transcript from the ugpb gene was not identified as differentially expressed in our cdna microarray analysis between the most and the least invasive cultures. therefore, under our experimental conditions, this omp seems not to be involved in the higher invasiveness of the late-log phase cultures. it is possible that the composition of the cell culture medium does not induce the expression of ugpb, or it is also possible that ugpb is constitutively expressed and/or act in concert with other factors. although genetic analysis reveals that ugpb may belong to an operon that encodes for a sn-glycerol-3-phosphate abc transporter, the experimental evidence does not support this hypothesis. a previous study showed that the product of ugpb in b. suis is indeed a surface-exposed protein with adhesion and invasion activity. in fact, in this study, three of the transcripts predicted to encode the transport system were highly up-regulated in late-log phase cultures, when compared to stationary phase cultures. in concordance with previous experimental evidence, our microarray data would support the finding of others that ugpb does not belong to an operon that encodes for a sn-glycerol-3-phosphate abc transporter. in addition, our results support growth-phase regulation of the sn-glycerol-3-phosphate abc transport system, which has been implicated in brucella pathogenesis. the ability of brucella to invade host cells is linked to its om properties. b. melitensis omp profile changes during culture growth, as gene expression is transcriptional regulated by environmental conditions. although the different expression pattern observed in gene-encoding cell envelope products between late-log and stationary phase cultures here may have been due to medium composition, it is more likely that they are due to growth phase regulation. in this study, we did not evaluate the role of the omp in internalization in epithelial cells and therefore their individual participation in increased invasiveness of late-log phase cultures could not be determined. only two differentially expressed genes encoding for o-chain and peptidoglycan layer biosynthesis from this study, were previously evaluated in brucella pathogenesis, although not in epithelial cells internalization. due to the importance that the cell envelope in initial host:pathogen interaction, the regulation and role of gene-encoding om products differentially expressed in this study should be addressed in future studies. rapid adaptive physiological response to multiple environmental and cellular signals in bacteria is mainly mediated by transcriptional regulators and two-component regulatory systems. prokaryotic genes putatively coding for transcriptional regulators are grouped in families based on sequence similarity and functional criteria. twenty-two transcripts, belonging to families of transcriptional regulators, were differentially expressed in our study. it was recently reported that b. melitensis mutants for of these transcriptional regulators were not attenuated after one-week of infection in mice. however, effects of these transcriptional regulators on internalization of b. melitensis by non-phagocytic cells have not been examined. their contribution to invasion therefore remains unknown. luxr is a well-known family of transcriptional activators that regulates various functions in microbes. there are two loci that encode transcripts belonging to this family of transcriptional regulators in the b. melitensis genome, and their expression is required for transcription of virulence factors such as virb operon and flagella. the transcriptional regulator vjbr was not differentially expressed in our study, but the other luxr homolog, was 221-fold up-regulated in the late-log phase of growth, compared to stationary phase cultures. the targets of blxr are currently unidentified, but regulatory effects on other transcriptional-regulatory proteins and proteins predicted to be involved in cell envelope biogenesis was observed. it may be possible that some of these gene products regulated by blxr positively influence b. melitensis invasion of hela cells. analysis of the invasive phenotype of a b. melitensis blxr deletion mutant in hela cells would be the first step in determining the importance of this transcriptional regulator during the initial host:brucella interaction, followed by the identification of the effector gene it regulates. bvrr/bvrs is a well characterized two-component regulatory system that controls the expression of genes essential for brucella abortus invasion to non-phagocytic cells. high level of identity is present between b. melitensis chvi/chvg and the b. abortus bvrr/bvrs proteins. in our study, no transcriptional change was observed in bmei2036/i orfs between the most and the least invasive cultures. likely, brucella maintain a basal expression level of the regulatory locus, as a change in the phosphorylation of the protein required for activity rather than transcription. twenty orfs dedicated to signal transduction were identified in b. melitensis genome. the importance of some of them in brucella virulence had been characterized lately, including blxr, vjbr, ftcr and bvrr/bvrs. however, their contribution to internalization in non-phagocytic cells is less known. recently, mutants with defective expression in two transcriptional regulators had an altered pattern in initial host:pathogen interaction due to surface modifications. future identification of the target genes of these regulators would clarify brucella physiology, metabolism and virulence regulation. several motility-related genes were more highly expressed at late-log phase compared to stationary phase, including kinesin-like protein, chemotaxis motd protein and genes related to flagellum apparatus synthesis and functions, e.g. flagellin itself. flagellin has been well-characterized as a contributor to bacterial virulence through chemotaxis and adhesion to and invasion of host cells. the extent to which flagellar machinery participates in the process of invasion seems to depend at least partly on the species of bacteria and/or the host cell type. for instance, flagellar-associated motility in salmonella is not required but accelerates invasion of caco- colonic epithelial cells, whereas the invasion of acanthamoeba astronyxis by burkholderia pseudomallei absolutely requires an intact flagellum apparatus. in the case of b. melitensis, a previous study demonstrated that expression of flagella is growth curve-dependent and required for persistent disease in a mouse model but not for invasion in cellular models. that study reported that a functional flagellum was assembled in early-log growth phase cultures but not at later time points. in our study, we did not analyze gene expression at early time points of the growth curve, but the results indicated that some flagellar genes were expressed more in late-log phase cultures as compared to stationary phase cultures. the differences in flagellum gene expression between the study of fretin et al. and ours could be attributed to evaluation of different steps of the process, different culture media used, or post-transcriptional regulatory mechanisms. the flagellar apparatus is built hierarchically under complex regulation. thirty-one flagellar genes distributed in three clusters on chromosome ii and along with three transcriptional regulators of flagellar system expression have been identified in b. melitensis. however, the order of flagella gene expression and the whole system regulation in brucellae has not been established. here, only five genes from two loci encoding different parts of the flagellar apparatus were differentially expressed in late-log phase cultures compared to stationary phase cultures. detection of expression of some but not all genes from an operon is not uncommon with microarray data, due to the inherent nature of microarrays that produces variation that leads to some false negatives. in a previous study, rambow-larsen et al. using a cdna microarray, also identified only of the flagellar genes, belonging to different flagellar loci and encoding for distinct parts of the flagellar apparatus, expressed under a putative quorum-sensing regulator blxr. similarly, microarray detected changes in expression of only some of the genes of the flagellar operon in salmonella enterica serovar typhimurium, which is transcribed with a polycistronic message, despite a 10-fold difference in some genes of each operon. two different functions, motility and protein secretion have been ascribed to flagella, but these roles have yet to be demonstrated in brucellae. we were not able to evaluate the role of b. melitensis flagellar gene expression in invasion under our experimental conditions, but undoubtedly, the presence of flagellar machinery and other adhesion/motility factors at late-log phase, and their exact contribution to the brucella invasion process warrant further studies. the virb operon has been reported to be essential for intracellular survival and multiplication of brucella, but its role in adherence and internalization is contradictory. in our study, three genes from the operon were up-regulated in late-log growth phase cultures compared to the stationary phase of growth. virb is transcribed as an operon, with no secondary promoters. it is maximally expressed in b. melitensis at the early exponential phase of the growth curve, and its expression decays as the bacteria reach the stationary phase. however, the half-lives of the individual segments of the virb transcript are not known. under our experimental conditions, it is possible that virb was expressed earlier in the growth curve, and the different rate of transcript degradation allowed the detection of expression of some genes of the operon in late-log phase but not in stationary phase cultures. distinct half-lives of the individual segments of the polycistronic mrna, results in the differential expression of each gene in the operon. this phenomenon has been well characterized in other bacteria, and is worthy to additional evaluation of b. melitensis virb operon. in addition, and similar to mention for flagellar genes, microarray could detect expression of some but not all genes from an operon, due to the inherent nature of the technique. further, our analysis method was particularly stringent in order to greatly reduce false positives at the risk of additional false negatives. thus, other genes in the virb operon were increased in expression such as virb virb virb virb and virb although not statistically significant because of the stringency of our statistical analysis. finally, genes with uncharacterized function that were differentially expressed at late-log phase compared with the stationary phase also deserve some special consideration. this group of "hidden genes" represents 22% of the differentially expressed genes identified in this study, and it may contain some of the heretofore unknown virulence factors utilized for b. melitensis to invade and infect the host, as was previously suggested. conversely, brucella internalization should not be disregarded as a product of synergistic action among several gene products in non-phagocytic cells.  our study reveals that b. melitensis grown in cell culture medium at late-log phase are more invasive in non-phagocytic cells than cultures grown at mid-log or stationary growth phases. cdna microarrays provide informative differential transcriptional profiles of the most and the least invasive b. melitensis cultures. we consider these data a platform for conducting further studies on the brucella:host initial interaction. since the roles of the majority of differentially expressed genes in this study are not well defined in brucella pathogenesis, future studies on brucella virulence can now be specifically focused to more precisely delineate the roles of candidate genes identified in this study. 
b. melitensis at the late logarithmic phase of growth are more invasive to hela cells than mid-logarithmic or stationary growth phases. microarray analysis of b. melitensis gene expression identified up- and down-regulated genes in late-log growth phase compared to the stationary growth phase. are the etiological agents of brucellosis, a zoonotic infectious disease that causes abortion in animals and chronic debilitating illness in humans. as expected, the majority of up-regulated genes in late-log phase cultures were those associated with growth, including dna replication, transcription, translation, intermediate metabolism, energy production and conversion, membrane transport, and biogenesis of the cell envelope and outer membrane; while the down-regulated genes were distributed among several functional categories. in this study, we characterized changes in genome-wide transcript abundance of the most and the least invasive growth phases of b. melitensis cultures to hela cells, as a preliminary approach for identifying candidate pathogen genes involved in invasion of epithelial cells. further characterization of some genes found differentially expressed in the most invasive culture will likely bring new insights into the initial molecular interactions between brucella and its host. natural brucella infections occur primarily through an incompletely defined mechanism of adhesion to and penetration of mucosal epithelium.
 brucella spp. are the etiological agents of brucellosis, a zoonotic infectious disease that causes abortion in animals and chronic debilitating illness in humans. natural brucella infections occur primarily through an incompletely defined mechanism of adhesion to and penetration of mucosal epithelium. in this study, we characterized changes in genome-wide transcript abundance of the most and the least invasive growth phases of b. melitensis cultures to hela cells, as a preliminary approach for identifying candidate pathogen genes involved in invasion of epithelial cells.  b. melitensis at the late logarithmic phase of growth are more invasive to hela cells than mid-logarithmic or stationary growth phases. microarray analysis of b. melitensis gene expression identified up- and down-regulated genes in late-log growth phase compared to the stationary growth phase. as expected, the majority of up-regulated genes in late-log phase cultures were those associated with growth, including dna replication, transcription, translation, intermediate metabolism, energy production and conversion, membrane transport, and biogenesis of the cell envelope and outer membrane; while the down-regulated genes were distributed among several functional categories.  this brucella global expression profile study provides novel information on growth phase-specific gene expression. further characterization of some genes found differentially expressed in the most invasive culture will likely bring new insights into the initial molecular interactions between brucella and its host. 