Welcome to Drosophila Cell Landscape
Drosophila melanogaster scRNA-seq
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
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.
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).
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).
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).
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).
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.
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.
Primer | Sequence |
---|---|
TSO_RT primer | AAGCAGTGGTATCAACGCAGAGTGAAT/rG/rG/iXNA_G |
TSO_PCR primer | AAGCAGTGGTATCAACGCAGAGT |
P5 primer | AATGATACGGCGACCACCGAGATCTACACGCCTGTCCGCGGAAGCAGTGGTATCAACGCAGAGT*A*C |
Read1CustSeq | GCCTGTCCGCGGAAGCAGTGGTATCAACGCAGAGTACGT |
We used single-cell RNA sequencing to determine the cell type composition of Drosophila and construct a basic scheme for the Drosophila cell landscape (DCL). To allow for public access of the resource, we created a Drosophila cell landscape (DCL) website at http://bis.zju.edu.cn/DCL/.
DCL website consists of seven web pages.
1. Home: This page contains three sections. About describes the functions and update of the DCL website, Microwell-seq protocol describes how Microwell-seq operates.
2. Landscape: We analysed >400,000 single cells from 7 stages (Day11, Day16, Day21, Day26, Day31, Day36, Day41). In a global view, the complete Drosophila dataset are grouped into 85 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-scDCL: You can upload the DGE matrix of RNA-seq data (scRNA-seq data or bulk RNA-seq data) and the scDCL 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.
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.
Raw Data: GSE198831, GSE178149
Figshare: https://figshare.com/s/1ab3c6d7648d12247eb2
On August, 2022: Article online: Construction of a cross-species cell landscape at single-cell level. Nucleic Acids Research, 2022.
On April 21, 2021: Add Drosophila tissue single-cell data in gallery.
On April 10, 2021: Add Drosophila 7 stage single-cell data in gallery.
On December 1, 2020: DCL Version 1.0 released.