# Combining laser-assisted microdissection (LAM) and RNA-seq allows to perform a comprehensive transcriptomic analysis of epidermal cells of Arabidopsis embryo

Kaori Sakai, Ludivine Taconnat, Néro Borrega, Jennifer Yansouni, Véronique Brunaud, Christine Paysant-Le Roux, Etienne Delannoy, Marie-Laure Martin Magniette, Loic Lepiniec, Jean Denis Faure, et al.

► To cite this version:

Kaori Sakai, Ludivine Taconnat, Néro Borrega, Jennifer Yansouni, Véronique Brunaud, et al. Combining laser-assisted microdissection (LAM) and RNA-seq allows to perform a comprehensive transcriptomic analysis of epidermal cells of Arabidopsis embryo. Plant Methods, 2018, 14, pp.10.10.1186/s13007-018-0275-x. http://hal-02625394

HAL Id: hal-02625394

https://hal.inrae.fr/hal-02625394

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Combining laser-assisted microdissection (LAM) and RNA-seq allows to perform a comprehensive transcriptomic analysis of epidermal cells of Arabidopsis embryo

Kaori Sakai<fnref n="1" />†, Ludivine Taconnat<fnref n="2" />,³†, Nero Borrega<fnref n="1" />, Jennifer Yansouni<fnref n="2" />,³, Véronique Brunaud<fnref n="2" />,³, Christine Paysant-Le Roux<fnref n="2" />,³, Etienne Delannoy<fnref n="2" />,³, Marie-Laure Martin Magniette<fnref n="2" />,³,⁴, Loïc Lepiniec<fnref n="1" />, Jean Denis Faure<fnref n="1" />, Sandrine Balzergue<fnref n="2" />,³,⁵* and Bertrand Dubreucq<fnref n="1" />* ID

**Abstract**

**Background:** Genome-wide characterization of tissue- or cell-specific gene expression is a recurrent bottleneck in biology. We have developed a sensitive approach based on ultra-low RNA sequencing coupled to laser assisted microdissection for analyzing different tissues of the small Arabidopsis embryo.

**Methods and results:** We first characterized the number of genes detected according to the quantity of tissue yield and total RNA extracted. Our results revealed that as low as 0.02 mm<fnref n="2" /> of tissue and 50 pg of total RNA can be used without compromising the number of genes detected. The optimised protocol was used to compare the epidermal versus mesophyll cell transcriptomes of cotyledons at the torpedo-shaped stage of embryo development. The approach was validated by the recovery of well-known epidermal genes such *AtML1* or *AtPDF2* and genes involved in flavonoid and cuticular waxes pathways. Moreover, the interest and sensitivity of this approach were highlighted by the characterization of several transcription factors preferentially expressed in epidermal cells.

**Conclusion:** This technical advance unlocks some current limitations of transcriptomic analyses and allows to investigate further and efficiently new biological questions for which only a very small amounts of cells need to be isolated. For instance, it paves the way to increasing the spatial accuracy of regulatory networks in developing small embryo of Arabidopsis or other plant tissues.

**Keywords:** Laser-assisted microdissection, Ultra-low-RNA-seq, Epidermis, Embryo cells, Arabidopsis

**Background**

Multicellular (higher) organisms like plants are characterized by cell-specific differentiation and tissue

formation during development. Over the last decades, many studies have addressed the question of gene expression during plant growth, under stress conditions, at precise stages or in specific genetic backgrounds. One aim of modern biology is to provide more quantitative and tissue- or cell-specific analyses with regard to the regulatory networks controlling biological processes. Nevertheless, the characterization and isolation of cell types organised in a three-dimensional space requires specific and tricky physical manipulations (for review see [13]). Different

<footnote>*Correspondence: sandrine.balzergue@inra.fr; bertrand.dubreucq@inra.fr</footnote>

<footnote>†Kaori Sakai and Ludivine Taconnat have contributed equally to this work</footnote>

<footnote n="1">¹Institut Jean-Pierre Bourgin (IJPB), INRA, AgroParisTech, CNRS, Université Paris-Saclay, RD10, 78026 Versailles Cedex, France</footnote>

<footnote n="2">²Institute of Plant Sciences Paris Saclay IPS2, CNRS, INRA, Université Paris-Sud, Université Evry, Université Paris-Saclay, Bâtiment 630, 91405 Orsay, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
approaches have been used including, the expression of reporter genes coupled to affinity purification system (e.g. INTACT, [10], fluorescent-activated cell sorting [2, 13, 35, 44], or manual dissection based on morphological characters). These techniques are limited in their application due to the need for protoplasting or genetic transformations with cell type-specific markers.

Laser Assisted Microdissection (LAM) allows precise recovery of specific tissues or cell types, based on their morphology or fluorescence staining when cell specific markers have been introduced. If based on paraffin embedded sections, the identification of specific zones can be facilitated by the use of a microscope, although decreasing the size of the section lowers the amount of RNA extracted. Similarly, very small areas of interest will have poor total RNA content, hampering comprehensive transcriptional analysis. The size of the zone of interest and its subsequent RNA content are therefore a crucial issue in the production of good-quality data. Coupling LAM to quantitative RT-PCR and then to DNA chips has increased the number of genes detected and thus, these approaches have been widely and successfully used on seed [4, 5, 37]. However, the development of NGS (Next Generation Sequencing) technologies offers many other advantages: sensitivity, ability to quantify expression in species for which no genome sequences are available (i.e. new species of interest), access to differentially spliced forms or to non-coding RNAs. Moreover, in a specific case of plant-pathogen interaction studies, the LAM provides also a cellular-level resolution to reveal the often low coverage of pathogen transcripts [14]. The NGS has an additional advantage compared to traditional chip hybridization, as it requires a smaller amount of total RNA for analysis. Some studies combined LAM, RNA ribosomal depletion and RNA-seq experiment in order to reveal the entire diversity of transcripts [16, 30]. In these cases, the amounts of total RNA coming from LAM is in the range of ng, requiring spending a lot of time to microdissection step, often at the expense of RNA quality or/and biological repetition number. Therefore, we combined ultra-low-RNA-seq sequencing with laser microdissection for undertaking precise and comprehensive analyses of the epidermal versus mesophyll cell transcriptions of the cotyledons, at the torpedo-shaped stage of the small Arabidopsis embryo.

The development of Arabidopsis embryo has been extensively described [17, 41]. Briefly, after fertilization, the zygotic cell gains polarity and develops following a precise pattern of divisions to give rise to specific cell types (epidermis, vascular bundles, cortex, and mesophyll). The epidermis, which originates from the differentiation of the protoderm at the dermatogen stage, marks the junction between the embryo and the external

environment [15]. This cell layer forms a hydrophobic barrier over the aerial surfaces of the plant [39]. Specific secondary metabolites accumulate in the epidermis in the form of cuticular fatty acids. The cutin and flavonoids then form protective compounds. Previous analyses have shown that several genes are specifically expressed in the epidermis during aerial organ development [39] or during embryo development [15, 27, 40].

In this paper, we investigated the effect of decreasing the total amount of tissue and RNA template using ultra-low RNA-seq methods and the technical limit of this approach. After validation, the optimized method was used to characterize genes that are differentially expressed in the epidermis versus mesophyll cells at the torpedo stage of Arabidopsis embryo development.

## Methods

### Plant material and growth conditions

*Arabidopsis thaliana* plants, accession Columbia (Col-0), were grown in a greenhouse under the following conditions: 13 h of light, 25 °C/17 °C day/night, and irrigated three times per week with mineral nutrient solution. To harvest seeds at defined developmental stages, individual flowers were tagged on the day of opening, and then opened flowers and developing siliques were counted daily. Siliques at 8 days after fertilization, corresponding to seeds containing embryos at linear stage were harvested under RNase free conditions: all materials and working surfaces were treated with RNase Zap (Ambion) and immediately fixed in 3:1 (vol/vol) ethanol:acetic acid at 4 °C on ice. Siliques were cut at the edge into ~ 1 cm segments before fixation to allow better penetration of the fixator. The seeds were fixed under vacuum for 1 h and left O/N in the fixator at 4 °C.

The plant material was dehydrated in a graded ethanol series (70% 1 h, 85% 1 h, 95% 1 h, 100% 1 h two times, 100% ethanol O/N), and infiltrated with histoclear (1:3 1 h, 1:1 1 h, 3:1 2h30 histoclear: ethanol). This was followed by 100% histoclear for 30 min three times. Samples were then incubated with 1:1 paraffin/histoclear for 1 h and paraffin 100% at 60 °C O/N. The paraffin was replaced twice over 1–2 days. Seeds were sectioned at 8 µm thickness using an automatic microtome (Microm HM 355S) and mounted on polyethylene napthalate (PEN)-membrane slides (Zeiss) in RNase-free conditions. Slides were dried with a hot plate set at 24 °C and deparaffinized twice in 100% histoclear for 1 min and dehydrated in 100% ethanol for 1 min. Laser capture microdissection was performed using a PALM Micro-Beam system (Zeiss). For pilot experiment, 20 whole cotyledons sections coming from 6 siliques on 3 individual plants were microdissected and captured. For the second experiment, each tissue type (i.e. Epidermis and
Mesophyll) of each seed was separately microdissected to minimize contamination from adjacent cell and tissue types (Fig. 2a). Four biological replicates harvested at four different dates were captured for each tissue type. Each biological replicate consisted of around 25 microdissected tissue sections from at least 5 siliques coming from one individual plant. All tissues were captured in a collection tube with adhesive cap (Zeiss) within 10 min to maximize the quality of total RNA extraction.

## RNA-seq experiments

Microdissected samples were harvested and incubated directly into RNA extraction buffer. Total RNA was extracted using the Arcturus Pico RNA extraction kit (Thermo Fisher Scientific, Inc) and then treated with RNase-free DNase (1:8 dilution of DNase I in RDD buffer; Qiagen). RNA quantity and quality were checked by microcapillary electrophoresis RNA 6000 Pico Chip (Agilent 2100 BioAnalyzer; Agilent Technologies Waldbroon, Germany). RIN (RNA Integrity Number, Agilent) obtained were around 6.3 (6.0 for epiderm and 6.5 mesophyll: range between 5.7 and 6.7).

For the pilot experiment, in order to focus on optimizing the quality of extracted RNA in relation to the time spent on microdissection, 6 dilutions of a same RNA sample of whole cotyledons were used (5 ng, 100 pg, 75 pg, 50 pg, 25 pg, and 10 pg). For the second experiment, the objective of which was to compare epidermis and mesophyll tissues, a dilution of 100 pg of total RNA was used. cDNA syntheses were performed using the SMARTer Ultra Low Input RNA Kit for Sequencing-v4 (Clontech Laboratories, Inc.) and libraries were prepared according to DNA Sample Preparation Illumina kit instructions with a different bar code for each sample (Illumina, Cat. Nos. FC-131-1024). Ultra-low RNA-seq libraries were checked for their quality on High-sensitivity DNA chip using Agilent 2100 bioanalyzer (Waldbroon, Germany) before Illumina sequencing (Illumina®, California, U.S.A.). The UltraLowRNA-seq samples were sequenced in Paired-End (PE) with a read length of 100 bases. For the pilot experiment, the six libraries were sequenced on Hiseq 2000 machine. Samples with a dilution of 100 pg and 5 ng were first sequenced to obtain around 50 million reads/sample. Samples with a dilution of 100 pg, 75 pg, 50 pg, 25 pg, and 10 pg were then sequenced at a second date to obtain around 30 million reads/sample. The 100 pg sample was thus sequenced twice and used to correct a sequencing bias effect. The value of Spearman's correlation coefficient between the two 100 pg samples was 0.61 before correction and 0.86 after. After correction, the results of the two sequencing for the 100 pg sample were similar and only the second one was kept for the analyses. For the second experiment,

the eight libraries were sequenced on Hiseq 2000 machine to obtain around 35 million reads/sample with a multiplexing rate of 4 libraries/lane.

## Bioinformatics and statistical analyses

The raw data (fastq) were trimmed for Phred Quality Score > 20, read length > 30 bases, and the ribosome sequences were removed with tool sortMeRNA [19]. The mapper Bowtie version 2 [20] was used to align reads against the *A. thaliana* transcriptome (with 'local' option and other default parameters). The 33602 annotated genes were extracted from TAIR (v10) database corresponding to the representative gene model (longest CDS) given by TAIR. The abundance of each gene was calculated by a home-made script which counts only paired-end reads for which both reads map unambiguously one gene, removing multi-hits.

All the statistical analyses were done with the R software using also EdgeR package version 3.8.6 [28]. For the pilot experiment, to compare the gene expression, the raw counts were normalized to take the difference of the library sizes into account with TMM method and a sequencing date effect. It was done with a negative binomial generalized linear model with one factor (sequencing date). Normalized counts equal to the raw counts divided by the scaling factor minus a date sequencing effect. The second experiment concerned the tissue comparison between epidermis and mesophyll. First, genes which did not have at least 1 read after a count per million normalization in at least one half of the samples, were discarded. Then, raw counts were normalized using TMM method and count distribution was modelled with a negative binomial generalized linear model where the tissue type and the harvest date were taken into account and where the dispersion is estimated by the edgeR method. A likelihood ratio test was performed to evaluate a tissue effect. Raw p-values were adjusted with the Benjamini-Hochberg procedure to control the False Discovery Rate (FDR). A gene was declared differentially expressed if its adjusted *p* value ≤ 0.05.

## Data deposition

RNA-Seq projects were deposited in the international repository GEO (Gene Expression Omnibus, Edgard et al. [11]): http://www.ncbi.nlm.nih.gov/geo/; accession no. GSEGSE NGS2014_01_MicroD GSE98176 163 and NGS2012_02_MicroD GSE98313, according to the MINSEQE 'minimum information about a high-throughput sequencing experiment'. All steps of the experiment, from growth conditions to bioinformatics analyses, were detailed in CATdb [12]: http://tools.ips2.u-psud.fr/fr/CATdb/; Project: NGS2012_02_microD and NGS2014_01.
## Results

We analysed the quality and quantity of total RNA extracted from microdissected tissues performed on paraffin-embedded seeds according to the protocol previously described. In order to optimize RNA quality taking into consideration the microdissection process itself, we first evaluated the quantity and quality of total RNA extracted from the entire embryo. Focusing on the relationship between RNA quality and quantity extracted, we observed that RIN values are very variable, between 0 and 7.8, for very low concentrations of total RNA (from 0 to 200 pg/µl) (Fig. 1a). At a surface of 20.000 µm², the quality appears to stabilize between RIN 5.5 and 6 for higher concentration (Fig. 1c). Although, as expected, the RNA quantity increases with the amount of microdissected surface (Fig. 1b) whereas no improvement of the RNA quality was obtained by increasing the microdissected

surface (Fig. 1c). This surface of 20.000 µm² was thus chosen as a reference to reach around 200 pg of total RNA. High quality Epidermis/Mesophyll microdissected total RNAs (RIN = 6.8) coming from torpedo stage embryos were obtained and used for RNA-seq (Fig. 2b).

We first investigated the number of genes detected according to the quantity of mRNA template by using the pilot experiment (see Methods section). The numbers of reads and detected genes obtained for each sample after mapping and gene allocations are shown in Table 1. The quality of sequencing was very good for each library, with around 90% of reads mapped and unambiguously associated to genes. The number of detected genes decreases from 19266 to 15936 in parallel with RNA quantities ranging from 100 to 10 pg. The percentage of genes having at least 1 read in the 5 ng RNA sample and conserved in other samples is around 90% for 75 and 100 pg samples, and around 85% for the other. Normalized count distribution is presented in Fig. 3. Median of the gene expression level seems to be lower for samples with a dilution of 10 and 25 pg. Moreover, the first axis of the Principal Component Analysis (PCA) made from the normalized counts states a clear separation of 10 pg sample from the others and the second axis a clear separation of 25 pg sample from the others. Both axes explain more than 50% of the variability (Additional file 1: Fig. S1). This discrepancy between 10 and 25 pg samples and the others is also observed on the dendrogram cluster graph (Additional file 2: Fig. S2). A scatter-plot matrix showed a correlation between samples increasing with the RNA quantity from 10 pg to 5 ng, and always greater than 0.83 from 5 ng (Fig. 4).

We then performed a microdissection experiment using a 100 pg RNA template and 4 biological replicates, with the aim to identify differentially expressed genes of epidermis versus mesophyll cells. Here, we used cotyledons at the torpedo stage of embryo development (see Methods section). To assess data quality, we performed a principal component analysis (Additional file 3: Fig. S3). The latter clearly showed no replicate bias and the two first axes, which represent 41.46% of the explained variance discriminate the two types of tissue. It is to be noted that epidermis samples are much more variable than mesophyll samples that may reflect a contamination by endosperm-specific expressed genes during the dissection process.

The analysis of differential gene expression based on four biological replicates allowed to characterize 870 Differentially Expressed Genes (DEGs) with an adjusted *p* value < 0.05 (Table 2). We compared this number with respect to the number of biological replicates used. As expected, the results showed that the overall number of DEGs increased with the number of biological replicates.

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Micro-dissected samples in this study. **a** Mesophyll or epiderm of linear staged embryos were microdissected at X40. Image shows the different steps of microdissection process: area selection, laser cutting, catapulting and capture of the sample. In red the epidermis, in blue the mesophyll. Bar = 30 µm. **b** Quality of the total RNA (Agilent Bioanalyzer profile) extracted after microdissection of the epidermis. 1: marker, 2: small RNA, 3: 18S rRNA, 4: 28S rRNA</figcaption></figure>

The significant log ratio beyond which a gene is declared differentially expressed did not change: an average of 0.65 for 2 replicates and 0.75 with 3 replicates, whereas a significant log-ratio of 0.6 was obtained with the 4 replicates (Table 2). Thus, in this experiment, increasing the number of biological replicates did not reduce the detected differential expression between two conditions, but increased the completeness of the information.

To validate the approach and confirm the biological significance of the DEGs identified, we selected two different metabolic pathways known to take place more specifically in the epidermal cells, wax and flavonoid synthesis, and checked the expression of the genes involved in our data. RNA-Seq analysis revealed a specific subset of fatty acid biosynthetic genes, whose expression is higher in the epidermal cell layers than in mesophyll cells (Fig. 5). One of the main physiological functions of the epidermis is to produce a protective hydrophobic barrier made of epicuticular waxes and cutin. It involves the synthesis of alkanes, alcohols aldehydes, free fatty acids, and wax esters that are all derived from long and very long chain acyl-CoA, elongated in the Endoplasmic

<table><caption>**Table 1** Read mapping and gene detection statistics based on starting RNA quantity</caption><thead><tr><th>RNA template</th><th>5 ng</th><th>100 pg</th><th>75 pg</th><th>50 pg</th><th>25 pg</th><th>10 pg</th></tr></thead><tbody><tr><td>Librairy size (million)</td><td>53.341628</td><td>38.317268</td><td>25.174298</td><td>31.510988</td><td>32.682558</td><td>29.791112</td></tr><tr><td>Mapped reads (million) (percentage)</td><td>51.127950 (96%)</td><td>30.780261 (81%)</td><td>22.402608 (89%)</td><td>28.652941 (91%)</td><td>29.123427 (89%)</td><td>26.701774 (90%)</td></tr><tr><td>Detected genes (level 0)</td><td>17446</td><td>19266</td><td>19539</td><td>18485</td><td>16829</td><td>15936</td></tr><tr><td>Common reads with 5 ng reference sample (percentage)</td><td></td><td>15704 (90%)</td><td>15880 (98.8%)</td><td>15317 (87.8%)</td><td>14312 (85%)</td><td>13748 (86.3%)</td></tr></tbody></table>

Both first lines are the number of sequenced PE reads in million(s), and number of PE reads kept after mapping and gene association. Detected genes are the number of genes with a 1 PE or more. The last line is the number of genes in common between 17446 genes detected with 5 ng compared to the other samples. The raw counts are already rid of rRNA reads, 3–6% have been removed for each sample

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Boxplot of normalized counts after Log2 + 1 transformation, for the 6 samples of the pilot experiment (RNA quantity from 5 ng to 10 pg)</figcaption></figure>
<figure><img src="image_7.png" /></figure>

Reticulum [29]. The list of genes detected is indicated in the Fig. 5 and as well as those strongly expressed in the epidermis identified by a red dot. A limited subset of 27 genes was found differentially expressed in embryo epidermis. Interestingly, these genes were mostly found either upstream or downstream of the wax and cutin biosynthetic pathways. Thirteen genes were involved in

acylCoA synthesis (ACCase, LACS and elongase), and ten were associated with transport/export of wax and cutin constituents of the extracellular space.

We then checked the expression of known flavonoid genes. This pathway is usually presented with a common part made of the “early biosynthetic genes” (EBGs) that leads to flavonoid precursors (i.e. dihydroflavonols)
<table><caption>**Table 2** Number of differentially expressed genes detected based on number of biological replicates</caption><thead><tr><th rowspan="2"></th><th colspan="3">Number of repeats</th></tr><tr><th>2</th><th>3</th><th>4</th></tr></thead><tbody><tr><td>Log ratio for DE genes</td><td>0.65</td><td>0.75</td><td>0.65</td></tr><tr><td>Number of DE genes</td><td>752</td><td>803</td><td>870</td></tr></tbody></table>

Top—comparison of gene expression levels needed (log ratio) to validate a differential expression (DE) according to number of repeats with 0.1 ng total microdissected input RNA. Bottom—number of genes differentially expressed according to organs and repeats

and more specific branches made of the late biosynthetic genes (LBG) involved in the synthesis of specific flavonoids (e.g. anthocyanins or tannins) [21]. The analysis showed that most of the EBGs are DEGs in the epidermis compared to the mesophyll (Fig. 6) that is fully consistent with epidermal biosynthesis of flavonoids. Interestingly, flavonoids transcriptional regulators *ENHANCER OF GLABRA 3* (*EGL3*, At1g63650) and *MYB111* (At5g49330), involved in anthocyanins and flavonols biosynthesis, respectively [43], were also shown to be differentially accumulated in the epidermis.

To go further in the transcriptional regulations occurring specifically in the epidermis, we focused on the transcription factors (TF) that are differentially expressed with a log ratio > 1 in the epidermis when compared to the mesophyll cell (Table 3). From the 870 DEGs genes, 571 had a log ratio > 1 and from these 571, 43 (9.49%) were transcription factors, preferentially expressed in epidermal tissues, that is slightly higher than the normal occurrence of TF in Arabidopsis genome (5.6% - [9], and 33 have been found in embryo (http://bar.utoronto.ca/efp2/Arabidopsis/Arabidopsis_eFPBrowser2.html) and, according to literature, 9 more specifically in epidermis (Table 3). The log ratios of TF DEGs with a log ratio > 1, ranges from 1.08 for *FUS3* to > 10 for *MYB56* with 7 TF showing a log ratio above 5.

## Discussion

In this paper we have undertaken a microdissection followed by RNA-Seq in developing Arabidopsis embryos. We showed that robust and comprehensive RNA-Seq experiment can be performed starting with very small amount of starting material and optimized microdissection time. Embedded sections usually produce lower quality RNA than fresh tissue. However, the quality of the section remains much higher in embedded tissues compared to cryo-sections, which shrivel quickly at room temperature under the LAM. A key advantage of embedded sections is to reduce damages to organ structure before and during micro-dissection. Thus, the combination of LAM and RNA-seq provide new powerful

approach for investigating the transcriptome of specific tissues or cell types. First, based on mapped reads and detected genes, we concluded that a low amount of total RNA (i.e. 10 pg), is sufficient to generate libraries and produce RNA-seq data. The number of genes detected was comprised between 15936 and 17446, which is fully comparable with embryo data obtained in previous studies based on microarrays technologies [7, 37]. Surprisingly, the number of genes detected in the 5 ng sample is lower than expected when compared to other samples (i.e. 100–50 pg), nevertheless we can observe a decrease of the number of genes detected according to the RNA quantity. The overlap between the genes detected in the reference sample (5 ng) and the other samples was good (≥ 90% with 75 and 100 pg), regardless of the variations of the library size. Based on these results and PCA/distribution profiles, we therefore recommend a minimum of 50 pg of total RNA to perform a comprehensive analysis of the expressed genome. Lower quantity of total RNA (up to 10 pg) can be used to raise libraries and sequence successfully, essentially in order to check for the presence of a given gene, but is not fully quantitative. Last, as expected, increasing the number of biological repeats positively impacted the number of DEGs detected (e.g. an additional repeat increases by around 10% the number of DEG detected).

Dissecting the epidermis for RNA-seq profiling at the torpedo stage of embryo development was a real challenge to tackle with Arabidopsis. Thus, we evaluated the approach on two well-known primary and secondary metabolic pathways, namely wax and flavonoids, both occurring preferentially in epidermal cells. Only 4 differentially genes that are supposed to be directly involved in wax and cutin biosynthesis were found: one *MAH-like* (At4g39150), *WAX2* (At5g57800), *CUT1* (At1g68530) and *HOTHEAD* (At1g72970). The well-characterized genes involved in epicuticular waxes and cutin like *CER1*, *CER3* or *CER4* are however not found in the dissected RNA dataset. It has to be noted that the precise biochemical functions of these 4 genes are still unknown. On the contrary, the function of acyl-CoA biosynthetic is better characterized. As expected, *ACC1* gene, which is involved in malonyl-CoA synthesis was found to be preferentially expressed in epidermal cells that is consistent with previous data [18]. Acyl-CoA is associated with epicuticular wax synthesis in developing stems [39]. Several of them were found to have also a differentially enhanced expression in the embryo epidermis. For instance, out of the 21 KCS genes coding the first enzyme of the acyl-CoA elongation complex), eleven showed an induced expression in stem epidermis, among which six were also induced in embryo epidermis [39]. Two out of the three remaining genes of the elongase complex, *KCR* and *ECR*
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Genes involved in VLCFA and WAXES are up regulated in the epidermis. Schematic representation of the very long chain fatty acid and waxes biosynthetic pathway. Adapted from http://aralip.plantbiology.msu.edu/pathways/fatty_acid_elongation_wax_biosynthesis. The red dots highlight the differentially expressed genes</figcaption></figure>
<figure><img src="image_9.png" /></figure>

were also found to be induced in the embryo epidermis. The remaining member HACD coded by *PASTICCINO* gene was strongly expressed in seedling epidermis and was found as the most-differentially expressed gene in the apical part of the globular embryo [8, 26]. Other epicuticular waxes and cutin genes like *CER1*, *CER3* or *CER4* were not isolated, suggesting that they are either not differentially expressed between the two tissues studied or expressed at low levels, below the detection threshold of the methods used. In conclusion, the embryo epidermis showed specific expression of fatty acid genes involved in the epicuticular waxes and cutin biosynthesis leading most probably to an increase in acyl-CoA elongation and in wax or cutin export.

For the flavonoid pathway, most of the EBGs involved in the biosynthesis of flavonoid precursors (i.e. dihydroflavonols) were found to be preferentially expressed in the epidermal cells. This is fully consistent with the simultaneous characterization of *MYB111*, which is involved in the regulation of these EBGs [38] and *ENHANCER OF GLABRA 3* (*EGL3*, At1g63650) acting both as a regulator of epidermal cell fate (see below) and flavonoid

biosynthesis [43]. Moreover, this is also consistent with the usual epidermal localization of flavonoids, which can participate to the cuticle layer, providing some protective barrier against biotic or abiotic stresses [23, 24].

Last, we focused on DE transcription factors (TF). Among the most differentially expressed genes in the epidermal cells compared to the mesophyll cell, we identified *PROTODERMAL FACTOR 2* (*PDF2*) and *ARABIDOPSIS THALIANA MERISTEM LAYER 1* (*AtML1*), two homologous HD-Zip expressed in epidermal cells [1, 27, 32], as well as *WEREWOLF* (*WER*) and *EGL3* that interact to control epidermal cell fate in root [33, 36]. Our results are fully consistent with the epidermal localisation and function of these TF and again, validate the approach. It has to be noticed that *MYB118* was identified as DE expressed with a high log ratio (> 5, Table 3). This transcription factor was functionally analysed and mRNA accumulation as well as promoter analysis showed that its expression is restricted to the endosperm [3]. This clearly suggests a probable contamination of the epidermis by endosperm tissues, probably when the microdissection itself was performed, a little bit outside the epidermis to protect the tissue. Thus the data we provide here are not strictly epidermis specific but represent an epidermis-enriched fraction of the genes expressed during embryogenesis. We could have narrow the dissection beam to the epidermis's cell wall but would probably have decreased the quality of the mRNA during the process. Nevertheless, our results also provide new putative epidermal-specific or preferentially-expressed TFs (Table 3). They pave the way for new functional analyses of these interesting candidates for the regulation of epidermal cell differentiation and metabolic pathways. Many other factors have been identified as DE without preliminary data about their possible localisation, some of them being highly expressed in the epidermis, suggesting major roles in the specificity or in the differentiation of the epidermis in the growing embryo. It is to be noticed that, for many of the TF, the count number was relatively high in both cell layers, suggesting a difference of expression but no cell specificity. For some TF the count number was very low, or even not detected at all, outside the epidermis, suggesting these genes are not only DE in the epidermis compared to mesophyll, but specific for the epidermis. To go further in the characterization of these DE genes, it would be also very interesting to study the promoter sequences of all these epidermal specific or preferentially expressed genes in order to characterize putative conserved epidermis specific cis-elements and then isolate new regulatory genes and build the regulatory networks involved. More broadly, identifying co-expressed genes in specific tissues or cell layers will allow a much more detailed analysis of these regulatory networks.
<table><caption>Table 3 Transcription factors differentially expressed in the epidermis in the developing embryo. TAF from http://datf.cbi.pku.edu.cn/download.php</caption><thead><tr><th>AGI</th><th>Name</th><th></th><th>LogDiff Epidermis/ mesophyll</th><th>Expressed in the embryo</th><th>Localisation</th><th>References</th></tr></thead><tbody><tr><td>At5g17800</td><td>AtMYB56</td><td>myb domain protein 56</td><td>&gt; 10</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G15510</td><td>ATNAC2</td><td>NAC domain containing protein 2</td><td>9.70</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G45980</td><td>WOX8</td><td>WUSCHEL-related homeobox 8</td><td>9.09</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4g31680</td><td></td><td>Transcriptional factor B3 family protein</td><td>6.5</td><td>Yes</td><td></td><td></td></tr><tr><td>AT2G20825</td><td>ULT2</td><td>ULTRAPETALA 2</td><td>5.23</td><td></td><td></td><td></td></tr><tr><td>AT1G49770</td><td></td><td>bHLH</td><td>5.08</td><td></td><td></td><td></td></tr><tr><td>AT3G27785</td><td>MYB118</td><td>myb domain protein 118</td><td>5.03</td><td>Yes</td><td></td><td></td></tr><tr><td>AT2G38470</td><td>#N/A</td><td></td><td>4.97</td><td></td><td>Epidermis</td><td>Suh et al. [39]</td></tr><tr><td>AT1G27730</td><td>STZ</td><td>SALT TOLERANCE ZINC FINGER</td><td>4.71</td><td></td><td></td><td></td></tr><tr><td>AT4G25490</td><td>CBF1</td><td>C-REPEAT/DRE BINDING FACTOR 1</td><td>4.58</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G62470</td><td>MYB96</td><td>myb domain protein 96</td><td>4.43</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G61250</td><td>AtMYB17</td><td>myb domain protein 17</td><td>4.23</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G47230</td><td>ERF5</td><td>ETHYLENE RESPONSIVE ELEMENT BINDING FACTOR 5</td><td>4.15</td><td></td><td></td><td></td></tr><tr><td>AT5G59820</td><td>RHL41</td><td>RESPONSIVE TO HIGH LIGHT 41</td><td>3.95</td><td></td><td></td><td></td></tr><tr><td>AT1G65620</td><td>AS2</td><td>ASYMMETRIC LEAVES 2</td><td>3.83</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G18270</td><td>ANAC087</td><td></td><td>3.81</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G47600</td><td>MYB94</td><td>myb domain protein 94</td><td>3.71</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G46880</td><td>HB-7</td><td>homeobox-7</td><td>3.65</td><td>Yes</td><td>Epidermis</td><td>Nakamura et al. [25]</td></tr><tr><td>AT2G27050</td><td>EIL1</td><td>ETHYLENE-INSENSITIVE3-LIKE 1</td><td>3.54</td><td>Yes</td><td></td><td></td></tr><tr><td>AT5G49330</td><td>AtMYB111</td><td>myb domain protein 111</td><td>3.23</td><td></td><td></td><td></td></tr><tr><td>AT2G36890</td><td>ATMYB38</td><td>myb domain protein 38</td><td>3.11</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4G38620</td><td>MYB4</td><td>myb domain protein 4</td><td>3.06</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4G04890</td><td>PDF2</td><td>PROTODERMAL FACTOR2</td><td>3.03</td><td>Yes</td><td>Epidermis</td><td>Abe et al. [1]</td></tr><tr><td>AT4G21750</td><td>ATML1</td><td>MERISTEM LAYER 1</td><td>2.86</td><td>Yes</td><td>Epidermis</td><td>Sessions et al. [34]</td></tr><tr><td>AT5G14750</td><td>WER</td><td>WEREWOLF 1</td><td>2.78</td><td>Yes</td><td>Epidermis</td><td>Ryu et al. [31]</td></tr><tr><td>AT1G32640</td><td>ATMYC2</td><td>JASMONATE INSENSITIVE 1</td><td>2.67</td><td></td><td></td><td></td></tr><tr><td>AT3G62610</td><td>AtMYB11</td><td>myb domain protein 11</td><td>2.64</td><td>Yes</td><td></td><td>Stracke et al. 2007</td></tr><tr><td>AT1G63650</td><td>EGL3</td><td>ENHANCER OF GLABRA3</td><td>2.51</td><td>Yes</td><td>Epidermis</td><td>Bernhardt et al. [6]</td></tr><tr><td>AT3G52910</td><td>AtGRF4</td><td>GROWTH-REGULATING FACTOR 4</td><td>2.37</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4G01250</td><td>WRKY22</td><td>WRKY DNA-binding protein 22</td><td>2.34</td><td>Yes</td><td></td><td></td></tr><tr><td>AT1G21970</td><td>LEC1</td><td>LEAFY COTYLEDON 1</td><td>2.31</td><td>Yes</td><td>Epidermis</td><td>Lotan et al. [22]</td></tr><tr><td>AT4G16780</td><td>ATHB-2</td><td>Homeobox-leucine zipper protein HAT4</td><td>2.29</td><td></td><td></td><td></td></tr><tr><td>AT4G25470</td><td>CBF2</td><td>FREEZING TOLERANCE QTL 4</td><td>2.18</td><td></td><td></td><td></td></tr><tr><td>AT1G14687</td><td>ATHB32</td><td></td><td>2.14</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G16770</td><td>ATEBP/RAP2.3</td><td></td><td>2.09</td><td>Yes</td><td></td><td></td></tr><tr><td>AT2G45190</td><td>AFO</td><td>ABNORMAL FLORAL ORGANS</td><td>1.99</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4G24240</td><td>WRKY7</td><td></td><td>1.88</td><td>Yes</td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 3 continued**</caption><thead><tr><th>AGI</th><th>Name</th><th></th><th>LogDiff Epidermis/mesophyll</th><th>Expressed in the embryo</th><th>Localisation</th><th>References</th></tr></thead><tbody><tr><td>AT4G31550</td><td>WRKY11</td><td></td><td>1.51</td><td>Yes</td><td></td><td></td></tr><tr><td>AT1G14440</td><td>ATHB31</td><td></td><td>1.46</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G01460</td><td>#N/A</td><td></td><td>1.40</td><td>Yes</td><td></td><td></td></tr><tr><td>AT4G36930</td><td>SPT</td><td>SPATULA</td><td>1.18</td><td>Yes</td><td></td><td></td></tr><tr><td>AT2G34710</td><td>PHB</td><td>PHABULOSA</td><td>1.17</td><td>Yes</td><td></td><td></td></tr><tr><td>AT3G26790</td><td>FUS3</td><td>FUSCA 3</td><td>1.08</td><td>Yes</td><td>Epidermis</td><td>Tsuchiya et al. [42]</td></tr><tr><td>AT3G62670</td><td>ARR20</td><td>ARABIDOPSIS RESPONSE REGULATOR</td><td>0.86</td><td>Yes</td><td></td><td></td></tr></tbody></table>

## Conclusion

We have set up a robust protocol for RNA seq analysis of ultra low quantity of template RNA obtained after laser-assisted-microdissection. We have shown that the protocol is still working below 50 pg of total RNA despite lower number of detected genes. This could be useful for some applications necessarily requiring a very small amount of starting material although not fully quantitative. We have successfully applied this method to the analysis of the epidermis in developing embryos, allowing detection of many differentially expressed transcription factors in this cell layer. This analysis will help for very precise transcriptomic analysis of specific tissues that could be of great use for a better understanding of gene regulation in specific organ or cell types.

## Additional files

**Additional file 1: Fig. S1.** Principal Components Analysis (PCA) on normalized counts of the samples of the pilot experiment (RNA quantity from 5 ng to 10 pg). First and second components are shown, along with the percentage of variance explained.

**Additional file 2: Fig. S2.** Sample clustering based on normalized counts of the 6 samples of the pilot experiment (RNA quantity from 5 ng to 10 pg) after a transformation of the counted reads data as moderated log-counts-per-million. A Euclidean distance is computed between samples, and the dendrogram is built upon the Ward criterion.

**Additional file 3: Fig. S3.** Principal Components Analysis (PCA) on normalized counts of the comparison between Epidermis vs Mesophyll tissues. First and second components are shown, along the percentage of the variance explained.

### Authors' contributions

BD and SB conceived the project, and designed experiments. KS and ND performed LAM experiment. LT and YJ performed RNA-seq libraries experiment. VB, CPL, ED and MLMM performed RNA-seq analysis. BD, LL and JDF analysed the transcriptomic data from Arabidopsis embryo. BD, SB, LL, JDF, VB and MLMM wrote the manuscript. All authors read and approved the final manuscript.

### Author details

### Acknowledgements

The authors also acknowledge the IG-CNS partner for the libraries sequencing. The IJPB and IPS2 institutes benefit from the support of the LabEx Saclay Plant Sciences-SPS (ANR-10-LABX-0040-SPS).

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and materials

The transcriptomic data have been submitted to GEO data depository NGS2014_01_MicroD GSE98176 and NGS2012_02_MicroD GSE98313 accessible through https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=obetg.uoslnifniz&acc=GSE98176 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=kszcacybtzmhwf&acc=GSE98313.

### Consent for publication

Not applicable.

### Ethics approval and consent to participate

Not applicable.

### Funding

This work has benefited from a French State grant (LabEx Saclay Plant Sciences-SPS, ref. ANR-10-LABX-0040-SPS), managed by the French National Research Agency under an "Investments for the Future" program (ref. ANR-11-IDEX-0003-02).

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 14 June 2017 Accepted: 15 January 2018

Published online: 03 February 2018

### References

1. Abe M, Katsumata H, Komeda Y, Takahashi T. Regulation of shoot epidermal cell differentiation by a pair of homeodomain proteins in Arabidopsis. Development. 2003;130:635–43.

2. Bargmann BO, Birnbaum KD. Fluorescence activated cell sorting of plant protoplasts. J Vis Exp. 2010 (36). https://doi.org/10.3791/1673.

3. Barthole G, To A, Marchive C, Brunaud V, Soubigou-Taconnat L, Berger N, Dubreucq B, Lepiniec L, Baud S. MYB118 represses endosperm maturation in seeds of Arabidopsis. Plant Cell. 2014;26:3519–37.

<footnote n="1">¹ Institut Jean-Pierre Bourgin (IJPB), INRA, AgroParisTech, CNRS, Université Paris-Saclay, RD10, 78026 Versailles Cedex, France. ² Institute of Plant Sciences Paris Saclay IPS2, CNRS, INRA, Université Paris-Sud, Université Evry, Université Paris-Saclay, Bâtiment 630, 91405 Orsay, France. ³ Institute of Plant</footnote>

<footnote>Sciences Paris-Saclay IPS2, Paris Diderot, Sorbonne Paris-Cité, Bâtiment 630, 91405 Orsay, France. ⁴ UMR MIA-Paris, AgroParisTech, INRA, Université Paris-Saclay, 75005 Paris, France. ⁵ Present Address: IRHS, Université d'Angers, INRA, AGROCAMPUS-Ouest, SFR4207 QUASAV, Université Bretagne Loire, 49045 Angers, France.</footnote>
4. Becker MG, Hsu SW, Harada JJ, Belmonte MF. Genomic dissection of the seed. Front Plant Sci. 2014;5:464.

5. Belmonte MF, Kirkbride RC, Stone SL, Pelletier JM, Bui AQ, Yeung EC, Hashimoto M, Fei J, Harada CM, Munoz MD, Le BH, Drews GN, Brady SM, Goldberg RB, Harada JJ. Comprehensive developmental profiles of gene activity in regions and subregions of the Arabidopsis seed. Proc Natl Acad Sci USA. 2013;110:E435–44.

6. Bernhardt C, et al. The bHLH genes GL3 and EGL3 participate in an intercellular regulatory circuit that controls cell patterning in the Arabidopsis root epidermis. Development. 2005;132(2):291–8.

7. Braybrook SA, Stone SL, Park S, Bui AQ, Le BH, Fischer RL, Goldberg RB, Harada JJ. Genes directly regulated by LEAFY COTYLEDON2 provide insight into the control of embryo maturation and somatic embryogenesis. Proc Natl Acad Sci USA. 2006;103:3468–73.

8. Casson S, Spencer M, Walker K, Lindsey K. Laser capture microdissection for the analysis of gene expression during embryogenesis of Arabidopsis. Plant J. 2005;42:111–23.

9. Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E. AGRIS: arabidopsis gene regulatory information server, an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinform. 2003;4:25.

10. Deal RB, Henikoff S. The INTACT method for cell type-specific gene expression and chromatin profiling in *Arabidopsis thaliana*. Nat Protoc. 2011;6:56–68.

11. Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.

12. Gagnot S, Tamby JP, Martin-Magniette ML, Bitton F, Taconnat L, Balzergue S, Aubourg S, Renou JP, Lecharny A, Brunaud V. CATdb: a public access to Arabidopsis transcriptome data from the URGV-CATMA platform. Nucl Acids Res. 2008;36:D986–90.

13. Galbraith DW, Birnbaum K. Global studies of cell type-specific gene expression in plants. Annu Rev Plant Biol. 2006;57:451–75.

14. Girard IJ, Mcloughlin AG, de Kievit TR, Fernando DW, Belmonte MF (2016) Integrating large-scale data and RNA technology to protect crops from fungal pathogens. Front Plant Sci. 2016;7:631.

15. Javelle M, Vernoud V, Rogowsky PM, Ingram GC. Epidermis: the formation and functions of a fundamental plant tissue. New Phytol. 2011;189:17–39.

16. Jardinaud MF, Boivin S, Rodde N, Catrice O, Kisiala A, Lepage A, Moreau S, Roux B, Cottret L, Sallet E, Brault M, Emery RJ, Gouzy J, Frugier F, Gamas P. A laser dissection-RNAseq analysis highlights the activation of cytokinin pathways by nod factors in the medicago truncatula root epidermis. Plant Physiol. 2016;171(3):2256–76. https://doi.org/10.1104/pp.16.00711.

17. Jurgens G. Pattern formation in the flowering plant embryo. Curr Opin Genet Dev. 1992;2:567–70.

18. Kajiwara T, Furutani M, Hibara K, Tasaka M. The GURKE gene encoding an acetyl-CoA carboxylase is required for partitioning the embryo apex into three subregions in Arabidopsis. Plant Cell Physiol. 2004;45:1122–8.

19. Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012;28(24):3211–7.

20. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

21. Lepiniec L, Debeaujon I, Routaboul J-M, Baudry A, Pourcel L, Nesi N, Caboche M. Genetics and biochemistry of seed flavonoids. Ann Rev Plant Biol. 2006;57:405–30.

22. Lotan T, Ohto M, Yee KM, West MA, Lo R, Kwong RW, Yamagishi K, Fischer RL, Goldberg RB, Harada JJ. Arabidopsis LEAFY COTYLEDON1 is sufficient to induce embryo development in vegetative cells. Cell. 1998;93(7):1195–205.

23. Mahroug S, Courdavault V, Thiersault M, St-Pierre B, Burlat V. Epidermis is a pivotal site of at least four secondary metabolic pathways in *Catharanthus roseus* aerial organs. Planta. 2006;223:1191–200.

24. Murata J, Roepke J, Gordon H, De Luca V. The leaf epidermome of *Catharanthus roseus* reveals its biochemical specialization. Plant Cell. 2008;20:524–42.

25. Nakamura M, et al. Characterization of the class IV homeodomain-Leucine Zipper gene family in Arabidopsis. Plant Physiol. 2006;141(4):1363–75.

26. Nobusawa A, Kim M, Kaira K, Miyashita G, Negishi A, Oriuchi N, Higuchi T, Tsushima Y, Kanai Y, Yokoo S, Oyama T. Diagnostic usefulness of (1) (8)-F-AMT PET and L-type amino acid transporter 1 (LAT1) expression in oral squamous cell carcinoma. Eur J Nucl Med Mol Imaging. 2013;40:1692–700.

27. Ogawa E, Yamada Y, Sezaki N, Kosaka S, Kondo H, Kamata N, Abe M, Komeda Y, Takahashi T. ATML1 and PDF2 play a redundant and essential role in Arabidopsis embryo development. Plant Cell Physiol. 2015;56:1183–92.

28. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40. https://doi.org/10.1093/bioinformatics/btp616.

29. Roudier F, Gissot L, Beaudoin F, Haslam R, Michaelson L, Marion J, Molino D, Lima A, Bach L, Morin H, Tellier F, Palauqui JC, Bellec Y, Renne C, Miquel M, Dacosta M, Vignard J, Rochat C, Markham JE, Moreau P, Napier J, Faure JD. Very-long-chain fatty acids are involved in polar auxin transport and developmental patterning in Arabidopsis. Plant Cell. 2010;22:364–75.

30. Roux B, Rodde N, Jardinaud MF, Timmers T, Sauviac L, Cottret L, Carrere S, Sallet E, Courcelle E, Moreau S, Debelle F, Capela D, de Carvalho-Niebel F, Gouzy J, Bruand C, Gamas P. An integrated analysis of plant and bacterial gene expression in symbiotic root nodules using laser-capture microdissection coupled to RNA sequencing. Plant J. 2014;77(6):817–37.

31. Ryu JS, et al. Phytochrome-specific type 5 phosphatase controls light signal flux by enhancing phytochrome stability and affinity for a signal transducer. Cell. 2005;120(3):395–406.

32. San-Bento R, Farcot E, Galletti R, Creff A, Ingram G. Epidermal identity is maintained by cell-cell communication via a universally active feedback loop in *Arabidopsis thaliana*. Plant J. 2014;77:46–58.

33. Schiefelbein J, Huang L, Zheng X. Regulation of epidermal cell fate in Arabidopsis roots: the importance of multiple feedback loops. Front Plant Sci. 2014;5:47.

34. Sessions A, Weigel D, Yanofsky MF. The Arabidopsis thaliana MERISTEM LAYER 1 promoter specifies epidermal expression in meristems and young primordia. Plant J. 1999;20(2):259–63.

35. Slane D, Kong J, Berendzen KW, Kilian J, Henschen A, Kolb M, Schmid M, Harter K, Mayer U, De Smet I, Bayer M, Jurgens G. Cell type-specific transcriptome analysis in the early *Arabidopsis thaliana* embryo. Development. 2014;141:4831–40.

36. Song SK, Kwak SH, Chang SC, Schiefelbein J, Lee MM. WEREWOLF and ENHANCER of GLABRA3 are interdependent regulators of the spatial expression pattern of GLABRA2 in Arabidopsis. Biochem Biophys Res Commun. 2015;467:94–100.

37. Spencer MW, Casson SA, Lindsey K. Transcriptional profiling of the Arabidopsis embryo. Plant Physiol. 2007;143:924–40.

38. Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F, Niehaus K, Weisshaar B. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the *Arabidopsis thaliana* seedling. Plant J. 2007;50:660–77.

39. Suh MC, Samuels AL, Jetter R, Kunst L, Pollard M, Ohlrogge J, Beisson F. Cuticular lipid composition, surface structure, and gene expression in Arabidopsis stem epidermis. Plant Physiol. 2005;139:1649–65.

40. Takada S, Jurgens G. Transcriptional regulation of epidermal cell fate in the Arabidopsis embryo. Development. 2007;134:1141–50.

41. ten Hove CA, Lu KJ, Weijers D. Building a plant: cell fate specification in the early Arabidopsis embryo. Development. 2015;142:420–30.

42. Tsuchiya Y, Nambara E, Naito S, McCourt P. The FUS3 transcription factor functions through the epidermal regulator TTG1 during embryogenesis in Arabidopsis. Plant J. 2004;37(1):73–81.

43. Xu W, Lepiniec L, Dubos C. New insights toward the transcriptional engineering of proanthocyanidin biosynthesis. Plant Signaling Behav. 2014;9(4). pii:e28736.

44. Zhang C, Barthelson RA, Lambert GM, Galbraith DW. Global characterization of cell-specific gene expression through fluorescence-activated sorting of nuclei. Plant Physiol. 2008;147:30–40.