BACKGROUND
stathmin  is a ubiquitous microtubule  destabilizing protein linked to cancer and cell health: stathmin is highly over-expressed in leukemias and a number of other cancers, where its expression level often correlates with cancer stage progression and prognosis for survival  <cit> .
stathmin is the founding member of a family of mt destabilizers known as the stathmin family of proteins, which includes scg <dig> , sclip , and rb <dig>   <cit> , each expressed from separate genes.
each of the four stathmin family proteins shares a homologous tubulin binding site that functions as both a mt destabilizer and tubulin heterodimer sequestering protein.
stathmin is expressed in a wide range of tissues and is present as a soluble cytosolic protein  <cit> , while scg <dig> , and rb <dig>  are neuron-specific homologues of stathmin localized to membranes in developing  and mature  nerve cells.
though it is well established that stathmin regulates mts, many have suggested alternative functions for stathmin.
stathmin has been called a cell survival factor because its level of overexpression correlates with cancer stage progression, invasion, and metastasis for many cancer types .
for example, knockdown of stathmin protein by sirna  <cit> , shrna  <cit> , or ribozymes  <cit>  leads to apoptosis of several cancer cell lines in culture.
stathmin has also been linked to cell motility and metastasis.
overexpression of stathmin stimulates motility of both gn- <dig> neurons  <cit>  and ht- <dig> fibrosarcoma cells  <cit> .
in fibrosarcoma cells, stathmin activity is regulated by p27kip <dig>  <cit> .
ng and coauthors  <cit>  have also proposed a role for stathmin in mouse embryonic fibroblast  cell migration, although their results indicate that stathmin inhibits, rather than promotes, cell migration.
in this study, stathmin activity was regulated by the transcription factor stat <dig>  <cit> .
the ability of stathmin to positively or negatively regulate motility may be context-specific, where stathmin promotes motility in 3d matrices, but not in 2d  <cit> .
although stathmin has a significant function in regulating the mt cytoskeleton, surprisingly stmn <dig> knockout  mice develop normally except for some minor age-onset axonopathies associated with stmn1-/-  <cit>  and a lack of learned or innate fear response  <cit> .
it is not known why stathmin is dispensable for normal development but required for survival of many cancer cell lines.
for example, it is not known whether compensatory changes occur in mice lacking stmn <dig> to permit normal development.
here we performed an unbiased screen for transcriptome level changes associated with genetic knockout of stmn <dig>  microarray studies of mefs, genetically designated stmn1+/+ or stmn1-/-, were used to compare global transcriptome level changes due to deletion of the stmn <dig> gene.
examination of genes encoding proteins either related to stathmin, or regulated by stathmin, revealed differential expression of the stathmin family proteins and tubulin isotypes.
additional classification of differential regulation was performed using gene ontology  and kegg pathway analyses, and demonstrated that expression of genes functioning in cell adhesion and extracellular matrix make up the majority of up and down regulated genes, supporting a role for stmn <dig> in cell motility.
methods
isolation of mefs
c57bl/ <dig> stmn1+/- male and female mice  were mated.
<dig>  days to  <dig>  days post coitus, pregnant females were sacrificed by cervical dislocation, and fibroblasts were isolated as described by tessarollo  <cit> .
fibroblast cells from individual embryos were plated and allowed to grow for 1– <dig> days prior to storage of aliquots in liquid nitrogen.
genotyping of mefs was performed as described by liedtke  <cit> .
briefly, dna was isolated from embryonic tissue using isoamyl alcohol/phenol extraction methods.
samples were amplified using pcr to identify embryos with intact stmn <dig> intron iii or the neomycin cassette used to disrupt the stmn <dig> gene  <cit> .
the pcr jump start® redtaq kit  was used for amplification; each sample contained  <dig>  mm mgcl <dig>  deoxynucleotide triphosphates , primers , and taq polymerase .
for the wild-type amplification,  <dig> temperature cycles were performed .
the following primer set was used as a probe for wild-type alleles: forward sense primer 5'-3' , corresponding to a region of intron iii deleted in the mutant allele; and a reverse anti-sense primer 5'-3'  also missing in the mutant, corresponding to a region of intron iii of the stmn <dig> gene that yields an amplification product  <dig> bp long.
a second set of reactions were performed using mutant-allele specific primers that anneal to the neomycin insert  <cit> ; forward sense primer 5'-3'  and a reverse anti-sense primer 5'-3' , corresponding to a region of the inserted neomycin cassette of the mutant allele that yields an amplification product  <dig> bp long.
cell culture
all cells were cultured at 37°c in a humidified atmosphere of 5% carbon dioxide.
cells were grown in dmem  supplemented with  <dig>  g/l d-glucose, l-glutamine ,  <dig> mm sodium bicarbonate, 1× antibiotic/antimycotic ,  <dig> mm sodium pyruvate , and 10% fetal bovine serum  .
cell cultures were passaged at ~80% confluency, about every 3– <dig> days.
cultures were discarded after passage  <dig>
rna isolation
total rna was isolated from cultured mef cell lines using trizol® reagent  following manufacturer's instructions, followed by dnase treatment.
rna quality was assessed by a260/a <dig> ratio  and further validated by the agilent  <dig> bioanalyzer.
the same mrna was used in both microarray and qrt-pcr experiments.
microarray
stmn+/+ and stmn1-/- mrna were fluorescently labeled with either cy <dig> or cy <dig> using the agilent low rna input fluorescent linear amplification kit .
hybridized two-color samples were prepared using agilent whole mouse genome  <dig> ×  <dig> k oligo microarray kit format  according to agilent instructions and using agilent reagents.
each array includes  <dig>  biological features with a total of  <dig>  probes of 60-mer controls and gene probes.
arrays were run as color-swap duplicates and scanned with an agilent microarray scanner system, which generated the tiff images of low and high intensity scans utilized by agilent feature extraction software .
feature extraction processing of fluorescent data corrected signals for background noise, foreground intensities, positive and negative spot controls, background subtraction, and signal normalization.
results were collected into tab delimited text files for each of the four experimental arrays for analysis with agilent technologies software genespring gx .
data were processed in genespring gx  by first filtering on flags for features either present or marginal in any  <dig> of  <dig> arrays, followed by error filtering to a coefficient of variance < 50%, then statistical analysis using the t-test against zero with a false discovery rate of <  <dig> .
the resultant list comprised  <dig> genes,  <dig> of which were up-regulated and  <dig> down-regulated.
further analyses were done using fold change cut-offs.
qrt-pcr
rna was used to synthesize cdna with superscript™ iii first-strand synthesis system for quantitative rt-pcr .
rna was isolated from at least two cell lines from each genotype for cdna preparation and qrt-pcr.
all stathmin family proteins and isotypes of mouse α- and β-tubulin were probed using primers listed in .
oligonucleotide probes were designed using primer express® software   and cdna sequences from the mouse genome  <cit> .
for each probe, we used at least three dilutions of cdna, each dilution run in triplicate.
commercial 96-well plate arrays for the focal adhesion and extracellular matrix pathways were also assayed for each mef genotype .
two plates were run for each genotype.
message levels were quantified for each genotype using oligonucleotide probes designed to amplify cdna fragments containing target gene exons.
results were normalized to the signal from gapdh for designed probes, or to the signal from five internal standards  for the sabiosciences arrays.
amplification was achieved using power sybr® green pcr master mix  and applied biosystems  <dig> real time pcr system with sds v <dig>  software.
rt-pcr amplicon specificity was checked by electrophoresis of rt-pcr products on 2% agarose gels .
following manufacturer instructions , the threshold for determining ct values was set to log scale  <dig>  and internal normalization of genotype results to gapdh was first calculated .
the ct values for gapdh were stable and consistent across genotypes.
one measure of the consistent gapdh amplification is the small standard deviations for the gapdh ct among replicates.
the standard deviation for gapdh was ≤  <dig>  –  <dig> % of the mean ct value.
Δct values were not calculated for target mrna samples that exceeded  <dig> cycles prior to crossing the threshold.
the relative abundance of target mrna between genotypes was calculated as  <dig> raised to the negative of the difference in Δct values
all statistical analyses were done using excel .
rt-pcr data are presented as fold changes relative to mef stmn1+/+ samples.
standard errors  were computed using Δct values transformed to se values of fold change using the formula:
significance was determined by anova of Δct values  <cit> .
amplification efficiency was determined to be 100% for gapdh for each genotype per manufacturer's instructions , and it was assumed that efficiency of all target genes was also 100% for our statistical analysis.
RESULTS
microarray analysis of gene regulation in stmn1-/- mefs
changes to stathmin protein levels regulate microtubule polymer in mefs  <cit>  and cancer cell metastatic potential and invasiveness  <cit> .
to probe whether stmn <dig> gene knockout resulted in a change to the mef gene expression profile, we identified the transcriptome differences between stmn1+/+ and stmn1-/- mefs using agilent 2-color microarrays run in quadruplicate with cy3/cy <dig> dye color-swap.
the initial  <dig> statistically significant probes were further parsed for differential gene expression fold change of ≥  <dig>  to encompass as many potentially relevant gene-level alterations as possible.
the inclusive list of genes totaled  <dig> probes, of which  <dig> were up-regulated and  <dig> were down-regulated in stmn1-/- compared to the stmn1+/+ controls.
of the  <dig> probes,  <dig> have known functions.
of these,  <dig> were up-regulated and  <dig> were down-regulated as shown in the heatmap of figure 1a.
more stringent analysis was performed on the list of  <dig> genes at a fold change threshold of ≥  <dig> .
at this level, comparison of stmn1+/+ and stmn1-/- cell lines resulted in expression changes for  <dig> genes , of which  <dig> were up-regulated and  <dig> were down-regulated .
expression of stathmin family and tubulin isotypes measured by microarray
stathmin is the founding member of a family of proteins, all of which contain a homologous tubulin binding site  <cit> .
the microarray did not show complete loss of stmn <dig> in the stmn-/- mefs because the probe design for stmn <dig> is flawed.
the agilent 60-mer was found to match both stmn <dig> on chromosome  <dig> and a region of chromosome  <dig>  a source of immunodominant mhc-associated peptides   <cit> .
however, by both qrt-pcr and western blotting, stmn <dig> mrna and protein are absent in stmn-/- mefs  <cit> .
we next looked for changes to stathmin family mrnas as a consequence of stmn <dig> knockout and found that stmn <dig> decreased by  <dig>  fold and stmn <dig> increased by  <dig>  fold .
stmn <dig> was not found to be differentially regulated on the microarray.
tubulins are expressed from  <dig> α-tubulin genes and  <dig> β-tubulin genes.
the cohort of tubulin isotypes present can assemble into mts with different dynamic properties in vitro  <cit> .
therefore, we examined whether stmn <dig> genotype influences the tubulin isotype expression profile of each cell line.
no significant changes in expression of α-tubulin isotypes were detected.
in contrast, small, but significant  changes in the expression levels of β-tubulin isotypes were observed in the microarray , with β-tubulin isotype classes iia/b increased  <dig> / <dig>  fold, class-iii increased  <dig>  fold, class ivb decreased  <dig>  fold, and class-v decreased  <dig>  fold.
magnitude and direction of gene expression changes found using microarray analysis compared to qrt-pcr with designed primers.
only those tubulin isotypes whose expression was significantly changed, as measured by microarray with fdr <  <dig> , are shown compared to qrt-pcr .
data are given as fold change in stmn-/- mefs relative to stmn+/+ mef.
confirmation of stathmin family and tubulin isotype expression by qrt-pcr
changes to stathmin family mrna and tubulin isotype levels were verified by qrt-pcr using designed primers .
for stmn1-/- cells relative to stmn+/+ cells, upregulation of stmn <dig>  was seen by both microarray and qrt-pcr.
in contrast, stmn <dig>  showed down regulation by microarray, but upregulation by qrt-pcr .
the reason for the difference in expression pattern by the two methods is not known but suggests that changes in stmn <dig> expression are either not significant or not consistently altered compared to the level in wt cells.
we also used qrt-pcr to examine the message levels for stathmin and related proteins in mefs heterozygous for stathmin .
the stmn1+/- cells showed a  <dig>  fold increase in stmn <dig>  similar to the increased expression in stmn-/- mefs measured by qrt-pcr.
for stmn <dig>  mrna was increased  <dig> fold in stmn+/- mefs compared to that in wild-type cells, consistent with the up regulation observed in stmn-/- cells.
the stmn <dig> mrna was not present at quantifiable levels in rt-pcr for any of the three genotypes, consistent with the microarray for wild-type and knockout genotypes.
similar changes in stathmin family expression were also observed in an additional cell line for each genotype.
protein products other than stathmin were undetectable with available antibodies, suggesting that protein products of stmn <dig> and stmn <dig> are present at low levels.
qrt-pcr was also used to confirm expression of tubulin isotypes for each mouse tubulin isotype.
mefs expressed five α- and seven β-tubulin isotypes at quantifiable levels by qrt-pcr .
for the α-tubulin isotypes, stmn-/- mefs expressed slightly higher levels of tuba1c and tuba <dig> mrnas compared to wild-type mefs  .
by microarray, the probes for tuba1c and tuba <dig> did not pass the stringency filters applied to the raw data set and were not flagged as significantly changed in expression.
examination of the unfiltered microarray data showed opposite regulation to that measured by rt-pcr for tuba1c and tuba <dig>  the discrepancies between rt-pcr and microarray may reflect difficulty detecting small changes in expression level.
as shown in table  <dig>  microarray and rt-pcr detected similar changes in expression of five β tubulins isotypes.
stmn-/- mefs showed higher expression of isotype classes iia, iib and iii, while isotype classes ivb and v were reduced in stmn1-/- mefs relative to wild-type cells .
by rt-pcr, we also detected upregulation of β-tubulin class i and down regulation of class iva in stmn-/- mefs compared to wild-type cells .
the microarray data for class i and iva β-tubulins did not pass the applied data filters.
although they did not pass the data filters, class iva tubulin showed reduced expression by both rt-pcr and microarray, but changes to class i tubulin were in opposite directions by microarray and rt-pcr.
therefore, with a few exceptions, the data from rt-pcr for tubulin isotypes generally confirmed the microarray results.
microarray results for stmn  <dig> and stmn <dig> were also confirmed by qrt-pcr, and in general, both stmn family and tubulin isotype expression confirmed the microarray results .
correlation between microarray and qrt-pcr can depend on a large number of factors  <cit> .
qualitatively, the correlations we found between microarray and qrt-pcr are consistent with those described by others  <cit> .
in addition, we also examined tubulin isotype expression in the heterozygous stmn1+/- cell line and observed changes in β-tubulin mrna levels similar to that measured in knockout mefs.
we also found that the stmn1+/- cell lines expressed lower mrna levels of most α-isotypes, but showed higher expression of tuba4a mrna compared to wild-type cells.
tubulin mrna levels do not always correlate with measured protein levels  <cit> .
therefore, we used western blots to estimate relative protein levels for those isotypes where specific antibodies are available .
the anti-α tubulin antibody dm1a recognizes the c-terminal amino acid sequence found in tuba1a and tuba1b.
immunoblots probed with dm1a indicated that the combined protein level of tuba1a and tuba1b was reduced in heterozygous cells, which is consistent with the decreased mrna level for each isotype relative to wild-type cells.
the combined protein level of tuba1a and b was higher in stmn1-/- cells compared to wild-type cells, while the corresponding mrna levels were slightly reduced in the stmn1-/- cells .
the protein levels for β-tubulin isotype classes ii and iii decreased with loss of stmn <dig>  while class iv β-tubulin increased in the stmn1-/- line.
for each mef cell line, the protein levels of β-tubulin isotypes were inversely related to their corresponding mrna level.
for example, in stmn1-/- cells, both β-tubulin classes ii and iii showed increased mrna levels and decreased protein levels relative to wild-type cells.
for class iv β-tubulins, mrnas decreased while protein level increased in stmn1-/- cells relative to wild-type cells.
stmn <dig> genotype and expression of mt regulators and motors
the microarray data were then examined for expression of mt regulators and motor proteins, the kinesin family and dyneins, which transport cargo along mts.
as shown in table  <dig>  we found that stmn <dig> knockout resulted in small fold-change differences to a number of mt stabilizing and destabilizing proteins, without a clear pattern of up-regulation of expression for proteins that might compensate for stathmin or for down-regulation of proteins that normally work in opposition to stathmin.
most notably, we did not see a change in the expression of map <dig>  a protein previously shown to function antagonistically to stathmin  <cit> .
in general, expression of genes for other proteins that destabilize mts, including mt severing proteins, were down-regulated in stmn-/- mefs relative to stmn1+/+ cells, making it unlikely that these mt destabilizers substitute for stmn <dig>  surprisingly, the largest group of genes whose expression in stmn1-/- mefs differed from that in wt mefs was the group encoding motor proteins.
we found differential expression of  <dig> kinesin genes and  <dig> genes encoding dynein light or intermediate chains.
two kinesins, kif 21b and kif26b showed >  <dig> fold up-regulation in stmn1-/- mefs.
these two kinesins have not been well studied and it is unclear what the consequences of their up regulation may be.
magnitude and direction of gene expression changes found using microarray analysis.
gene names are in m. musculus format, fold changes are from microarray data.
stmn <dig> loss impact categorized by kegg pathways and gene ontology
the analyses described above involved directed searches for genes encoding proteins functioning as part of the mt cytoskeleton.
to expand our analysis and take advantage of the non-biased nature of microarrays, we performed pathway analysis on those genes showing a ≥  <dig>  fold change using the online resource kegg pathways  <cit>  for mus musculus.
top pathway hits include focal adhesion , extra-cellular matrix  receptor interactions , cell communication , and regulation of actin cytoskeleton  .
an alternate tool to kegg pathway analysis is the newly developed kegg spider  <cit> .
kegg spider was used to better characterize the interconnectedness of the differentially expressed genes.
the output from kegg spider includes enriched gene ontology  categories based on an interconnection network algorithm that also provides p-value significance measures.
go characterization of  <dig> genes with fold changes ≥  <dig>  was performed within kegg spider with a p-value cutoff of  <dig> , yielding  <dig> individual categories .
the most highly represented go terms included cell adhesion and extracellular matrix, and the sum of all cell motility-related genes comprised about 50% of all go terms represented .
given the changes to genes implicated in cell adhesion and cell motility, we used a pathway specific  rt-pcr kit to further confirm the microarray results.
relative transcript levels in stmn1+/+ and stmn1-/- mefs were measured for  <dig> focal adhesion and ecm mrnas.
data that met a fold change threshold of ≥  <dig>  and p-value of <  <dig>  in the pcr array were considered significant for comparison to microarray data .
there were  <dig> differentially regulated genes that matched in the direction  of fold change, while  <dig> genes were oppositely regulated compared to the microarray .
the majority  of mrnas for cell adhesion and ecm were up or down regulated in the same direction  for both microarray and qrt-pcr, providing general support for the microarray results.
the most extreme differentially expressed entities in stmn1-/- compared to stmn1+/+ measured by the qrt-pcr array were catenin β  <dig>  which was down-regulated  <dig>  fold in support of enhanced migration, and p-selectin  which was up-regulated by  <dig>  fold.
the qrt-pcr measurements were made from a commercial 96-well plate format.
data for qrt-pcr is shown as the mean of duplicate assays.
the top  <dig> entries show mrnas with matching directionality of expression fold change between rt-pcr and microarray.
the last  <dig> entries list mrnas in which expression fold change measured by rt-pcr is in the opposite direction of the microarray data.
discussion
stathmin is an integral protein involved in control of mt polymer level during interphase of the cell cycle  <cit> .
many cancers express high levels of stathmin  and increased stathmin level is correlated with reduced patient survival  <cit> .
cancer cells over-expressing stathmin protein display invasive and metastatic behaviors  <cit> , and alternative stathmin binding partners have recently been revealed  <cit> .
these roles for stathmin in regulation of the mt cytoskeleton and cell motility, and stathmin's as yet undefined role in cancer, led us to use microarray technologies as an unbiased screen for transcriptome level changes due to genetic knockout of the stmn <dig> gene to provide clues to stathmin function beyond the mt cytoskeleton.
loss of stmn <dig> is correlated with expression changes to genes encoding proteins functioning in cell motility pathways
several recent studies have indicated a role for stathmin in regulation of cell motility  <cit> .
here, we show that genetic knockout of stmn <dig> has transcriptome level impact on genes broadly classified in cell motility pathways, including integrins, ecm components, and cell adhesion , which comprised 50% of all go categories derived from our microarray data.
additionally, qrt-pcr using a pathway specific array gave strong supporting evidence of cell adhesion and ecm protein transcript changes as a result of stmn <dig> loss .
of the  <dig> differentially regulated genes that match results from the  <dig> microarray replicates,  <dig> integrins are represented, thus linking changes in stmn <dig> expression to change in focal adhesion proteins that also link to the actin cytoskeleton.
although we cannot predict whether the measured changes in gene expression would promote or inhibit cell migration, our data point to a role for stathmin in regulating the expression level of a wide range of proteins required for cell adhesion and locomotion.
these data indicate that stmn <dig> knockout has consequences beyond simply regulating assembly of the mt cytoskeleton by regulating the expression level of a cohort of proteins required for cell motility.
whether these expression changes are a downstream consequence of altered mt assembly dynamics is not known.
loss of stmn <dig> impacts the mt system
cells differing in stmn <dig> genotype expressed different mrna levels for the stathmin related protein stmn <dig> , with stmn-/- mefs showing upregulation of stmn <dig>  by either microarray or rt-pcr analysis.
we cannot detect stmn <dig> expression in mefs with antibodies currently available, indicating that stmn <dig> protein is likely present at a low level in these cells.
we note that the expression of other mt destabilizers are reduced in stmn1-/- mefs, making it unlikely that cells compensate for loss of stmn <dig> through expression of stmn <dig> or other mt destabilizers .
several previous reports have also noted that stathmin family proteins are upregulated in stmn <dig> knockout mice.
liedtke et al.
<cit>  demonstrated that stmn <dig> expression is higher in the brains of aged knock-out mice and yoshie et al.
<cit>  showed increased expression of stmn <dig> in the peri-implantation uteri of knockout mice.
it is possible that stmn <dig> may compensate for loss of stmn <dig> in these tissues.
we did not detect changes in stmn <dig> expression in stmn1-/- mefs, either by microarray or qrt-pcr, indicating that stmn <dig> does not compensate for stmn <dig> in these cells.
in addition to changing expression levels of stmn <dig>  stmn1-/- mefs also expressed different levels of tubulin isotypes compared to wild-type cells .
tubulin isotype mrna levels do not necessarily correspond to protein levels  <cit> , as we also found in mefs.
for example, stmn1-/- mefs expressed higher levels of β-tubulin class iii mrna but lower levels of the corresponding protein .
alterations to β-tubulin protein isotype ratio, specifically the levels of classes ii, iii, and iv proteins, can impact rates of mt growth and shortening in vitro  <cit>  although this is not always observed in cells  <cit> .
changes to β-tubulin isotype protein expression have also been linked to paclitaxel resistance in cancers  <cit>  via increased mt dynamic instability.
in mefs, the amount of β-tubulin class ii and iii proteins fell dramatically with loss of stmn <dig> while class iv protein increased.
class ii and iv form more stable mts, while class iii forms more dynamic mts in vitro  <cit> , and therefore, mefs do not display a clear pattern of tubulin isotype up or down regulation to compensate for loss of stathmin.
CONCLUSIONS
our microarray and qrt-pcr results showed that deletion of both copies of the stmn <dig> gene resulted in changes to the expression level of many genes, including a large cohort of genes related to cell adhesion and motility.
it is possible that stathmin depletion has additional, or different, consequences in those cells where it is expressed at higher levels than that in embryonic fibroblasts.
overall, our data support a role for stmn <dig> in cell migration and cross-talk between the mt and actin components of the cytoskeleton  <cit> , as suggested by baldassarre et al.
<cit>  and belletti et al.
<cit> .
authors' contributions
dnr participated in experimental design, performed experiments and analyses; assembled tables and figures and drafted the manuscript.
lc participated in experimental design, assembly of tables and figures and helped to draft the manuscript.
all authors read and approved the final version of the manuscript.
authors' information
dnr is a graduate student in the chemistry department, lehigh university.
the manuscript describes a portion of her ph.d. research project in lc's laboratory in biological sciences, lehigh university.
supplementary material
additional file 1
primers designed for stathmin family and tubulin isotype qrt-pcr.
all primers were designed as described in methods.
bold nucleotides represent exon junctions in mrnas.
click here for file
additional file 2
genes differentially expressed in stmn1-/- mefs compared to stmn1+/+ mefs.
differentially expressed genes listed with accession numbers, gene symbols, and average fold change over quadruplicate microarrays with color-swap.
click here for file
additional file 3
tubulin isotype mouse versus human nomenclature and antibody specificities.
tubulin isotypes for α- and β-tubulins are listed with accession numbers, current gene name, old gene name, new classification scheme , c-terminal sequences of each, and antibodies that recognize those sequences if known.
click here for file
additional file 4
kegg pathways represented by genes in additional table  <dig>  data entered into the kegg pathways website  yielded a number of pathway hits.
shown here are the top  <dig> pathways and the genes represented in each.
click here for file
additional file 5
go terms represented by genes in additional file  <dig>  data entered into kegg spider  used to generate the pie chart in figure  <dig>
click here for file
acknowledgements
the authors would like to thank dr. jutta marzillier, genomics facility at lehigh university, for assistance with microarray; dr. ping-shi wu, mathematics department, lehigh university, for help with statistical analyses; dr. bob skibbens for helpful discussions; and dr. gleb shumyatsky, rutgers university, for heterozygous mice.
supported by nih gm <dig> to lc.

 stathmin is a ubiquitous microtubule destabilizing protein linked to cancer and cell health: stathmin is highly over-expressed in leukemias and a number of other cancers, where its expression level often correlates with cancer stage progression and prognosis for survival. stathmin is the founding member of a family of mt destabilizers known as the stathmin family of proteins, which includes scg, sclip, and rb, each expressed from separate genes. each of the four stathmin family proteins shares a homologous tubulin binding site that functions as both a mt destabilizer and tubulin heterodimer sequestering protein. stathmin is expressed in a wide range of tissues and is present as a soluble cytosolic protein, while scg, and rb are neuron-specific homologues of stathmin localized to membranes in developing and mature nerve cells. though it is well established that stathmin regulates mts, many have suggested alternative functions for stathmin. stathmin has been called a cell survival factor because its level of overexpression correlates with cancer stage progression, invasion, and metastasis for many cancer types. for example, knockdown of stathmin protein by sirna, shrna, or ribozymes leads to apoptosis of several cancer cell lines in culture. stathmin has also been linked to cell motility and metastasis. overexpression of stathmin stimulates motility of both gn- neurons and ht- fibrosarcoma cells. in fibrosarcoma cells, stathmin activity is regulated by p27kip. ng and coauthors have also proposed a role for stathmin in mouse embryonic fibroblast cell migration, although their results indicate that stathmin inhibits, rather than promotes, cell migration. in this study, stathmin activity was regulated by the transcription factor stat. the ability of stathmin to positively or negatively regulate motility may be context-specific, where stathmin promotes motility in 3d matrices, but not in 2d. although stathmin has a significant function in regulating the mt cytoskeleton, surprisingly stmn knockout mice develop normally except for some minor age-onset axonopathies associated with stmn1-/- and a lack of learned or innate fear response. it is not known why stathmin is dispensable for normal development but required for survival of many cancer cell lines. for example, it is not known whether compensatory changes occur in mice lacking stmn to permit normal development. here we performed an unbiased screen for transcriptome level changes associated with genetic knockout of stmn microarray studies of mefs, genetically designated stmn1+/+ or stmn1-/-, were used to compare global transcriptome level changes due to deletion of the stmn gene. examination of genes encoding proteins either related to stathmin, or regulated by stathmin, revealed differential expression of the stathmin family proteins and tubulin isotypes. additional classification of differential regulation was performed using gene ontology and kegg pathway analyses, and demonstrated that expression of genes functioning in cell adhesion and extracellular matrix make up the majority of up and down regulated genes, supporting a role for stmn in cell motility. methods isolation of mefs c57bl/ stmn1+/- male and female mice were mated. days to days post coitus, pregnant females were sacrificed by cervical dislocation, and fibroblasts were isolated as described by tessarollo. fibroblast cells from individual embryos were plated and allowed to grow for 1– days prior to storage of aliquots in liquid nitrogen. genotyping of mefs was performed as described by liedtke. briefly, dna was isolated from embryonic tissue using isoamyl alcohol/phenol extraction methods. samples were amplified using pcr to identify embryos with intact stmn intron iii or the neomycin cassette used to disrupt the stmn gene. the pcr jump start® redtaq kit was used for amplification; each sample contained mm mgcl deoxynucleotide triphosphates, primers, and taq polymerase. for the wild-type amplification, temperature cycles were performed. the following primer set was used as a probe for wild-type alleles: forward sense primer 5'-3', corresponding to a region of intron iii deleted in the mutant allele; and a reverse anti-sense primer 5'-3' also missing in the mutant, corresponding to a region of intron iii of the stmn gene that yields an amplification product bp long. a second set of reactions were performed using mutant-allele specific primers that anneal to the neomycin insert; forward sense primer 5'-3' and a reverse anti-sense primer 5'-3', corresponding to a region of the inserted neomycin cassette of the mutant allele that yields an amplification product bp long. cell culture all cells were cultured at 37°c in a humidified atmosphere of 5% carbon dioxide. cells were grown in dmem supplemented with g/l d-glucose, l-glutamine, mm sodium bicarbonate, 1× antibiotic/antimycotic, mm sodium pyruvate, and 10% fetal bovine serum. cell cultures were passaged at ~80% confluency, about every 3– days. cultures were discarded after passage rna isolation total rna was isolated from cultured mef cell lines using trizol® reagent following manufacturer's instructions, followed by dnase treatment. rna quality was assessed by a260/a ratio and further validated by the agilent bioanalyzer. the same mrna was used in both microarray and qrt-pcr experiments. microarray stmn+/+ and stmn1-/- mrna were fluorescently labeled with either cy or cy using the agilent low rna input fluorescent linear amplification kit. hybridized two-color samples were prepared using agilent whole mouse genome × k oligo microarray kit format according to agilent instructions and using agilent reagents. each array includes biological features with a total of probes of 60-mer controls and gene probes. arrays were run as color-swap duplicates and scanned with an agilent microarray scanner system, which generated the tiff images of low and high intensity scans utilized by agilent feature extraction software. feature extraction processing of fluorescent data corrected signals for background noise, foreground intensities, positive and negative spot controls, background subtraction, and signal normalization. results were collected into tab delimited text files for each of the four experimental arrays for analysis with agilent technologies software genespring gx. data were processed in genespring gx by first filtering on flags for features either present or marginal in any of arrays, followed by error filtering to a coefficient of variance < 50%, then statistical analysis using the t-test against zero with a false discovery rate of <. the resultant list comprised genes, of which were up-regulated and down-regulated. further analyses were done using fold change cut-offs. qrt-pcr rna was used to synthesize cdna with superscript™ iii first-strand synthesis system for quantitative rt-pcr. rna was isolated from at least two cell lines from each genotype for cdna preparation and qrt-pcr. all stathmin family proteins and isotypes of mouse α- and β-tubulin were probed using primers listed in. oligonucleotide probes were designed using primer express® software and cdna sequences from the mouse genome. for each probe, we used at least three dilutions of cdna, each dilution run in triplicate. commercial 96-well plate arrays for the focal adhesion and extracellular matrix pathways were also assayed for each mef genotype. two plates were run for each genotype. message levels were quantified for each genotype using oligonucleotide probes designed to amplify cdna fragments containing target gene exons. results were normalized to the signal from gapdh for designed probes, or to the signal from five internal standards for the sabiosciences arrays. amplification was achieved using power sybr® green pcr master mix and applied biosystems real time pcr system with sds v software. rt-pcr amplicon specificity was checked by electrophoresis of rt-pcr products on 2% agarose gels. following manufacturer instructions, the threshold for determining ct values was set to log scale and internal normalization of genotype results to gapdh was first calculated. the ct values for gapdh were stable and consistent across genotypes. one measure of the consistent gapdh amplification is the small standard deviations for the gapdh ct among replicates. the standard deviation for gapdh was ≤ – % of the mean ct value. Δct values were not calculated for target mrna samples that exceeded cycles prior to crossing the threshold. the relative abundance of target mrna between genotypes was calculated as raised to the negative of the difference in Δct values all statistical analyses were done using excel. rt-pcr data are presented as fold changes relative to mef stmn1+/+ samples. standard errors were computed using Δct values transformed to se values of fold change using the formula: significance was determined by anova of Δct values. amplification efficiency was determined to be 100% for gapdh for each genotype per manufacturer's instructions, and it was assumed that efficiency of all target genes was also 100% for our statistical analysis.  microarray analysis of gene regulation in stmn1-/- mefs changes to stathmin protein levels regulate microtubule polymer in mefs and cancer cell metastatic potential and invasiveness. to probe whether stmn gene knockout resulted in a change to the mef gene expression profile, we identified the transcriptome differences between stmn1+/+ and stmn1-/- mefs using agilent 2-color microarrays run in quadruplicate with cy3/cy dye color-swap. the initial statistically significant probes were further parsed for differential gene expression fold change of ≥ to encompass as many potentially relevant gene-level alterations as possible. the inclusive list of genes totaled probes, of which were up-regulated and were down-regulated in stmn1-/- compared to the stmn1+/+ controls. of the probes, have known functions. of these, were up-regulated and were down-regulated as shown in the heatmap of figure 1a. more stringent analysis was performed on the list of genes at a fold change threshold of ≥. at this level, comparison of stmn1+/+ and stmn1-/- cell lines resulted in expression changes for genes, of which were up-regulated and were down-regulated. expression of stathmin family and tubulin isotypes measured by microarray stathmin is the founding member of a family of proteins, all of which contain a homologous tubulin binding site. the microarray did not show complete loss of stmn in the stmn-/- mefs because the probe design for stmn is flawed. the agilent 60-mer was found to match both stmn on chromosome and a region of chromosome a source of immunodominant mhc-associated peptides. however, by both qrt-pcr and western blotting, stmn mrna and protein are absent in stmn-/- mefs. we next looked for changes to stathmin family mrnas as a consequence of stmn knockout and found that stmn decreased by fold and stmn increased by fold. stmn was not found to be differentially regulated on the microarray. tubulins are expressed from α-tubulin genes and β-tubulin genes. the cohort of tubulin isotypes present can assemble into mts with different dynamic properties in vitro. therefore, we examined whether stmn genotype influences the tubulin isotype expression profile of each cell line. no significant changes in expression of α-tubulin isotypes were detected. in contrast, small, but significant changes in the expression levels of β-tubulin isotypes were observed in the microarray, with β-tubulin isotype classes iia/b increased / fold, class-iii increased fold, class ivb decreased fold, and class-v decreased fold. magnitude and direction of gene expression changes found using microarray analysis compared to qrt-pcr with designed primers. only those tubulin isotypes whose expression was significantly changed, as measured by microarray with fdr <, are shown compared to qrt-pcr. data are given as fold change in stmn-/- mefs relative to stmn+/+ mef. confirmation of stathmin family and tubulin isotype expression by qrt-pcr changes to stathmin family mrna and tubulin isotype levels were verified by qrt-pcr using designed primers. for stmn1-/- cells relative to stmn+/+ cells, upregulation of stmn was seen by both microarray and qrt-pcr. in contrast, stmn showed down regulation by microarray, but upregulation by qrt-pcr. the reason for the difference in expression pattern by the two methods is not known but suggests that changes in stmn expression are either not significant or not consistently altered compared to the level in wt cells. we also used qrt-pcr to examine the message levels for stathmin and related proteins in mefs heterozygous for stathmin. the stmn1+/- cells showed a fold increase in stmn similar to the increased expression in stmn-/- mefs measured by qrt-pcr. for stmn mrna was increased fold in stmn+/- mefs compared to that in wild-type cells, consistent with the up regulation observed in stmn-/- cells. the stmn mrna was not present at quantifiable levels in rt-pcr for any of the three genotypes, consistent with the microarray for wild-type and knockout genotypes. similar changes in stathmin family expression were also observed in an additional cell line for each genotype. protein products other than stathmin were undetectable with available antibodies, suggesting that protein products of stmn and stmn are present at low levels. qrt-pcr was also used to confirm expression of tubulin isotypes for each mouse tubulin isotype. mefs expressed five α- and seven β-tubulin isotypes at quantifiable levels by qrt-pcr. for the α-tubulin isotypes, stmn-/- mefs expressed slightly higher levels of tuba1c and tuba mrnas compared to wild-type mefs. by microarray, the probes for tuba1c and tuba did not pass the stringency filters applied to the raw data set and were not flagged as significantly changed in expression. examination of the unfiltered microarray data showed opposite regulation to that measured by rt-pcr for tuba1c and tuba the discrepancies between rt-pcr and microarray may reflect difficulty detecting small changes in expression level. as shown in table microarray and rt-pcr detected similar changes in expression of five β tubulins isotypes. stmn-/- mefs showed higher expression of isotype classes iia, iib and iii, while isotype classes ivb and v were reduced in stmn1-/- mefs relative to wild-type cells. by rt-pcr, we also detected upregulation of β-tubulin class i and down regulation of class iva in stmn-/- mefs compared to wild-type cells. the microarray data for class i and iva β-tubulins did not pass the applied data filters. although they did not pass the data filters, class iva tubulin showed reduced expression by both rt-pcr and microarray, but changes to class i tubulin were in opposite directions by microarray and rt-pcr. therefore, with a few exceptions, the data from rt-pcr for tubulin isotypes generally confirmed the microarray results. microarray results for stmn and stmn were also confirmed by qrt-pcr, and in general, both stmn family and tubulin isotype expression confirmed the microarray results. correlation between microarray and qrt-pcr can depend on a large number of factors. qualitatively, the correlations we found between microarray and qrt-pcr are consistent with those described by others. in addition, we also examined tubulin isotype expression in the heterozygous stmn1+/- cell line and observed changes in β-tubulin mrna levels similar to that measured in knockout mefs. we also found that the stmn1+/- cell lines expressed lower mrna levels of most α-isotypes, but showed higher expression of tuba4a mrna compared to wild-type cells. tubulin mrna levels do not always correlate with measured protein levels. therefore, we used western blots to estimate relative protein levels for those isotypes where specific antibodies are available. the anti-α tubulin antibody dm1a recognizes the c-terminal amino acid sequence found in tuba1a and tuba1b. immunoblots probed with dm1a indicated that the combined protein level of tuba1a and tuba1b was reduced in heterozygous cells, which is consistent with the decreased mrna level for each isotype relative to wild-type cells. the combined protein level of tuba1a and b was higher in stmn1-/- cells compared to wild-type cells, while the corresponding mrna levels were slightly reduced in the stmn1-/- cells. the protein levels for β-tubulin isotype classes ii and iii decreased with loss of stmn while class iv β-tubulin increased in the stmn1-/- line. for each mef cell line, the protein levels of β-tubulin isotypes were inversely related to their corresponding mrna level. for example, in stmn1-/- cells, both β-tubulin classes ii and iii showed increased mrna levels and decreased protein levels relative to wild-type cells. for class iv β-tubulins, mrnas decreased while protein level increased in stmn1-/- cells relative to wild-type cells. stmn genotype and expression of mt regulators and motors the microarray data were then examined for expression of mt regulators and motor proteins, the kinesin family and dyneins, which transport cargo along mts. as shown in table we found that stmn knockout resulted in small fold-change differences to a number of mt stabilizing and destabilizing proteins, without a clear pattern of up-regulation of expression for proteins that might compensate for stathmin or for down-regulation of proteins that normally work in opposition to stathmin. most notably, we did not see a change in the expression of map a protein previously shown to function antagonistically to stathmin. in general, expression of genes for other proteins that destabilize mts, including mt severing proteins, were down-regulated in stmn-/- mefs relative to stmn1+/+ cells, making it unlikely that these mt destabilizers substitute for stmn surprisingly, the largest group of genes whose expression in stmn1-/- mefs differed from that in wt mefs was the group encoding motor proteins. we found differential expression of kinesin genes and genes encoding dynein light or intermediate chains. two kinesins, kif 21b and kif26b showed > fold up-regulation in stmn1-/- mefs. these two kinesins have not been well studied and it is unclear what the consequences of their up regulation may be. magnitude and direction of gene expression changes found using microarray analysis. gene names are in m. musculus format, fold changes are from microarray data. stmn loss impact categorized by kegg pathways and gene ontology the analyses described above involved directed searches for genes encoding proteins functioning as part of the mt cytoskeleton. to expand our analysis and take advantage of the non-biased nature of microarrays, we performed pathway analysis on those genes showing a ≥ fold change using the online resource kegg pathways for mus musculus. top pathway hits include focal adhesion, extra-cellular matrix receptor interactions, cell communication, and regulation of actin cytoskeleton. an alternate tool to kegg pathway analysis is the newly developed kegg spider. kegg spider was used to better characterize the interconnectedness of the differentially expressed genes. the output from kegg spider includes enriched gene ontology categories based on an interconnection network algorithm that also provides p-value significance measures. go characterization of genes with fold changes ≥ was performed within kegg spider with a p-value cutoff of, yielding individual categories. the most highly represented go terms included cell adhesion and extracellular matrix, and the sum of all cell motility-related genes comprised about 50% of all go terms represented. given the changes to genes implicated in cell adhesion and cell motility, we used a pathway specific rt-pcr kit to further confirm the microarray results. relative transcript levels in stmn1+/+ and stmn1-/- mefs were measured for focal adhesion and ecm mrnas. data that met a fold change threshold of ≥ and p-value of < in the pcr array were considered significant for comparison to microarray data. there were differentially regulated genes that matched in the direction of fold change, while genes were oppositely regulated compared to the microarray. the majority of mrnas for cell adhesion and ecm were up or down regulated in the same direction for both microarray and qrt-pcr, providing general support for the microarray results. the most extreme differentially expressed entities in stmn1-/- compared to stmn1+/+ measured by the qrt-pcr array were catenin β which was down-regulated fold in support of enhanced migration, and p-selectin which was up-regulated by fold. the qrt-pcr measurements were made from a commercial 96-well plate format. data for qrt-pcr is shown as the mean of duplicate assays. the top entries show mrnas with matching directionality of expression fold change between rt-pcr and microarray. the last entries list mrnas in which expression fold change measured by rt-pcr is in the opposite direction of the microarray data. discussion stathmin is an integral protein involved in control of mt polymer level during interphase of the cell cycle. many cancers express high levels of stathmin and increased stathmin level is correlated with reduced patient survival. cancer cells over-expressing stathmin protein display invasive and metastatic behaviors, and alternative stathmin binding partners have recently been revealed. these roles for stathmin in regulation of the mt cytoskeleton and cell motility, and stathmin's as yet undefined role in cancer, led us to use microarray technologies as an unbiased screen for transcriptome level changes due to genetic knockout of the stmn gene to provide clues to stathmin function beyond the mt cytoskeleton. loss of stmn is correlated with expression changes to genes encoding proteins functioning in cell motility pathways several recent studies have indicated a role for stathmin in regulation of cell motility. here, we show that genetic knockout of stmn has transcriptome level impact on genes broadly classified in cell motility pathways, including integrins, ecm components, and cell adhesion, which comprised 50% of all go categories derived from our microarray data. additionally, qrt-pcr using a pathway specific array gave strong supporting evidence of cell adhesion and ecm protein transcript changes as a result of stmn loss. of the differentially regulated genes that match results from the microarray replicates, integrins are represented, thus linking changes in stmn expression to change in focal adhesion proteins that also link to the actin cytoskeleton. although we cannot predict whether the measured changes in gene expression would promote or inhibit cell migration, our data point to a role for stathmin in regulating the expression level of a wide range of proteins required for cell adhesion and locomotion. these data indicate that stmn knockout has consequences beyond simply regulating assembly of the mt cytoskeleton by regulating the expression level of a cohort of proteins required for cell motility. whether these expression changes are a downstream consequence of altered mt assembly dynamics is not known. loss of stmn impacts the mt system cells differing in stmn genotype expressed different mrna levels for the stathmin related protein stmn, with stmn-/- mefs showing upregulation of stmn by either microarray or rt-pcr analysis. we cannot detect stmn expression in mefs with antibodies currently available, indicating that stmn protein is likely present at a low level in these cells. we note that the expression of other mt destabilizers are reduced in stmn1-/- mefs, making it unlikely that cells compensate for loss of stmn through expression of stmn or other mt destabilizers. several previous reports have also noted that stathmin family proteins are upregulated in stmn knockout mice. liedtke et al. demonstrated that stmn expression is higher in the brains of aged knock-out mice and yoshie et al. showed increased expression of stmn in the peri-implantation uteri of knockout mice. it is possible that stmn may compensate for loss of stmn in these tissues. we did not detect changes in stmn expression in stmn1-/- mefs, either by microarray or qrt-pcr, indicating that stmn does not compensate for stmn in these cells. in addition to changing expression levels of stmn stmn1-/- mefs also expressed different levels of tubulin isotypes compared to wild-type cells. tubulin isotype mrna levels do not necessarily correspond to protein levels, as we also found in mefs. for example, stmn1-/- mefs expressed higher levels of β-tubulin class iii mrna but lower levels of the corresponding protein. alterations to β-tubulin protein isotype ratio, specifically the levels of classes ii, iii, and iv proteins, can impact rates of mt growth and shortening in vitro although this is not always observed in cells. changes to β-tubulin isotype protein expression have also been linked to paclitaxel resistance in cancers via increased mt dynamic instability. in mefs, the amount of β-tubulin class ii and iii proteins fell dramatically with loss of stmn while class iv protein increased. class ii and iv form more stable mts, while class iii forms more dynamic mts in vitro, and therefore, mefs do not display a clear pattern of tubulin isotype up or down regulation to compensate for loss of stathmin.  our microarray and qrt-pcr results showed that deletion of both copies of the stmn gene resulted in changes to the expression level of many genes, including a large cohort of genes related to cell adhesion and motility. it is possible that stathmin depletion has additional, or different, consequences in those cells where it is expressed at higher levels than that in embryonic fibroblasts. overall, our data support a role for stmn in cell migration and cross-talk between the mt and actin components of the cytoskeleton, as suggested by baldassarre et al. and belletti et al.. authors' contributions dnr participated in experimental design, performed experiments and analyses; assembled tables and figures and drafted the manuscript. lc participated in experimental design, assembly of tables and figures and helped to draft the manuscript. all authors read and approved the final version of the manuscript. authors' information dnr is a graduate student in the chemistry department, lehigh university. the manuscript describes a portion of her ph.d. research project in lc's laboratory in biological sciences, lehigh university. supplementary material additional file 1 primers designed for stathmin family and tubulin isotype qrt-pcr. all primers were designed as described in methods. bold nucleotides represent exon junctions in mrnas. click here for file additional file 2 genes differentially expressed in stmn1-/- mefs compared to stmn1+/+ mefs. differentially expressed genes listed with accession numbers, gene symbols, and average fold change over quadruplicate microarrays with color-swap. click here for file additional file 3 tubulin isotype mouse versus human nomenclature and antibody specificities. tubulin isotypes for α- and β-tubulins are listed with accession numbers, current gene name, old gene name, new classification scheme, c-terminal sequences of each, and antibodies that recognize those sequences if known. click here for file additional file 4 kegg pathways represented by genes in additional table data entered into the kegg pathways website yielded a number of pathway hits. shown here are the top pathways and the genes represented in each. click here for file additional file 5 go terms represented by genes in additional file data entered into kegg spider used to generate the pie chart in figure click here for file acknowledgements the authors would like to thank dr. jutta marzillier, genomics facility at lehigh university, for assistance with microarray; dr. ping-shi wu, mathematics department, lehigh university, for help with statistical analyses; dr. bob skibbens for helpful discussions; and dr. gleb shumyatsky, rutgers university, for heterozygous mice. supported by nih gm to lc.
to expand our understanding of stathmin function in normal cells, we compared gene expression profiles, measured by microarray and qrt-pcr, of mouse embryo fibroblasts isolated from stmn1+/+ and stmn1-/- mice to determine the transcriptome level changes present in the genetic knock-out of stathmin. microarray analysis of stmn loss at a fold change threshold of ≥ revealed expression changes for genes, of which were up-regulated and were down-regulated. stathmin over-expression has been correlated with cancer stage progression, while stathmin depletion leads to death of some cancer cell lines in culture. these results support models indicating a role for stathmin in regulating cell locomotion, but also suggest that this functional activity may involve changes to the cohort of proteins expressed in the cell, rather than as a direct consequence of stathmin-dependent regulation of the microtubule cytoskeleton. microarray data and qrt-pcr analysis of mrna expression demonstrated changes in the message levels for stmn encoding rb a protein related to stathmin, and in alterations to many tubulin isotype mrnas. kegg pathway analysis of the microarray data indicated changes to cell motility-related genes, and qrt-pcr plates specific for focal adhesion and ecm proteins generally confirmed the microarray data. approximately 50% of genes up or down regulated in stmn1-/- mouse embryo fibroblasts function broadly in cell adhesion and motility. several microtubule assembly regulators and motors were also differentially regulated in stmn1-/- cells, but these changes should not compensate for loss of stathmin. several stathmin binding partners, in addition to tubulin, have been shown to affect cell motility in culture.
 stathmin protein functions to regulate assembly of the microtubule cytoskeleton by destabilizing microtubule polymers. stathmin over-expression has been correlated with cancer stage progression, while stathmin depletion leads to death of some cancer cell lines in culture. in contrast, stathmin-null mice are viable with minor axonopathies and loss of innate fear response. several stathmin binding partners, in addition to tubulin, have been shown to affect cell motility in culture. to expand our understanding of stathmin function in normal cells, we compared gene expression profiles, measured by microarray and qrt-pcr, of mouse embryo fibroblasts isolated from stmn1+/+ and stmn1-/- mice to determine the transcriptome level changes present in the genetic knock-out of stathmin.  microarray analysis of stmn loss at a fold change threshold of ≥ revealed expression changes for genes, of which were up-regulated and were down-regulated. microarray data and qrt-pcr analysis of mrna expression demonstrated changes in the message levels for stmn encoding rb a protein related to stathmin, and in alterations to many tubulin isotype mrnas. kegg pathway analysis of the microarray data indicated changes to cell motility-related genes, and qrt-pcr plates specific for focal adhesion and ecm proteins generally confirmed the microarray data. several microtubule assembly regulators and motors were also differentially regulated in stmn1-/- cells, but these changes should not compensate for loss of stathmin.  approximately 50% of genes up or down regulated in stmn1-/- mouse embryo fibroblasts function broadly in cell adhesion and motility. these results support models indicating a role for stathmin in regulating cell locomotion, but also suggest that this functional activity may involve changes to the cohort of proteins expressed in the cell, rather than as a direct consequence of stathmin-dependent regulation of the microtubule cytoskeleton. 