Single-cell analysis for bioprocessing

Cell population heterogeneity has attracted great interest for understanding the individual cellular performances in their response to external stimuli and in the production of targeted products. Physical characterization of single cells and analysis of dynamic gene expression, synthesized proteins, and cellular metabolites from one single cell are reviewed. Advanced techniques have been developed to achieve high-throughput and ultrahigh resolution or sensitivity. Single cell capture methods are discussed as well. How to make use of cellular heterogeneities for maximizing cellular productivity is still in the infant stage, and control strategies will be formulated after the causes for heterogeneity have been elucidated.

Single-cell functional analysis of parathyroid adenomas reveals distinct classes of calcium sensing behaviour in primary hyperparathyroidism.

Primary hyperparathyroidism (PHPT) is a common endocrine neoplastic disorder caused by a failure of calcium sensing secondary to tumour development in one or more of the parathyroid glands. Parathyroid adenomas are comprised of distinct cellular subpopulations of variable clonal status that exhibit differing degrees of calcium responsiveness. To gain a clearer understanding of the relationship among cellular identity, tumour composition and clinical biochemistry in PHPT, we developed a novel single cell platform for quantitative evaluation of calcium sensing behaviour in freshly resected human parathyroid tumour cells. Live-cell intracellular calcium flux was visualized through Fluo-4-AM epifluorescence, followed by in situ immunofluorescence detection of the calcium sensing receptor (CASR), a central component in the extracellular calcium signalling pathway. The reactivity of individual parathyroid tumour cells to extracellular calcium stimulus was highly variable, with discrete kinetic response patterns observed both between and among parathyroid tumour samples. CASR abundance was not an obligate determinant of calcium responsiveness. Calcium EC50 values from a series of parathyroid adenomas revealed that the tumours segregated into two distinct categories. One group manifested a mean EC50 of 2.40 mM (95% CI: 2.37-2.41), closely aligned to the established normal range. The second group was less responsive to calcium stimulus, with a mean EC50 of 3.61 mM (95% CI: 3.45-3.95). This binary distribution indicates the existence of a previously unappreciated biochemical sub-classification of PHPT tumours, possibly reflecting distinct etiological mechanisms. Recognition of quantitative differences in calcium sensing could have important implications for the clinical management of PHPT.

Single-Cell Transcriptome Analysis of Olfactory Sensory Neurons

Olfactory sensory neurons (OSNs), which detect a myriad of odorants, are known to express one allele of one olfactory receptor (OR) gene (Olfr) from the largest gene family in the mammalian genome. The OSNs expressing the same OR project their axons to the main olfactory bulb where they converge to form glomeruli. This “One neuron-one receptor rule” makes the olfactory epithelium (OE), which consists of a vast number of OSNs expressing unique ORs, one of the most heterogeneous cell populations. However, the mechanism of how the single OR allele is chosen remains unclear along with the question of whether one OSN only expresses a single OR gene, a hypothesis that has not been rigorously verified while we performed the experiments. Moreover, failure of axonal targeting to single glomerulus was observed in MeCP2 deficient OSNs where delayed development was proposed as an explanation for the phenotype. How Mecp2 mutation caused this aberrant targeting is not entirely understood.

In this dissertation, we explored the transcriptomes of single and mature OSNs by single-cell RNA-Seq to reveal their heterogeneity and further studied the OR gene expression from these isolated OSNs. The singularity of sequenced OSNs was ensured by the observation of monoallelic expression of X-linked genes from the hybrid samples from crosses between mice of different strains where strain-specific polymorphisms could be used to track the allelic origins of SNP-containing reads. The clustering of expression profiles from triplicates that originated from the same cell assured that the transcriptomic identities of OSNs were maintained through the experimental process. The average gene expression profiles of sequenced OSNs correlated well to the conventional transcriptome data of FACS-sorted Omp-positive cells, and the top-ranked expression of OR was conceded in the single-OSN transcriptomes. While exploring cellular diversity, in addition to OR genes, we revealed nearly 200 differentially expressed genes among the sequenced OSNs in this study. Among the 36 sequenced OSNs, eight cells (22.2%) showed multiple OR gene expression and the presences of additional ORs were not restricted to the neighbor loci that shared the transcriptional effect of the primary OR expression, suggesting that the “One neuron-one receptor rule” might not be strictly true at the transcription level. All of the inferable ORs, including additional co-expressed ORs, were shown to be monoallelic. Our sequencing of 21 Mecp2308 mutant OSNs, of which 62% expressed more than one OR genes, and the expression levels of the additional ORs were significantly higher than those in the wild-type, suggested that MeCP2 plays a role in the regulation of singular OR gene expression. Dual label in situ hybridization along with the sequence data revealed that dorsal and ventral ORs were co-expressed in the same Mecp2 mutant OSN, further implying that MeCP2 might be involved in regulation of OR territories in the OE. Our results suggested a new role of MeCP2 in OR gene choice and ratified that this multiple-OR expression caused by Mecp2 mutation did not accompany delayed OSN development that has been observed in the previous studies on the Mecp2 mutants.

Microfluidics-based single-cell functional proteomics microchip for portraying protein signal transduction networks within the framework of physicochemical principles, with applications in fundamental and translational cancer research

Single-cell functional proteomics assays can connect genomic information to biological function through quantitative and multiplex protein measurements. Tools for single-cell proteomics have developed rapidly over the past 5 years and are providing unique opportunities. This thesis describes an emerging microfluidics-based toolkit for single cell functional proteomics, focusing on the development of the single cell barcode chips (SCBCs) with applications in fundamental and translational cancer research.

The microchip designed to simultaneously quantify a panel of secreted, cytoplasmic and membrane proteins from single cells will be discussed at the beginning, which is the prototype for subsequent proteomic microchips with more sophisticated design in preclinical cancer research or clinical applications. The SCBCs are a highly versatile and information rich tool for single-cell functional proteomics. They are based upon isolating individual cells, or defined number of cells, within microchambers, each of which is equipped with a large antibody microarray (the barcode), with between a few hundred to ten thousand microchambers included within a single microchip. Functional proteomics assays at single-cell resolution yield unique pieces of information that significantly shape the way of thinking on cancer research. An in-depth discussion about analysis and interpretation of the unique information such as functional protein fluctuations and protein-protein correlative interactions will follow.

The SCBC is a powerful tool to resolve the functional heterogeneity of cancer cells. It has the capacity to extract a comprehensive picture of the signal transduction network from single tumor cells and thus provides insight into the effect of targeted therapies on protein signaling networks. We will demonstrate this point through applying the SCBCs to investigate three isogenic cell lines of glioblastoma multiforme (GBM).

The cancer cell population is highly heterogeneous with high-amplitude fluctuation at the single cell level, which in turn grants the robustness of the entire population. The concept that a stable population existing in the presence of random fluctuations is reminiscent of many physical systems that are successfully understood using statistical physics. Thus, tools derived from that field can probably be applied to using fluctuations to determine the nature of signaling networks. In the second part of the thesis, we will focus on such a case to use thermodynamics-motivated principles to understand cancer cell hypoxia, where single cell proteomics assays coupled with a quantitative version of Le Chatelier's principle derived from statistical mechanics yield detailed and surprising predictions, which were found to be correct in both cell line and primary tumor model.

The third part of the thesis demonstrates the application of this technology in the preclinical cancer research to study the GBM cancer cell resistance to molecular targeted therapy. Physical approaches to anticipate therapy resistance and to identify effective therapy combinations will be discussed in detail. Our approach is based upon elucidating the signaling coordination within the phosphoprotein signaling pathways that are hyperactivated in human GBMs, and interrogating how that coordination responds to the perturbation of targeted inhibitor. Strongly coupled protein-protein interactions constitute most signaling cascades. A physical analogy of such a system is the strongly coupled atom-atom interactions in a crystal lattice. Similar to decomposing the atomic interactions into a series of independent normal vibrational modes, a simplified picture of signaling network coordination can also be achieved by diagonalizing protein-protein correlation or covariance matrices to decompose the pairwise correlative interactions into a set of distinct linear combinations of signaling proteins (i.e. independent signaling modes). By doing so, two independent signaling modes – one associated with mTOR signaling and a second associated with ERK/Src signaling have been resolved, which in turn allow us to anticipate resistance, and to design combination therapies that are effective, as well as identify those therapies and therapy combinations that will be ineffective. We validated our predictions in mouse tumor models and all predictions were borne out.

In the last part, some preliminary results about the clinical translation of single-cell proteomics chips will be presented. The successful demonstration of our work on human-derived xenografts provides the rationale to extend our current work into the clinic. It will enable us to interrogate GBM tumor samples in a way that could potentially yield a straightforward, rapid interpretation so that we can give therapeutic guidance to the attending physicians within a clinical relevant time scale. The technical challenges of the clinical translation will be presented and our solutions to address the challenges will be discussed as well. A clinical case study will then follow, where some preliminary data collected from a pediatric GBM patient bearing an EGFR amplified tumor will be presented to demonstrate the general protocol and the workflow of the proposed clinical studies.

Metastability and dynamics of stem cells: from direct observations to inference at the single cell level

Organismal development, homeostasis, and pathology are rooted in inherently probabilistic events. From gene expression to cellular differentiation, rates and likelihoods shape the form and function of biology. Processes ranging from growth to cancer homeostasis to reprogramming of stem cells all require transitions between distinct phenotypic states, and these occur at defined rates. Therefore, measuring the fidelity and dynamics with which such transitions occur is central to understanding natural biological phenomena and is critical for therapeutic interventions.

While these processes may produce robust population-level behaviors, decisions are made by individual cells. In certain circumstances, these minuscule computing units effectively roll dice to determine their fate. And while the 'omics' era has provided vast amounts of data on what these populations are doing en masse, the behaviors of the underlying units of these processes get washed out in averages.

Therefore, in order to understand the behavior of a sample of cells, it is critical to reveal how its underlying components, or mixture of cells in distinct states, each contribute to the overall phenotype. As such, we must first define what states exist in the population, determine what controls the stability of these states, and measure in high dimensionality the dynamics with which these cells transition between states.

To address a specific example of this general problem, we investigate the heterogeneity and dynamics of mouse embryonic stem cells (mESCs). While a number of reports have identified particular genes in ES cells that switch between 'high' and 'low' metastable expression states in culture, it remains unclear how levels of many of these regulators combine to form states in transcriptional space. Using a method called single molecule mRNA fluorescent in situ hybridization (smFISH), we quantitatively measure and fit distributions of core pluripotency regulators in single cells, identifying a wide range of variabilities between genes, but each explained by a simple model of bursty transcription. From this data, we also observed that strongly bimodal genes appear to be co-expressed, effectively limiting the occupancy of transcriptional space to two primary states across genes studied here. However, these states also appear punctuated by the conditional expression of the most highly variable genes, potentially defining smaller substates of pluripotency.

Having defined the transcriptional states, we next asked what might control their stability or persistence. Surprisingly, we found that DNA methylation, a mark normally associated with irreversible developmental progression, was itself differentially regulated between these two primary states. Furthermore, both acute or chronic inhibition of DNA methyltransferase activity led to reduced heterogeneity among the population, suggesting that metastability can be modulated by this strong epigenetic mark.

Finally, because understanding the dynamics of state transitions is fundamental to a variety of biological problems, we sought to develop a high-throughput method for the identification of cellular trajectories without the need for cell-line engineering. We achieved this by combining cell-lineage information gathered from time-lapse microscopy with endpoint smFISH for measurements of final expression states. Applying a simple mathematical framework to these lineage-tree associated expression states enables the inference of dynamic transitions. We apply our novel approach in order to infer temporal sequences of events, quantitative switching rates, and network topology among a set of ESC states.

Taken together, we identify distinct expression states in ES cells, gain fundamental insight into how a strong epigenetic modifier enforces the stability of these states, and develop and apply a new method for the identification of cellular trajectories using scalable in situ readouts of cellular state.

Purification, biochemical characterization and localization at single cell resolution of Matrix Gla Protein (MGP) and Bone Gla Protein (BGP) in the teleost fish Argyrosomus regius

Tese de dout. em Química, Faculdade de Ciências do Mar e do Ambiente, Univ. do Algarve, 2002

Phenotypic and molecular investigation of circulating tumor cells (CTCs) isolated from breast cancer patients at single cell resolution

Despite the paramount advances in cancer research, breast cancer (BC) still ranks one of the leading causes of cancer-related death worldwide. Thanks to the screening campaign started in developed countries, BC is often diagnosed at early stages (non-metastatic BC, nmBC), but disease relapse occurrence even after decades and at distant sites is not an uncommon phenomenon. Conversely, metastatic BC (mBC) is considered an incurable disease. The major perpetrators of tumor spread to secondary organs are circulating tumor cells (CTCs), a rare population of cells detectable in the peripheral blood of oncologic patients. In this study, CTCs from patients diagnosed with luminal nmBC and mBC (hormone receptor positive, Human Epidermal Growth Factor Receptor 2 (HER2) negative) were characterized at both phenotypic and molecular levels. To better understand the molecular mechanisms underlying their biology and their metastatic potential, next-generation sequencing (NGS) analyses were performed at single-cell resolution to assess copy number aberrations (CNAs), single nucleotide variants (SNVs) and gene expression profiling. The findings of this study arise hints in CTC detection, and pave the way to new application in CTC research.

Increased natural killer activity does not prevent progression of experimental Kala-azar

Kala-azar is the visceral form of leishmaniasis and it is caused by intracellular parasites from the complex Leishmania donovani. Golden hamster (Mesocricetus auratus) infected with Leishmania donovani develop a disease very similar to human Kala-azar. There is conspicuous hipergammaglobulinaemia and their T cells do not respond to stimulation with parasite antigens. We used this experimental model to evaluate the natural killer (NK) activity during the initial phase of the disease. Outbred hamsters infected by intravenous route with 5.106 amastigotes of L. donovani 1S showed a concurrent increase in the spleen weight and in the spleen cell number. Using the single cell assay we detected a significant increase in the percentage of NK effector cells on the 4th day of infection. Imprints from spleen and liver showed at days 14 and 28 a significant increase in the parasite burden . These results show that the increased NK activity in the beginning of the infection was not able to restrain the progression of the disease in this experimental model.

Studies of natural-killer-cells in pregnancy-induced hypertension

The number and activity of natural killer (NK) cells were studied in 20 patients with pregnancy-induced hypertension (PIH), 15 uncomplicated pregnant women and 16 healthy non-pregnant women, All the pregnant women were primigravidae and were evaluated during the third trimester of gestation. Peripheral blood NK cells were detected with monoclonal antibodies by indirect immunofluorescence and cytotoxic activity was measured using a single-cell assay against K562 target cells. Hypertensive pregnant women had an increased number of circulating NK cells associated with a significant decrease of NK activity, the cytotoxic activity was significantly lower in normal pregnant and PIH women when compared with non-pregnant controls. The onset of immature NK cells in peripheral blood and the impairment of their cytotoxic activity in PIH patients may be associated with hormones and immunosuppressive substances produced by tissues occurring at the maternal-fetal interface.

Immunological alterations in patients with primary tumors in central nervous system

Natural killer (NK) cells play an important role in immune surveillance against tumors. The present work aimed to study the cytotoxic activity of NK cells and T cell subsets in peripheral blood of 13 patients with primary tumors in central nervous system (CNS). As controls 29 healthy subjects with the age range equivalent to the patients were studied. The methods employed were: a) determination of cytotoxic activity of NK cells towards K562 target cells, evaluated by single cell-assay; b) enumeration of CD3+ lymphocytes and their CD4+ and CD8+ subsets defined by monoclonal antibodies; c) the identification of tumors were done by histologic and immunochemistry studies. The results indicated that adults and children with tumor in CNS display reduced percentage of total T cells, helper/inducer subset and low helper/suppressor ratio. The cytotoxic activity of NK cells was decreased in patients with CNS tumors due mainly to a decrease in the proportion of target-binding lymphocytes. These results suggest that cytotoxic activity of NK cells may be affected by the immunoregulatory disturbances observed in patients with primary tumors in CNS.

Increased natural killer activity does not prevent progression of experimental kala-azar

Kala-azar is the visceral form of leishmaniasis and it is caused by intracellular parasites from the complex Leishmania donovani. Golden hamster (Mesocricetus auratus) infected with Leishmania donovani develop a disease very similar to human Kala-azar. There is conspicuous hipergammaglobulinaemia and their T cells do not respond to stimulation with parasite antigens. We used this experimental model to evaluate the natural killer (NK) activity during the initial phase of the disease. Outbred hamsters infected by intravenous route with 5.106 amastigotes of L. donovani 1S showed a concurrent increase in the spleen weight and in the spleen cell number. Using the single cell assay we detected a significant increase in the percentage of NK effector cells on the 4th day of infection. Imprints from spleen and liver showed at days 14 and 28 a significant increase in the parasite burden. These results show that the increased NK activity in the beginning of the infection was not able to restrain the progression of the disease in this experimental model.

Studies of natural killer cells in patients with paracoccidioidomycosis

The number and activity of natural killer (NK) cells were studied in 34 untreated patients with paracoccidioidomycosis, 20 with the chronic form of the disease and 14 with the acute form. NK cells were detected with monoclonal antibody Leu-11c and the cytotoxic activity was measured using a single cell assay against K562 target cells. Both groups of patients had an increased number of circulating NK cells, their cytotoxic activity being significantly lower than in the healthy controls. These findings may be of importance in the immunological disturbances associated with paracoccidioidomycosis since NK cells exert important immune effector functions and may play a role in resistance against Paracoccidioides brasiliensis.

DNA damage by ochratoxin A in rat kidney assessed by the alkaline comet assay

There are few studies of ochratoxin A (OTA) genotoxicity in experimental animals and the results obtained with cell cultures are inconsistent, although the carcinogenic potential of OTA for the kidney of experimental animals has been well established. We studied the genotoxic potential of OTA in the kidney of adult female Wistar rats (5 in each group) treated intraperitoneally with OTA (0.5 mg kg body weight-1 day-1 for 7, 14, and 21 days) measuring DNA mobility on agarose gel stained with ethidium-bromide using standard alkaline single-cell gel electrophoresis (comet assay). Negative control animals were treated with solvent (Tris buffer, 1.0 mg/kg) and positive control animals were treated with methyl methanesulfonate (40 mg/kg) according to the same schedule. OTA concentrations in plasma and kidney homogenates in 7-, 14-, and 21-day treated animals were 4.86 ± 0.53, 7.52 ± 3.32, 7.85 ± 2.24 µg/mL, and 0.87 ± 0.09, 0.99 ± 0.06, 1.09 ± 0.15 µg/g, respectively. In all OTA-treated groups, the tail length, tail intensity, and tail moment in kidney tissue were significantly higher than in controls (P < 0.05). The tail length and tail moment were higher after 14 days than after 7 days of treatment (P < 0.05), and still higher after 21 days (P < 0.05). The highest tail intensity was observed in animals treated for 21 days, and it differed significantly from animals treated for 7 and 14 days (P < 0.05). OTA concentrations in plasma and kidney tissue increased steadily and OTA concentration in kidney tissue strongly correlated with tail intensity and tail moment values. These results confirm the genotoxic potential of OTA, and show that the severity of DNA lesions in kidney correlates with OTA concentration.

Transcriptional regulation of effector CD8+ T cell differentiation and molecular dysfunction during HIV-1 infection

Les cellules T CD8+ jouent un rôle primordial dans le contrôle des infections virales en limitant la dissémination des cellules infectées. Lors de l’infection chronique par le virus HIV, les cellules T CD8+ HIV-spécifiques ne se différencient pas en cellules effectrices fonctionnelles capables de tuer les cellules infectées par le virus ; ces cellules ne sont plus capables de proliférer ou de produire l’ IL-2. Ces cellules expriment PD-1 et l’engagement de PD-1, par son ligand, aboutit a plusieurs de ces déficits fonctionnels des cellules T . Le rôle de PD-1 dans la régulation d'évènements transcriptionnels contrôlant la différentiation et l'obtention des fonction effectrices des cellules T CD8+ reste à démontrer. Id2 joue un rôle central dans la différenciation des cellules T CD8+ effectrices. Nous avons émis l’hypothèse que le défaut de maturation observé chez les cellules T CD8+ PD-1 high HIV-spécifiques (CD8+PD-1hi) au cours de l’infection chronique par le virus HIV pouvait être lié à la diminution d’expression du régulateur Id2. Nous avons ainsi démontré que l'engagement de PD-1 contribuait à une diminution d'expression de Id2 et de ses cibles transcriptionnelles. La surexpression de Id2 de ces cellules a permis de restaurer l'expression de marqueurs tels que Granzyme B et Bcl-2 et diminuir l’expression du marqueur de maturation de CD27. La famille des cytokines à chaine gamma joue un rôle clef dans la survie et l’homéostasie des cellules T. Dans ce travail, nous avons démontré que l’IL-15 était unique grâce à ses capacités de stimulation de l’expression d’Id2 et ses propriétés favorisant la survie ainsi que la différenciation des cellules T CD8+ effectrices. l’IL-15 induit la prolifération de toutes les populations de cellules T mémoires provenant de donneurs sains. L’addition de cette cytokine aux sous-populations cellulaires Ttm et Tem a permis leur différenciation en cellules effectrices capables de produire Granzyme B alors que la stimulation par l’IL-15 des cellules Tcm ne favorise pas leur différenciation. Un test de cytotoxicitié par cytométrie en flux nous a permis de confirmer que la stimulation de cellules T CD8+ HIV spécifiques par l’IL-15 favorisait l’expression de Id2 et restaurait les fonctions cytotoxiques des cellules T CD8+ HIV spécifiques. En conclusion, nous avons pour la première fois dans cette thèse défini les mécanismes moléculaires impliqués dans la modulation de l’expression du régulateur transcriptionnel Id2 par l’IL-15. Nous avons également révélé comment l’engagement de PD-1 conduisait a une altération de l’expression et de la fonction d’Id2 et favorisait la diminution des fonctions effectrices des cellules T CD8-HIV spécifiques. Une perspective de traitement avec des agents tels que l’IL-15 ou le bloquage de PD-1, en combinaison avec les traitements conventionnels, pourrait contribuer à une meilleure stimulation des réponses immunes favorisant ainsi la réactivation des cellules T CD8+ et permettant la destruction de cellules T CD4+ infectées de manière latente.