Experimental Design for Laser Microdissection RNA-Seq: Lessons from an Analysis of Maize Leaf Development

Podsumowanie

Many developmentally important genes have cell- or tissue-specific expression patterns. This paper describes LM RNA-seq experiments to identify genes that are differentially expressed at the maize leaf blade-sheath boundary and in lg1-R mutants compared to wild-type. The experimental considerations discussed here apply to transcriptomic analyses of other developmental phenomena.

Streszczenie

Genes with important roles in development frequently have spatially and/or temporally restricted expression patterns. Often these gene transcripts are not detected or are not identified as differentially expressed (DE) in transcriptomic analyses of whole plant organs. Laser Microdissection RNA-Seq (LM RNA-Seq) is a powerful tool to identify genes that are DE in specific developmental domains. However, the choice of cellular domains to microdissect and compare, and the accuracy of the microdissections are crucial to the success of the experiments. Here, two examples illustrate design considerations for transcriptomics experiments; a LM RNA-seq analysis to identify genes that are DE along the maize leaf proximal-distal axis, and a second experiment to identify genes that are DE in liguleless1-R (lg1-R) mutants compared to wild-type. Key elements that contributed to the success of these experiments were detailed histological and in situ hybridization analyses of the region to be analyzed, selection of leaf primordia at equivalent developmental stages, the use of morphological landmarks to select regions for microdissection, and microdissection of precisely measured domains. This paper provides a detailed protocol for the analysis of developmental domains by LM RNA-Seq. The data presented here illustrate how the region selected for microdissection will affect the results obtained.

Wprowadzenie

The maize leaf is an ideal model to study the formation of developmental fields during morphogenesis, as it has a distinct boundary between the blade and sheath that is amenable to genetic dissection (Figure 1A). During the early stages of leaf development, a linear band of smaller cells, the preligule band (PLB), subdivides the leaf primordium into pre-blade and pre-sheath domains. A fringe-like ligule and triangular auricles develop from the PLB (Figure 1A, C, D). Genetic screens have identified mutations that disrupt the blade-sheath boundary. For example, recessive liguleless1 (lg1) mutations delete the ligule and auricles^1,2,3,4 (Figure 1B). In situ hybridization revealed that lg1 transcript accumulates at the PLB and emerging ligule, making it an excellent marker for ligule development^5,6 (Figure 1E).

Figure 1: Wild-type and liguleless1-R maize leaves. (A) Blade-sheath boundary region of mature wild-type leaf showing ligule and auricle structures. (B) Blade-sheath boundary region of mature liguleless1-R leaf showing absence of ligule and auricle structures. Leaves in A and B have been cut in half along the midrib. (C) Longitudinal section through wild-type leaf primordium. Sample has been processed and stained for histological analysis. The initiating ligule is apparent as a bump protruding from the plane of the leaf (arrowhead). (D) Longitudinal section through wild-type leaf primordium. Sample has been processed for LM as described in the text. Arrowhead indicates initiating ligule. (E) lg1 in situ hybridization of shoot apex lateral longitudinal section. Asterisks indicate lg1 transcript accumulation at the PLB of the P6 leaf primordium. Arrows indicate base of P6 primordium. Bar indicates measurement from the base of the primordium to the PLB. Scale bars in A and B = 20 mm. Scale bars in C-E = 100 µm. This figure has been modified from reference⁶(Copyright American Society of Plant Biologists). Please click here to view a larger version of this figure.

In this study, LM RNA-Seq was employed to identify a suite of genes that are differentially expressed (DE) at the blade-sheath boundary relative to other parts of the leaf primordium and to identify genes that are DE in lg1-R mutants relative to wild-type siblings. LM RNA-Seq is a method of quantifying transcript accumulation in specific cells or cellular domains⁷. LM systems combine a laser and a microscope with a digital camera. Sectioned tissue is mounted on slides and viewed through the microscope. The LM software typically includes drawing tools that allow the user to outline any selected region for microdissection. The laser cuts along the line, and the selected tissue is catapulted off the slide and into a tube suspended above the slide. LM allows the user to microdissect precise domains, including specific cell layers and even single cells^8,9. RNA can then be extracted from the microdissected tissue. Subsequently, the RNA-Seq component utilizes next-generation sequencing to sequence cDNA libraries generated from the extracted RNA^10,11.

Key advantages of LM RNA-seq are the ability to quantify transcript accumulation in precisely defined domains and the capacity to profile the entire transcriptome simultaneously⁷. The technique is particularly suited to probing early developmental events where the region of interest is often microscopic. Previous studies have utilized LM combined with microarray technology to study developmental processes in plants^9,12,13. RNA-Seq has the advantage of quantifying transcripts across a broad dynamic range, including low-expressed genes, and prior sequence information is not required^10,11. Moreover, LM RNA-Seq has the potential to highlight developmentally important genes that may be missed in mutagenesis screens due to genetic redundancy or to lethality of the loss-of-function mutant.

Developmentally important genes, such as narrow sheath1 (ns1) and cup-shaped cotyledon2 (cuc2), often have specific expression patterns of just one or a few cells^17,18,19,20. Many are expressed only during early developmental stages and not in the mature organ. When whole organs or large domains are analyzed, these cell-specific transcripts are diluted and may not be detected in more conventional analyses. By permitting analyses of precisely defined domains, LM RNA-Seq enables these tissue-specific genes to be identified and quantified.

Crucial factors in the success of the experiments described here were a thorough histological analysis that guided selection of the appropriate developmental stage and domain for analysis, and precise measurement of cell-tissue domains for LM. To ensure that equivalent domains were sampled for all replicates, tissue was collected from leaf primordia at the same developmental stage and the microdissected domains were measured relative to morphological landmarks such as the emerging ligule (Figure 2). It is known that some genes are expressed in a gradient from the tip to the base of the leaf. By measuring precise domains, variation due to sampling from different locations along the leaf proximal-distal axis was kept to a minimum (Figure 3A). By microdissecting domains of the same size, variation due to differential dilution of cell-specific transcripts was also reduced (Figure 3B). Lateral longitudinal sections of the shoot apex were used for all microdissections. These are sections that are perpendicular to the midrib-margin axis (Figure 4). Using only sections that include the SAM ensures that equivalent lateral regions of leaf primordia are analyzed.

In samples processed and sectioned for LM, the first morphological sign of ligule outgrowth is a bump on the adaxial side due to periclinal cell divisions in the adaxial epidermis (Figure 1D, Figure 2). It was determined that the emerging ligule could be reliably identified at plastochron 7 stage leaf primordia. We were interested in genes expressed in the entire ligule region, including the emerging ligule and the cells immediately distal that will form the auricle. In order to ensure that equivalent tissue selections were made, the ligule bump was used as a morphological landmark and a 100 µm rectangle centered on the ligule bump was selected for LM (Figure 2A, 2B). Equivalent sized rectangles of pre-blade and pre-sheath were selected from the same leaf primordia.

Analyses of liguleless mutant plants presented a different challenge; lg1-R mutants do not form a ligule, therefore this morphological feature could not be used to select the region for LM. Instead, the domain of lg1 transcript accumulation in wild-type leaf primordia was determined, and a region that would encompass this domain was defined. These preliminary analyses were performed on seedlings from the same planting as were used for the final analysis, since previous work has shown that the location of the PLB varies depending on growth conditions. In situ hybridization indicated that lg1 transcripts accumulate in the PLB of P6 leaf primordia (Figure 1E). We selected a domain 400-900 µm from the base of the leaf primordia that encompassed the domain of lg1 expression (purple rectangles, Figure 2A) and captured these equivalent regions from wild-type and lg1-R plants. To minimize variation in genetic background and growth conditions when comparing transcript accumulation in lg1-R and wild-type plants, segregating families of mutants and wild-type siblings were used.

Protokół

NOTE: Fix tissue for histological analysis at the same time that tissue is fixed for LM. Examine stained sections for morphological features that will guide later LM. When comparing mutant to wild-type, perform in situ hybridization or immunolocalization to define the domain where the gene of interest is expressed (in this case lg1).

1. Tissue Fixation and Processing

1. Grow flats of maize seedlings to two weeks old under standard conditions⁶.
2. Dissect shoot apices for lateral sections (Figure 4).
   1. Excise seedling just below the soil line.
   2. Using a razor blade, remove thin slices from the base of the stem (cuts 1, Figure 4A) until an oval of culm encircled by one or two mature leaves is visible (Figure 4B).
   3. Make another cut approximately 10 mm above the base (cut 2, Figure 4A). This 10 mm segment will contain the SAM and young leaf primordia.
   4. Turn the 10 mm segment so the base is facing up. Make two cuts parallel to the lateral axis so that a slice of tissue 2-3 mm thick is obtained (cuts 3, and 4, Figure 4B). Discard the outer two portions and retain the central slice for fixation and embedding.
      NOTE: Outer leaves may be trimmed and discarded.
3. Fix tissue and process for embedding.
   1. Ensure that all materials to be used in subsequent steps are RNase free. Treat solutions with diethyl pyrocarbonate (DEPC) (1 ml DEPC per liter of solution. Incubate overnight with occasional shaking, and autoclave). Bake glassware in an oven at 200 °C or higher for at least 6 hr and treat plastic ware with RNase decontamination solution.
   2. Day 1: Immerse tissue slices in ~10 ml of Farmer's fix (3:1 Ethanol: Acetic acid) in glass vial on ice. After all samples have been dissected, apply vacuum to remove air bubbles and aid penetration of fixative. Hold under vacuum for 10 min then release the vacuum slowly. Replace fixative and incubate at 4 °C overnight with gentle shaking.
   3. Day 2: Incubate in the following series of solutions, ~10 ml each, 1 hr each, all with gentle shaking; 85% ethanol at 4 °C, 95% ethanol at 4 °C, 100% ethanol at 4 °C, 100% ethanol at 4 °C, 100% ethanol at 4 °C, 1:1 ethanol: xylenes at room temperature, 100% xylenes at room temperature, 100% xylenes at room temperature.
      NOTE: Xylenes are toxic by contact and inhalation. Work in a fume hood and use appropriate gloves.
   4. Add paraffin tissue embedding medium pellets to approximately half volume of xylenes and incubate overnight at room temperature with gentle shaking.
   5. Day 3: Transfer the vial to 60 °C oven until the pellets melt. Pour off solution and replace with fresh melted tissue embedding medium. Change the medium two more times during the day.
   6. Day 4: Change tissue embedding medium once in the morning. Return to 60 °C oven until afternoon.
4. Cast blocks
   1. Place embedding molds on hot plate of tissue embedding station. Use forceps to transfer the tissue samples to the embedding molds with the cut surface facing down. Top up the mold with melted paraffin and place the embedding ring on top of the mold. Transfer to a cold plate until paraffin has solidified. Store the paraffin blocks at 4 °C in an airtight container with silica gel.

2. Sectioning and Slide Preparation

1. Cut 10 µm sections on a microtome²⁵.
2. Examine ribbons and choose median sections. Median sections are those that include the SAM, which appears as a dome of cells surrounded by leaf primordia.
3. Mount sections on slides.
   1. Place slides that are suitable for LM (either RNase free or baked) on 42 °C slide warmer and apply several drops of 50% ethanol solution to cover the slide.
   2. Float sections on ethanol solution until the sections have expanded.
      NOTE: Floating sections on ethanol solution rather than water keeps RNA in a precipitated state reducing RNA degradation.
   3. Tilt slide and remove excess ethanol solution by aspiration with a disposable transfer pipette. Use lint-free wipes to wick away any additional ethanol solution.
   4. Dry slides at 42 °C for several hours or overnight. Store slides at 4 °C in an airtight container with silica gel.
4. Deparaffinize slides on the day of use.
   1. Prepare three glass Coplin jars containing; 100% xylenes (xylenes I), 100% xylenes (xylenes II), and 100% ethanol (~50 ml of each solution).
   2. Using clean forceps to transfer slides, immerse slides in xylenes I for 2 min, xylenes II for 2 min, and 100% ethanol for 1 min.
   3. Drain slides on lint-free wipes and air dry at room temperature.

3. Microdissection of Blade, Ligule and Sheath Samples from Plastochron 7 Leaf Primordia

1. Secure the slides on stage of LM microscope. Use five or six slides for each replicate, utilizing five sections per slide.
   NOTE: Tissue pooling for a single replicate is illustrated in Figure 5.
2. Examine slides and identify the five most median sections on each slide, using the SAM apex as the central reference point.
   NOTE: This can be done at low magnification, usually a 5X objective is sufficient.
3. Using 10X or 20X objective, identify the position of the ligule on the plastochron 7 leaf primordium of each section. The ligule will be visible as a bump protruding from the adaxial surface of the leaf primordium. Mark this position using the drawing tool of the LM software; select the pencil icon, move the cursor to the appropriate position and click and drag the mouse to draw.
   NOTE: A 10X or 20X objective is appropriate for this and subsequent steps. When using lateral sections, the two sides of each leaf primordium will be present in each section (Figure 2A).
4. Using the ruler tool and the rectangular drawing tool, measure 100 µm high rectangles centered on the ligule of each section (red rectangles, Figure 2A, 2B). These will be the "Ligule" sample.
   1. To use the ruler tool; select the ruler icon, move the cursor to one end of the object to be measured, click and drag to measure the object. The length of the ruler will be shown on the screen.
   2. To draw a rectangle; select the rectangle icon, move the cursor to a point that will be one corner of the rectangle, click and drag to draw a rectangle of the appropriate size. Alternatively, select the straight line drawing tool and draw four straight lines.
   3. Measure 100 µm rectangles positioned 50 µm above and below the "Ligule" rectangle.
      NOTE: These will be the "Blade" and "Sheath" samples, respectively (green and blue rectangles, Figure 2A, 2B). Based on our histological data, a 100 µm rectangle encompasses the entire ligule region. Equivalent sized portions of blade and sheath were chosen to ensure that similar amounts of tissue were collected for each. Spacers of 50 µm were used to ensure that no ligule region tissue is inadvertently included in the blade or sheath microdissections.
5. Microdissect measured rectangles (Figure 2D-2F)^7,8,9, collecting Ligule, Blade and Sheath samples in separate tubes. Use the laser cut function to cut through the tissue section along the outline of the selected domain. Use the catapult function to propel the rectangle of tissue off the slide and into the lid of the tube (Figure 2D-2F).

4. Microdissection of Blade, Ligule and Sheath Adaxial Epidermal Samples from Plastochron 7 Leaf Primordia

1. Select sections and use the ruler tool to measure 100 µm high segments centered on the plastochron 7 ligule, as described in section 3 (above).
   1. Select only the adaxial epidermal cells of each 100 µm high "Blade" and "Sheath" segment (green and blue selections, Figure 2C) by outlining with the drawing tool. The epidermis is the outer cell layer; the adaxial side is the one closest to the SAM.
   2. For the "Ligule" sample, select only the cells of the emerging ligule bump as described in Section 3.3 (red selection, Figure 2C).
2. Microdissect selected regions, collecting Blade, Ligule and Sheath epidermis samples in separate tubes, as described in section 3.5.

5. Microdissection of Plastochron 6 Leaf Primordia from lg1-R and Wild-type Siblings

1. Grow segregating families of mutant (lg1-R) and wild-type plants.
2. Fix and process shoot apices for LM, as described in sections 1.2-1.4. Fix wild-type and mutant shoot apices in separate vials in order to keep them separate. Samples from the same planting should be fixed and processed for in situ hybridization.
3. Determine where lg1 is transcribed in wild-type siblings, by performing lg1 in situ hybridization^6,26,27. Measure position of lg1 transcript accumulation from base of leaf primordium in multiple samples (Figure 1E).
4. Based on in situ hybridization data, choose portion of leaf primordium that encompasses the region where lg1 is transcribed. In this case, 400-900 µm from the base of plastochron 6 leaf primordia (purple rectangles, Figure 2A).
5. Microdissect selected portion of leaf primordia, as described in section 3.5, collecting lg1-R and wild-type samples in separate tubes.

6. Apply RNA Extraction Buffer

1. Apply 50 µl RNA extraction buffer to microdissected tissue and proceed with RNA extraction. Continue with RNA extraction, RNA amplification, library construction, sequencing and bioinformatics analysis as described in reference⁶.

Wyniki

Using the LM scheme outlined in Figure 2, approximately 1,000,000-1,500,000 µm² of tissue was collected for each replicate in the all-cell-layers LM (Figure 5), and 200,000 µm² per replicate for the adaxial epidermis LM. Approximately 2,500,000 µm² of tissue was collected for each replicate in the LM of lg1-R and wild-type leaf primordia. Two rounds of linear RNA amplification yielded microgram quan...

Dyskusje

Experimental design is a critical factor in RNA-seq experiments. Key considerations are the precise domain(s) and developmental stage(s) to be analyzed, and what comparisons will be made. It is crucial to think in terms of comparisons, since the output is typically a list of genes that are DE between two or more conditions. As with all experiments, it is important to alter only one variable at a time. For example, when comparing different leaf domains, leaves of the same age and developmental stage, grown under the same ...

Materiały

Name	Company	Catalog Number	Comments
Diethyl pyrocarbonate	Sigma-Aldrich	159220	Used for RNase treatment of solutions
Razor blades	Electron Microscopy Sciences	72000
RNase Zap	Sigma-Aldrich	R2020-250ML	RNase decontamination solution
Ethanol absolute 200 proof	Fisher Scientific	BP28184
Acetic acid, glacial	Sigma-Aldrich	A6283
Glass vials - 22 ml	VWR	470206-384
Xylenes, histological grade	Sigma-Aldrich	534056-4L
Paraplast plus	Sigma-Aldrich	P3683-1KG
Disposable base molds, 15 mm x 15 mm x 5 mm	VWR	15154-072
Embedding rings	VWR	15154-303
Silica gel packets	Electron Microscopy Sciences	71206-01	Desiccant for storage of paraffin blocks
Oven	Fisher Scientific	15-103-0503	Oven must maintain temperature of 60 °C
Paraffin embedding station	Leica	EG1160
Microtome	Leica	RM2235
Slide warmer	Electron Microscopy Sciences	71317-10
Coplin jars	Electron Microscopy Sciences	70316-02
Laser microdissector	Zeiss
KIMWIPES™ Delicate Task Wipers	Kimberly-Clark Professional	34120	Lint-free wipes for wicking excess solutions from microscope slides
Membrane Slide 1.0 PEN	Zeiss	415190-9041-000	Slides for laser microdissection
Adhesive Cap 200 opaque	Zeiss	415190-9181-000	Tubes for laser microdissection
PicoPure RNA Isolation Kit	ThermoFisher Scientific	KIT0204