Tamar Hashimshony

CEL-Seq: single-cell RNA-Seq by multiplexed linear amplification.

High-throughput sequencing has allowed for unprecedented detail in gene expression analyses, yet its efficient application to single cells is challenged by the small starting amounts of RNA. We have developed CEL-Seq, a method for overcoming this limitation by barcoding and pooling samples before linearly amplifying mRNA with the use of one round of in vitro transcription. We show that CEL-Seq gives more reproducible, linear, and sensitive results than a PCR-based amplification method. We demonstrate the power of this method by studying early C. elegans embryonic development at single-cell resolution. Differential distribution of transcripts between sister cells is seen as early as the two-cell stage embryo, and zygotic expression in the somatic cell lineages is enriched for transcription factors. The robust transcriptome quantifications enabled by CEL-Seq will be useful for transcriptomic analyses of complex tissues containing populations of diverse cell types.

The role of DNA methylation in setting up chromatin structure during development

DNA methylation inhibits gene expression in animal cells, probably by affecting chromatin structure. Biochemical studies suggest that this process may be mediated by methyl-specific binding proteins that recruit enzymatic machinery capable of locally altering histone modification. To test whether DNA methylation actually has a role in the assembly of chromatin during normal development, we used cell transfection and a transgene construct genetically programmed to be either methylated or unmethylated in all cell types of the mouse. Chromatin immunoprecipitation (ChIP) analysis shows that the presence of DNA methylation brings about the deacetylation of histone H4 and methylation of Lys9 of histone H3 (H3 Lys9) and prevents methylation of Lys4 of histone H3 (H3 Lys4), thus generating a structure identical to that of methylated sequences in the genome. These results indicate that the methylation pattern established in early embryogenesis is profoundly important in setting up the structural profile of the genome.

DNA methylation models histone acetylation

One of the main determinants of chromatin structure is histone acetylation. Local chromosomal acetylation can be regulated dynamically, both through the involvement of transactivating factors with intrinsic histone acetylase activity, and through the recruitment of deacetylase complexes that repress gene expression. Histone acetylation status is transiently modified from one state to another in response to physiological changes operating in the cell. It is not yet known, however, how the basic histone acetylation profiles on tissue-specific and housekeeping gene sequences are established during normal development. Here we show that DNA methylation takes part in this process by inducing decreased levels of chromatin acetylation.

Establishment of transcriptional competence in early and late S phase

In animal cells, the process of DNA replication takes place in a programmed manner, with each gene region designated to replicate at a fixed time slot in S phase. Housekeeping genes undergo replication in the first half of S phase in all cell types, whereas the replication of many tissue specific genes is developmentally controlled, being late in most tissues but early in the tissue of expression. Here we employ nuclear DNA injection as an experimental system to test whether this phenomenon is due to differences in the ability to set up transcriptional competence during S phase. Our results show that, regardless of sequence, exogenous genes are a better template for transcription when injected into nuclei of cells in early as opposed to late S phase, and this expression state, once initiated, is preserved after cell division. DNA injected in late S phase is apparently repressed because it is packaged into chromatin containing deacetylated histones, and the same is true for late replicating chromosomal DNA. These findings suggest a mechanistic connection between replication timing and gene expression that might help to explain how epigenetic states can be maintained in vivo.

Developmental Milestones Punctuate Gene Expression in the Caenorhabditis Embryo

A fundamental question in developmental biology relates to the connection between morphological stages and their underlying molecular activity. Here we demonstrate that, at the molecular level, embryonic development in five Caenorhabditis species proceeds through two distinct milestones in which the transcriptome is resistant to differences in species-specific developmental timings. By comparing the complete protein-coding transcriptomes of individually timed embryos across ten morphological markers, we found that developmental milestones can be characterized by their expression dynamics and activation of key developmental regulators. This approach led us to discover the nematode phylotypic stage and to show that in chordates and arthropods it is represented as two distinct stages, suggesting that animal body plans might evolve by uncoupling and elaboration on formerly synchronous processes.

CEL-Seq2: sensitive highly-multiplexed single-cell RNA-Seq.

Single-cell transcriptomics requires a method that is sensitive, accurate, and reproducible. Here, we present CEL-Seq2, a modified version of our CEL-Seq method, with threefold higher sensitivity, lower costs, and less hands-on time. We implemented CEL-Seq2 on Fluidigm’s C1 system, providing its first single-cell, on-chip barcoding method, and we detected gene expression changes accompanying the progression through the cell cycle in mouse fibroblast cells. We also compare with Smart-Seq to demonstrate CEL-Seq2’s increased sensitivity relative to other available methods. Collectively, the improvements make CEL-Seq2 uniquely suited to single-cell RNA-Seq analysis in terms of economics, resolution, and ease of use.

An upstream repressor element plays a role in Igf2 imprinting

The imprinted Igf2 gene is associated with a small upstream region that is differentially methylated on the active paternal allele. We have identified a repressor element within this sequence and shown that repression is probably mediated through a trans‐ acting factor, GCF2. DNA methylation of this site abrogates both protein binding and repressor activity. Targeting experiments demonstrate that this element plays a role in the repression of the maternal Igf2 gene in vivo.

Lymphatic vessels arise from specialized angioblasts within a venous niche

How cells acquire their fate is a fundamental question in developmental and regenerative biology. Multipotent progenitors undergo cell-fate restriction in response to cues from the microenvironment, the nature of which is poorly understood. In the case of the lymphatic system, venous cells from the cardinal vein are thought to generate lymphatic vessels through trans-differentiation. Here we show that in zebrafish, lymphatic progenitors arise from a previously uncharacterized niche of specialized angioblasts within the cardinal vein, which also generates arterial and venous fates. We further identify Wnt5b as a novel lymphatic inductive signal and show that it also promotes the ‘angioblast-to-lymphatic’ transition in human embryonic stem cells, suggesting that this process is evolutionarily conserved. Our results uncover a novel mechanism of lymphatic specification, and provide the first characterization of the lymphatic inductive niche. More broadly, our findings highlight the cardinal vein as a heterogeneous structure, analogous to the haematopoietic niche in the aortic floor.

Spatiotemporal transcriptomics reveals the evolutionary history of the endoderm germ layer

The concept of germ layers has been one of the foremost organizing principles in developmental biology, classification, systematics and evolution for 150 years (refs 1, 2, 3). Of the three germ layers, the mesoderm is found in bilaterian animals but is absent in species in the phyla Cnidaria and Ctenophora, which has been taken as evidence that the mesoderm was the final germ layer to evolve. The origin of the ectoderm and endoderm germ layers, however, remains unclear, with models supporting the antecedence of each as well as a simultaneous origin. Here we determine the temporal and spatial components of gene expression spanning embryonic development for all Caenorhabditis elegans genes and use it to determine the evolutionary ages of the germ layers. The gene expression program of the mesoderm is induced after those of the ectoderm and endoderm, thus making it the last germ layer both to evolve and to develop. Strikingly, the C. elegans endoderm and ectoderm expression programs do not co-induce; rather the endoderm activates earlier, and this is also observed in the expression of endoderm orthologues during the embryology of the frog Xenopus tropicalis, the sea anemone Nematostella vectensis and the sponge Amphimedon queenslandica. Querying the phylogenetic ages of specifically expressed genes reveals that the endoderm comprises older genes. Taken together, we propose that the endoderm program dates back to the origin of multicellularity, whereas the ectoderm originated as a secondary germ layer freed from ancestral feeding functions.

The mid-developmental transition and the evolution of animal body plans

Animals are grouped into ~35 ‘phyla’ based upon the notion of distinct body plans. Morphological and molecular analyses have revealed that a stage in the middle of development—known as the phylotypic period—is conserved among species within some phyla. Although these analyses provide evidence for their existence, phyla have also been criticized as lacking an objective definition, and consequently based on arbitrary groupings of animals. Here we compare the developmental transcriptomes of ten species, each annotated to a different phylum, with a wide range of life histories and embryonic forms. We find that in all ten species, development comprises the coupling of early and late phases of conserved gene expression. These phases are linked by a divergent ‘mid-developmental transition’ that uses species-specific suites of signalling pathways and transcription factors. This mid-developmental transition overlaps with the phylotypic period that has been defined previously for three of the ten phyla, suggesting that transcriptional circuits and signalling mechanisms active during this transition are crucial for defining the phyletic body plan and that the mid-developmental transition may be used to define phylotypic periods in other phyla. Placing these observations alongside the reported conservation of mid-development within phyla, we propose that a phylum may be defined as a collection of species whose gene expression at the mid-developmental transition is both highly conserved among them, yet divergent relative to other species.