Transcriptome profiling of fruit development and maturation in Chinese white pear (Pyrus bretschneideri Rehd)

Transcriptome Analysis Reveals Fruit Quality Formation in Actinidia eriantha Benth

Actinidia chinensis Planch. is a fruit tree originating from China that is abundant in the wild. Actinidia eriantha Benth. is a type of A. chinensis that has emerged in recent years. The shape of A. eriantha is an elongated oval, and the skin is covered with dense, non-shedding milk-white hairs. The mature fruit has flesh that is bright green in colour, and the fruit has a strong flavour and a grass-like smell. It is appreciated for its rich nutrient content and unique flavour. Vitamin C, sugar, and organic acids are key factors in the quality and flavour composition of A. eriantha but have not yet been systematically analysed. Therefore, we sequenced the transcriptome of A. eriantha at three developmental stages and labelled them S1, S2, and S3, and comparisons of S1 vs. S2, S1 vs. S3, and S2 vs. S3 revealed 1218, 4019, and 3759 upregulated differentially expressed genes and 1823, 3415, and 2226 downregulated differentially expressed genes, respectively. Furthermore, the upregulated differentially expressed genes included 213 core genes, and Gene Ontology enrichment analysis showed that they were enriched in hormones, sugars, organic acids, and many organic metabolic pathways. The downregulated differentially expressed genes included 207 core genes, which were enriched in the light signalling pathway. We further constructed the metabolic pathways of sugars, organic acids, and vitamin C in A. eriantha and identified the genes involved in vitamin C, sugar, and organic acid synthesis in A. eriantha fruits at different stages. During fruit development, the vitamin C content decreased, the carbohydrate compound content increased, and the organic acid content decreased. The gene expression patterns were closely related to the accumulation patterns of vitamin C, sugars, and organic acids in A. eriantha. The above results lay the foundation for the accumulation of vitamin C, sugars, and organic acids in A. eriantha and for understanding flavour formation in A. eriantha.

In-silico mining and characterization of MYB family genes in wilt-resistant hybrid guava (Psidium guajava × Psidium molle)

BACKGROUND

The MYB family is one of the most significant groups of transcription factors in plants. However, several MYBs have been linked to secondary metabolism and are important for determining the color of fruit's peel and pulp. Despite being a substantial fruit crop in tropical and subtropical areas of the world, wilt-resistant hybrid guava (Psidium guajava × Psidium molle; PGPM) has not yet been the subject of a thorough examination. This study's goal was to assess the expression of MYB in guava fruit pulp, roots, and seeds to predict its function by in silico analysis of the guava root transcriptome data.

RESULTS

In the current study, we have mined the MYBs family of MYB genes from the transcriptome of the PGPM guava root. We have mined 15 distinct MYB transcription factor genes/transcripts viz MYB3, MYB4, MYB23, MYB86, MYB90, MYB308, MYB5, MYB82, MYB114, MYB6, MYB305, MYB44, MYB51, MYB46, and MYB330. From the analyses, it was found that R2-MYB and R3-MYB domains are conserved in all known guava MYB proteins. The expression of six different MYB TFs was examined using semi-quantitative RT-PCR in "Shweta" pulp (white colour pulp), "Lalit" pulp (red color pulp), "Lalit" root, and "Lalit" seed.

CONCLUSION

There were 15 MYB family members observed in guava. They were unequally distributed across the chromosomes, most likely as a result of gene duplication. Additionally, the expression patterns of the particular MYBs showed that MYB may be involved in the control of wilt, fruit ripening, seed development, and root development. Our results allow for a more thorough functional characterization of the guava MYB family genes and open the door to additional research into one essential MYB transcription factor family of genes and its involvement in the growth and ripening of guava fruit.

Metabolome and transcriptome profiling revealed the enhanced synthesis of volatile esters in Korla pear

BACKGROUND

Flavor contributes to the sensory quality of fruits, including taste and aroma aspects. The quality of foods is related to their flavor-associated compounds. Pear fruits have a fruity sense of smell, and esters are the main contributor of the aroma. Korla pear are well known due to its unique aroma, but the mechanism and genes related to volatile synthesis have not been fully investigated.

RESULTS

Flavor-associated compounds, including 18 primary metabolites and 144 volatiles, were characterized in maturity fruits of ten pear cultivars from five species, respectively. Based on the varied metabolites profiles, the cultivars could be grouped into species, respectively, by using orthogonal partial least squares discrimination analysis (OPLS-DA). Simultaneously, 14 volatiles were selected as biomarkers to discriminate Korla pear (Pyrus sinkiangensis) from others. Correlation network analysis further revealed the biosynthetic pathways of the compounds in pear cultivars. Furthermore, the volatile profile in Korla pear throughout fruit development was investigated. Aldehydes were the most abundant volatiles, while numerous esters consistently accumulated especially at the maturity stages. Combined with transcriptomic and metabolic analysis, Ps5LOXL, PsADHL, and PsAATL were screened out as the key genes in ester synthesis.

CONCLUSION

Pear species can be distinguished by their metabolic profiles. The most diversified volatiles as well as esters was found in Korla pear, in which the enhancement of lipoxygenase pathway may lead to the high level of volatile esters at maturity stages. The study will benefit the fully usage of pear germplasm resources to serve fruit flavor breeding goals.

Advances in genomics for diversity studies and trait improvement in temperate fruit and nut crops under changing climatic scenarios

Genetic improvement of temperate fruit and nut crops through conventional breeding methods is not sufficient alone due to its extreme time-consuming, cost-intensive, and hard-to-handle approach. Again, few other constraints that are associated with these species, viz., their long juvenile period, high heterozygosity, sterility, presence of sexual incompatibility, polyploidy, etc., make their selection and improvement process more complicated. Therefore, to promote precise and accurate selection of plants based on their genotypes, supplement of advanced biotechnological tools, viz., molecular marker approaches along with traditional breeding methods, is highly required in these species. Different markers, especially the molecular ones, enable direct selection of genomic regions governing the trait of interest such as high quality, yield, and resistance to abiotic and biotic stresses instead of the trait itself, thus saving the overall time and space and helping screen fruit quality and other related desired traits at early stages. The availability of molecular markers like SNP (single-nucleotide polymorphism), DArT (Diversity Arrays Technology) markers, and dense molecular genetic maps in crop plants, including fruit and nut crops, led to a revelation of facts from genetic markers, thus assisting in precise line selection. This review highlighted several aspects of the molecular marker approach that opens up tremendous possibilities to reveal valuable information about genetic diversity and phylogeny to boost the efficacy of selection in temperate fruit crops through genome sequencing and thus cultivar improvement with respect to adaptability and biotic and abiotic stress resistance in temperate fruit and nut species.

Analysis of the β-Glucosidase Family Reveals Genes Involved in the Lignification of Stone Cells in Chinese White Pear (Pyrus bretschneideri Rehd.)

BGLU β-glucosidases in glycoside hydrolase family 1 (GH1) are involved in many processes of plant secondary metabolism. In particular, its de-glycosylation function plays an important role in the transport of lignin monolignols. No comprehensive study of the BGLU family in Chinese pear (Pyrus bretschneideri Rehd.) has been reported yet. In this study, the 50 BGLU family members from Chinese white pear were identified. Three candidate genes, PbBGLU1, PbBGLU15, and PbBGLU16, that may be involved in lignin synthesis were screened by bioinformatics analysis and qRT-PCR. Subcellular localization showed that all three of these candidate genes were expressed in the extracellular region. Then, we analyzed the functions of PbBGLU1 and PbBGLU16. In situ hybridization analysis showed that PbBGLU1 transcripts were not only localized to some pulp cell walls, lignin deposition, and stone cell areas of a pear fruit, but that was also a small amount of enrichment in normal pear flesh cells. PbBGLU16 transcripts were only enriched in lignin deposition and stone cell areas of pear fruit. Enzyme activity analysis revealed that GST-PbBGLU1 and GST-PbBGLU16 had a stronger activity and higher catalytic efficiency for coniferin than syringin. In addition, GST-PbBGLU16 exhibited the higher activity and catalytic efficiency for the two substrates compared with GST-PbBGLU1. The transformation of PbBGLU1 and PbBGLU16 into Arabidopsis identified that the lignin contents of Arabidopsis BGLU-45 mutant, PbBGLU1-RE, and PbBGLU16-RE were not changed than that of wild-type. However, compared with wild-type Arabidopsis, the overexpression of the plant's lignin increased in varying degrees. The effect of PbBGLU16 on the lignin increment was greater than that of PbBGLU1 in Arabidopsis. In pear fruits, with transient overexpression of PbBGLU1, the contents of lignin and stone cells were significantly higher (0.01 < P < 0.05) than those with empty vector injection pear fruits. After transient expression of PbBGLU16, lignin in pear fruit increased significantly (0.01 < P < 0.05) and stone cells showed a very significant difference (P < 0.01) compared with the control group. However, RNA interference silenced these two genes in pear fruit, which seemed to have no impression on lignin and stone cells. This study provides a molecular biological basis for improving pear fruit quality at the molecular level.

MYB transcription factor family in sweet cherry (Prunus avium L.): genome-wide investigation, evolution, structure, characterization and expression patterns

BACK GROUND

MYB Transcription factors (TFs) are most imperative and largest gene family in plants, which participate in development, metabolism, defense, differentiation and stress response. The MYB TFs has been studied in various plant species. However, comprehensive studies of MYB gene family in the sweet cherry (Prunus avium L.) are still unknown.

RESULTS

In the current study, a total of 69 MYB genes were investigated from sweet cherry genome and classified into 28 subfamilies (C1-C28 based on phylogenetic and structural analysis). Microcollinearity analysis revealed that dispersed duplication (DSD) events might play an important role in the MYB genes family expansion. Chromosomal localization, the synonymous (Ks) and nonsynonymous (Ka) analysis, molecular characteristics (pI, weight and length of amino acids) and subcellular localization were accomplished using several bioinformatics tools. Furthermore, the members of distinct subfamilies have diverse cis-acting regions, conserved motifs, and intron-exon architectures, indicating functional heterogeneity in the MYB family. Moreover, the transcriptomic data exposed that MYB genes might play vital role in bud dormancy. The quantitative real-time qRT-PCR was carried out and the expression pattern indicated that MYB genes significantly expressed in floral bud as compared to flower and fruit.

CONCLUSION

Our comprehensive findings provide supportive insights into the evolutions, expansion complexity and functionality of PavMYB genes. These PavMYB genes should be further investigated as they seem to be brilliant candidates for dormancy manipulation in sweet cherry.

Quantitative proteomic analysis of pear (Pyrus pyrifolia cv. "Hosui") flesh provides novel insights about development and quality characteristics of fruit

A total of 6763 proteins were identified in the developing pear flesh, which were further screened for differentially expressed proteins related to fruit quality and ATP-binding cassette transporters. To obtain further details on changes in protein levels during fruit ripening and to identify and evaluate changes in various metabolic pathways that affect fruit quality, a proteomic method using tandem mass tags was implemented at three developmental stages in Pyrus pyrifolia cv. "Hosui" that identified 6763 proteins. Subcellular localization and Gene Ontology enrichment analysis revealed major functions of all identified proteins. Kyoto Encyclopedia of Genes and Genomes pathway analysis suggested that all metabolic processes are reflected in the up- and downregulation of differentially expressed proteins during fruit development, which play predominant roles in cell division, cell expansion, and fruit ripening. Among the examined differentially expressed proteins, 160 related to fruit quality, and 14 ATP-binding cassette transporters related to fruit development were identified and analyzed. The quantitative data were validated by parallel reaction monitoring, which confirmed the reliability of the experimental results.

Evidence for the Involvement of Vernalization-related Genes in the Regulation of Cold-induced Ripening in 'D'Anjou' and 'Bartlett' Pear Fruit

European pear (Pyrus communis L.) cultivars require a genetically pre-determined duration of cold-temperature exposure to induce autocatalytic system 2 ethylene biosynthesis and subsequent fruit ripening. The physiological responses of pear to cold-temperature-induced ripening have been well characterized, but the molecular mechanisms underlying this phenomenon continue to be elucidated. This study employed previously established cold temperature conditioning treatments for ripening of two pear cultivars, 'D'Anjou' and 'Bartlett'. Using a time-course transcriptomics approach, global gene expression responses of each cultivar were assessed at four stages of developmental during the cold conditioning process. Differential expression, functional annotation, and gene ontology enrichment analyses were performed. Interestingly, evidence for the involvement of cold-induced, vernalization-related genes and repressors of endodormancy release was found. These genes have not previously been described to play a role in fruit during the ripening transition. The resulting data provide insight into cultivar-specific mechanisms of cold-induced transcriptional regulation of ripening in European pear, as well as a unique comparative analysis of the two cultivars with very different cold conditioning requirements.

Transcriptomic and metabolomic analysis provides insights into anthocyanin and procyanidin accumulation in pear

BACKGROUND

Pear is one of the most important fruit crops worldwide. Anthocyanins and procyanidins (PAs) are important secondary metabolites that affect the appearance and nutritive quality of pear. However, few studies have focused on the molecular mechanism underlying anthocyanin and PA accumulation in pear.

RESULTS

We conducted metabolome and transcriptome analyses to identify candidate genes involved in anthocyanin and PA accumulation in young fruits of the pear cultivar 'Clapp Favorite' (CF) and its red mutation cultivar 'Red Clapp Favorite' (RCF). Gene-metabolite correlation analyses revealed a 'core set' of 20 genes that were strongly correlated with 10 anthocyanin and seven PA metabolites. Of these, PcGSTF12 was confirmed to be involved in anthocyanin and PA accumulation by complementation of the tt19-7 Arabidopsis mutant. Interestingly, PcGSTF12 was found to be responsible for the accumulation of procyanidin A3, but not petunidin 3, 5-diglucoside, opposite to the function of AtGSTs in Arabidopsis. Transformation with PcGSTF12 greatly promoted or repressed genes involved in anthocyanin and PA biosynthesis, regulation, and transport. Electrophoretic mobility shift and luciferase reporter assays confirmed positive regulation of PcGSTF12 by PcMYB114.

CONCLUSION

These findings identify a core set of genes for anthocyanin and PA accumulation in pear. Of these, PcGSTF12, was confirmed to be involved in anthocyanin and PA accumulation. Our results also identified an important anthocyanin and PA regulation node comprising two core genes, PcGSTF12 and PcMYB114. These results provide novel insights into anthocyanin and PA accumulation in pear and represent a valuable data set to guide future functional studies and pear breeding.

Comparative genomic analysis of the PAL genes in five Rosaceae species and functional identification of Chinese white pear

Phenylalanine ammonia lyase (PAL) plays an important role in the biosynthesis of secondary metabolites regulating plant growth response. To date, the evolutionary history of the PAL family in Rosaceae plants remains unclear. In this study, we identified 16 PAL homologous genes in five Rosaceae plants (Pyrus bretschneideri, Fragaria vesca, Prunus mume, Prunus persica, and Malus × domestica). We classified these PALs into three categories based on phylogenetic analysis, and all PALs were distributed on 13 chromosomes. We tracked gene duplication events and performed sliding window analysis. These results revealed the evolution of PALs in five Rosaceae plants. We predicted the promoter of the PbPALs by PLANT CARE online software, and found that the promoter region of both PbPAL1 and PbPAL3 have at least one AC element. The results of qRT-PCR analysis found that PbPAL1 and PbPAL2 were highly expressed in the stems and roots, while expression level of PbPAL3 was relatively low in different tissues. The expression of PbPAL1 and PbPAL2 increased firstly and then decreased at different developmental periods of pear fruit. Among them, the expression of PbPAL1 reached the highest level 55 days after flowering. Three PbPALs were induced by abiotic stress to varying degrees. We transfected PbPAL1 and PbPAL2 into Arabidopsis thaliana, which resulted in an increase in lignin content and thickening of the cell walls of intervascular fibres and xylem cells. In summary, this research laid a foundation for better understanding the molecular evolution of PALs in five Rosaceae plants. Furthermore, the present study revealed the role of PbPALs in lignin synthesis, and provided basic data for regulating lignin synthesis and stone cells development in pear plants.

Transcriptome profiling of 'Kyoho' grape at different stages of berry development following 5-azaC treatment

Background

5-Azacytidine (5-azaC) promotes the development of ‘Kyoho’ grape berry but the associated changes in gene expression have not been reported. In this study, we performed transcriptome analysis of grape berry at five developmental stages after 5-azaC treatment to elucidate the gene expression networks controlling berry ripening.

Results

The expression patterns of most genes across the time series were similar between the 5-azaC treatment and control groups. The number of differentially expressed genes (DEGs) at a given developmental stage ranged from 9 (A3_C3) to 690 (A5_C5). The results indicated that 5-azaC treatment had not very great influences on the expressions of most genes. Functional annotation of the DEGs revealed that they were mainly related to fruit softening, photosynthesis, protein phosphorylation, and heat stress. Eight modules showed high correlation with specific developmental stages and hub genes such as PEROXIDASE 4, CAFFEIC ACID 3-O-METHYLTRANSFERASE 1, and HISTONE-LYSINE N-METHYLTRANSFERASE EZA1 were identified by weighted gene correlation network analysis.

Conclusions

5-AzaC treatment alters the transcriptional profile of grape berry at different stages of development, which may involve changes in DNA methylation.

Isolation and Detection Methods of Plant miRNAs

Small RNAs (sRNAs) are RNAs of low abundance in organisms. Among sRNAs, miRNAs are included and represent approximately 10% of the total number of sRNAs. The isolation of sRNAs is critical for miRNA detection and analysis. The precipitation of low-molecular-weight (LMW) RNAs from total RNA extracts has allowed enrichment of sRNAs. Here, we describe a simple method to isolate sRNAs from different plant species. The main advantage of this method is that it does not need first an extraction of total RNA and it is not based on TRIzol^® reagent. This method has been successfully used for miRNA analyses by Northern blot assay and RT-qPCR (these techniques are as well described in this chapter), as well as sRNA library preparation.

Sequencing of Euscaphis konishii Endocarp Transcriptome Points to Molecular Mechanisms of Endocarp Coloration

Flower and fruit colors are of vital importance to the ecology and economic market value of plants. The mechanisms of flower and fruit coloration have been well studied, especially among ornamental flower plants and cultivated fruits. As people pay more attention to exocarp coloration, the endocarp coloration in some species has often been ignored. Here, we report on the molecular mechanism of endocarp coloration in three development stages of Euscaphis konishii. The results show that endocarp reddening is closely related to anthocyanin accumulation, and a total of 86,120 unigenes were assembled, with a mean length of 893 bp (N50 length of 1642 bp). We identified a large number of differentially expressed genes associated with endocarp coloration, including anthocyanin biosynthesis, carotenoid biosynthesis, and chlorophyll breakdown. The genes participating in each step of the anthocyanin biosynthesis were found in the transcriptome dataset, but a few genes were found in the carotenoid biosynthesis and chlorophyll breakdown. In addition, the candidate R2R3-MYB transcription factors and candidate glutathione S-transferase transport genes, which likely regulate the anthocyanin biosynthesis, were identified. This study offers a platform for E. konishii functional genomic research and provides a reference for revealing the regulatory mechanisms of endocarp reddening.

Expression of circular RNAs during C2C12 myoblast differentiation and prediction of coding potential based on the number of open reading frames and N6-methyladenosine motifs

The importance of circular RNAs (circRNAs) as regulators of muscle development and muscle-associated disorders is becoming increasingly apparent. To explore potential regulators of muscle differentiation, we determined the expression profiles of circRNAs of skeletal muscle C2C12 myoblasts and myotubes using microarray analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to explore circRNA functions. We also established competing endogenous RNA (ceRNA) networks using bioinformatics methods and predicted the coding potential of differentially expressed circRNAs. We found that 581 circRNAs were differentially regulated between C2C12 myoblasts and myotubes. Bioinformatics analysis suggested that the primary functions of the linear transcripts of the circRNAs were linked with organization of the cytoskeleton, calcium signaling, cell cycle, and metabolic pathways. ceRNA networks showed that the myogenic-specific genes myogenin, myocyte enhancer factor 2a, myosin heavy chain (Myh)-1, Myh7, and Myh7b could combine with 91 miRNAs and the top 30 upregulated circRNAs, forming 239 edges. According to the number of open reading frames and N6-methyladenosine motifs, we identified 224 circRNAs with coding potential, and performed GO and KEGG analyses based on the linear counterparts of 75 circRNAs. We determined that the 75 circRNAs were related to regulation of the actin cytoskeleton and metabolic pathways. We established expression profiles of circRNAs during C2C12 myoblast differentiation and predicted the function of differentially expressed circRNAs, which might be involved in skeletal muscle development. Our study offers new insight into the functions of circRNAs in skeletal muscle growth and development.

Systematic analysis and comparison of the PHD-Finger gene family in Chinese pear (Pyrus bretschneideri) and its role in fruit development

PHD-finger proteins, which belongs to the type of zinc finger family, and that play an important role in the regulation of both transcription and the chromatin state in eukaryotes. Currently, PHD-finger proteins have been well studied in animals, while few studies have been carried out on their function in plants. In the present study, 129 non-redundant PHD-finger genes were identified from 5 Rosaceae species (pear, apple, strawberry, mei, and peach); among them, 31 genes were identified in pear. Subsequently, we carried out a bioinformatics analysis of the PHD-finger genes. Thirty-one PbPHD genes were divided into 7 subfamilies based on the phylogenetic analysis, which are consistent with the intron-exon and conserved motif analyses. In addition, we identified five segmental duplication events, implying that the segmental duplications might be a crucial role in the expansion of the PHD-finger gene family in pear. The microsynteny analysis of five Rosaceae species showed that there were independent duplication events in addition to the genome-wide duplication of the pear genome. Subsequently, ten expressed PHD-finger genes of pear fruit were identified using qRT-PCR, and one of these genes, PbPHD10, was identified as an important candidate gene for the regulation of lignin synthesis. Our research provides useful information for the further analysis of the function of PHD-finger gene family in pear.

A Novel LncRNA, MuLnc1, Associated With Environmental Stress in Mulberry (Morus multicaulis)

Environmental stresses are major constraints that limit the leaf productivity and quality of mulberry. LncRNAs have emerged as important regulators in response to biotic and abiotic stresses in plants. However, the functions and mechanisms of most lncRNAs remain largely unknown. A novel lncRNA designated as MuLnc1 was found to be cleaved by mul-miR3954 and produce secondary siRNAs in a 21 nt phase in mulberry. It was demonstrated that one of the siRNAs produced, si161579, can silence the expression of the calmodulin-like protein gene CML27 of mulberry (MuCML27). When MuCML27 was heterologously expressed in Arabidopsis, the transgenic plants exhibited enhanced resistance to Botrytis cinerea and Pseudomonas syringae pv tomato DC3000. In addition, the transgenic MuCML27-overexpressing Arabidopsis plants are more tolerant to salt and drought stresses. Furthermore, the network of mul-miR3954-MuLnc1-siRNAs-mRNAs was modeled to elucidate the interaction between lncRNAs and sRNAs with mRNAs. All of these, taken together, suggest that MuLnc1 was associated with environmental stress in mulberry and may be considered as a potential genetic improvement target gene of mulberry. The information provided may shed light on the complicated gene expression regulatory mechanisms in mulberry stress responses.

Using Next-Generation Sequencing to Detect Differential Expression Genes in Bradysia odoriphaga after Exposure to Insecticides

Bradysia odoriphaga (Diptera: Sciaridae) is the most important pest of Chinese chive. Insecticides are used widely and frequently to control B. odoriphaga in China. However, the performance of the insecticides chlorpyrifos and clothianidin in controlling the Chinese chive maggot is quite different. Using next generation sequencing technology, different expression unigenes (DEUs) in B. odoriphaga were detected after treatment with chlorpyrifos and clothianidin for 6 and 48 h in comparison with control. The number of DEUs ranged between 703 and 1161 after insecticide treatment. In these DEUs, 370–863 unigenes can be classified into 41–46 categories of gene ontology (GO), and 354–658 DEUs can be mapped into 987–1623 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The expressions of DEUs related to insecticide-metabolism-related genes were analyzed. The cytochrome P450-like unigene group was the largest group in DEUs. Most glutathione S-transferase-like unigenes were down-regulated and most sodium channel-like unigenes were up-regulated after insecticide treatment. Finally, 14 insecticide-metabolism-related unigenes were chosen to confirm the relative expression in each treatment by quantitative Real Time Polymerase Chain Reaction (qRT-PCR). The results of qRT-PCR and RNA Sequencing (RNA-Seq) are fairly well-established. Our results demonstrate that a next-generation sequencing tool facilitates the identification of insecticide-metabolism-related genes and the illustration of the insecticide mechanisms of chlorpyrifos and clothianidin.

Transcriptome changes between compatible and incompatible graft combination of Litchi chinensis by digital gene expression profile

Plant grafting has been practiced widely in horticulture and proved as a useful tool in science. However, the mechanisms of graft healing or graft incompatibility remain poorly understood. In this study, Litchi chinensis cv. 'Jingganghongnuo' homograft ('J/J') and 'Jingganghongnuo'/'zhuangyuanhong' heterograft ('J/Z') as compatible and incompatible combination, respectively, was used to study transcriptional changes between incompatible and compatible graft during graft union formation. Anatomical observation indicated that three stages (2 h, 14 d and 21 d after grafting) were critical for graft union formation and selected for high-throughput sequencing. Results indicated 6060 DEGs were differentially expressed in the compatible combination and 5267 DEGs exhibiting in the incompatible one. KEGG pathway enrichment analysis revealed that DEGs were involved in metabolism, wound response, phenylpropanoid biosynthesis and plant hormone signal transduction. The expression of 9 DEGs annotated in auxin pathway was up-regulated in compatible combination than that in incompatible combination. The IAA concentration confirmed that the IAA might promote the graft compatibility. In addition, 13 DEGs related to lignin biosynthesis were differentially expressed during graft healing process. Overall, our results provide abundant sequence resources for studying mechanisms underlying graft compatibility and establish a platform for further studies of litchi and other evergreen fruit trees.

Transcriptome profiling of litchi leaves in response to low temperature reveals candidate regulatory genes and key metabolic events during floral induction

Background

Litchi (Litchi chinensis Sonn.) is an economically important evergreen fruit tree widely cultivated in subtropical areas. Low temperature is absolutely required for floral induction of litchi, but its molecular mechanism is not fully understood. Leaves of litchi played a key role during floral induction and could be the site of low temperature perception. Therefore, leaves were treated under different temperature (15 °C/25 °C), and high-throughput RNA sequencing (RNA-Seq) performed with leaf samples for the de novo assembly and digital gene expression (DGE) profiling analyses to investigate low temperature-induced gene expression changes.

Results

83,107 RNA-Seq unigenes were de novo assembled with a mean length of 1221 bp and approximately 61% of these unigenes (50,345) were annotated against public protein databases. Differentially-expressed genes (DEGs) under low temperature treatment in comparison with the control group were the main focus of our study. Hierarchical clustering analysis arranged 2755 DEGs into eight groups with three significant expression clusters (p-value ≤ 0.05) during floral induction. With the increasing contents of sugars and starch, the expression of genes involved in metabolism of sugars increased dramatically after low temperature induction. One FT gene (Unigene0025396, LcFT1) which produces a protein called ‘florigen’ was also detected among DEGs of litchi. LcFT1 exhibited an apparent specific tissue and its expression was highly increased after low temperature induction, GUS staining results also showed GUS activity driven by LcFT1 gene promoter can be induced by low temperature, which indicated LcFT1 probably played a pivotal role in the floral induction of litchi under low temperature.

Conclusions

Our study provides a global survey of transcriptomes to better understand the molecular mechanisms underlying changes of leaves in response to low temperature induction in litchi. The analyses of transcriptome profiles and physiological indicators will help us study the complicated metabolism of floral induction in the subtropic evergreen plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3747-x) contains supplementary material, which is available to authorized users.

Transcriptional changes during ovule development in two genotypes of litchi (Litchi chinensis Sonn.) with contrast in seed size

Litchi chinensis is a subtropical fruit crop, popular for its nutritional value and taste. Fruits with small seed size and thick aril are desirable in litchi. To gain molecular insight into gene expression that leads to the reduction in the size of seed in Litchi chinensis, transcriptomes of two genetically closely related genotypes, with contrasting seed size were compared in developing ovules. The cDNA library constructed from early developmental stages of ovules (0, 6, and 14 days after anthesis) of bold- and small-seeded litchi genotypes yielded 303,778,968 high quality paired-end reads. These were de-novo assembled into 1,19,939 transcripts with an average length of 865 bp. A total of 10,186 transcripts with contrast in expression were identified in developing ovules between the small- and large- seeded genotypes. A majority of these differences were present in ovules before anthesis, thus suggesting the role of maternal factors in seed development. A number of transcripts indicative of metabolic stress, expressed at higher level in the small seeded genotype. Several differentially expressed transcripts identified in such ovules showed homology with Arabidopsis genes associated with different stages of ovule development and embryogenesis.

Quantitative Proteomics-Based Reconstruction and Identification of Metabolic Pathways and Membrane Transport Proteins Related to Sugar Accumulation in Developing Fruits of Pear (Pyrus communis)

During their 6 month development, pear (Pyrus communis) fruits undergo drastic changes in their morphology and their chemical composition. To gain a better understanding of the metabolic pathways and transport processes active during fruit development, we performed a time-course analysis using mass spectrometry (MS)-based protein identification and quantification of fruit flesh tissues. After pre-fractionation of the samples, 2,841 proteins were identified. A principal component analysis (PCA) separated the samples from seven developmental stages into three distinct clusters representing the early, mid and late developmental phase. Over-representation analysis of proteins characteristic of each developmental phase revealed both expected and novel biological processes relevant at each phase. A high abundance of aquaporins was detected in samples from fruits in the cell expansion stage. We were able quantitatively to reconstruct basic metabolic pathways such as the tricarboxylic acid (TCA) cycle, which indicates sufficient coverage to reconstruct other metabolic pathways. Most of the enzymes that presumably contribute to sugar accumulation in pear fruits could be identified. Our data indicate that invertases do not play a major role in the sugar conversions in developing pear fruits. Rather, sucrose might be broken down by sucrose synthases. Further focusing on sugar transporters, we identified several putative sugar transporters from diverse families which showed developmental regulation. In conclusion, our data set comprehensively describes the proteome of developing pear fruits and provides novel insights about sugar accumulation as well as candidate genes for key reactions and transport steps.

Identification and expression analysis of genes related to calyx persistence in Korla fragrant pear

BACKGROUND

The objective of this study was to increase understanding about genetic mechanisms affecting calyx persistence in Korla fragrant pear (Pyrus brestschneideri Rehd). Flowers were collected at early bloom, full bloom, and late bloom. The RNA was extracted from the flowers and then combined according to calyx type. Transcriptome and digital gene expression (DGE) profiles of flowers, ovaries, and sepals with persistent calyx (SC_hua, SC_ep, and SC_zf, respectively) were compared with those of flowers, ovaries, and sepals with deciduous calyx (TL_hua, TL_ep, and TL_zf, respectively). Temporal changes in the expression of selected genes in floral organs with either persistent or deciduous calyx were compared using real-time quantitative PCR (qRT-PCR).

RESULTS

Comparison of the transcriptome sequences for SC_hua and TL_hua indicated 26 differentially expressed genes (DEGs) with known relationship to abscission and 10 DEGs with unknown function. We identified 98 MYB and 21 SPL genes from the assembled unigenes. From SC_zf vs TL_zf, we identified 21 DEGs with known relationship to abscission and 18 DEGs with unknown function. From SC_ep vs TL_ep, 12 DEGs with known relationship to abscission were identified along with 11 DEGs with unknown function. Ten DEGs were identified by both transcriptome sequencing and DGE sequencing.

CONCLUSIONS

More than 50 DEGs were observed that were related to calyx persistence in Korla fragrant pear. Some of the genes were related to cell wall degradation, plant hormone signal transduction, and stress response. Other DEGs were identified as zinc finger protein genes and lipid transfer protein genes. Further analysis showed that calyx persistence in Korla fragment pear was a metabolic process regulated by many genes related to cell wall degradation and plant hormones.

MYB Transcription Factors in Chinese Pear (Pyrus bretschneideri Rehd.): Genome-Wide Identification, Classification, and Expression Profiling during Fruit Development

The MYB family is one of the largest families of transcription factors in plants. Although, some MYBs were reported to play roles in secondary metabolism, no comprehensive study of the MYB family in Chinese pear (Pyrus bretschneideri Rehd.) has been reported. In the present study, we performed genome-wide analysis of MYB genes in Chinese pear, designated as PbMYBs, including analyses of their phylogenic relationships, structures, chromosomal locations, promoter regions, GO annotations, and collinearity. A total of 129 PbMYB genes were identified in the pear genome and were divided into 31 subgroups based on phylogenetic analysis. These PbMYBs were unevenly distributed among 16 chromosomes (total of 17 chromosomes). The occurrence of gene duplication events indicated that whole-genome duplication and segmental duplication likely played key roles in expansion of the PbMYB gene family. Ka/Ks analysis suggested that the duplicated PbMYBs mainly experienced purifying selection with restrictive functional divergence after the duplication events. Interspecies microsynteny analysis revealed maximum orthology between pear and peach, followed by plum and strawberry. Subsequently, the expression patterns of 20 PbMYB genes that may be involved in lignin biosynthesis according to their phylogenetic relationships were examined throughout fruit development. Among the 20 genes examined, PbMYB25 and PbMYB52 exhibited expression patterns consistent with the typical variations in the lignin content previously reported. Moreover, sub-cellular localization analysis revealed that two proteins PbMYB25 and PbMYB52 were localized to the nucleus. All together, PbMYB25 and PbMYB52 were inferred to be candidate genes involved in the regulation of lignin biosynthesis during the development of pear fruit. This study provides useful information for further functional analysis of the MYB gene family in pear.

Omics studies of citrus, grape and rosaceae fruit trees

Recent advance of bioinformatics and analytical apparatuses such as next generation DNA sequencer (NGS) and mass spectrometer (MS) has brought a big wave of comprehensive study to biology. Comprehensive study targeting all genes, transcripts (RNAs), proteins, metabolites, hormones, ions or phenotypes is called genomics, transcriptomics, proteomics, metabolomics, hormonomics, ionomics or phenomics, respectively. These omics are powerful approaches to identify key genes for important traits, to clarify events of physiological mechanisms and to reveal unknown metabolic pathways in crops. Recently, the use of omics approach has increased dramatically in fruit tree research. Although the most reported omics studies on fruit trees are transcriptomics, proteomics and metabolomics, and a few is reported on hormonomics and ionomics. In this article, we reviewed recent omics studies of major fruit trees, i.e. citrus, grapevine and rosaceae fruit trees. The effectiveness and prospects of omics in fruit tree research will as well be highlighted.

De novo assembly and functional annotation of Myrciaria dubia fruit transcriptome reveals multiple metabolic pathways for L-ascorbic acid biosynthesis

Background

Myrciaria dubia is an Amazonian fruit shrub that produces numerous bioactive phytochemicals, but is best known by its high L-ascorbic acid (AsA) content in fruits. Pronounced variation in AsA content has been observed both within and among individuals, but the genetic factors responsible for this variation are largely unknown. The goals of this research, therefore, were to assemble, characterize, and annotate the fruit transcriptome of M. dubia in order to reconstruct metabolic pathways and determine if multiple pathways contribute to AsA biosynthesis.

Results

In total 24,551,882 high-quality sequence reads were de novo assembled into 70,048 unigenes (mean length = 1150 bp, N50 = 1775 bp). Assembled sequences were annotated using BLASTX against public databases such as TAIR, GR-protein, FB, MGI, RGD, ZFIN, SGN, WB, TIGR_CMR, and JCVI-CMR with 75.2 % of unigenes having annotations. Of the three core GO annotation categories, biological processes comprised 53.6 % of the total assigned annotations, whereas cellular components and molecular functions comprised 23.3 and 23.1 %, respectively. Based on the KEGG pathway assignment of the functionally annotated transcripts, five metabolic pathways for AsA biosynthesis were identified: animal-like pathway, myo-inositol pathway, L-gulose pathway, D-mannose/L-galactose pathway, and uronic acid pathway. All transcripts coding enzymes involved in the ascorbate-glutathione cycle were also identified. Finally, we used the assembly to identified 6314 genic microsatellites and 23,481 high quality SNPs.

Conclusions

This study describes the first next-generation sequencing effort and transcriptome annotation of a non-model Amazonian plant that is relevant for AsA production and other bioactive phytochemicals. Genes encoding key enzymes were successfully identified and metabolic pathways involved in biosynthesis of AsA, anthocyanins, and other metabolic pathways have been reconstructed. The identification of these genes and pathways is in agreement with the empirically observed capability of M. dubia to synthesize and accumulate AsA and other important molecules, and adds to our current knowledge of the molecular biology and biochemistry of their production in plants. By providing insights into the mechanisms underpinning these metabolic processes, these results can be used to direct efforts to genetically manipulate this organism in order to enhance the production of these bioactive phytochemicals.

The accumulation of AsA precursor and discovery of genes associated with their biosynthesis and metabolism in M. dubia is intriguing and worthy of further investigation. The sequences and pathways produced here present the genetic framework required for further studies. Quantitative transcriptomics in concert with studies of the genome, proteome, and metabolome under conditions that stimulate production and accumulation of AsA and their precursors are needed to provide a more comprehensive view of how these pathways for AsA metabolism are regulated and linked in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2225-6) contains supplementary material, which is available to authorized users.

Transcriptome sequencing and analysis of major genes involved in calcium signaling pathways in pear plants (Pyrus calleryana Decne.)

Background

Pears (Pyrus spp. L.) are an important genus of trees that produce one of the world’s oldest fruit crops. Salinity stress is a common limiting factor for plant productivity that significantly affects the flavor and nutritional quality of pear fruits. Much research has shown that calcium signaling pathways, mediated by Calcineurin B-like proteins (CBLs) and their interacting kinases (CIPKs), are closely associated with responses to stresses, including salt. However, little is known about the molecular mechanisms that govern the relationship between salt stress and calcium signaling pathways in pear plants. The available genomic information for pears has promoted much functional genomic analysis and molecular breeding of the genus. This provided an ample foundation for characterizing the transcriptome of pear under salt stress.

Results

A high-throughput Illumina RNA-seq technology was used to identify a total of 78,695 unigenes that were successfully annotated by BLASTX analysis, using the publicly available protein database. Additionally, 2,855 novel transcripts, 218,167 SNPs, 23,248 indels and 18,322 alternative splicing events occurred. Assembled unique sequences were annotated and classified with Gene Ontology (GO), Clusters of Orthologous Group (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, which revealed that the main activated genes in pear are predominately involved in functions such as basic physiological processes, metabolic pathways, operation of cellular components, signal transduction mechanisms, and other molecular activities. Through targeted searches of the annotations, the majority of the genes involved in calcium signaling pathways were identified, among which, four genes were validated by molecular cloning, while 11 were validated by RT-qPCR expression profiles under salt stress treatment.

Conclusions

These results facilitate a better understanding of the molecular genetics and functional genomic mechanisms of salt stress in pear plants. Furthermore, they provide a valuable foundation for additional research on the molecular biology and functional genomics of pear and related species.

Electronic supplementary material

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

Genome-Wide Analysis of Sorbitol Dehydrogenase (SDH) Genes and Their Differential Expression in Two Sand Pear (Pyrus pyrifolia) Fruits

Through RNA-seq of a mixed fruit sample, fourteen expressed sorbitol dehydrogenase (SDH) genes have been identified from sand pear (Pyrus pyrifolia Nakai). Comparative phylogenetic analysis of these PpySDHs with those from other plants supported the closest relationship of sand pear with Chinese white pear (P. bretschneideri). The expression levels varied greatly among members, and the strongest six (PpySDH2, PpySDH4, PpySDH8, PpySDH12, PpySDH13 and PpySDH14) accounted for 96% of total transcript abundance of PpySDHs. Tissue-specific expression of these six members was observed in nine tissues or organs of sand pear, with the greatest abundance found in functional leaf petioles, followed by the flesh of young fruit. Expression patterns of these six PpySDH genes during fruit development were analyzed in two sand pear cultivars, “Cuiguan” and “Cuiyu”. Overall, expression of PpySDHs peaked twice, first at the fruitlet stage and again at or near harvest. The transcript abundance of PpySDHs was higher in “Cuiguan” than in “Cuiyu”, accompanied by a higher content of sugars and higher ratio of fructose to sorbitol maintained in the former cultivar at harvest. In conclusion, it was suggested that multiple members of the SDH gene family are possibly involved in sand pear fruit development and sugar accumulation and may affect both the sugar amount and sugar composition.

Proteome analysis of pear reveals key genes associated with fruit development and quality

Comparative and association analyses of the proteome and transcriptome for pear fruit development were conducted for the first time in this study. Pear fruit development involves complex physiological and biochemical processes, but there is still little knowledge available at proteomic and transcriptomic levels, which would be helpful for understanding the molecular mechanisms of fruit development and quality in pear. In our study, three important stages, including early development (S4-22), middle development (S6-27), and near ripening (S8-30), were investigated in 'Dangshansuli' by isobaric tags for relative and absolute quantitation (iTRAQ) labeling technology, identifying a total of 1,810 proteins during pear fruit development. The association analysis of proteins and transcript expression revealed 1,724, 1,722, and 1,718 associated proteins identified in stages S4-22, S6-27, and S8-30, respectively. A total of 237, 318, and 425 unique proteins were identified as differentially expressed during S4-22 vs S6-27, S6-27 vs S8-30, S4-22 vs S8-30, respectively, and the corresponding correlation coefficients of the overall differentially expressed proteins and transcripts data were 0.6336, 0.4113, and 0.7049. The phenylpropanoid biosynthesis pathway, which is related to lignin formation of pear fruit, was identified as a significantly enriched pathway during early stages of fruit development. Finally, a total of 35 important differentially expressed proteins related to fruit quality were identified, including three proteins related to sugar formation, seven proteins related to aroma synthesis, and sixteen proteins related to the formation of lignin. In addition, qRT-PCR verification provided further evidence to support differentially expressed gene selection. This study is the first to reveal protein and associated mRNA variations in pear during fruit development and quality conformation, and identify key genes and proteins helpful for future functional genomics studies, and provides gene resources for improvement of pear quality.

Transcriptomic analysis of Litchi chinensis pericarp during maturation with a focus on chlorophyll degradation and flavonoid biosynthesis

Background

The fruit of litchi (Litchi chinensis) comprises a white translucent edible aril surrounded by a pericarp. The pericarp of litchi has been the focus of studies associated with fruit size, coloration, cracking and shelf life. However, research at the molecular level has been limited by the lack of genomic and transcriptomic information. In this study, an analysis of the transcriptome of litchi pericarp was performed to obtain information regarding the molecular mechanisms underlying the physiological changes in the pericarp, including those leading to fruit surface coloration.

Results

Coincident with the rapid break down of chlorophyll, but substantial increase of anthocyanins in litchi pericarp as fruit developed, two major physiological changes, degreening and pigmentation were visually apparent. In this study, a cDNA library of litchi pericarp with three different coloration stages was constructed. A total of 4.7 Gb of raw RNA-Seq data was generated and this was then de novo assembled into 51,089 unigenes with a mean length of 737 bp. Approximately 70% of the unigenes (34,705) could be annotated based on public protein databases and, of these, 3,649 genes were significantly differentially expressed between any two coloration stages, while 156 genes were differentially expressed among all three stages. Genes encoding enzymes involved in chlorophyll degradation and flavonoid biosynthesis were identified in the transcriptome dataset. The transcript expression patterns of the Stay Green (SGR) protein suggested a key role in chlorophyll degradation in the litchi pericarp, and this conclusion was supported by the result of an assay over-expressing LcSGR protein in tobacco leaves. We also found that the expression levels of most genes especially late anthocyanin biosynthesis genes were co-ordinated up-regulated coincident with the accumulation of anthocyanins, and that candidate MYB transcription factors that likely regulate flavonoid biosynthesis were identified.

Conclusions

This study provides a large collection of transcripts and expression profiles associated with litchi fruit maturation processes, including coloration. Since most of the unigenes were annotated, they provide a platform for litchi functional genomic research within this species.

Electronic supplementary material

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