Dissection of the style's response to pollination using transcriptome profiling in self-compatible (Solanum pimpinellifolium) and self-incompatible (Solanum chilense) tomato species

The influence of the pollination compatibility type on the pistil S-RNase expression in European pear (Pyrus communis)

The S-RNase gene plays an essential role in the gametophytic self-incompatibility (GSI) system of Pyrus. It codes for the stylar-expressed S-RNase protein which inhibits the growth of incompatible pollen tubes through cytotoxicity and the induction of programmed cell death in the pollen tube. While research on the Pyrus GSI system has primarily focused on the S-RNase gene, there is still a lack of insight into its spatiotemporal expression profile and the factors that regulate it. Previous studies have suggested that S-RNase expression in the style is influenced by pollination and is dependent on the compatibility type. We here continue on this basic hypothesis by analyzing the spatiotemporal expression of the S-RNase alleles in Pyrus communis "Conference" styles in response to different types of pollination; namely, upon full- and semi-compatible pollination and upon incompatible selfing. The results revealed that temporal dynamics of S-RNase expression are influenced by the pollen's compatibility type, indicating the presence of a signaling mechanism between pollen and style to control S-RNase production during pollen tube growth. In our experiment, S-RNase expression continuously decreased after cross-pollination and in the unpollinated control. However, after a fully incompatible pollination, S-RNase expression remained constant. Finally, semi-compatible pollination showed a initially constant S-RNase expression for both alleles followed by a strong decrease in expression. Based on these results and previous findings, we propose a regulatory mechanism to explain the effect of pollination and the associated compatibility type on S-RNase expression in the style. This proposed mechanism could be used as a starting point for future research.

Comparative stigmatic transcriptomics reveals self and cross pollination responses to heteromorphic incompatibility in Plumbago auriculata Lam

"Heteromorphic self-incompatibility" (HetSI) in plants is a mechanism of defense to avoid self-pollination and promote outcrossing. However, the molecular mechanism underlying HetSI remains largely unknown. In this study, RNA-seq was conducted to explore the molecular mechanisms underlying self-compatible (SC, "T × P" and "P × T") and self-incompatible (SI, "T × T" and "P × P") pollination in the two types of flowers of Plumbago auriculata Lam. which is a representative HetSI plant. By comparing "T × P" vs. "T × T", 3773 (1407 upregulated and 2366 downregulated) differentially expressed genes (DEGs) were identified, 1261 DEGs between "P × T" and "P × P" (502 upregulated and 759 downregulated). The processes in which these DEGs were significantly enriched were "MAPK (Mitogen-Activated Protein Kinases-plant) signaling pathway", "plant-pathogen interaction","plant hormone signal transduction", and "pentose and glucuronate interconversion" pathways. Surprisingly, we discovered that under various pollination conditions, multiple notable genes that may be involved in HetSI exhibited distinct regulation. We can infer that the HetSI strategy might be unique in P. auriculata. It was similar to "sporophytic self-incompatibility" (SSI) but the HetSI mechanisms in pin and thrum flowers are diverse. In this study, new hypotheses and inferences were proposed, which can provide a reference for crop production and breeding.

Genome-wide transcriptome analysis uncovers gene networks regulating fruit quality and volatile compounds in mango cultivar 'Tainong' during postharvest

Mango is one of the most economically important fruit; however, the gene regulatory mechanism associated with ripening and quality changes during storage remains largely unclear. This study explored the relationship between transcriptome changes and postharvest mango quality. Fruit quality patterns and volatile components were obtained using headspace gas chromatography and ion-mobility spectrometry (HS-GC-IMS). The changes in mango peel and pulp transcriptome were analyzed during four stages (pre-harvesting, harvesting, maturity, and overripe stages). Based on the temporal analysis, multiple genes involved in the biosynthesis of secondary metabolites were upregulated in both the peel and pulp during the mango ripening process. Moreover, cysteine and methionine metabolism related to ethylene synthesis were upregulated in the pulp over time. Weighted gene co-expression network analysis (WGCNA) further showed that the pathways of pyruvate metabolism, citrate cycle, propionate metabolism, autophagy, and SNARE interactions in vesicular transport were positively correlated with the ripening process. Finally, a regulatory network of important pathways from pulp to peel was constructed during the postharvest storage of mango fruit. The above findings provide a global insight into the molecular regulation mechanisms of postharvest mango quality and flavor changes.

Gene coexpression analysis reveals key pathways and hub genes related to late-acting self-incompatibility in Camellia oleifera

INTRODUCTION

Self-incompatibility (SI) is an important strategy for plants to maintain abundant variation to enhance their adaptability to the environment. Camellia oleifera is one of the most important woody oil plants and is widely cultivated in China. Late acting self-incompatibility (LSI) in C. oleifera results in a relatively poor fruit yield in the natural state, and understanding of the LSI mechanism remains limited.

METHODS

To better understand the molecular expression and gene coexpression network in the LSI reaction in C. oleifera, we conducted self- and cross-pollination experiments at two different flower bud developmental stages (3-4 d before flowering and 1 d before flowering), and cytological observation, fruit setting rate (FSR) investigation and RNA-Seq analysis were performed to investigate the mechanism of the male -female interaction and identify hub genes responsible for the LSI in C. oleifera.

RESULTS

Based on the 21 ovary transcriptomes, a total of 7669 DEGs were identified after filtering out low-expression genes. Weighted gene coexpression network analysis (WGCNA) divided the DEGs into 15 modules. Genes in the blue module (1163 genes) were positively correlated with FSR, and genes in the pink module (339 genes) were negatively correlated with FSR. KEGG analysis indicated that flavonoid biosynthesis, plant MAPK signaling pathways, ubiquitin-mediated proteolysis, and plant-pathogen interaction were the crucial pathways for the LSI reaction. Fifty four transcription factors (TFs) were obtained in the two key modules, and WRKY and MYB were potentially involved in the LSI reaction in C. oleifera. Network establishment indicated that genes encoding G-type lectin S-receptor-like serine (lecRLK), isoflavone 3'-hydroxylase-like (CYP81Q32), cytochrome P450 87A3-like (CYP87A3), and probable calcium-binding protein (CML41) were the hub genes that positively responded to the LSI reaction. The other DEGs inside the two modules, including protein RALF-like 10 (RALF), F-box and pectin acetylesterase (MTERF5), might also play vital roles in the LSI reaction in C. oleifera.

DISCUSSION

Overall, our study provides a meaningful resource for gene network studies of the LSI reaction process and subsequent analyses of pollen-pistil interactions and TF roles in the LSI reaction, and it also provides new insights for exploring the mechanisms of the LSI response.

Combined transcriptomic and proteomic analysis reveals multiple pathways involved in self-pollen tube development and the potential roles of FviYABBY1 in self-incompatibility in Fragaria viridis

Fragaria viridis exhibits S-RNase-based gametophytic self-incompatibility, in which S-RNase is the major factor inhibiting pollen tube growth. However, the pathways involved in and the immediate causes of the inhibition of pollen tube growth remain unknown. Here, interactive RNA sequencing and proteome analysis revealed changes in the transcriptomic and proteomic profiles of F. viridis styles harvested at 0 and 24 h after self-pollination. A total of 2,181 differentially expressed genes and 200 differentially abundant proteins were identified during the pollen development stage of self-pollination. Differentially expressed genes and differentially abundant proteins associated with self-incompatible pollination were further mined, and multiple pathways were found to be involved. Interestingly, the expression pattern of the transcription factor FviYABBY1, which is linked to polar growth, differed from those of other genes within the same family. Specifically, FviYABBY1 expression was extremely high in pollen, and its expression trend in self-pollinated styles was consistent with that of S-RNase. Furthermore, FviYABBY1 interacted with S-RNase in a non-S haplotype way. Therefore, FviYABBY1 affects the expression of polar growth-related genes in self-pollen tubes and is positively regulated by S-RNase.

Stigmatic Transcriptome Analysis of Self-Incompatible and Compatible Pollination in Corylus heterophylla Fisch. × Corylus avellana L

Self-incompatibility (SI) protects plants from inbreeding depression due to self-pollination and promotes the outcrossing process to maintain a high degree of heterozygosity during evolution. Corylus is an important woody oil and nut species that shows sporophytic SI (SSI). Yet the molecular mechanism of SI in Corylus remains largely unknown. Here we conducted self- ("Dawei" × "Dawei") and cross-pollination ("Dawei" × "Liaozhen No. 7") experiments and then performed an RNA-Seq analysis to investigate the mechanism of pollen-stigma interactions and identify those genes that may be responsible for SSI in Corylus. We uncovered 19,163 up- and 13,314 downregulated genes from the comparison of different pollination treatments. These differentially expressed genes (DEGs) were significantly enriched in plant-pathogen interaction, plant hormone signal transduction, and MAPK signaling pathway-plant. We found many notable genes potentially involved in pollen-stigma interactions and SSI mechanisms, including genes encoding receptor-like protein kinases (RLK), calcium-related genes, disease-resistance genes, and WRKY transcription factors. Four upregulated and five downregulated DEGs were consistently identified in those comparison groups involving self-incompatible pollination, suggesting they had important roles in pollen-pistil interactions. We further identified the S-locus region of the Corylus heterophylla genome based on molecular marker location. This predicted S-locus contains 38 genes, of which 8 share the same functional annotation as the S-locus genes of Corylus avellana: two PIX7 homologous genes (EVM0002129 and EVM0025536), three MIK2 homologous genes (EVM0002422, EVM0005666, and EVM0009820), one aldose 1-epimerase (EVM0002095), one 3-dehydroquinate synthase II (EVM0021283), and one At3g28850 homologous gene (EVM0016149). By characterizing the pistil process during the early postpollination phase via transcriptomic analysis, this study provides new knowledge and lays the foundation for subsequent analyses of pollen-pistil interactions.

Tomato SlIDA has a critical role in tomato fertilization by modifying reactive oxygen species homeostasis

Anther development and pollen tube elongation are key steps for pollination and fertilization. The timing and spatial distribution of reactive oxygen species (ROS) and programmed cell death are central to these processes, but the regulatory mechanism of ROS production is not well understood. Inflorescence deficient in abscission (IDA) is implicated in many plant development and responses to environmental stimuli. However, their role in reproductive development is still unknown. We generated tomato knockout lines (CR-slida) of an IDA homolog (SlIDA), which is expressed in the tapetum, septum and pollen tube, and observed a severe defect in male gametes. Further analysis indicated that there was a programmed cell death defect in the tapetum and septum and a failure of anther dehiscence in the CR-slida lines, likely related to insufficient ROS signal. Liquid chromatography-tandem mass spectrometry identified mature SlIDA as a 14-mer EPIP peptide, which was shown to be secreted, and a complementation experiment showed that application of a synthetic 14-mer EPIP peptide rescued the CR-slida defect and enhanced the ROS signal. Moreover, the application of the ROS scavengers diphenyleneiodonium or Mn-TMPP suppressed peptide function. Collectively, our results revealed that SlIDA plays an essential role in pollen development and pollen tube elongation by modulating ROS homeostasis.

Pollination Type Recognition from a Distance by the Ovary Is Revealed Through a Global Transcriptomic Analysis

Sexual reproduction in flowering plants involves intimate contact and continuous interactions between the growing pollen tube and the female reproductive structures. These interactions can trigger responses in distal regions of the flower well ahead of fertilization. While pollination-induced petal senescence has been studied extensively, less is known about how pollination is perceived at a distance in the ovary, and how specific this response is to various pollen genotypes. To address this question, we performed a global transcriptomic analysis in the ovary of a wild potato species, Solanum chacoense, at various time points following compatible, incompatible, and heterospecific pollinations. In all cases, pollen tube penetration in the stigma was initially perceived as a wounding aggression. Then, as the pollen tubes grew in the style, a growing number of genes became specific to each pollen genotype. Functional classification analyses revealed sharp differences in the response to compatible and heterospecific pollinations. For instance, the former induced reactive oxygen species (ROS)-related genes while the latter affected genes associated to ethylene signaling. In contrast, incompatible pollination remained more akin to a wound response. Our analysis reveals that every pollination type produces a specific molecular signature generating diversified and specific responses at a distance in the ovary in preparation for fertilization.

Transcriptomic Analysis Reveals a Comprehensive Calcium- and Phytohormone-Dominated Signaling Response in Leymus chinensis Self-Incompatibility

Sheepgrass (Leymus chinensis (Trin.) Tzvel.) is an economically and ecologically important forage in the grass family. Self-incompatibility (SI) limits its seed production due to the low seed-setting rate after self-pollination. However, investigations into the molecular mechanisms of sheepgrass SI are lacking. Therefore, microscopic observation of pollen germination and pollen tube growth, as well as transcriptomic analyses of pistils after self- and cross-pollination, were performed. The results indicated that pollen tube growth was rapidly inhibited from 10 to 30 min after self-pollination and subsequently stopped but preceded normally after cross-pollination. Time course comparative transcriptomics revealed different transcriptome dynamics between self- and cross-pollination. A pool of SI-related signaling genes and pathways was generated, including genes related to calcium (Ca2+) signaling, protein phosphorylation, plant hormone, reactive oxygen species (ROS), nitric oxide (NO), cytoskeleton, and programmed cell death (PCD). A putative SI response molecular model in sheepgrass was presented. The model shows that SI may trigger a comprehensive calcium- and phytohormone-dominated signaling cascade and activate PCD, which may explain the rapid inhibition of self-pollen tube growth as observed by cytological analyses. These results provided new insight into the molecular mechanisms of sheepgrass (grass family) SI.

Transcriptome and phytohormone analysis reveals a comprehensive phytohormone and pathogen defence response in pear self-/cross-pollination

Candidate genes were identified and the role of phytohormones such as JA-Me and ABA in the synthesis of S-RNase was emphasized in pear self-incompatibility. Self-incompatibility (SI) occurs widely in flowering plants as an intraspecific reproductive barrier. This phenomenon promotes variation within species, but for some species such as Pyrus, SI is a nuisance rather than a benefit in agricultural production. Although many studies have been conducted on SI in pears, its mechanism remains unclear. In this study, high-throughput Illumina RNA sequencing (RNA-seq) was used to identify SI-related genes in pear styles. Using transcriptome comparisons, differentially expressed genes of unpollinated (UP), cross-pollinated (CP), and self-pollinated (SP) styles were identified after 48 h. A total of 1796 and 1890 genes were identified in DSC (UP vs. CP) and DSI (UP vs. SP), respectively. KEGG analysis revealed that genes involved in the "plant hormone signal transduction pathway" and "plant-pathogen interaction pathway" were significantly enriched in DSI (UP vs. SP) compared to those in DSC (UP vs. CP). The expression level of S-glycoprotein ribonuclease (S-RNase) was dramatically reduced in cross-pollinated (CP) styles. To better understand the relationship between the expression patterns of S-RNase and two major KEGG pathways, the concentrations of phytohormones were measured, and the expression pattern of S-RNase was analysed using qRT-PCR. Our results demonstrate that methyl jasmonate and abscisic acid may enhance the expression level of S-RNase, and pollination can affect the synthesis of methyl jasmonate and abscisic acid in pear styles. Overall, this study is a global transcriptome analysis of SI in pear. A relationship between self-rejection, plant hormones, and pathogen defence was shown in pear.

Dissecting stylar responses to self-pollination in wild tomato self-compatible and self-incompatible species using comparative proteomics

Self-incompatibility (SI), a phenomenon that is widespread among flowering plants (angiosperms), promotes outbreeding, resulting in increased genetic diversity and species survival. SI is also important in establishing intra- or interspecies reproductive barriers, such as those that are evident in the tomato clade, Solanum section Lycopersicon, where they limit the use of wild species inbreeding programs to improve cultivated tomato. However, the molecular mechanisms underlying SI are poorly understood in the tomato clade. In this study, an SI (Solanum chilense, LA0130) and a self-compatible (SC, Solanum pimpinellifolium, LA1585) tomato species were chosen to dissect the mechanism of SI formation using a comparative proteomics approach. A total of 635 and 627 protein spots were detected in two-dimensional electrophoresis (2-DE) maps of proteins from the SI and SC species, respectively. In the SC species, 22 differently expressed proteins (DEPs) were detected in SCP versus SCUP (self-pollination versus non-pollination in SC species). Of these, 3 and 18 showed an up-or down-regulated expression in the SCP protein sample, respectively, while only one DEP (MSRA, Solyc03g111720) was exclusively expressed in the SCP sample. In the SI species, 14 DEPs were found between SIP/SIUP, and 5 of these showed higher expression in SIP, whereas two DEPs (MLP-like protein 423-like, gene ID, 460386008 and (ATP synthase subunit alpha, gene ID, Solyc00g042130) were exclusively expressed in SIP or SIUP, respectively. Finally, two S-RNases (gene IDs, 313247946 and 157377662) were exclusively expressed in the SI species. Sequence homology analysis and a gene ontology tool were used to assign the DEPs to the 'metabolism', 'energy', 'cytoskeleton dynamics', 'protein degradation', 'signal transduction', 'defence/stress responses', 'self-incompatibility' and 'unknown' protein categories. We discuss the putative functions of the DEPs in different biological processes and how these might be associated with the regulation of SI formation in the tomato clade.

Transcriptome analysis reveals self-incompatibility in the tea plant (Camellia sinensis) might be under gametophytic control

BACKGROUND

Self-incompatibility (SI) is under genetic control and prevents inbreeding depression in angiosperms. SI mechanisms are quite complicated and still poorly understood in many plants. Tea (Camellia sinensis L.) belonging to the family of Theaceae, exhibits high levels of SI and high heterozygosity. Uncovering the molecular basis of SI of the tea plant may enhance breeding and simplify genomics research for the whole family.

RESULTS

The growth of pollen tubes following selfing and crossing was observed using fluorescence microscopy. Self-pollen tubes grew slower than cross treatments from 24 h to 72 h after pollination. RNA-seq was employed to explore the molecular mechanisms of SI and to identify SI-related genes in C. sinensis. Self and cross-pollinated styles were collected at 24 h, 48 h and 72 h after pollination. Six RNA-seq libraries (SP24, SP48, SP72, CP24 CP48 and CP72; SP = self-pollinated, CP = cross-pollinated) were constructed and separately sequenced. In total, 299.327 million raw reads were generated. Following assembly, 63,762 unigenes were identified, and 27,264 (42.76 %) unigenes were annotated in five public databases: NR, KOG, KEGG, Swiss-Port and GO. To identify SI-related genes, the fragments per kb per million mapped reads (FPKM) values of each unigene were evaluated. Comparisons of CP24 vs. SP24, CP48 vs. SP48 and CP72 vs. SP72 revealed differential expression of 3,182, 3,575 and 3,709 genes, respectively. Consequently, several ubiquitin-mediated proteolysis, Ca(2+) signaling, apoptosis and defense-associated genes were obtained. The temporal expression pattern of genes following CP and SP was analyzed; 6 peroxidase, 1 polyphenol oxidase and 7 salicylic acid biosynthetic process-related genes were identified. The RNA-seq data were validated by qRT-PCR of 15 unigenes. Finally, a unigene (CL25983Contig1) with strong homology to the S-RNase was analyzed. It was mainly expressed in styles, with dramatically higher expression in self-pollinated versus cross-pollinated tissues at 24 h post-pollination.

CONCLUSIONS

The present study reports the transcriptome of styles after cross- and self-pollination in tea and offers novel insights into the molecular mechanism behind SI in C. sinensis. We believe that this RNA-seq dataset will be useful for improvement in C. sinensis as well as other plants in the Theaceae family.

Gene Expression Profiles in Rice Developing Ovules Provided Evidence for the Role of Sporophytic Tissue in Female Gametophyte Development

The development of ovule in rice (Oryza sativa) is vital during its life cycle. To gain more understanding of the molecular events associated with the ovule development, we used RNA sequencing approach to perform transcriptome-profiling analysis of the leaf and ovules at four developmental stages. In total, 25,401, 23,343, 23,647 and 23,806 genes were identified from the four developmental stages of the ovule, respectively. We identified a number of differently expressed genes (DEGs) from three adjacent stage comparisons, which may play crucial roles in ovule development. The DEGs were then conducted functional annotations and Kyoto encyclopedia of genes and genomes (KEGG) pathway analyses. Genes related to cellular component biogenesis, membrane-bounded organelles and reproductive regulation were identified to be highly expressed during the ovule development. Different expression levels of auxin-related and cytokinin-related genes were also identified at various stages, providing evidence for the role of sporophytic ovule tissue in female gametophyte development from the aspect of gene expression. Generally, an overall transcriptome analysis for rice ovule development has been conducted. These results increased our knowledge of the complex molecular and cellular events that occur during the development of rice ovule and provided foundation for further studies on rice ovule development.

Dissection of the style's response to pollination using transcriptome profiling in self-compatible (Solanum pimpinellifolium) and self-incompatible (Solanum chilense) tomato species

BACKGROUND

Tomato (Solanum lycopersicum) self-compatibility (SC) is defined as self-pollen tubes that can penetrate their own stigma, elongate in the style and fertilize their own ovules. Self-incompatibility (SI) is defined as self-pollen tubes that are prevented from developing in the style. To determine the influence of gene expression on style self-pollination, a transcriptome-wide comparative analysis of SC and SI tomato unpollinated/pollinated styles was performed using RNA-sequencing (RNA-seq) data.

RESULTS

Transcriptome profiles of 24-h unpollination (UP) and self-pollination (P) styles from SC and SI tomato species were generated using high-throughput next generation sequencing. From the comparison of SC self-pollinated and unpollinated styles, 1341 differentially expressed genes (DEGs) were identified, of which 753 were downregulated and 588 were upregulated. From the comparison of SI self-pollinated and unpollinated styles, 804 DEGs were identified, of which 215 were downregulated and 589 were upregulated. Nine gene ontology (GO) terms were enriched significantly in SC and 78 GO terms were enriched significantly in SI. A total of 105 enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified in SC and 80 enriched KEGG pathways were identified in SI, among which "Cysteine and methionine metabolism pathway" and "Plant hormone signal transduction pathway" were significantly enriched in SI.

CONCLUSIONS

This study is the first global transcriptome-wide comparative analysis of SC and SI tomato unpollinated/pollinated styles. Advanced bioinformatic analysis of DEGs uncovered the pathways of "Cysteine and methionine metabolism" and "Plant hormone signal transduction", which are likely to play important roles in the control of pollen tubes growth in SI species.