Oncotarget

Research Papers:

FGFR2 is overexpressed in myxoid liposarcoma and inhibition of FGFR signaling impairs tumor growth in vitro

PDF |  HTML  |  Supplementary Files  |  How to cite

Oncotarget. 2015; 6:20215-20230. https://doi.org/10.18632/oncotarget.4046

Metrics: PDF 2867 views  |   HTML 4261 views  |   ?  

Helen Künstlinger _, Jana Fassunke, Hans-Ulrich Schildhaus, Benedikt Brors, Carina Heydt, Michaela Angelika Ihle, Gunhild Mechtersheimer, Eva Wardelmann, Reinhard Büttner and Sabine Merkelbach-Bruse

Abstract

Helen Künstlinger1, Jana Fassunke1, Hans-Ulrich Schildhaus2, Benedikt Brors3, Carina Heydt1, Michaela Angelika Ihle1, Gunhild Mechtersheimer4, Eva Wardelmann5, Reinhard Büttner1 and Sabine Merkelbach-Bruse1

1 Institute of Pathology, University Hospital Cologne, Cologne, Germany

2 Institute of Pathology, University Hospital Göttingen, Göttingen, Germany

3 Computational Oncology, German Cancer Research Center (DKFZ), Heidelberg, Germany

4 Institute of Pathology, University Hospital Heidelberg, Heidelberg, Germany

5 Gerhard-Domagk-Institute of Pathology, University Hospital Münster, Münster, Germany

Correspondence to:

Helen Künstlinger, email:

Keywords: myxoid liposarcoma, whole-genome microarray, gene expression, fibroblast growth factor, FGFR inhibition

Received: November 20, 2014 Accepted: April 22, 2015 Published: May 08, 2015

Abstract

Myxoid liposarcomas account for more than one third of liposarcomas and about 10% of all adult soft tissue sarcomas. The tumors are characterized by specific chromosomal translocations leading to the chimeric oncogenes FUS-DDIT3 or EWS1R-DDIT3. The encoded fusion proteins act as aberrant transcription factors. Therefore, we implemented comparative expression analyses using whole-genome microarrays in tumor and fat tissue samples. We aimed at identifying differentially expressed genes which may serve as diagnostic or prognostic biomarkers or as therapeutic targets. Microarray analyses revealed overexpression of FGFR2 and other members of the FGF/FGFR family. Overexpression of FGFR2 was validated by qPCR, immunohistochemistry and western blot analysis in primary tumor samples. Treatment of the myxoid liposarcoma cell lines MLS 402 and MLS 1765 with the FGFR inhibitors PD173074, TKI258 (dovitinib) and BGJ398 as well as specific siRNAs reduced cell proliferation, induced apoptosis and delayed cell migration. Combination of FGFR inhibitors with trabectedin further increased the effect. Our study demonstrates overexpression of FGFR2 and a functional role of FGFR signaling in myxoid liposarcoma. As FGFR inhibition showed effects on proliferation and cell migration and induced apoptosis in vitro, our data indicate the potential use of FGFR inhibitors as a targeted therapy for these tumors.


Introduction

Liposarcomas are rare malignant tumors that preferentially occur in deep soft tissue. However, representing 10-15% of all soft tissue sarcomas they are one of the most frequent sarcoma subtypes [1]. Based on histological and cytogenetic properties, liposarcomas can be distinguished into three different entities, i.e. well-differentiated/dedifferentiated (WD/DDLPS), myxoid/round cell (MRLPS) as well as pleomorphic liposarcomas (PLPS) [2, 3]. WD/DDLPS are genetically defined by giant marker or ring chromosomes with an amplification on chromosome 12 affecting among others the genes MDM2 and CDK4 [4]. Likewise, myxoid and round cell liposarcomas are considered as a common tumor entity which is characterized by a reciprocal translocation of the DDIT3 gene with either FUS (>90%) or EWS1R [5-8]. The translocations t(12;16)(q13;p11) and accordingly t(12;22)(q13;p12) are specific for this tumor entity and absent in other myxoid look-alikes such as myxofibrosarcoma [9]. The translocation leads to the fusion of the involved genes and the formation of a chimeric protein.

Besides the initial translocation only little is known about tumorigenic pathways deregulated by the chimeric protein [10]. As the DDIT3 fusion proteins are most likely to act as aberrant transcription factors, the transcriptional control of many genes may be altered. Such differentially expressed genes could be diagnostic or prognostic biomarkers as well as therapeutic targets. In order to identify multiple differentially expressed genes at the same time, cDNA microarrays are particularly suitable. They have already been performed in different sarcoma entities and revealed subtype specific expression signatures as well as basis for novel therapeutic approaches [11-15]. Nevertheless, the identification of candidate target structures based on gene expression profiling alone does not provide reliable evidence for the implementation of new therapeutic strategies. Rather, results need to be carefully validated and functional studies have to confirm the suitability of identified candidates.

Treatment of liposarcomas with the new chemotherapeutic compound trabectedin (ET-743) revealed promising results [16, 17], but molecularly targeted therapies are not yet available. For other soft tissue sarcomas some success with targeted therapies has been achieved in specific subtypes, but their implementation remains far behind treatment regimes in carcinomas [18, 19, 20]. An interesting family of receptor tyrosine kinases that can successfully be targeted and whose role in tumorigenesis has been revealed for different sarcoma entities are fibroblast growth factor receptors (FGFRs). FGFR1 amplification has been described in osteosarcomas and rhabdomyosarcomas were identified to carry activating mutations in FGFR4 [21, 22]. Furthermore, activation of FGFR signaling through amplification of the adaptor FRS2 has recently been described in high-grade liposarcomas [23].

The implementation of therapies with defined molecular targets requires the identification of new key molecules. Therefore, we implemented comparative whole-genome microarray analyses in primary myxoid liposarcomas and fat tissue samples. FGFR2, together with other members of the FGF/FGFR family, showed overexpression. FGFR2 expression was further analyzed in primary tumors and myxoid liposarcoma cell lines were treated with FGFR silencing siRNAs and tyrosine kinase inhibitors. Inhibitors were additionally combined with trabectedin. Our study investigates a potential role of FGFR signaling in myxoid liposarcomas and the use of FGFR inhibitors as a novel targeted treatment approach for these tumors.

Results

Microarray analyses

In order to identify new key molecules in the pathogenesis of myxoid liposarcoma, whole-genome microarray analyses were performed with seven tumor samples and an RNA pool of eight normal fat tissue samples with excellent RNA quality (RNA integrity number (RIN) value ≥ 7.0, Table 1 and Supplemental Figure S1 A).

Table 1: Cohort of tumor and fat control samples. Samples with RIN values printed in bold were used for microarray analysis.

Table-1

A principal component analysis (PCA) was carried out to analyze the relative similarity between the expression profiles. The PCA plot (Figure 1A) illustrates a clear difference between the fat control pool and the tumor samples, which show a high degree of similarity among each other. Due to this similarity the tumor samples were compared to the control as a whole. By means of this comparison a multitude of differentially expressed genes could be identified. We detected 7,946 genes with significant twofold or more up- or down-regulation (p-value < 0.01). Maximal measured fold changes were ±100. In the range of +80 to +100 and –80 to –100 were 26 genes, respectively. A heat map reflecting the differential expression of several genes in MLS tumor samples compared to normal fat tissue is given in Supplemental Figure S1B.

All MLS samples included in the study were proven to carry a DDIT3 translocation involving either FUS or EWS1R using FISH analysis (Supplemental Figure S1C). Additionally, in samples used for microarray analysis and MLS cell lines the exact type of DDIT3 translocation was determined by RT-PCR and sequencing (Supplemental Table S1). In two tumor samples our analysis revealed so far undescribed variants. Case MLS 12 carries both fusion transcripts of type I and type III. For MLS 1 a variant of type I lacking the last 24 codons of FUS exon 7 was identified. A correlation between the differential expression of genes and the particular type of fusion transcript was not detected.

Whole-genome microarray analysis and validation.

Figure 1: Whole-genome microarray analysis and validation. A. Three dimensional principal component analysis (PCA) of the gene expression profiles. On X, Y and Z axis arbitrary units for the three different principal components are indicated. Each profile is condensed to a single data point in the three dimensional PCA. The data points’ size and distance reflect the three dimensional position within the grid and thus indicate the similarity among the different profiles. B. qPCR validation of selected candidate genes in cryo-conserved (top) and FFPE (bottom) primary tumor samples. Red or green colored bars represent genes with confirmed up- or down-regulation (same expression changes as discovered by microarray analysis previously). Data are represented as analysis output of REST [48] (mean fold change of triplicate measurements). *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001.

Microarray validation and identification of suitable reference genes

We analyzed the expression stability of 16 candidate reference genes by means of qPCR using TaqMan® Array Human Endogenous Control Panels in cryo-conserved and formalin-fixed tumor tissue as well as in cryo-conserved fat tissue samples. Data were analyzed using GeNorm and NormFinder software tools, which both detected comparable results. Expression stabilities generated with GeNorm are shown in Supplemental Figure S2. Suitable reference genes should exhibit stability values below 1.5, which was seen in cryo-conserved tissue samples for all genes apart from PGK1. In formalin-fixed tissue, stability values for the candidate reference genes turned out to be much higher. As the qPCR validation was also performed in a cohort of formalin-fixed samples, these data were considered preferentially for selection of reference genes. IPO8 was among the two most stable reference genes in all three tissue types and was therefore selected as one reference gene. B2M was chosen as the second reference gene due to its stability in formalin-fixed as well as in cryo-conserved tissue samples. The evaluation of expression data with NormFinder revealed likewise stability values below 1.5 for IPO8 and B2M in all three tissue types.

Results of microarray analyses were validated for eleven genes in the whole tumor group by SYBR Green qPCR. The analyzed genes were selected from those genes with at least twenty fold significantly differential expression (Table 2). RNA pools of cryo-conserved or post-fixed fat tissue were used as controls. Data obtained by qPCR experiments were evaluated using REST (Relative Expression Software Tool) and graphically visualized using GraphPad Prism (Figure 1B). Validation in cryo-conserved tumor samples revealed the same differential expression of all tested genes as measured with microarray analysis. In formalin-fixed tumor tissue ten of eleven genes showed the same expression changes as discovered by microarray analyses previously. Thus, results of microarray analyses could be reproduced by qPCR.

Table 2: Candidate genes for qPCR validation. Differential expression of selected candidate genes as detected with whole-genome microarrays. Genes in the upper section were up-regulated whereas genes in the bottom section were down-regulated in MLS tumor samples compared to fat tissue as control.

Table-2

FGFR expression in myxoid liposarcoma

By evaluating the microarray analyses multiple genes expressed at significantly different levels were identified. One of these genes was FGFR2, which was highly overexpressed in myxoid liposarcoma. In addition to FGFR2 also other members of the FGF/FGFR family showed an overexpression in the microarray analyses (Figure 2A) further reinforcing the potential role of FGFR signaling in myxoid liposarcomas. In qPCR analyses a clear up-regulation of FGFR2 gene expression was detected in cryo-conserved as well as in formalin-fixed tumor samples (Figure 2B) confirming the overexpression identified by microarray analyses in the whole tumor cohort.

FGFR2 expression in primary myxoid liposarcomas was also shown on protein level by immunohistochemistry (Figure 2C) and western blot analysis (Figure 2D). Vascular endothelial cells as well as the tumor cells were positive for FGFR2. FGFR2 protein expression was also present in the myxoid liposarcoma cell lines MLS 402 and MLS 1765 (Figures 2C and 2D). In contrast, fat tissue samples showed only weak or no FGFR2 protein expression (Figure 2D).

FGFR expression.

Figure 2: FGFR expression. A. Differential expression of members of the FGF/FGFR family as detected with the microarray analysis. B. Differential expression of FGFR2 in myxoid liposarcomas detected with qPCR in the whole tumor cohort. Data are represented as analysis output of REST [48] (mean fold change of triplicate measurements). *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001. C. FGFR2 immunohistochemistry in MLS primary tumor tissue (left) and cell lines (right). D. FGFR2 protein expression in cryo-conserved tumor and fat tissue samples as well as in MLS cell lines detected by western blot analysis. Kato-III cells served as positive control for FGFR2 expression.

FGFR silencing with specific siRNAs

In order to analyze the functional role of FGFR signaling in the myxoid liposarcoma cell lines MLS 402 and MLS 1765, cells were transiently transfected with specific siRNAs. Silencing efficiency of the used siRNAs was shown by qPCR 48 h after transfection (Supplemental Figure S3). As both cell lines express all four FGFRs (Figure 3A), the receptors were knocked-down separately as well as in combination. Figure 3B shows apoptosis induction upon FGFR silencing, either of FGFR1, 2, 3 and 4 alone or by combined knock-down of several receptors. In MLS 402 single receptor silencing of FGFR2 induced apoptosis, whereas knock-down of the other receptors alone showed no effect (Figure 3B, left upper graph). Accordingly, the combined knock-down of more than one FGF receptor induced apoptosis only when FGFR2 siRNA was included (Figure 3B, left bottom graph). In MLS 1765 cells knock-down of FGFRs induced apoptosis as well, but to a lesser extent than in MLS 402. Silencing of FGFR2 alone did not cause apoptosis, but apoptosis induction was detected upon knock-down of either FGFR1 or FGFR3 (Figure 3B, right upper graph). In line with this all combinations of different siRNAs including those against FGFR1 and/or FGFR3 induced apoptosis in MLS 1765 cells (Figure 3B, right bottom graph).

In MLS 402 cells the combined silencing of the four FGF receptors caused attenuation of ERK mediated downstream signaling, shown by the reduced phosphorylation of ERK1/2 (Figure 3C). In MLS 1765 cells even single receptor silencing of either FGFR1 or FGFR3 lead to a reduction of ERK1/2 phosphorylation (Figure 3C). This result is in concordance with apoptosis induction in MLS 1765 cells after silencing of FGFR1 or FGFR3.

Effects of FGFR siRNA silencing.

Figure 3: Effects of FGFR siRNA silencing. A. FGFR1, FGFR2, FGFR3 and FGFR4 protein expression in MLS cell lines detected by western blot analysis. B. Apoptosis induction upon FGFR siRNA silencing in MLS cell lines, either of the different receptors alone (top) or combined knock-down of several FGFRs (bottom). As negative control cells were transfected with Stealth RNAi™ Negative Control Duplexes with the corresponding GC content (Contr. F1, F2, F3, F4). AllStars Hs Cell Death siRNA served as transfection control and as positive control for apoptosis induction (Pos.Contr.). Data of quintuplicates are represented as mean +/- SEM. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001. F1:FGFR1; F2:FGFR2; F3: FGFR3; F4:FGFR4. C. Effects of FGFR silencing on downstream signaling in MLS cell lines 48 h after transfection with specific siRNAs. As negative control cells were transfected with Stealth RNAi™ Negative Control Duplexes with the corresponding GC content.

FGFR inhibition with small molecules

As silencing of FGF receptors with specific siRNAs induced apoptosis in the myxoid liposarcoma cell lines MLS 402 and MLS 1765, also the functional effects of small molecules directed against FGFRs were examined. Cells were treated with the in vitro compound PD173074 (purchased from Sigma Aldrich) as well as with the two clinically applicable FGFR inhibitors dovitinib (TKI258) and BGJ398 (both from Novartis). As shown in Figure 4A, viability of both cell lines was reduced through FGFR inhibition by all three compounds in a concentration dependent manner. IC50 values were in the micromolar range and inhibitory effects could be seen under serum-reduced as well as under full-serum conditions. The detected effects upon FGFR inhibition by TKI258 and BGJ398 on the viability of the MLS cell lines were compared to those in a known FGFR inhibitor sensitive cell line. For that purpose the FGFR2 amplified gastric cancer cell line Kato-III was used. This cell line was markedly more sensitive to the highly FGFR specific compound BGJ398, whereas the sensitivity towards the more broad-spectrum tyrosine kinase inhibitor TKI258 was similar in MLS and Kato-III cells (Supplemental Figure S4A).

Effects of FGFR inhibitors on the viability of myxoid liposarcoma cell lines were further differentiated using the ApoTox-Glo™ Triplex Assay (Promega). Inhibitor concentrations reducing the viability of the cells (1 – 10 µM) were proven to be not cell toxic, but specifically inducing apoptosis (Figure 4B). BGJ398 and TKI258 already induced apoptosis after 48 hours of treatment in both cell lines, whereas apoptosis induction upon treatment with PD173074 in MLS 402 was detected only after 72 hours (Supplemental Figure S4B).

Effects of FGFR inhibitors on cell viability and apoptosis induction.

Figure 4: Effects of FGFR inhibitors on cell viability and apoptosis induction. A. Effects of FGFR inhibition with PD173074, TKI258 and BGJ398 on the viability of MLS cell lines. Cells were treated with different inhibitor concentrations under full serum (10% FBS) or serum reduced (2% FBS) conditions and analyzed using MTT assay after 48 h. Data of quintuplicates are represented as mean +/- SEM. B. Apoptosis induction in MLS cell lines 48 h after treatment with FGFR inhibitors. As negative control cells were treated with 0.1% DMSO. Camptothecin treated cells served as positive control for apoptosis induction (Pos.Contr.). FGFRi: FGFR inhibitors. Data of quintuplicates are represented as mean +/- SEM. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001.

FGFR inhibition in myxoid liposarcoma cell lines reduced kinase activity of the receptors, as detected by reduced receptor phosphorylation (Supplemental Figure S5A). Furthermore, the downstream signaling was affected, as shown exemplarily for ERK1/2 (Figure 5A). All three FGFR inhibitors reduced the phosphorylation of ERK1/2, being one of the main downstream targets of FGFRs.

Additionally, the effects of FGFR inhibition on the migratory activity of MLS 402 and MLS 1765 cells were examined by scratch assay. The FGFR inhibitors PD173074, BGJ398 and TKI258 were applied in a final concentration of 0.1 µM, a concentration which does not affect cell viability in the regarded timeframe of 24 hours. Supplemental Figure S5B exemplarily shows the effects of BGJ398 on the migration of MLS 1765 cells. Scratch closure was examined in FGFR inhibitor treated and DMSO treated control cells and the migratory activity was calculated as relative scratch closure compared to control. All three FGFR inhibitors reduced the migratory activity of both myxoid liposarcoma cell lines (Figure 5B). However, PD173074 and BGJ398 showed much stronger effects than TKI258, which had only little impact on the migration of MLS 402 cells. In MLS 1765 the effect of TKI258 on cell migration was indeed more obvious but still not significant.

Effects of FGFR inhibitors on downstream signaling and cell migration.

Figure 5: Effects of FGFR inhibitors on downstream signaling and cell migration. A. Effects of FGFR inhibitors [1 µM] on downstream signaling through phosphorylation of ERK1/2. Protein lysates were collected 5 h after treatment with the respective inhibitor. Lysates of cells treated with 0.1% DMSO served as controls. B. Effects of FGFR inhibitors [0.1 µM each] on the migratory activity of MLS cell lines. Control cells were treated with 0.1% DMSO. Data of triplicates are represented as mean +/- SEM. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001.

Combination of FGFR inhibitors with trabectedin

Additionally to single agent treatment, FGFR inhibitors were applied in combination with trabectedin (ET-743). The corresponding combination scheme is given in Figure 6A. As shown in Figure 6B for MLS 1765, the supplementary administration of FGFR inhibitors enhanced the effects of trabectedin in comparison to single agent treatment. Cell viability was reduced and apoptosis was induced by all three inhibitors in combination with trabectedin in a concentration below its single agent efficacy. The same holds true for the migratory activity of myxoid liposarcoma cell lines. In MLS 402 as well as in MLS 1765 cells migration was delayed by the additional administration of FGFR inhibitors when compared to cells treated with trabectedin alone (Figure 6C). As seen in the single agent migration assay, PD173074 and BGJ398 displayed a greater impact on the migration of myxoid liposarcoma cells than TKI258.

Combined treatment with FGFR inhibitors and trabectedin.

Figure 6: Combined treatment with FGFR inhibitors and trabectedin. A. Combination scheme of FGFR inhibitors and trabectedin. B. Effects of combined treatment with FGFR inhibitors together with trabectedin on MLS 1765 cells in comparison to trabectedin alone. Cells were treated with different combinations of compound concentrations (as shown in A) and effects on cell viability as well as induction of apoptosis were determined after 48 h. As negative control cells were treated with 0.2% DMSO. Digitonin and camptothecin treated cells served as positive controls for reduction of viability and induction of apoptosis, respectively. Data of quintuplicates are represented as mean +/- SEM. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001. C. Migratory activity of MLS 402 and MLS 1765 cells under treatment with trabectedin [0.1 nM] alone or combined with FGFR inhibitors [0.1 µM each]. Control cells were treated with 0.2% DMSO. Data of triplicates are represented as mean +/- SEM. *: p ≤ 0.05; **: p ≤ 0.01; ***: p ≤ 0.001. ET-743 / ET: trabectedin, FGFRi: FGFR inhibitors, PD: PD173074, BGJ: BGJ398, TKI: dovitinib (TKI258).

Discussion

The treatment of soft tissue sarcomas in general and of myxoid liposarcomas in particular is so far restricted to unspecific chemotherapy and preferably total surgical resection. The implementation of new specific therapies requires the identification of suitable targetable molecular structures. An approach to reveal such target structures is comparative gene expression profiling as applied in this study by whole-genome microarray analyses on seven myxoid liposarcoma samples and an RNA pool of eight fat tissue samples. That way differential expression of multiple genes was identified. A principal component analysis demonstrated a joint clustering of the tumor samples apart from the control sample. We detected nearly 8,000 genes with a significant fold change of at least ±2. Compared to other gene expression profiling studies in sarcomas the number of differentially expressed genes identified in this study is rather high, but previous studies mainly compared different sarcoma entities to each other and did not include whole-genome expression analyses [11, 12, 14]. Furthermore, the high number of differentially expressed genes demonstrates the strong homogeneity in our tumor cohort. This feature of myxoid liposarcomas became already apparent in previous gene expression and methylation profiling studies that observed tight clusters for myxoid liposarcomas clearly distinguishable from other entities or corresponding normal tissue [12-14].

In order to confirm the reliability of microarray data by qPCR validation, it is important to determine suitable references because a universal reference for all kinds of tissues and experimental conditions does not exist [24-27]. Our data revealed that widely used reference genes, such as GAPDH, are not self-evidently appropriate for formalin fixed tissue with intense degradation of RNA. We could identify B2M and IPO8 as suitable reference genes for myxoid liposarcoma and fat control samples in cryo-conserved as well as in formalin-fixed tissue. The herein presented whole-genome expression analyses could be reproduced very well by means of qPCR even by using FFPE tissue and served therefore as method to identify differentially expressed candidate genes.

Receptor tyrosine kinases belong to the most promising target structures as their constitutive activation has been shown in various tumor entities and they can successfully be targeted [19, 28, 29]. A constitutively active receptor tyrosine kinase may be due to DNA mutations or gene amplification but as well overexpression can lead to permanent receptor signaling [30, 31]. The present study has revealed the receptor tyrosine kinase fibroblast growth factor receptor 2 (FGFR2) as a gene of interest being highly overexpressed in myxoid liposarcoma.

FGF receptors were already reported to play a role in various cancers and their inhibition is effectively applied in targeted therapies with specific small molecule inhibitors [32-38]. Even in sarcomas a role of FGF receptors was described. Girnita et al. showed the importance of the bFGF pathway for the maintenance of a malignant phenotype of Ewing’s sarcoma cells [39]. The functional role and the effective inhibition of FGF receptors in synovial sarcoma were shown by Ishibe et al. [40]. A recent study in canine and human sarcomas detected co-expression of various FGFRs and demonstrated growth inhibitory effects by FGFR inhibition in vitro [41]. An expression study in different sarcomas by Baird et al. revealed an association of FGFR2, FGFR4 and FGF18 gene expression with liposarcomas [11]. This association became also apparent in the present study where we could detect overexpression not only of FGFR2 but also of FGFR4 as well as of the ligands FGF5, FGF11 and FGF18. The parallel overexpression of several FGF/FGFR family members reinforces a potential role of FGFR signaling in the pathogenesis of myxoid liposarcoma. Furthermore, FGFR signaling is described as being essential in preadipocyte differentiation providing another possible connection to the pathogenesis of liposarcomas that are characterized by the presence of premature adipocytes [42].

In vitro experiments in myxoid liposarcoma cell lines confirmed a functional role of FGFR signaling. The administration of specific siRNAs directed against FGFRs induced apoptosis in MLS 402 and MLS 1765 cells whereby the two cell lines reacted differently to selective FGFR silencing. While cell survival in MLS 402 cells seems to be dependent on FGFR2, in MLS 1765 rather FGFR1 and FGFR3 are important. Thus, in myxoid liposarcoma cells there is obviously not a single FGF receptor acting as a strong tumor driver. Accordingly, small molecules inhibiting the whole group of FGF receptors reduced cell viability, induced apoptosis and delayed migration in both myxoid liposarcoma cell lines.

In the examined myxoid liposarcoma cell lines the impact of FGFR signaling on survival and migration seemed to outweigh that on proliferation. FGFR inhibitor concentrations needed to achieve the observed effects differed in at least one decimal power. Furthermore, the FGFR specific inhibitors PD173074 and BGJ398 had a much greater impact on cell migration than the rather broad-spectrum tyrosine kinase inhibitor TKI258 also sustaining the hypothesis that in myxoid liposarcoma the role of FGFR signaling is especially prominent for cell migration. The herein used FGFR inhibitors TKI258 and BGJ398 are already applied in clinical trials of other tumor entities, such as FGFR3 mutated bladder cancer or squamous cell lung cancer with FGFR1 amplification (e.g. NCT01004224 or NCT01831726). We now demonstrated FGFR overexpression in myxoid liposarcoma and the effective application of FGFR inhibition in vitro indicating a novel targeted therapy approach. This approach might involve the combination of FGFR inhibitors with chemotherapeutic agents like trabectedin. In this study we showed that the additional administration of FGFR inhibitors can improve the impact on myxoid liposarcoma cells in comparison to trabectedin alone.

In order to examine whether FGFR inhibitors alone or in combination with chemotherapy are efficient as a potent therapy in myxoid liposarcoma, further in vivo studies are needed. In an in vivo setting also the role of FGFR signaling in angiogenesis will surely become relevant. Several FGFs are described to promote neoangiogenesis, among them FGF2 but also FGF5 and FGF18, for which an overexpression in myxoid liposarcomas was detected in this study [43-45]. Given the fact that myxoid liposarcomas are characterized by a typical branching capillary network, an inhibitor blocking also angiogenic pathways might be particularly effective in this entity.

In conclusion, our study revealed overexpression of FGFR2 as well as a functional role of FGFR signaling in myxoid liposarcoma and provides basis for the use of FGFR inhibitors as a novel targeted treatment approach for these tumors. Further in vivo studies are warranted to affirm the efficacy of FGFR inhibition either alone or in combination with chemotherapy.

materials and Methods

Human tissue samples

A cohort of 16 human primary myxoid liposarcomas and 15 fat tissue samples was assorted (Table 1). Tumor samples were provided by the Biobank of the Center for Integrated Oncology (CIO) Cologne Bonn and by the Institute of Pathology in Heidelberg. Fat tissue samples were obtained during intra-operative diagnostics in the scope of the CIO Cologne Bonn. Cryo-conserved material of eight fat tissue samples was additionally post-fixed in formalin and paraffin-embedded (FFPE) to obtain suitable controls for FFPE tumor samples. The study was approved by the ethics committees of the medical schools of the universities Bonn and Cologne. All cases were histologically characterized by an experienced pathologist (HUS, GM, EW or RB) and DDIT3-translocation was proven with fluorescence in-situ hybridization (FISH) and reverse transcription (RT)-PCR. FISH probes for DDIT3, FUS and EWS1R were from Abbott Molecular (Abbott Park IL, U.S.A.). Primer sequences for RT-PCR reactions are given in Supplemental Table S1.

RNA extraction, quality control and cDNA synthesis

Total RNA extraction was performed with miRNeasy mini or miRNeasy FFPE Kit (both Qiagen, Hilden, Germany). RNA was quantified using Nanodrop 1000 spectrophotometer (Thermo Scientific, Waltham MA, U.S.A.). Assessment of RNA degradation was performed via Agilent 2100 Bioanalyzer and RNA 6000 Nano Assay (Agilent Technologies, Santa Clara CA, U.S.A.). cDNA was synthesized using random hexamer primers and Omniscript Reverse Transcription Kit (Qiagen). Two independent cDNA syntheses were pooled for each sample. In order to detect FGFR knock-down after transfection of specific siRNAs the High-Capacity RNA-to-cDNA™ Kit (Life Technologies, Carlsbad CA, U.S.A.) was used for cDNA synthesis.

Microarray analyses

RNA from seven tumor samples and an RNA pool of eight fat tissue samples with a minimal RIN value of 7.0 were selected for whole-genome microarray analyses with Agilent Whole Human Genome Oligo Microarrays 4x44K. Microarray analyses were carried out by the genomic service of Miltenyi Biotec (Bergisch Gladbach, Germany). Raw data of microarray hybridizations have been deposited in NCBI’s Gene Expression Omnibus [46] and are accessible through GEO Series accession number GSE62747 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62747).

Quantitative real-time PCR (qPCR)

qPCR reactions were performed using Applied Biosystems 7900HT Fast Real-Time PCR System (Life Technologies) or LightCycler® 480 II system (Roche Diagnostics, Rotkreuz, Switzerland). Expression stability of 16 potential reference genes was analyzed with TaqMan® Array Human Endogenous Control Panels (Life Technologies). Obtained raw Cq (cycle of quantification) values were converted into relative quantities and analyzed using GeNorm [47] and NormFinder (MDL, Aarhus, Denmark). Expression of candidate genes was analyzed using Power SYBR® Green PCR Kit (Life Technologies). Corresponding primer pairs are given in Supplemental Table S2. Gene expression data were analyzed using the Relative Expression Software Tool (REST) [48] and graphically visualized using GraphPad Prism (GraphPad Software, La Jolla CA, U.S.A.). In order to detect FGFR knock-down in MLS cell lines after siRNA transfection TaqMan® Gene Expression Assays for FGFR1, 2, 3, and 4 as well as the reference genes IPO8 and B2M were used along with the TaqMan® Universal Master Mix II (all from Life Technologies).

FGFR2 expression in tissue samples

Gene expression of FGFR2 in primary tumor samples was detected by qPCR with Power SYBR® Green PCR Kit (Life Technologies) and specific primers (for 5’-GTGAAACTTGGTACTTCATGGTGA, rev 5’-GAGATGGCATTCTTGTTGTTACTG). Expression data were normalized using IPO8 and B2M as reference genes and relative gene expression was calculated with REST. Immunohistochemistry for FGFR2 was performed using rabbit anti-FGFR2 antibody from Zytomed Systems (Berlin, Germany). FGFR2 protein expression in primary tumor and fat tissue samples was analyzed with western blot using rabbit anti-FGFR2 antibody from Sigma Aldrich (St. Louis MO, U.S.A.) in a dilution of 1:500. For detection of reference protein expression rabbit anti-HPRT antibody from abcam (Cambridge, GB) in a dilution of 1:1000 was used.

Cell culture and FGFR inhibition in vitro

Myxoid liposarcoma cell lines MLS 402 and MLS 1765 were provided by Prof. Pierre Åman, University of Gothenburg, Sweden and cultured in RPMI 1640 medium supplemented with 2 mM L-glutamine and 2% or 10% FBS under humidified conditions with 5% CO2. Cell lines’ entity was verified by proof of the specific FUS-DDIT3 translocation. Furthermore, uniqueness of the cell lines was confirmed by STR profiling performed at DSMZ, Braunschweig, Germany. Gastric cancer cell line Kato-III was purchased from Cell Lines Service (Eppelheim, Germany) and cultured in F12 medium supplemented with 2 mM L-glutamine and 10% FBS under humidified conditions with 5% CO2. FGFR2 amplification in Kato-III cells was proven with FISH analysis using Zytolight SPEC FGFR2/CEN 10 Dual Color Probe from ZytoVision (Bremerhaven, Germany).

FGFR inhibitors PD173074 (purchased from Sigma Aldrich), BGJ398 and TKI258 (both provided by Novartis, Basel, Switzerland) as well as trabectedin (provided by PharmaMar, Madrid, Spain) were solved and prediluted in DMSO. Final DMSO concentration during assays was 0.1% for single agent or 0.2% for combined treatment. For FGFR silencing cells were transiently transfected with 100 nM siRNA using Lipofectamine™ RNAiMax (Life Technologies) and Stealth RNAi™ siRNAs against FGFR1, 2, 3 and 4 (set of 3 siRNAs each) as well as Stealth RNAi™ Negative Control Duplexes (all from Life Technologies). AllStars Hs Cell Death siRNA (Qiagen) served as transfection control.

Detection of cell viability, toxicity and apoptosis

Cell viability was measured in quintuplicate using MTT assay. 10 µl MTT staining solution (5 mg/ml in PBS, sterile filtered) were added to each well of a 96well plate and reaction was stopped after 5 h with 100 µl MTT solvent (10% SDS in 0.01 M HCl). Formed crystals were lyzed overnight at 37°C and absorbance was detected at 550 nm deducting background at 690 nm. To further differentiate observed effects cells were analyzed using ApoTox-Glo™ Triplex Assay (Promega, Madison WI, U.S.A.).

Western blot

Total protein lysates were extracted from untreated cell lines as well as 48 h after siRNA transfection or 2 or 5 h after inhibitor treatment. For western blot analysis equal amounts of total protein were loaded on NuPAGE® Bis-Tris Gels (Life Technologies). Electrophoretically separated proteins were transferred to nitrocellulose membrane and unspecific binding was blocked with 5% milk/PBS-Tween. Membranes were incubated with primary antibodies overnight at 4°C (FGFR2 1:500; ACTB 1:5000, AC-15, both Sigma Aldrich; FGFR1 1:500; FGFR4 1:500, AM11076PU-N, both from Acris Antibodies, San Diego CA, U.S.A.; FGFR3 1:500, D2G7E; phospho-FGFR 1:500, 55H2; ERK1/2 1:1000; phospho-ERK1/2 1:1000, 20G11, all from CellSignaling Technology, Danvers MA, U.S.A.; HPRT 1:1000, abcam). After incubation with HRP-conjugated secondary antibody SuperSignal West Pico or Femto Chemiluminescent Substrate (Thermo Scientific) was used for detection.

Scratch assay

Cells were grown to 100% confluence in medium containing 10% FBS in 24 well culture plates. In each well a scratch was set crosswise with a 100 µl pipet tip, debris was washed away using PBS. Cells were incubated for 24 h with medium containing 2% FBS and the respective treatment. Pictures of each well with centered scratch were taken before and after treatment period and uncovered areas were determined using AxioVision 4.8 software (Carl Zeiss, Oberkochen, Germany). Scratch closure was calculated as the proportion of previously uncovered area closed by cell migration after the considered period. Migratory activity was defined as relative scratch closure compared to control.

Acknowledgments

This project was funded by the Wilhelm Sander Foundation for cancer research. We thank Prof. Pierre Åman for kindly providing the myxoid liposarcoma cell lines MLS 402 and MLS 1765. We thank the companies Novartis and PharmaMar for the complimentary supply with the FGFR inhibitors BGJ398 and TKI258 and trabectedin, respectively.

References

1. Wibmer C, Leithner A, Zielonke N, Sperl M and Windhager R. Increasing incidence rates of soft tissue sarcomas? A population-based epidemiologic study and literature review. Ann Oncol. 2010; 21:1106-1111.

2. Hameed M. Pathology and genetics of adipocytic tumors. Cytogenet Genome Res. 2007; 118:138-147.

3. Coindre JM, Pedeutour F and Aurias A. Well-differentiated and dedifferentiated liposarcomas. Virchows Arch. 2010; 456:167-179.

4. Muller CR, Paulsen EB, Noordhuis P, Pedeutour F, Saeter G and Myklebost O. Potential for treatment of liposarcomas with the MDM2 antagonist Nutlin-3A. Int J Cancer. 2007; 121:199-205.

5. Antonescu CR, Tschernyavsky SJ, Decuseara R, Leung DH, Woodruff JM, Brennan MF, Bridge JA, Neff JR, Goldblum JR and Ladanyi M. Prognostic impact of P53 status, TLS-CHOP fusion transcript structure, and histological grade in myxoid liposarcoma: a molecular and clinicopathologic study of 82 cases. Clin Cancer Res. 2001; 7:3977-3987.

6. Crozat A, Aman P, Mandahl N and Ron D. Fusion of CHOP to a novel RNA-binding protein in human myxoid liposarcoma. Nature. 1993; 363:640-644.

7. Panagopoulos I, Hoglund M, Mertens F, Mandahl N, Mitelman F and Aman P. Fusion of the EWS and CHOP genes in myxoid liposarcoma. Oncogene. 1996; 12:489-494.

8. Rabbitts TH, Forster A, Larson R and Nathan P. Fusion of the dominant negative transcription regulator CHOP with a novel gene FUS by translocation t(12;16) in malignant liposarcoma. Nat Genet. 1993; 4:175-180.

9. Antonescu CR, Elahi A, Humphrey M, Lui MY, Healey JH, Brennan MF, Woodruff JM, Jhanwar SC and Ladanyi M. Specificity of TLS-CHOP rearrangement for classic myxoid/round cell liposarcoma: absence in predominantly myxoid well-differentiated liposarcomas. J Mol Diagn. 2000; 2:132-138.

10. Beck AH, West RB and van de Rijn M. Gene expression profiling for the investigation of soft tissue sarcoma pathogenesis and the identification of diagnostic, prognostic, and predictive biomarkers. Virchows Arch. 2010; 456:141-151.

11. Baird K, Davis S, Antonescu CR, Harper UL, Walker RL, Chen Y, Glatfelter AA, Duray PH and Meltzer PS. Gene expression profiling of human sarcomas: insights into sarcoma biology. Cancer Res. 2005; 65:9226-9235.

12. Francis P, Namlos HM, Muller C, Eden P, Fernebro J, Berner JM, Bjerkehagen B, Akerman M, Bendahl PO, Isinger A, Rydholm A, Myklebost O and Nilbert M. Diagnostic and prognostic gene expression signatures in 177 soft tissue sarcomas: hypoxia-induced transcription profile signifies metastatic potential. BMC Genomics. 2007; 8:73.

13. Renner M, Wolf T, Meyer H, Hartmann W, Penzel R, Ulrich A, Lehner B, Hovestadt V, Czwan E, Egerer G, Schmitt T, Alldinger I, Renker EK, et al. Integrative DNA methylation and gene expression analysis in high-grade soft tissue sarcomas. Genome Biol. 2013; 14:r137.

14. Singer S, Socci ND, Ambrosini G, Sambol E, Decarolis P, Wu Y, O‘Connor R, Maki R, Viale A, Sander C, Schwartz GK and Antonescu CR. Gene expression profiling of liposarcoma identifies distinct biological types/subtypes and potential therapeutic targets in well-differentiated and dedifferentiated liposarcoma. Cancer Res. 2007; 67:6626-6636.

15. Edris B, Fletcher JA, West RB, van de Rijn M and Beck AH. Comparative gene expression profiling of benign and malignant lesions reveals candidate therapeutic compounds for leiomyosarcoma. Sarcoma. 2012; 2012:805614.

16. Casali PG, Sanfilippo R and D‘Incalci M. Trabectedin therapy for sarcomas. Curr Opin Oncol. 2010; 22:342-346.

17. Germano G, Frapolli R, Simone M, Tavecchio M, Erba E, Pesce S, Pasqualini F, Grosso F, Sanfilippo R, Casali PG, Gronchi A, Virdis E, Tarantino E, et al. Antitumor and anti-inflammatory effects of trabectedin on human myxoid liposarcoma cells. Cancer Res. 2010; 70:2235-2244.

18. Davis IJ, McFadden AW, Zhang Y, Coxon A, Burgess TL, Wagner AJ and Fisher DE. Identification of the receptor tyrosine kinase c-Met and its ligand, hepatocyte growth factor, as therapeutic targets in clear cell sarcoma. Cancer Res. 2010; 70:639-645.

19. Butrynski JE, D‘Adamo DR, Hornick JL, Dal Cin P, Antonescu CR, Jhanwar SC, Ladanyi M, Capelletti M, Rodig SJ, Ramaiya N, Kwak EL, Clark JW, Wilner KD, et al. Crizotinib in ALK-rearranged inflammatory myofibroblastic tumor. N Engl J Med. 2010; 363:1727-1733.

20. Judson I. Targeted therapies in soft tissue sarcomas. Ann Oncol. 2010; 21 Suppl 7:vii277-280.

21. Taylor JGt, Cheuk AT, Tsang PS, Chung JY, Song YK, Desai K, Yu Y, Chen QR, Shah K, Youngblood V, Fang J, Kim SY, Yeung C, et al. Identification of FGFR4-activating mutations in human rhabdomyosarcomas that promote metastasis in xenotransplanted models. J Clin Invest. 2009; 119:3395-3407.

22. Fernanda Amary M, Ye H, Berisha F, Khatri B, Forbes G, Lehovsky K, Frezza AM, Behjati S, Tarpey P, Pillay N, Campbell PJ, Tirabosco R, Presneau N, et al. Fibroblastic growth factor receptor 1 amplification in osteosarcoma is associated with poor response to neo-adjuvant chemotherapy. Cancer medicine. 2014; 3:980-987.

23. Zhang K, Chu K, Wu X, Gao H, Wang J, Yuan YC, Loera S, Ho K, Wang Y, Chow W, Un F, Chu P and Yen Y. Amplification of FRS2 and activation of FGFR/FRS2 signaling pathway in high-grade liposarcoma. Cancer Res. 2013; 73:1298-1307.

24. Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GA and Zumla A. The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal Biochem. 2005; 344:141-143.

25. Fassunke J, Blum MC, Schildhaus HU, Zapatka M, Brors B, Kunstlinger H, Buttner R, Wardelmann E and Merkelbach-Bruse S. qPCR in gastrointestinal stromal tumors: Evaluation of reference genes and expression analysis of KIT and the alternative receptor tyrosine kinases FLT3, CSF1-R, PDGFRB, MET and AXL. BMC Mol Biol. 2010; 11:100.

26. Huggett J, Dheda K, Bustin S and Zumla A. Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005; 6:279-284.

27. Bustin SA, Benes V, Nolan T and Pfaffl MW. Quantitative real-time RT-PCR—a perspective. J Mol Endocrinol. 2005; 34:597-601.

28. Pao W, Miller V, Zakowski M, Doherty J, Politi K, Sarkaria I, Singh B, Heelan R, Rusch V, Fulton L, Mardis E, Kupfer D, Wilson R, et al. EGF receptor gene mutations are common in lung cancers from “never smokers” and are associated with sensitivity of tumors to gefitinib and erlotinib. Proc Natl Acad Sci U S A. 2004; 101:13306-13311.

29. Tuveson DA, Willis NA, Jacks T, Griffin JD, Singer S, Fletcher CD, Fletcher JA and Demetri GD. STI571 inactivation of the gastrointestinal stromal tumor c-KIT oncoprotein: biological and clinical implications. Oncogene. 2001; 20:5054-5058.

30. Turner N and Grose R. Fibroblast growth factor signaling: from development to cancer. Nat Rev Cancer. 2010; 10:116-129.

31. Di Martino E, Tomlinson DC and Knowles MA. A Decade of FGF Receptor Research in Bladder Cancer: Past, Present, and Future Challenges. Adv Urol. 2012; 2012:Article ID 429213.

32. Byron SA, Gartside MG, Wellens CL, Goodfellow PJ, Birrer MJ, Campbell IG and Pollock PM. FGFR2 mutations are rare across histologic subtypes of ovarian cancer. Gynecol Oncol. 2010; 117:125-129.

33. Kunii K, Davis L, Gorenstein J, Hatch H, Yashiro M, Di Bacco A, Elbi C and Lutterbach B. FGFR2-amplified gastric cancer cell lines require FGFR2 and Erbb3 signaling for growth and survival. Cancer Res. 2008; 68:2340-2348.

34. Pardo OE, Latigo J, Jeffery RE, Nye E, Poulsom R, Spencer-Dene B, Lemoine NR, Stamp GW, Aboagye EO and Seckl MJ. The fibroblast growth factor receptor inhibitor PD173074 blocks small cell lung cancer growth in vitro and in vivo. Cancer Res. 2009; 69:8645-8651.

35. Tomlinson DC, Lamont FR, Shnyder SD and Knowles MA. Fibroblast growth factor receptor 1 promotes proliferation and survival via activation of the mitogen-activated protein kinase pathway in bladder cancer. Cancer Res. 2009; 69:4613-4620.

36. Weiss J, Sos ML, Seidel D, Peifer M, Zander T, Heuckmann JM, Ullrich RT, Menon R, Maier S, Soltermann A, Moch H, Wagener P, Fischer F, et al. Frequent and focal FGFR1 amplification associates with therapeutically tractable FGFR1 dependency in squamous cell lung cancer. Sci Transl Med. 2010; 2:62ra93.

37. Guagnano V, Furet P, Spanka C, Bordas V, Le Douget M, Stamm C, Brueggen J, Jensen MR, Schnell C, Schmid H, Wartmann M, Berghausen J, Drueckes P, et al. Discovery of 3-(2,6-dichloro-3,5-dimethoxy-phenyl)-1-{6-[4-(4-ethyl-piperazin-1-yl)-phenylamin o]-pyrimidin-4-yl}-1-methyl-urea (NVP-BGJ398), a potent and selective inhibitor of the fibroblast growth factor receptor family of receptor tyrosine kinase. J Med Chem. 2011; 54:7066-7083.

38. Gavine PR, Mooney L, Kilgour E, Thomas AP, Al-Kadhimi K, Beck S, Rooney C, Coleman T, Baker D, Mellor MJ, Brooks AN and Klinowska T. AZD4547: an orally bioavailable, potent, and selective inhibitor of the fibroblast growth factor receptor tyrosine kinase family. Cancer Res. 2012; 72:2045-2056.

39. Girnita L, Girnita A, Wang M, Meis-Kindblom JM, Kindblom LG and Larsson O. A link between basic fibroblast growth factor (bFGF) and EWS/FLI-1 in Ewing’s sarcoma cells. Oncogene. 2000; 19:4298-4301.

40. Ishibe T, Nakayama T, Okamoto T, Aoyama T, Nishijo K, Shibata KR, Shima Y, Nagayama S, Katagiri T, Nakamura Y, Nakamura T and Toguchida J. Disruption of fibroblast growth factor signal pathway inhibits the growth of synovial sarcomas: potential application of signal inhibitors to molecular target therapy. Clin Cancer Res. 2005; 11:2702-2712.

41. Schweiger N, Hauck M, Steinhoff H, Sampl S, Reifinger M, Walter I, Kreilmeier T, Marian B, Grusch M, Berger W, Holzmann K and Kleiter M. Canine and human sarcomas exhibit predominant FGFR1 expression and impaired viability after inhibition of signaling. Molecular carcinogenesis. 2014:In press.

42. Patel NG, Kumar S and Eggo MC. Essential role of fibroblast growth factor signaling in preadipoctye differentiation. J Clin Endocrinol Metab. 2005; 90:1226-1232.

43. Kilvaer TK, Valkov A, Sorbye SW, Smeland E, Bremnes RM, Busund LT and Donnem T. Fibroblast growth factor 2 orchestrates angiogenic networking in non-GIST STS patients. J Transl Med. 2011; 9:104.

44. Allerstorfer S, Sonvilla G, Fischer H, Spiegl-Kreinecker S, Gauglhofer C, Setinek U, Czech T, Marosi C, Buchroithner J, Pichler J, Silye R, Mohr T, Holzmann K, et al. FGF5 as an oncogenic factor in human glioblastoma multiforme: autocrine and paracrine activities. Oncogene. 2008; 27:4180-4190.

45. Antoine M, Wirz W, Tag CG, Mavituna M, Emans N, Korff T, Stoldt V, Gressner AM and Kiefer P. Expression pattern of fibroblast growth factors (FGFs), their receptors and antagonists in primary endothelial cells and vascular smooth muscle cells. Growth Factors. 2005; 23:87-95.

46. Edgar R, Domrachev M and Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002; 30:207-210.

47. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A and Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002; 3:RESEARCH0034.

48. Pfaffl MW, Horgan GW and Dempfle L. Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002; 30:e36.


Creative Commons License All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.
PII: 4046