Model-based extension of high-throughput to high-content data
- Andrea C Pfeifer†^{1, 2},
- Daniel Kaschek†^{3},
- Julie Bachmann^{1},
- Ursula Klingmüller^{1, 2} and
- Jens Timmer^{3, 4}Email author
https://doi.org/10.1186/1752-0509-4-106
© Pfeifer et al; licensee BioMed Central Ltd. 2010
Received: 30 April 2010
Accepted: 5 August 2010
Published: 5 August 2010
Abstract
Background
High-quality quantitative data is a major limitation in systems biology. The experimental data used in systems biology can be assigned to one of the following categories: assays yielding average data of a cell population, high-content single cell measurements and high-throughput techniques generating single cell data for large cell populations. For modeling purposes, a combination of data from different categories is highly desirable in order to increase the number of observable species and processes and thereby maximize the identifiability of parameters.
Results
In this article we present a method that combines the power of high-content single cell measurements with the efficiency of high-throughput techniques. A calibration on the basis of identical cell populations measured by both approaches connects the two techniques. We develop a mathematical model to relate quantities exclusively observable by high-content single cell techniques to those measurable with high-content as well as high-throughput methods. The latter are defined as free variables, while the variables measurable with only one technique are described in dependence of those. It is the combination of data calibration and model into a single method that makes it possible to determine quantities only accessible by single cell assays but using high-throughput techniques. As an example, we apply our approach to the nucleocytoplasmic transport of STAT5B in eukaryotic cells.
Conclusions
The presented procedure can be generally applied to systems that allow for dividing observables into sets of free quantities, which are easily measurable, and variables dependent on those. Hence, it extends the information content of high-throughput methods by incorporating data from high-content measurements.
Background
In systems biology, a wide range of experimental data is used for mathematical modeling. Qualitative data mostly serves as a basis for determining network structures, whereas dynamic pathway modeling relies on high-quality quantitative data. In general, experimental data describing biological systems can be divided into three groups. Firstly, data generated from large cell populations yields an average information of the whole population behavior. However, cell population assays such as biochemical measurements or microarray studies can be misleading as large cell-to-cell variations are often observed, even in seemingly uniform populations. This stochasticity can be caused by asynchronous cell cycles, differences in cell sizes and varying protein states or expression levels [1–3]. Secondly, single cell data with high-content information from a limited number of cells result in a stochastic distribution of measured quantities. Many single cell approaches are based on microscopy, but other technologies are under development to investigate for example gene expression or proteins in single cells [4–6]. The third group covers a small range of experimental techniques that generate single cell data from large cell populations in a high-throughput format. Most common among those is flow cytometry, which however is limited to measurements from cells in suspension. Moreover, in contrast to microscopy, standard flow cytometry can only detect average whole cell fluorescence intensities lacking spatially resolved information. Currently, high-throughput imaging techniques as well as imaging flow cytometers digitally imaging cells directly in flow are being developed, with the goal to gather high-content information from a large number of single cells [7, 8]. This will increase the number of parameters that can be determined in parallel by high-throughput and high-content techniques.
For modeling purposes it is essential to link data from different types of experiments in order to include as many details of the system as possible in the modeling process and to avoid non-identifiabilities during the parameter estimation. However, some of the components can only be measured by time consuming high-content techniques. For models describing entire cell populations, high-content data for large cell numbers is necessary but often impossible to provide. In contrast, high-throughput techniques can generate these large data sets, despite a lack in detailed single cell information.
A signaling pathway that has been extensively investigated by dynamic pathway modeling is the JAK-STAT pathway [9]. Upon binding of an extracellular ligand to the respective receptor latent signal transducers and activators of transcription (STATs) are activated by Janus kinases (JAK) leading to rapid nucleocytoplasmic cycling of STATs. In addition, constitutive nucleocytoplasmic cycling of unphosphorylated STAT has been shown for several STAT proteins by biochemical and microscopic experiments [10–15]. It has been proposed that import of STAT is enhanced upon activation [16], while export of activated STAT is slowed down either through retention in the nucleus by DNA binding [17] or possibly a different export mechanism [15]. Previously, rapid nucleocytoplasmic cycling of activated STAT5 has been identified as the step most sensitive to perturbation within the core module of the JAK2/STAT5 pathway by mathematical modeling based on biochemical data [18], but import and export rates could not be measured experimentally. These transport steps are crucial as important decisions regulating cell fate are made by the nuclear reactions of STATs.
A method to determine the rates for nuclear import and export of STAT5 is fluorescence recovery after photobleaching (FRAP). FRAP is a single cell fluorescence microscopy method routinely used to measure the kinetics of transport processes between cell compartments as well as diffusion and dynamic binding reactions [19, 20]. One prerequisite for a quantitative FRAP experiment is that the investigated system is in a steady state on the time scale of the experiment otherwise a mathematical description of the data is difficult to obtain. The JAK-STAT system is only in a steady state in unstimulated cells, ligand stimulation induces phosphorylation of STATs and thereby perturbs the steady state. Therefore, we focused on the nuclear import and export rates of unphosphorylated STAT5 with the goal to generate rates for the steady state in unstimulated cells that can be set to a fixed value in a larger pathway model. Biochemical data describing the phosphorylation dynamics of the pathway components after stimulation in combination with mathematical modeling can then serve to indirectly determine nuclear import and export rates for phosphorylated STAT5.
Here, we present a model for extracting the import and export rates from FRAP experiments of STAT5B-GFP in the steady state of unstimulated NIH3T3-EpoR cells. Furthermore, the dependence of these rates on STAT5B-GFP concentration and cell size is shown. To be able to combine this information with biochemical data from cell populations expressing STAT5B-GFP, cell size distribution and STAT5B-GFP concentration distribution within the respective cell population are additionally measured by flow cytometry. Cell size as well as STAT5B-GFP concentration are estimated directly from flow cytometry data after calibration of these data to microscopy data.
The calibration procedure can be generally applied to link data from powerful high-content techniques and fast, efficient high-throughput methods. In combination with the mathematical model, it provides a novel rationale to determine formerly inaccessible information for large cell populations by less time-consuming high-throughput measurements.
Results and Discussion
Data calibration links high-content with high-throughput data
i.e. X_{ T }(Y ) depends linearly on the high-content quantity X_{ C }(Y ). The slope $\frac{m}{{m}^{\prime}}$ and intercept $\frac{d-{d}^{\prime}}{{m}^{\prime}}$ of eq. (4) need to be determined in order to translate X_{ C }into X_{ T }and vice versa. For this purpose, the distribution quantiles of X_{ C }and X_{ T }are used.
Assuming that N_{ C }and N_{ T }> N_{ C }measurements have been performed for the high-content and high-throughput techniques, respectively, the ordered set of measurements {X_{ C,i }}_{i = 1},...,N_{ C }is an estimate of the N_{ C }equally spaced quantiles ${\left\{{\stackrel{~}{X}}_{C,i}\right\}}_{i=1,...,{N}_{C}}$ of the theoretical distribution of X_{ C }. In the same way, the sample quantiles ${\left\{{\stackrel{~}{X}}_{T,i}\right\}}_{i=1,...,{N}_{C}}$of {X_{ T,i }}_{i = 1},...,N_{ T }estimate the N_{ C }theoretical quantiles of X_{ T }. According to eq. (4) the distributions of X_{ T }and X_{ C }belong to the same location-scale family. Consequently, the QQ-plot of X_{ T }versus X_{ C }is supposed to follow a straight line with intercept $\frac{d-{d}^{\prime}}{{m}^{\prime}}$ and slope $\frac{m}{{m}^{\prime}}$. A least squares fit of the QQ-plot gives asymptotically unbiased estimates of slope and intercept for a large class of theoretical distributions. The convergence of the sample quantiles to the theoretical quantiles as well as the convergence of the least squares estimator is well known and is carried out rigorously in [21].
An implementation of this calibration procedure is provided by the R script [Additional file 1] in the supplement. A sample configuration is given by [Additional file 2].
Nucleocytoplasmic cycling of STAT5B is modeled as saturatable pump
As an example, the analysis of nuclear import and export of the transcription factor STAT5B was chosen. Nucleocytoplasmic cycling is only measurable by single cell microscopy, namely FRAP, whereas other crucial features such as the dynamic changes of the phosphorylation state of the proteins are accessed by biochemical measurements from cell populations. In eukaryotic cells, the nucleus is separated from the cytoplasm by the nuclear envelope. Molecules can only migrate between those two compartments through nuclear pores forming small holes in the membrane. Small molecules (<20-40 kDa) can diffuse freely through nuclear pores whereas larger molecules require active transport aided by soluble transport proteins that interact with the cargo molecule as well as the nuclear pore. Active nuclear import and export are regulated by different mechanisms. In order to be imported into the nucleus, proteins usually carry a nuclear localization signal (NLS) to which importins can bind and enable nuclear translocation. Similarly, a nuclear export signal (NES) within the cargo protein structure is recognized by an exportin. For most proteins of the STAT family, the respective importins have been identified (reviewed in [22]). In the case of STAT5B however, so far no importins could be identified that directly interact with the transcription factor [23]. Instead, import of STAT5B has been suggested to require additional factors acting as chaperones between the importins and STAT5B [24]. Active nuclear export of STATs is generally mediated by the exportin CRM1.
For large concentrations (c ≫ γ) I saturates with saturation value β. For small concentrations (c ≪ γ) I depends linearly on c with slope $\frac{\beta}{\gamma}$.
Equation (5) is the resulting current for all pores of a cell. The constants β and γ may still vary within a population, i.e. from cell to cell. In a next step the saturation value β = _{ k }K as well as the slope $\frac{\beta}{\gamma}={\kappa}^{\prime}K$ are assumed to depend linearly on a quantity K which is the product of the abundance of transport factors and the number of nuclear pores. This is appropriate for two reasons:
- 1.
K is dominated by the number of nuclear pores which have a similar density throughout the cell population. Hence, K depends linearly on the nuclear surface area A _{ nuc }.
- 2.
K is dominated by the number of cytoplasmic transport factors with the same concentration in all cells which is proportional to the cytoplasmic volume V _{ cyt }.
- 3.
K is dominated by the number of nuclear transport factors with the same concentration in all cells which is proportional to the nuclear volume V _{ nuc }.
demonstrates the concentration dependency of the normalized transport currents under the assumption that the parameters β_{ l }and γ_{ l }are constant throughout the population. This assumption is necessary for a valid formula describing the import and export currents within a population. The second formulation with α_{ l }(c) follows the idea of a linearly increasing current for small concentrations and will also be used.
Import and export current distribution for STAT5B
Import and export currents depend on STAT5B concentration and cell size
To determine the import and export rates α_{ imp }and α_{ exp }39 FRAP data sets generated from cells expressing varying concentrations of STAT5B-GFP were fitted with eq. (23) described in the Methods section. Variable protein levels were achieved by a tightly regulatable expression system that we developed based on a Tet-inducible promoter.
Estimated Parameters
β[mol/s] | γ[μ M] | |
---|---|---|
I _{ imp } | 18.04 ± 6.55 | 3.56 ± 2.63 |
I _{ exp } | 19.01 ± 7.06 | 6.49 ± 4.71 |
β_{ A }[mM/s] | γ [μ M] | |
J _{ imp,A } | 16.13 ± 3.88 | 3.12 ± 1.63 |
J _{ exp,A } | 74.72 ± 22.21 | 7.07 ± 3.96 |
β_{ V }[mM μ m/s] | γ [μ M] | |
${J}_{imp,{V}_{nuc}}$ | 86.37 ± 18.82 | 2.98 ± 1.44 |
${J}_{exp,{V}_{nuc}}$ | 40.13 ± 10.54 | 6.81 ± 3.42 |
${J}_{imp,{V}_{cyt}}$ | 22.14 ± 5.30 | 3.43 ± 1.71 |
${J}_{exp,{V}_{cyt}}$ | 8.28 ± 2.37 | 5.31 ± 3.25 |
The significance of the χ^{2} reduction has been tested with a bootstrap method: from the 39 data points 39 points have been drawn randomly with replacement. Then for all models, i.e. without normalization, K ∝ A, K ∝ V_{ nuc }and K ∝ V_{ cyt }pairwise differences ${\delta}_{ij}={\chi}_{i}^{2}-{\chi}_{j}^{2}$of the χ^{2} values have been computed leading to 6 difference values for the import and 6 difference values for the export models. This procedure has been repeated 10^{4} times resulting in 2 × 6 distributions ${\delta}_{ij}^{(exp)}$ and ${\delta}_{ij}^{(imp)}$of χ^{2} difference values. The position of zero with respect to such a distribution decides whether one of the compared models is superior to the other. More precisely:
For the export distributions model 3 - normalization by nucleus volume - is superior to all other models at a 99% confidence level (3σ). For the import data the situation is not so clear. Models 3 and 4 cannot be discriminated and seem to describe the data equally well. Both models are clearly superior to model 1 and exceed model 2 at a 1σ level.
We decided to follow the hypothesis of normalization by the volumes of the originating compartments, i.e. export model 3 and import model 4. From a biological point of view this seems to be the most reasonable hypothesis. From a practical point of view, models 3 and 4 describe the import equally well and cannot be distinguished given the data at hand.
Calibration of flow cytometry data to microscopy data yields comparable quantities
To yield comparable distributions, it has to be considered that different quantities are measured by the two techniques. Microscopy data directly result in absolute numbers for compartment volumes and protein concentrations, whereas flow cytometry data from the fluorescence intensity channel F_{2} are logarithmic due to the amplification of the signal by the instrument and the forward scatter of a flow cytometer using light scattering is an approximate measure of the cell cross-section area (see [27]). Therefore, values from either flow cytometry or microscopy measurements had to be transformed to yield comparable quantities. For practical reasons, the microscopy fluorescence intensities I_{ micro }were logarithmized yielding X_{ M }= log I_{ micro }. Similarly, cell volumes V_{ micro }determined by microscopy were converted to cross-section areas assuming a spherical shape of the cell as is the case for flow cytometry samples. This yields ${X}_{M}=\sqrt[3]{36\pi}\xb7{V}_{micro}^{\frac{2}{3}}$. For flow cytometry, X_{ F }= F_{2} and X_{ F }= F_{0} for fluorescence intensity and cross section area respectively. X_{ M }and X_{ F }defined like this build the basis for the calibration method described above.
where p ∈ [0, 1] is the fraction of the population that is dropped and n = 2 for quadratic penalization. The penalization is chosen on purpose to fulfill ⟨(n + 1)p^{ n })⟩ = 1 for uniformly distributed p. This guarantees that the penalization is of the same magnitude as χ^{2}. The resulting X^{2} curve for the size distribution indicated that only the complete flow cytometry data set lead to the best accordance, while a local minimum existed for the fluorescence intensity distribution (fig. 4B).
After choosing the optimal cut-off, a least squares regression was applied to the QQ-plot. The linearity of the data points confirmed that the shapes of the two distributions are the same. However, even after two cuts there were deviations for the border points that result from a small population of cells which is detected differently by flow cytometry and by microscopy. To exclude biased fit parameters the least squares regression was restricted to the inner 66% region of points (fig. 4C).
Thus, data preprocessing and subsequent least squares regression of the QQ-plot lead to comparable quantities obtained by different experimental techniques. All functions for preprocessing the flow cytometry data and for calibration of flow cytometry to microscopy data are included in the R script 3.1 provided in the supplement.
Distributions of transport currents for an exemplary cell population are calculated
In order to compute the distribution of currents for a sample flow cytometry measurement, the calibration was combined with the formula describing the currents (eq. (8)). As has been shown above (fig. 2D), the rates α_{ in }and α_{ out }depend on the STAT5B concentrations C_{ cyt }in the cytoplasm and C_{ nuc }in the nucleus as well as the compartment volumes V_{ nuc }and V_{ cyt }for nuclear export and import, respectively.
Since the individual cell compartments cannot be distinguished by flow cytometry, an average ratio of the cytoplasmic to nuclear quantities had to be estimated from microscopy data. For every FRAP data set, the fractions ${f}_{V}=\frac{{V}_{cyt}}{{V}_{nuc}}$ and ${f}_{c}=\frac{{c}_{cyt}}{{c}_{nuc}}$ were determined and averaged. In addition, we tested if the fraction f_{ V }and the cell volume V or f_{ C }and the total concentration c of STAT5B-GFP are correlated. A large correlation value would indicate that an additional model for describing the dependency of the compartment quantities on the overall quantities would be necessary. The data lead to f_{ V }= 4.27 ± 0.11 and f_{ c }= 0.645 ± 0.015 and the correlation test revealed cor(f_{ V }, V ) = 0.15 ± 0.28 and cor(f_{ c }, c) = -0.37 ± 0.26. Thus, the assumption of a correlation for the cell volume would not lead to a better estimate of f_{ V }Even for f_{ c }considering the correlation would have a minor effect.
Conclusion
In this study, transport rates for unphosphorylated STAT5B were determined in single cells by FRAP and found to follow saturation kinetics dependent on both STAT5B-GFP expression level as well as size of the originating compartment. This reflects a saturation of cofactors necessary for active transport of STAT5B through the nuclear pore complex. The parameters for concentration and volume dependency of the cycling currents were estimated. To predict transport currents with the saturation model for large cell populations, STAT5B-GFP concentration and cell size distribution were measured by flow cytometry. As flow cytometry only yields relative values for cell size and total cell fluorescence, a calibration to absolute numbers generated by single cell microscopy is required. For calibration, the concentration of the transcription factor STAT5B as well as cell size were determined by confocal microscopy and flow cytometry from the identical cell population, resulting in a linear calibration curve. Subsequently, absolute cell size and STAT5B-GFP concentration distributions were computed from flow cytometry data using the calibration curve. Finally, transport current distributions and thereby cell-to-cell variation were predicted using the saturation model.
In recent years, other members of the STAT protein family have been studied by FRAP to investigate nucleocytoplasmic cycling [12, 15, 28], but the data have so far only been interpreted qualitatively. Our results provide a procedure to link directly measured import and export rates of unphosphorylated STAT5B with data indirectly describing the nucleocytoplasmic cycling of activated STAT5B generated by biochemical experiments. Furthermore, by using an inducible expression system for STAT5B-GFP, we identify a saturation-like behavior of STAT5B nuclear import and export, indicating a limitation in transport factors. The nature of these factors remains to be identified.
The proposed method is generic and is applicable as long as two conditions are fulfilled. First, the quantities that are measured by a certain high-content method have to be functionally related. This allows for expressing a subset of the quantities, defined as the dependent variables, as a function of the remaining, free variables. Second, the free variables have to be part of the quantities that can be measured by a given high-throughput method. If these conditions hold, it is possible to determine the function connecting free and dependent variables by setting up a mathematical model and estimating its parameters. Furthermore, it is possible to calibrate the two measurement techniques against each other as the high-throughput quantities are in particular part of the high-content quantities. This means that a high-throughput measurement can be translated into the ambit of a high-content measurement. Applying the fixed parameter model then leads to a prediction of the dependent variables' distributions representing an indirect determination of these variables for every cell of the population. The method is especially useful if there is a great discrepancy in accessibility between free and dependent variables. This combination of two experimental approaches results in a higher degree of measured variables suitable for mathematical modeling and a reduction of non-identifiabilities in the parameter estimation.
Methods
Experimental Procedures
The retroviral expression vector pMOWS containing the cDNA for murine HA-EpoR was introduced into NIH3T3 cells (ATCC) and a single cell clone stably expressing HA-EpoR was obtained by selection with G418. pMOWSIN-TREt-STAT5B-GFP was cotransduced into NIH3T3-EpoR cells together with the cDNA for the transactivator protein contained in pMOWS-rtTAM2. A single cell clone stably expressing murine STAT5B-GFP was obtained by selection with puromycin. Expression of STAT5BGFP was regulated by a Tet-inducible promoter included in pMOWSIN-TREt. pMOWSIN-TREt was generated by digesting pTRE-tight (Clontech) and inserting TREt into the self-inactivating (SIN) retroviral vector pMOWSIN. pMOWS-rtTAM2 was generated by introducing cDNA of rtTAM2 from pUHrT-62-1 (H. Bujard, Heidelberg, Germany) into pMOWS using BamHI/EcoRI restriction sites [29]. To simplify identification of the nuclei cells used for FRAP experiments also were transduced with pMOWS-H2B-mCherry. All cells were maintained in DMEM supplemented with 10% calf serum and 1% PenStrep.
To determine import and export currents, STAT5B-GFP was photobleached in the entire nuclear region with 100% laser power (488 nm). For analysis of the transport dynamics 10 prebleach and approximately 240 postbleach images of the whole cell were acquired for 30-40 min after bleaching.
Flow cytometry analysis of STAT5B-GFP expression level and measurement of the approximate cell size were performed on a BD FACSCalibur system with the software package CellQuest. Cells were grown in 60 mm cell culture dishes and were treated as for microscopy. Cells were detached from the dishes by 0.05% trypsin/EDTA and washed once in PBS/0.3% BSA. For each cell population 20 000 cells were measured. Forward and side scatter were detected linearly, for fluorescence intensity detection the signal was logarithmically amplified. NIH3T3-EpoR cells were used as control for cellular autofluorescence and cell size. Raw data was extracted from CellQuest files with FCSExtract [31]. Fluorescence intensity values were directly used for analysis. Values for the forward scatter were assumed to be approximately proportional to the cross-section area of the cell [27]. Cell shape was assumed to be roughly spherical for detached cells and therefore the relation between cross-section area and volume is known.
For the calibration measurement, cells from one cell population were seeded in 60 mm dishes as well as Labtek chambers 20 hours before the experiment. STAT5B-GFP expression was induced with either 10 or 250 ng/ml doxycycline 16 hours prior to serum-starvation. Flow cytometry analysis was performed as described above. z-stacks of 100 tiled frames were acquired by confocal microscopy. For each doxycyline treatment the cell volume and the total amount of STAT5B-GFP per cell were determined for 200 cells as described for FRAP experiments above.
Mathematical Model
Import and export currents from FRAP data
As in the previous section K = A, V_{ cyt }, V_{ nuc }accounts for the normalization. The associated index l is omitted as an index of j. The ε-terms describe the continuous bleaching due to constant laser exposition during postbleach image acquisition.
with ${\lambda}_{2}=\frac{{\alpha}_{imp}}{{V}_{cyt}}+\frac{{\alpha}_{exp}}{{V}_{nuc}}$.
Consequently, the experimentally accessible quantities cyt and nuc are directly associated to the concentrations appearing in the ODE system. Note that the exponential decrease of the signal (due to continuous bleaching) and the proportionality factor between the signal S and the concentration c_{ L }drop out. This is even true if the proportionality factor is time dependent.
The resulting rate function α_{ l }(c) or equivalently j_{ l }(c) = c α_{ l }(c) can be used for microscopy data: cell images are analyzed for the quantity K and for fluorescence intensities which allow calculating the protein concentrations of interest. Plugging these values in the formula for j_{ l }(c) yields an estimate for the current between nucleus and cytoplasm of the investigated cell without measuring it explicitly.
Notes
Declarations
Acknowledgements
The authors thank Hermann Bujard for the generous gift of the cDNA for rtTAM2. We also thank Clemens Kreutz and Verena Becker for valuable discussions and critically reading the manuscript.
This work was supported by the German Federal Ministry of Education and Research (BMBF) grants FORSYS-ViroQuant (#0313923), LungSys (FKZ0315415E) and the Excellence Initiative of the German Federal and State Governments.
Authors’ Affiliations
References
- Sigal A, Milo R, Cohen A, Geva-Zatorsky N, Klein Y, Liron Y, Rosenfeld N, Danon T, Perzov N, Alon U: Variability and memory of protein levels in human cells. Nature. 2006, 444 (7119): 643-646. 10.1038/nature05316View ArticlePubMedGoogle Scholar
- Spencer SL, Gaudet S, Albeck JG, Burke JM, Sorger PK: Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis. Nature. 2009, 459 (7245): 428-432. 10.1038/nature08012PubMed CentralView ArticlePubMedGoogle Scholar
- Halter M, Elliott JT, Hubbard JB, Tona A, Plant AL: Cell volume distributions reveal cell growth rates and division times. J Theor Biol. 2009, 257: 124-130. 10.1016/j.jtbi.2008.10.031View ArticlePubMedGoogle Scholar
- Diercks A, Kostner H, Ozinsky A: Resolving cell population heterogeneity: real-time PCR for simultaneous multiplexed gene detection in multiple single-cell samples. PLoS One. 2009, 4 (7): e6326- 10.1371/journal.pone.0006326PubMed CentralView ArticlePubMedGoogle Scholar
- Masujima T: Live single-cell mass spectrometry. Anal Sci. 2009, 25 (8): 953-960. 10.2116/analsci.25.953View ArticlePubMedGoogle Scholar
- Lapainis T, Rubakhin SS, Sweedler JV: Capillary electrophoresis with electrospray ionization mass spectrometric detection for single-cell metabolomics. Anal Chem. 2009, 81 (14): 5858-5864. 10.1021/ac900936gPubMed CentralView ArticlePubMedGoogle Scholar
- Pepperkok R, Ellenberg J: High-throughput fluorescence microscopy for systems biology. Nat Rev Mol Cell Biol. 2006, 7 (9): 690-696. 10.1038/nrm1979View ArticlePubMedGoogle Scholar
- George T, Morrissey P, Cui C, Singh S, Fitzgerald-Bocarsly P: Measurement of cytoplasmic to nuclear translocation. Curr Protoc Cytom. 2009Google Scholar
- Pfeifer AC, Timmer J, Klingmüller U: Systems biology of JAK/STAT signalling. Essays Biochem. 2008, 45: 109-120. 10.1042/BSE0450109View ArticlePubMedGoogle Scholar
- Meyer T, Begitt A, Lödige I, van Rossum M, Vinkemeier U: Constitutive and IFN-gamma-induced nuclear import of STAT1 proceed through independent pathways. EMBO J. 2002, 21 (3): 344-354. 10.1093/emboj/21.3.344PubMed CentralView ArticlePubMedGoogle Scholar
- Bhattacharya S, Schindler C: Regulation of Stat3 nuclear export. J Clin Invest. 2003, 111 (4): 553-559.PubMed CentralView ArticlePubMedGoogle Scholar
- Pranada AL, Metz S, Herrmann A, Heinrich PC, Müller-Newen G: Real time analysis of STAT3 nucleocytoplasmic shuttling. J Biol Chem. 2004, 279 (15): 15114-15123. 10.1074/jbc.M312530200View ArticlePubMedGoogle Scholar
- Liu L, McBride KM, Reich NC: STAT3 nuclear import is independent of tyrosine phosphorylation and mediated by importin-alpha3. Proc Natl Acad Sci USA. 2005, 102 (23): 8150-8155. 10.1073/pnas.0501643102PubMed CentralView ArticlePubMedGoogle Scholar
- Zeng R, Aoki Y, Yoshida M, ichi Arai K, Watanabe S: Stat5B shuttles between cytoplasm and nucleus in a cytokine-dependent and -independent manner. J Immunol. 2002, 168 (9): 4567-4575.View ArticlePubMedGoogle Scholar
- Iyer J, Reich NC: Constitutive nuclear import of latent and activated STAT5a by its coiled coil domain. FASEB J. 2008, 22 (2): 391-400. 10.1096/fj.07-8965comView ArticlePubMedGoogle Scholar
- Marg A, Shan Y, Meyer T, Meissner T, Brandenburg M, Vinkemeier U: Nucleocytoplasmic shuttling by nucleoporins Nup153 and Nup214 and CRM1-dependent nuclear export control the subcellular distribution of latent Stat1. J Cell Biol. 2004, 165 (6): 823-833. 10.1083/jcb.200403057PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer T, Marg A, Lemke P, Wiesner B, Vinkemeier U: DNA binding controls inactivation and nuclear accumulation of the transcription factor Stat1. Genes Dev. 2003, 17 (16): 1992-2005. 10.1101/gad.268003PubMed CentralView ArticlePubMedGoogle Scholar
- Swameye I, Müller TG, Timmer J, Sandra O, Klingmüller U: Identification of nucleocytoplasmic cycling as a remote sensor in cellular signaling by databased modeling. Proc Natl Acad Sci USA. 2003, 100 (3): 1028-1033. 10.1073/pnas.0237333100PubMed CentralView ArticlePubMedGoogle Scholar
- Axelrod D, Koppel DE, Schlessinger J, Elson E, Webb WW: Mobility measurement by analysis of fluorescence photobleaching recovery kinetics. Bio-phys J. 1976, 16 (9): 1055-69.Google Scholar
- Sprague BL, McNally JG: FRAP analysis of binding: proper and fitting. Trends Cell Biol. 2005, 15 (2): 84-91. 10.1016/j.tcb.2004.12.001View ArticlePubMedGoogle Scholar
- Das B, Resnick SI: QQ - Plots, Random Sets and Data from a Heavy Tailed Distribution. Stochastic Models. 2008, 24: 103-10.1080/15326340701828308.View ArticleGoogle Scholar
- Reich NC, Liu L: Tracking STAT nuclear traffic. Nat Rev Immunol. 2006, 6 (8): 602-612. 10.1038/nri1885View ArticlePubMedGoogle Scholar
- Ma J, Cao X: Regulation of Stat3 nuclear import by importin alpha5 and importin alpha7 via two different functional sequence elements. Cell Signal. 2006, 18 (8): 1117-1126. 10.1016/j.cellsig.2005.06.016View ArticlePubMedGoogle Scholar
- Kawashima T, Bao YC, Minoshima Y, Nomura Y, Hatori T, Hori T, Fukagawa T, Fukada T, Takahashi N, Nosaka T, Inoue M, Sato T, Kukimoto-Niino M, Shirouzu M, Yokoyama S, Kitamura T: A Rac GTPase-activating protein, MgcRacGAP, is a nuclear localizing signal-containing nuclear chaperone in the activation of STAT transcription factors. Mol Cell Biol. 2009, 29 (7): 1796-1813. 10.1128/MCB.01423-08PubMed CentralView ArticlePubMedGoogle Scholar
- Timney BL, Tetenbaum-Novatt J, Agate DS, Williams R, Zhang W, Chait BT, Rout MP: Simple kinetic relationships and nonspecific competition govern nuclear import rates in vivo. J Cell Biol. 2006, 175 (4): 579-593. 10.1083/jcb.200608141PubMed CentralView ArticlePubMedGoogle Scholar
- Kopito RB, Elbaum M: Reversibility in nucleocytoplasmic transport. Proceedings of the National Academy of Sciences. 2007, 104 (31): 12743-12748. 10.1073/pnas.0702690104.View ArticleGoogle Scholar
- Hawley TS, Hawley RG: Flow Cytometry Protocols (Methods in Moldecular Biology). 2004, Humana PressView ArticleGoogle Scholar
- Köster M, Frahm T, Hauser H: Nucleocytoplasmic shuttling revealed by FRAP and FLIP technologies. Curr Opin Biotechnol. 2005, 16: 28-34. 10.1016/j.copbio.2004.11.002View ArticlePubMedGoogle Scholar
- Ketteler R, Glaser S, Sandra O, Martens UM, Klingmüller U: Enhanced transgene expression in primitive hematopoietic progenitor cells and embryonic stem cells efficiently transduced by optimized retroviral hybrid vectors. Gene Ther. 2002, 9 (8): 477-487. 10.1038/sj.gt.3301653View ArticlePubMedGoogle Scholar
- Rasband W: ImageJ 1997 - 2009., http://rsb.info.nih.gov/ij
- Glynn E: FCSExtract. 2006, http://research.stowers-institute.org/efg/ScientificSoftware/Utility/FCSExtract/index.htmGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.