James M. Piret

High-throughput microfluidic single-cell RT-qPCR

A long-sought milestone in microfluidics research has been the development of integrated technology for scalable analysis of transcription in single cells. Here we present a fully integrated microfluidic device capable of performing high-precision RT-qPCR measurements of gene expression from hundreds of single cells per run. Our device executes all steps of single-cell processing, including cell capture, cell lysis, reverse transcription, and quantitative PCR. In addition to higher throughput and reduced cost, we show that nanoliter volume processing reduced measurement noise, increased sensitivity, and provided single nucleotide specificity. We apply this technology to 3,300 single-cell measurements of (i) miRNA expression in K562 cells, (ii) coregulation of a miRNA and one of its target transcripts during differentiation in embryonic stem cells, and (iii) single nucleotide variant detection in primary lobular breast cancer cells. The core functionality established here provides the foundation from which a variety of on-chip single-cell transcription analyses will be developed.

High-throughput analysis of single hematopoietic stem cell proliferation in microfluidic cell culture arrays

Heterogeneity in cell populations poses a major obstacle to understanding complex biological processes. Here we present a microfluidic platform containing thousands of nanoliter-scale chambers suitable for live-cell imaging studies of clonal cultures of nonadherent cells with precise control of the conditions, capabilities for in situ immunostaining and recovery of viable cells. We show that this platform mimics conventional cultures in reproducing the responses of various types of primitive mouse hematopoietic cells with retention of their functional properties, as demonstrated by subsequent in vitro and in vivo (transplantation) assays of recovered cells. The automated medium exchange of this system made it possible to define when Steel factor stimulation is first required by adult hematopoietic stem cells in vitro as the point of exit from quiescence. This technology will offer many new avenues to interrogate otherwise inaccessible mechanisms governing mammalian cell growth and fate decisions.

Correlation of Murine Embryonic Stem Cell Gene Expression Profiles with Functional Measures of Pluripotency

Global gene expression profiling was performed on murine embryonic stem cells (ESCs) induced to differentiate by removal of leukemia inhibitory factor (LIF) to identify genes whose change in expression correlates with loss of pluripotency. To identify appropriate time points for the gene expression analysis, the dynamics of loss of pluripotency were investigated using three functional assays: chimeric mouse formation, embryoid body generation, and colony‐forming ability. A rapid loss of pluripotency was detected within 24 hours, with very low residual activity in all assays by 72 hours. Gene expression profiles of undifferentiated ESCs and ESCs cultured for 18 and 72 hours in the absence of LIF were determined using the Affymetrix GeneChip U74v2. In total, 473 genes were identified as significantly differentially expressed, with approximately one third having unknown biological function. Among the 275 genes whose expression decreased with ESC differentiation were several factors previously identified as important for, or markers of, ESC pluripotency, including Stat3, Rex1, Sox2, Gbx2, and Bmp4. A significant number of the decreased genes also overlap with previously published mouse and human ESC data. Furthermore, several membrane proteins were among the 48 decreased genes correlating most closely with the functional assays, including the stem cell factor receptor c‐Kit. Through identification of genes whose expression closely follows functional properties of ESCs during early differentiation, this study lays the foundation for further elucidating the molecular mechanisms regulating the maintenance of ESC pluripotency and facilitates the identification of more reliable molecular markers of the undifferentiated state.

Comprehensive microRNA expression profiling of the hematopoietic hierarchy

The hematopoietic system produces a large number of highly specialized cell types that are derived through a hierarchical differentiation process from a common stem cell population. miRNAs are critical players in orchestrating this differentiation. Here, we report the development and application of a high-throughput microfluidic real-time quantitative PCR (RT-qPCR) approach for generating global miRNA profiles for 27 phenotypically distinct cell populations isolated from normal adult mouse hematopoietic tissues. A total of 80,000 RT-qPCR assays were used to map the landscape of miRNA expression across the hematopoietic hierarchy, including rare progenitor and stem cell populations. We show that miRNA profiles allow for the direct inference of cell lineage relations and functional similarity. Our analysis reveals a close relatedness of the miRNA expression patterns in multipotent progenitors and stem cells, followed by a major reprogramming upon restriction of differentiation potential to a single lineage. The analysis of miRNA expression in single hematopoietic cells further demonstrates that miRNA expression is very tightly regulated within highly purified populations, underscoring the potential of single-cell miRNA profiling for assessing compartment heterogeneity.

Assessing differentiation status of human embryonic stem cells noninvasively using Raman microspectroscopy.

Raman microspectroscopy is an attractive approach for chemical imaging of biological specimens, including live cells, without the need for chemi-selective stains. Using a microspectrometer, near-infrared Raman spectra throughout the range 663 cm(-1) to 1220 cm(-1) were obtained from colonies of CA1 human embryonic stem cells (hESCs) and CA1 cells that had been stimulated to differentiate for 3 weeks by 10% fetal bovine serum on gelatin. Distributions and intensities of spectral bands attributed to proteins varied significantly between undifferentiated and differentiated cells. Importantly, compared to proteins and lipids, the band intensities of nucleic acids were dominant in undifferentiated cells with a dominance-reversal in differentiated cells. Thus, we could identify intensity ratios of particular protein-related bands (e.g., 757 cm(-1) tryptophan) to nucleic acid bands (784 cm(-1) DNA/RNA composite) that were effective in discriminating between spectra of undifferentiated and differentiated cells. We observed no discernible negative effects due to the laser exposure in terms of morphology, proliferation, or pluripotency of the stem cells. We conclude that Raman microscopy and complementary data processing procedures provide a rapid, noninvasive approach that can distinguish hESCs from differentiated cells. This is the first report to identify specific Raman markers for the differentiation status of hESCs.

Metabolic flux analysis of CHO cells in perfusion culture by metabolite balancing and 2D [13C, 1H] COSY NMR spectroscopy.

The physiological state of CHO cells in perfusion culture was quantified by determining fluxes through the bioreaction network using (13)C glucose and 2D-NMR spectroscopy. CHO cells were cultivated in a 2.5L perfusion bioreactor with glucose and glutamine as the primary carbon and energy sources. The reactor was inoculated at a cell density of 8 x 10(6)cells/mL and operated at approximately 10 x 10(6)cells/mL using unlabeled glucose for the first 13 days. The second phase lasted 12 days and the medium consisted of 10% [U-(13)C]glucose, 40% labeled [1-(13)C]glucose with the balance unlabeled. After the culture attained isotopic steady state, biomass samples from the last 3 days of cultivation were considered representative and used for flux estimation. They were hydrolyzed and analyzed by 2D [(13)C, (1)H] COSY measurements using the heteronuclear single quantum correlation sequence with gradients for artifacts suppression. Metabolic fluxes were determined using the 13C-Flux software package by minimizing the residuals between the experimental and the simulated NMR data. Normalized residuals exhibited a Gaussian distribution indicating good model fit to experimental data. The glucose consumption rate was 5-fold higher than that of glutamine with 41% of glucose channeled through the pentose phosphate pathway. The fluxes at the pyruvate branch point were almost equally distributed between lactate and the TCA cycle (55% and 45%, respectively). The anaplerotic conversion of pyruvate to oxaloacetate by pyruvate carboxylase accounted for 10% of the pyruvate flux with the remaining 90% entering the TCA cycle through acetyl-CoA. The conversion of malate to pyruvate catalyzed by the malic enzyme was 70% higher than that for the anaplerotic reaction catalyzed by pyruvate carboxylase. Most amino acid catabolic and biosynthetic fluxes were significantly lower than the glycolytic and TCA cycle fluxes. Metabolic flux data from NMR analysis validated a simplified model where metabolite balancing was used for flux estimation. In this reduced flux space, estimates from these two methods were in good agreement. This simplified model can routinely be used in bioprocess development experiments to estimate metabolic fluxes with much reduced analytical investment. The high resolution flux information from 2D-NMR spectroscopy coupled with the capability to validate a simplified metabolite balancing based model for routine use make (13)C-isotopomer analysis an attractive bioprocess development tool for mammalian cell cultures.

Mammalian cell retention devices for stirred perfusion bioreactors

Within the spectrum of current applications for cell culture technologies, efficient large-scale mammalian cell production processes are typically carried out in stirred fed-batch or perfusion bioreactors. The specific aspects of each individual process that can be considered when determining the method of choice are presented. A major challenge for perfusion reactor design and operation is the reliability of the cell retention device. Current retention systems include cross-flow membrane filters, spin-filters, inclined settlers, continuous centrifuges and ultrasonic separators. The relative merits and limitations of these technologies for cell retention and their suitability for large-scale perfusion are discussed.

Limitations to the amplification and stability of human tissue-type plasminogen activator expression by Chinese hamster ovary cells

Chinese hamster ovary cell production of recombinant tissue-type plasminogen activator (t-PA) was increased by amplification of cotransfected dihydrofolate reductase cDNA using stepwise adaptation to increasing methotrexate (MTX) concentrations. The highest producing clones were isolated at 5 microM MTX and yielded 26,000 U/10(6) cells/day t-PA (43 microgram/10(6) cells/day). Above 25 microM MTX, cell specific t-PA production rates became increasingly variable and the cDNA copynumbers decreased. No apparent correlation between the cell specific t-PA production rate and the growth rate was observed upon subcloning of the amplified cells. When MTX selection was removed, the t-PA production rate decreased up to tenfold within 40 days; this was accompanied by an up to 60% drop in cDNA copynumber. Subclones isolated after 108 days of culture in the absence of MTX were, on average, sixfold more stable than their parental cells. In culture without MTX, the maximum stable t-PA production rate obtained (over 250 days) was 7000 +/- 750 U/10(6) cells/day (approximately 12 microgram/10(6) cells/day), approximately threefold lower than the maximum unstable levels of production reached under selective pressure. Taken together, these results define a wide range of the highest t-PA expression rates obtained under MTX selection, for which stable expression without selection has not been reported.

Characterization of polyhormonal insulin-producing cells derived in vitro from human embryonic stem cells.

Human embryonic stem cells (hESCs) were used as a model system of human pancreas development to study characteristics of the polyhormonal cells that arise during fetal pancreas development. HESCs were differentiated into fetal-like pancreatic cells in vitro using a 33-day, 7-stage protocol. Cultures were ~90-95% PDX1-positive by day (d) 11 and 70-75% NKX6.1-positive by d17. Polyhormonal cells were scattered at d17, but developed into islet-like clusters that expressed key transcription factors by d33. Human C-peptide and glucagon secretion were first detected at d17 and increased thereafter in parallel with INS and GCG transcript levels. HESC-derived cells were responsive to KCl and arginine, but not glucose in perifusion studies. Compared to adult human islets, hESC-derived cells expressed ~10-fold higher levels of glucose transporter 1 (GLUT1) mRNA, but similar levels of glucokinase (GCK). In situ hybridization confirmed the presence of GLUT1 transcript within endocrine cells. However, GLUT1 protein was excluded from this population and was instead observed predominantly in non-endocrine cells, whereas GCK was co-expressed in insulin-positive cells. In rubidium efflux assays, hESC-derived cells displayed mild potassium channel activity, but no responsiveness to glucose, metabolic inhibitors or glibenclamide. Western blotting experiments revealed that the higher molecular weight SUR1 band was absent in hESC-derived cells, suggesting a lack of functional KATP channels at the cell surface. In addition, KATP channel subunit transcript levels were not at a 1:1 ratio, as would be expected (SUR1 levels were ~5-fold lower than KIR6.2). Various ratios of SUR1:KIR6.2 plasmids were transfected into COSM6 cells and rubidium efflux was found to be particularly sensitive to a reduction in SUR1. These data suggest that an impaired ratio of SUR1:KIR6.2 may contribute to the observed KATP channel defects in hESC-derived islet endocrine cells, and along with lack of GLUT1, may explain the absence of glucose-stimulated insulin secretion.

Logistic Equations Effectively Model Mammalian Cell Batch and Fed‐Batch Kinetics by Logically Constraining the Fit

A four‐parameter logistic equation was used to fit batch and fed‐batch time profiles of viable cell density in order to estimate net growth rates from the inoculation through the cell death phase. Reduced three‐parameter forms were used for nutrient uptake and metabolite/product formation rate calculations. These logistic equations constrained the fits to expected general concentration trends, either increasing followed by decreasing (four‐parameter) or monotonic (three‐parameter). The applicability of this approach was first verified for Chinese hamster ovary (CHO) cells cultivated in 15‐L batch bioreactors. Cell density, metabolite, and nutrient concentrations were monitored over time and used to estimate the logistic parameters by nonlinear least squares. The logistic models fit the experimental data well, supporting the validity of this approach. Further evidence to this effect was obtained by applying the technique to three previously published batch studies for baby hamster kidney (BHK) and hybridoma cells in bioreactors ranging from 100 mL to 300 L. In 27 of the 30 batch data sets examined, the logistic models provided a statistically superior description of the experimental data than polynomial fitting. Two fed‐batch experiments with hybridoma and CHO cells in benchtop bioreactors were also examined, and the logistic fits provided good representations of the experimental data in all 25 data sets. From a computational standpoint, this approach was simpler than classical approaches involving Monod‐type kinetics. Since the logistic equations were analytically differentiable, specific rates could be readily estimated. Overall, the advantages of the logistic modeling approach should make it an attractive option for effectively estimating specific rates from batch and fed‐batch cultures.