BACKGROUND
gastrointestinal stromal tumors  are the most common mesenchymal tumors of the gastrointestinal tract and are characterized by the expression of the kit receptor  and to a lesser extent of pdgfra , representing two closely related receptor tyrosine kinases   <cit> .
the majority of gist shows oncogenic mutations either in kit or pdgfra  <cit> .
mainly, mutations in exon  <dig> or  <dig> of the kit gene or in exon  <dig> of pdgfra lead to ligand independent, constitutive activation of the kinase function  <cit> .
about 60% of all gist carry an exon  <dig> mutation of kit which encodes the juxtamembrane domain of the receptor possessing an autoinhibitory function  <cit> .
less common mutations in pdgfra  are detected in gist that often display gastric location and epithelioid morphology  <cit> .
in a minority of cases  no mutations in the known kit or pdgfra hot spots are detected although these tumors express the kit protein.
this subgroup is called wild type gist  and comprises tumors in pediatric patients, in patients affected by the carney triad, neurofibromatosis type  <dig>  associated gist and a subset of sporadic adult gist  <cit> .
the pathogenetic mechanisms underlying wt-gist are poorly understood and there is limited benefit of imatinib therapy in these patients  <cit> .
therefore the identification of additional genetic factors contributing to the pathogenesis of gist may help to find new concepts of individualized therapy.
recently, the braf mutation p.v600e was found in 4-13% of wt-gist  <cit> .
for another subgroup of wt-gist including pediatric tumors, a strong igf1r expression combined partly with gene amplification was described  <cit> .
two other alternative rtk probably involved in the pathogenesis of gist are axl and met.
both kinases have been shown to be upregulated in gist resistant to treatment  <cit> .
axl is a member of the ufo/axl subfamily and activates the same signaling pathway as kit.
the tyrosine kinase domain of met is mutated in sporadic papillary renal carcinomas.
some mutations in the met gene are located in codons homologous to those in kit and it is suggested that these missense mutations lead to constitutive activation of the met protein  <cit> .
to develop additional therapy approaches it would be of interest to know whether this rtk also plays a role in wt-gist.
besides kit and pdgfra, csf1r , flt <dig>  and pdgfrb  belong to the same family of type iii rtk.
these five tyrosine kinases show a homologous structure and a comparable function in activation, proliferation and suppressing apoptosis  <cit> .
aberrant expression and mutations in either csf1r, flt <dig> and pdgfrb or their ligands have been described in several malignant diseases.
overexpression of csf1r is found in epithelial tumors such as breast and ovarian cancer  <cit> .
the translocation t of its ligand csf <dig> leads to the development of the tenosynovial giant cell tumor  <cit> .
aberrantly expressed flt <dig> is observed at high levels in a spectrum of hematologic malignancies  <cit> .
additionally, in aml an internal tandem duplication in the transmembrane domain of flt <dig> was identified which leads to constitutive activation of its kinase domain.
it seems that this mutation is not present or very rare in gist  <cit> .
pdgfrb is overexpressed in malignant peripheral nerve sheath tumors  and chordomas  <cit> .
to analyse expression profiles of certain tumors for research and diagnostic purposes, qpcr  is frequently applied because of its reproducibility and high sensitivity  <cit> .
this method is based on the normalization of the target gene expression on stably expressed internal reference genes.
a major challenge is the application of suitable reference genes which have to be tested and verified under defined experimental conditions  <cit> .
ideal reference genes have to be non-regulated, stable and not affected by biological or experimental conditions.
the target gene is amplified together with the reference gene in order to minimize experimental variability concerning reverse transcription enzymatic efficiencies, pcr efficiency, amount of starting material and differences between human tissues.
the reference gene and the target gene should have very robust and stable expression profiles to ensure accurate normalization and interpretation of results.
the most stable expressed gene from a set of genes can be identified by genorm  <cit> , a software program which additionally provides the number of genes required to calculate a robust normalization factor based on the geometric mean of these genes.
typical reference genes regulate basic and ubiquitous cellular functions and are responsible for the cellular maintenance, e.g.
gapdh or ß-actin.
however, these commonly used reference genes vary considerably in different tissue types or under different experimental conditions  <cit> .
there is no standard reference gene for all kinds of tissue types.
to our best knowledge only few studies investigated the alteration of stability of reference genes in different mesenchymal tumour entities.
in the present study we attempted to identify suitable reference genes in gastrointestinal stromal tumors by using a set of sixteen reference genes which are currently applied in qpcr procedures.
furthermore we evaluated the gene expression of kit as well as the alternative rtk flt <dig>  csf1-r, pdgfrb, axl and met in mutated and non-mutated gastrointestinal stromal tumors by qpcr using the identified reference genes.
the study was focussed on wt-gist compared to samples with mutations in kit exon  <dig> and  <dig> and pdgfra exon  <dig> in order to evaluate whether overexpression of these alternative rtk might contribute to the pathogenesis of wt-gist.
methods
samples
a total of  <dig> samples were included into this study.
all specimens were obtained in the years  <dig> and  <dig> under approved ethical protocols and with informed consent from each patient.
all samples were fixed in neutral-buffered formalin prior to paraffin embedding.
<dig> samples from normal tissue  as control group and  <dig> gist representing different mutational subgroups were evalutated.
<dig> samples of wt-gist,  <dig> samples of wt-gist associated with neurofibromatosis type  <dig>   <dig> samples with exon  <dig> mutation in kit,  <dig> samples with exon  <dig> mutation in kit and  <dig> samples with exon  <dig> mutation in pdgfra).
sequence analysis of kit  and pdgfra  was carried out as described earlier  <cit> .
additionally, in all samples the wild type status of the braf gene was ascertained.
all gist samples were stained immunohistochemically for cd <dig>  four normal tissues from the gastrointestinal region were available as fresh frozen and formalin-fixed, paraffin-embedded samples.
rna extraction and cdna synthesis
prior to rna extraction, paraffin-embedded tissues were cut into  <dig> μm sections and mounted on glass-slides.
six slides of each specimen were used for rna extraction.
the sections were deparaffinized by extracting twice in xylene for  <dig> min at room temperature.
rehydration was done in 100% ethanol, 90% ethanol, 80% ethanol and 70% ethanol made with depc-water for  <dig> min each.
tumor areas previously marked on a h&e slide were scraped from the sections with a sterile scalpel.
sections were transferred into a sterile  <dig>  ml tube.
extraction and purification of rna was done using the rneasy ffpe kit  according to the manufacturer's recommendations.
after tumor localization by h&e staining fresh frozen material was rasped into  <dig> μm thick pieces and rna purification was carried out using the rneasy kit .
both rneasy kits contain a step of dnase treatment.
finally rna from fresh and from fixed tissue was eluted in water.
the quantification was done spectrophotometrically .
<dig> ng of rna from each sample was reverse-transcribed using a random-hexamer primer and avian myoblastosis virus reverse transcriptase  according to the manufacturer's protocol .
cdna from those four fresh frozen and corresponding formalin-fixed, paraffin-embedded  control samples was generated in only one reverse transcription and loaded on the reference low density arrays .
further cdna samples were generated from the four ffpe control tissues, pooled within one patient and loaded onto a second independent reference low density array.
identification of reference genes
the reference genes used here were preselected because of their constitutive, non-regulated stable expression over a wide spectrum of tissues.
but nevertheless the preselected reference genes are not suitable for every kind of tissue and therefore need to be analyzed prior to use in a certain study.
the detection of suitable reference genes was carried out using taqman low density arrays .
these arrays are prefabricated 384-well cards where gene-specific primer and probe sets are spotted in small reaction chambers during manufacturing.
the cards have  <dig> separate loading ports leading into  <dig> wells each.
in this study, arrays with triplicates of  <dig> putative reference genes were used, so  <dig> different samples could be analyzed.
the list of assays is given in table  <dig>  for each sample,  <dig> ng cdna was mixed with 2× taqman universal pcr master mix .
<dig> μl of this mixture was loaded into each port and distributed into the reaction chambers by centrifugation.
the card was sealed and the quantitative pcr  was performed on an abi prism ht  <dig>  sequence detection system.
after pre-incubation for  <dig> min at 50°c and  <dig> min at 95°c, the pcr reaction was performed .
the fluorescent signal was measured in each cycle.
analyses of gene expression by qpcr
qpcr analysis was performed using the assays-on-demand products  listed in table  <dig>  these gene-specific qpcr assays consist of a pair of unlabeled pcr primers and a fam labeled specific probe.
according to the manufacturer of these assays, probe and primer sets that would amplify pseudogenes are excluded in the process of development.
reactions were carried out in a reaction volume containing  <dig> μl pcr master mix ,  <dig>  μl forward and reverse primer mix,  <dig> ng cdna ad  <dig> μl a.dest.
triplicate reactions were carried out for each transcript.
control reactions were performed using a minus rt preparation and a sample with a.dest instead of rna.
pcr conditions were the same as for the taqman low density arrays.
statistical analyses
the genorm applet for microsoft excel was used to determine the most stable genes among the sixteen candidate reference genes.
raw cq values were converted into relative quantities for analysis with genorm, where the highest relative quantity for each gene is set to  <dig>  the program selects from a panel of candidate reference genes the two most stable genes or a combination of multiple stable genes for normalization.
the gene expression stability  value is based on the combined estimate of intra- and intergroup expression variations of the genes studied and takes the pcr efficiency into account.
the limited m-value is  <dig> .
the most stable genes are stepwise selected from the investigated gene panel to estimate how many reference genes should be used.
the normalization factors define the optimal number of reference genes required for a precised normalization design.
the analysis of the expression data of tumor samples compared to control samples was performed with the rest software   <cit> .
the software normalizes the measured cq-values of the target genes with those of the reference genes and compares the expression data of tumor and control samples by considering the pcr efficiency and the mean crossing point deviation.
RESULTS
rna quality
rna quality of the samples was inspected on a 1% agarose gel.
the concentration and purity of the rna was characterized by the mean a260/ <dig> ratio and was on average  <dig>  for fresh frozen as well as for formalin fixed, paraffin-embedded samples and reflected pure and protein-free rna.
expression variability of cdna synthesis
four normal control ffpe tissue samples from the gastric and bowel wall were used to determine suitable reference genes.
one tlda was done with cdna from only one reverse transcription for each sample.
as it is assumed that cdna synthesis from ffpe tissue varies strongly in efficiency, pooled cdna synthesis samples from the same ffpe tissue were used for a second tlda.
the pooled cdna preparation showed an essentially lower variance in contrast to the cdna samples from only one reverse transcription .
the median of the m-values ranges from  <dig>  to  <dig> , thus the variability of the technical replicates is low.
within a single specimen the range is also very small.
the median of the cdna samples from a single reverse transcription ranges from  <dig>  to  <dig> .
however, each cdna sample from one reverse non-pooled transcription showed a high variance of the values .
therefore the variance of the samples can be lowered by cdna pooling.
the median of the cdna samples from one reverse transcription from fresh frozen tissue ranged from  <dig>  to  <dig>  and can be compared to the median of the pooled cdna from the ffpe tissue.
within a single specimen the range is higher than in the pooled cdna samples, but smaller than in the cdna from only one reverse transcription from the ffpe samples.
identification of suitable reference genes for normalization
all  <dig> putative reference genes had a high expression stability and the 'm'-value  was clearly below the 'm'-cutoff-value of  <dig> .
by stepwise exclusion of genes, the expression stability value 'm' is calculated  and the expression stability rises with the exclusion of further genes.
samples generated from only one reverse transcription, showed polr2a, tfrc, rplpo and gapdh as the most stable genes .
the average stability value was between  <dig>  and  <dig> .
the range of the values was  <dig> - <dig> .
the most stable genes for pooled cdna samples were pgk <dig>  ppia, rplpo and ipo <dig> with a median of  <dig> - <dig>  .
the range of the values was lower and varied between  <dig>  and  <dig> .
all in all the gene expression variability of the pooled cdna samples is much lower compared to the single reverse transcription cdna synthesis.
the gene ubc is the most stable one in fresh frozen tissue with the smallest value range , but rplpo and ppia belong also to the most stable genes in fresh frozen tissue with a median of  <dig> - <dig>  and the range of the values was  <dig> - <dig>  .
tbp was not expressed in the tissue.
by combining the lowest variability values of fixed and fresh tissue, the genes polr2a, ppia, rplpo and tfrc were detected .
hence these four genes were used for further analysis with the gist samples.
expression profiles of kit, flt <dig>  csf1-r, pdgfrb, axl and met
six patient groups  were used to analyze the gene expression of kit, flt <dig>  csf1-r, pdgfrb, axl and met.
here, the rest analysis with polr2a is shown exemplarily.
the rest analysis of the target genes kit, csf1r, flt <dig>  pdgfrb, axl and met with the reference genes rplpo, ppia and trfc showed the same significant results for the differential expression or at least the same trend in mutated and non-mutated gist .
as suspected, we found a significant overexpression of kit in exon  <dig> and exon  <dig> mutated gist in comparison with normal tissue.
also a significantly lowered expression of pdgfrb in both groups compared to normal tissue was shown.
the same effect was observed in pdgfra exon  <dig> mutated gist compared to normal tissue .
expression analysis of kit, csf1r, flt <dig>  pdgfrb, axl and met was carried out in six gist cohorts .
significant variances are in bold.
the expression data were calculated using the rest software with polr2a as reference gene.
nt: normal tissue; p <  <dig> ; +/- indicates up- or downregulation
on closer inspection of normal tissue and wt-gist a clear overexpression of kit in the tissue of wild type tumors was shown.
the same effect, but even stronger, could be detected by comparing nf1-associated wt-gist with normal tissue.
concerning kit expression, the tumors without mutation behave like kit mutated gist .
this was in accordance with immunohistochemical staining .
the other target genes, flt <dig>  csf1r, pdgfrb, axl and met showed lower expression levels in the tumor tissue compared to normal tissue except axl in nf <dig> associated gist.
met was even significantly reduced in wt-gist, kit exon  <dig>   <dig> and pdgfra exon  <dig> mutated gist.
no changes in expression levels were detected comparing normal tissue and wt-gist with kit exon  <dig> mutated gist, indicating that both groups possess similar expression levels for the analyzed receptors.
none of the groups showed significant expression alterations for csf1r and flt <dig>  an exception is the group of kit exon  <dig> mutated gist.
a lower expression of csf1r, flt <dig>  pdgfrb and met compared to normal tissue was identified in contrast to kit, which showed a significant overexpression.
in summary, overexpression of kit was detected in each gist subgroup compared to the corresponding normal tissue by using the preassigned reference gene polr2a .
gist with a mutation in exon  <dig> of pdgfra did not show a significant upregulation of kit compared to normal tissue.
gist with a mutation in kit exon  <dig>   <dig> and in the pdgfra exon  <dig> showed significant pdgfrb downregulation.
discussion
most gastrointestinal stromal tumors exhibit mutations in exon  <dig> of the kit gene.
this exon encodes for the juxtamembrane domain of the receptor which possesses an autoinhibitory function.
in wt-gist, activating mutations are found neither in the kit nor in the pdgfra gene.
clinical treatment with the tyrosine kinase inhibitor imatinib targets the atp binding site in the kinase domain of the kit receptor.
in wt-gist, response to treatment is often poor  <cit>  although most of them overexpress the kit protein.
a screen for activating mutations in the juxtamembrane domain of the alternative receptor tyrosine kinases csf1r, flt <dig>  pdgfrb in  <dig> wt-gist samples yielded only wild type sequences  in agreement with previously published results  <cit> .
therefore we investigated in this study whether the expression of alternative receptor tyrosine kinases may contribute to the pathogenesis of wt-gist and thus may help to identify wt-gist subgroups with different response to imatinib and elucidate novel therapeutic targets.
the pathogenesis of several malignant tumors is associated with overexpression of csf1r and pdgfrb  <cit> .
imatinib inhibits cell invasion in malignant peripheral nerve sheath tumors by blocking pdgfrb  <cit>  and it has been found to have antitumor activity in patients with chordoma  <cit> .
in our study, the gene-expression level of csf1r, flt <dig> and pdgfrb was determined in a cohort of  <dig> gist samples.
furthermore, we assessed the expression of axl and met, two receptor tyrosine kinases which were found to be alternatively activated in therapy resistant gist  <cit> .
gene expression analysis by qpcr requires suitable reference genes.
the expression of reference genes like gapdh or beta-actin is regulated differentially depending on the tissue type.
therefore they are not suited as univocal reference genes  <cit> .
the determination of reference genes with stable expression in the experimental system used is essential to ensure accurate normalization and interpretation of results.
whereas fresh tissue is frequently not available for genetic analysis, ffpe material is the standard.
after formalin-fixation and paraffin-embedding of tissue, the isolated rna is often heavily fragmented.
in our study we used fresh as well as fixed material from the gastrointestinal tract to determine suitable reference genes and to analyze whether their expression levels are comparable.
we then validated the reference genes by qpcr in our gist cohort.
as reported also by others  <cit>  the variability of expression levels for the reference genes was very diverse comparing fresh frozen and formalin-fixed tissue.
therefore, we decided to validate separate reference genes for each tissue type.
the genes trfc, polr2a, ppia and rplpo were validated as appropriate reference genes for ffpe tissue.
for fresh frozen tissue, ppia and rplpo were also found to be suitable.
additionally, ubc is suited as a reference gene in fresh frozen tissue.
the variability of the reference genes for fresh frozen tissue was lower than for fixed tissue.
to overcome the problem of high variability in ffpe tissue, we pooled two independent cdna syntheses from one sample as recommended in the miqe guidelines  <cit> .
the miqe guidelines give considerations for a consistent application of the qpcr technology including experimental details, data analysis and reporting principles.
unequal efficiency of cdna synthesis might be a reason for deviation.
additionally, we tried to select consistent patient material for our cohort by choosing paraffin blocks having the same age.
it was shown by bibikova  <cit> , that cq-values in qpcr experiments depend on the age and condition of the tissue blocks.
because the expression value depends also on amplicon length  <cit> , only primer sets generating amplicons of about  <dig> bp were chosen for qpcr.
this corresponds to the fragment length of degraded rna between  <dig> and  <dig> bp.
the application of the miqe guidelines results in a minimum variability for reference genes.
in summary, considering the miqe guidelines ffpe material can be used reliably for expression analysis in gist, but the use of separate reference genes for ffpe tissue is indispensable.
our qpcr analysis included gist with wild type sequences in the hot-spot regions of kit and pdgfra , kit exon  <dig> and exon  <dig> mutated samples, pdgfra exon  <dig> mutated samples, nf-1-associated wt-gist and normal tissue controls.
the wt-gist and the samples with the two different kit mutations showed a significantly increased expression of kit in contrast to the normal tissue.
the results correspond to the immunohistochemical stainings of our samples and were in agreement with published data  <cit> .
it was shown immunohistochemically that pdgfra mutated gist have only a slightly increased protein expression of kit  <cit> .
our data revealed the same trend on rna level compared to the normal tissue.
additionally, we could show in all groups of mutated gist a significantly reduced expression of pdgfrb compared to the normal tissue.
this could be due to the concomitant increase of kit expression.
when comparing nf-1-associated gist without kit mutation with wt-gist without nf- <dig> association the latter showed a lower expression of kit.
thus, the results of the two wt-gist groups give a heterogenous profile, which suggests that different genomic events may be responsible for the development of these tumors.
comparing our sample groups with each other, no significant difference in the gene expression levels of flt <dig>  csf1r and axl were determined.
an exception was seen in the sample group with kit exon  <dig> mutation.
here, a significantly reduced expression of csf1r, flt <dig> and pdgfrb compared to the normal tissue was detected.
the results lead us to the assumption that kit exon  <dig> mutated gist play a special role compared to gist carrying other mutations.
interestingly, kit exon  <dig> mutated gist need a double daily dose of the tyrosine kinase inhibitor imatinib to be effectively treated  <cit> .
furthermore, they develop preferentially in the small intestine but only rarely in the stomach where the majority of gist are detected  <cit> .
all qpcr data were calculated four times with trfc, polr2a, ppia and rplpo as reference genes using the rest software.
concordant results with the four reference genes are based on the extensive and complex preselection of our cohort and the preparation of cdna synthesis according to the miqe guidelines.
CONCLUSIONS
in summary, we conclude that none of the alternative receptor tyrosine kinases analyzed here are associated with the pathogenesis of wild type or mutated gist.
it remains to be clarified whether an increased expression of receptor tyrosine kinase ligands is responsible for tumorigenesis of wt-gist as it is described for dermatofibrosarcoma protuberans  and tenosynovial giant cell tumor   <cit> .
further studies are needed to elucidate the role of ligand-driven pathogenesis in wt-gist.
authors' contributions
mcb, jf and hk carried out the molecular genetic studies and jf wrote and drafted the manuscript.
smb and hus developed the design of the study.
smb coordinated the study together with ew and hus.
bb and mz performed the statistical analysis of the expression data.
ew and hus participated in the diagnosis and selection of tumor material.
rb revised the manuscript for important intellectual content and approved the final manuscript.
all authors read and approved the manuscript.
acknowledgements and funding
this study was supported by a grant of the bonfor programme of the bonn medical centre to hus.

 gastrointestinal stromal tumors are the most common mesenchymal tumors of the gastrointestinal tract and are characterized by the expression of the kit receptor and to a lesser extent of pdgfra, representing two closely related receptor tyrosine kinases. the majority of gist shows oncogenic mutations either in kit or pdgfra. mainly, mutations in exon or of the kit gene or in exon of pdgfra lead to ligand independent, constitutive activation of the kinase function. about 60% of all gist carry an exon mutation of kit which encodes the juxtamembrane domain of the receptor possessing an autoinhibitory function. less common mutations in pdgfra are detected in gist that often display gastric location and epithelioid morphology. in a minority of cases no mutations in the known kit or pdgfra hot spots are detected although these tumors express the kit protein. this subgroup is called wild type gist and comprises tumors in pediatric patients, in patients affected by the carney triad, neurofibromatosis type associated gist and a subset of sporadic adult gist. the pathogenetic mechanisms underlying wt-gist are poorly understood and there is limited benefit of imatinib therapy in these patients. therefore the identification of additional genetic factors contributing to the pathogenesis of gist may help to find new concepts of individualized therapy. recently, the braf mutation p.v600e was found in 4-13% of wt-gist. for another subgroup of wt-gist including pediatric tumors, a strong igf1r expression combined partly with gene amplification was described. two other alternative rtk probably involved in the pathogenesis of gist are axl and met. both kinases have been shown to be upregulated in gist resistant to treatment. axl is a member of the ufo/axl subfamily and activates the same signaling pathway as kit. the tyrosine kinase domain of met is mutated in sporadic papillary renal carcinomas. some mutations in the met gene are located in codons homologous to those in kit and it is suggested that these missense mutations lead to constitutive activation of the met protein. to develop additional therapy approaches it would be of interest to know whether this rtk also plays a role in wt-gist. besides kit and pdgfra, csf1r, flt and pdgfrb belong to the same family of type iii rtk. these five tyrosine kinases show a homologous structure and a comparable function in activation, proliferation and suppressing apoptosis. aberrant expression and mutations in either csf1r, flt and pdgfrb or their ligands have been described in several malignant diseases. overexpression of csf1r is found in epithelial tumors such as breast and ovarian cancer. the translocation t of its ligand csf leads to the development of the tenosynovial giant cell tumor. aberrantly expressed flt is observed at high levels in a spectrum of hematologic malignancies. additionally, in aml an internal tandem duplication in the transmembrane domain of flt was identified which leads to constitutive activation of its kinase domain. it seems that this mutation is not present or very rare in gist. pdgfrb is overexpressed in malignant peripheral nerve sheath tumors and chordomas. to analyse expression profiles of certain tumors for research and diagnostic purposes, qpcr is frequently applied because of its reproducibility and high sensitivity. this method is based on the normalization of the target gene expression on stably expressed internal reference genes. a major challenge is the application of suitable reference genes which have to be tested and verified under defined experimental conditions. ideal reference genes have to be non-regulated, stable and not affected by biological or experimental conditions. the target gene is amplified together with the reference gene in order to minimize experimental variability concerning reverse transcription enzymatic efficiencies, pcr efficiency, amount of starting material and differences between human tissues. the reference gene and the target gene should have very robust and stable expression profiles to ensure accurate normalization and interpretation of results. the most stable expressed gene from a set of genes can be identified by genorm, a software program which additionally provides the number of genes required to calculate a robust normalization factor based on the geometric mean of these genes. typical reference genes regulate basic and ubiquitous cellular functions and are responsible for the cellular maintenance, e.g. gapdh or ß-actin. however, these commonly used reference genes vary considerably in different tissue types or under different experimental conditions. there is no standard reference gene for all kinds of tissue types. to our best knowledge only few studies investigated the alteration of stability of reference genes in different mesenchymal tumour entities. in the present study we attempted to identify suitable reference genes in gastrointestinal stromal tumors by using a set of sixteen reference genes which are currently applied in qpcr procedures. furthermore we evaluated the gene expression of kit as well as the alternative rtk flt csf1-r, pdgfrb, axl and met in mutated and non-mutated gastrointestinal stromal tumors by qpcr using the identified reference genes. the study was focussed on wt-gist compared to samples with mutations in kit exon and and pdgfra exon in order to evaluate whether overexpression of these alternative rtk might contribute to the pathogenesis of wt-gist. methods samples a total of samples were included into this study. all specimens were obtained in the years and under approved ethical protocols and with informed consent from each patient. all samples were fixed in neutral-buffered formalin prior to paraffin embedding. samples from normal tissue as control group and gist representing different mutational subgroups were evalutated. samples of wt-gist, samples of wt-gist associated with neurofibromatosis type samples with exon mutation in kit, samples with exon mutation in kit and samples with exon mutation in pdgfra). sequence analysis of kit and pdgfra was carried out as described earlier. additionally, in all samples the wild type status of the braf gene was ascertained. all gist samples were stained immunohistochemically for cd four normal tissues from the gastrointestinal region were available as fresh frozen and formalin-fixed, paraffin-embedded samples. rna extraction and cdna synthesis prior to rna extraction, paraffin-embedded tissues were cut into μm sections and mounted on glass-slides. six slides of each specimen were used for rna extraction. the sections were deparaffinized by extracting twice in xylene for min at room temperature. rehydration was done in 100% ethanol, 90% ethanol, 80% ethanol and 70% ethanol made with depc-water for min each. tumor areas previously marked on a h&e slide were scraped from the sections with a sterile scalpel. sections were transferred into a sterile ml tube. extraction and purification of rna was done using the rneasy ffpe kit according to the manufacturer's recommendations. after tumor localization by h&e staining fresh frozen material was rasped into μm thick pieces and rna purification was carried out using the rneasy kit. both rneasy kits contain a step of dnase treatment. finally rna from fresh and from fixed tissue was eluted in water. the quantification was done spectrophotometrically. ng of rna from each sample was reverse-transcribed using a random-hexamer primer and avian myoblastosis virus reverse transcriptase according to the manufacturer's protocol. cdna from those four fresh frozen and corresponding formalin-fixed, paraffin-embedded control samples was generated in only one reverse transcription and loaded on the reference low density arrays. further cdna samples were generated from the four ffpe control tissues, pooled within one patient and loaded onto a second independent reference low density array. identification of reference genes the reference genes used here were preselected because of their constitutive, non-regulated stable expression over a wide spectrum of tissues. but nevertheless the preselected reference genes are not suitable for every kind of tissue and therefore need to be analyzed prior to use in a certain study. the detection of suitable reference genes was carried out using taqman low density arrays. these arrays are prefabricated 384-well cards where gene-specific primer and probe sets are spotted in small reaction chambers during manufacturing. the cards have separate loading ports leading into wells each. in this study, arrays with triplicates of putative reference genes were used, so different samples could be analyzed. the list of assays is given in table for each sample, ng cdna was mixed with 2× taqman universal pcr master mix. μl of this mixture was loaded into each port and distributed into the reaction chambers by centrifugation. the card was sealed and the quantitative pcr was performed on an abi prism ht sequence detection system. after pre-incubation for min at 50°c and min at 95°c, the pcr reaction was performed. the fluorescent signal was measured in each cycle. analyses of gene expression by qpcr qpcr analysis was performed using the assays-on-demand products listed in table these gene-specific qpcr assays consist of a pair of unlabeled pcr primers and a fam labeled specific probe. according to the manufacturer of these assays, probe and primer sets that would amplify pseudogenes are excluded in the process of development. reactions were carried out in a reaction volume containing μl pcr master mix, μl forward and reverse primer mix, ng cdna ad μl a.dest. triplicate reactions were carried out for each transcript. control reactions were performed using a minus rt preparation and a sample with a.dest instead of rna. pcr conditions were the same as for the taqman low density arrays. statistical analyses the genorm applet for microsoft excel was used to determine the most stable genes among the sixteen candidate reference genes. raw cq values were converted into relative quantities for analysis with genorm, where the highest relative quantity for each gene is set to the program selects from a panel of candidate reference genes the two most stable genes or a combination of multiple stable genes for normalization. the gene expression stability value is based on the combined estimate of intra- and intergroup expression variations of the genes studied and takes the pcr efficiency into account. the limited m-value is. the most stable genes are stepwise selected from the investigated gene panel to estimate how many reference genes should be used. the normalization factors define the optimal number of reference genes required for a precised normalization design. the analysis of the expression data of tumor samples compared to control samples was performed with the rest software. the software normalizes the measured cq-values of the target genes with those of the reference genes and compares the expression data of tumor and control samples by considering the pcr efficiency and the mean crossing point deviation.  rna quality rna quality of the samples was inspected on a 1% agarose gel. the concentration and purity of the rna was characterized by the mean a260/ ratio and was on average for fresh frozen as well as for formalin fixed, paraffin-embedded samples and reflected pure and protein-free rna. expression variability of cdna synthesis four normal control ffpe tissue samples from the gastric and bowel wall were used to determine suitable reference genes. one tlda was done with cdna from only one reverse transcription for each sample. as it is assumed that cdna synthesis from ffpe tissue varies strongly in efficiency, pooled cdna synthesis samples from the same ffpe tissue were used for a second tlda. the pooled cdna preparation showed an essentially lower variance in contrast to the cdna samples from only one reverse transcription. the median of the m-values ranges from to, thus the variability of the technical replicates is low. within a single specimen the range is also very small. the median of the cdna samples from a single reverse transcription ranges from to. however, each cdna sample from one reverse non-pooled transcription showed a high variance of the values. therefore the variance of the samples can be lowered by cdna pooling. the median of the cdna samples from one reverse transcription from fresh frozen tissue ranged from to and can be compared to the median of the pooled cdna from the ffpe tissue. within a single specimen the range is higher than in the pooled cdna samples, but smaller than in the cdna from only one reverse transcription from the ffpe samples. identification of suitable reference genes for normalization all putative reference genes had a high expression stability and the'm'-value was clearly below the'm'-cutoff-value of. by stepwise exclusion of genes, the expression stability value'm' is calculated and the expression stability rises with the exclusion of further genes. samples generated from only one reverse transcription, showed polr2a, tfrc, rplpo and gapdh as the most stable genes. the average stability value was between and. the range of the values was -. the most stable genes for pooled cdna samples were pgk ppia, rplpo and ipo with a median of -. the range of the values was lower and varied between and. all in all the gene expression variability of the pooled cdna samples is much lower compared to the single reverse transcription cdna synthesis. the gene ubc is the most stable one in fresh frozen tissue with the smallest value range, but rplpo and ppia belong also to the most stable genes in fresh frozen tissue with a median of - and the range of the values was -. tbp was not expressed in the tissue. by combining the lowest variability values of fixed and fresh tissue, the genes polr2a, ppia, rplpo and tfrc were detected. hence these four genes were used for further analysis with the gist samples. expression profiles of kit, flt csf1-r, pdgfrb, axl and met six patient groups were used to analyze the gene expression of kit, flt csf1-r, pdgfrb, axl and met. here, the rest analysis with polr2a is shown exemplarily. the rest analysis of the target genes kit, csf1r, flt pdgfrb, axl and met with the reference genes rplpo, ppia and trfc showed the same significant results for the differential expression or at least the same trend in mutated and non-mutated gist. as suspected, we found a significant overexpression of kit in exon and exon mutated gist in comparison with normal tissue. also a significantly lowered expression of pdgfrb in both groups compared to normal tissue was shown. the same effect was observed in pdgfra exon mutated gist compared to normal tissue. expression analysis of kit, csf1r, flt pdgfrb, axl and met was carried out in six gist cohorts. significant variances are in bold. the expression data were calculated using the rest software with polr2a as reference gene. nt: normal tissue; p <; +/- indicates up- or downregulation on closer inspection of normal tissue and wt-gist a clear overexpression of kit in the tissue of wild type tumors was shown. the same effect, but even stronger, could be detected by comparing nf1-associated wt-gist with normal tissue. concerning kit expression, the tumors without mutation behave like kit mutated gist. this was in accordance with immunohistochemical staining. the other target genes, flt csf1r, pdgfrb, axl and met showed lower expression levels in the tumor tissue compared to normal tissue except axl in nf associated gist. met was even significantly reduced in wt-gist, kit exon and pdgfra exon mutated gist. no changes in expression levels were detected comparing normal tissue and wt-gist with kit exon mutated gist, indicating that both groups possess similar expression levels for the analyzed receptors. none of the groups showed significant expression alterations for csf1r and flt an exception is the group of kit exon mutated gist. a lower expression of csf1r, flt pdgfrb and met compared to normal tissue was identified in contrast to kit, which showed a significant overexpression. in summary, overexpression of kit was detected in each gist subgroup compared to the corresponding normal tissue by using the preassigned reference gene polr2a. gist with a mutation in exon of pdgfra did not show a significant upregulation of kit compared to normal tissue. gist with a mutation in kit exon and in the pdgfra exon showed significant pdgfrb downregulation. discussion most gastrointestinal stromal tumors exhibit mutations in exon of the kit gene. this exon encodes for the juxtamembrane domain of the receptor which possesses an autoinhibitory function. in wt-gist, activating mutations are found neither in the kit nor in the pdgfra gene. clinical treatment with the tyrosine kinase inhibitor imatinib targets the atp binding site in the kinase domain of the kit receptor. in wt-gist, response to treatment is often poor although most of them overexpress the kit protein. a screen for activating mutations in the juxtamembrane domain of the alternative receptor tyrosine kinases csf1r, flt pdgfrb in wt-gist samples yielded only wild type sequences in agreement with previously published results. therefore we investigated in this study whether the expression of alternative receptor tyrosine kinases may contribute to the pathogenesis of wt-gist and thus may help to identify wt-gist subgroups with different response to imatinib and elucidate novel therapeutic targets. the pathogenesis of several malignant tumors is associated with overexpression of csf1r and pdgfrb. imatinib inhibits cell invasion in malignant peripheral nerve sheath tumors by blocking pdgfrb and it has been found to have antitumor activity in patients with chordoma. in our study, the gene-expression level of csf1r, flt and pdgfrb was determined in a cohort of gist samples. furthermore, we assessed the expression of axl and met, two receptor tyrosine kinases which were found to be alternatively activated in therapy resistant gist. gene expression analysis by qpcr requires suitable reference genes. the expression of reference genes like gapdh or beta-actin is regulated differentially depending on the tissue type. therefore they are not suited as univocal reference genes. the determination of reference genes with stable expression in the experimental system used is essential to ensure accurate normalization and interpretation of results. whereas fresh tissue is frequently not available for genetic analysis, ffpe material is the standard. after formalin-fixation and paraffin-embedding of tissue, the isolated rna is often heavily fragmented. in our study we used fresh as well as fixed material from the gastrointestinal tract to determine suitable reference genes and to analyze whether their expression levels are comparable. we then validated the reference genes by qpcr in our gist cohort. as reported also by others the variability of expression levels for the reference genes was very diverse comparing fresh frozen and formalin-fixed tissue. therefore, we decided to validate separate reference genes for each tissue type. the genes trfc, polr2a, ppia and rplpo were validated as appropriate reference genes for ffpe tissue. for fresh frozen tissue, ppia and rplpo were also found to be suitable. additionally, ubc is suited as a reference gene in fresh frozen tissue. the variability of the reference genes for fresh frozen tissue was lower than for fixed tissue. to overcome the problem of high variability in ffpe tissue, we pooled two independent cdna syntheses from one sample as recommended in the miqe guidelines. the miqe guidelines give considerations for a consistent application of the qpcr technology including experimental details, data analysis and reporting principles. unequal efficiency of cdna synthesis might be a reason for deviation. additionally, we tried to select consistent patient material for our cohort by choosing paraffin blocks having the same age. it was shown by bibikova, that cq-values in qpcr experiments depend on the age and condition of the tissue blocks. because the expression value depends also on amplicon length, only primer sets generating amplicons of about bp were chosen for qpcr. this corresponds to the fragment length of degraded rna between and bp. the application of the miqe guidelines results in a minimum variability for reference genes. in summary, considering the miqe guidelines ffpe material can be used reliably for expression analysis in gist, but the use of separate reference genes for ffpe tissue is indispensable. our qpcr analysis included gist with wild type sequences in the hot-spot regions of kit and pdgfra, kit exon and exon mutated samples, pdgfra exon mutated samples, nf-1-associated wt-gist and normal tissue controls. the wt-gist and the samples with the two different kit mutations showed a significantly increased expression of kit in contrast to the normal tissue. the results correspond to the immunohistochemical stainings of our samples and were in agreement with published data. it was shown immunohistochemically that pdgfra mutated gist have only a slightly increased protein expression of kit. our data revealed the same trend on rna level compared to the normal tissue. additionally, we could show in all groups of mutated gist a significantly reduced expression of pdgfrb compared to the normal tissue. this could be due to the concomitant increase of kit expression. when comparing nf-1-associated gist without kit mutation with wt-gist without nf- association the latter showed a lower expression of kit. thus, the results of the two wt-gist groups give a heterogenous profile, which suggests that different genomic events may be responsible for the development of these tumors. comparing our sample groups with each other, no significant difference in the gene expression levels of flt csf1r and axl were determined. an exception was seen in the sample group with kit exon mutation. here, a significantly reduced expression of csf1r, flt and pdgfrb compared to the normal tissue was detected. the results lead us to the assumption that kit exon mutated gist play a special role compared to gist carrying other mutations. interestingly, kit exon mutated gist need a double daily dose of the tyrosine kinase inhibitor imatinib to be effectively treated. furthermore, they develop preferentially in the small intestine but only rarely in the stomach where the majority of gist are detected. all qpcr data were calculated four times with trfc, polr2a, ppia and rplpo as reference genes using the rest software. concordant results with the four reference genes are based on the extensive and complex preselection of our cohort and the preparation of cdna synthesis according to the miqe guidelines.  in summary, we conclude that none of the alternative receptor tyrosine kinases analyzed here are associated with the pathogenesis of wild type or mutated gist. it remains to be clarified whether an increased expression of receptor tyrosine kinase ligands is responsible for tumorigenesis of wt-gist as it is described for dermatofibrosarcoma protuberans and tenosynovial giant cell tumor. further studies are needed to elucidate the role of ligand-driven pathogenesis in wt-gist. authors' contributions mcb, jf and hk carried out the molecular genetic studies and jf wrote and drafted the manuscript. smb and hus developed the design of the study. smb coordinated the study together with ew and hus. bb and mz performed the statistical analysis of the expression data. ew and hus participated in the diagnosis and selection of tumor material. rb revised the manuscript for important intellectual content and approved the final manuscript. all authors read and approved the manuscript. acknowledgements and funding this study was supported by a grant of the bonfor programme of the bonn medical centre to hus.
by combining the lowest variability values of fixed and fresh tissue, the genes polr2a, ppia, rplpo and tfrc were chosen for further analysis of the gist samples. gene expression variability of the pooled cdna samples is much lower than the single reverse transcription cdna synthesis. approximately 10% of gist are so-called wild type gist without mutations in the hot spots. gastrointestinal stromal tumors represent the most common mesenchymal tumors of the gastrointestinal tract. overexpression of kit compared to the corresponding normal tissue was detected in each gist subgroup except in gist with pdgfra exon mutation. in the present study we evaluated appropriate reference genes for the expression analysis of formalin-fixed, paraffin-embedded and fresh frozen samples from gastrointestinal stromal tumors. comparing our sample groups, no significant differences in the gene expression levels of flt csf1r and axl were determined. about 85% carry an activating mutation in the kit or pdgfra gene. none of the alternative receptor tyrosine kinases analyzed is associated with the pathogenesis of wild-type or mutated gist. we evaluated the gene expression of kit as well as of the alternative receptor tyrosine kinase genes flt csf1-r, pdgfrb, axl and met by qpcr. an exception was the sample group with kit exon mutation. as the variability of expression levels for the reference genes is very high comparing fresh frozen and formalin-fixed tissue there is a strong need for validation in each tissue type. a significantly reduced expression of csf1r, flt and pdgfrb compared to the normal tissue was detected. wt-gist were compared to samples with mutations in kit exon and and pdgfra exon in order to evaluate whether overexpression of these alternative rtk might contribute to the pathogenesis of wt-gist. it remains to be clarified whether an autocrine or paracrine mechanism by overexpression of receptor tyrosine kinase ligands is responsible for the tumorigenesis of wt-gist. gist with mutations in kit exon and and in pdgfra exon showed a significant pdgfrb downregulation.
 gastrointestinal stromal tumors represent the most common mesenchymal tumors of the gastrointestinal tract. about 85% carry an activating mutation in the kit or pdgfra gene. approximately 10% of gist are so-called wild type gist without mutations in the hot spots. in the present study we evaluated appropriate reference genes for the expression analysis of formalin-fixed, paraffin-embedded and fresh frozen samples from gastrointestinal stromal tumors. we evaluated the gene expression of kit as well as of the alternative receptor tyrosine kinase genes flt csf1-r, pdgfrb, axl and met by qpcr. wt-gist were compared to samples with mutations in kit exon and and pdgfra exon in order to evaluate whether overexpression of these alternative rtk might contribute to the pathogenesis of wt-gist.  gene expression variability of the pooled cdna samples is much lower than the single reverse transcription cdna synthesis. by combining the lowest variability values of fixed and fresh tissue, the genes polr2a, ppia, rplpo and tfrc were chosen for further analysis of the gist samples. overexpression of kit compared to the corresponding normal tissue was detected in each gist subgroup except in gist with pdgfra exon mutation. comparing our sample groups, no significant differences in the gene expression levels of flt csf1r and axl were determined. an exception was the sample group with kit exon mutation. a significantly reduced expression of csf1r, flt and pdgfrb compared to the normal tissue was detected. gist with mutations in kit exon and and in pdgfra exon showed a significant pdgfrb downregulation.  as the variability of expression levels for the reference genes is very high comparing fresh frozen and formalin-fixed tissue there is a strong need for validation in each tissue type. none of the alternative receptor tyrosine kinases analyzed is associated with the pathogenesis of wild-type or mutated gist. it remains to be clarified whether an autocrine or paracrine mechanism by overexpression of receptor tyrosine kinase ligands is responsible for the tumorigenesis of wt-gist. 