Rapid changes in chromatin structure during dedifferentiation of primary hepatocytes in vitro.

Primary hepatocytes are widely used in the pharmaceutical industry to screen drug candidates for hepatotoxicity, but hepatocytes quickly dedifferentiate and lose their mature metabolic function in culture. Attempts have been made to better recapitulate the in vivo liver environment in culture, but the full spectrum of signals required to maintain hepatocyte function ex vivo remains elusive. To elucidate molecular changes that accompany, and may contribute to dedifferentiation of hepatocytes ex vivo, we performed lineage tracing and comprehensive profiling of alterations in their gene expression profiles and chromatin landscape during culture. First, using genetically tagged hepatocytes we demonstrate that expression of the fetal gene alpha-fetoprotein in cultured hepatocytes comes from cells that previously expressed the mature gene albumin, and not from a population of albumin-negative precursor cells, proving mature hepatocytes undergo true dedifferentiation in culture. Next we studied the dedifferentiation process in detail through bulk RNA-sequencing of hepatocytes cultured over an extended period. We identified three distinct phases of dedifferentiation: an early phase, where mature hepatocyte genes are rapidly downregulated in a matter of hours; a middle phase, where fetal genes are activated; and a late phase, where initially rare contaminating non-parenchymal cells proliferate, taking over the culture. Lastly, to better understand the signaling events that result in the rapid downregulation of mature genes in hepatocytes, we examined changes in chromatin accessibility in these cells during the first 24 h of culture using Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). We find that drastic and rapid changes in chromatin accessibility occur immediately upon the start of culture. Using binding motif analysis of the areas of open chromatin sharing similar temporal profiles, we identify several candidate transcription factors potentially involved in the dedifferentiation of primary hepatocytes in culture.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization.

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Interspecies chimeric conditions affect the developmental rate of human pluripotent stem cells.

Human pluripotent stem cells hold significant promise for regenerative medicine. However, long differentiation protocols and immature characteristics of stem cell-derived cell types remain challenges to the development of many therapeutic applications. In contrast to the slow differentiation of human stem cells in vitro that mirrors a nine-month gestation period, mouse stem cells develop according to a much faster three-week gestation timeline. Here, we tested if co-differentiation with mouse pluripotent stem cells could accelerate the differentiation speed of human embryonic stem cells. Following a six-week RNA-sequencing time course of neural differentiation, we identified 929 human genes that were upregulated earlier and 535 genes that exhibited earlier peaked expression profiles in chimeric cell cultures than in human cell cultures alone. Genes with accelerated upregulation were significantly enriched in Gene Ontology terms associated with neurogenesis, neuron differentiation and maturation, and synapse signaling. Moreover, chimeric mixed samples correlated with in utero human embryonic samples earlier than human cells alone, and acceleration was dose-dependent on human-mouse co-culture ratios. The altered gene expression patterns and developmental rates described in this report have implications for accelerating human stem cell differentiation and the use of interspecies chimeric embryos in developing human organs for transplantation.

Reproducibility across single-cell RNA-seq protocols for spatial ordering analysis.

As newer single-cell protocols generate increasingly more cells at reduced sequencing depths, the value of a higher read depth may be overlooked. Using data from three different single-cell RNA-seq protocols that lend themselves to having either higher read depth (Smart-seq) or many cells (MARS-seq and 10X), we evaluate their ability to recapitulate biological signals in the context of spatial reconstruction. Overall, we find gene expression profiles after spatial reconstruction analysis are highly reproducible between datasets despite being generated by different protocols and using different computational algorithms. While UMI-based protocols such as 10X and MARS-seq allow for capturing more cells, Smart-seq's higher sensitivity and read-depth allow for analysis of lower expressed genes and isoforms. Additionally, we evaluate trade-offs for each protocol by performing subsampling analyses and find that optimizing the balance between sequencing depth and number of cells within a protocol is necessary for efficient use of resources. Our analysis emphasizes the importance of selecting a protocol based on the biological questions and features of interest.

Automated minute scale RNA-seq of pluripotent stem cell differentiation reveals early divergence of human and mouse gene expression kinetics.

Pluripotent stem cells retain the developmental timing of their species of origin in vitro, an observation that suggests the existence of a cell-intrinsic developmental clock, yet the nature and machinery of the clock remain a mystery. We hypothesize that one possible component may lie in species-specific differences in the kinetics of transcriptional responses to differentiation signals. Using a liquid-handling robot, mouse and human pluripotent stem cells were exposed to identical neural differentiation conditions and sampled for RNA-sequencing at high frequency, every 4 or 10 minutes, for the first 10 hours of differentiation to test for differences in transcriptomic response rates. The majority of initial transcriptional responses occurred within a rapid window in the first minutes of differentiation for both human and mouse stem cells. Despite similarly early onsets of gene expression changes, we observed shortened and condensed gene expression patterns in mouse pluripotent stem cells compared to protracted trends in human pluripotent stem cells. Moreover, the speed at which individual genes were upregulated, as measured by the slopes of gene expression changes over time, was significantly faster in mouse compared to human cells. These results suggest that downstream transcriptomic response kinetics to signaling cues are faster in mouse versus human cells, and may offer a partial account for the vast differences in developmental rates across species.

RNA-Seq of Human Neural Progenitor Cells Exposed to Lead (Pb) Reveals Transcriptome Dynamics, Splicing Alterations and Disease Risk Associations.

Lead (Pb) is a well-known toxicant, especially for the developing nervous system, albeit the mechanism is largely unknown. In this study, we use time series RNA-seq to conduct a genome-wide survey of the transcriptome response of human embryonic stem cell-derived neural progenitor cells to lead treatment. Using a dynamic time warping algorithm coupled with statistical tests, we find that lead can either accelerate or decelerate the expression of specific genes during the time series. We further show that lead disrupts a neuron- and brain-specific splicing factor NOVA1 regulated splicing network. Using lead induced transcriptome change signatures, we predict several known and novel disease risks under lead exposure. The findings in this study will allow a better understanding of the mechanism of lead toxicity, facilitate the development of diagnostic biomarkers and treatment for lead exposure, and comprise a highly valuable resource for environmental toxicology. Our study also demonstrates that a human (embryonic stem) cell-derived system can be used for studying the mechanism of toxicants, which can be useful for drug or compound toxicity screens and safety assessment.

Uniform neural tissue models produced on synthetic hydrogels using standard culture techniques.

The aim of the present study was to test sample reproducibility for model neural tissues formed on synthetic hydrogels. Human embryonic stem (ES) cell-derived precursor cells were cultured on synthetic poly(ethylene glycol) (PEG) hydrogels to promote differentiation and self-organization into model neural tissue constructs. Neural progenitor, vascular, and microglial precursor cells were combined on PEG hydrogels to mimic developmental timing, which produced multicomponent neural constructs with 3D neuronal and glial organization, organized vascular networks, and microglia with ramified morphologies. Spearman's rank correlation analysis of global gene expression profiles and a comparison of coefficient of variation for expressed genes demonstrated that replicate neural constructs were highly uniform to at least day 21 for samples from independent experiments. We also demonstrate that model neural tissues formed on PEG hydrogels using a simplified neural differentiation protocol correlated more strongly to in vivo brain development than samples cultured on tissue culture polystyrene surfaces alone. These results provide a proof-of-concept demonstration that 3D cellular models that mimic aspects of human brain development can be produced from human pluripotent stem cells with high sample uniformity between experiments by using standard culture techniques, cryopreserved cell stocks, and a synthetic extracellular matrix. Impact statement Pluripotent stem (PS) cells have been characterized by an inherent ability to self-organize into 3D "organoids" resembling stomach, intestine, liver, kidney, and brain tissues, offering a potentially powerful tool for modeling human development and disease. However, organoid formation must be quantitatively reproducible for applications such as drug and toxicity screening. Here, we report a strategy to produce uniform neural tissue constructs with reproducible global gene expression profiles for replicate samples from multiple experiments.

Species-specific developmental timing is maintained by pluripotent stem cells ex utero.

How species-specific developmental timing is controlled is largely unknown. By following human embryonic stem (ES) cell and mouse epiblast stem (EpiS) cell differentiation through detailed RNA-sequencing time courses, here we show that pluripotent stem cells closely retain in vivo species-specific developmental timing in vitro. In identical neural differentiation conditions in vitro, gene expression profiles are accelerated in mouse EpiS cells compared to human ES cells with relative rates of differentiation closely reflecting the rates of progression through the Carnegie stages in utero. Dynamic Time Warping analysis identified 3389 genes that were regulated more quickly in mouse EpiS cells and identified none that were regulated more quickly in human ES cells. Interestingly, we also find that human ES cells differentiated in teratomas maintain the same rate of differentiation observed in vitro in spite of being grown in a mouse host. These results suggest the existence of a cell autonomous, species-specific developmental clock that pluripotent stem cells maintain even out of context of an intact embryo.

Human pluripotent stem cell-derived neural constructs for predicting neural toxicity.

Human pluripotent stem cell-based in vitro models that reflect human physiology have the potential to reduce the number of drug failures in clinical trials and offer a cost-effective approach for assessing chemical safety. Here, human embryonic stem (ES) cell-derived neural progenitor cells, endothelial cells, mesenchymal stem cells, and microglia/macrophage precursors were combined on chemically defined polyethylene glycol hydrogels and cultured in serum-free medium to model cellular interactions within the developing brain. The precursors self-assembled into 3D neural constructs with diverse neuronal and glial populations, interconnected vascular networks, and ramified microglia. Replicate constructs were reproducible by RNA sequencing (RNA-Seq) and expressed neurogenesis, vasculature development, and microglia genes. Linear support vector machines were used to construct a predictive model from RNA-Seq data for 240 neural constructs treated with 34 toxic and 26 nontoxic chemicals. The predictive model was evaluated using two standard hold-out testing methods: a nearly unbiased leave-one-out cross-validation for the 60 training compounds and an unbiased blinded trial using a single hold-out set of 10 additional chemicals. The linear support vector produced an estimate for future data of 0.91 in the cross-validation experiment and correctly classified 9 of 10 chemicals in the blinded trial.

A cost-effective RNA sequencing protocol for large-scale gene expression studies.

RNA sequencing has increasingly become an indispensable tool for biological research. While sequencing costs have fallen dramatically in recent years, the current cost of RNA sequencing, nonetheless, remains a barrier to even more widespread adoption. Here, we present a simple RNA sequencing protocol with substantially reduced costs. This protocol uses as little as 10 ng of total RNA, allows multiplex sequencing of up to 96 samples per lane, and is strand specific. Extensive validation using human embryonic stem cells showed high consistency between technical replicates at various multiplexing levels.

Comparative RNA-seq analysis in the unsequenced axolotl: the oncogene burst highlights early gene expression in the blastema.

The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.

Single-round, multiplexed antibody mimetic design through mRNA display.

In a single round: By combining the high-efficiency enrichment through the continuous-flow magnetic separation (CFMS) technique with the analytical power of next-generation sequencing, the generation of antibody mimetics with a single round of mRNA display is made possible. This approach eliminates iterative selection cycles and provides a path to fully automated ligand generation (see picture).

Single read and paired end mRNA-Seq Illumina libraries from 10 nanograms total RNA.

Whole transcriptome sequencing by mRNA-Seq is now used extensively to perform global gene expression, mutation, allele-specific expression and other genome-wide analyses. mRNA-Seq even opens the gate for gene expression analysis of non-sequenced genomes. mRNA-Seq offers high sensitivity, a large dynamic range and allows measurement of transcript copy numbers in a sample. Illumina's genome analyzer performs sequencing of a large number (> 10(7)) of relatively short sequence reads (< 150 bp).The "paired end" approach, wherein a single long read is sequenced at both its ends, allows for tracking alternate splice junctions, insertions and deletions, and is useful for de novo transcriptome assembly. One of the major challenges faced by researchers is a limited amount of starting material. For example, in experiments where cells are harvested by laser micro-dissection, available starting total RNA may measure in nanograms. Preparation of mRNA-Seq libraries from such samples have been described(1, 2) but involves significant PCR amplification that may introduce bias. Other RNA-Seq library construction procedures with minimal PCR amplification have been published(3, 4) but require microgram amounts of starting total RNA. Here we describe a protocol for the Illumina Genome Analyzer II platform for mRNA-Seq sequencing for library preparation that avoids significant PCR amplification and requires only 10 nanograms of total RNA. While this protocol has been described previously and validated for single-end sequencing(5), where it was shown to produce directional libraries without introducing significant amplification bias, here we validate it further for use as a paired end protocol. We selectively amplify polyadenylated messenger RNAs from starting total RNA using the T7 based Eberwine linear amplification method, coined "T7LA" (T7 linear amplification). The amplified poly-A mRNAs are fragmented, reverse transcribed and adapter ligated to produce the final sequencing library. For both single read and paired end runs, sequences are mapped to the human transcriptome(6) and normalized so that data from multiple runs can be compared. We report the gene expression measurement in units of transcripts per million (TPM), which is a superior measure to RPKM when comparing samples(7).

Proteomic and phosphoproteomic comparison of human ES and iPS cells.

Combining high-mass-accuracy mass spectrometry, isobaric tagging and software for multiplexed, large-scale protein quantification, we report deep proteomic coverage of four human embryonic stem cell and four induced pluripotent stem cell lines in biological triplicate. This 24-sample comparison resulted in a very large set of identified proteins and phosphorylation sites in pluripotent cells. The statistical analysis afforded by our approach revealed subtle but reproducible differences in protein expression and protein phosphorylation between embryonic stem cells and induced pluripotent cells. Merging these results with RNA-seq analysis data, we found functionally related differences across each tier of regulation. We also introduce the Stem Cell-Omics Repository (SCOR), a resource to collate and display quantitative information across multiple planes of measurement, including mRNA, protein and post-translational modifications.

Chemically defined conditions for human iPSC derivation and culture.

We re-examine the individual components for human embryonic stem cell (ESC) and induced pluripotent stem cell (iPSC) culture and formulate a cell culture system in which all protein reagents for liquid media, attachment surfaces and splitting are chemically defined. A major improvement is the lack of a serum albumin component, as variations in either animal- or human-sourced albumin batches have previously plagued human ESC and iPSC culture with inconsistencies. Using this new medium (E8) and vitronectin-coated surfaces, we demonstrate improved derivation efficiencies of vector-free human iPSCs with an episomal approach. This simplified E8 medium should facilitate both the research use and clinical applications of human ESCs and iPSCs and their derivatives, and should be applicable to other reprogramming methods.