This video shows the basic steps for performing whole transcriptome analysis on dissected chick embryonic spinal cord samples after transfection with in ovo electroporation.
In ovo electroporation of the chick neural tube is a fast and inexpensive method for identification of gene function during neural development. Genome wide analysis of differentially expressed transcripts after such an experimental manipulation has the potential to uncover an almost complete picture of the downstream effects caused by the transfected construct. This work describes a simple method for comparing transcriptomes from samples of transfected embryonic spinal cords comprising all steps between electroporation and identification of differentially expressed transcripts. The first stage consists of guidelines for electroporation and instructions for dissection of transfected spinal cord halves from HH23 embryos in ribonuclease-free environment and extraction of high-quality RNA samples suitable for transcriptome sequencing. The next stage is that of bioinformatic analysis with general guidelines for filtering and comparison of RNA-Seq datasets in the Galaxy public server, which eliminates the need of a local computational structure for small to medium scale experiments. The representative results show that the dissection methods generate high quality RNA samples and that the transcriptomes obtained from two control samples are essentially the same, an important requirement for detection of differential expression genes in experimental samples. Furthermore, one example is provided where experimental overexpression of a DNA construct can be visually verified after comparison with control samples. The application of this method may be a powerful tool to facilitate new discoveries on the function of neural factors involved in spinal cord early development.
Genetic studies on live organisms frequently use the chicken embryo as a model because in ovo electroporation represents a fast and inexpensive way to partially transfect embryonic structures in vivo with DNA constructs1-5. Bicistronic expression vectors that encode one transcript containing a fluorescent reporter together with the gene of interest, such as pCIG6 or pMES7, allow quick verification of transfection quality in the desired region under a stereomicroscope before processing embryos for downstream analysis.
With recent advances in DNA sequencing, it is now possible to obtain digital whole transcriptome expression profiles from RNA samples using RNA-Seq8,9. Therefore, instead of time-consuming methods that enable analysis of only a few gene products in parallel, such as in situ hybridization, immunohistochemistry and qPCR, preparation of cDNA libraries for high-throughput sequencing can provide information of the whole transcriptome. Observation of the experimental effects on the expression levels of every single gene may in turn provide powerful insights on the pathways altered by the construct used. This is particularly easy if a genomic assembly for the organism used is available, thus the chick embryo again is a suitable model.
If the user has access to a reliable genome sequencing center, the main issue is obtaining high quality RNA samples. This work demonstrates a method for obtaining RNA samples from transfected neural tubes suitable for RNA-Seq, as well as the downstream steps for in silico analysis of the resulting datasets. After electroporation at HH12-13 followed by 48 hr incubation, transfected trunk spinal cord halves were dissected in ribonuclease-free conditions, immediately lysed and further protected from degradation by ultra-low freezing until RNA purification and library preparation.
The Galaxy public server10 provides access to free resources for analyzing high-throughput sequencing data. The amount of space offered for registered users is enough for small experiments, thus eliminating the need for local computational infrastructure. RNA-Seq data analysis includes filtering, alignment and quantitation/comparison of gene expression. Although there are general guidelines for RNA-Seq data analysis, filtering parameters will largely depend on the quality of sequencing and should be determined after quality control analysis of the raw data.
This protocol describes the steps to obtain and compare RNA-Seq profiles from chicken embryonic spinal cords transfected in vivo. As a representative result, comparison of datasets obtained from two independent control samples transfected with an empty vector showed that the method is reproducible and should allow identification of differentially expressed transcripts in experimental samples. Therefore, application of this method has the potential to greatly increase the number of discoveries after one single experiment and thus provide a large set of data for subsequent investigation.
1. Electroporate the Construct of Interest in the Neural Tubes of HH12-13 Chicken Embryos In ovo
Use a fluorescent reporter, preferably in a bicistronic transcript or fused to the protein of interest with an intercalating viral 2a peptide11, to enable quick identification of individuals with satisfactory levels of transfection. NOTE: Typical incubation time for this stage is around 48 hr at 37.8 °C.
2. Harvest Successfully Transfected Embryos and Dissect the Electroporated Spinal Cord Halves
NOTE: This procedure lasts from 2-3 hr for each 10 embryos and yields 1-2 µg total RNA per embryo. The dissection of neural tubes requires fine movements and some practice sessions may be necessary before applying it to experimental samples.
3. Purify RNA and Prepare Libraries for Illumina Sequencing
4. Compare RNA-Seq Datasets in the Galaxy Public Server10.
NOTE: Server storage space is limited and experiments comparing many samples may require deletion of intermediate files generated in the filtering steps.
To validate the method described, two control samples were generated from embryos electroporated with the empty vector pMES7 and one sample from embryos transfected with the same vector containing the insert construct SCRT2-ZNF, which encodes the zinc-finger domains of chicken SCRT218. Samples yielding 11-22 µg total RNA each were obtained from pools of 8 to 12 embryos. All three samples scored an optimal RIN value of 10 in a quality check performed with the Bioanalyzer instrument (Figure 1), demonstrating that high quality RNA can be obtained with this protocol.
After sequencing, the described filtering process increased the alignment rate from 83 to 91% of reads (data not shown). As expected, comparison of the transcriptomes from two control samples showed no clear difference when globally compared in a scatter plot (Figure 2) and the Pearson correlation coefficient was 0.99. Additionally, display of coverage for all the three datasets along the locus that contains the portion of the SCRT2 sequence showed a clear increase in the number of reads mapped to this region in the experimental sample (Figure 3), indicating that the overexpression was successful and that differences between control and experimental samples should be detectable with the method described here.
Figure 1: RNA quality control. Bioanalyzer electropherograms showing high level of integrity for the two control total RNA samples indicated by the 18S and 28S ribosomal RNA peaks. Insets represent a virtual gel for the run.
Figure 2: Reproducibility of RNA-Seq quantitation results. Scatter plot showing that RPKM for most genes is the same in both control samples. Virtually all genes that do not follow this pattern have very low RPKM values and thus are usually ignored from transcriptome comparisons.
Figure 3: Visualization of differences in an experimental sample. UCSC Genome Browser display of the SCRT2 locus in the chicken genome (galGal4) showing a visible increase in coverage of the region that codes for the protein domain (SCRT2-ZnF) ectopically expressed in the experimental sample.
Here we provide guidelines for analyzing effects after electroporation of the chicken spinal cord. Although electroporation of DNA vectors is more frequently used to overexpress a gene of interest, one can also use constructs encoding dominant negatives, quimeric proteins or precursors for siRNAs to generate gene function knock-down conditions19-21. In fact, comparison of transcriptome profiles resulting from both gain- and loss-of-function analysis may point out genes differentially expressed in opposite ways in each condition, thus having the potential to uncover more specific effects.
As any method applied downstream to RNA extraction, synthesis of RNA-Seq libraries requires a minimum initial amount of starting material. With this protocol, 8 embryos should be sufficient to yield around 10 µg of total RNA. Starting with a lower number of individuals without compromising an adequate working concentration is probably possible with the use of RNA extraction systems that allow final dilution in smaller volumes.
Proper segregation of the neural tube from all mesoderm should be practiced a few times before applying to important samples, as this is a delicate step that requires some familiarity. To avoid excessive RNA degradation during the dissection step, it is also important to complete the whole procedure for all neural tubes in about 1-2 hr. Again, practice will be important to improve dissection time. One important observation is that the connection between mesoderm tissues is stronger than the connection between mesoderm and neural tube, thus pulling mesodermal structures, specially the notochord, is more efficient than trying to cut unwanted tissue. Also, researchers should feel free to develop their own technique optimized to their own dissection conditions. However, care should be taken for samples compared in a same experiment to be dissected similarly so that level of contamination with mesodermal tissues is as homogeneous as possible between samples.
Our protocol demonstrates the comparison of transcriptome profiles from stage HH23 chicken embryonic spinal cords transfected at stage HH12-13. However, small adjustments should allow its application to embryos harvested at earlier stages. On the other hand, if later stages are desired for analysis of the end result, the bicistronic vectors mentioned may not be suitable due to a continuous reduction in the number of molecules per cell as development proceeds. Therefore, especially if the intended end point is achieved after more than 72 hr after transfection, it may be necessary to use vectors encoding transposon mediated genomic integration of the construct being used22.
Another issue to be considered is the heterogeneity of the electroporation efficiency. The number of non-transfected cells ranges from a few to a considerable amount, especially the ones in the ventral neural tube. Since after dissection all cells are subjected to downstream analysis, a low electroporation efficiency will reduce the assay’s sensitivity. An approach to enrich the sample with electroporated cells is the use of fluorescence activated cell sorting (FACS) to reduce the population of non-transfected cells from the final sample23. If such an approach is used, however, a greater number of embryos will be necessary to achieve sufficient yield for RNA-Seq library preparation.
Illumina sequencing is recommended as this is more widely offered by sequencing centers and in silico RNA-Seq data analysis is simplified when compared to other available platforms24. However, other options are still adequate with modifications in data filtering steps more suitable to the platform chosen. RNA-Seq data analysis guidelines are still not completely established and comparison of outputs from different available tools, some of which are already included in the Galaxy public server, is highly recommended. Nevertheless, the in silico analysis described here should be suitable for obtaining good comparison results from high quality datasets. As in all other areas in Biology, high-throughput sequencing increases enormously the amount of information that can be obtained from a single experiment, and the method described here is one way of applying it to the study of spinal cord development.
The authors have nothing to disclose.
CYIY is supported by FAPESP (2012/14421-5) and FMV is supported by a fellowship from FAPESP (2009/53695-0). We thank the DNA Technologies Core at University of California, Davis for preparation and sequencing of the RNA-Seq libraries used in this work, the Galaxy team for providing an excellent and free interface for high-throughput sequencing data analysis tools and Dr. Marianne Bronner for allowing us to film in her laboratory space.
Name of Material/ Equipment | Company | Catalog Number | Comments/Description |
Indian ìnk | Any supplier | ||
18G x 1 1/2 needle | BD | 305196 | |
3 ml syringe | BD | 309657 | |
Tris | MP Biomedicals | 819620 | |
F D & C Blue 1 | Spectra Colors Corporation | 5.FC.0010P0 | Food Dye used to visualize DNA injection |
Ringer's solution (1 L) | |||
– 7.2 g NaCl | Sigma-Aldrich | S-9888 | |
– 0.37 g KCl | Sigma-Aldrich | P-4504 | |
– 0.225 g CaCl2.2H2O | Fisher Scientific | 10043-52-4 | |
– 0.217 g Na2HPO4.7H2O | Sigma-Aldrich | S-9390 | |
– 0.02 g KH2PO4 | Sigma-Aldrich | P-5379 | |
– ddH2O to 800 ml and adjust pH to 7.4 | |||
– ddH2O to 1 L | |||
– Filter and autoclave | |||
PBS (Phosphate Buffered Solution) (1 L) | |||
– 8 g NaCl | Sigma-Aldrich | S-9888 | |
– 0.2 g KCl | Sigma-Aldrich | P-4504 | |
– 1.15 g Na2HPO4.7H2O | Sigma-Aldrich | S-9390 | |
– 0.2 g KH2PO4 | Sigma-Aldrich | P-5379 | |
– ddH2O to 800 ml and adjust pH to 7.2 | |||
– ddH2O to 1 L and filter | |||
– 1 ml DEPC (Diethylpyrocarbonate) | Sigma-Aldrich | D-5758 | |
– Shake vigorously and let stand overnight at room temperature | |||
– Autoclave | |||
Sieved spoon | Any supplier | ||
Sterile 60 x 15 mm polystirene Petri dishes | Corning Life Sciences | 351007 | |
RNaseZap RNase decontamination solution | Life Technologies | AM9780 | |
Fine point surgical scissors | Any supplier | ||
Straight fine point tweezers | Any supplier | ||
Pulled glass needle made from 1.1 x 75 mm glass capillary tubes | Kimble Chase | 40A502 | |
9'' glass Pasteur pipette | Any supplier | ||
Manual pipette pump | Any supplier | ||
Clear 1.5 ml microcentrifuge polypropilene tubes | Corning Life Sciences | MCT-150-C | |
Microcentrifuge | Any supplier | ||
RNAlater solution for RNA stabilization | Life Technologies | AM7020 | |
RNAqueous total RNA isolation kit | Life Technologies | AM1912 | |
Molecular Biology Grade Ethanol | Any supplier | ||
RNA 6000 Nano Kit | Agilent Technologies | 5067-1511 | |
2100 Electrophoresis Bioanalyzer Instrument | Agilent Technologies | G2939AA | |
Truseq RNA Sample Preparation Kit v2 | Illumina | RS-122-2001/RS-122-2002 |