Click
here to close Hello! We notice that
you are using Internet Explorer, which is not supported by Echinobase
and may cause the site to display incorrectly. We suggest using a
current version of Chrome,
FireFox,
or Safari.
Genome Res
2014 May 01;245:860-8. doi: 10.1101/gr.167668.113.
Show Gene links
Show Anatomy links
General approach for in vivo recovery of cell type-specific effector gene sets.
Barsi JC
,
Tu Q
,
Davidson EH
.
???displayArticle.abstract???
Differentially expressed, cell type-specific effector gene sets hold the key to multiple important problems in biology, from theoretical aspects of developmental gene regulatory networks (GRNs) to various practical applications. Although individual cell types of interest have been recovered by various methods and analyzed, systematic recovery of multiple cell type-specific gene sets from whole developing organisms has remained problematic. Here we describe a general methodology using the sea urchin embryo, a material of choice because of the large-scale GRNs already solved for this model system. This method utilizes the regulatory states expressed by given cells of the embryo to define cell type and includes a fluorescence activated cell sorting (FACS) procedure that results in no perturbation of transcript representation. We have extensively validated the method by spatial and qualitative analyses of the transcriptome expressed in isolated embryonic skeletogenic cells and as a consequence, generated a prototypical cell type-specific transcriptome database.
???displayArticle.pubmedLink???
24604781
???displayArticle.pmcLink???PMC4009615 ???displayArticle.link???Genome Res ???displayArticle.grants???[+]
Figure 1. Steps required for attaining cell type-specific transcriptomes. (1) Microinject a group of zygotes with an artificial chromosome genetically engineered to express a fluorescent reporter in the cell type of interest; (2) culture transgenic zygotes until they have reached the desired developmental stage, then disaggregate the embryos to individual cells; (3) run the cells through flow cytometry in order to segregate populations of cells by way of FACS; and (4) isolate, amplify, sequence, and quantify the pool of mRNA from each population in order to calculate relative enrichment for each gene. (AU) Arbitrary units; (BACs) bacterial artificial chromosomes.
Figure 2. Assemblage of control experiments. (A) nCounter analysis system measurements quantifying the expression level of 181 regulatory mRNAs used to assess the gene expression profile of disaggregated embryos. Red line represents a twofold change in gene expression from moderate abundance onward. (A1) Cells from disaggregated embryos collected by flow cytometry show no change in gene expression when compared to undisturbed whole embryos. (A2) Cells from disaggregated embryos collected by flow cytometry show no change in gene expression among biological replicates. (B) Fixed Sp blastula with PMCs labeled in purple. Labeling was achieved by whole mount RNA in situ hybridization of cah10l mRNA. (C) Live transgenic Sp blastula carrying an artificial chromosome that contains GFP under control of the cis-regulatory apparatus governing tbr expression. (C1) DIC image of transgenic blastula. (C2) Epifluorescent image of the same specimen shown in C1 reveals that GFP is exclusively expressed in PMCs. (C3) Composite image generated by merging C1 with C2. (D,E) Flow cytometry data in the form of 2% probability contour plots. Events falling outside of the lowest contour have been depicted as dots within the graph. A vertical red line demarcates the value beneath which an event likely reflects a technical artifact. (D) The exclusion of nonviable cells required that threshold parameters be calibrated for 7-AAD fluorescence. (D1) Data reflects the variation in size among the cells utilized for 7-AAD calibration. (D2) Untreated cells were used in determining baseline fluorescence, depicted as a horizontal red line. (D3) Following 7-AAD treatment, cells observed to fluoresce above baseline were excluded from further analysis. (E) Segregation of cell populations by FACS required that threshold parameters be calibrated for GFP fluorescence. (E1) Data reflects the variation in size among the cells utilized for GFP calibration. (E2) Cells from uninjected embryos were used in determining baseline fluorescence, depicted as a horizontal red line. (E3) Cells from transgenic embryos were then segregated by FACS, according to their fluorescence relative to baseline. (AU) Arbitrary units.
Figure 3. Flow cytometry gating strategy across replicates. (A1–C3) Flow cytometry data in the form of 2% probability contour plots. In all instances, events falling outside the lowest contour have been depicted as dots within the graph. (A,B,C) Each series reflects data from one of the three replicates. (A) Replicate #1 obtained by using an alx1:GFP BAC. Replicate #2 (B) and Replicate #3 (C) obtained by using a tbr:GFP BAC. (A1,B1,C1) Each graph depicts all the events detected by flow cytometry. Events enclosed within the polygonal red line were visually corroborated to constitute individual cells, hence promoted to a second round of analysis. (A2,B2,C2) Each graph reflects the fraction of cells that have incorporated 7-AAD. Consequently, these were excluded from further study. A polygonal red line encloses the cell population promoted to a third and final round of analysis. (A3,B3,C3) Each graph reflects the fraction of viable cells segregated by GFP FACS. Data points shown above the horizontal red line represent PMCs, whereas those below represent a heterogeneous population containing all cell types. Percentage displayed in red at the corner of each quadrant. A vertical red line demarcates the value beneath which an event likely reflects a technical artifact. (AU) Arbitrary units.
Figure 4. Cell type-specific transcriptome. (A–D) Scatterplots compare transcriptome data unique to PMCs (GFP+) with that of all cell types (GFP−). Each data point reflects the relative mRNA abundance of a gene as estimated by the number of sequencing reads that map to its locus. The values in the graph reflect measurements taken across three independent replicates. Data points that fall above the diagonal indicate augmented levels of expression among PMCs, relative to other cell types. Genes statistically determined to be differentially expressed are shown as data points colored either red or blue to indicate an associated P-value of less than 0.01 or 0.05, respectively. Data points corresponding to genes of interest have been encircled in order to distinguish them from the rest of the data set. In select cases, a label revealing gene identity accompanies encircled data points. (A) The cohort of transcription factors enriched among PMCs is in accordance with, and expands upon, the known PMC regulatory state. (B) Genes directly involved in biomineralization are overrepresented in the catalog of PMC enriched transcripts, illuminating the skeletogenic fate these cells acquire. (C) Marked data points reflect the subset of PMC enriched transcripts for which spatial expression has been corroborated by RNA in situ hybridization. (D) Marked data points reflect genes reported in the literature to be enriched within PMCs.
Figure 5. Corroboration of cell type specificity. (A1–A4) Four of the genes tested for cell type-specific expression. Fixed Sp blastula with PMCs labeled in purple. Labeling was achieved by whole mount RNA in situ hybridization of cah10l (A1), msp130r1 (A2), p16 (A3), and 3apcol (A4). (B) Scatterplot identical to those described in Figure 4. Data points corresponding to the genes shown in A1–A4 have been encircled in order to distinguish them from other differentially expressed genes and labeled accordingly. Supplemental Table 1 lists all genes corroborated in a similar manner and the supporting data are shown in Supplemental Figures 4–12.
Anders,
Differential expression analysis for sequence count data.
2010, Pubmed
Anders,
Differential expression analysis for sequence count data.
2010,
Pubmed
Ben-Tabou de-Leon,
Gene regulatory control in the sea urchin aboral ectoderm: spatial initiation, signaling inputs, and cell fate lockdown.
2013,
Pubmed
,
Echinobase
Cannon,
Global analysis of the haematopoietic and endothelial transcriptome during zebrafish development.
2013,
Pubmed
Damle,
Precise cis-regulatory control of spatial and temporal expression of the alx-1 gene in the skeletogenic lineage of s. purpuratus.
2011,
Pubmed
,
Echinobase
Dobin,
STAR: ultrafast universal RNA-seq aligner.
2013,
Pubmed
Duloquin,
Localized VEGF signaling from ectoderm to mesenchyme cells controls morphogenesis of the sea urchin embryo skeleton.
2007,
Pubmed
,
Echinobase
Estrada,
An integrated strategy for analyzing the unique developmental programs of different myoblast subtypes.
2006,
Pubmed
Geiss,
Direct multiplexed measurement of gene expression with color-coded probe pairs.
2008,
Pubmed
Haenni,
Analysis of C. elegans intestinal gene expression and polyadenylation by fluorescence-activated nuclei sorting and 3'-end-seq.
2012,
Pubmed
Knapp,
Recombinant sea urchin vascular endothelial growth factor directs single-crystal growth and branching in vitro.
2012,
Pubmed
,
Echinobase
Lee,
A highly efficient Escherichia coli-based chromosome engineering system adapted for recombinogenic targeting and subcloning of BAC DNA.
2001,
Pubmed
Li,
Encoding regulatory state boundaries in the pregastrular oral ectoderm of the sea urchin embryo.
2014,
Pubmed
,
Echinobase
Li,
Direct and indirect control of oral ectoderm regulatory gene expression by Nodal signaling in the sea urchin embryo.
2012,
Pubmed
,
Echinobase
Li,
New regulatory circuit controlling spatial and temporal gene expression in the sea urchin embryo oral ectoderm GRN.
2013,
Pubmed
,
Echinobase
Livingston,
A genome-wide analysis of biomineralization-related proteins in the sea urchin Strongylocentrotus purpuratus.
2006,
Pubmed
,
Echinobase
Longabaugh,
BioTapestry: a tool to visualize the dynamic properties of gene regulatory networks.
2012,
Pubmed
Materna,
Diversification of oral and aboral mesodermal regulatory states in pregastrular sea urchin embryos.
2013,
Pubmed
,
Echinobase
Oliveri,
Global regulatory logic for specification of an embryonic cell lineage.
2008,
Pubmed
,
Echinobase
Peter,
Modularity and design principles in the sea urchin embryo gene regulatory network.
2009,
Pubmed
,
Echinobase
Peter,
A gene regulatory network controlling the embryonic specification of endoderm.
2011,
Pubmed
,
Echinobase
Peter,
Predictive computation of genomic logic processing functions in embryonic development.
2012,
Pubmed
,
Echinobase
Poustka,
On the origin of the chordate central nervous system: expression of onecut in the sea urchin embryo.
2004,
Pubmed
,
Echinobase
Rafiq,
The genomic regulatory control of skeletal morphogenesis in the sea urchin.
2012,
Pubmed
,
Echinobase
Ransick,
Detection of mRNA by in situ hybridization and RT-PCR.
2004,
Pubmed
,
Echinobase
Spencer,
A spatial and temporal map of C. elegans gene expression.
2011,
Pubmed
Tu,
Gene structure in the sea urchin Strongylocentrotus purpuratus based on transcriptome analysis.
2012,
Pubmed
,
Echinobase
Tu,
Quantitative developmental transcriptomes of the sea urchin Strongylocentrotus purpuratus.
2014,
Pubmed
,
Echinobase
Wahl,
The cis-regulatory system of the tbrain gene: Alternative use of multiple modules to promote skeletogenic expression in the sea urchin embryo.
2009,
Pubmed
,
Echinobase
Wu,
Ingression of primary mesenchyme cells of the sea urchin embryo: a precisely timed epithelial mesenchymal transition.
2007,
Pubmed
,
Echinobase
Zhu,
A large-scale analysis of mRNAs expressed by primary mesenchyme cells of the sea urchin embryo.
2001,
Pubmed
,
Echinobase
de-Leon,
Information processing at the foxa node of the sea urchin endomesoderm specification network.
2010,
Pubmed
,
Echinobase