such as deoxynivalenol (DON), which remain in the food chain and constitute a threat to .... sion (10,000 macroconidia/mL), DON (5 g/L in water), or mock. The.
Aug 12, 2015 - Asian soybean rust (ASR) caused by Phakopsora pachyrhizi is a devastating foliar disease affecting soybean production worldwide. Understanding nonhost resistance against ASR may provide an avenue to engineer soybean to confer durable r
cally, such as low-density lipoprotein (LDL)-cholesterol, lack of .... zymatic colorimetric reactions using commercial kits (Spinreact, Girona,. Spain; values in ...
NADH/NAD+ ratio and also affected the photosynthesis .... They all belong to the NAD(P)H: quinone ...... spray IIIsource (AB SCIEX, Concord, ON) and a pulled.
Jan 26, 2018 - PORA, CHLM, CHLH, PsaK, PsaD, and GUN4 are potential target genes of GLK in .... growth chamber (16-h light/8-h dark, 27â¦C/19â¦C). Three lines were verified as ... User Manual (Clontech, Palo Alto, CA, USA). The domain-.
Aug 13, 2014 - This paper is available online free of all access charges (see ..... formed using the CLUSTER program (http://bonsai.ims.u-tokyo.
May 27, 2016 - analysis revealed that melatonin acted through activation of ROS-scavenging and ROS signal transduction pathways, including antioxidant ...
Aug 17, 2017 - evidence that apoptosis plays a central role during infection with this virus in vitro and in mice (Parquet et al., 2001; Yang et al.,. 2002, 2008 ...
Feb 7, 2017 - Ca2+, protein kinase, ROS, and hormones (auxin, abscisic acid (ABA), and gibberellin (GA)) play important roles in cold-induced gene expression19,20. ...... 11A (Leici Instrument Factory, Shanghai, China) and a last conductivity reading
mitosis/meiosis transition and the meiosis/post-meiotic ..... testes sections from adult CD-1 mice were fixed in Bouin's solution ..... 0.01 based the t test.
Aug 28, 2015 - organic acids, amino acids, and nucleotides), is a modern technique often used for ... Human Development Brain and Tissue Bank for Developmental .... To test for correlations between lutein concentrations and age, data ..... HPC: hippo
Mar 24, 2017 - sequencing and investigated 4 subgroups of transcriptomic profiles, which ...... of mammalian brain RNA-Seq data utilizing WGCNA and Mantel ...
Jun 19, 2017 - ... National University, Seoul, South Korea, 2 College of Pharmacy and .... lanosterol synthase; MEP, 2-C-methyl-D-erythritol 4-phosphate; MVK, ...
in RPTEC/TERT1 cells in response to Cyclosporine A. Anja Wilmes0, Lydia ... full omics analysis of the effects of CsA in cultured renal epi- thelial cells. Hepatic ...
Mar 12, 2008 - Background: Our knowledge on zebrafish reproduction is very limited. ... We found that five members of the GTP-binding septin gene family, from which only one .... female embryonic gonads , microarray analyses to identify genes ...
Feb 17, 2017 - Intracellular ionic and redox homeostases are essential for the ... DEGs in calcium homeostasis and signaling pathway (CDC31, CMK2, LCB5, ...
Sep 21, 2015 - Salvia miltiorrhiza is a well-known traditional Chinese medicine, and .... In all, we extracted 132 and 125 independent analytes that represented putative ..... The RNA samples > 300 ng Î¼ Lâ1 with high purity (OD260/280 > 1.8, ...
Feb 10, 2016 - Fermentation of Penicillium expansum KACC 40815. .... which consisted of an LC pump (Varian 212), a photodiode array detector (ProStar 335 ...
Feb 4, 2015 - RNA-binding protein that binds the 30S ribosomal subunit, required for ...... Tanaka, Y.; Tsumoto, K. Crystal structure of the enzyme CapF of.
Sep 14, 2015 - Citation: Westrop GD, Williams RAM, Wang L,. Zhang T, Watson DG, Silva AM, et al. (2015). Metabolomic Analyses of Leishmania Reveal Multiple. Species Differences and Large Differences in Amino. Acid Metabolism. PLoS ONE 10(9): e0136891
Sep 2, 2011 - Pnpla8. Patatin-like phospholipase domain containing 8. Ð2.72. 24. Bile acid .... lipases Pnpla3, Pnpla8, Cel, Lpl, and Pnliprp1 (Boverhof et al.,.
Aug 15, 2005 - Immunoglobulin deposition, redox imbalance, and impaired .... of the aortic root (measured as averaged lesion per cross-section of the aortic sinus). .... chain reactions. ... the pathogenetic outcome is determined by the balance ... M
Dec 9, 2016 - were edited in Jalview  and Adobe Illustrator CS6. We retrieved 188 sequences from the InterProt database (http://www.ebi.ac.uk/interpro/), GenBank database (http://www.ncbi.nlm.nih.gov/) or the available literature corresponding to
International Journal of
Molecular Sciences Article
Transcriptomic and GC-MS Metabolomic Analyses Reveal the Sink Strength Changes during Petunia Anther Development Yuanzheng Yue 1,2 , Shaoze Tian 1 , Yu Wang 1 , Hui Ma 1 , Siyu Liu 1 , Yuqiao Wang 1 and Huirong Hu 1, * 1
Key Laboratory of Urban Agriculture in Central China, Ministry of Agriculture, Key Laboratory of Horticultural Plant Biology, Ministry of Education, College of Horticulture and Forestry Sciences, Huazhong Agricultural University, Wuhan 430070, China; [email protected] (Y.Y.); [email protected] (S.T.); [email protected] (Y.W.); [email protected] (H.M.); [email protected] (S.L.); [email protected] (Y.W.) Key Laboratory of Landscape Architecture, Jiangsu Province, College of Landscape Architecture, Nanjing Forestry University, Nanjing 210037, China Correspondence: [email protected]; Tel.: +86-158-2742-1300
Received: 6 February 2018; Accepted: 18 March 2018; Published: 23 March 2018
Abstract: Petunia, which has been prevalently cultivated in landscaping, is a dicotyledonous herbaceous flower of high ornamental value. Annually, there is a massive worldwide market demand for petunia seeds. The normal development of anther is the necessary prerequisite for the plants to generate seeds. However, the knowledge of petunia anther development processes is still limited. To better understand the mechanisms of petunia anther development, the transcriptomes and metabolomes of petunia anthers at three typical development stages were constructed and then used to detect the gene expression patterns and primary metabolite profiles during the anther development processes. Results suggested that there were many differentially-expressed genes (DEGs) that mainly participated in photosynthesis and starch and sucrose metabolism when DEGs were compared between the different development stages of anthers. In this study, fructose and glucose, which were involved in starch and sucrose metabolism, were taken as the most important metabolites by partial least-squares discriminate analysis (PLS-DA). Additionally, the qRT-PCR analysis of the photosynthetic-related genes all showed decreased expression trends along with the anther development. These pieces of evidence indicated that the activities of energy and carbohydrate metabolic pathways were gradually reduced during all the development stages of anther, which affects the sink strength. Overall, this work provides a novel and comprehensive understanding of the metabolic processes in petunia anthers. Keywords: Petunia hybrida; anther development; transcriptome; metabolome
1. Introduction Anther is the key component of plant reproductive organs, for which any abnormal development processes will definitely influence the male plant fertility, and many studies have demonstrated that there are plenty of genes participating in anther development . Petunia is a typical representative of the Solanaceae, which are located on a crucial node in plant phylogeny and have excellent ornamental characteristics and economic value in horticulture. Now, it has been widely used in plant molecular biology studies, as well as in gardens [2,3]. Different from the Arabidopsis thaliana, which has been taken as a model plant in fertility-related research, petunia has larger anthers and could be easily used
Int. J. Mol. Sci. 2018, 19, 955; doi:10.3390/ijms19040955
Int. J. Mol. Sci. 2018, 19, 955
2 of 15
in biochemical and histological analysis . Hence, it has great potential to be taken as an ideal model species of dicots in studying anther development. Recently, the functions of many key genes, including a series of transcription factors, participating in the anther/pollen development have been identified in model and crop plants . However, the gene networks of anther development are very complex, and the understanding of their regulation still needs to be broadened . Transcriptomic sequencing technology could be used to discover a great number of genes’ information, which has been widely utilized in anther transcriptome analysis of many species, such as in rice , wheat , maize , lily  and tobacco . In petunia, the microarray and cDNA-amplified fragment length polymorphism (cDNA-AFLP) have shown that many anther/pollen-specific/predominant genes are expressed in the different development stages of anthers [12–14]. However, due to the limitations of these technologies, the transcriptomic data during petunia anther development are still fragmentary and call for more studies. The gas chromatography-mass spectrometer (GC-MS) technique has been widely employed in metabolic profiling studies of plant tuber growth , fruit maturation , response to environmental changes  and disease [18,19]. The metabolomics of petunia floral scent production have also been well analyzed by this technique, in which benzaldehyde, phenylacetaldehyde, methyl benzoate, phenylethyl alcohol, iso-eugenol and benzyl benzoate were identified as the most abundant benzenoid compounds with the circadian emission rhythms . Recently, two lines of Petunia axillaris (a strongly-scented line and a weakly-scented line) were used to profile the concentrations of metabolites and explore the regulation mechanism of the nocturnal emission rhythm of floral scent concentrations in the scent biosynthetic pathway . However, this technique has not been reported to be used in the anthers, and the metabolic profiling of petunia anther development is still unclear. Omics based on high-throughput technologies have been a useful tool in exploring the complex profile of gene regulatory networks and analyzing the physiological and biochemical changes of metabolic pathways. Recently, many biological phenotypes have been deeply analyzed using the transcriptomic and metabolomic integration strategy, such as the accelerated fruit senescence mechanism after cold storage in litchi , the adaptive responses to ultraviolet radiation of grapevine , the sugar and organic acid metabolism and puffing disorder of Citrus fruit , the copper stress acclimation of brown algae (Ectocarpus siliculosus)  and the transcriptional regulation of the metabolism of sulfate starvation and resupply in Arabidopsis . In order to broaden the knowledge of gene regulation networks and the metabolite profiling during petunia anther development, herein, the transcriptomic and metabolomic of petunia anthers at three critical development stages were analyzed using RNA-seq and GC-MS technologies, in which many significantly differentially-expressed unigenes (DEGs) and metabolic pathways, as well as critical metabolites were identified. This is the first time that comparative multi-omics systems have been used to explore the development process of petunia anthers, and we hope it will be helpful in the following anther development research. 2. Results 2.1. The Chlorophyll Content Is Reduced during the Anther Maturation Process Bulleted in this study, the chlorophyll content was successfully tested in pollen mother cell stage (Stage 1), microsporogenesis stage (Stage 2) and pollen grains stage (Stage 3) of petunia anthers (Figure 1). The content of chlorophyll has a significant decreasing tendency with the development of anthers (Figure 2). When compared with leaf, the concentration of chlorophyll is much less, though it has been proven that chlorophyll indeed exists in anthers, and chlorophyll could play important roles in anther development.
Semi-thin sections of petunia “Fantasy Stage 1, Stage and Stage Figure Figure 1.1. Semi-thin Semi-thin sections sections of of petunia petunia “Fantasy “Fantasy Red” Red” anthers. anthers. Stage Stage 1, 1, Stage Stage 222 and and Stage Stage 33 mean mean for pollen mother cell, microsporogenesis and pollen grain stages, respectively, which the pollen mother cell, microsporogenesis and pollen grain stages, respectively, for which theanthers anthersare are from bud length (without sepal) of 0.5 mm, 15 ±± 0.5 and 35 ±± 0.5 from flower flower buds budswith withaaabud budlength length(without (withoutsepal) sepal)of of333±±±0.5 0.5mm, mm,15 15± 0.5 mm mmand and35 35± 0.5 mm mm individually. 100 μm. individually.Scale Scalebars barsrepresent represent100 100µm. μm.
Figure 2. chlorophyll petunia Figure 2. 2. Contents Contents of of chlorophyll chlorophyll in in petunia petunia anthers anthers and and leaves. leaves. Stage Stage 111 means means the the anthers anthers from from the the Figure Contents of in anthers and leaves. Stage means the anthers from the 33±±0.5 mm flower buds; Stage 2 means the anthers from the 15 ± 0.5 mm flower buds; Stage 3 means means the the anthers anthers from fromthe the15 15± ± 0.5 mm flower buds; buds; Stage Stage 33 means means 3 ± 0.5 mm flower buds; Stage 2 means the from the 35 ±± 0.5 the anthers anthers from from the the 35 35± mm flower flower buds. buds. Different each column column stand stand for for aaa the anthers 0.5 mm mm flower buds. Different capital capital letters letters in in each each column stand for significant difference at 0.01 level. significantdifference differenceat at0.01 0.01level. level. significant
2.2. Sequencing 2.2. Transcriptome TranscriptomeSequencing Sequencing 2.2. Transcriptome Three samples isolated from petunia anthers of different development stages Three total total RNA RNA samples samples were were isolated isolated from from petunia petunia anthers anthers of of different different development development stages stages Three total RNA were −1 (Stage 1, Stage 2, Stage 3). These RNA samples were at a concentration of about 2000 ng·μL with −1 − 1 (Stage 1, 1,Stage Stage2,2,Stage Stage3). 3).These TheseRNA RNAsamples sampleswere wereatataaconcentration concentrationofofabout about2000 2000ngng·μL with (Stage ·µL with OD 260/280 ≥ 1.8, and the RNA integrity numbers (RINs) of 9.1–9.6 were used for cDNA libraries’ OD 260/280 ≥ 1.8, and the RNA integrity numbers (RINs) of 9.1–9.6 were used for cDNA libraries’ OD260/280 ≥ 1.8, and the RNA integrity numbers (RINs) of 9.1–9.6 were used forTMcDNA libraries’ construction. sequenced using the Illumina HiSeq TM 2000 construction. The The three three cDNA cDNA libraries libraries were were sequenced sequenced using using the theIllumina IlluminaHiSeq HiSeqTM 2000 platform. platform. construction. The three cDNA libraries were 2000 platform. A total of about 162 million raw sequencing reads was generated, and after discarding the Atotal totalof of about million raw sequencing reads was generated, and after the discarding the A about 162 162 million raw sequencing reads was generated, and after discarding low-quality low-quality reads, we obtained about 158 million clean reads, 97.53% of the raw reads. For all three low-quality reads, we obtained about 158 million clean reads, 97.53% of the raw reads. For all three reads, we obtained about 158 million clean reads, 97.53% of the raw reads. For all three samples, samples, the quality score above 20 was around 98.50%, and the were samples, the quality score 20 (Q20) (Q20) was98.50%, aroundand 98.50%, andpercentages the GC GC percentages percentages were 43.48%, 43.48%, the quality score above 20 above (Q20) was around the GC were 43.48%, 43.26%, 43.26%, 43.46%, respectively. We achieved 89.7%, 90.2% and 89.8% genome alignment ratios 43.26%, respectively. 43.46%, respectively. We 89.7%, achieved 89.7%, and 89.8% genomeratios alignment ratios of of 43.46%, We achieved 90.2% and 90.2% 89.8% genome alignment of unigenes in unigenes in Stage 1, Stage 2 and Stage 3 of petunia anthers, respectively (Table S1). Approximately unigenes in Stage 1, Stage 2 and Stage 3 of petunia anthers, respectively (Table S1). Approximately Stage 1, Stage 2 and Stage 3 of petunia anthers, respectively (Table S1). Approximately 82% (27,152) of 82% (27,152) of the were detected to be in the three stages development 82%genes (27,152) ofdetected the genes genes detected in tothe be expressed expressed inof the threedevelopment stages of of anther anther development the were towere be expressed three stages anther compared to the compared to the total predicted genes (32,928) of the petunia genome. At last, a total of 4118 putative compared to the total predicted genes (32,928) of the petunia genome. At last, a total of 4118 putative total predicted genes (32,928) of the petunia genome. At last, a total of 4118 putative transcription transcription factors (TFs), belonged to TF families, identified from the transcription factorsbelonged (TFs), which which to 55 55 TF identified families, were were from the database database of of factors (TFs), which to 55belonged TF families, were from identified the database of transcriptomes transcriptomes (Figure 3). transcriptomes (Figure 3). (Figure 3).
2.3. Unigenes between Anther Transcriptomes 2.3.Identify IdentifyDifferentially-Expressed Differentially-ExpressedUnigenes Unigenesbetween betweenAnther AntherTranscriptomes Transcriptomes 2.3. Identify Differentially-Expressed The ofofof anther maturation progress are shown ininFigure 4.4.The TheDEGs DEGsin inpairs pairsofof two stages anther maturation progress areare shown Figure The The DEGs in oftwo twostages stages anther maturation progress shown in Figure 4. transcriptome of the pollen grain Stage 2 in comparison to the microsporogenesis Stage 1 transcriptome of the pollen grain Stage The transcriptome of the pollen grain Stage2 2inincomparison comparisonto to the the microsporogenesis Stage 11 is isis characterized by 6065 DEGs, which 3220 are upregulated unigenes 2845 downregulated. As characterizedby by6065 6065DEGs, DEGs, which 3220 are upregulated unigenes and 2845 downregulated. As characterized ofofof which 3220 are upregulated unigenes andand 2845 downregulated. As the the anthers develop from the microspores of Stage 2 to the mature pollen of Stage 3, 6934 unigenes the anthers develop from the microspores of Stage 2 to the mature pollen of Stage 3, 6934 unigenes anthers develop from the microspores of Stage 2 to the mature pollen of Stage 3, 6934 unigenes are are ofofwhich are and Remarkably, compared areDEGs, DEGs, which 1619upregulated areupregulated upregulated and5315 5315downregulated. downregulated. Remarkably, when compared DEGs, of which 16191619 are and 5315 downregulated. Remarkably, whenwhen compared with with 1,1,3Stage 33 had 3094 unigenes and unigenes. The with Stage Stage3094 had 3094 upregulated upregulated unigenes and 7413 7413 downregulated downregulated unigenes. The Stage 1,Stage Stage had upregulated unigenes and 7413 downregulated unigenes. The comparison comparison between Stage 1 and the other anther developmental stages (2 and 3) implied that the comparison between Stage 1 and the other anther developmental stages (2 and 3) implied that between Stage 1 and the other anther developmental stages (2 and 3) implied that the number the of number unigenes increased with the numberofofdownregulated downregulated unigenessignificantly significantly increased with thepollen pollenformation. formation. downregulated unigenes significantly increased with the pollen formation.
2.4. 2.4.KEGG KEGGPathways PathwaysEnrichment Enrichment 2.4. KEGG Pathways Enrichment According KEGG annotation, anther transcriptomes According to to KEGG KEGG annotation, annotation, the the DEGs DEGs among among the the anther anther transcriptomes transcriptomes were were mainly mainly According to the DEGs among the were mainly assigned 127 pathways. KEGG assigned to 127pathways. pathways. The The upregulated upregulatedand and downregulated downregulatedgenes genes were wereused usedin in the theKEGG KEGG assigned toto127 The upregulated and downregulated genes were used in the pathway enrichment analysis, respectively. Based on the KEGG pathway enrichment analysis of pathway enrichment analysis, respectively. Based on the KEGG pathway enrichment analysis the pathway enrichment analysis, respectively. Based on the KEGG pathway enrichment analysis ofofthe the upregulated expressing genes, the “photosynthesis” and “starch and sucrose metabolism” pathways upregulated expressing genes, the “photosynthesis” and “starch and sucrose metabolism” pathways upregulated expressing genes, the “photosynthesis” and “starch and sucrose metabolism” pathways were significantly different compared weresignificantly significantlydifferent differentas ascompared comparedwithin withineach eachof oftwo twotranscriptomes transcriptomes(Stage (Stage111with withStage Stage2, were as within each of two transcriptomes (Stage with Stage 2,2, Stage Stage11with withStage Stage3,3,Stage Stage22with withStage Stage3). 3).Interestingly, Interestingly,the theKEGG KEGGpathway pathwayenrichment enrichmentanalysis analysisofof
Int. J. Mol. Sci. 2018, 19, 955
5 of 15
Stage 1 with Stage 3, Stage 2 with Stage 3). Interestingly, the KEGG pathway enrichment analysis of the downregulated genes between each of two transcriptomes also obtained the same results. The statistically-enriched pathways between each two transcriptome are shown in Figures S1–S3. 2.5. Validation of the Gene Expression Profiles by qRT-PCR In order to confirm the transcription profile revealed by RNA-seq data, 18 transcription factors, 14 starch and sucrose metabolism pathways and 6 photosynthesis-related genes were selected to design specific primers for qRT-PCR analysis (Figure 5A–C). The overall correlation coefficient (R = 0.89) between RNA-seq and qRT-PCR data was obtained by linear regression analysis, which showed a good correlation between these two data (Figure 5D), indicating that the three transcriptomics data were reliable. 2.6. The Expression Patterns of Photosynthesis-Related Genes In this study, two key enzyme genes RBCL and RBCS, which play important roles in controlling the photosynthesis rate, were used to explore their expression profiles during the different anther development processes. The results showed that their expression levels were very high in the initial development stage of anthers and then decreased gradually with the maturing of anthers (Figure 6A). Then, 11 enzyme genes related to chlorophyll biosynthesis were selected to test their expression patterns. Interestingly, they had an overall downregulated expression trend, which was very similar to RBCL and RBCS (Figure 6A,B). Lastly, the results of the three chlorophyll degradation genes showed that, from the mother cell stage until the microsporogenesis stage, the expression levels of RCCR remained high, and NYC1 increased gradually, while PAO remained low until a sudden increase occurred in the microsporogenesis stage. After that, expressions of all three genes significantly decreased in the 2.5-cm stage and remained low at the pollen grain stage (Figure 6C). 2.7. Metabolome Analysis of Petunia Anther Development by GC-MS To investigate the metabolic profile changes of anthers in the maturation processes, we obtained the GC-MS total ion current (TIC) chromatograms for 18 petunia anther samples from three typical development stages. The obvious differences of chromatographic peaks were observed between sample groups, and the retention times were fairly consistent and reproducible (Figure 7A). In this research, a total of 37 metabolites out of 126 total peaks could be identified in our sample libraries across all samples, and most were either amino acids or carbohydrates. 2.8. Metabolic Profiling and Core Differential Metabolites’ Identification These identified metabolites were visualized in HeatMap, in which they were generally classified into four clusters according to their expression data (Figure 7B). Clusters I, II, III and IV featured the high enrichment of metabolites at anther Stage(s) 1/and 2, only 2, 2 and 3 and only 4, respectively. In Stage 1, the most abundant metabolites were organic acids (ethanedioic acid, hexadecanoic acid, octadecanoic acid, butanoic acid and butanedioic acid) and alcohols (myo-inositol and glycerol), and the abundance of amino acids and sugars (fructose, glucose and galactose) was very low. The greatest abundance of sugars (fructose, glucose, galactose and sucrose) was in Stage 2, then their contents were significantly reduced in Stage 3; while other organic acids and nearly all of the amino acids had the greatest enrichment in Stage 3 (Figure 7B).
Int. J. Mol. Sci. 2018, 19, 955 Int. J. Mol. Sci. 2018, 19, x FOR PEER REVIEW
6 of 15 6 of 15
Figure5.5.The TheqRT-PCR qRT-PCRvalidation validationofofDEGs. DEGs.The Therelative relative expression expression levels levels of of 18 18 transcription transcription factors Figure factors (A), 14 starch and sucrose metabolism pathways (B) and 6 photosynthesis related genes (C). The leftY (A), 14 starch and sucrose metabolism pathways (B) and 6 photosynthesis related genes (C). The left Y axis represents the relative transcript amount obtained by qRT-PCR. The right Y axis represents axis represents the relative transcript amount obtained by qRT-PCR. The right Y axis represents the the fragmentsper per million fragments (FPKM) value of each using RNA-Seq analysis. fragments kbkb perper million fragments (FPKM) value of each genegene using RNA-Seq analysis. ErrorError bars bars indicate the standard errors. (D) Correlation of the gene expression ratios qRT-PCR between indicate the standard errors. (D) Correlation analysisanalysis of the gene expression ratios between qRT-PCR and RNA-seq. and RNA-seq.
Int. J. Mol. Sci. 2018, 19, 955 Int. J. Mol. Sci. 2018, 19, x FOR PEER REVIEW
7 of 15 7 of 15
Figure 6. Expression analyses of the photosynthesis-related genes during petunia anther development. The expression patterns ofanalyses 2 photosynthesis genes (A), 11 chlorophyll genes (B) and Figure 6. Expression of the ratio photosynthesis-related genes biosynthesis during petunia anther 3development. chlorophyll degradation genes patterns (C) in allof the2petunia anther development sample of The expression photosynthesis ratio genes stages. (A), 11The chlorophyll 0.3anther is at genes the pollen stage; thedegradation samples of 0.5anther, 1.5anther at the biosynthesis (B) mother and 3 cell chlorophyll genes (C)1.0anther, in all the petuniaareanther development stages. The sample of 0.3anther is at theare pollen microspore stage; the samples of 2.5anther and 3.5anther at themother pollen cell grainstage; stage.the samples of 0.5anther, 1.0anther, 1.5anther are at the microspore stage; the samples of 2.5anther and 3.5anther are at the pollen grain stage.
Int. J. Mol. Sci. 2018, 19, 955 Int. J. Mol. Sci. 2018, 19, x FOR PEER REVIEW
8 of 15 8 of 15
Figure7.7.GC-MS GC-MSmetabolomic metabolomic analysis analysis of of petunia petunia anthers. (A) GC-MS Figure GC-MS total total ion ion chromatography chromatographyof of anthers under Stage 1, Stage 2 and Stage 3; (B) clustering of the 37 identified metabolites analyzed by anthers under Stage 1, Stage 2 and Stage 3; the 37 identified metabolites analyzed by GC-MS;(C) (C)the thePLS-DA PLS-DAanalysis analysisof ofmetabolites metabolites in in different different anther anther development stages of GC-MS; of petunia. petunia.
To assess systematically the metabolic profile changes during anther development, the PLS-DA To assess systematically the metabolic profile changes during anther development, the PLS-DA plot was generated from the GC-MS metabolite data of Stage 1, Stage 2 and Stage 3 anthers and plot was generated from the GC-MS metabolite data of Stage 1, Stage 2 and Stage 3 anthers and showed showed clear metabolic differences between two stages. Remarkably, Stage 1, Stage 2 and Stage 3 clear metabolic differences between two stages. Remarkably, Stage 1, Stage 2 and Stage 3 anthers anthers could be completely separated sufficiently by use of two principal components. The first could be completely separated sufficiently by use of two principal components. The first principle principle component (PC1) accounting for 48.76% of the variation in the data could separate all three component (PC1) accounting for 48.76% of the variation in the data could separate all three types types of anthers with no outliers, and the second component (PC2) accounting for 47.43% of the of anthers with no outliers, and the second component (PC2) accounting for 47.43% of the variation variation could result in clear separation of Stage 2 anthers from the others (Figure 7C). The could result in clear separation of Stage 2 anthers from the others (Figure 7C). The contribution of each contribution of each variable to PC1 and PC2 was also calculated by giving each variable a weight variable to PC1 and PC2 was also calculated by giving each variable a weight value. The top 10 core value. The top 10 core differential metabolites of PC1 and PC2 discrimination are illustrated in differential metabolites of PC1 and PC2 discrimination are illustrated in Table S2. Interestingly, the two Table S2. Interestingly, the two most important core metabolites were fructose (relative levels of 1.00, most important core metabolites were fructose (relative levels of 1.00, 5.66 and 2.20 in Stage 1, Stage 2 5.66 and 2.20 in Stage 1, Stage 2 and Stage 3 samples, respectively) and glucose (relative levels of and Stage 3 samples, respectively) and glucose (relative levels of 1.00, 6.24 and 1.95 in Stage 1, Stage 2 1.00, 6.24 and 1.95 in Stage 1, Stage 2 and Stage 3 samples, respectively). Both of their contents were significantly increased from Stage 1 to Stage 2 and then obviously decreased in Stage 3 (Figure 7B).
Int. J. Mol. Sci. 2018, 19, 955
9 of 15
and Stage 3 samples, respectively). Both of their contents were significantly increased from Stage 1 to Stage 2 and then obviously decreased in Stage 3 (Figure 7B). 3. Discussion Petunia hybrida, which derived from crosses between two wild parents P. axillaris and P. inflata, is the most popular bedding plant in the world. Furthermore, it has also been taken as a model species for research on biological phenomena. Recently, the high-quality assembled reference genomes of these two wild parents have been accomplished, and this will directly accelerate further research into petunia’s unique characterization . As one of the most important reproductive organs, the development of anther is a complex process, and any disturbing of it could significantly influence the pollen’s fertility, which would bring great losses in breeding. The previous results implied that petunia anther development is a complex process, and some anther-/pollen-specific genes/promoters were obtained [7,12,28,29]. However, due to the different concerns of experiment designs and technical limitation, the knowledge related to petunia anther development still needs to be refined. In this study, three typical development stages of petunia anther, which are the pollen mother cell, microsporogenesis and mature pollen grains stages , were selected to sequence and identify gene expression at the whole genome level using Illumina sequencing technology. It has been known that the basic helix-loop-helix protein (bHLH) transcription factors that participate in various development processes constitute one of the largest families in plants . Here, a total of 55 TF families were identified from the database of the transcriptomes, of which the largest absolute number of TFs (460) was from the members of the bHLH family (Figure 3). These data indicated that the bHLH TFs could have key functions in controlling petunia anther development. When comparing to the abundance of DEGs between two different anther development stages, the number of downregulated genes was significantly increased during anther maturation, and on the contrary, the number of upregulated genes had a decreasing trend (Figure 4), implying that the genes’ expression and the activity of anther in the initial development stage are more vigorous than the later stages. KEGG pathway enrichment analysis was employed to further analyze the biological functions of DEGs. The “photosynthesis” and “starch and sucrose metabolism” pathways relevant to energy metabolism were enriched in both groups (Figures S2–S4). This result suggested that these pathways could play critical roles in petunia anther development. In this research, the whole anthers, which contain a high percentage of biomass from non-pollen tissues, were taken as samples to implement the transcriptome analysis. Nowadays, laser microdissection technology has made it possible to exactly isolated target tissues of pollen. Without the once green anther cuticle, neither “photosynthesis”—nor “starch and sucrose metabolism”—related genes/pathways had been found important in pollen development [32,33]. There are many enzymes participating in chlorophyll biosynthesis, and any disruption of these would seriously restrain the synthesis of chlorophyll, which could significantly influence the photosynthesis rate . In this research, the expression of 11 enzyme genes related to chlorophyll biosynthesis featured sustained declines (Figure 6B), implying that the ability of chlorophyll biosynthesis was gradually weakened during anther development. It has been elucidated that chlorophyll degradation is usually remarkably initiated at the final stage of a plant organ development process, such as fruit ripening and leaf senescence [35,36]. However, it has not been reported in plant reproductive organs, and the molecular mechanism of chlorophyll degradation remains unclear. Through transcriptomic analysis, a series of chlorophyll degradation-related genes was identified in petunia anthers. Dramatically, the expression patterns of three typical chlorophyll degradation pathway-related genes (RCCR, NYC1 and PAOs) were not strictly negative correlated with the chlorophyll biosynthesis-related genes, and all of them were significantly reduced before the microsporogenesis stage (Figure 6C). These results implied that the molecular mechanism of chlorophyll degradation is fairly complex in petunia anthers, and further exploration of this metabolic pathway during petunia anther development is important. Interestingly, chlorophyll was also detected in anthers, and the content decreased with anther growth (Figure 2). The RBCL and RBCS genes, which
Int. J. Mol. Sci. 2018, 19, 955
10 of 15
encode ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO), could control the photosynthesis rate by CO2 fixing . Here, the expression levels of these genes also showed a decreased trend during the maturing processes of anthers (Figure 6A). These results suggested that the photosynthesis pathway existed in petunia anthers and was significantly active in the early development stages. It has been reported that the INVERTASE (INV) gene could catalyze the hydrolysis of sucrose to glucose and fructose, and the HEXOKINASE (HK) and FRUCTOKINASE (FK) genes participate in the glycolysis process of glucose and fructose [38,39]. Here, the expression patterns between INV with HK and FK were completely opposite (Figure S4). This result indicated that the production of glucose and fructose was most active at the beginning of petunia anther development, and then, the glycolysis ratios of glucose and fructose were significantly enhanced, causing the increasing consumptions of glucose and fructose in the following development stages. This conclusion also has been confirmed by the metabolomic analysis, in which the contents of glucose and fructose were firstly increased and then decreased with the development of anthers (Figure 7B). The first component (PC1) and the second component (PC2) of the variation in the data could separate all three stages of anthers with no outliers (Figure 7C), indicating that our metabolomic analysis was reliable and could sufficiently reflect the metabolic profile changes of anthers. A previous study has proven that the sugar signaling pathway might be a central signaling network of cotton anthers, and the conversion of sucrose to fructose and glucose could be suppressed under high temperature condition . Interestingly, the top two core different metabolites during the anther development processes were confirmed as fructose and glucose (Table S2), which were two critical materials in the starch and sucrose metabolism pathways, indicating that this pathway could be critically important in petunia anther development. This finding is in accordance with our above KEGG pathway enrichment result. These evidence suggested that the sink strength was altered in the different development stages of petunia anthers. It is worth mentioning that many peaks detected through GC-MS analysis are still unidentified, and these unknown metabolites would be the subject of our future studies. 4. Materials and Methods 4.1. Plant Materials Plants of Petunia hybrida “Fantasy Red” are a multi- and small-flowered cultivar that was grown in a mix of peat moss, perlite and vermiculite at a volume ratio of 2-1-1 in the experimental greenhouse of Huazhong Agricultural University. Anthers were selected for the pending analysis at immature (pollen mother cell), intermediate (microsporogenesis) and mature (pollen grain) stages, which represent three typical stages of petunia anther development, measured as flower bud lengths (sepals were not included) of 3 ± 0.5 mm (Stage 1), 15 ± 0.5 mm (Stage 2) and 35 ± 0.5 mm (Stage 3), respectively (Figure 1). In order to satisfy the basic sample weight of every experiment, anthers of three plants were mixed with equal amounts and taken as one biological replicate. The fresh samples were taken and immediately frozen in liquid nitrogen and stored at −80 ◦ C until analysis. 4.2. Chlorophylls and Carotenoids Content Assay Briefly, 100 mg of mature leaf and different anther (Stage 1, Stage 2 and Stage 3) tissues of each replicate plant were extracted with 5 mL 96% (v/v) ethanol to determine the content of chlorophyll. Extraction was performed in the dark for 24 h at 4 ◦ C and then analyzed by spectrophotometry using a Hitachi U-2000 spectrophotometer. The absorbance readings were performed at 665 nm and 649 nm for chlorophyll a and chlorophyll b, respectively. Contents of chlorophyll a (Ca ), chlorophyll b (Cb ) and chlorophylls (Pchlorophylls ) were calculated according to the following equation: Ca (mg/mL) = 1.95 × A665 − 6.88 × A649 , Cb (mg/mL) = 24.96 × A649 − 7.32 × A665 and Pchlorophylls (mg/g) = 50 × (Ca + Cb ), where A represents absorbance at the specified wavelength. Every experiment was repeated with three biological replicates.
Int. J. Mol. Sci. 2018, 19, 955
11 of 15
4.3. RNA Extraction Total RNA was extracted from the anther samples above using RNAiso Reagent (Takara, Tokyo, Japan) according to the previously described method . The quality of total RNA samples was assessed using 1% agarose gels. The concentration of RNA samples was tested using NanoDrop (Thermo Scientific, Waltham, MA, USA). The RIN value of RNA was determined by the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). 4.4. RNA-Seq and Reads Mapping Three petunia anther cDNA libraries of three different developmental stages (Stage 1, Stage 2 and Stage 3) were prepared using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. After the purification and fragmentation of RNA, cDNA synthesis, end repair, adapter ligation and PCR amplification, the cDNA library products (ranging from 200–700 bp) were sequenced on the Illumina HiSeqTM 2000 instrument using paired-end sequencing technology by the staff at Beijing Genome Institute (BGI) (Shenzhen, China). The datasets involving three different anther stages were deposited in the NCBI Sequence Read Archive (SRA) with Accession Number SRP126902 under BioProject Number PRJNA422657. Before assembly, the raw reads were first filtered to obtain high-quality clean reads by removing adapter sequences, reads with more than 20% low-quality bases (quality value < 20) or ambiguous nucleotides (denoted with an “N” in the sequence trace). After the purifying process, the clean reads were mapped on the Petunia axillaris reference genome (https://solgenomics.net/organism/Petunia_axillaris/genome), which is more closely related to P. hybrida “Fantasy Red” than other sequenced species Petunia inflata. The SOAP2 software (http://soap.genomics.org.cn/soapaligner.html) was used to map the clean reads under strict default parameters (mismatch number ≤ 5). The reference genes mapped by at least one read, in at least one sample, were defined as expression genes and selected for further analysis. The Pfam database of the petunia anther transcriptome was obtained from the reference annotation of the P. axillaris genome. 4.5. Gene Expression Difference Analysis The unigene expression level was calculated using the FPKM method (fragments per kb per million fragments) . The significance of unigenes’ expression difference was determined by using “FDR ≤ 0.001 and the absolute value of log2 Ratio ≥ 2” as the threshold [42,43]. All the upregulated or downregulated genes were mapped to the KEGG pathway database, and the numbers of unigenes for every KEGG Orthology (KO) term were calculated. Significantly enriched KO terms from the set of upregulated or downregulated genes were identified using the formula of the hypergeometric test , when comparing these upregulated or downregulated genes with the whole petunia anther transcriptome background. 4.6. qRT-PCR Analysis A TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix kit (Transgene, Shenzhen, China) was used to synthesize the first-strand cDNA on the basis of the manufacturer’s instructions. The qRT-PCR was carried out using the SYBR Premix Ex Taq™ IIkit (Takara, Tokyo, Japan) in the ABI 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) according to the protocol described previously . The gene expression level was calculated by the 2−∆∆Ct method, with the beta-actin gene of petunia as the reference gene . Each experiment was repeated with three independent replications, and data are shown as the mean values ± SE (standard error). All primer pairs were designed by Primer Premier 5 software (Premier Biosoft Ltd., Palo Alto, CA, USA) (Table S3), the specificity of which was confirmed by sequencing of qRT-PCR products. The beta-actin gene was used as the internal control.
Int. J. Mol. Sci. 2018, 19, 955
12 of 15
4.7. Metabolite Extraction and Derivatization For metabolite extraction, the petunia anthers (Stage 1, Stage 2 and Stage 3) stored at −80 ◦ C were ground with a mortar and pestle in liquid nitrogen and transferred to 1.5-mL centrifuge tubes. Metabolites were extracted from the samples (ca. 50 mg) in 750 µL 100% methanol, and the samples were vortexed for 30 s. Then, 31.5 µL ribitol (0.2 mg mL−1 ) were added into each sample tube as an internal quantitative standard. The mixed samples were incubated for 15 min at 28 ◦ C in an incubator shaker and then centrifuged for 10 min at 11,000× g. The supernatant (535 µL) of each sample was transferred into a 2-mL centrifuge tube, then 536 µL and 1071 µL ddH2 O (double-distilled H2 O) were added into each 2-mL tube. After being centrifuged (15 min, 2200× g), 400 µL supernatant of each sample were transferred to a new Eppendorf tube and immediately dried in a vacuum. The residue was dissolved in 90 µL 20 mg·mL−1 methoxyamine pyridine solution and incubated for 120 min at 37 ◦ C in an incubator shaker. Finally, the mixture was treated with 60 µL MSTFA reagent (containing 1% TMCS) and incubated for 30 min at 37 ◦ C. Every experiment was repeated by six biological replicates. 4.8. GC-MS Analysis A volume of 1 µL analyte for each sample was absorbed with a split ratio of 10:1 and injected into a DB-5 MS capillary column (30 m × 0.25 mm inner diameter, 0.25 µm film thickness; J&W Scientific, Folsom, CA, USA). The injection and transfer line temperatures were 280 and 270 ◦ C. Helium was selected as the carrier gas, and the gas flow rate was 1 mL·min−1 . The column initial temperature was retained at 70 ◦ C for 5 min and finally increased to 300 ◦ C at a heating rate of 5 ◦ C·min−1 for 3 min. Following the same program of samples, the retention times of C10–C40 alkanes were obtained and used to generate retention indices (RIs) of sample peaks that have similarities of more than 700 with the spectral catalogs in NIST05 (National Institute of Standards). The peak area value of ribitol was taken as the internal standard to calculate the metabolite abundance. The SIMCA-P 11.5 software (Umetrics AB, Umea, Sweden) was used to examine the differences in the metabolite levels of petunia anthers of different development stages by partial least-squares discriminate analysis (PLS-DA). In order to confirm the significant differential metabolites at least between two different stages, the SPSS 13.0 software (SPSS Inc., Chicago, IL, USA) was used for variance analysis (p < 0.05). 4.9. Statistical Analysis Statistical analysis was performed using the SAS 8.1 software (SAS Institute, Cary, NC, USA) followed by Tukey’s multiple range tests (p < 0.01). 5. Conclusions Based on the sequence data and mass chromatographic peaks obtained by GC-MS and RNA-Seq analysis at three typical development stages of petunia anthers, many critical unigenes and metabolites were identified, which showed that anther development is a complex process. The unigenes’ expression profiles and metabolite content changes through three petunia development stages were analyzed further. When combining the transcriptomic with GC-MS metabolomic results, the “starch and sucrose metabolism” pathway was indicated as the most important metabolic pathway in the petunia anther ripening processes. We found that the chlorophyll content, as well as the expression levels of the photosynthesis rate and chlorophyll biosynthetic genes all showed a decreasing trend in petunia anther ripening. The content of the two most important core different metabolites (fructose and glucose) all reached the peak in the microsporogenesis stage. These pieces of evidence suggested that the anther could be taken as a metabolic source, as well as a strong metabolic sink in the pollen mother cell stage, then it became only a weak sink in the following development stages, and the photosynthesis-related pathways could play important roles in the changes of the sink strength. This study offers new insights
Int. J. Mol. Sci. 2018, 19, 955
13 of 15
into the molecular and metabolic mechanisms in petunia anther development characteristics and could provide excellent platforms and study contents for future reproductive organ development research. Supplementary Materials: Supplementary materials can be found at http://www.mdpi.com/1422-0067/19/4/955/s1. Acknowledgments: This research work was funded by the National Natural Science Foundation of China (Grant No. 31601785) and the Fundamental Research Funds for the Central Universities (Projects 2662016PY042, 2662017PY075). Author Contributions: Yuanzheng Yue and Huirong Hu designed the study. Yuanzheng Yue performed most of the experiments and drafted the manuscript. Shaoze Tian and Yu Wang reviewed and edited the article. Hui Ma and Siyu Liu arranged the data. Yuanzheng Yue and Yuqiao Wang arranged the figures and tables. Conflicts of Interest: The authors declare no conflict of interest.
References 1. 2. 3.
11. 12. 13.
14. 15. 16.
McCormick, S. Molecular analysis of male gametogenesis in plants. Trends Genet. 1991, 7, 298–303. [CrossRef] Gerats, T.; Vandenbussche, M. A model system for comparative research: Petunia. Trends Plant Sci. 2005, 10, 251–256. [CrossRef] [PubMed] Zenoni, S.; D’Agostino, N.; Tornielli, G.B.; Quattrocchio, F.; Chiusano, M.L.; Koes, R.; Zethof, J.; Guzzo, F.; Delledonne, M.; Frusciante, L.; et al. Revealing impaired pathways in the an11 mutant by high-throughput characterization of Petunia axillaris and Petunia inflata transcriptomes. Plant J. 2011, 68, 11–27. [CrossRef] [PubMed] Vandenbussche, M.; Chambrier, P.; Rodrigues Bento, S.; Morel, P. Petunia, Your Next Supermodel? Front. Plant Sci. 2016, 7, 72. [CrossRef] [PubMed] Yue, L.; Twell, D.; Kuang, Y.F.; Liao, J.P.; Zhou, X.Q. Transcriptome analysis of Hamelia patens (Rubiaceae) anthers reveals candidate genes for tapetum and pollen wall development. Front. Plant Sci. 2016, 7, 1991. [CrossRef] [PubMed] Deveshwar, P.; Bovill, W.D.; Sharma, R.; Able, J.A.; Kapoor, S. Analysis of anther transcriptomes to identify genes contributing to meiosis and male gametophyte development in rice. BMC Plant Biol. 2011, 11, 78. [CrossRef] [PubMed] Hobo, T.; Suwabe, K.; Aya, K.; Suzuki, G.; Yano, K.; Ishimizu, T.; Fujita, M.; Kikuchi, S.; Hamada, K.; Miyano, M.; et al. Various spatiotemporal expression profiles of anther-expressed genes in rice. Plant Cell Physiol. 2008, 49, 1417–1428. [CrossRef] [PubMed] Golkari, S.; Gilbert, J.; Prashar, S.; Procunier, J.D. Microarray analysis of Fusarium graminearum-induced wheat genes: Identification of organ-specific and differentially expressed genes. Plant Biotechnol. J. 2007, 5, 38–49. [CrossRef] [PubMed] Ma, J.; Duncan, D.; Morrow, D.J.; Fernandes, J.; Walbot, V. Transcriptome profiling of maize anthers using genetic ablation to analyze pre-meiotic and tapetal cell types. Plant J. 2007, 50, 637–648. [CrossRef] [PubMed] Hsu, Y.F.; Tzeng, J.D.; Liu, M.C.; Yei, F.L.; Chung, M.C.; Wang, C.S. Identification of anther-specific/predominant genes regulated by gibberellin during development of lily anthers. J. Plant Physiol. 2008, 165, 553–563. [CrossRef] [PubMed] Bokvaj, P.; Hafidh, S.; Honys, D. Transcriptome profiling of male gametophyte development in Nicotiana tabacum. Genomics Data 2015, 3, 106–111. [CrossRef] [PubMed] Cnudde, F.; Moretti, C.; Porceddu, A.; Pezzotti, M.; Gerats, T. Transcript profiling on developing Petunia hybrida floral organs. Sex. Plant Reprod. 2003, 16, 77–85. [CrossRef] Cnudde, F.L.; Hedatale, V.; de Jong, H.; Pierson, E.S.; Rainey, D.Y.; Zabeau, M.; Weterings, K.; Gerats, T.; Peters, J.L. Changes in gene expression during male meiosis in Petunia hybrida. Chromosome Res. 2006, 14, 919–932. [CrossRef] [PubMed] Shimamura, K.; Ishimizu, T.; Nishimura, K.; Matsubara, K.; Kodama, H.; Watanabe, H.; Hase, S.; Ando, T. Analysis of expressed sequence tags from Petunia flowers. Plant Sci. 2007, 173, 495–500. [CrossRef] Roessner, U.; Wagner, C.; Kopka, J.; Trethewey, R.N.; Willmitzer, L. Simultaneous analysis of metabolites in potato tuber by gas chromatography–mass spectrometry. Plant J. 2000, 23, 131–142. [CrossRef] [PubMed] Smith, E.D.; Whiting, M.D.; Rudell, D.R. Metabolic profiling of ethephon-treated sweet cherry (Prunus avium L.). Metabolomics 2011, 7, 126–133. [CrossRef]
Int. J. Mol. Sci. 2018, 19, 955
33. 34. 35.
14 of 15
Jänkänpää, H.J.; Mishra, Y.; Schröder, W.P.; Jansson, S. Metabolic profiling reveals metabolic shifts in Arabidopsis plants grown under different light conditions. Plant Cell Environ. 2012, 35, 1824–1836. [CrossRef] [PubMed] Gai, Y.P.; Han, X.J.; Li, Y.Q.; Yuan, C.Z.; Mo, Y.Y.; Guo, F.Y.; Liu, Q.X.; Ji, X.L. Metabolomic analysis reveals the potential metabolites and pathogenesis involved in mulberry yellow dwarf disease. Plant Cell Environ. 2014, 37, 1474–1490. [CrossRef] [PubMed] Cevalloscevallos, J.M.; Futch, D.B.; Shilts, T.; Folimonova, S.Y.; Reyesdecorcuera, J.I. GC-MS metabolomic differentiation of selected citrus varieties with different sensitivity to citrus huanglongbing. Plant Physiol. Biochem. 2012, 53, 69–76. [CrossRef] [PubMed] Verdonk, J.C.; Ch, R.D.V.; Verhoeven, H.A.; Haring, M.A.; van Tunen, A.J.; Schuurink, R.C. Regulation of floral scent production in petunia revealed by targeted metabolomics. Phytochemistry 2003, 62, 997–1008. [CrossRef] Oyama-Okubo, N.; Sakai, T.; Ando, T.; Nakayama, M.; Soga, T. Metabolome profiling of floral scent production in Petunia axillaris. Phytochemistry 2013, 90, 37. [CrossRef] [PubMed] Yun, Z.; Qu, H.X.; Wang, H.; Zhu, F.; Zhang, Z.K.; Duan, X.W.; Yang, B.; Cheng, Y.J.; Jiang, Y.M. Comparative transcriptome and metabolome provides new insights into the regulatory mechanisms of accelerated senescence in litchi fruit after cold storage. Sci. Rep. 2016, 6, 19356. [CrossRef] [PubMed] Matus, J.T. Transcriptomic and metabolomic networks in the grape berry illustrate that it takes more than flavonoids to fight against ultraviolet radiation. Front. Plant Sci. 2016, 7, 1337. [CrossRef] [PubMed] Ibanez, A.M.; Martinelli, F.; Reagan, R.L.; Uratsu, S.L.; Vo, A.; Tinoco, M.A.; Phu, M.L.; Chen, Y.; Rocke, D.M.; Dandekar, A.M. Transcriptome and metabolome analysis of citrus fruit to elucidate puffing disorder. Plant Sci. 2014, 218, 87–98. [CrossRef] [PubMed] Ritter, A.; Dittami, S.M.; Goulitquer, S.; Correa, J.A.; Boyen, C.; Potin, P.; Tonon, T. Transcriptomic and metabolomic analysis of copper stress acclimation in Ectocarpus siliculosus highlights signaling and tolerance mechanisms in brown algae. BMC Plant Biol. 2014, 14, 116. [CrossRef] [PubMed] Bielecka, M.; Watanabe, M.; Morcuende, R.; Scheible, W.R.; Hawkesford, M.J.; Hesse, H.; Hoefgen, R. Transcriptome and metabolome analysis of plant sulfate starvation and resupply provides novel information on transcriptional regulation of metabolism associated with sulfur, nitrogen and phosphorus nutritional responses in Arabidopsis. Front. Plant Sci. 2015, 5, 805. [CrossRef] [PubMed] Bombarely, A.; Moser, M.; Amrad, A.; Bao, M.Z.; Bapaume, L.; Barry, C.S.; Bliek, M.; Boersma, M.R.; Borghi, L.; Bruggmann, R.; et al. Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida. Nat. Plants 2016, 2, 16074. [CrossRef] [PubMed] Yue, Y.Z.; Yin, C.Q.; Guo, R.; Peng, H.; Yang, Z.N.; Liu, G.F.; Bao, M.Z.; Hu, H.R. An anther-specific gene PhGRP is regulated by PhMYC2 and causes male sterility when overexpressed in petunia anthers. Plant Cell Rep. 2017, 36, 1401–1415. [CrossRef] [PubMed] Yue, Y.Z.; Yang, D.; Sun, J.; Peng, H.; Yin, C.Q.; Guo, R.; Ning, G.G.; Hu, H.R. A novel PhLRR gene promoter is sufficient for engineering male sterility in petunia. Plant Mol. Biol. Rep. 2016, 34, 970–977. [CrossRef] Porceddu, A.; Reale, L.; Lanfaloni, L.; Moretti, C.; Sorbolini, S.; Tedeschini, E.; Ferranti, F.; Pezzotti, M. Cloning and expression analysis of a Petunia hybrida flower specific mitotic-like cyclin. FEBS Lett. 1999, 462, 211–215. [CrossRef] Zhang, J.; Wu, Y.; Li, D.; Wang, G.; Li, X.; Xia, Y.P. Transcriptomic analysis of the underground renewal buds during dormancy transition and release in “Hangbaishao” peony (Paeonia lactiflora). PLoS ONE 2015, 10, e0119118. [CrossRef] [PubMed] Tang, X.; Zhang, Z.Y.; Zhang, W.J.; Zhao, X.M.; Li, X.; Zhang, D.; Liu, Q.Q.; Tang, W.H. Global gene profiling of laser-captured pollen mother cells indicates molecular pathways and gene subfamilies involved in rice meiosis. Plant Physiol. 2010, 154, 1855–1870. [CrossRef] [PubMed] Yuan, T.L.; Huang, W.J.; He, J.; Zhang, D.; Tang, W.H. Stage-specific gene profiling of germinal cells helps delineate the mitosis/meiosis transition. Plant Physiol. 2018, 176, 1610–1626. [CrossRef] [PubMed] Eckhardt, U.; Grimm, B.; Hortensteiner, S. Recent advances in chlorophyll biosynthesis and breakdown in higher plants. Plant Mol. Biol. 2004, 56, 1–14. [CrossRef] [PubMed] Schelbert, S.; Aubry, S.; Burla, B.; Agne, B.; Kessler, F.; Krupinska, K.; Hörtensteiner, S. Pheophytin pheophorbide hydrolase (Pheophytinase) is involved in chlorophyll breakdown during leaf senescence in Arabidopsis. Plant Cell 2009, 21, 767–785. [CrossRef] [PubMed]