Next Article in Journal
Genome Editing in Cowpea Vigna unguiculata Using CRISPR-Cas9
Next Article in Special Issue
Quantitative Proteomic Analysis Reveals Novel Insights into Intracellular Silicate Stress-Responsive Mechanisms in the Diatom Skeletonema dohrnii
Previous Article in Journal
Alterations of Fatty Acid Profile May Contribute to Dyslipidemia in Chronic Kidney Disease by Influencing Hepatocyte Metabolism
Previous Article in Special Issue
Metabolome and Transcriptome Association Analysis Reveals Dynamic Regulation of Purine Metabolism and Flavonoid Synthesis in Transdifferentiation during Somatic Embryogenesis in Cotton
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of Leaf Sheath Maturation in Maize

1
College of Agronomy, Synergetic Innovation Centre of Henan Grain Crops and National Key Laboratory of Wheat and Maize Crop Science, Henan Agricultural University, Zhengzhou 450002, China
2
State Key Laboratory of Crop Biology, College of Agronomic Sciences, Shandong Agricultural University, Tai’an 271018, China
3
The Institute of Tropical Bioscience and Biotechnology (ITBB), Chinese Academy of Tropical Agricultural Sciences (CATAS), Haikou 571101, China
4
Department of Statistics, Iowa State University, Ames, IA 50011, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2019, 20(10), 2472; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102472
Submission received: 2 April 2019 / Revised: 11 May 2019 / Accepted: 17 May 2019 / Published: 19 May 2019

Abstract

:
The morphological development of the leaf greatly influences plant architecture and crop yields. The maize leaf is composed of a leaf blade, ligule and sheath. Although extensive transcriptional profiling of the tissues along the longitudinal axis of the developing maize leaf blade has been conducted, little is known about the transcriptional dynamics in sheath tissues, which play important roles in supporting the leaf blade. Using a comprehensive transcriptome dataset, we demonstrated that the leaf sheath transcriptome dynamically changes during maturation, with the construction of basic cellular structures at the earliest stages of sheath maturation with a transition to cell wall biosynthesis and modifications. The transcriptome again changes with photosynthesis and lignin biosynthesis at the last stage of sheath tissue maturation. The different tissues of the maize leaf are highly specialized in their biological functions and we identified 15 genes expressed at significantly higher levels in the leaf sheath compared with their expression in the leaf blade, including the BOP2 homologs GRMZM2G026556 and GRMZM2G022606, DOGT1 (GRMZM2G403740) and transcription factors from the B3 domain, C2H2 zinc finger and homeobox gene families, implicating these genes in sheath maturation and organ specialization.

1. Introduction

Maize (Zea mays L.) is one of the most important grain crops globally. It is not only used as a grain food and feed but is also preserved as corn silage and as a bioenergy source. Regardless of its final use, the accumulation of maize biomass is important. The leaf is the primary source of photosynthate; thus, its morphological development greatly influences plant biomass and yield. A typical mature maize leaf is made up of three parts: The leaf blade, the ligular region and the leaf sheath. The blade is the distal part of the maize leaf and the major organ for photosynthesis, the ligular region is a wedge-shaped structure connecting the leaf blade and leaf sheath that acts as a hinge to project the leaf blade away from the stem, and the sheath wraps around the stem and provides strength for the growth and development of the leaf blade. The shape and development of these three regions control the architecture of the maize leaf, which is important for photosynthesis and ultimately maize yields [1,2,3,4,5].
The development of grass leaves proceeds basipetally, with the distal cells differentiating and maturing first while the basal cells divide and expand. Several recent studies have portrayed a very dynamic biochemical differentiation process along the maize leaf blade. Here, basic cellular functions—such as DNA synthesis and cell wall synthesis—become enriched in the basal region of the leaf sink, which transitions to secondary cell wall biosynthesis and the establishment of the photosynthetic machinery in the source-sink transition zone and the dominant photosynthetic reactions in the distal part of the leaf [6,7]. Also, several mutants have been identified that influence the morphogenesis of the leaf blade in maize, e.g., ns (narrow sheath) and rld1 (rolled leaf 1) [8,9], and produce abnormal leaf blade development in maize. Several genes that control ligule development in maize have been identified. These include LIGULELESS1 (LG1), which encodes a Squamosa binding protein (SBP) family transcription factor; LIGULELESS2 (LG2), encoding a basic leucine zipper (bZIP) family transcription factor; and LIGULELESS NARROW (LGN), a putative serine-threonine kinase that controls early ligule formation in maize leaves [10,11,12,13,14]. Compared to studies on the leaf blade and ligule tissues, very few studies have focused on the leaf sheath, and our understanding of sheath development is limited.
The leaf sheath wraps around the stem in maize and is believed to provide strength for the growth and development of the leaf blade. Hatfield and colleagues compared changes in the cell wall component in maize leaf blades [15], midribs and sheaths from nodes 9 to 14 and concluded that sheath and midrib tissues always accumulate more neutral sugars, lignin and total phenolics than blade tissues. These metabolite measurements may explain the higher mechanical strength of leaf sheaths; however, information at the molecular level is still limited.
Most mutations that affect the blade-sheath boundary appear to extend the sheath tissue into the blade, but the functions of genes that regulate the development of the sheath tissues are poorly understood, and almost no genes have been identified that distinguish the leaf sheath from the blade and ligule tissues. To investigate the transcriptional dynamics associated with the maturation of maize sheath tissues, we conducted RNA sequencing (RNA-seq) profiling experiments to capture the progressive stages of sheath maturation. These data reveal a dynamic transcriptional process during leaf sheath maturation. In particular, we have defined transcriptional regulators that may specify unique activities in sheath tissue relative to those in the leaf blade. This dataset serves as a foundation for future studies of maize sheath development and maturation.

2. Results

2.1. Defining the Leaf Sheath Transcriptome

In the mature maize leaf, the proximal sheath and distal blade tissues are separated by the leaf ligule (Figure 1A,D). Relative to the leaf blade, epidermal cells in the sheath tissue are smaller and contain fewer cell files, as determined by counting the number of cell files between the same veins across the sheath and blade (Figure 1A,B). To investigate developmental and maturation changes in the sheath tissue, we planted maize seeds every day for 13 days and then harvested the leaf blades and sheaths from the third leaf on the same day, when seedlings ranged in age from 9 to 13 days after planting in order to minimize the environmental variation at the time of harvest. Two millimeter-wide sections of sheath and blade tissues just below or above the ligule were harvested from the third leaf at 10 (stage 1, S1), 11 (stage 2, S2), 12 (stage 3, S3), and 13 (stage 4, S4) days after planting. The sheath tissue was poorly defined on the third leaf at 9 days after planting (stage 0), so we could not harvest tissue from this stage (Figure 1C,D). Although the specification of sheath cell fate is complete in the leaf primordium [16], the sheath tissues we harvested from stage 1 to 4 were still undergoing differentiation and maturation (Figure 1C,D). In total, 221.9 million high-quality reads were obtained from the RNA-seq analysis, of which 164.0 million (73.9%) were uniquely mapped to the maize reference genome. We detected 24,704, 25,225, 25,224, and 24,498 expressed genes from stages 1 to 4, respectively, and 23,557 genes were expressed in all four stages during sheath maturation.

2.2. Maturation Gradient in the Leaf Sheath of Maize

We identified 7918 differentially expressed (DE) genes in the leaf sheath at four maturation stages with a false discovery rate (FDR) controlled at 0.001 (Supplemental Table S1), which accounted for approximately 31% of the sheath transcriptome. We then assigned genes to functional categories and grouped the genes by developmental dynamics using the k-means clustering algorithm. We identified nine clusters (K1–K9; Figure 2; Supplemental Figure S1 and Supplemental Table S1) from the leaf sheath and eight main clusters (K1–K8) that accounted for approximately 99% of the DE genes in four stages of development and maturation. As shown in Figure 2A,B, genes encoding proteins involved in cellular organization (e.g., actin and tubulin), vesicular transport, (e.g., syntaxin and snare), and DNA synthesis/chromatin structure were greatly enriched in clusters K1 and K2 and represented genes required for the installation of the basic cellular infrastructure. These genes are expressed at their highest levels in the early stage of sheath maturation (stage 1). Differences also existed between the K1 and K2 clusters. Namely, genes encoding fatty acid synthesis and elongation-related products, such as 3-ketoacyl-CoA synthase 4, 3-ketoacyl-CoA synthase 6, 3-ketoacyl-CoA synthase 9, 3-ketoacyl-CoA synthase 11 and 3-ketoacyl-CoA synthase 20 were differentially expressed [17]. The expression of genes related to respiration (including glycolysis, the tricarboxylic acid cycle and mitochondrial electron transport); signaling genes, especially G-proteins and LRR receptor kinases; and some transcription factors were greatly enriched in cluster K1 and slightly enriched in cluster K2, and genes involved in protein synthesis and targeting were only enriched in K2. The genes in clusters K3 and K4 were expressed at their highest levels in stage 2 and mainly participate in cell wall biosynthesis/modification and secondary metabolism. The genes that were highly expressed in stage 3 were included in clusters K5–K7 and are mainly involved in pathways related to photosynthesis (K6 and K7), hormone metabolism (K6 and K7), and redox and transcription regulation (K6 and K7), indicating that the sheath tissues at stage 3 begin building photosynthetic machinery. The genes in the K8 cluster were enriched in the functional categories of photosynthesis, sulfate assimilation, and secondary metabolism, e.g., lignin biosynthesis to strengthen the sheath, as well as cell organization, especially actin and fibrillin family proteins that are required for plastoglobule development. This highlighted the functions of both photosynthesis and plant strengthening in the mature sheath of maize. Taken together, these data revealed metabolic changes during the maturation of the sheath tissue in the maize leaf.

2.3. Cell Wall and Lignin Synthesis during Leaf Sheath Maturation in Maize

The cell wall plays an important role in shaping and strengthening cells. We noticed a dramatic change in the expression of cell wall-related genes during sheath maturation (Figure 3A). As shown in Figure 3A, cellulose synthases (CESAs)—the key genes of cellulose biosynthesis in primary and secondary cell walls—exhibited the highest expression levels in stages 1 or 2. Most of them exhibited dramatically decreased expression in stages 3 and 4, during late sheath maturation. Consistent with the expression patterns of the CESAs, most genes encoding CSL (cellulose synthase-like) proteins and cell wall proteins—for example, arabinogalactan proteins (AGP), leucine-rich repeat family proteins (LRR), and reversibly glycosylated polypeptide (RGP)—also showed decreased expression patterns in the late stages (S3 and S4) of sheath maturation. This indicated that the cell wall might be built in stages 1 and 2.
As a main component of the secondary cell wall, lignin plays an important role in providing strength and rigidity to support the cells and the plant body. As shown in Figure 3B (Supplemental Table S2), key enzymes related to lignin biosynthesis [18,19]—e.g., 4-coumarate-CoA ligase (4CL), cinnamoyl CoA reductase (CCR), O-methyltransferase (COMT), and cinnamyl alcohol dehydrogenase (CAD)—exhibited increased expression during the maturation of the sheath, which indicated the accumulation of lignin with the maturation of the leaf sheath tissues.

2.4. Changes in Transcription Factors during Sheath Maturation in Maize

As key factors that regulate gene expression, transcription factors (TFs) play important roles in plant growth, development and the response to various environmental stress. We detected 456 differentially expressed TFs (q < 0.001) during sheath maturation and grouped them into six groups (G1 to G6, Figure 4A and Supplemental Table S3) using a hierarchical clustering (HCL) program. As shown in Figure 4A, groups G3 to G6 included most of the TFs (91%). The TFs in group G3 exhibited their highest expression levels in stage 4 during sheath maturation, and those in G4 were expressed at their peak level in stage 3. The goladen2-like (GLK) family of TFs, including GLK1, which regulates chloroplast development in maize [20], were exclusively enriched in groups G3 and G4 (Figure 4B). The DNA binding with one finger (DOF) family of TFs, such as CDF3, likely regulate photoperiod gene expression [21]. CO-like family members, such as COL3 [22,23], regulate gene expression in photomorphogenesis and during lateral root development and function as a day length-sensitive regulator of shoot branching. The MYB family of TFs, such as MYB4 [24,25], which respond to UV-B and LHY (involved in circadian rhythm [26]), were all mainly enriched in groups G3 and G4. As such, they indicated an increase in photosynthesis in stage 3 and 4 sheath tissues. Groups G5 and G6 included TFs that were highly expressed in stage 1, the earliest maturation stage that we harvested. The trihelix family members were highly enriched in both groups G5 and G6, and the basic helix-loop-helix (bHLH) family (e.g., MUTE) controlled meristemoid differentiation during the early stage of stomatal development [27]. The changes in the expression of TFs suggested the transcriptional regulation of gene expression during sheath maturation.

2.5. Identification of Genes Expressed at High Levels in the Leaf Sheath

To determine whether some genes were specifically expressed in the leaf sheath, we compared the sheath transcriptome with that of the blade at each stage during sheath maturation. In total, 167 genes were expressed at higher levels and 362 genes were expressed at lower levels in the sheath, than that in the blade, across the four maturation stages, with the FDR controlled at 0.001 (Supplemental Table S4).
As shown in Figure 5, of the 362 genes that were expressed at lower levels in the sheath than in the blade, 22% of them participated in the photosynthetic pathway. In contrast, of the 167 genes that were expressed at higher levels in the sheath than in the blade, none participated in photosynthesis. In addition to the photosynthetic pathway, genes involved in tetrapyrrole synthesis and redox scavenging were also specifically expressed at high levels in the leaf blade, consistent with the function of the green leaf blade to harvest light for photosynthesis, and redox reactive species generated by the light need to be controlled by enzymes such as ascorbate peroxidase. Interestingly, for genes in the functional categories of DNA synthesis/chromatin structure, cell wall development, lipid metabolism and cell organization, the enriched genes were expressed at much higher levels in the leaf sheath than in the blade. This may reflect a prolonged phase of cell growth and division in sheath tissues relative to the blade. As almost no sheath-specific genes were reported in a previous study on maize, we focused on identifying genes that were uniquely expressed at high levels in the leaf sheath.
As shown in Table 1, we identified 15 genes that were expressed at very low levels or not significant levels in the leaf blade, while their expression was high in the leaf sheath across the four developmental stages sampled. Among them were two genes (GRMZM2G026556 and GRMZM2G022606) that are homologous to BLADE ON PETIOLE2 (BOP2) in Arabidopsis and rice. BOP2 acts in cells adjacent to lateral organ boundaries to repress genes that confer meristem cell fate and induce genes that promote lateral organ fate and polarity in Arabidopsis, and plays a primary role in rice to regulate the proximal-distal axis [28,29]. DON-glucosyltransferase 1 (DOGT1, GRMZM2G403740), encodes an enzyme that presumably regulates BR activity in Arabidopsis [30,31], and a homolog of SMALLER WITH VARIABLE BRANCHES (SVB, GRMZM2G131409), a protein with a conserved domain of unknown function (DUF538) that influences trichome development in Arabidopsis [32] were also expressed at higher levels in sheath relative to blade. In addition, glycerol-3-phosphate acyltransferase (GPAT2, GRMZM2G033767) and GDSL-like Lipase/Acylhydrolase superfamily protein (GDSL-like Lipase, GRMZM5G862317), which likely participates in lipid metabolism [33,34]; carboxyesterase (CXE18, GRMZM2G104141), which may play a role in cell elongation [35] and senescence-associated gene (SAG12, GRMZM2G061879), which controls nitrogen allocation during senescence in Arabidopsis [36], were all highly expressed in the sheath tissue. This indicated the possible role of these genes in sheath development and maturation.
Transcription factors may also play important roles in distinguishing sheath tissues from blade tissues. The B3 domain-containing transcription factor NGA1, which negatively regulates cell proliferation in Brassica rapa [37,38], is expressed at a high level in the sheath (GRMZM2G082227) tissue only, which may suggest its role in prohibiting the development of sheath tissue into blade tissue. The maize NTT (NO TRANSMITTING TRACT) homologs GRMZM2G071101 and GRMZM2G445684, two members of the C2H2-type zinc finger family that determine the distal cell fate in the root and influence transmitting tract development and pollen tube growth in Arabidopsis [39,40,41], as well as the homeobox family member GRMZM2G034113, the homolog of which (the AtHB7 gene) is involved in ABA signaling during water stress, were all highly expressed in the sheath tissue and exhibited quite low expression in the blade [42,43]. This suggested that these TFs may play a role in the development or tissue-specific function of the leaf sheath. To confirm the reliability of the RNA-seq data, the expression levels of these four transcription factors at the four blades and sheath maturation points discussed above were also quantified by qRT-PCR (Supplemental Figure S2, Supplemental Table S5). The consistency between the results of the two methods verified the reliability of the enriched expression of these TFs in the leaf sheath.

3. Discussion

Cereal crops, such as maize, are a primary calorie source for human and animal diets, and the most important organs for the conversion of photosynthetic energy into carbon are the leaves. Therefore, understanding leaf development in maize is important for agriculture to further increase maize yields. Several studies have uncovered photosynthetic development in the maize leaf blade [6,7]. However, very few studies have focused on leaf sheath tissues. Consequently, very few sheath-specific genes in maize have been identified and our understanding of sheath development is also limited.
Using a comprehensive transcriptome dataset, we captured the dynamic transcriptome changes that occurred during leaf sheath tissue maturation. In the earliest stage of maturation (stage 1), basic cellular structures were constructed, as the genes associated with cell organization, vesicle transport and DNA synthesis/chromatin structure were highly expressed. In stage 2, as the sheath tissue continued to expand, there was an increase in the expression of the genes required for cell wall biosynthesis, modification, and secondary metabolism. By stage 3, the sheath cells were well primed for photosynthesis, as they expressed genes required for photosynthesis and redox reactions at high levels. In addition to genes required for photosynthesis, by stage 4, there was an increased expression of genes participating in secondary metabolism and lignin biosynthesis, which suggests an increase in the strength of the sheath tissues. Therefore, these four stages of leaf sheath differentiation define distinct steps in the sheath maturation process.
We inferred from the RNA-seq data that different tissues in the maize leaf are highly specialized in their biological functions. By comparing the transcriptome of the leaf blade and sheath and identifying the genes expressed at high or very low levels in the sheath tissue across four stages of maturation, we found that increased expression of photosynthesis-related genes in the leaf blade corresponded with its exclusive function in photosynthesis. We also found increased expression of genes that participate in DNA synthesis/chromatin structure, cell wall maintenance, lipid metabolism and cell organization in the sheath tissue, which are involved in building and strengthening the tissue. Interestingly, several genes were highly expressed, specifically in the sheath tissue (e.g., the BOP2 homologs GRMZM2G026556 and GRMZM2G022606). The BOP gene is required for the patterning of the leaf petiole and blade in Arabidopsis [28], and a recent study demonstrated that it activates proximal sheath differentiation and suppresses distal blade differentiation in rice [29], indicating its important role in sheath development. More research needs to be performed to test the functions of DOGT1 (GRMZM2G403740) and SVB (GRMZM2G131409) in maize, even though their Arabidopsis homologs affect hormone activity and influence plant development [30,31,32].
In our survey, we found that TFs may play important roles in both sheath development and functional specialization. During the early stage of sheath maturation, trihelix family and bHLH family members may participate in cell development, and in the late stage of maturation. When the sheath starts to carry out photosynthesis, the members of the GLK, DOF and CO-like families may play roles in light signaling and chloroplast development. Interestingly, TFs from the B3 domain (GRMZM2G082227), C2H2-type zinc finger (GRMZM2G071101 and GRMZM2G445684) and homeobox (GRMZM2G034113) families were expressed at a higher level in the sheath tissue than in the blade, which may suggest a role for these transcription factors in sheath development and functional specialization. Taken together, our results provide detailed insight into dynamic changes during maize leaf sheath maturation and serve as a valuable resource for maize functional genomics study.

4. Materials and Methods

4.1. Plant Material and Sampling

Maize inbred B73 were planted in the growth chamber every day for 13 days with a light intensity of 500 µmol/m2/s, 12:12 L/D, 31 °C L/22 °C D, and 50% relative humidity. The first planting was 13 days, and the blade and sheath of the third leaves were harvested at 2 mm sections above and below the leaf ligule from the seedlings on the 10th (stage 1), 11th (stage 2), 12th (stage 3) and 13th (stage 4) days after planting. Tissue from 10 leaves was harvested as a pool and frozen in liquid nitrogen for RNA extraction and library construction, and 3 biological replicates were harvested at each stage.

4.2. SEM Observation

Scanning electron microscopy was used to observe the cell structure in the maize leaf blade and sheath tissues. Stage 4 seedlings, including the leaf blade, ligule and sheath were fixed in 4% glutaraldehyde in 25 mM sodium phosphate buffer (pH 6.8) at 4  °C for at least 24  h. Then, the samples were rinsed by 25 mM sodium phosphate (pH 6.8) and dehydrated in a graded ethanol series at 4 °C. Thereafter, the samples were immersed in tertiary butyl alcohol and stored in a refrigerator at 4 °C. After sublimation of the frozen tertiary butyl alcohol in a vacuum, the samples were mounted on SEM stubs with double-sided tape, sputter-coated with gold, and examined under a scanning electron microscope (JSM-6610LV, JEOL, Tokyo, Japan).

4.3. RNA-Seq Library Construction and Sequencing

Total RNA was extracted from each sample using the TRIZOL reagent (Invitrogen, Carlsbad, CA, USA). The concentration of the RNA was determined using a DeNovix Spectrophotometer (DeNovix, DS-11, Wilmington, DE, USA), and the RNA quality was determined by 1% agarose gel electrophoresis. The RNA-Seq library construction method of Wang et al. was used [44]. Briefly, the total RNA was purified using the TURBO DNA-free Kit (Ambion, Austin, TX, USA) to completely remove genomic DNA contamination, after that, the poly(A) RNA was isolated from the purified total RNA using poly(T) oligonucleotide-attached magnetic beads (Invitrogen, Carlsbad, CA, USA). Following purification, the mRNA was fragmented into small pieces using divalent cations under elevated temperatures, and the cleaved RNA fragments were reverse-transcribed into first-strand cDNA using reverse transcriptase and random primers. Second-strand cDNA synthesis was performed using DNA polymerase I and RNase H, and the cDNA fragments were processed for end repair, a single adenine base was added, and sequences were ligated to the adapters. These products were then purified and enriched by PCR to create the final cDNA libraries and sequenced on the Illumina Hi-Seq 2500 for single-ended sequencing. Sequencing data have been deposited in the NCBI Sequence Read Archive under the accession number SRP133466.

4.4. Mapping of Reads

The FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit/) was used to remove the adapters of the raw reads with the fastx_clipper setting a parameter of “-a ADAPTER”. The sequence quality was examined by FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) setting “-t 6 –noextract”, and low quality reads were filtered by the FASTX-toolkit setting parameters as “q20p80” (i.e., for each read kept, 80% of bases must have sequence quality greater than 20, which indicates 1% sequencing error rate). The clean reads were mapped to the Maize B73 genome (version 3) obtained from Phytozome using Tophat v2.0.10 (http://tophat.cbcb.umd.edu/) setting -read-mismatches 2 --num-threads 6’. Count data were generated by Cuffdiff embed in Cufflinks pipeline v2.1.1 setting “--num-threads 6 -frag-bias-correct” [45]. The gene expression level was normalized as FPKM, and genes with FPKM > 1 were considered to be expressed.

4.5. Defining Differentially Expressed Genes and Cluster Analysis

Differentially expressed (DE) genes were identified by DESeq v1.34.1 [46] in Bioconductor (http://www.bioconductor.org/), based on a comparison across all the sheath samples (pairwise comparisons with stage 1 as control) or between the blade and sheath in each stage with a false discovery rate (FDR) controlled by the Benjamini and Hochberg method set at 0.001. The MapMan program (http://mapman.gabipd.org) was used to assign genes into functional categories. Clustering analysis of KMC (k-means clustering)—based on the Kendall’s Tau distance metric with maximum iterations set as 50—and HCL support trees, based on Kendall’s Tau distance metric with average linkage method, was performed through the MEV software (http://www.tm4.org/). Functional enrichment analysis was performed using Fisher’s exact test, according to the method by Li et al. [6]. The ID of the DE genes were used to generate a text file. Then, the R software was used to complete the enrichment analysis.

4.6. Quantitative RT-PCR Analysis

cDNA was prepared using the EasyScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TRAN, Beijing, China) and qRT–PCR analyses were conducted using TransStart Tip Green qPCR SuperMix (TRAN, Beijing, China) on a Step One System (Applied Biosystems, Branchburg, NJ, USA). The quantification method (2−∆∆Ct) was used and the variation in expression was estimated using three biological replicates [47]. The maize Ubi2 (UniProtKB/TrEMBL, Q42415, https://www.uniprot.org/statistics/TrEMBL) gene was used as an internal control to normalize the data. The PCR conditions consisted of an initial denaturation step at 94 °C for 30 s, followed by 40 cycles at 94 °C for 5 s and 60 °C for 30 s.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/10/2472/s1. Figure S1. The original graphs showing the K-mean clusters of the differentially expressed genes (DEGs); Figure S2. qRT-PCR and RNA-seq comparison; Table S1. K-mean clusters of DEGs during sheath maturation; Table S2. Expression of genes involved in cell wall and lignin synthesis; Table S3. Hierarchical clustering (HCL) of transcription factors related to sheath maturation in maize; Table S4. Genes expressed at higher or lower levels in the sheath tissue compared with their expression in the blade tissue in maize across four stages of maturation (FDR < 0.001); Table S5. Primers used for qRT-PCR.

Author Contributions

Formal analysis, L.D., L.Q., Z.D. and R.B.; Software, X.D. and P.L.; Supervision, P.L.; Visualization, X.W.; Writing—original draft, P.L.; Writing—review & editing, Y.C., T.P.B. and P.L. All authors reviewed the manuscript.

Acknowledgments

This work is supported by National Key Research and Development Program of China 2016YFD0101003, NSFC 91435108, as well as the Taishan Pandeng program.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Duvick, D. Genetic progress in yield of United States maize (Zea mays L.). Maydica 2005, 50, 193. [Google Scholar]
  2. Kong, F.; Zhang, T.; Liu, J.; Heng, S.; Shi, Q.; Zhang, H.; Wang, Z.; Ge, L.; Li, P.; Lu, X.; et al. Regulation of Leaf Angle by Auricle Development in Maize. Mol. Plant 2017, 10, 516–519. [Google Scholar] [CrossRef] [PubMed]
  3. Lewis, M.W.; Bolduc, N.; Hake, K.; Htike, Y.; Hay, A.; Candela, H.; Hake, S. Gene regulatory interactions at lateral organ boundaries in maize. Development 2014, 141, 4590–4597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lewis, M.W.; Hake, S. Keep on growing: Building and patterning leaves in the grasses. Curr. Opin. Plant Biol. 2016, 29, 80–86. [Google Scholar] [CrossRef] [PubMed]
  5. Tian, F.; Bradbury, P.J.; Brown, P.J.; Hung, H.; Sun, Q.; Flint-Garcia, S.; Rocheford, T.R.; McMullen, M.D.; Holland, J.B.; Buckler, E.S. Genome-wide association study of leaf architecture in the maize nested association mapping population. Nat. Genet. 2011, 43, 159–162. [Google Scholar] [CrossRef] [PubMed]
  6. Li, P.; Ponnala, L.; Gandotra, N.; Wang, L.; Si, Y.; Tausta, S.L.; Kebrom, T.H.; Provart, N.; Patel, R.; Myers, C.R.; et al. The developmental dynamics of the maize leaf transcriptome. Nat. Genet. 2010, 42, 1060–1067. [Google Scholar] [CrossRef] [PubMed]
  7. Wang, L.; Czedik-Eysenberg, A.; Mertz, R.A.; Si, Y.; Tohge, T.; Nunes-Nesi, A.; Arrivault, S.; Dedow, L.K.; Bryant, D.W.; Zhou, W.; et al. Comparative analyses of C(4) and C(3) photosynthesis in developing leaves of maize and rice. Nat. Biotechnol. 2014, 32, 1158–1165. [Google Scholar] [CrossRef]
  8. Juarez, M.T.; Kui, J.S.; Thomas, J.; Heller, B.A.; Timmermans, M.C.P. microRNA-mediated repression of rolled leaf1 specifies maize leaf polarity. Nature 2004, 428, 84–88. [Google Scholar] [CrossRef]
  9. Nardmann, J. The maize duplicate genes narrow sheath1 and narrow sheath2 encode a conserved homeobox gene function in a lateral domain of shoot apical meristems. Development 2004, 131, 2827–2839. [Google Scholar] [CrossRef] [Green Version]
  10. Moreno, M.A.; Harper, L.C.; Krueger, R.W.; Dellaporta, S.L.; Freeling, M. liguleless1 encodes a nuclear-localized protein required for induction of ligules and auricles during maize leaf organogenesis. Genes Dev. 1997, 11, 616. [Google Scholar]
  11. Becraft, P.W.; Bongard-Pierce, D.K.; Sylvester, A.W.; Poethig, R.S.; Freeling, M. The liguleless-1 gene acts tissue specifically in maize leaf development. Dev. Biol. 1990, 141, 220–232. [Google Scholar] [CrossRef]
  12. Sylvester, A.W.; Cande, W.Z.; Freeling, M. Division and differentiation during normal and liguleless-1 maize leaf development. Development 1990, 110, 985–1000. [Google Scholar]
  13. Walsh, J.; Waters, C.A.; Freeling, M. The maize gene liguleless2 encodes a basic leucine zipper protein involved in the establishment of the leaf blade-sheath boundary. Genes Dev. 1998, 12, 208. [Google Scholar]
  14. Moon, J.; Candela, H.; Hake, S. The Liguleless narrow mutation affects proximal-distal signaling and leaf growth. Development 2013, 140, 405–412. [Google Scholar] [CrossRef] [PubMed]
  15. Hatfield, R.D.; Marita, J.M. Maize development: Cell wall changes in leaves and sheaths. Am. J. Plant Sci. 2017, 8, 1248. [Google Scholar] [CrossRef]
  16. Johnston, R.; Wang, M.; Sun, Q.; Sylvester, A.W.; Hake, S.; Scanlon, M.J. Transcriptomic analyses indicate that maize ligule development recapitulates gene expression patterns that occur during lateral organ initiation. Plant Cell 2014, 26, 4718–4732. [Google Scholar] [CrossRef]
  17. Chen, Y.; Kelly, E.E.; Masluk, R.P.; Nelson, C.L.; Cantu, D.C.; Reilly, P.J. Structural classification and properties of ketoacyl synthases. Protein Sci. 2011, 20, 1659–1667. [Google Scholar] [CrossRef] [Green Version]
  18. Boerjan, W.; Ralph, J.; Baucher, M. Lignin biosynthesis. Annu. Rev. Plant Biol. 2003, 54, 519–546. [Google Scholar] [CrossRef]
  19. Barros, J.; Serk, H.; Granlund, I.; Pesquet, E. The cell biology of lignification in higher plants. Ann. Bot. 2015, 115, 1053. [Google Scholar] [CrossRef]
  20. Waters, M.T.; Moylan, E.C.; Langdale, J.A. GLK transcription factors regulate chloroplast development in a cell-autonomous manner. Plant J. Cell Mol. Biol. 2008, 56, 432–444. [Google Scholar] [CrossRef] [Green Version]
  21. Song, Y.H.; Ito, S.; Imaizumi, T. Flowering time regulation: Photoperiod- and temperature-sensing in leaves. Trends Plant Sci. 2013, 18, 575–583. [Google Scholar] [CrossRef] [PubMed]
  22. Datta, S.; Hettiarachchi, G.H.C.M.; Deng, X.-W.; Holm, M. Arabidopsis CONSTANS-LIKE3 Is a Positive Regulator of Red Light Signaling and Root Growth. Plant Cell 2006, 18, 70–84. [Google Scholar] [CrossRef] [PubMed]
  23. Tripathi, P.; Carvallo, M.; Hamilton, E.E.; Preuss, S.; Kay, S.A. Arabidopsis B-BOX32 interacts with CONSTANS-LIKE3 to regulate flowering. Proc. Natl. Acad. Sci. USA 2017, 114, 172–177. [Google Scholar] [CrossRef] [PubMed]
  24. Jin, H.; Cominelli, E.; Bailey, P.; Parr, A.; Mehrtens, F.; Jones, J.; Tonelli, C.; Weisshaar, B.; Martin, C. Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. Embo J. 2014, 19, 6150–6161. [Google Scholar] [CrossRef] [PubMed]
  25. Zhao, J.; Zhang, W.; Zhao, Y.; Gong, X.; Guo, L.; Zhu, G.; Wang, X.; Gong, Z.; Schumaker, K.S.; Guo, Y. SAD2, an importin -like protein, is required for UV-B response in Arabidopsis by mediating MYB4 nuclear trafficking. Plant Cell 2007, 19, 3805–3818. [Google Scholar] [CrossRef]
  26. Shalit-Kaneh, A.; Kumimoto, R.W.; Filkov, V.; Harmer, S.L. Multiple feedback loops of the Arabidopsis circadian clock provide rhythmic robustness across environmental conditions. Proc. Natl. Acad. Sci. USA 2018, 115, 7147–7152. [Google Scholar] [CrossRef] [PubMed]
  27. Pillitteri, L.J.; Bogenschutz, N.L.; Torii, K.U. The bHLH protein, MUTE, controls differentiation of stomata and the hydathode pore in Arabidopsis. Plant Cell Physiol. 2008, 49, 934–943. [Google Scholar]
  28. Ha, C.M.; Jun, J.H.; Nam, H.G.; Fletcher, J.C. BLADE-ON-PETIOLE 1 and 2 control Arabidopsis lateral organ fate through regulation of LOB domain and adaxial-abaxial polarity genes. Plant Cell 2007, 19, 1809–1825. [Google Scholar] [CrossRef]
  29. Toriba, T.; Tokunaga, H.; Shiga, T.; Nie, F.; Naramoto, S.; Honda, E.; Tanaka, K.; Taji, T.; Itoh, J.I.; Kyozuka, J. BLADE-ON-PETIOLE genes temporally and developmentally regulate the sheath to blade ratio of rice leaves. Nat. Commun. 2019, 10, 619. [Google Scholar] [CrossRef] [PubMed]
  30. Poppenberger, B.; Fujioka, S.; Soeno, K.; George, G.L.; Vaistij, F.E.; Hiranuma, S.; Seto, H.; Takatsuto, S.; Adam, G.; Yoshida, S.; et al. The UGT73C5 of Arabidopsis thaliana glucosylates brassinosteroids. Proc. Natl. Acad. Sci. USA 2005, 102, 15253–15258. [Google Scholar] [CrossRef] [Green Version]
  31. Husar, S.; Berthiller, F.; Fujioka, S.; Rozhon, W.; Khan, M.; Kalaivanan, F.; Elias, L.; Higgins, G.S.; Li, Y.; Schuhmacher, R.; et al. Overexpression of the UGT73C6 alters brassinosteroid glucoside formation in Arabidopsis thaliana. BMC Plant Biol. 2011, 11, 51. [Google Scholar] [CrossRef] [PubMed]
  32. Marks, M.D.; Wenger, J.P.; Gilding, E.; Jilk, R.; Dixon, R.A. Transcriptome analysis of Arabidopsis wild-type and gl3-sst sim trichomes identifies four additional genes required for trichome development. Mol. Plant 2009, 2, 803–822. [Google Scholar] [CrossRef] [PubMed]
  33. Beisson, F.; Li, Y.; Bonaventure, G.; Pollard, M.; Ohlrogge, J.B. The acyltransferase GPAT5 is required for the synthesis of suberin in seed coat and root of Arabidopsis. Plant Cell 2007, 19, 351–368. [Google Scholar] [CrossRef] [PubMed]
  34. Lai, C.P.; Huang, L.M.; Chen, L.O.; Chan, M.T.; Shaw, J.F. Genome-wide analysis of GDSL-type esterases/lipases in Arabidopsis. Plant Mol. Biol. 2017, 95, 181–197. [Google Scholar] [CrossRef] [PubMed]
  35. Qin, Y.; Leydon, A.R.; Manziello, A.; Pandey, R.; Mount, D.; Denic, S.; Vasic, B.; Johnson, M.A.; Palanivelu, R. Penetration of the stigma and style elicits a novel transcriptome in pollen tubes, pointing to genes critical for growth in a pistil. PLoS Genet. 2009, 5, e1000621. [Google Scholar] [CrossRef] [PubMed]
  36. James, M.; Poret, M.; Masclaux-Daubresse, C.; Marmagne, A.; Coquet, L.; Jouenne, T.; Chan, P.; Trouverie, J.; Etienne, P. SAG12, a Major Cysteine Protease Involved in Nitrogen Allocation during Senescence for Seed Production in Arabidopsis thaliana. Plant Cell Physiol. 2018, 59, 2052–2063. [Google Scholar] [CrossRef] [PubMed]
  37. Kwon, S.H.; Chang, S.C.; Ko, J.-H.; Song, J.T.; Kim, J.H. Overexpression of Brassica rapa NGATHA1 Gene Confers De-Etiolation Phenotype and Cytokinin Resistance on Arabidopsis thaliana. J. Plant Biol. 2011, 54, 119–125. [Google Scholar] [CrossRef]
  38. Lee, B.H.; Kwon, S.H.; Lee, S.J.; Park, S.K.; Song, J.T.; Lee, S.; Lee, M.M.; Hwang, Y.S.; Kim, J.H. The Arabidopsis thaliana NGATHA transcription factors negatively regulate cell proliferation of lateral organs. Plant Mol. Biol. 2015, 89, 529–538. [Google Scholar] [CrossRef] [PubMed]
  39. Crawford, B.C.; Ditta, G.; Yanofsky, M.F. The NTT gene is required for transmitting-tract development in carpels of Arabidopsis thaliana. Curr. Biol. 2007, 17, 1101–1108. [Google Scholar] [CrossRef]
  40. Crawford, B.C.; Sewell, J.; Golembeski, G.; Roshan, C.; Long, J.A.; Yanofsky, M.F. Plant development. Genetic control of distal stem cell fate within root and embryonic meristems. Science 2015, 347, 655–659. [Google Scholar]
  41. Marsch-Martinez, N.; Zuniga-Mayo, V.M.; Herrera-Ubaldo, H.; Ouwerkerk, P.B.; Pablo-Villa, J.; Lozano-Sotomayor, P.; Greco, R.; Ballester, P.; Balanza, V.; Kuijt, S.J.; et al. The NTT transcription factor promotes replum development in Arabidopsis fruits. Plant J. 2014, 80, 69–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Ré, D.A.; Capella, M.; Bonaventure, G.; Chan, R.L. Arabidopsis AtHB7 and AtHB12 evolved divergently to fine tune processes associated with growth and responses to water stress. BMC Plant Biol. 2014, 14, 1–14. [Google Scholar] [CrossRef]
  43. Olsson, A.; Engström, P.; Söderman, E. The homeobox genes ATHB12 and ATHB7 encode potential regulators of growth in response to water deficit in Arabidopsis. Plant Mol. Biol. 2004, 55, 663–677. [Google Scholar] [CrossRef] [PubMed]
  44. Lin, W.; Si, Y.; Dedow, L.K.; Ying, S.; Peng, L.; Brutnell, T.P. A Low-Cost Library Construction Protocol and Data Analysis Pipeline for Illumina-Based Strand-Specific Multiplex RNA-Seq. PLoS ONE 2011, 6, e26426. [Google Scholar]
  45. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  47. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Observation of the leaf sheath and blade tissues. (A) Scanning electron micrograph of a mature maize leaf. B73 seedlings were grown for 13 days (stage 4). The red line indicates the ligular region, the tissues above the line are the blade and the tissues below the line are sheath tissues. The black line indicates the veins across the sheath, ligule and blade tissues of maize. (B) Cell file numbers and the ratio between the cell length and width were calculated in the leaf blade and sheath tissue between two lateral veins. (C) The blade and sheath tissues above or below the ligule of the third leaf from seedlings 10 (stage 1), 11 (stage 2), 12 (stage 3) and 13 (stage 4) days after planting were harvested at the same time. The arrowheads point to the leaf ligule. (D) Microscopic view of the blade (B), ligule (L) and sheath tissues (S) from stage 0 to stage 4 plants. The arrowheads point to the leaf ligule.
Figure 1. Observation of the leaf sheath and blade tissues. (A) Scanning electron micrograph of a mature maize leaf. B73 seedlings were grown for 13 days (stage 4). The red line indicates the ligular region, the tissues above the line are the blade and the tissues below the line are sheath tissues. The black line indicates the veins across the sheath, ligule and blade tissues of maize. (B) Cell file numbers and the ratio between the cell length and width were calculated in the leaf blade and sheath tissue between two lateral veins. (C) The blade and sheath tissues above or below the ligule of the third leaf from seedlings 10 (stage 1), 11 (stage 2), 12 (stage 3) and 13 (stage 4) days after planting were harvested at the same time. The arrowheads point to the leaf ligule. (D) Microscopic view of the blade (B), ligule (L) and sheath tissues (S) from stage 0 to stage 4 plants. The arrowheads point to the leaf ligule.
Ijms 20 02472 g001
Figure 2. Transcriptome profiling during sheath maturation in maize. (A) K-means clustering showing the expression profile of the maize sheath transcriptome. Eight major clusters (K1–K8) were identified along the four developmental stages (stage 1 to 4) in sheath tissues from 7918 differentially expressed genes. (B) Functional category enrichment was calculated using the MapMan binning method among the eight major clusters in A. The shade of red represents the significance level of log 10 transformed p-values calculated from Fisher’s exact test. Gray blocks indicate functional category enrichment was absent.
Figure 2. Transcriptome profiling during sheath maturation in maize. (A) K-means clustering showing the expression profile of the maize sheath transcriptome. Eight major clusters (K1–K8) were identified along the four developmental stages (stage 1 to 4) in sheath tissues from 7918 differentially expressed genes. (B) Functional category enrichment was calculated using the MapMan binning method among the eight major clusters in A. The shade of red represents the significance level of log 10 transformed p-values calculated from Fisher’s exact test. Gray blocks indicate functional category enrichment was absent.
Ijms 20 02472 g002
Figure 3. Heatmaps showing the expression profiles of genes related to cell wall (A) and lignin biosynthesis (B). Relative gene expression was calculated from the maximum fragments per kilobase million (FPKM) values among the maturation zones. The detailed expression patterns and identities of the genes in each of these biosynthetic pathways are shown in Supplementary Table S2.
Figure 3. Heatmaps showing the expression profiles of genes related to cell wall (A) and lignin biosynthesis (B). Relative gene expression was calculated from the maximum fragments per kilobase million (FPKM) values among the maturation zones. The detailed expression patterns and identities of the genes in each of these biosynthetic pathways are shown in Supplementary Table S2.
Ijms 20 02472 g003
Figure 4. Dynamics of accumulated transcription factor profiles. (A) Hierarchical clustering (HCL) of the transcription factors related to sheath maturation in maize, six lineages (G1 to G6) were identified. The bar represents the expression of transcription factors normalized by row across four stages during sheath maturation in the HCL program. (B) Distribution of transcription factor families in groups G1 to G6 in the leaf sheath.
Figure 4. Dynamics of accumulated transcription factor profiles. (A) Hierarchical clustering (HCL) of the transcription factors related to sheath maturation in maize, six lineages (G1 to G6) were identified. The bar represents the expression of transcription factors normalized by row across four stages during sheath maturation in the HCL program. (B) Distribution of transcription factor families in groups G1 to G6 in the leaf sheath.
Ijms 20 02472 g004
Figure 5. Functional category distribution of genes expressed at high or low levels in the sheath tissue.
Figure 5. Functional category distribution of genes expressed at high or low levels in the sheath tissue.
Ijms 20 02472 g005
Table 1. Genes expressed at high levels in the sheath tissue.
Table 1. Genes expressed at high levels in the sheath tissue.
BladeSheath
IDShort DescriptionSymbolStageStage
12341234
GRMZM2G026556BTB/POZ domain-containing proteinBOP20.30.281.050.739.9118.4731.5638.95
GRMZM2G022606BTB/POZ domain-containing proteinBOP20.40.30.550.39.0910.5814.6420.82
GRMZM2G403740don-glucosyltransferase 1DOGT10.050.060.550.110.372.123.8913.76
GRMZM2G131409Encodes Smaller with Variable BranchesSVB0.580.520.3300.4915.3916.428.07
GRMZM2G033767glycerol-3-phosphate acyltransferase 2GPAT200.551.250.60.8363.1105.9133.77
GRMZM5G862317GDSL-like Lipase/Acylhydrolase superfamily proteinGDSL-like Lipase0.512.793.183.5657.1942.2550.5739.53
GRMZM2G104141carboxyesterase 18CXE1800.722.110.87021.0673.5739.17
GRMZM2G061879senescence-associated gene 12SAG121.416.927.681.381.5851.99136.22196.46
GRMZM2G082227AP2/B3-like transcriptional factor family protein NGATHA1NGA11.512.082.482.47.136.8114.3816.34
GRMZM2G445684C2H2-type zinc finger family protein No Transmitting TractNTT0.220.170.130.1617.6815.7939.7942.81
GRMZM2G071101C2H2-type zinc finger family protein No Transmitting TractNTT0.311.770.750.3410.0819.0945.6842.21
GRMZM2G034113homeobox 7HB-71.670.911.976.7914.2417.743.0843.7
GRMZM5G842695MATE efflux family protein-0.110.60.380.120.123.075.096.33
GRMZM2G343291Protein of unknown function-4.4713.774.610.572.0669.06117.9594.52
GRMZM2G136571alpha/beta-Hydrolases superfamily protein-0.830.771.530.550.637.0334.1225.2
GRMZM2G412436Protein of unknown function-0.982.391.310.172.2518.5615.164.33
Note: The expression of genes was represented by the FPKM value in the table.

Share and Cite

MDPI and ACS Style

Dong, L.; Qin, L.; Dai, X.; Ding, Z.; Bi, R.; Liu, P.; Chen, Y.; Brutnell, T.P.; Wang, X.; Li, P. Transcriptomic Analysis of Leaf Sheath Maturation in Maize. Int. J. Mol. Sci. 2019, 20, 2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102472

AMA Style

Dong L, Qin L, Dai X, Ding Z, Bi R, Liu P, Chen Y, Brutnell TP, Wang X, Li P. Transcriptomic Analysis of Leaf Sheath Maturation in Maize. International Journal of Molecular Sciences. 2019; 20(10):2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102472

Chicago/Turabian Style

Dong, Lei, Lei Qin, Xiuru Dai, Zehong Ding, Ran Bi, Peng Liu, Yanhui Chen, Thomas P. Brutnell, Xianglan Wang, and Pinghua Li. 2019. "Transcriptomic Analysis of Leaf Sheath Maturation in Maize" International Journal of Molecular Sciences 20, no. 10: 2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102472

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop