Population Transcriptomics Research Articles

Population transcriptomics with single‐cell resolution: A new field made possible by microfluidics

Tissues contain complex populations of cells. Like countries, which are comprised of mixed populations of people, tissues are not homogeneous. Gene expression studies that analyze entire populations of cells from tissues as a mixture are blind to this diversity. Thus, critical information is lost when studying samples rich in specialized but diverse cells such as tumors, iPS colonies, or brain tissue. High throughput methods are needed to address, model and understand the constitutive and stochastic differences between individual cells. Here, we describe microfluidics technologies that utilize a combination of molecular biology and miniaturized labs on chips to study gene expression at the single cell level. We discuss how the characterization of the transcriptome of each cell in a sample will open a new field in gene expression analysis, population transcriptomics, that will change the academic and biomedical analysis of complex samples by defining them as quantified populations of single cells.

Gene transcription reflects poor health status of resident European eel chronically exposed to environmental pollutants

Understanding the effects of chronic exposure to pollutants on the genome and transcriptome of diadromous fish populations is crucial for their resilience under combined anthropogenic and environmental selective pressures. The catadromous European eel (Anguilla anguilla L.) has suffered a dramatic decline in recruitment for three decades, necessitating a thorough assessment of the transcriptional effects of environmental pollutants on resident and migrating eels in natural systems. We investigated the relationship between muscular bioaccumulation levels of metals (Hg, Cd, Pb, Cu, Zn, Ni, Cr, As and Se), PCBs and organochlorine pesticides (DDTs), the health status (condition factor and lipid reserves) and the associated transcriptional response in liver and gill tissues for genes involved in metal detoxification (metallothionein, MT) and oxidative metabolism (cytochrome P4501A, CYP1A) of xenobiotic compounds. In total 84 resident eels originating from three Belgian river basins (Scheldt, Meuse and Yzer) were analyzed along with five unpolluted aquaculture samples as control group. There was a large spatial variation in individual contaminant intensity and profile, while tissue pollution levels were strongly and negatively associated with condition indices, suggesting an important impact of pollution on the health of sub-adult resident eels. Gene transcription patterns revealed a complex response mechanism to a cocktail of pollutants, with a high variation at low pollution levels, but strongly down-regulated hepatic and gill gene transcription in highly polluted eels. Resident eels clearly experience a high pollution burden and seem to show a dysfunctional gene transcription regulation of detoxification genes at higher pollutant levels, correlated with low energy reserves and condition. To fully understand the evolutionary implications of pollutants on eel reproductive fitness, analyses of mature migrating eels and the characterization of their transcriptome-wide gene transcription response would be appropriate to unveil the complex responses associated with multiple interacting stressors and the long-term consequences at the entire species level. In the meanwhile, jointly monitoring environmental and tissue pollution levels at a European scale should be initiated, while preserving high quality habitats to increase the recovery chance of European eel in the future.

The Transcriptional Profile of Mesenchymal Stem Cell Populations in Primary Osteoporosis Is Distinct and Shows Overexpression of Osteogenic Inhibitors

Primary osteoporosis is an age-related disease characterized by an imbalance in bone homeostasis. While the resorptive aspect of the disease has been studied intensely, less is known about the anabolic part of the syndrome or presumptive deficiencies in bone regeneration. Multipotent mesenchymal stem cells (MSC) are the primary source of osteogenic regeneration. In the present study we aimed to unravel whether MSC biology is directly involved in the pathophysiology of the disease and therefore performed microarray analyses of hMSC of elderly patients (79–94 years old) suffering from osteoporosis (hMSC-OP). In comparison to age-matched controls we detected profound changes in the transcriptome in hMSC-OP, e.g. enhanced mRNA expression of known osteoporosis-associated genes (LRP5, RUNX2, COL1A1) and of genes involved in osteoclastogenesis (CSF1, PTH1R), but most notably of genes coding for inhibitors of WNT and BMP signaling, such as Sclerostin and MAB21L2. These candidate genes indicate intrinsic deficiencies in self-renewal and differentiation potential in osteoporotic stem cells. We also compared both hMSC-OP and non-osteoporotic hMSC-old of elderly donors to hMSC of ∼30 years younger donors and found that the transcriptional changes acquired between the sixth and the ninth decade of life differed widely between osteoporotic and non-osteoporotic stem cells. In addition, we compared the osteoporotic transcriptome to long term-cultivated, senescent hMSC and detected some signs for pre-senescence in hMSC-OP.Our results suggest that in primary osteoporosis the transcriptomes of hMSC populations show distinct signatures and little overlap with non-osteoporotic aging, although we detected some hints for senescence-associated changes. While there are remarkable inter-individual variations as expected for polygenetic diseases, we could identify many susceptibility genes for osteoporosis known from genetic studies. We also found new candidates, e.g. MAB21L2, a novel repressor of BMP-induced transcription. Such transcriptional changes may reflect epigenetic changes, which are part of a specific osteoporosis-associated aging process.

Activation of two distinct Sox9-EGFP-expressing intestinal stem cell populations during crypt regeneration after irradiation

Recent identification of intestinal epithelial stem cell (ISC) markers and development of ISC reporter mice permit visualization and isolation of regenerating ISCs after radiation to define their functional and molecular phenotypes. Previous studies in uninjured intestine of Sox9-EGFP reporter mice demonstrate that ISCs express low levels of Sox9-EGFP (Sox9-EGFP Low), whereas enteroendocrine cells (EEC) express high levels of Sox9-EGFP (Sox9-EGFP High). We hypothesized that Sox9-EGFP Low ISCs would expand after radiation, exhibit enhanced proliferative capacities, and adopt a distinct gene expression profile associated with rapid proliferation. Sox9-EGFP mice were given 14 Gy abdominal radiation and studied between days 3 and 9 postradiation. Radiation-induced changes in number, growth, and transcriptome of the different Sox9-EGFP cell populations were determined by histology, flow cytometry, in vitro culture assays, and microarray. Microarray confirmed that nonirradiated Sox9-EGFP Low cells are enriched for Lgr5 mRNA and mRNAs enriched in Lgr5-ISCs and identified additional putative ISC markers. Sox9-EGFP High cells were enriched for EEC markers, as well as Bmi1 and Hopx, which are putative markers of quiescent ISCs. Irradiation caused complete crypt loss, followed by expansion and hyperproliferation of Sox9-EGFP Low cells. From nonirradiated intestine, only Sox9-EGFP Low cells exhibited ISC characteristics of forming organoids in culture, whereas during regeneration both Sox9-EGFP Low and High cells formed organoids. Microarray demonstrated that regenerating Sox9-EGFP High cells exhibited transcriptomic changes linked to p53-signaling and ISC-like functions including DNA repair and reduced oxidative metabolism. These findings support a model in which Sox9-EGFP Low cells represent active ISCs, Sox9-EGFP High cells contain radiation-activatable cells with ISC characteristics, and both participate in crypt regeneration.

Expression QTL mapping in regulatory and helper T cells from the BXD family of strains reveals novel cell-specific genes, gene-gene interactions and candidate genes for auto-immune disease

BackgroundRegulatory T cells (Tregs) play an essential role in the control of the immune response. Treg cells represent important targets for therapeutic interventions of the immune system. Therefore, it will be very important to understand in more detail which genes are specifically activated in Treg cells versus T helper (Th) cells, and which gene regulatory circuits may be involved in specifying and maintaining Treg cell homeostasis.ResultsWe isolated Treg and Th cells from a genetically diverse family of 31 BXD type recombinant inbred strains and the fully inbred parental strains of this family--C57BL/6J and DBA/2J. Subsequently genome-wide gene expression studies were performed from the isolated Treg and Th cells. A comparative analysis of the transcriptomes of these cell populations allowed us to identify many novel differentially expressed genes. Analysis of cis- and trans-expression Quantitative Trait Loci (eQTLs) highlighted common and unique regulatory mechanisms that are active in the two cell types. Trans-eQTL regions were found for the Treg functional genes Nrp1, Stat3 and Ikzf4. Analyses of the respective QTL intervals suggested several candidate genes that may be involved in regulating these genes in Treg cells. Similarly, possible candidate genes were found which may regulate the expression of F2rl1, Ctla4, Klrb1f. In addition, we identified a focused group of candidate genes that may be important for the maintenance of self-tolerance and the prevention of allergy.ConclusionsVariation of expression across the strains allowed us to find many novel gene-interaction networks in both T cell subsets. In addition, these two data sets enabled us to identify many differentially expressed genes and to nominate candidate genes that may have important functions for the maintenance of self-tolerance and the prevention of allergy.

Interindividual Variability in Gene Expression Profiles in Human Hepatocytes and Comparison with HepaRG Cells

Interindividual variations in functions other than drug metabolism activity, remain poorly elucidated in human liver. In the present study, the whole transcriptome of several human hepatocyte populations and the differentiated human HepaRG cell line have been analyzed and compared, using oligonucleotide pangenomic microarrays. We show that, although the variation in the percentages of expressed genes did not exceed 14% among the primary human hepatocyte populations, huge interindividual differences in the transcript levels of many genes were observed. These genes were related to various functions; in addition to drug metabolism, they mainly concerned carbohydrate, amino acid, and lipid metabolism. HepaRG cells expressed from 81 to 92% of the genes active in human hepatocytes and, in addition, a specific gene subset mainly related to their transformed status, some chromosomal abnormalities, and the presence of primitive biliary epithelial cells. Of interest, a relationship was evidenced between abnormal basal expression levels of some target genes and their corresponding previously reported fold changes in one of four human hepatocyte populations treated with the hepatotoxic drug troglitazone and not with other nonhepatotoxic peroxisome proliferator-activated receptor agonists (PLoS One 6:e18816, 2011). Taken together, our results support the view that HepaRG cells express most of the genes active in primary human hepatocytes and show that expression of most human hepatic genes can quantitatively greatly vary among individuals, thereby contributing to explain the huge interindividual variability in susceptibility to drugs and other environmental factors.

Comprehensive transcriptome and immunophenotype analysis of renal and cardiac MSC-like populations supports strong congruence with bone marrow MSC despite maintenance of distinct identities

Cells resembling bone marrow mesenchymal stem cells (MSC) have been isolated from many organs but their functional relationships have not been thoroughly examined. Here we compared the immunophenotype, gene expression, multipotency and immunosuppressive potential of MSC-like colony-forming cells from adult murine bone marrow (bmMSC), kidney (kCFU-F) and heart (cCFU-F), cultured under uniform conditions. All populations showed classic MSC morphology and in vitro mesodermal multipotency. Of the two solid organ-specific CFU-F, only kCFU-F displayed suppression of T-cell alloreactivity in vitro, albeit to a lesser extent than bmMSC. Quantitative immunophenotyping using 81 phycoerythrin-conjugated CD antibodies demonstrated that all populations contained high percentages of cells expressing diagnostic MSC surface markers (Sca1, CD90.2, CD29, CD44), as well as others noted previously on murine MSC (CD24, CD49e, CD51, CD80, CD81, CD105). Illumina microarray expression profiling and bioinformatic analysis indicated a correlation of gene expression of 0.88–0.92 between pairwise comparisons. All populations expressed approximately 66% of genes in the pluripotency network (Plurinet), presumably reflecting their stem-like character. Furthermore, all populations expressed genes involved in immunomodulation, homing and tissue repair, suggesting these as conserved functions for MSC-like cells in solid organs. Despite this molecular congruence, strong biases in gene and protein expression and pathway activity were seen, suggesting organ-specific functions. Hence, tissue-derived MSC may also retain unique properties potentially rendering them more appropriate as cellular therapeutic agents for their organ of origin.

Population transcriptomics: insights from Drosophila simulans, Drosophila sechellia and their hybrids

Sequence differentiation has been widely studied between populations and species, whereas interest in expression divergence is relatively recent. Using microarrays, we compared four geographically distinct populations of Drosophila simulans and a population of Drosophila sechellia, and interspecific hybrids. We observed few differences between populations, suggesting a slight population structure in D. simulans. This structure was observed in direct population comparisons, as well as in interspecific comparisons (hybrids vs. parents, D. sechellia vs. D. simulans). Expression variance is higher in the French and Zimbabwean populations than in the populations from the ancestral range of D. simulans (Kenya and Seychelles). This suggests a large scale phenomenon of decanalization following the invasion of a new environment. Comparing D. simulans and D. sechellia, we revealed 304 consistently differentially expressed genes, with striking overrepresentation of genes of the cytochrome P450 family, which could be related to their role in detoxification as well as in hormone regulation. We also revealed differences in genes involved in Juvenile hormone and Dopamine differentiation. We finally observed very few differentially expressed genes between hybrids and parental populations, with an overrepresentation of X-linked genes.

Equivalence testing in microarray analysis: similarities in the transcriptome of human atherosclerotic and nonatherosclerotic macrophages

We focus on similarities in the transcriptome of human Kupffer cells and alveolar, splenic, and atherosclerotic plaque-residing macrophages. We hypothesized that these macrophages share a common expression signature. We performed microarray analysis on mRNA from these subsets (4 patients) and developed a novel statistical method to identify genes with significantly similar expression levels. Phenotypic and functional diversity between macrophage subpopulations reflects their plasticity to respond to microenvironmental signals. Apart from detecting differences in expression profiles, the comparison of the transcriptomes of different macrophage populations may also allow the definition of molecular similarities between these subsets. This new method calculates the maximum difference in gene expression level, based on the estimated confidence interval on that gene's expression variance. We listed the genes by equivalence ranking relative to expression level. FDR estimation was used to determine significance. We identified 500 genes with significantly equivalent expression levels in the macrophage subsets at 5.5% FDR using a confidence level of α = 0.05 for equivalence. Among these are the established macrophage marker CD68, IL1 receptor antagonist, and MHC-related CD1C. These 500 genes were submitted to IPA and GO clustering using DAVID. Additionally, hierarchical clustering of these genes in the Novartis human gene expression atlas revealed a subset of 200 genes specifically expressed in macrophages. Equivalently expressed genes, identified by this new method, may not only help to dissect common molecular mechanisms, but also to identify cell- or condition-specific sets of marker genes that can be used for drug targeting and molecular imaging.

Stage-Specific Role for Shh in Dopaminergic Differentiation of Human Embryonic Stem Cells Induced by Stromal Cells

Stromal cells have been used to induce dopaminergic differentiation of mouse, primate, and human embryonic stem cells (hESCs), but the mechanism that governs this induction is unknown. In this manuscript, we show that medium conditioned by the stromal cell line PA6 (PA6-CM) can induce dopaminergic differentiation in neural stem cells (NSCs) derived from hESCs but not directly from hESCs, indicating that soluble factors produced by PA6 cells act at the NSC stage to specify a dopaminergic fate. To identify such soluble factors, we analyzed the transcriptomes of PA6 cells, NSCs, and dopaminergic populations induced by PA6-CM from hESC-derived NSCs. We focused our analysis on growth factors expressed by PA6 and receptors expressed by NSCs, and generated a list of growth factors and receptors that are differentially expressed. Some of the growth factor/receptor pairs are categorized into the Shh, Wnt5A, TGFbeta, and IGF pathways. The expression of genes activated by these pathways in dopaminergic populations was analyzed to confirm that these signals were likely candidates for specifying dopaminergic fate. Results were verified for Shh by using perturbation agents such as cyclopamine to show that Shh is indeed one of the active agents in PA6-CM, and by showing that Shh and FGF8 can substitute for PA6-CM at the NSC induction stage. We conclude that PA6-CM can induce dopaminergic differentiation in hESCs in a stage-specific manner. Shh is likely an important soluble dopaminergic inducing factor secreted by stromal cells and acts after the neural fate determination.

Neural Differentiation of Embryonic Stem Cells In Vitro: A Road Map to Neurogenesis in the Embryo

BackgroundThe in vitro generation of neurons from embryonic stem (ES) cells is a promising approach to produce cells suitable for neural tissue repair and cell-based replacement therapies of the nervous system. Available methods to promote ES cell differentiation towards neural lineages attempt to replicate, in different ways, the multistep process of embryonic neural development. However, to achieve this aim in an efficient and reproducible way, a better knowledge of the cellular and molecular events that are involved in the process, from the initial specification of neuroepithelial progenitors to their terminal differentiation into neurons and glial cells, is required.Methodology/Principal FindingsIn this work, we characterize the main stages and transitions that occur when ES cells are driven into a neural fate, using an adherent monolayer culture system. We established improved conditions to routinely produce highly homogeneous cultures of neuroepithelial progenitors, which organize into neural tube-like rosettes when they acquire competence for neuronal production. Within rosettes, neuroepithelial progenitors display morphological and functional characteristics of their embryonic counterparts, namely, apico-basal polarity, active Notch signalling, and proper timing of production of neurons and glia. In order to characterize the global gene activity correlated with each particular stage of neural development, the full transcriptome of different cell populations that arise during the in vitro differentiation protocol was determined by microarray analysis. By using embryo-oriented criteria to cluster the differentially expressed genes, we define five gene expression signatures that correlate with successive stages in the path from ES cells to neurons. These include a gene signature for a primitive ectoderm-like stage that appears after ES cells enter differentiation, and three gene signatures for subsequent stages of neural progenitor development, from an early stage that follows neural induction to a final stage preceding terminal differentiation.Conclusions/SignificanceOverall, our work confirms and extends the cellular and molecular parallels between monolayer ES cell neural differentiation and embryonic neural development, revealing in addition novel aspects of the genetic network underlying the multistep process that leads from uncommitted cells to differentiated neurons.

Differentially Expressed and Novel Transcripts in Highly Purified Chronic Phase CML Stem Cells

The chronic phase of CML is sustained by rare BCR-ABL+ stem cells. These cells share many properties with normal pluripotent hematopoietic stem cells, but also differ in critical ways that alter their growth, drug responsiveness and genome stability. Understanding the molecular mechanisms underlying the biological differences between normal and CML stem cells is key to the development of more effective CML therapies. To obtain new insights into these mechanisms, we generated Long Serial Analysis of Gene Expression (SAGE) libraries from paired isolates of highly purified lin-CD34+CD45RA-CD36- CD71-CD7-CD38+ and lin-CD34+CD45RA-CD36-CD71-CD7-CD38- cells from 3 chronic phase CML patients (all with predominantly Ph+/BCR-ABL+ cells in both subsets) and from 3 control samples: a pool of 10 normal bone marrows (BMs), a single normal BM and a pool of G-CSF-mobilized blood cells from 9 donors. In vitro bioassays showed the CD34+CD38+ cells were enriched in CFCs (CML: 3–20% pure; normal: 4–19% pure) and the CD34+CD38- cells were enriched in LTC-ICs (CML: 0.2–26% pure; normal: 12–52% pure). Each of the 12 libraries was then sequenced to a depth of ~200,000 tags and tags from libraries prepared from like phenotypes were compared between genotypes using DiscoverySpace software and hierarchical clustering. 1687 (355 with clustering) and 1258 (316 with clustering) transcripts were thus identified as differentially expressed in the CML vs control CD34+CD38− and CD34+CD38+ subsets, respectively. 266 of these transcripts (11 with clustering) were differentially expressed in both subsets. The differential expression of 5 genes (GAS2, IGF2BP2, IL1R1, DUSP1 &SELL) was confirmed by real-time PCR analysis of lin-CD34+ cells isolated from an additional 5 normal BMs and 11 CMLs, and lin-CD34+CD38− cells from an additional 2 normal BMs and 2 CMLs (with dominant Ph+ cells). GAS2 and IL1R1 transcript levels were correlated with BCR-ABL transcript levels in both primitive subsets, and predicted differences in expression of IL1R1 and SELL were apparent within 3 days in CD34+ cord blood cells transduced with a lenti-BCR-ABL-IRES-GFP vs a control lenti-GFP vector (n=3). These findings support a direct role of BCR-ABL in perturbing the expression of these 3 genes. Further comparison of the meta CD34+CD38− and CD34+CD38+ CML cell libraries with most publicly accessible SAGE data revealed 69 novel tags in the CD34+ CML cells that correspond to unique but conserved genomic sequences. Nine of these were recovered by 5′- and 3′- RACE applied to cDNAs pooled from several human leukemic cell lines. These results illustrate the power of SAGE to reveal key components of the transcriptomes of rare human CML stem cell populations including transcripts of genes not previously known to exist. Continuing investigation of their biological roles in primary CML cells and primitive BCR-ABL-transduced human cells offer important strategies for delineating their potential as therapeutic targets.

Population transcriptomics of life‐history variation in the genusSalmo

In this study, we demonstrate the power of applying complementary DNA (cDNA) microarray technology to identifying candidate loci that exhibit subtle differences in expression levels associated with a complex trait in natural populations of a nonmodel organism. Using a highly replicated experimental design involving 180 cDNA microarray experiments, we measured gene-expression levels from 1098 transcript probes in 90 individuals originating from six brown trout (Salmo trutta) and one Atlantic salmon (Salmo salar) population, which follow either a migratory or a sedentary life history. We identified several candidate genes associated with preparatory adaptations to different life histories in salmonids, including genes encoding for transaldolase 1, constitutive heat-shock protein HSC70-1 and endozepine. Some of these genes clustered into functional groups, providing insight into the physiological pathways potentially involved in the expression of life-history related phenotypic differences. Such differences included the down-regulation of genes involved in the respiratory system of future migratory individuals. In addition, we used linear discriminant analysis to identify a set of 12 genes that correctly classified immature individuals as migratory or sedentary with high accuracy. Using the expression levels of these 12 genes, 17 out of 18 individuals used for cross-validation were correctly assigned to their respective life-history phenotype. Finally, we found various candidate genes associated with physiological changes that are likely to be involved in preadaptations to seawater in anadromous populations of the genus Salmo, one of which was identified to encode for nucleophosmin 1. Our findings thus provide new molecular insights into salmonid life-history variation, opening new perspectives in the study of this complex trait.

Comparative Transcriptome Analysis of Different Subsets of CD34+ Normal and Chronic Myeloid Leukemia Cells Identifies Novel Perturbations in the CML Stem Cell Population.

Chronic myeloid leukemia (CML) arises from a hematopoietic stem cell that acquires a BCR-ABL fusion gene. During the chronic phase of the disease, this cell produces an expanding multi-lineage clone that usually comes to dominate all terminal stages of myelopoiesis. However, eventually, further mutations are acquired which cause progression to a rapidly fatal acute leukemia. To obtain new insights into the molecular perturbations that cause CML stem cells to initiate chronic phase disease, we generated Long Serial Analysis of Gene Expression (LongSAGE) libraries (∼200,000 tags/library) from extracts of highly purified lin−CD34+CD45RA−CD36−CD71−CD7−CD38+ and lin−CD34+CD45RA−CD36−CD71−CD7−CD38− normal bone marrow (BM) and G-CSF-mobilized peripheral blood (G-mPB) cells as well as cells from 3 chronic phase CML patients with predominantly Ph+/BCR-ABL+ cells in both of these very primitive cell subsets. Long term culture-initiating cell (LTC-IC) and direct colony-forming cell (CFC) assays performed on an aliquot of each of these cell populations showed the frequency of LTC-IC were 4 to 130-fold higher in the 34+38− cells than in the matching 34+38+ cells with the opposite trend for the CFCs. Comparison of the tags present in the pooled CML and pooled normal BM and G-mPB libraries revealed many differentially expressed genes. Real-time RT-PCR analysis of lin−CD34+ cells from 14 chronic phase CML patients and 3 normal BMs confirmed the differential expression of 14 candidate transcripts identified by SAGE (changes ranging from 3-fold lower to 80-fold higher in the CML cells, p<0.05). The altered levels of expression of 5 of these genes (i.e., beta-catenin, MLLT3, IL1R1, LY6E and GAS2) were highly correlated with the relative levels of BCR-ABL transcripts in the same cells (r ≥0.6). 5 of the 14 genes (IL1R1, vWF, SOX4, SELL and RHOB) were found to be differentially expressed in the 3-day post-transduction progeny of CD34+ cord blood cells exposed to a BCR/ABL-IRES-GFP vs a control GFP- lentivirus preparation (n=2). 3 (GAS2, DUSP1 and TP53BP2) were upregulated (5 to 11-fold) in imatinib-treated K562 cells (as compared to untreated K562 cells) but their expression remained unchanged in similarly treated KG1 cells (a primitive BCR-ABL-negative human AML cell line) providing further evidence that their deregulated expression is secondary to the kinase activity of p210BCR-ABL. In addition, from a comparison of the meta-library for the 34+38− CML cells with the meta-libraries for both the normal 34+38− and 34+38+ cells, the meta-library for the 34+38+ CML cells and most publicly accessible SAGE data, we were able to identify 65 novel tags in the 34+38− CML cells. 32 of these unique tags were located within conserved genomic regions and >2 kb away from known transcripts, and of these 32, 3 were confirmed to represent novel transcripts using a PCR approach. These results illustrate the potential of SAGE to reveal novel as well as known components in the transcriptomes of rare normal and cancer stem cell populations. Investigation of their roles in primitive human cells transduced with BCR-ABL and BCR-ABL+cell lines indicates the utility of these models for further delineation of the complex effects of BCR-ABL expression in chronic phase CML stem cells.

Chromosome instability, chromosome transcriptome, and clonal evolution of tumor cell populations

Chromosome instability and aneuploidy are hallmarks of cancer, but it is not clear how changes in the chromosomal content of a cell contribute to the malignant phenotype. Previously we have shown that we can readily isolate highly proliferative tumor cells and their revertants from highly invasive tumor cell populations, indicating how phenotypic shifting can contribute to malignant progression. Here we show that chromosome instability and changes in chromosome content occur with phenotypic switching. Further, we show that changes in the copy number of each chromosome quantitatively impose a proportional change in the chromosome transcriptome ratio. This correlation also applies to subchromosomal regions of derivative chromosomes. Importantly, we show that the changes in chromosome content and the transcriptome favor the expression of a large number of genes appropriate for the specific tumor phenotype. We conclude that chromosome instability generates the necessary chromosome diversity in the tumor cell populations and, therefore, the transcriptome diversity to allow for environment-facilitated clonal expansion and clonal evolution of tumor cell populations.

Presumptive pre‐sertoli cells express genes involved in cell proliferation and cell signalling during a critical window in early testis differentiation

In mammals, the pre-Sertoli cell of the male genital ridge is the first cell type to display sex specific differentiation and differential gene expression. The genetic cascade driving the differentiation of pre-Sertoli cells and ultimately testis formation is beginning to be unravelled, but many questions remain. A better understanding of the transcriptome of pre-Sertoli cells immediately after sex determination is essential in order to further understand this differentiation process. A mouse model expressing Red Fluorescent Protein (RFP) under the control of a hybrid mouse/pig SRY promoter (HybSRYp-RFP) was used to purify cells from embryonic day 12.0 (e12.0) male genital ridges. To compare the transcriptomes of HybSRYp-RFP cell populations versus age matched whole female genital ridges, RNA was extracted and used to generate molecular probes that were hybridized onto Affymetrix Mouse Genome 430 2.0 micro-arrays. The expression of genes considered markers for pre-Sertoli cells, including Sox9, Mis, Dhh and Fgf9 were identified within the HybSRYp-RFP expressing cell population, while markers for germ cells (Oct4, SSEA-1) and endothelial cells (Ntrk3) were not identified. In contrast, markers for ovarian somatic cell expression, including Fst and Bmp2, were identified as overexpressed within the ovarian cell population. In a general fashion, genes identified as 2.5-fold over expressed in HybSRYp-RFP expressing cells coded notably for cell signalling and extra cellular proteins. The expression of Sox10, Stc2, Fgf18, Fgf13 and Wnt6 were further characterized via whole mount in situ hybridization (WISH) on male and female genital ridges between e11.5 and e14.5. Sox10, Fgf18, Fgf13 and Stc2 gene expression was detected within the male genital ridges while Wnt6 was found diffusely within both the male and female genital ridges. These data represent the earliest comprehensive microarray expression analysis of purified presumptive pre-Sertoli cells available to date.

The Transcriptional Program of Terminal Granulocytic Differentiation in Man.

To characterize the transcriptional program that governs terminal granulocytic differentiation in vivo, we performed comprehensive microarray analyses of human bone marrow (BM) populations highly enriched in promyelocytes (PMs), myelocytes/metamyelocytes (MYs), and neutrophils (bm-PMNs). These analyses identified 11,310 genes involved in differentiation, of which 6,700 were differentially regulated including granule proteins and surface receptors previously not identified in neutrophils. Functional clustering demonstrated that differentially expressed genes were assigened to all of the major gene categories including apoptosis, cell cycle, chaperones, enzymes, immunity proteins, kinases, motor proteins, signal transducers, structural proteins, and transcription factors. In an attempt to define the developmental distance between PM, MY, and bm-PMN populations, we assessed the Pearsons correlation coefficient (γ) by correlating the transcriptomes (=all 44.760 probe sets, Affymetrix) of BM populations. The correlation coefficients among replicates of BM populations were in the range of 0.97–0.99. On the other hand the correlation coefficients between BM populations were 0.81 (γPM - MY), 0.79 (γMY - bm-PMN), and 0.52 (γPM - bm-PMN), and thus, provided a useful quantitative measure reflecting the hierarchical relationship between the three BM populations. Additional correlation of bm-PMN and peripheral blood-PMN transcriptomes revealed a high similarity of both populations (γbm-PMN - pb-PMN=0.95) indicating a termination of granulocytic differentiation in the BM microenvironment during steady state hematopoesis. Differentiation of PMs towards MYs was accompanied by a marked decline of proliferative and general cellular activity as defined by downregulation of E2F target genes, cyclin dependent kinases 2/4/6, and various metabolic, proteasomal, and mitochondrial genes. Expression patterns of apoptosis genes indicated death control by the p53-pathway in PMs and by death receptor pathways in bm-PMNs. Effector proteins critical for host defense were expressed successively throughout granulocytic differentiation, whereas receptors and receptor ligands essential for the activation of the host defense program were terminally upregulated in bm-PMNs. The upregulation of ligand-receptor pairs, which are defined inducers as well as target genes of NF-kB, suggests a constitutive activation of NF-kB in bm-PMNs by autocrine loops. Overall, these results define a granulocytic differentiation model governed by a highly coordinated fail-safe program, which promotes completion of differentiation before cells gain responsiveness towards activating stimuli that accompany infections.

Darwin's Naturalization Hypothesis Revisited

In The Origin of Species, Darwin (1859) drew attention to observations by Alphonse de Candolle (1855) that floras gain by naturalization far more species belonging to new genera than species belonging to native genera. Darwin (1859, p. 86) goes on to give a specific example: “In the last edition of Dr. Asa Gray’s ‘Manual of the Flora of the United States’ ... out of the 162 naturalised genera, no less than 100 genera are not there indigenous.” Darwin used these data to support his theory of intense competition between congeners, described only a few pages earlier: “As the species of the same genus usually have, though by no means invariably, much similarity in habits and constitution, and always in structure, the struggle will generally be more severe between them” (1859, p. 60). Darwin’s intriguing observations have recently attracted renewed interest, as comprehensive lists of naturalized plants have become available for various regions of the world. Two studies (Mack 1996; Rejmanek 1996, 1998) have concluded that naturalized floras provide some support for Darwin’s hypothesis, but only one of these studies used statistical tests. Analyses of additional floras are needed to test the generality of Darwin’s naturalization hypothesis. Mack (1996) tabulated data from six regional floras within the United States and noted that naturalized species more often belong to alien genera than native genera, with the curious exception of one region (New York). In addition to the possibility of strong competition between native and introduced congeners, Mack (1996) proposed that specialist native herbivores, or pathogens, may be

How Does Immigration Influence Local Adaptation? A Reexamination of a Familiar Paradigm

How Does Immigration Influence Local Adaptation? A Reexamination of a Familiar Paradigm