Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit. Alphonso is known as the “King of mangos” due to its unique flavor, attractive color, low fiber pulp and long shelf life. We analyzed the transcriptome of Alphonso mango through Illumina sequencing from seven stages of fruit development and ripening as well as flower.

1. Changes During Fruit Ripening
2. Stages Of Fruit Ripening

Total transcriptome data from these stages ranged between 65 and 143 Mb. Importantly, 20,755 unique transcripts were annotated and 4,611 were assigned enzyme commission numbers, which encoded 142 biological pathways. These included ethylene and flavor related secondary metabolite biosynthesis pathways, as well as those involved in metabolism of starch, sucrose, amino acids and fatty acids. Differential regulation (p-value ≤ 0.05) of thousands of transcripts was evident in various stages of fruit development and ripening. Novel transcripts for biosynthesis of mono-terpenes, sesqui-terpenes, di-terpenes, lactones and furanones involved in flavor formation were identified. Large number of transcripts encoding cell wall modifying enzymes was found to be steady in their expression, while few were differentially regulated through these stages. Novel 79 transcripts of inhibitors of cell wall modifying enzymes were simultaneously detected throughout Alphonso fruit development and ripening, suggesting controlled activity of these enzymes involved in fruit softening.

Introduction Mango ( Mangifera indica L.) is one of the most popular and highly favored fruit. Global mango production was reported to be 43.3 million metric tons in 2013 preceding banana, apple, grape and orange. There are thousands of mango cultivars worldwide, among which Alphonso, Keitt, Kent, Lilli, Zill, Osteen, Haden, Kesar, Pairi, Dashehari, Langra and Banganapalli are well known. These varieties vary in their fruit color, size, shape, flavor, taste and ripening period and pattern. To understand the composition and biosynthesis of their unique flavor and complex ripening process various studies have been carried out at metabolic –, proteomic –, genetic – and post-harvest processing – levels. Whole genome sequencing and RNA sequencing (RNAseq) are the two important high throughput technologies recently adopted to understand the complex cellular and physiological processes in fruits such as citrus, tomato and strawberry as well as domestication and diseases tolerance in citrus,. Although mango genome sequence is yet not available, few recent studies have described the transcriptomic analysis from various tissues of few mango cultivars.

The first report from Zill mango provided extensive transcriptomic and proteomic profiling from pulp and skin tissues of four fruit developing stages using pooled RNA but not stage specific and differentially expressed transcripts. Another study of leaf transcriptome and chloroplast genome sequencing from cultivar Langra provided information about the production of several bioactive compounds.

Transcriptome analysis from two (raw and ripe) and three (raw, mid ripe and ripe) stages of Kent and Dashehari fruit pulp, respectively provided important insights into the ripening process and flavor biogenesis in these mango cultivars. India is the largest producer and exporter of mango with 40.6% share in international mango market. Among the Indian mango cultivars Alphonso is globally favored and highly exported mango due to its unique and attractive flavor, low fiber containing pulp and high carotene content,. The ripening duration of Alphonso mango is 15 days from harvest, which is the highest among all mango cultivars; for example, the ripening duration for Kent and Dashehari mango fruit is 10 and 6 days, respectively,. Fruit ripening in Alphonso mango progresses from skin towards the stone leading to attractive skin color and easy monitoring of ripening progress.

On the other hand various mango varieties, viz. Haden, Keitt, Kent, Tommy Atkins (National Mango Board, USA; ) and Dashehari show polarity of their ripening from fruit stone to skin, making it difficult to identify ripened fruits. Longer ripening duration and shelf life of Alphonso mango provides sufficient time for its transportation globally. The mechanisms underlying these unique properties of Alphonso mango need to be explored in depth at spatial and temporal level of fruit development and ripening using transcriptomic, proteomic and metabolomic approaches as they can be correlated to the specific phenotype.

In the present study, we performed transcriptome analysis of Alphonso mango through eight different tissues such as flower, whole fruit at 30 and 60 DAP (Days After Pollination), pulp and skin of 90 DAP fruit (mature raw fruit) and pulp from three fruit ripening stages i.e. 5 DAH (Days After Harvest): table green stage; 10 DAH: mid ripe stage and 15 DAH: ripe stage to analyze various fruit developing and ripening processes in Alphonso mango. Alphonso mango transcriptome Alphonso mango transcriptome was screened through eight tissue samples. To map differentially expressed transcripts a merged assembly was generated from the reads of all the tissues, which reflected upon overall Alphonso mango flower and fruit transcriptome.

For each tissue read numbers were more than 100 million, which were assembled using k-mers 67, 75 and 83 separately and then merged for individual tissue. Average number of unique transcripts post assembly was 76,043 and with N50 and N80 values as 1,835 and 1,008 bp, respectively (Table ). The minimum and maximum lengths of transcript from these assemblies were 100 and 17,342 bp, respectively (Table ). Average number of transcripts was 11,925 upon filtering for redundancy and identifying candidate coding regions with maximum 90% identity and minimum 70% coverage (Table ). Number of non redundant (NR) transcripts/sample. Unique transcripts from each assembled and filtered set were subjected to BLASTx against the non-redundant dataset from NCBI.

From a total of 20,755 unique transcripts from the merged assembly, 92.22% transcripts were annotated while 954 transcripts (4.59%) encoded hypothetical proteins and 661 (3.17%) remained unidentified. BLASTx statistics revealed maximum hits from Citrus sinensis and C. Clementina followed by Theobroma cacao, Jatropha curcas, Vitis venifera and Ricinus communis (Fig. ). Total 74,330 GO terms related to various biological processes (BP), molecular functions (MF) and cellular components (CC) were assigned to these 20,755 unique transcripts. Metabolic, cellular and single organism processes followed by biological regulation and localization were the most abundant terms under the BP category; binding and catalytic activity related terms under MF category; while, organelle, membrane and macromolecular complex terms under the CC category, respectively (Fig., Supplementary Figure ).

Blast statistics showing ( a) distribution of top hit species, ( b) distribution of top gene ontologies from BP: biological processes, MF: molecular functions and CC: cellular components and ( c) Enzyme code distribution: number of transcripts (Y-axis) encoding. Unique transcripts were assigned enzyme commission (EC) number to determine the involvement of these transcripts in various BPs. Total 4,611 ECs were assigned from oxidoreductase, transferase, hydrolase, ligase, lyase and isomerase classes; wherein, transferases were the most abundant followed by hydrolases in all the eight tissues (Fig. ). These assigned ECs represented 142 known pathways from the KEGG database –, which are potentially functional in Alphonso mango fruit development and ripening. Most of these pathways were saturated with higher number of annotations from transcriptome data e.g.

Metabolism of starch, sucrose and various amino acids including methionine and biosynthesis pathway of ethylene, phenylpropanoids and flavonoids (Supplementary Figure ). Differentially expressed transcripts Comparison between adjacent tissue stages was carried out to identify differentially expressed transcripts at each stage of fruit development and ripening (Table, Supplementary Table ). Flower to fruit of 30DAP, 524 transcripts were down regulated while 181 were up-regulated.

Among the down-regulated transcripts alpha-amylase and subtilisin inhibitor-like (contig6593), carbonic anhydrase (contig5377), chitinase (contig4907) and maternal effect embryo arrest 59 (contig4949) exhibited the highest fold change (2-fold); while homeodomain-like protein (contig12912 and 12913), cytochrome p450 - cyp72a219-like (contig4083), heat shock cognate 70 kDa protein (contig5948) and inositol-3-phosphate synthase (contig9620) were highly up-regulated (2-fold). Number of differentially expressed transcripts. Transition from 30 DAP to 60 DAP resulted in up-regulation of 73 and down-regulation of seven transcripts. Important down-regulated transcripts were, n-acetyltransferase (contig1261), 9- cis-epoxycarotenoid dioxygenase (contig1680), protein reversion-to-ethylene sensitivity (contig9164) and ethylene receptor 2 (contig6147).

While important up-regulated transcripts were nucleotide sugar transporter family protein (contig6337), beta-xylosyltransferase (contig4150), various cellulose synthase catalytic subunits and laccases. Comparison of 60 DAP fruit tissue with 90 DAP pulp and skin tissue, respectively revealed down-regulation of beta-xylosyltransferase, beta-d-xylosidase, cellulose synthase and galacturonosyl transferase; whereas spx and exs domain-containing protein was up-regulated in both 90 DAP pulp and skin.

Up-regulation of homeobox protein sbh1, gdsl esterase lipase, caffeoyl shikimate esterase, hydroperoxide lyase, udp-rhamnose:rhamnosyl transferase and pectinesterase inhibitor was evinced from 90 DAP pulp tissue compared to that in 60 DAP fruit tissue. Evaluation of differentially expressed transcripts between pulp and skin of 90 DAP revealed very few transcripts down-regulated in skin with none of them showing change more than 2-fold, while 54 transcripts were up-regulated in 90 DAP skin tissue including the important ones as DNA mismatch repair protein msh5, acyl carrier protein, amino-acid permease, calcium-transporting ATPase, isoflavone reductase and various disease resistance proteins. During ripening of Alphonso pulp from 90 DAP to 5 DAH, 42 transcripts were up-regulated (2-fold) including methyltransferase, amino-acid permease, chloroplastic 9- cis-epoxycarotenoid dioxygenase, beta-galactosidase and protein phosphatase mainly, whereas 52 transcripts were down-regulated, important ones being few disease resistance proteins, peroxidase, sucrose synthase and glycerol-3-phosphate dehydrogenase. The highest level of differential expression was evinced through transition from 5 DAH to 10 DAH among all the ripening tissues wherein 418 transcripts were down and 191 were up-regulated. These stages are known for non-climacteric to climacteric transition during Alphonso fruit ripening based on metabolite analysis,. Prominently down regulated transcripts were phospholipase-A LCAT3, sugar phosphate exchanger, auxin-responsive protein IAA9, abscisate β-glucosyltransferase and membrane-associated kinase regulator with more than 3-fold change, while up-regulated transcripts were aspartic proteinase nepenthesin, UDP-glucose 6-dehydrogenase, 1-aminocyclopropane-1-carboxylate oxidase, peroxidase, bidirectional sugar transporter sweet1-like, omega-6 fatty acid desaturase and squalene monooxygenase. During the transition from 10 DAH to 15 DAH stage, increased phosphate metabolism was evident.

Total 30 transcripts were up-regulated mainly including phospholipase-D, inorganic pyrophosphatase, phosphate transporter, transcription factor glk2 and DNA translocase, whereas 12 transcripts were down-regulated and important ones were aspartic protease, plastocyanin-like domain protein, annexin d2-like and glucuronosyltransferase (2-fold). Idiosyncratic transcripts for the stage Various transcripts unique to each of these stages (Fig., Table, Supplementary Table ) representing multiple stage specific processes were identified during this comparative analysis. A total of 388 transcripts were identified as unique to flower tissue, which mostly included various transcription and translation factors, late embryogenesis abundant proteins, stress-sensitive and dehydration-responsive element-binding protein-1b and various ribosomal proteins. Distinct transcripts.

Similarly, distinct transcripts specific to fruit developing and ripening stages were identified, wherein 1,090 and 27 transcripts were idiosyncratic to the fruit developing and ripening stages, respectively (Supplementary Table ). Various auxin and gibberellin induced and regulated proteins as well as several proteins responsible for various vacuole activities, multiple disease resistance proteins and various terpene synthases were distinct to the developing stages. Transcripts encoding for multiple transcription and translation factors during Alphonso fruit development were also detected. Interestingly multiple ethylene responsive transcription factors along with the protein reversion-to-ethylene sensitivity were uniquely revealed in the developing stages.

Whereas, NIN-like protein, respiratory burst oxidase homolog protein d-like, rhamnogalacturonate lyase b-like, lectin receptor kinase, gag protein and methionine Υ-lyase were exclusive to the Alphonso ripening stages. Similarly, WRKY transcription factor 43, b3 domain-containing val3 and ap2 ERF domain-containing transcription factors were uniquely identified from the ripening stages. Several hypothetical and uncharacterized proteins were also found to be Alphonso ripening specific and their characterization might help to reveal the ripening process in Alphonso. Gene ontology (GO) enrichment Fisher’s exact test was performed to understand over- and down-represented GOs (p-value ≤ 0.001) during the transition, which gave an overall picture of the Alphonso mango development and ripening (Fig. ).

During the flower to 30 DAP fruit transition certain GOs were over-expressed such as post-embryonic developmental and anatomical structure morphogenesis process, response to abiotic stimulus and various plastid and thylakoid processes. At the same time, hydrolysis of O-glycosyl compounds and UDP-glycosyltransferase activity along with the methyltransferase, sucrose metabolic and lipid biosynthetic activities coding GOs were down-represented. 30 DAP to 60 DAP transition described enriched GOs for starch and sucrose metabolic process, fatty acid, cellulose and vitamin biosynthetic processes, protein glycosylation, response to oxidative stress and heat along with HSP binding. During the same event cell differentiation, growth and cell death in addition to the anatomical structure morphogenesis, secondary metabolic process and response to biotic stimulus coding GOs were observed to be decreased.

Comparison between 60 DAP fruit and 90 DAP pulp revealed over-represented GOs for cell differentiation and cell growth, response to biotic stimulus, secondary metabolic process, tropism and regulation of gene expression. While, starch and sucrose metabolic processes, O-glycosyl hydrolase UDP-glucosyltransferase activities, fatty acid biosynthesis and catabolic processes along with the lipid transport, various amino acid and terpenoid biosynthesis and cell wall biogenesis were down-represented. 90 DAP pulp and skin tissues showed over-expression of various ion, ADP and sequence-specific DNA binding activities along with the response to hormone, chemical and organic substances in the skin tissue; whereas only cytoplasm and cytoplasmic part related GOs were found to be down-expressed in the skin compared to that in the pulp at 90 DAP. Over and down expressed gene ontologies (GO) between stages of development and ripening. Transition of 90 DAP pulp to 5 DAH pulp revealed only over-representation of GOs such as biosynthesis of RNA, aromatic compounds, and nucleobase containing compounds; binding of various ions; regulation of primary, cellular, nitrogenous compounds, RNA biosynthesis processes and gene expression and response to hormone, chemical and organic substance. Surprisingly, none of the GOs showed significant (P-value ≤ 0.001) enrichment during the transition from 5 DAH to 10 DAH and reflected as the stationary phase of Alphonso mango ripening.

The 10 DAH (mid ripe stage) to 15 DAH (ripe stage) transition showed over-expression of many GOs such as binding of various ions; phospholipid, glycolipid, amino acid, co-enzyme and water soluble vitamin biosynthetic processes; methyltransferase, UDP-glycosyltransferase and O-glycosyl hydrolase activities; starch, sucrose metabolic processes, gluconeogenesis and glycolytic process; protein phosphorylation, dephosphorylation, folding, protein modification and small protein conjugation along with the regulation and response to oxidative stress, oxidoreductase activity and redox homeostasis. While the GOs related to plastid and cytosol; post embryonic development; response to abiotic and endogenous stimulus; anatomical structure morphogenesis, plasma membrane, cell differentiation and communication were down-represented. Spatial changes in transcriptome at 90 DAP In case of Alphonso mango, 90 DAP stage is a mature raw stage of the fruit and is considered as the right stage of fruit harvest (0 DAH) for further artificial ripening of the fruit.

Hence transcriptome analysis of skin and pulp were separately carried out at this stage. Overall 90 DAP skin was found to be metabolically more active compared to the 90 DAP pulp with respect to differentially expressed genes (Supplementary Table ), unique genes (Supplementary Table ) and enriched GOs (Fig. ). In the skin, 54 transcripts were up-regulated whereas only four were down-regulated compared to the pulp. Among the up-regulated transcripts (2-fold), isoflavone reductase, transcription and translation regulatory proteins, hydrolases, methyltransferase etc. Were more prominent; whereas four down-regulated transcripts were membrane and cytoplasm related GOs. Unique transcripts upon comparison between 90 DAP pulp and skin showed carotene and xanthophylls biosynthesis related contigs, beta-carotene hydroxylase (contig5406) and anthocyanidin 3-o-glucosyltransferase, flavor related various terpene synthases and ripening related contigs such as ethylene-responsive transcription factors, pectate lyase, pectin esterase and cellulase as unique to skin compared to the pulp.

Interestingly, in comparison of 5DAH stage with 90 DAP pulp few of these transcripts such as isoflavone reductase (contig5337) and beta-carotene hydroxylase (contig5406) were observed to be idiosyncratic to 5 DAH stage as seen in 90 DAP skin stage. Similarly, transcripts encoding enzymes involved in pectin degradation (contigs9578, 9579, 18096 and 3603) were found to be upregulated (2-fold) in the comparison between 10DAH vs 90 DAP pulp and 15 DAH vs 90 DAP pulp. These finding suggest progression of ripening related molecular processes in the pulp of 5 DAH and onward stages which were initiated at 90 DAP skin and highlight initiation of Alphonso ripening process from skin and its probable progress towards fruit stone. Genes involved in the flavor biogenesis in Alphonso mango Quantitatively mono-terpenes are abundant in Alphonso followed by sesqui-terpenes,. The present data revealed six contigs encoding mono-terpene synthases ( limonene synthase1, limonene synthase2, beta ocimine synthase1, beta ocimine synthase2, isoprene synthase1 and isoprene synthase2), five contigs encoding sesquiterpene synthases ( germacreneD synthase1, germacreneD synthase2, nerolidol synthase1, nerolidol synthase2 and alpha farnesene synthase) and three contigs coding for di-terpene synthases ( ent-kaurene synthase and casbene synthase (E,E)-geranillinalool synthase).

Changes During Fruit Ripening

Phylogenetic analysis (Fig. ) of these genes along with other plant terpene synthases (TPS) showed distribution of these genes in to the TPS-a, TPS-b, TPS-e and TPS-f clades, respectively. Cladogram representing phylogenetic analysis of proteins encoded by various terpene synthases identified in the present study (blue color) along with the terpene synthases from other Angiosperm plants (black color) using the neighbor-joining method. Furaneol and mesifuran are the two furanones from Alphonso mango and their synthesis by enone oxidoreductase ( EO) and O-methyltransferase ( OMTS), respectively have been described earlier,. Multiple contigs coding for quinone oxidoreductase and O-methyltransferases were detected in the present analysis. Phylogenetic analysis of these contigs with the characterized genes revealed another similar transcript variant for the MiEO (Supplementary Figure ), whereas none of the contigs showed similarity to the MiOMTS (Supplementary Figure ). Green grassy aroma of unripe fruits is due to the C6 volatiles formed during the lipoxygenase and hydroperoxide lyase (HPL) pathways. In the present study a single contig encoding the hydroperoxide lyase and six contigs encoding 13-lipoxygenase were detected.

Involvement of 9-lipoxygenase ( Mi9LOX) and epoxide hydrolase 2 ( MiEH2) in the biogenesis of lactones from Alphonso mango was confirmed in our recent study. One more transcript encoding 9-lipoxygenase similar to that of the characterized Mi9LOX (Supplementary Figure ) and three contigs coding for epoxide hydrolase 2 grouping with the MiEH2 (Supplementary Figure ) were additionally detected in the present study. Three more novel contigs coding for epoxide hydrolase were also detected but neither grouped with the other EH1or EH2 from different plant species (Supplementary Figure ). We further analyzed the differential expression of all these flavor related genes in terms of transcript abundance in flower and fruit developing and ripening stages (Supplementary Figure ). Among these the contigs encoding various TPS were abundant in flower and during early developing stages.

Interestingly, many contigs encoding EH (contigs 8280, 3904, 3901, 14123 and 8281), LOX (contigs 18105, 12748, 12747 and 12385) and EO (contigs 8026, 5618, 15594, 14137 and 13600) were found to be ripening specific and might be playing crucial role in generating unique aroma volatiles during Alphonso fruit ripening. Glycosidases and cell wall degrading enzymes from Alphonso mango Glycosidases are involved in a variety of functions such as hydrolysis of complex carbohydrates (storage and structural) to mono-saccharides, removal of sugars from various glycans including glycosidically bound aroma volatiles, which serve as storage pool for the aroma compounds. In the present study, many glycosidases acting on various sugars i.e. Glucose, galactose, mannose, fructose, xylose, fucose and rhamnose, were detected. Among these, class glucosidase with the highest number of contigs (51 contigs) was observed to contain 28 and 21contigs coding for glucan β-glucosidase and general β-glucosidase, respectively.

Two contigs coding for α-glucosidases were identified of which one coded for glucan α-glucosidase and the other for general α-glucosidase. Among these, at 30 DAP stage contig7442 and contig7857 were found to be down-regulated (3-fold) compared to those in flower stage. Contig16888 was down-regulated whereas contig17138 was up-regulated at 10 DAH than those at 5 DAH. At 15 DAH contig9072 was down-regulated (1.5-fold) compared to that at 10 DAH. Among the galactosidase class, nine contigs encoding α-galactosidase and 21 contigs encoding β-galactosidase were detected. Two of these were down-regulated (contig1095 and contig1096) at 30 DAP compared to those in flower. Contig3844 was down regulated in 90 DAP pulp compared to 60 DAP.

Transition from 90 DAP to 5 DAH reflected in to up-regulation of contig3844 (2.72-fold) whereas, contig1525 coding for α-galactosidase was down-regulated in 10 DAH compared to that at 5 DAH. In the mannosidase class, five α-mannosidase and eight endo-β-mannosidase coding contigs were identified out of which two encoding for endo-beta-mannosidase (contig15554 and contig15558) were down-regulated at 10 DAH compared to 5 DAH, rest didn’t show differential regulation. Among two α-xylosidase and five β-xylosidase from Alphonso mango only contig1633 showed differential regulation which was up-regulated at 60 DAP compared to 30 DAP and was further down-regulated in both the 90 DAP tissues. Two transcripts coding for the acid beta-fructofuranosidase (contig12501 and contig12502) were detected but did not show any differential regulation in various tissues analyzed. Also, seven contigs coding for the alpha-l-fucosidase were identified from Alphonso mango but did not show differential regulation. Degradation of plant cell wall components, namely cellulose, hemi-cellulose and pectin by cellulases and glucanases in ripening fruits is responsible for the fruit softness.

Cellulases encoding five contigs were identified among which one coded for acidic cellulase and found to be abundant during flower and early fruit developing stages. Total 18 contigs encoding glucanase were detected of which contig4148 was up-regulated in flower compared to 30 DAP, whereas, four transcripts (contig9268, contig19283, contig9267 and contig17145) were down-regulated (2-fold) at 10 DAH compared to 5 DAH stage. Pectin is another component of fruit cell wall and is degraded by a set of enzymes viz. Pectate lyase (PL), pectin esterase (PE), polygalacturonase (PG) and rhamnogalacturonate lyase. In the present study 18 PL, 24 PG and 10 PE coding transcripts were detected. However, only few were differentially expressed, for example, only PL contig7696 was down-regulated in 30 DAP fruit compared to flower and PL contig9578 was up-regulated in 10 DAH fruit than 5 DAH fruit.

PG non-catalytic subunit jp650 coding contig9895 was down-regulated in 5 DAH pulp than in 90 DAP pulp, whereas PG contig3471 was down regulated and PG contig1614 was up-regulated in 10 DAH compared to 5 DAH stage. Likewise, PE contig7997 was down-regulated and contig9162 was up-regulated at 30 DAP stage compared to flower.

Similarly, eight transcripts encoding rhamnogalacturonate lyase were detected of which two (contig7745 and contig7746) were distinct to the ripening stages (Supplementary Table ). Discussion Recently transcriptome studies on mango cultivars namely, Zill, Langra, Kent and Dashehari have put forth important information regarding the fruit and leaf physiology of mango. These studies have identified genes encoding multiple enzymes involved in various pathways of primary and secondary metabolism such as, citrate cycle, glycolysis and gluconeogenesis from carbohydrate metabolism; fatty acid biosynthesis, beta oxidation and salicylic acid biosynthesis from fatty acid metabolism; biosynthesis and degradation of various amino acids as well as ethylene biosynthesis from methionine.

Genes involved in the flavonoid biosynthesis, vitamin biosynthesis (β-carotene and α-tocopherols) as well as terpenoid backbone synthesis (mevalonate pathway) have also been well explored. In the present study genes involved in all these pathways were identified and their differential expression was also evident through various stages of the Alphonso mango fruit development and ripening (Supplementary Table ). In addition the present study revealed some novel findings highlighting better understanding of various processes involved in mango fruit development and ripening and some were unique to the most favoured Alphonso mango fruit, which are discussed below. Quantitative abundance of terpenes in mangos is well known, Transcriptome and gene expression studies in mango have explored the terpene biosynthesis pathway till GPP, FPP and GGPP synthesis, –. In this study, we identified six, five and three genes encoding mono-terpene synthases (MTPS), sesqui- terpene synthases (STPS) and di-terpene synthases (DTPS), respectively that are involved in biosynthesis of specific terpene molecules (Supplementary Figure ). These genes were abundant in the flower tissue followed by 30 DAP stage.

Further, the transcript abundance of many of these terpene synthase genes in the present study has been depicted to be idiosyncratic to the developing stages, leading to their least expression in the ripening stages of Alphonso fruit (Figs and, Supplementary Figure and Supplementary Table ). Previous aroma volatile analysis from Alphonso mango supports this observation wherein flower had the highest concentration of mono-terpenes, oxygenated mono-terpenes and sesqui-terpenes which decreased through the fruit development. Quantitative reverse-transcriptase PCR validation of various transcripts obtained through RNAseq from mono-terpene metabolism through various tissues (flower and fruit developing and ripening stages). Vertical bars at each data point represent. Another flavor related pathway is lipoxygenase (LOX) followed by HPL pathway –, which produces C6 GLVs and lactones through peroxygenase pathway. The transcriptomes analyzed from Kent and Dashehari mangos reported the presence of six and five genes coding for the LOX family, respectively,.

Here, we report detailed annotation of Alphonso mango LOX genes wherein, two genes encode the 9-LOX and six genes encode the 13-LOX. Likewise, Peroxygenase and epoxide hydrolase ( EH) genes have been well-studied for biosynthesis of cutin biopolymer and defense related compounds,; while our recent study has shown involvement of these genes in the production of lactones. In spite of their biological significance, none of the previous transcriptome studies on mango could determine the presence of peroxygenase and epoxide hydrolase genes. On the other hand in the current study various contigs encoding novel EH (three contigs), EH2 (four contigs) and peroxygenase (three contigs) were detected. Similarly, multiple transcripts encoding enone oxidoreductase and O-methyltransferase having role in furanone biosynthesis, were also identified and expression profiles of many of them were shown to be ripening specific (Supplementary Figure ).

Abundance and ripening related expression of large number of these unique flavor related genes in Alphonso mango (Supplementary Figure ) signifies their role in synthesis of diverse aroma volatiles and their unique blend giving sweet and fruity flavor in Alphonso as shown by our previous studies,. Fruit ripening is a complex physiological process and can be characterized by means of fruit softening due to changes in the cell wall structure –, increased sugar content by polysaccharide hydrolysis and changes in the aroma volatiles. Starch and pectin are the major storage and structural polysaccharides in the mango fruit, respectively. Various hydrolases and lyases are known to carry out polysaccharide and cell wall hydrolysis, –, while amylases degrade the starch in to soluble sugars during ripening.

We identified four α-amylase, three isoamylase and 13 β-amylase coding transcripts. Only one transcript coding for α-amylase (contig1439) was down-regulated through flower to 30 DAP fruit transition while others were expressed throughout the fruit developing and ripening stages. On the other hand, transcriptome analysis of Kent mango revealed identification of four β-amylase and three α-amylase transcripts out of which only two β-amylase coding transcripts were found to be up-regulated in ripe tissue, none of the other showed differential expression. In addition, we also detected large number of PL, PG and PE encoding transcripts known to be responsible for degradation of complex hetero-polysaccharide, pectin were detected (Supplementary Table ). Most of these had stable expression and only few were differentially expressed in Alphonso mango. These results are in contrast to the observations from Kent and Dashehari mango transcriptomic data in terms of the number of unigenes detected and their differential regulation.

In case of Dashehari mango four PL and none of the PE or PG encoding transcripts were up-regulated; whereas in Kent, four PL, six PE and six PG encoding unigenes were up-regulated. These results signify controlled and steady activity of pectin degradation leading to slow and balanced transitions during Alphonso fruit ripening and this might be one of the reasons for its longer shelf life. Another interesting observation in our study was, the identification of 79 contigs encoding 20 different inhibitor classes from Alphonso mango (Supplementary Figure ). Previously, only three unigenes coding for the cysteine proteinase inhibitors were reported from Langra leaf transcriptome, while no such inhibitors were reported from Zill, Kent and Dashehari mango transcriptomes.

The present study suggests Alphonso mango transcriptome to be rich in these inhibitors throughout the fruit development and ripening stages. Thus, the presence of large numbers of amylase, PL, PE and PG transcripts with very few of them differentially regulated, perpetual expression of most of the starch and cell wall hydrolyzing enzymes along with the persistent presence of inhibitors for amylase, pectinesterase, polygalacturonase and ethylene sensitivity can be cumulatively suggested to play a crucial role in controlled and slow ripening and longer fruit shelf life of Alphonso mango.

As a representative example the balance between pectinesterases and their inhibitors during Alphonso fruit development and ripening has been demonstrated in Fig. As heatmap and area chart of expression of PE and PEI. We further observed that phosphate metabolism related gene ontologies were enriched at complete ripe stage (15 DAH) of Alphonso mango. Among these some transcripts were uniquely found at 15 DAH stage (Supplementary Table ) suggesting accelerated primary metabolism. Moreover, contig 6946 encoding cytokinin riboside 5 -monophosphate phosphoribohydrolase was uniquely observed at 15 DAH stage. This enzyme plays important role in converting endogenous inactive cytokinin nucleotides to the biologically active free cytokinin responsible for delayed ripening,.

On the other hand contig6770 encoding phospholipase-D (PLD) involved in degradation of important cell membrane component i.e. Phospholipids was observed to be upregulated (3.4fold). A key role of PLD in softening was well explained in tomato plants transformed with antisense PLD cDNA which showed delayed ripening and increased firmness in tomato fruits.

This counter play at molecular level might be responsible for fine tuning of ripening process in Alphonso mango providing longer shelf life to the fruits. ( a) Heatmap representing expression profiles of contigs encoding pectinesterase (PE) and pectinesterase inhibitor (PEI) from RNAseq data. ( b) Representation of percent co-expression (Y-axis) of 10 contigs encoding pectinesterase and 12 contigs encoding. Additionally oxidative burst, oxidoreductase activities and oxidative stress related gene ontologies were observed during ripening. These factors can generate reactive oxygen species and lead to cell death and fruit damage. Such reactive oxygen species induced cell death is suppressed by Bax inhibitor, which was well studied in Arabidopsis thaliana. Interestingly, four contigs coding for Bax inhibitor in Alphonso mango showed ripening specific expression probably responsible for preventing cell death in ripening fruits and thus longer shelf life.

Further detailed study on these inhibitors might help to understand jelly formation in Dashehari mangos due to excessive ripening and spongy tissue formation in Alphonso mangos due to uneven ripening,. It is known that ripened fruits are prone to attack by various pathogens. A well distinct defense mechanism was observed in Alphonso mango wherein various defense related proteins (227 contigs) and chitinases (19 contigs) acting on fungal cell wall were differentially regulated (Supplementary Table ). Chitinases were found to be accumulated in the flower and in the early fruit developing stages. Insect driven pollination has the risk of fungal infection to flower and further spore accumulation around ovary causing internal infection to the fruit. This might be restricted by the presence of various chitinases in Alphonso mango. Similarly, various disease resistance proteins might play role during fruit development and ripening process to defend infections.

Conclusions The transcriptome of Alphonso mango analyzed through eight stages of flower to fruit development and ripening transitions revealed various differentially regulated and stage specific genes. Unique transcript profiles probably responsible for distinct and favorable characteristics of Alphonso mango fruit such as flavor, color, ripening duration, skin to stone ripening pattern and longer shelf life were identified and analyzed. This study provides large data sets for further functional validation of fruit ripening process. Plant material Flower and fruits of cv. Alphonso were collected from three independent trees as biological replicates from the Mango Research Sub Centre, Deogad (16.528336 N, 73.344790 E) affiliated to Dr. Balasaheb Sawant Konkan Agricultural University, Dapoli, Maharashtra, India.

Flowers from inflorescence were collected and snap frozen. Fruits from developing stages were collected at 30, 60 and 90 days after pollination (DAP). Fruits from 30 and 60 DAP were analyzed as whole fruit in the present study after removing fruit stone; whereas at mature raw stage (90 DAP) pulp (mesocarp) and skin (exocarp) were separated, snap frozen in liquid nitrogen and stored at −80 °C until further analysis. A set of fruits were additionally harvested at their mature raw stage and kept in the hay containing boxes at ambient temperature for ripening and only pulp tissue for ripening stages as table green, mid ripe and ripe were collected at 5, 10 and 15 days after harvest (DAH), respectively. At each ripening stage fruits were removed from the box, pulp and skin were separated and pulp was frozen in liquid nitrogen and stored at −80 °C till further use. RNA isolation and cDNA synthesis Total RNA isolation was carried out for all the tissues sampled for current study using RNeasy Plus mini kit (Quiagen, Hilden, Germany). RNA quality as 260 nm/ 280 nm ratio was evaluated using Nanodrop 1000 (Thermo Fisher Scientific, Waltham, Massachusetts, USA) and RNA integrity was checked using Bioanalyzer 2100 (Agilent Technologies, Santa Clara, USA).

Two microgram of total RNA was used to carry out reverse transcription for synthesis of cDNA using High Capacity cDNA reverse transcription kit (Applied Biosystem, Carlsbad, CA, USA). Library preparation and sequencing High quality RNAs from seven progressive stages of fruit development and ripening as well as flower tissue from single representative biological replicate were sent to the Next Generation Genomics Facility (NGGF) at the Centre for Cellular and Molecular Platforms (C-CAMP), Bangalore for performing transcriptome sequencing. Briefly, 1 ug of total RNA from each sample was used to prepare eight individual libraries and mRNA was purified using polyT oligo beads. The purified mRNA was fragmented in the range of 100 to 140 bases with optimum at around 120 bases and the cDNA was synthesized. End repair, A-Tailing, Adapter ligation and the library preparation were performed using Tru Seq RNA sample preparation kit v2 (Illumina, San Diego, USA) as per manufacturer’s instructions.

PCR enrichment was performed for 15 cycles and the sample was validated on the Bioanalyzer 2100. Libraries were sequenced in a Paired End 100 base run, using TruSeq SBS Kit v3-HS (Catalog No.: FC-401-3001) for sequencing on the Illumina HiSeq. 1000 platform according to the manufacturer’s recommended protocols. Bioinformatics data analysis Paired end RNA sequencing was performed using Illumina Hiseq2000. Quality check on raw data was performed using FastQC.

Adapter free, good quality reads (Q = 30; min read length = 85) were obtained using Cutadapt. Alphonso transcriptome for each stage was assembled using Velvet-Oases with K-mers 67, 75, 83 and merging them at 27 k-mer. Additionally, a merged transcriptome was also generated using k-mer 55. Further, for all the merged assembled transcripts, we used Transdecoder to extract potential candidate coding regions within transcripts. Partial cds were discarded and only those transcripts with start and stop codon were made non redundant based upon sequence identity cut-off of 90% using CD-HIT-est and used for downstream analyses. Further, merged full length transcripts (from merged assembly) were used as reference to map back all the raw reads from each stage using default parameters of Bowtie. DESeq2 was used to identify differentially expressed transcripts and were filtered based upon p-value ≤ 0.05 and expression value 0.

All those transcripts having mapping count zero were excluded from further analysis such as in identification of uniquely expressed transcripts in a particular stage or a specific set of stages (e.g. Developing and ripening stages). Unique and common list of transcripts were represented using Venny and provided in Supplementary Table. Full length transcripts (from merged assembly) were used as reference for differential expression analysis using DESeq2. Annotation, enzyme code distribution and GO mapping and InterProScan analyses were carried using BLAST2GO 3.1.3 workbench (Biobam Bioinformatics S.L., Valencia, Spain) as described in the user manual. GO enrichment analysis was carried out in given test and reference sets by Fisher’s exact test in BLAST2GO with p-value (0.001) and FDR filters. Quantitative reverse-transcriptase PCR Quantitative reverse-transcriptase PCR was performed using the Fast Start Universal SYBR Green master mix (Roche Inc.

Stages Of Fruit Ripening

Indianapolis, Indiana, USA) and elongation factor 1α ( EF1α) as an endogenous reference gene employing the primers reported earlier. Various transcripts selected from transcriptome data were amplified using gene specific primers (Supplementary Table ) and quantification was performed using 7900HT Fast Real-Time PCR System (Applied Biosystems, California, USA) having thermal cycle program of initial denaturation at 95 °C for 10 min with subsequent 40 cycles of 95 °C for 3 sec and 60 °C for 30 sec followed by a dissociation curve analysis of transcripts.

Relative quantification (ΔΔCT method) and statistical analysis was carried out using DataAssist™ v3.01 software (Applied Biosystems, California, USA). Eight tissue samples as described above from three independent biological replicates were used for this analysis. For individual transcript the highest expression at a particular fruit development or ripening stage was considered as 1 and expression in other tissues including flower was normalized to that for better graphical presentation. The authors would like to thank Dr. Aarti Desai for helpful discussion and Dr.

Abhay Jere (Persistent Systems Ltd., Pune, India) for his critical suggestions, guidance and support. The authors acknowledge the Next Generation Genomics Facility (NGGF) at the Centre for Cellular and Molecular Platforms (C-CAMP), Bangalore (Project ID: BT/PR3481/INF/22/140/2011) for performing transcriptome sequencing of the eight mango tissues. Authors would like to thank Dr. Sanjiv Galande, Dr. Ram Kulkarni and Mr. Saurabh Pradhan (Indian Institute of Science Education and Research, Pune, India) for their help in evaluation of RNA quality by Bioanalyzer 2100.

This research was funded by the Council of Scientific and Industrial Research, New Delhi, India under the project CSC0133 (FUNHEALTH) to CSIR-NCL and 21(0997)/16/EMR-II grant as Emeritus Scientist scheme to VSG. Author Contributions A.B.D.

Designed, standardized and performed majority of the experiments, analyzed the data and wrote the manuscript. Carried out bioinformatics analysis.

Helped with high quality RNA isolation, cDNA synthesis and qRT PCR studies. Contributed by providing the experimental tissue. Contributed with suggestions for the work and for manuscript planning, while NYK contributed in bioinformatic analysis and manuscript writing. Designed the study, participated in all parts of the work and contributed to the preparation of the manuscript.

All authors read and approved the final manuscript.