De novo Transcriptome Assembly of Floral Buds of Pineapple and Identification of Differentially Expressed Genes in Response to Ethephon Induction

Comparative transcriptome analysis reveals candidate genes for cold stress response and early flowering in pineapple

Pineapple originates from tropical regions in South America and is therefore significantly impacted by cold stress. Periodic cold events in the equatorial regions where pineapple is grown may induce early flowering, also known as precocious flowering, resulting in monetary losses due to small fruit size and the need to make multiple passes for harvesting a single field. Currently, pineapple is one of the most important tropical fruits in the world in terms of consumption, and production losses caused by weather can have major impacts on worldwide exportation potential and economics. To further our understanding of and identify mechanisms for low-temperature tolerance in pineapple, and to identify the relationship between low-temperature stress and flowering time, we report here a transcriptomic analysis of two pineapple genotypes in response to low-temperature stress. Using meristem tissue collected from precocious flowering-susceptible MD2 and precocious flowering-tolerant Dole-17, we performed pairwise comparisons and weighted gene co-expression network analysis (WGCNA) to identify cold stress, genotype, and floral organ development-specific modules. Dole-17 had a greater increase in expression of genes that confer cold tolerance. The results suggested that low temperature stress in Dole-17 plants induces transcriptional changes to adapt and maintain homeostasis. Comparative transcriptomic analysis revealed differences in cuticular wax biosynthesis, carbohydrate accumulation, and vernalization-related gene expression between genotypes. Cold stress induced changes in ethylene and abscisic acid-mediated pathways differentially between genotypes, suggesting that MD2 may be more susceptible to hormone-mediated early flowering. The differentially expressed genes and module hub genes identified in this study are potential candidates for engineering cold tolerance in pineapple to develop new varieties capable of maintaining normal reproduction cycles under cold stress. In addition, a total of 461 core genes involved in the development of reproductive tissues in pineapple were also identified in this study. This research provides an important genomic resource for understanding molecular networks underlying cold stress response and how cold stress affects flowering time in pineapple.

Differential gene expression during floral transition in pineapple

Pineapple (Ananas comosus var. comosus) and ornamental bromeliads are commercially induced to flower by treatment with ethylene or its analogs. The apex is transformed from a vegetative to a floral meristem and shows morphological changes in 8 to 10 days, with flowers developing 8 to 10 weeks later. During eight sampling stages ranging from 6 h to 8 days after treatment, 7961 genes were found to exhibit differential expression (DE) after the application of ethylene. In the first 3 days after treatment, there was little change in ethylene synthesis or in the early stages of the ethylene response. Subsequently, three ethylene response transcription factors (ERTF) were up-regulated and the potential gene targets were predicted to be the positive flowering regulator CONSTANS-like 3 (CO), a WUSCHEL gene, two APETALA1/FRUITFULL (AP1/FUL) genes, an epidermal patterning gene, and a jasmonic acid synthesis gene. We confirm that pineapple has lost the flowering repressor FLOWERING LOCUS C. At the initial stages, the SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) was not significantly involved in this transition. Another WUSCHEL gene and a PHD homeobox transcription factor, though not apparent direct targets of ERTF, were up-regulated within a day of treatment, their predicted targets being the up-regulated CO, auxin response factors, SQUAMOSA, and histone H3 genes with suppression of abscisic acid response genes. The FLOWERING LOCUS T (FT), TERMINAL FLOWER (TFL), AGAMOUS-like APETELAR (AP2), and SEPETALA (SEP) increased rapidly within 2 to 3 days after ethylene treatment. Two FT genes were up-regulated at the apex and not at the leaf bases after treatment, suggesting that transport did not occur. These results indicated that the ethylene response in pineapple and possibly most bromeliads act directly to promote the vegetative to flower transition via APETALA1/FRUITFULL (AP1/FUL) and its interaction with SPL, FT, TFL, SEP, and AP2. A model based on AP2/ERTF DE and predicted DE target genes was developed to give focus to future research. The identified candidate genes are potential targets for genetic manipulation to determine their molecular role in flower transition.

Effect of plant growth regulators DA-6 and COS on drought tolerance of pineapple through bromelain and oxidative stress

BACKGROUND

Due to global warming, drought climates frequently occur on land, and despite being drought resistant, pineapples are still subjected to varying degrees of drought stress. Plant growth regulators can regulate the stress tolerance of plants through hormonal effects. This experiment aims to investigate the regulatory effects of different plant growth regulators on Tainong- 16 and MD-2 Pineapple when subjected to drought stress.

RESULTS

In this experiment, we examined the regulatory effects of two different plant growth regulators, sprayed on two pineapple varieties: MD-2 Pineapple and Tainong-16. The main component of T1 was diethyl aminoethyl hexanoate (DA-6) and that of T2 is chitosan oligosaccharide (COS). An environment similar to a natural drought was simulated in the drought stress treatments. Then, pineapples at different periods were sampled and a series of indicators were measured. The experimental results showed that the drought treatments treated with T1 and T2 plant growth regulators had a decrease in malondialdehyde, an increase in bromelain and antioxidant enzyme indicators, and an increase in phenotypic and yield indicators.

CONCLUSION

This experiment demonstrated that DA-6 and COS can enhance the drought resistance of pineapple plants to a certain extent through bromelain and oxidative stress. Therefore, DA-6 and COS have potential applications and this experiment lays the foundation for further research.

Functional diversification and molecular mechanisms of FLOWERING LOCUS T/TERMINAL FLOWER 1 family genes in horticultural plants

Flowering is an important process in higher plants and is regulated by a variety of factors, including light, temperature, and phytohormones. Flowering restriction has a considerable impact on the commodity value and production cost of many horticultural crops. In Arabidopsis, the FT/TFL1 gene family has been shown to integrate signals from various flowering pathways and to play a key role in the transition from flower production to seed development. Studies in several plant species of the FT/TFL1 gene family have revealed it harbors functional diversity in the regulation of flowering. Here, we review the functional evolution of the FT/TFL1 gene family in horticulture plants and its unique regulatory mechanisms; in addition, the FT/TFL1 family of genes as an important potential breeding target is explored.

A Complex Gene Network Mediated by Ethylene Signal Transduction TFs Defines the Flower Induction and Differentiation in Olea europaea L

The olive tree (Olea europaea L.) is a typical Mediterranean crop, important for olive and oil production. The high tendency to bear fruits in an uneven manner, defined as irregular or alternate bearing, results in a significant economic impact for the high losses in olives and oil production. Buds from heavy loaded (‘ON’) and unloaded (‘OFF’) branches of a unique olive tree were collected in July and the next March to compare the transcriptomic profiles and get deep insight into the molecular mechanisms regulating floral induction and differentiation. A wide set of DEGs related to ethylene TFs and to hormonal, sugar, and phenylpropanoid pathways was identified in buds collected from ‘OFF’ branches. These genes could directly and indirectly modulate different pathways, suggesting their key role during the lateral bud transition to flowering stage. Interestingly, several genes related to the flowering process appeared as over-expressed in buds from March ‘OFF’ branches and they could address the buds towards flower differentiation. By this approach, interesting candidate genes related to the switch from vegetative to reproductive stages were detected and analyzed. The functional analysis of these genes will provide tools for developing breeding programs to obtain olive trees characterized by more constant productivity over the years.

Floral transcriptomes reveal gene networks in pineapple floral growth and fruit development

Proper flower development is essential for sexual reproductive success and the setting of fruits and seeds. The availability of a high quality genome sequence for pineapple makes it an excellent model for studying fruit and floral organ development. In this study, we sequenced 27 different pineapple floral samples and integrated nine published RNA-seq datasets to generate tissue- and stage-specific transcriptomic profiles. Pairwise comparisons and weighted gene co-expression network analysis successfully identified ovule-, stamen-, petal- and fruit-specific modules as well as hub genes involved in ovule, fruit and petal development. In situ hybridization confirmed the enriched expression of six genes in developing ovules and stamens. Mutant characterization and complementation analysis revealed the important role of the subtilase gene AcSBT1.8 in petal development. This work provides an important genomic resource for functional analysis of pineapple floral organ growth and fruit development and sheds light on molecular networks underlying pineapple reproductive organ growth.

Small RNA sequencing revealed various microRNAs involved in ethylene-triggered flowering process in Aechmea fasciata

Ethylene-triggered flowering is a common phenomenon in plants of the family Bromeliaceae, but its molecular mechanism remains unclear. As a classical group of small RNAs, microRNAs play an essential role in the regulation of flowering. In this study, we found that various miRNAs participate in the ethylene-triggered flowering process in Aechmea fasciata via small RNA sequencing using juvenile and adult plants treated with ethylene for 24 hours. Finally, 63 known miRNAs, 52 novel miRNAs and 1721 target genes were identified or predicted. Expression changes of specific miRNAs were validated by qRT-PCR and northern blotting. Some predicted targets, including SPL, GAMYB and ARF, were verified in RLM-RACE experiments. Gene Ontology (GO) and KEGG analysis showed that numerous developmental and RNA-related processes were enriched. Integrated analysis of the transcriptomic data with small RNA sequencing revealed that numerous miRNAs and targets involved in ethylene-triggered flowering in A. fasciata. Our study is helpful for illuminating the molecular basis of the ethylene-triggered flowering phenomenon in Bromeliaceae.

Novel findings to the biosynthetic pathway of magnoflorine and taspine through transcriptomic and metabolomic analysis of Croton draco (Euphorbiaceae)

BACKGROUND

Croton draco is an arboreal species and its latex as well as some other parts of the plant, are traditionally used in the treatment of a wide range of ailments and diseases. Alkaloids, such as magnoflorine, prevent early atherosclerosis progression while taspine, an abundant constituent of latex, has been described as a wound-healer and antitumor-agent. Despite the great interest for these and other secondary metabolites, no omics resources existed for the species and the biosynthetic pathways of these alkaloids remain largely unknown.

RESULTS

To gain insights into the pathways involved in magnoflorine and taspine biosynthesis by C. draco and identify the key enzymes in these processes, we performed an integrated analysis of the transcriptome and metabolome in the major organs (roots, stem, leaves, inflorescences, and flowers) of this species. Transcript profiles were generated through high-throughput RNA-sequencing analysis while targeted and high resolution untargeted metabolomic profiling was also performed. The biosynthesis of these compounds appears to occur in the plant organs examined, but intermediaries may be translocated from the cells in which they are produced to other cells in which they accumulate.

CONCLUSIONS

Our results provide a framework to better understand magnoflorine and taspine biosynthesis in C. draco. In addition, we demonstrate the potential of multi-omics approaches to identify candidate genes involved in the biosynthetic pathways of interest.

Transcriptional foliar profile of the C3-CAM bromeliad Guzmania monostachia

Guzmania monostachia is an epiphytic tank bromeliad that displays the inducible CAM photosynthesis under stressful conditions and had the highest stomata density in the leaf apex, while the base portion has the highest density of trichomes, which are specialized structures used to acquire water and nutrients from the tank solution. In order to correlate the genetic factors behind these morpho-physiological characteristics along the leaf blade of G. monostachia, a comparative transcriptome analysis was performed to identify the functional enriched pathways and unigenes that could play a role in the apical, middle and basal leaf portions. A total of 653 million reads were used for de novo transcriptome assembly, resulting in 48,051 annotated unigenes. Analysis of differentially expressed genes (DEGs) among distinct leaf regions revealed that 806 DEGs were upregulated in the apex compared to the middle portion, while 9685 DEGs were upregulated in the apex and 9784 DEGs were upregulated in the middle portions compared to the base. Our outcomes correlated some DEGs and identified unigenes with their physiological functions, mainly suggesting that the leaf apex was related to the regulation of stomatal movement, production of chlorophyll, cellular response to stress, and H2O2 catabolic process. In contrast, the middle portion showed DEGs associated with the transport of amino acids. Furthermore, DEGs from the leaf base were mainly correlated with responses to nutrients and nitrogen compounds, regulation of potassium ion import, response to water deprivation, and trichome branching, indicating that, at least in part, this leaf portion can replace some of the root functions of terrestrial plants. Therefore, possibly candidate unigenes and enriched pathways presented here could be prospected in future experimental work, opening new possibilities to bioengineer non-inducible CAM plants and/or improve the fertilization use efficiency by increasing leaf nutrient acquisition of crop plants.

Comparative Transcriptome Analysis between Fertile and CMS Flower Buds in Wucai (Brassica campestris L.)

Background

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a variant of nonheading Chinese cabbage (Brassica campestris L.), which is one of the major vegetables in China. Cytoplasmic male sterility (CMS) has been used for Wucai breeding in recent years. However, the underlying molecular mechanism of Wucai CMS remains unclear. In this study, the phenotypic and cytological features of Wucai CMS were observed by anatomical analysis, and a comparative transcriptome analysis was carried out to identify genes related to male sterility using Illumina RNA sequencing technology (RNA-Seq).

Results

Microscopic observation demonstrated that tapetum development was abnormal in the CMS line, which failed to produce fertile pollen. Bioinformatics analysis detected 4430 differentially expressed genes (DEGs) between the fertile and sterile flower buds. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to better understand the functions of these DEGs. Among the DEGs, 35 genes (53 DEGS) were implicated in anther and pollen development, and 11 genes were involved in pollen cell wall formation and modification; most of these showed downregulated expression in sterile buds. In addition, several genes related to tapetum development (A6, AMS, MS1, MYB39, and TSM1) and a few genes annotated to flowering (CO, AP3, VIN3, FLC, FT, and AGL) were detected and confirmed by qRT-PCR as being expressed at the meiosis, tetrad, and uninucleate microspore stages, thus implying possible roles in specifying or determining the fate and development of the tapetum, male gametophyte and stamen. Moreover, the top four largest transcription factor families (MYB, bHLH, NAC and WRKY) were analyzed, and most showed reduced expression in sterile buds. These differentially expressed transcription factors might result in abortion of pollen development in Wucai.

Conclusion

The present comparative transcriptome analysis suggested that many key genes and transcription factors involved in anther development show reduced gene expression patterns in the CMS line, which might contribute to male sterility in Wucai. This study provides valuable information for a better understanding of CMS molecular mechanisms and functional genome studies in Wucai.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5331-4) contains supplementary material, which is available to authorized users.

Comprehensive tissue-specific transcriptome profiling of pineapple (Ananas comosus) and building an eFP-browser for further study

Pineapple is one of the most economically important tropical or subtropical fruit trees. However, few studies focus on the development of its unique collective fruit. In this study, we generated a genome-wide developmental transcriptomic profile of 14 different tissues of the collective fruit of the pineapple covering each of the three major fruit developmental stages. In total, 273 tissue-specific and 1,051 constitutively expressed genes were detected. We also performed gene co-expression analysis and 18 gene modules were classified. Among these, we found three interesting gene modules; one was preferentially expressed in bracts and sepals and was likely involved in plant defense; one was highly expressed at the beginning of fruit expansion and faded afterward and was probably involved in endocytosis; Another gene module increased expression level with pineapple fruit development and was involved in terpenoid and polyketide metabolism. In addition, we built a pineapple electronic fluorescent pictograph (eFP) browser to facilitate exploration of gene expression during pineapple fruit development. With this tool, users can visualize expression data in this study in an intuitive way. Together, the transcriptome profile generated in this work and the corresponding eFP browser will facilitate further study of fruit development in pineapple.

Integrated DNA methylome and transcriptome analysis reveals the ethylene-induced flowering pathway genes in pineapple

Ethylene has long been used to promote flowering in pineapple production. Ethylene-induced flowering is dose dependent, with a critical threshold level of ethylene response factors needed to trigger flowering. The mechanism of ethylene-induced flowering is still unclear. Here, we integrated isoform sequencing (iso-seq), Illumina short-reads sequencing and whole-genome bisulfite sequencing (WGBS) to explore the early changes of transcriptomic and DNA methylation in pineapple following high-concentration ethylene (HE) and low-concentration ethylene (LE) treatment. Iso-seq produced 122,338 transcripts, including 26,893 alternative splicing isoforms, 8,090 novel transcripts and 12,536 candidate long non-coding RNAs. The WGBS results suggested a decrease in CG methylation and increase in CHH methylation following HE treatment. The LE and HE treatments induced drastic changes in transcriptome and DNA methylome, with LE inducing the initial response to flower induction and HE inducing the subsequent response. The dose-dependent induction of FLOWERING LOCUS T-like genes (FTLs) may have contributed to dose-dependent flowering induction in pineapple by ethylene. Alterations in DNA methylation, lncRNAs and multiple genes may be involved in the regulation of FTLs. Our data provided a landscape of the transcriptome and DNA methylome and revealed a candidate network that regulates flowering time in pineapple, which may promote further studies.

Transcriptomic and Hormonal Analyses Reveal that YUC-Mediated Auxin Biogenesis Is Involved in Shoot Regeneration from Rhizome in Cymbidium

Cymbidium, one of the most important orchid genera in horticulture, can be classified into epiphytic and terrestrial species. Generally, epiphytic Cymbidium seedlings can be easily propagated by tissue culture, but terrestrial seedlings are difficult to propagate. To date, the molecular mechanisms underlying the differences in the ease with which terrestrial and epiphytic cymbidiums can be propagated are largely unknown. Using RNA-sequencing, quantitative reverse transcription PCR and enzyme-linked immunosorbent assay, Cymbidium 'Xiaofeng' (CXF), which can be efficiently micropropagated, and terrestrial Cymbidium sinense 'Qijianbaimo' (CSQ), which has a low regeneration ability, were used to explore the molecular mechanisms underlying the micropropagation ability of Cymbidium species. To this end, 447 million clean short reads were generated, and 31,264 annotated unigenes were obtained from 10 cDNA libraries. A total of 1,290 differentially expressed genes (DEGs) were identified between CXF and CSQ during shoot induction. Gene ontology (GO) enrichment analysis indicated that the DEGs were significantly enriched in auxin pathway-related GO terms. Further analysis demonstrated that YUC and GH3 family genes, which play crucial roles in the regulation of auxin/IAA (indole-3-acetic acid) metabolism, acted quickly in response to shoot induction culture in vitro and were closely correlated with variation in shoot regeneration between CXF and CSQ. In addition, the study showed that IAA accumulated rapidly and significantly during shoot induction in CXF compared to that in CSQ; in contrast, no significant changes in other hormones were observed between CXF and CSQ. Furthermore, shoot regeneration in CXF was inhibited by a yucasin-auxin biosynthesis inhibitor, indicating that increased IAA level is required for high-frequency shoot regeneration in CXF. In conclusion, our study revealed that YUC-mediated auxin biogenesis is involved in shoot regeneration from rhizome in Cymbidium.

Transcriptome Analysis and Identification of Genes Associated with Floral Transition and Flower Development in Sugar Apple (Annona squamosa L.)

Sugar apple (Annona squamosa L.) is a semi-deciduous subtropical tree that progressively sheds its leaves in the spring. However, little information is available on the mechanism involved in flower developmental pattern. To gain a global perspective on the floral transition and flower development of sugar apple, cDNA libraries were prepared independently from inflorescent meristem and three flowering stages. Illumina sequencing generated 107,197,488 high quality reads that were assembled into 71,948 unigenes, with an average sequence length of 825.40 bp. Among the unigenes, various transcription factor families involved in floral transition and flower development were elucidated. Furthermore, a Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that unigenes exhibiting differential expressions were involved in various phytohormone signal transduction events and circadian rhythms. In addition, 147 unigenes exhibiting sequence similarities to known flowering-related genes from other plants were differentially expressed during flower development. The expression patterns of 20 selected genes were validated using quantitative-PCR. The expression data presented in our study is the most comprehensive dataset available for sugar apple so far and will serve as a resource for investigating the genetics of the flowering process in sugar apple and other Annona species.