Welcome to Zebrafish Cell Development Landscape

Danio rerio Version2 scRNA-seq

Tissues
Zebrafish Cell Landscape View in gallery

Microwell-seq protocol

Materials

We taught a local company (G-Bio Biotechnology Co. Ltd) to provide Microwell-seq Kit/Service http://www.igeneseq.com/research_detail/dxbcp/microwell.html

Preparation of beads

Beads were washed twice with 1 ml of 100% EtOH and twice with 1 ml of TE/TW (10 mM Tris pH 8.0, 1 mM EDTA, 0.01% Tween). The bead pellet was resuspended in 1ml TE/TW for long-term storage at 4 ºC. For Microwell- Seq, 50 µl beads suspension was removed from the stock tube, washed twice with 500 µl of DPBS/EDTA (20 mM EDTA), then resuspended in 500 µl of DPBS/EDTA for a bead concentration (1:10, stock bead: bead for use). 500 µl diluted bead suspension was used for one experiment.

Fabrication of microwell array

PDMS washed once by 100% EtOH, twice by H2O and air drying. Agarose (5%) microwell arrays were cast from the PDMS before each experiment. Agarose microwell arrays could with DPBS/EDTA for short-term storage at 4ºC (in 10 days).

Protocol

Cell collection and lysis

Cell concentration should be carefully controlled in Microwell-seq. Both cell and bead concentrations were estimated by hemocytometer. The proper cell concentration is ~ 170,000 cells/ml (10% wells will be occupied by single cells). Evenly distributed cell suspension were pipetted onto microwell array (5 min~10 min) and extra cells were washed away. To eliminate cell doublets, the microwell plate was inspected under a microscope and microwell with overfull cell doublets can be discarded. Bead suspension(500 µl) was then loaded on the microwell plate which was placed on a magnet. Excess beads were washed away slowly. Cold lysis buffer (0.1 M Tris-HCL pH 7.5, 0.5 M LiCl, 1% LiSDS, 10 mM EDTA, and 5 mM dithiothreitol) was pipetted over the surface of the plate, and removed away after 10 min incubation. Then beads were collected and transferred to an RNase-free tube, washed twice with 200 µl 6X SSC, and once with 200 µl of 5x RT buffer (or 50 mM Tris buffer).

Reverse transcription

Then 20 µl RT mix was added. The RT mix contained dNTP mix (1mM), 100U SuperScript II reverse transcriptase, 1x Superscript II first-strand buffer (Takara, cat #2690A), 10 U Rnase Inhibitor (Lucigen, cat #30281-2 or other brands), 1M Betaine (Sigma, cat #B0300-1VL) and 6 mM MgCl2 (Ambion, cat #AM9539G), 1 µM TSO_RT_primer (see table 1). The beads were incubated at 42ºC for 90 minutes (gently shake beads every 6 min). The beads were then washed once with 200 µl TE-SDS, once with 200 µl TE-TW, and once with 1X exonuclease I buffer (or 10 mM Tris pH 7.5).

Exonuclease I treatment

After washing by 1X exonuclease I buffer, beads were suspended in 200 µl exonuclease I mix (containing 1x exonuclease I buffer and 1 U/µl exonuclease I (NEB, cat #B0293S)), and incubated at 37ºC for 45minutes to remove oligonucleotides which did not capture mRNA. Beads were pooled and washed once with 200 µl TE-SDS, once with 200 µl TE-TW and once with 200 µl 10 mM Tris pH 7.5 (or ultrapure distilled water).

cDNA amplification

Beads were then distributed to 5 tubes. Every tube was added with 20 µl PCR mix, which includes 1x HiFi HotStart Readymix (Kapa Biosystems, cat #KK2602) and 0.4 µM TSO_PCR primer (see Table 1). The PCR program was as follows: 98ºC for 3 min, 14 cycles of (98ºC 20sec, 67ºC 15sec, 72ºC 6 min), 72ºC 5 min and 4ºC hold. Shake the PCR tube to mix the beads (first cDNA template) at the end of 98ºC in the first four cycles (important!). After pooling all PCR products, AMPure XP beads (Beckman Coulter, cat #A63881) were used to purify the cDNA library.

Transposase fragmentation and selective PCR

Prepare 3'-end cDNA fragments for sequencing by Nextera XT tagmentation reactions (Illumina, cat#FC-131-1096). Fragmentation reaction was done following the instructions of manufacturer. We replace the index 2 primers (N5xx) in the kit with our P5 primer (see Table 1) to specially amplify fragments which contain the 3' end of transcript. Other fragments will form self-loops, which impede their binding to PCR primers. PCR program is as follows: 72ºC 3 min, 98ºC 30 sec, 13-15 cycles of (98ºC 15 sec, 60ºC 30 sec, 72ºC 3 min), 72ºC 5 min and 4ºC hold. To eliminate primer dimer and large size fragments, 0.6 x AMPure XP beads were then used to purify cDNA library. Then size distribution was analyzed on an Agilent 2100 bioanalyzer, a peak at around 400~700bp range should be observed. Finally, the libraries were sequenced on the Illumina NextSeq 500 using 3 µl of 100 µM Read1CustSeq (see Table 1) for priming of read 1.

Table 1 Primer used in Microwell-seq

Primer Sequence
TSO_RT primer AAGCAGTGGTATCAACGCAGAGTGAAT/rG/rG/iXNA_G
TSO_PCR primer AAGCAGTGGTATCAACGCAGAGT
P5 primer AATGATACGGCGACCACCGAGATCTACACGCCTGTCCGCGGAAGCAGTGGTATCAACGCAGAGT*A*C
Read1CustSeq GCCTGTCCGCGGAAGCAGTGGTATCAACGCAGAGTACGT
Download protocol for making microwells and barcoded beads

About ZCL & ZCDL

We used single-cell RNA sequencing to determine the cell type composition of major zebrafish organs and construct a basic scheme for the Zebrafish Cell Landscape (ZCL). To allow for public access of the resource, we created a Zebrafish Cell Landscape (ZCL) website at http://bis.zju.edu.cn/ZCL/. By integrating five different zebrafish development stages with ZCL, we developed the Zebrafish Cell Development Landscape (ZCDL), also known as ZCL 2.0.

ZCL website consists of seven web pages.
  1. Home: This page contains three sections. About describes the functions and update of the ZCL website, zebrafish tissue info gives a brief description of each organization of the sample, Microwell-seq protocol describes how Microwell-seq operates.

  2. Landscape: For ZCL 1.0, We analysed >280,000 single cells from >16 zebrafish tissues (2-4 replicates per tissue in general) and cultures. In a global view, the complete zebrafish tissue dataset are grouped into 55 major clusters. For ZCL 2.0, We analysed >1,360,000 single cells from 5 zebrafish stages. In a global view, the complete zebrafish tissue dataset are grouped into 135 major clusters. We created this page to achieve better visualization for global view. You can view each group, each tissue and each gene by choosing the optional box on the left. Landscape view provides global view on single cell level, Marker list provides marker genes for each cluster.

  3. Gallery: In the Gallery, you can download the single-cell digital gene expression (DGE) matrix and get information of the cell number and sample source for each data. The results of the clustering analysis with marker genes are also shared on the Gallery.

  4. Search: You can search genes in the sample you selected, the bar chart of the gene will be returned.

  5. Run-scZCL: You can upload the DGE matrix of RNA-seq data (scRNA-seq data or bulk RNA-seq data) and the scZCL can help you to identify cell types in your data. The results are showed by interactive heatmap and table in csv format. In heatmap, the row means the defined cell type, the column means the query sample, and the color of block indicates the strength of the correlation. In the result table, Pearson correlation coefficients between query samples (to be identified) and reference cell types are provided.

  6. FAQ: Some frequently asked questions are listed on this page.

  7. Contact: Download links for data and contact information for data providers are placed in this page. We also provide message board on this page.

Citation

Renying Wang†, Peijing Zhang†,*, Jingjing Wang†, Lifeng Ma†, Weigao E†, Shengbao Suo†, Mengmeng Jiang†, Jiaqi Li†, Haide Chen, Huiyu Sun, Lijiang Fei, Ziming Zhou, Yincong Zhou, Yao Chen, Weiqi Zhang, Xinru Wang, Yuqing Mei, Zhongyi Sun, Chengxuan Yu, Jikai Shao, Yuting Fu, Yanyu Xiao, Fang Ye, Xing Fang, Hanyu Wu, Qile Guo, Xiunan Fang, Xia Li, Xianzhi Gao, Dan Wang, Peng-Fei Xu, Rui Zeng, Gang Xu, Lijun Zhu, Lie Wang, Jing Qu, Dan Zhang, Hongwei Ouyang, He Huang, Ming Chen, Shyh-Chang NG*, Guang-Hui Liu*, Guo-Cheng Yuan*, Guoji Guo* and Xiaoping Han*. Constrction of a cross-species cell landscape at single-cell level. Nucleic Acids Research, 2022. DOI: 10.1093/nar/gkac633.

Mengmeng Jiang†, Yanyu Xiao†, Weigao E†, Lifeng Ma†, Jingjing Wang, Haide Chen, Ce Gao, Yuan Liao, Qile Guo, Jinrong Peng*, Xiaoping Han* and Guoji Guo*. Characterization of the Zebrafish Cell Landscape at Single-Cell Resolution. Frontiers in Cell and Developmental Biology, 2021; 9: 743421. DOI: 10.3389/fcell.2021.743421.

Downloads

For ZCL 2.0

RAW DGE: GSE198571, GSE178150
Figshare: https://figshare.com/s/1ab3c6d7648d12247eb2

For ZCL 1.0

RAW DGE: RAW DGE, GSE130487
Figshare: https://figshare.com/s/44b3fbaea9722e12304c
Tissue Annotations: Tissue Annotations
Fig1 H5ad: Fig1 H5ad
Fig1 Metadata: Fig1 Metadata
Fig1 Anno&marker: Fig1 Anno&marker

News

On August, 2022: Article online: Construction of a cross-species cell landscape at single-cell level. Nucleic Acids Research, 2022.
On March 24, 2022: Add zebrafish 72hpf, 3month and 22month data.
On March 24, 2022: Zebrafish Cell Development Landscape (ZCL Version 2.0) released.
On Oct 1, 2021: Article online: Jiang, M. et al. Characterization of the Zebrafish Cell Landscape at Single-Cell Resolution. Front Cell Dev Biol, 2021 .
On April 20, 2021: Add zebrafish embryo single-cell data (4hpf-24hpf).
On Dec 1, 2020: Whole zebrafish single-cell data released.
On June 30, 2020: ZCL Version 1.0 released.