Impact Factor 3.677

The world's most-cited Plant Sciences journal

Frontiers in Plant Science

Plant Genetics and Genomics

This article is part of the Research Topic

Melatonin in plants

Original Research ARTICLE

Front. Plant Sci., 31 August 2016 | https://doi.org/10.3389/fpls.2016.01289

RNA-seq Analysis of Overexpressing Ovine AANAT Gene of Melatonin Biosynthesis in Switchgrass

Shan Yuan1, Yanhua Huang1,2, Sijia Liu1, Cong Guan1, Xin Cui1, Danyang Tian1, Yunwei Zhang1,3,4* and Fuyu Yang1,5*
  • 1College of Animal Science and Technology, China Agricultural University, Beijing, China
  • 2College of Agriculture, China Agricultural University, Beijing, China
  • 3Beijing Key Laboratory for Grassland Science, China Agricultural University, Beijing, China
  • 4National Energy R&D Center for Biomass, Beijing, China
  • 5Beijing Sure Academy of Biosciences, Beijing, China

Melatonin serves important functions in the promotion of growth and anti-stress regulation by efficient radical scavenging and regulation of antioxidant enzyme activity in various plants. To investigate its regulatory roles and metabolism pathways, the transcriptomic profile of overexpressing the ovine arylalkylamine N-acetyltransferase (oAANAT) gene, encoding the penultimate enzyme in melatonin biosynthesis, was compared with empty vector control using RNA-seq in switchgrass, a model plant of cellulosic ethanol conversion. The 85.22 million high quality reads that were assembled into 135,684 unigenes were generated by Illumina sequencing for transgenic oAANAT switchgrass with an average sequence length of 716 bp. A total of 946 differentially expression genes in transgenic line comparing to control switchgrass, including 737 up-regulated and 209 down-regulated genes, were mainly enriched with two main functional patterns of melatonin identifying by gene ontology analysis: the growth regulator and stress tolerance. Furthermore, KEGG maps indicated that the biosynthetic pathways of secondary metabolite (phenylpropanoids, flavonoids, steroids, stilbenoid, diarylheptanoid, and gingerol) and signaling pathways (MAPK signaling pathway, estrogen signaling pathway) were involved in melatonin metabolism. This study substantially expands the transcriptome information for switchgrass and provides valuable clues for identifying candidate genes involved in melatonin biosynthesis and elucidating the mechanism of melatonin metabolism.

Introduction

Since melatonin (N-acetyl-5-methoxytryptamine) was first discovered in plant in Dubbels et al. (1995) and Hattori et al. (1995), varying concentrations of melatonin have been found in many other plant species (Murch et al., 1997; Manchester et al., 2000; Simopoulos et al., 2005). The fundamental issues of melatonin biosynthetic pathways and physiological functions still need to be deciphered to utilize in plants (Tan et al., 2003, 2012; Arnao and Hernández-Ruiz, 2009; Arnao, 2014; Wei et al., 2016). Both exogenous melatonin treatments (Kolár et al., 2003; Arnao and Hernández-Ruiz, 2007; Sarropoulou et al., 2012; Zhang et al., 2013) and melatonin-rich transgenic plants (Kang et al., 2010; Park et al., 2012; Byeon and Back, 2014; Wang et al., 2014) have been carried on in order to determine its potential functional roles. The physiological role of melatonin in plants involved growth regulation, scavenging reactive oxygen species and increases of antioxidant enzyme activities (Paredes et al., 2008; Arnao and Hernández-Ruiz, 2014; Hardeland, 2016). Through the encoding genes for the catalytic reactions during melatonin biosynthesis have been identified and cloned recently, its biosynthetic pathways still not clearly deciphered yet in plants (Kang et al., 2011, 2013; Byeon et al., 2013a; Hardeland, 2016). The classic pathways of melatonin biosynthesis consist of the four catalytic reactions from tryptophan: decarboxylation by tryptophan decarboxylase (TDC), hydroxylation by tryptamine 5-hydroxylase (T5H) to serotonin, N-acetylation by serotonin N-acetyltransferase (SNAT) and the final O-methylation to melatonin by N-acetylserotonin O-methyltransferase (ASMT). The penultimate step catalyzes the same reaction as the non-homologous AANAT of vertebrate, which catalyzes conversion of serotonin into N-acetylserotonin (Byeon et al., 2013a). Overexpression of AANAT gene can promote the ability of biosynthesis of melatonin, and significantly improve the melatonin content in plants (Kang et al., 2010; Wang et al., 2014; Zhang et al., 2014).

Switchgrass (Panicum virgatum L.) is a perennial C4 grass native to North America, and is well-researched on germplasm collection, cultivation, genetic breeding as a model plant for cellulosic bioethanol production over the past several decades (Sanderson et al., 1996; Parrish and Fike, 2005). Aiming at the production demand of large biomass, the traditional breeding techniques have been challenged by its nature of outcrossing and polyploidy, as well as the general infertile environments for the biomass grass cultivation. Fortunately, modification of the functional genes that related to the growth and resistance could effectively assist the breeding process of switchgrass (Fu et al., 2011; Poovaiah et al., 2014; Baxter et al., 2015; Wuddineh et al., 2015). Several transgenic plants expressing the key genes for melatonin biosynthesis enzymes have been demonstrated to increase the resistance to environmental stresses (Kang et al., 2010; Park et al., 2012).

In this study, the transgenic melatonin-rich switchgrass (overexpressing of the oAANAT gene) and the transgenic control with empty vector (EV) were used to conduct a RNA-seq and analyze the effects of melatonin on gene expression. This research intends to provide differentially expressed genes (DEGs) in the melatonin-rich switchgrass and mechanism information on the biosynthesis, regulation, metabolism of melatonin for further investigations in plants.

Materials and Methods

Plant Materials and Cultivation

The switchgrass (Panicum virgatum L. var. Alamo) plants expressing sheep AANAT gene (ID: 25120) were used for this study. The transgenic plants were grown under 16 h light (26°C, 120 μmol/m2/s) and 8 h dark (18°C) conditions with watered every other day. Fully matured plants were chosen from each genotype for molecular characterization and transcriptome sequencing, which were all phenotypically identical with respect to overall size, tiller number and leaf pattern. Three replicates of relative high melatonin contents for each transgenic stem were frozen in liquid nitrogen and stored at 80°C until further analysis.

Characterization of Growth and Development of Transgenic Switchgrass

Tiller number, plant height, internode number, internode length, internode diameter, leaf blade length, leaf blade width, root number, root length, root diameter, and spike length were measured at the transgenic reproductive third (R3) stage (Figure 1), a stage with fully emerged spikelets and an emerged peduncle (Hardin et al., 2013). Internode 3 (I3) was used for measuring internode diameter. The leaves of I3 were used to measure leaf blade length and leaf blade width. Twelve replicates were measured for each transgenic line.

FIGURE 1
www.frontiersin.org

FIGURE 1. The phenotypes of transgenic oAANAT line and the control switchgrass in reproductive third stage.

RNA Isolation and Qualification

Total RNA was extracted from the stems using the TRIzol reagent method (Invitrogen, Carlsbad, CA, USA) and was treated with RNase-free DNase or 30 min at 37°C to remove residual DNA. RNA degradation and contamination was monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer spectrophotometer (Implen, Westlake Village, CA, USA). RNA concentration was determined by Qubit RNA Assay Kit in Qubit 2.0 Flurometer. RNA integrity was assessed with the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA).

Transcriptome Sample Preparation for Sequencing

The total amount of 1.5 μg RNA per sample was prepared for the RNA-seq. Sequencing libraries were generated using NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, USA). The cDNA fragments of 150∼200 bp were preferentially selected from the library by purification with AMPure XP system (Beckman Coulter, Beverly, MA, USA).

Clustering and Sequencing

The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia). After cluster generation, paired-end reads were generated by sequencing with the library preparations on an Illumina Hiseq platform.

Validation of RNA-seq Data by Real-Time Quantitative

Real-time quantitative PCR validation of RNA-seq data for 12 random genes was performed in 20 μl of reaction mixture (Supplementary Table S1). The amount of the amplified DNA was monitored by fluorescence at the end of each cycle using 7500 Real-Time PCR System (Applied Biosystems). Each plate was repeated three times in independent runs for all reference and selected genes. Gene expression was evaluated by the 2-ΔΔCt method (Livak and Schmittgen, 2001).

Data Analysis

Quality Control

Clean data (clean reads) were obtained by removing reads containing adapter, reads containing ploy-N and low quality reads from raw data. Q20, Q30, GC-contents and sequence duplication level of the clean data were calculated. Our raw data have been uploaded to NCBI1, and the accession number is PRJNA322585.

Transcriptome Assembly

The complete genome of switchgrass has not been released through the draft version was released in 2012 and updated to V3. De novo transcriptome assembly was accomplished using Trinity (Grabherr et al., 2011) with min_kmer_cov set to 2 by default and all other parameters set default.

Gene Functional Annotation

Gene function was annotated based on the following databases: Nr (NCBI non-redundant protein sequences); Nt (NCBI non-redundant nucleotide sequences); Pfam (Protein family); KOG/COG (Clusters of Orthologous Groups of proteins); Swiss-Prot (A manually annotated and reviewed protein sequence database); KO (KEGG Ortholog database); GO (Gene Ontology).

Differential Expression Analysis

Gene expression levels were estimated by RSEM (Li and Dewey, 2011). Differential expression analysis of two conditions/groups was performed using the DESeq R package (1.10.1). Genes with an adjusted P-value < 0.05 found by DESeq were assigned as differentially expressed using the Benjamini and Hochberg’s approach.

GO Enrichment Analysis and KEGG Pathway Enrichment Analysis

GOseq R packages were implemented to define the DEGs based Wallenius non-central hyper-geometric distribution (Young et al., 2010). We used KOBAS (Mao et al., 2005) software to test the statistical enrichment of differential expression genes in KEGG pathways2 (Kanehisa et al., 2008).

Results

Morphological Characterization of Transgenic Plants

The transgenic switchgrass exhibited the largely enhanced growth condition (shoot, leaf, and root) compared with the EV (Figure 1; Table 1). Average tillers, shoot height, stem nodes, third internode length, and stem diameter in transgenic line were 45.5, 80.2, 48.3, 63.3, and 22.99% higher than those of EV, respectively (P < 0.05). Leaf blade length was 23.4% longer in transgenic switchgrass (52.72 cm) than that of EV (43.53 cm, P < 0.05). There was no significant difference on leaf width between the two groups (P > 0.05). Average roots, root length, and root diameter were 24.2, 25.6, and 41.5% higher in melatonin-rich switchgrass than those of control groups. Mean spike length in AANAT transgenic switchgrass (23.17 cm) was more than fourfold of the control (5.4 cm, P < 0.05; Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Morphological characterization of transgenic switchgrass plants in reproductive third stage.

Illumina Paired-End Sequencing and Assembly

Total RNA was extracted from transgenic and control switchgrass in order to sequence using Illumina paired-end sequencing technology. In this study, exceeded 74G clean bases were acquired, and the average GC-rich content and the Q30 level of the six samples was 58.45 and 93.64%, respectively (Supplementary Table S2). After removal of adaptors and low-quality reads, average of 85,227,460 and 81,104,143 clean reads (95.9 and 94.2% of the raw data) of transgenic line and control were obtained, respectively (Supplementary Table S2). These reads were assembled into 264,869 transcripts with an average length of 1,052 bp and an N50 of 1,792 bp. After compared the different transcripts representing one unigene, the longest length transcript for each unigene was extracted. A total of 135,684 unigenes were obtained. The average length was 716 bp, and transcripts with lengths of more than 500 bp accounted for about 37.65% of all transcripts (Supplementary Table S3).

Annotation of All Non-redundant Unigenes

For the validation and annotation of the assembled unigenes, all assembled unigenes were submitted to a BLASTx search with an E value threshold of 1e-5 against the following databases: Nr (NCBI non-redundant protein sequences), Nt (NCBI non-redundant nucleotide sequences), Pfam (Protein family), Swiss-Prot (a manually annotated and reviewed protein sequence database), GO (Gene Ontology), KOG (eukaryotic orthologous groups) and KEGG (Kyoto Encyclopedia of Genes and Genomes). The unigenes were subjected to public databases for similarity searching. Among 135,684 unigenes, 56,968 (41.98%), 605,126 (44.59%) and 38,266 (28.2%) unigenes showed homology with sequences in the NCBI Nr, Nt and SwissProt databases, respectively (Supplementary Table S4; Supplementary Figure S1).

Functional Classification by GO and KOG

Gene ontology (GO, a standardized classification system for gene function) was assigned to classify the functions of predicted switchgrass unigenes. In total, 1,931 functional GO terms were assigned among 41,494 unigenes with BLAST matching to known proteins (Supplementary Table S4). The majority of the unigenes were assigned to the categories of biological processes (100,887, 47.61%), followed by cellular components (63,207, 29.83%) and molecular functions (47,821, 22.57%; Supplementary Figure S2). Under the category of biological processes, cellular processes (22,564, 22.37%) and metabolic processes (21,358, 21.17%) were prominently represented. Under the classification of molecular functions, the binding (22,620, 47.3%) and catalytic activities (17,760, 37.1%) represented the two largest categories, while other categories, such as those for transporter activities, structural molecule activity, molecular function regulator, and others, together contained only 7,441 unigenes representing 15.56% of the total number of unigenes. As for the cellular component, two categories, pertaining to cells and cell parts, accounted for approximately 40.38% of the cellular components that were identified; the organelle category accounted for approximately 13.53% of the cellular component unigenes, and the membrane and membrane part categories accounted for 19.43%.

In order to predict and classify possible functions, all unigenes were aligned to the euKaryotic Ortholog Groups (KOG) database in which orthologous gene products were classified. Out of 56,968 unigenes with significant similarity to nr proteins in this study, 21,165 sequences were assigned to KOG classifications. Among the 26 KOG categories, the cluster related to general function prediction (3,674, 17.36%) was the largest group, followed by those for post-translational modification (2,670, 12.62%); translation, translation, ribosomal structure, and biogenesis (1,693, 8.00%); and signal transduction mechanisms (1,605, 7.58%).

Functional Classifications Using KEGG Pathways

All the unigenes were analyzed with respect to the KEGG pathway database to further examine the transcriptome of transgenic switchgrass. Out of the 135,684 identified unigenes, 16,222 (11.96%) with significant matches were assigned to five main categories that included 131 KEGG pathways (Figure 2; Supplementary file-KEGG). Among the five main categories that were identified, metabolism held the greatest number of unigenes (9,257, 57.06%), followed by genetic information (4,484, 27.64%), cellular processes (1,000, 6.17%), organismal systems (816, 5.03%), and environmental information processing (665, 4.10%). These results indicate that active metabolic processes were occurring in transgenic switchgrass. As shown in Supplementary file-KEGG, the metabolism category contained 19 sub-categories, including environmental adaptation, nucleotide metabolism, metabolism of terpenoids and polyketides, energy metabolism, carbohydrate metabolism, folding, sorting and degradation, membrane transport, biosynthesis involved in other secondary metabolism, amino acid metabolism, transport and catabolism, signal transduction and (Figure 2).

FIGURE 2
www.frontiersin.org

FIGURE 2. Pathway assignment based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. (A) Classification based on cellular process categories, (B) classification based on environmental information processing categories, (C) classification based on genetic information processing categories, (D) classification based on metabolism categories, and (E) classification based on organismal systems categories.

Differential Gene Expression in Melatonin-Rich Transgenic Switchgrass

To reveal the molecular events occurring during the transgenic process, the digital gene expression libraries were constructed using RNA from the pools of control and the transgenic RNA samples and sequenced using Illumina technology. After Illumina sequencing and the removal of adaptors and low-quality reads, approximately 76,976,922,81,777,632, and 84,557,876 reads were obtained for the three control replicates, and 76,912,922, 88,415,412, and 90,354,046 reads were obtain for the three replicates for the transgenic lines. We then mapped the clean reads to the transcriptome reference data, and a total of 54,349, 56,895 and 73,352 unigene sequences were identified for the control replicates, and 53,595, 56,080 and 57,158 unigene sequences were identified for the drought replicates. After calculating the expression level of each mapped unigene, a total of 946 unigenes were detected that had levels of expression that were significantly different between the transgenic and control libraries, including 737 up-regulated and 209 down-regulated unigenes (Figure 3). The significantly higher average FPKM of AANAT gene of transgenic lines (648.53) than that of control group (0.11).

FIGURE 3
www.frontiersin.org

FIGURE 3. Volcano plot of the DEGs between the transgenic AANAT and control switchgrass.

The FPKM values were used to normalize gene expression levels and to compare their differences among transgenic lines and control. The percentage of high level expressed genes (FPKM beyond 15) in transgenic lines was slightly greater than in control samples, while that of low level expressed genes (FPKMs in the interval 0.1–3.57) in transgenic lines was relatively smaller than control samples (Supplementary Table S5).

Functional Annotation and Classification of the DEGs by Gene Ontology Analysis

Among the differentially expressed unigenes, 558 genes were significantly up-regulated in the melatonin-rich switchgrass by more than twofold of the levels in the EV (P < 0.05). These genes included F-box protein (a gene for the controlled degradation of cellular protein); disease resistance protein (a resistance protein guard the plant against pathogens); abscisic stress-ripening protein 3 (an abscisic acid-, stress-, and ripening-induced protein); heat shock protein 83 (a gene for promoting the maturation), structural maintenance and proper regulation of specific target proteins involved for instance in cell cycle control and signal transduction.

To identify the genes that are differentially expressed in transgenic lines, a functional categorization was carried out by GO analysis. A total of 511 unigenes, including 415 up-regulated genes and 96 down-regulated genes, were functionally assigned to the three categories of the GO database (Figure 4). The 2012 GO terms were functionally annotated with GO terms in transgenic and control groups. The GO terms of ‘organonitrogen compound biosynthetic process,’ ‘cellular amide metabolic process,’ ‘amide biosynthetic process,’ ‘aminoglycan metabolic process,’ and ‘peptide biosynthetic process’ in biological process were highly enriched in the DEGs (corrected P-value < 0.05), further suggesting the efficiency of the melatonin biosynthesis and the reliability of the gene expression data. Other terms, such as ‘cytoplasmic part,’ ‘non-membrane-bounded organelle,’ ‘ribonucleoprotein complex’ in cellular component, ‘structural molecule activity’ and ‘chitin binding’ in molecular function, were also significantly overrepresented (corrected P-value < 0.05). By comparing transgenic and control groups, 1687 up-regulated and 1094 down-regulated GO items.

FIGURE 4
www.frontiersin.org

FIGURE 4. Gene ontology (GO) classifications of DEGs across three comparisons. The Y-axis represents the number of DEGs in a category. The results of transgenic A line vs. EV (control) are summarized in three main categories: biological process (BP), cellular component (CC), and molecular function (MF).

KEGG Pathway Analysis of the Melatonin-Related Genes

To determine whether the melatonin-related genes engaged in specific pathways, the DEGs were used as objects to search against the KEGG pathway database. The top 20 obviously enriched pathways are shown in Figure 5. By comparing transgenic with control, ‘ribosome’ pathway enriched the most DEGs (Figure 6), followed by ‘oxidative phosphorylation,’ ‘Glutathione metabolism,’ ‘MAPK signaling pathway,’ ‘photosynthesis – antenna proteins’ and other pathways, suggesting that these pathways and processes might participate in the melatonin synthesis and metabolism (Supplementary file-KEGG).

FIGURE 5
www.frontiersin.org

FIGURE 5. Kyoto Encyclopedia of Genes and Genomes enrichments of the annotated DEGs across three comparisons. The left Y-axis indicates the KEGG pathway. The X-axis indicates the Rich factor. A high q value is represented by blue, and a low q value is represented by red. (a) transgenic A line vs. EV (control); (b) transgenic H line vs. EV (control); (c) transgenic A vs. H line.

FIGURE 6
www.frontiersin.org

FIGURE 6. The ribosome pathway enriched by KEGG analysis.

Transcription Factors

In this study, a total of 2214 transcription factors (TFs) were identified and classified into 79 different families, and the largest group of TFs was the MYB family (147, 6.64%), followed by NAC (108, 4.79%), Orphans (108, 4.88%), C2H2 (102, 4.62%), GRAS (101, 4.56%), C3H (95, 4.29%), bZIP (94, 4.29%), and WRKY (103, 6.71%; Table 2). These results further suggested that 33 DEGs encoding known or putative TFs were changed among transgenic and control groups, including the MYB, NAC, SNF2, and FAR1 TFs (Table 2).

TABLE 2
www.frontiersin.org

TABLE 2. Selected genes about transcription factor with altered expression (P < 0.05) in the two groups.

Validation of Gene Expression Profiles Using RT-qPCR

The 12 DEGs were randomly selected for qRT-PCR. Histograms were generated by comparing the FPKM determined by transcriptome analysis and qRT-PCR. A highly significant correlation (R2 = 0.703, P < 0.01) was found between the qPCR and RNA-Seq, indicating reproducibility and credible RNA-seq data (Supplementary Figure S3).

Discussion

Although, a growing body of molecular and gene expression research regarding melatonin has been documented on model plants (Kang et al., 2010; Park et al., 2012; Byeon et al., 2013b; Weeda et al., 2014), the expression profiling of genes on melatonin biosynthesis and metabolism in forage grass has not yet been investigated, which was potentially valuable for the molecular breeding of the switchgrass as an superior bioethanol grass. In this study, to offer an initial insight into the related genes of melatonin in plants, we conducted a RNA-seq analysis to generate the global transcriptomic profile in transgenic switchgrass plants overexpressing the sheep AANAT gene.

Overexpression of the oAANAT gene, encoding the key enzyme during melatonin biosynthesis process, significantly enhanced melatonin contents (over threefold of the EV from our previous work, P < 0.01) and promoted growth (shoot height, root length, stem diameter, leaf size) and reproductive (spike length) processes in transgenic switchgrass (Table 1). Through the classic view supports the last HIOMT/ASMT enzyme is rate-limiting during the biosynthesis of melatonin, the promotional roles of AANAT/SNAT were also reported in other species, e.g., the oAANAT transgenic ‘micro-tom’ tomato have higher melatonin levels than control (Wang et al., 2014). The positive regulatory role of melatonin on growth and development was confirmed in other transgenic melatonin-rich plants (Byeon and Back, 2014; Wang et al., 2014). For example, the transgenic oAANAT rice exhibited altered height, biomass and panicle numbers per plant, suggesting that melatonin took part in plant growth and reproduction (Byeon and Back, 2014). Moreover, exogenous melatonin treatments also enhanced root regeneration, photosynthetic pigments, total carbohydrates and biomass in the Chenopodium rubrum (Kolár et al., 2003), cherry rootstock (Sarropoulou et al., 2012); cucumber (Zhang et al., 2013), and soybean (Wei et al., 2015). The low concentrations melatonin increases photosynthetic activity by inducing porphyrin and chlorophyll biosynthesis, in contrast, high concentrations of melatonin induces the synthesis of proline and carbohydrate, which are beneficial for the osmoregulation of plants under stresses (Sarropoulou et al., 2012). Moreover, the stimulation of root generation and vitality and addition of the root/shoot ratio under melatonin treatment supported its effect on strengthening cucumber roots (Zhang et al., 2013). Furthermore, we detected the earlier flowering of transgenic oAANAT switchgrass than EV. Several flower-specific genes, FLC (Flowering Locus C, c52347_g1), AP2 (APETALA2, c34036_g2, c48942_g1), DELLA (c47591_g2), differentially expressed between transgenic oAANAT and EV plants (P < 0.05). However, more flowering related genes did not showed significantly regulated by overexpression of oAANAT genes (LFY, CO, TFL, AG, AP3). This might be attributed by the discrepancies of the gene expressions between plant stems and inflorescences, considering that the large number of putative floral-specific transcripts were detected in flowers of sweet potato using RNA-seq (Tao et al., 2013). The striking promotions on development patterns of stems, leaves and roots and alteration of flowering suggested the possible involvement of the melatonin in these physiological actions (Park, 2011; Arnao and Hernández-Ruiz, 2014; Hardeland, 2016).

Transcriptomic analysis exhibited that the large number of 946 DEGs in the melatonin-rich switchgrass comparing to the control (Figure 3). Among the DEGs, functional patterns of melatonin were divided into two main aspects identifying by gene ontology analysis: the growth regulator and stress tolerance. The growth related function of melatonin was supported by the GO terms of plant development (‘post-embryonic root development,’ ‘seed germination,’ ‘seedling development,’ and ‘pollen tube development’) and biosynthetic processes (‘organonitrogen compound biosynthetic process,’ ‘amide biosynthetic process,’ and ‘peptide biosynthetic process’). The genes encoding for signaling regulation and TFs such as MYB domain-containing protein, NAC and C2H2 type protein, leucine-rich repeat and zinc-finger genes were involved in plant growth and metabolisms (Table 2). These differentially expressed TFs were also detected in the other transcriptome analysis of growth related phenotypes in switchgrass (Palmer et al., 2012; Wang et al., 2012). The DREB subfamily, AP2 TFs, several different classes of zinc finger TFs (C2H2), auxin responsive TFs (ARF), Myb family TFs, and NAC were identified in the Arabidopsis dormancy related gene set (Wang et al., 2012). The APETALA2/ethylene response factor (AP2/ERF) superfamily of TFs plays essential roles in the regulation of various growth and developmental programs including stress responses (Wuddineh et al., 2015). This may be closely associated with the gene-regulation capacity of melatonin as a plant growth regulator (Park, 2011; Arnao and Hernández-Ruiz, 2014).

Furthermore, the stress tolerance function of melatonin was evidenced by the GO terms of ‘response to oxidative stress’ and ‘cellular response to stress.’ The F-box protein (the controlled degradation of cellular protein); disease resistance protein (a resistance protein guard the plant against pathogens); abscisic stress-ripening protein 3 (an abscisic acid-, stress-, and ripening-induced protein); heat shock protein 83 (a gene for promoting the maturation) involving for cell cycle control and signal transduction were differentially expressed in transgenic oAANAT switchgrass (Table 2). Similarly, the oxidative resistances were also detected in transgenic rice (Kang et al., 2010) and agrochemicals (Park et al., 2012) owing to the increased levels of endogenous melatonin. The promotion of the maturation, structural maintenance and proper regulation of specific target proteins was consistent with the result of a previous study in which a microarray analysis on the AANAT overexpressed rice, suggesting that melatonin is involved with stress responses (Byeon et al., 2013b). In addition, exogenous melatonin enhances abiotic tolerance (water, cold, salt stresses) in various plants (Zhang et al., 2013; Fan et al., 2015; Wei et al., 2015). The inhibitory effects of salt stress on gene expressions related to binding, oxidoreductase activity/process, and secondary metabolic processes under salt stress were alleviated by melatonin in soybean (Wei et al., 2015). Melatonin achieves its promotional roles by influencing gene expressions involved in metabolic activities: nitrogen metabolism, major carbohydrate metabolism, hormone metabolism, and secondary metabolism (Paredes et al., 2008; Nawaz et al., 2015).

Kyoto Encyclopedia of Genes and Genomes maps provided much information on exploration of metabolic pathways involved in melatonin and deconstruction of its biological functions. Many differentially expressed unigenes were enriched into the biosynthetic pathways of secondary metabolite (phenylpropanoids, flavonoids, steroids, stilbenoid, diarylheptanoid, and gingerol). Plant steroid is used as signaling molecules for physiological and developmental regulation and offers exciting potentials for enhancing crop yield (Vriet et al., 2012). In addition, signaling pathways (MAPK signaling pathway, estrogen signaling pathway) were significantly unregulated as well (Figure 5). Taking the MAPK (Mitogen-activated protein kinase) pathway as example, responses to various biotic and abiotic stresses as an integral component of cellular signaling during mitogenesis and differentiation of mitotic cells (Atkins et al., 1998; Group et al., 2002; Supplementary file-KEGG). The decades of secondary metabolic process and signal pathways related genes were differentially expressed in both endogenous melatonin-rich transgenic rice (Byeon et al., 2013b) and exogenous melatonin treatment (Zhang et al., 2014). These significantly altered expressions of metabolic pathways were not completely consistent between endogenous and exogenous melatonin treatment: the method of spraying leaves or coating seeds with melatonin merely drives short-term changes in morphological and physiological traits, while overexpression of endogenous melatonin gene enables sustainable synthesis of melatonin during the lifespan (Park, 2011). Besides, overexpression genes of the melatonin biosynthesis key enzyme participates in the de novo synthesis and is possibly involved in much more metabolism pathways than those from exogenous treatment, thus providing valuable information for the elucidation the biosynthesis pathway and functional mechanism of melatonin in plants (Kang et al., 2011, 2013; Byeon et al., 2013a,b; Zhang et al., 2016).

Previous studies have revealed that melatonin had significant effect in regulating hormones (such as ABA, GA4) in plant response to salinity and drought stress (Li et al., 2012; Zhang et al., 2016). Melatonin shared the common substrate (tryptophan) with IAA, however, the relationship of biosynthesis between melatonin and IAA is still in suspense and the independent relationship has been proposed in Arabidopsis root system architecture (Pelagio-Flores et al., 2012; Arnao and Hernández-Ruiz, 2014; Wei et al., 2016). In our study, the auxin-responsive protein IAA genes differentially expressed in the melatonin-rich switchgrass comparing with the control, and possibly resulted from the competitive relation for the same biosynthesis substrate. Thus, genome-wide expression analysis of melatonin-related genes in response to plant hormones supported the crosstalk between melatonin and IAA biosynthesis processes, and the intern mechanism of the relative expression needs more clues from other direct analysis of hormones treatments (Zhang et al., 2014; Wei et al., 2016).

Conclusion

This study is the first report of the transcriptomic profile of endogenous melatonin effects on in bioenergy crop switchgrass using RNA-seq technology. Switchgrass is an identified model species for bioethanol but the genetic background and genome information are not well-established, and these transcriptomic datasets will provide fundamental information and serve as new tools to genetically dissect melatonin-mediated pathways in other common grasses. The analysis based on DEGs reveals broad roles of melatonin in regulating plant growth, development and defense systems. Furthermore, the expression of many genes involved in signaling regulation such as MAPK signaling pathway, was also altered in response to the transgenic switchgrass of overexpression of melatonin biosynthesis gene. Taken together, studies on sequencing of transgenic melatonin-rich switchgrass suggests that melatonin plays a critical role in promotion of plant growth and may facilitate identification of melatonin’s functions in plants.

Author Contributions

Conceived and designed the experiments: YZ, FY. Performed the experiments: SY, YH, CG, SL, DT. Analyzed the data: SY, CG, XC. Wrote the paper: SY, YH, YZ. All authors reviewed and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer Y-DG declared a shared affiliation, though no other collaboration, with the authors to the handling Editor, who ensured that the process nevertheless met the standards of a fair and objective review.

Acknowledgments

Funding for this work was provided by the Ministry of Science and Technology, China (2014BAD23B03), National Natural Science Foundation of China (31272493), Natural Science Foundation of Beijing (6162016) and China Agricultural University (2014FG062).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01289

Footnotes

  1. ^http://www.ncbi.nlm.nih.gov/bioproject
  2. ^http://www.genome.jp/kegg/

References

Arnao, M. B. (2014). Phytomelatonin: discovery, content, and role in plants. Adv. Bot. 2014:815769. doi: 10.1155/2014/815769

CrossRef Full Text | Google Scholar

Arnao, M. B., and Hernández-Ruiz, J. (2007). Melatonin promotes adventitious- and lateral root regeneration in etiolated hypocotyls of Lupinus albus, L. J. Pineal Res. 42, 147–152. doi: 10.1111/j.1600-079X.2006.00396.x

CrossRef Full Text | Google Scholar

Arnao, M. B., and Hernández-Ruiz, J. (2009). Assessment of different sample processing procedures applied to the determination of melatonin in plants. Phytochem. Anal. 20, 14–18. doi: 10.1002/pca.1083

CrossRef Full Text | Google Scholar

Arnao, M. B., and Hernández-Ruiz, J. (2014). Melatonin: plant growth regulator and/or biostimulator during stress? Trends Plant Sci. 19, 789–797. doi: 10.1016/j.tplants.2014.07.006

CrossRef Full Text | Google Scholar

Atkins, C. M., Selcher, J. C., Petraitis, J. J., Trzaskos, J. M., and Sweatt, J. D. (1998). The MAPK cascade is required for mammalian associative learning. Nat. Neurosci. 1, 602–609. doi: 10.1038/2836

CrossRef Full Text | Google Scholar

Baxter, H. L., Poovaiah, C. R., Yee, K. L., Mazarei, M., Rodriguez, M., Thompson, O., et al. (2015). Field evaluation of transgenic switchgrass plants overexpressing PvMYB4 for reduced biomass recalcitrance. Bioenergy Res. 8, 1–12. doi: 10.1007/s12155-014-9570-1

CrossRef Full Text | Google Scholar

Byeon, Y., and Back, K. (2014). An increase in melatonin in transgenic rice causes pleiotropic phenotypes, including enhanced seedling growth, delayed flowering, and low grain yield. J. Pineal Res. 56, 408–414. doi: 10.1111/jpi.12129

CrossRef Full Text | Google Scholar

Byeon, Y., Lee, K., Park, Y. I., Park, S., and Back, K. (2013a). Molecular cloning and functional analysis of serotonin N-acetyltransferase from the cyanobacterium Synechocystis sp. pcc 6803. J. Pineal Res. 55, 371–376. doi: 10.1111/jpi.12080

CrossRef Full Text | Google Scholar

Byeon, Y., Park, S., Kim, Y. S., and Back, K. (2013b). Microarray analysis of genes differentially expressed in melatonin-rich transgenic rice expressing a sheep serotonin N-acetyltransferase. J. Pineal Res. 55, 357–363. doi: 10.1111/jpi.12077

CrossRef Full Text | Google Scholar

Dubbels, R., Reiter, R. J., Klenke, E., Goebel, A., Schnakenberg, E., Ehlers, C., et al. (1995). Melatonin in edible plants identified by radioimmunoassay and by high performance liquid chromatography-mass spectrometry. J. Pineal Res. 18, 28–31. doi: 10.1111/j.1600-079X.1995.tb00136.x

CrossRef Full Text | Google Scholar

Fan, J., Hu, Z., Xie, Y., Chan, Z., Chen, K., Amombo, E., et al. (2015). Alleviation of cold damage to photosystem II and metabolisms by melatonin in Bermudagrass. Front. Plant Sci. 6:925. doi: 10.3389/fpls.2015.00925

CrossRef Full Text | Google Scholar

Fu, C., Mielenz, J. R., Xiao, X. R., Ge, Y. X., Hamilton, C. Y., Rodriguez, M. R., et al. (2011). Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc. Natl. Acad. Sci. U.S.A. 108, 3803–3808. doi: 10.1073/pnas.1100310108

CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

CrossRef Full Text | Google Scholar

Group, M., Ichimura, K., Shinozaki, K., Tena, G., Sheen, J., Henry, Y., et al. (2002). Mitogen-activated protein kinase cascades in plants: a new nomenclature. Trends Plant Sci. 7, 301–308. doi: 10.1016/S1360-1385(02)02302-6

CrossRef Full Text | Google Scholar

Hardeland, R. (2016). Melatonin in plants – diversity of levels and multiplicity of functions. Front. Plant Sci. 7:198. doi: 10.3389/fpls.2016.00198

CrossRef Full Text | Google Scholar

Hardin, C. F., Fu, C., Hisano, H., Xiao, X., Shen, H., Stewart, C. N., et al. (2013). Standardization of switchgrass sample collection for cell wall and biomass trait analysis. Bioenergy Res. 6, 755–762. doi: 10.1007/s12155-012-9292-1

CrossRef Full Text | Google Scholar

Hattori, A., Migitaka, H., Iigo, M., Itoh, M., Yamamoto, K., Ohtani-Kaneko, R., et al. (1995). Identification of melatonin in plants and its effects on plasma melatonin levels and binding to melatonin receptors in vertebrates. Biochem. Mol. Biol. Int. 35, 627–634.

Google Scholar

Kanehisa, M., Araki, M., Goto, S., Hattori, M., Hirakawa, M., Itoh, M., et al. (2008). KEGG for linking genomes to life and the environment. Nucleic Acids Res. 36, 480–484. doi: 10.1093/nar/gkm882

CrossRef Full Text | Google Scholar

Kang, K., Kong, K., Park, S., Natsagdorj, U., Kim, Y. S., and Back, K. (2011). Molecular cloning of a plant N-acetylserotonin methyltransferase and its expression characteristics in rice. J. Pineal Res. 50, 304–309. doi: 10.1111/j.1600-079X.2010.00841.x

CrossRef Full Text | Google Scholar

Kang, K., Kyungjin, L., Sangkyu, P., Yeong, B., and Kyoungwhan, B. (2013). Molecular cloning of rice serotonin N-acetyltransferase, the penultimate gene in plant melatonin biosynthesis. J. Pineal Res. 55, 7–13. doi: 10.1111/jpi.12011

CrossRef Full Text | Google Scholar

Kang, K., Lee, K., Park, S., Kim, Y. S., and Back, K. (2010). Enhanced production of melatonin by ectopic overexpression of human serotonin N-acetyltransferase plays a role in cold resistance in transgenic rice seedlings. J. Pineal Res. 49, 176–182. doi: 10.1111/j.1600-079X.2010.00783.x

CrossRef Full Text | Google Scholar

Kolár, J., Johnson, H., and Machácková, I. (2003). Exogenously applied melatonin (N-acetyl-5-methoxytryptamine) affects flowering of the short-day plant Chenopodium rubrum. Physiol. Plant. 118, 605–612. doi: 10.1034/j.1399-3054.2003.00114.x

CrossRef Full Text | Google Scholar

Li, B., and Dewey, C. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12:323. doi: 10.1186/1471-2105-12-323

CrossRef Full Text | Google Scholar

Li, C., Wang, P., Wei, Z., Liang, D., Liu, C., Yin, L., et al. (2012). The mitigation effects of exogenous melatonin on salinity-induced stress in Malus hupehensis. J. Pineal Res. 53, 298–306. doi: 10.1111/j.1600-079X.2012.00999.x

CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408. doi: 10.1016/j.synthmet.2005.07.344

CrossRef Full Text | Google Scholar

Manchester, L. C., Tan, D. X., Reiter, R. J., Park, W., Monis, K., and Qi, W. (2000). High levels of melatonin in the seeds of edible plants: possible function in germ tissue protection. Life Sci. 67, 3023–3029. doi: 10.1016/S0024-3205(00)00896-1

CrossRef Full Text | Google Scholar

Mao, X., Cai, T., Olyarchuk, J. G., and Wei, L. (2005). Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 21, 3787–3793. doi: 10.1093/bioinformatics/bti430

CrossRef Full Text | Google Scholar

Murch, S. J., Simmons, C. B., and Saxena, P. K. (1997). Melatonin in feverfew and other medicinal plants. Lancet 350, 1598–1599. doi: 10.1016/S0140-6736(05)64014-7

CrossRef Full Text | Google Scholar

Nawaz, M. A., Huang, Y., Bie, Z., Ahmed, W., Reiter, R. J., Niu, M., et al. (2015). Melatonin: current status and future perspectives in plant science. Front. Plant Sci. 6:1230. doi: 10.3389/fpls.2015.01230

CrossRef Full Text | Google Scholar

Palmer, N. A., Saathoff, A. J., Kim, J., Benson, A., Tobias, C. M., Twigg, P., et al. (2012). Next-generation sequencing of crown and rhizome transcriptome from an upland, tetraploid switchgrass. Bioenergy Res. 5, 649–661. doi: 10.1007/s12155-011-9171-1

CrossRef Full Text | Google Scholar

Paredes, S. D., Korkmaz, A., Manchester, L. C., Tan, D. X., and Reiter, R. J. (2008). Phytomelatonin: a review. J. Exp. Bot. 60, 57–69. doi: 10.1093/jxb/ern284

CrossRef Full Text | Google Scholar

Park, S., Lee, D. E., Jang, H., Byeon, Y., Kim, Y. S., and Back, K. (2012). Melatonin-rich transgenic rice plants exhibit resistance to herbicide-induced oxidative stress. J. Pineal Res. 54, 258–263. doi: 10.1111/j.1600-079X.2012.01029.x

CrossRef Full Text | Google Scholar

Park, W. J. (2011). Melatonin as an endogenous plant regulatory signal: debates and perspectives. J. Plant Biol. 54, 143–149. doi: 10.1007/s12374-011-9159-6

CrossRef Full Text | Google Scholar

Parrish, D. J., and Fike, J. H. (2005). The biology and agronomy of switchgrass for biofuels. Crit. Rev. Plant Sci. 24, 423–459. doi: 10.1080/07352680500316433

CrossRef Full Text | Google Scholar

Pelagio-Flores, R., Muñoz-Parra, E., Ortiz-Castro, R., and López-Bucio, J. (2012). Melatonin regulates Arabidopsis root system architecture likely acting independently of auxin signaling. J. Pineal Res. 53, 279–288. doi: 10.1111/j.1600-079X.2012.00996.x

CrossRef Full Text | Google Scholar

Poovaiah, C. R., Mazarei, M., Decker, S. R., Turner, G. B., Sykes, R. W., Davis, M. F., et al. (2014). Transgenic switchgrass (Panicum virgatum, L.) biomass is increased by overexpression of switchgrass sucrose synthase (PvSUS1). Biotechnol. J. 10, 552–563. doi: 10.1002/biot.201400499

CrossRef Full Text | Google Scholar

Sanderson, M. A., Reed, R. L., McLaughlin, S. B., Wullschleger, S. D., Conger, B. V., Parrish, D. J., et al. (1996). Switchgrass as a sustainable bioenergy crop. Bioresour. Technol. 56, 83–93. doi: 10.1016/j.fusengdes.2006.04.059

CrossRef Full Text | Google Scholar

Sarropoulou, V., Dimassi-Theriou, K., Therios, I., and Koukourikou-Petridou, M. (2012). Melatonin enhances root regeneration, photosynthetic pigments, biomass, total carbohydrates and proline content in the cherry rootstock PHL-C (Prunus avium × Prunus cerasus). Plant Physiol. Biochem. 61, 162–168. doi: 10.1016/j.plaphy.2012.10.001

CrossRef Full Text | Google Scholar

Simopoulos, A. P., Tan, D., Manchester, L. C., and Reiter, R. J. (2005). Purslane: a plant source of omega-3 fatty acids and melatonin. J. Pineal Res. 39, 331–332. doi: 10.1111/j.1600-079X.2005.00269.x

CrossRef Full Text | Google Scholar

Tan, D. X., Hardeland, R., Manchester, L. C., Korkmaz, A., Ma, S., Rosales-Corral, S., et al. (2012). Functional roles of melatonin in plants, and perspectives in nutritional and agricultural science. J. Exp. Bot. 63, 577–597. doi: 10.1093/jxb/err256

CrossRef Full Text | Google Scholar

Tan, D. X., Manchester, L. C., Hardeland, R., Lopez-Burillo, S., Mayo, J. C., Sainz, R. M., et al. (2003). Melatonin: a hormone, a tissue factor, an autocoid, a paracoid, and an antioxidant vitamin. J. Pineal Res. 34, 75–78. doi: 10.1034/j.1600-079X.2003.02111.x

CrossRef Full Text | Google Scholar

Tao, X., Gu, Y. H., Jiang, Y. S., Zhang, Y. Z., and Wang, H. Y. (2013). Transcriptome analysis to identify putative floral-specific genes and flowering regulatory-related genes of sweet potato. Biosci. Biotechnol. Biochem. 77, 2169–2174. doi: 10.1271/bbb.130218

CrossRef Full Text | Google Scholar

Vriet, C., Russinova, E., and Reuzeau, C. (2012). Boosting crop yields with plant steroids. Plant Cell 24, 842–857. doi: 10.1105/tpc.111.094912

CrossRef Full Text | Google Scholar

Wang, L., Zhao, Y., Reiter, R. J., He, C., Liu, G., Lei, Q., et al. (2014). Changes in melatonin levels in transgenic ‘micro-tom’ tomato overexpressing ovine AANAT, and ovine HIOMT genes. J. Pineal Res. 56, 134–142. doi: 10.1111/jpi.12105

CrossRef Full Text | Google Scholar

Wang, Y., Zeng, X., Iyer, N. J., Bryant, D. W., Mockler, T. C., and Mahalingam, R. (2012). Exploring the switchgrass transcriptome using second-generation sequencing technology. PLoS ONE 7:e34225. doi: 10.1371/journal.pone.0034225

CrossRef Full Text | Google Scholar

Weeda, S., Zhang, N., Zhao, X., Ndip, G., Guo, Y., Buck, G. A., et al. (2014). Arabidopsis transcriptome analysis reveals key roles of melatonin in plant defense systems. PLoS ONE 9:e93462. doi: 10.1371/journal.pone.0093462

CrossRef Full Text | Google Scholar

Wei, W., Li, Q. T., Chu, Y. N., Reiter, R. J., Yu, X. M., Zhu, D. H., et al. (2015). Melatonin enhances plant growth and abiotic stress tolerance in soybean plants. J. Exp. Bot. 66, 695–707. doi: 10.1093/jxb/eru392

CrossRef Full Text | Google Scholar

Wei, Y., Zeng, H., Wei, H., Chen, L., He, C., and Shi, H. (2016). Comparative transcriptional profiling of melatonin synthesis and catabolic genes indicates the possible role of melatonin in developmental and stress responses in rice. Front. Plant Sci. 7:676. doi: 10.3389/fpls.2016.00676

CrossRef Full Text | Google Scholar

Wuddineh, W. A., Mazarei, M., Turner, G. B., Sykes, R. W., Decker, S. R., Davis, M. F., et al. (2015). Identification and molecular characterization of the switchgrass AP2/ERF transcription factor superfamily, and overexpression of PvERF001 for improvement of biomass characteristics for biofuel. Front. Bioeng. Biotechnol. 3:101. doi: 10.3389/fbioe.2015.00101

CrossRef Full Text | Google Scholar

Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, 1–14. doi: 10.1186/gb-2010-11-2-r14

CrossRef Full Text | Google Scholar

Zhang, N., Sun, Q., Li, H., Li, X., Cao, Y., Zhang, H., et al. (2016). Melatonin improved anthocyanin accumulation by regulating gene expressions and resulted in high reactive oxygen species scavenging capacity in cabbage. Front. Plant Sci. 7:197. doi: 10.3389/fpls.2016.00197

CrossRef Full Text | Google Scholar

Zhang, N., Zhang, H. J., Zhao, B., Sun, Q. Q., Cao, Y. Y., Li, R., et al. (2014). The RNA-seq approach to discriminate gene expression profiles in response to melatonin on cucumber lateral root formation. J. Pineal Res. 56, 39–50. doi: 10.1111/jpi.12095

CrossRef Full Text | Google Scholar

Zhang, N., Zhao, B., Zhang, H. J., Weeda, S., Yang, C., Yang, Z. C., et al. (2013). Melatonin promotes water-stress tolerance, lateral root formation, and seed germination in cucumber (Cucumis sativus, L.). J. Pineal Res. 54, 15–23. doi: 10.1111/j.1600-079X.2012.01015.x

CrossRef Full Text | Google Scholar

Keywords: switchgrass, melatonin, overexpression, oAANAT, RNA-seq

Citation: Yuan S, Huang Y, Liu S, Guan C, Cui X, Tian D, Zhang Y and Yang F (2016) RNA-seq Analysis of Overexpressing Ovine AANAT Gene of Melatonin Biosynthesis in Switchgrass. Front. Plant Sci. 7:1289. doi: 10.3389/fpls.2016.01289

Received: 15 July 2016; Accepted: 12 August 2016;
Published: 31 August 2016.

Edited by:

Haitao Shi, Hainan University, China

Reviewed by:

Ruediger Hardeland, University of Göttingen, Germany
Yang-Dong Guo, China Agricultural University, China

Copyright © 2016 Yuan, Huang, Liu, Guan, Cui, Tian, Zhang and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yunwei Zhang, zywei@126.com Fuyu Yang, yfuyu@cau.edu.cn