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

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.

The PPO family in Nicotiana tabacum is an important regulator to participate in pollination

Polyphenol oxidases (PPOs) are type-3 copper enzymes and are involved in many biological processes. However, the potential functions of PPOs in pollination are not fully understood. In this work, we have screened 13 PPO members in Nicotiana. tabacum (named NtPPO1-13, NtPPOs) to explore their characteristics and functions in pollination. The results show that NtPPOs are closely related to PPOs in Solanaceae and share conserved domains except NtPPO4. Generally, NtPPOs are diversely expressed in different tissues and are distributed in pistil and male gametes. Specifically, NtPPO9 and NtPPO10 are highly expressed in the pistil and mature anther. In addition, the expression levels and enzyme activities of NtPPOs are increased after N. tabacum self-pollination. Knockdown of NtPPOs would affect pollen growth after pollination, and the purines and flavonoid compounds are accumulated in self-pollinated pistil. Altogether, our findings demonstrate that NtPPOs potentially play a role in the pollen tube growth after pollination through purines and flavonoid compounds, and will provide new insights into the role of PPOs in plant reproduction.

Self S-RNase reduces the expression of two pollen-specific COBRA genes to inhibit pollen tube growth in pear

Due to self-incompatibility (SI) prevents self-fertilization, natural or artificial cross-pollination has been conducted in many orchards to stabilize fruit yield. However, it is still puzzled which routes of self S-RNase arresting pollen tube growth. Herein, 17 COBRA genes were isolated from pear genome. Of these genes, the pollen-specifically expressed PbCOB.A.1 and PbCOB.A.2 positively mediates pollen tube growth. The promoters of PbCOB.A.1 and/or PbCOB.A.2 were bound and activated by PbABF.E.2 (an ABRE-binding factor) and PbC2H2.K16.2 (a C2H2-type zinc finger protein). Notably, the expressions of PbCOB.A.1, PbCOB.A.2, and PbC2H2.K16.2 were repressed by self S-RNase, suggesting that self S-RNase reduces the expression of PbCOB.A.1 and PbCOB.A.2 by decreasing the expression of their upstream factors, such as PbC2H2.K16.2, to arrest pollen tube growth. PbCOB.A.1 or PbCOB.A.2 accelerates the growth of pollen tubes treated by self S-RNase, but can hardly affect level of reactive oxygen species and deploymerization of actin cytoskeleton in pollen tubes and cannot physically interact with any reported proteins involved in SI. These results indicate that PbCOB.A.1 and PbCOB.A.2 may not relieve S-RNase toxicity in incompatible pollen tube. The information provides a new route to elucidate the arresting pollen tube growth during SI reaction.

Specific functions of single pistil S-RNases in S-gene homozygous Pyrus germplasm

Gametophytic self-incompatibility (SI) is regulated by S-allele recognition; that is, pollen in a style with the same S-genotype will undergo programmed cell death and stop growing so that it is unable to complete double fertilization, ultimately resulting in the SI response. S-RNase is the female determinant of SI in pear (Pyrus). In the Pyrus genome, there are two different S-RNase alleles at the S-locus, which generate two different S-RNase products in the pistil. The extracted S-glycoprotein is actually a protein complex. In this study, artificial self-pollination was conducted at the bud stage to overcome SI in 'Huanghua' (S1S2) pear. Seven plants homozygous for S1-RNase and four homozygous for S2-RNase were selected from the selfed progeny of 'Huanghua' by S-gene molecular identification biotechnology. We investigated the function of single S-RNases isolated from the pistils of S-gene homozygous Pyrus germplasm. The pollen of 'Huanghua' could smoothly pass through the style of the S-gene homozygous germplasm and complete fertilization. S-RNases were extracted from flower styles of different genotypes and used to treat different types of pollen. The S-RNase from 'Huanghua' completely inhibited the growth of S1S2, S1S1, and S2S2 pollen, while the S-RNase from homozygous germplasm allowed some S1S2 pollen and different single genotypes of pollen to continue growing. These results further validate the core events of SI including cytoskeleton depolymerization and programmed cell death. By iTRAQ-based proteomic analysis of style proteins, a total of 13 S-RNase-related proteins were identified. In summary, we have created reliable S-RNase gene homozygous germplasm, which will play a crucial role in further research on SI in pear and in the development of the pear industry.

Comparative physiochemical and transcriptomic analysis reveals the influences of cross-pollination on ovary and fruit development in pummelo (Citrus maxima)

'Shuijingmiyou' pummelo (SJ), one of the most popular fruits in Yunnan province of China, is of relatively low fruit shape (FS) quality. In this study, we compared the FS promoting effects of cross pollinations using pollens from seven pummelo varieties, and found that 'Guanximiyou' pummelo (GX) cross-pollination showed the best FS promoting effects on SJ fruits by shortening its fruit neck. To explore the underlying mechanism, physiochemical and transcriptomic differences between self- and cross-pollinated SJ ovaries (SJO and GXO) were investigated. Higher salicylic acid, gibberellin and indole acetic acid contents and superoxide dismutase, peroxidase and catalase activities, and lower polyphenol oxidase activity were determined in GXO compared with SJO. Enrichment analysis of the identified 578 differentially expressed genes (123 up-regulated and 455 down-regulated) in GXO showed that genes involved in solute transport, RNA biosynthesis, phytohormone action and cell wall organization were significantly enriched. The results obtained in this study will be helpful in understanding the influences of cross-pollination on pummelo ovary and fruit development, and can provide the basis for clarifying the underlying mechanism of cross-pollination improved fruit quality.

PbrBZR1 interacts with PbrARI2.3 to mediate brassinosteroid-regulated pollen tube growth during self-incompatibility signaling in pear

S-RNase-mediated self-incompatibility (SI) prevents self-fertilization and promotes outbreeding to ensure genetic diversity in many flowering plants, including pear (Pyrus sp.). Brassinosteroids (BRs) have well-documented functions in cell elongation, but their molecular mechanisms in pollen tube growth, especially in the SI response, remain elusive. Here, exogenously applied brassinolide (BL), an active BR, countered incompatible pollen tube growth inhibition during the SI response in pear. Antisense repression of BRASSINAZOLE-RESISTANT1 (PbrBZR1), a critical component of BR signaling, blocked the positive effect of BL on pollen tube elongation. Further analyses revealed that PbrBZR1 binds to the promoter of EXPANSIN-LIKE A3 (PbrEXLA3) to activate its expression. PbrEXLA3 encodes an expansin that promotes pollen tube elongation in pear. The stability of dephosphorylated PbrBZR1 was substantially reduced in incompatible pollen tubes, where it is targeted by ARIADNE2.3 (PbrARI2.3), an E3 ubiquitin ligase that is strongly expressed in pollen. Our results show that during the SI response, PbrARI2.3 accumulates and negatively regulates pollen tube growth by accelerating the degradation of PbrBZR1 via the 26S proteasome pathway. Together, our results show that an ubiquitin-mediated modification participates in BR signaling in pollen and reveal the molecular mechanism by which BRs regulate S-RNase-based SI.

Transcriptome Analysis of the Late-Acting Self-Incompatibility Associated with RNase T2 Family in Camellia oleifera

The Camellia oil tree (Camellia oleifera Abel.) is an important nonwood forest species in China, and the majority of its cultivars are late-acting self-incompatibility (LSI) types. Although several studies have examined the mechanism of LSI, the process is quite complicated and unclear. In this study, pollen tube growth and fruit setting of two Camellia oil tree cultivars Huashuo (HS) and Huajin (HJ) were investigated after non and self-pollination, and transcriptomic analysis of the ovaries was performed 48 h after self-pollination to identify the potential genes implicated in the LSI of Camellia oil trees. The results showed that the fruit set of HS was significantly higher than that of HJ after self-pollination. Transcriptomic analysis revealed that plant hormone signal transduction, the phosphatidylinositol signaling system, ATP-binding cassette (ABC) transporters, reactive oxygen species (ROS) metabolism, and Ca2+ signaling were mainly contributed in the LSI of reaction of Camellia oil tree. Moreover, nine RNase T2 genes were identified from the transcriptome analysis, which also showed that CoRNase7 participated in the self-incompatibility reaction in HS. Based on phylogenetic analysis, CoRNase6 was closely related to S-RNase from coffee, and CoRNase7 and CoRNase8 were closely related to S-RNase from Camellia sinensis. The 9 RNase T2 genes successfully produced proteins in prokaryotes. Subcellular localization indicated that CoRNase1 and CoRNase5 were cytoplasmic proteins, while CoRNase7 was a plasma membrane protein. These results screened the main metabolic pathways closely related to LSI in Camellia oil tree, and SI signal transduction might be regulated by a large molecular regulatory network. The discovery of T2 RNases provided evidence that Camellia oil tree might be under RNase-based gametophytic self-incompatibility.

Integration of semi-in vivo assays and multi-omics data reveals the effect of galloylated catechins on self-pollen tube inhibition in Camellia oleifera

Camellia oil extracted from the seeds of Camellia oleifera Abel. is a popular and high-quality edible oil, but its yield is limited by seed setting, which is mainly caused by self-incompatibility (SI). One of the obvious biological features of SI plants is the inhibition of self-pollen tubes; however, the underlying mechanism of this inhibition in C. oleifera is poorly understood. In this study, we constructed a semi-in vivo pollen tube growth test (SIV-PGT) system that can screen for substances that inhibit self-pollen tubes without interference from the genetic background. Combined with multi-omics analysis, the results revealed the important role of galloylated catechins in self-pollen tube inhibition, and a possible molecular regulatory network mediated by UDP-glycosyltransferase (UGT) and serine carboxypeptidase-like (SCPL) was proposed. In summary, galloylation of catechins and high levels of galloylated catechins are specifically involved in pollen tube inhibition under self-pollination rather than cross-pollination, which provides a new understanding of SI in C. oleifera. These results will contribute to sexual reproduction research on C. oleifera and provide theoretical support for improving Camellia oil yield in production.

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.

A semi in vivo pollination technique to assess the level of gametophytic self-incompatibility and pollen tube growth in pear (Pyrus communis L.)

We describe a semi in vivo pollination technique to determine the compatibility relation between different pear cultivars. This assay provides a valuable addition to existing tools in GSI research. The gametophytic self-incompatibility (GSI) system in Pyrus inhibits fertilization by pollen that shares one of the two S-alleles of the style. Depending on their S-locus genotype, two pear cultivars therefore either show a cross-compatible, semi-compatible or incompatible interaction. Because GSI greatly influences seed and fruit set, accurate knowledge of the compatibility type of a cultivar is key for both pear fruit production and breeding. Currently, compatibility relations between different pear cultivars are generally assessed via S-genotyping. However, this approach is restricted to the currently known S-alleles in pear, and does not provide functional assessment of the level of (self-)incompatibility. We here present an optimized semi in vivo pollination assay, that enables quantitative analysis of (self-)incompatibility in pear, and that can also serve useful for more fundamental studies on pollen tube development and pollen-style interactions. This assay involves in vitro incubation of cut pollinated styles followed by microscopic counting of emerging pollen tubes at a specific time interval. The validity and selectivity of this method to determine compatibility interactions in pear is demonstrated in the cultivars "Celina" and "Packham's Triumph." Overall, this technique constitutes a valuable tool for quantitatively determining in vivo pollen tube growth and (cross-)compatibility in pear.

CoNPR1 and CoNPR3.1 are involved in SA- and MeSA- mediated growth of the pollen tube in Camellia oleifera

Salicylic acid (SA) is a plant hormone involved in a series of growth, development, and stress responses in plants. Nonexpressor of pathogenesis-related genes 1 (NPR1) is the core regulatory gene in the process of SA-mediated systemic acquired resistance (SAR). Whether NPR1 is involved in pollen tube growth mediated by SA and its derivative MeSA (methyl salicylate) remains to be explored. Here, we found that the contents of endogenous SA and MeSA in self- or cross-pollinated pistils changed significantly, and exogenous SA and MeSA significantly promoted pollen germination and pollen tube elongation of Camellia oleifera at lower concentrations. CoNPR1, CoNPR3.1, CoNPR3.2, and CoNPR5 were identified, and they were all located in the nucleus. A high level of consistency was observed across the phylogenetic relationships, gene structures, and functional domains, indicating a clear division of function, as observed in other species. The expression levels of CoNPR1, CoNPR3.1, CoNPR3.2, and CoNPR5 in self- and cross-pollinated pistils had certain regularity. Furthermore, they exhibited tissue-specific expression pattern. CoNPR1 and CoNPR3.1 were expressed in pollen tubes, whose expression was regulated by SA or MeSA, and their expression patterns were basically consistent with the trend of pollen germination. These results indicate that SA and MeSA are involved in the pollen tube growth of C. oleifera, and CoNPRs may play an important role therein.

Molecular characteristics of S-RNase alleles as the determinant of self-incompatibility in the style of Fragaria viridis

Strawberry (Fragaria spp.) is a member of the Rosoideae subfamily in the family Rosaceae. The self-incompatibility (SI) of some diploid species is a key agronomic trait that acts as a basic pollination barrier; however, the genetic mechanism underlying SI control in strawberry remains unclear. Two candidate S-RNases (S_a- and S_b-RNase) identified in the transcriptome of the styles of the self-incompatible Fragaria viridis 42 were confirmed to be SI determinants at the S locus following genotype identification and intraspecific hybridization using selfing progenies. Whole-genome collinearity and RNase T2 family analysis revealed that only an S locus exists in Fragaria; however, none of the compatible species contained S-RNase. Although the results of interspecific hybridization experiments showed that F. viridis (SI) styles could accept pollen from F. mandshurica (self-compatible), the reciprocal cross was incompatible. S_a and S_b-RNase contain large introns, and their noncoding sequences (promotors and introns) can be transcribed into long noncoding RNAs (lncRNAs). Overall, the genus Fragaria exhibits S-RNase-based gametophytic SI, and S-RNase loss occurs at the S locus of compatible germplasms. In addition, a type of SI-independent unilateral incompatibility exists between compatible and incompatible Fragaria species. Furthermore, the large introns and neighboring lncRNAs in S-RNase in Fragaria could offer clues about S-RNase expression strategies.

Genome-wide analysis of long noncoding RNAs affecting floral bud dormancy in pears in response to cold stress

The versatile role of long noncoding RNAs (lncRNAs) in plant growth and development has been established, but a systematic identification and analysis of lncRNAs in the pear has not been reported. Bud dormancy is a crucial and complicated protective mechanism for plants in winter. The roles of lncRNAs in the dormancy process remain largely unclear. In this study, we induced pear floral buds to enter into different dormant statuses by simulating four different chilling accumulation conditions. Then, a time series of RNA-seq analysis was performed and we identified 7594 lncRNAs in Pyrus pyrifolia (Burm. F.) Nakai that have not been identified. The sequence and expression of the lncRNAs were confirmed by PCR analysis. In total, 6253 lncRNAs were predicted to target protein-coding genes including 692 cis-regulated pairs (596 lncRNAs) and 13,158 trans-regulated pairs (6181 lncRNAs). Gene Ontology analysis revealed that most of lncRNAs' target genes were involved in catalytic activity, metabolic processes and cellular processes. In the trend analysis, 124 long-term cold response lncRNAs and 80 short-term cold response lncRNAs were predicted. Regarding the lncRNA-miRNA regulatory networks, 59 lncRNAs were identified as potential precursors for miRNA members of 20 families, 586 lncRNAs were targets of 261 pear miRNAs and 53 lncRNAs were endogenous target mimics for 26 miRNAs. In addition, three cold response lncRNAs, two miRNAs and their target genes were selected for expression confirmed. The trend of their expression was consistent with the predicted relationships among them and suggested possible roles of lncRNAs in ABA metabolic pathway. Our findings not only suggest the potential roles of lncRNAs in regulating the dormancy of pear floral buds but also provide new insights into the lncRNA-miRNA-mRNA regulatory network in plants.

Caspase-like proteases and the phytohormone cytokinin as determinants of S-RNAse-based self-incompatibility-induced PCD in Petunia hybrida L

S-RNAse-based self-incompatibility (SI) in petunia (Petunia hybrida L.) is a self-/non-self-recognition system underlying the pistil rejection of self-pollen. Using different methods, including a TUNEL assay, we have recently shown that programmed cell death (PCD) is a factor of the SI in petunia. Here, we show that the growth of self-incompatible pollen tubes in the style tissues during 4 h after pollination is accompanied by five-sixfold increase in a caspase-like protease (CLP) activity. Exogenous cytokinin (CK) inhibits the pollen tube growth and stimulates the CLP activity in compatible pollen tubes. The actin depolymerization with latrunculin B induces a sharp drop in the CLP activity in self-incompatible pollen tubes and its increase in compatible pollen tubes. Altogether, our results suggest that a CLP is involved in the SI-induced PCD and that CK is a putative activator of the CLP. We assume that CK provokes acidification of the cytosol and thus promotes the activation of a CLP. Thus, our results suggest that CK and CLP are involved in the S-RNAse-based SI-induced PCD in petunia. Potential relations between these components in PCD signaling are discussed.

Using RNA-seq to characterize pollen-stigma interactions for pollination studies

Insects are essential for the reproduction of pollinator-dependent crops and contribute to the pollination of 87% of wild plants and 75% of the world’s food crops. Understanding pollen flow dynamics between plants and pollinators is thus essential to manage and conserve wild plants and ensure yields are maximized in food crops. However, the determination of pollen transfer in the field is complex and laborious. We developed a field experiment in a pollinator-dependent crop and used high throughput RNA sequencing (RNA-seq) to quantify pollen flow by measuring changes in gene expression between pollination treatments across different apple (Malus domestica Borkh.) cultivars. We tested three potential molecular indicators of successful pollination and validated these results with field data by observing single and multiple visits by honey bees (Apis mellifera) to apple flowers and measured fruit set in a commercial apple orchard. The first indicator of successful outcrossing was revealed via differential gene expression in the cross-pollination treatments after 6 h. The second indicator of successful outcrossing was revealed by the expression of specific genes related to pollen tube formation and defense response at three different time intervals in the stigma and the style following cross-pollination (i.e. after 6, 24, and 48 h). Finally, genotyping variants specific to donor pollen could be detected in cross-pollination treatments, providing a third indicator of successful outcrossing. Field data indicated that one or five flower visits by honey bees were insufficient and at least 10 honey bee flower visits were required to achieve a 25% probability of fruit set under orchard conditions. By combining the genotyping data, the differential expression analysis, and the traditional fruit set field experiments, it was possible to evaluate the pollination effectiveness of honey bee visits under orchards conditions. This is the first time that pollen-stigma-style mRNA expression analysis has been conducted after a pollinator visit (honey bee) to a plant (in vivo apple flowers). This study provides evidence that mRNA sequencing can be used to address complex questions related to stigma–pollen interactions over time in pollination ecology.

Transcriptome Analysis of Pyrus betulaefolia Seedling Root Responses to Short-Term Potassium Deficiency

Potassium (K) plays a crucial role in multiple physiological and developmental processes in plants. Its deficiency is a common abiotic stress that inhibits plant growth and reduces crop productivity. A better understanding of the mechanisms involved in plant responses to low K could help to improve the efficiency of K use in plants. However, such responses remain poorly characterized in fruit tree species such as pears (Pyrus sp). We analyzed the physiological and transcriptome responses of a commonly used pear rootstock, Pyrus betulaefolia, to K-deficiency stress (0 mM). Potassium deprivation resulted in apparent changes in root morphology, with short-term low-K stress resulting in rapidly enhanced root growth. Transcriptome analyses indicated that the root transcriptome was coordinately altered within 6 h after K deprivation, a process that continued until 15 d after treatment. Potassium deprivation resulted in the enhanced expression (up to 5-fold) of a putative high-affinity K+ transporter, PbHAK5 (Pbr037826.1), suggesting the up-regulation of mechanisms associated with K+ acquisition. The enhanced root growth in response to K-deficiency stress was associated with a rapid and sustained decrease in the expression of a transcription factor, PbMYB44 (Pbr015309.1), potentially involved in mediating auxin responses, and the increased expression of multiple genes associated with regulating root growth. The concentrations of several phytohormones including indoleacetic acid (IAA), ABA, ETH, gibberellin (GA3), and jasmonic acid (JA) were higher in response to K deprivation. Furthermore, genes coding for enzymes associated with carbon metabolism such as SORBITOL DEHYDROGENASE (SDH) and SUCROSE SYNTHASE (SUS) displayed greatly enhanced expression in the roots under K deprivation, presumably indicating enhanced metabolism to meet the increased energy demands for growth and K+ acquisition. Together, these data suggest that K deprivation in P. betulaefolia results in the rapid re-programming of the transcriptome to enhance root growth and K+ acquisition. These data provide key insights into the molecular basis for understanding low-K-tolerance mechanisms in pears and in other related fruit trees and identifying potential candidates that warrant further analyses.

PbrPOE21 inhibits pear pollen tube growth in vitro by altering apical reactive oxygen species content

Genome-wide identification, tissue-specific expression analysis and functional characterization of selected genes containing the pear Pollen Olea europaea I domain reveal their roles in pollen tube growth. Genes containing the Pollen Olea europaea I (POE) domain play crucial roles in diverse growth and developmental processes. Nevertheless, the specific functions of POE family members in progression of pollen tube growth (PTG) remain uncharacterized. We identified 45 PbrPOE genes in the pear (Pyrus bretschneideri) genome, clustered into seven subclasses. PbrPOE genes contained 1 to 11 exons and 0 to 10 introns, with exon/intron structure mostly conserved within each subclass. Whole-genome duplication has mainly contributed to the duplication pattern of PbrPOE genes in pear. Expression profiles of 45 PbrPOE genes in 12 different pear tissues revealed that six PbrPOE genes (PbrPOE6, 12, 21, 29, 35 and 41) of subclass B were highly expressed during the growth of the pear pollen tube in vitro. PbrPOE21 was selected for further functional analysis on the basis of its high and differential expression pattern in pollen. Antisense oligodeoxynucleotide assays demonstrated that PTG was augmented in vitro when PbrPOE21 expression was significantly inhibited. Moreover, pollen tube length in vitro was reduced when PbrPOE21 was transitorily over-expressed using particle bombardment technology. Exogenous PbrPOE21 recombinant protein inhibited PTG in vitro at an optimum concentration of 1.8 µM. PbrPOE21 also affected reactive oxygen species content in the pear pollen tube apex. We suggest that PbrPOE21 inhibits PTG in vitro by altering apical reactive oxygen species content.

Comparative phosphoproteomic analysis of compatible and incompatible pollination in Brassica napus L

Self-incompatibility (SI) promotes outbreeding and prevents self-fertilization to promote genetic diversity in angiosperms. Several studies have been carried to investigate SI signaling in plants; however, protein phosphorylation and dephosphorylation in the fine-tuning of the SI response remain insufficiently understood. Here, we performed a phosphoproteomic analysis to identify the phosphoproteins in the stigma of self-compatible 'Westar' and self-incompatible 'W-3' Brassica napus lines. A total of 4109 phosphopeptides representing 1978 unique protein groups were identified. Moreover, 405 and 248 phosphoproteins were significantly changed in response to SI and self-compatibility, respectively. Casein kinase II (CK II) phosphorylation motifs were enriched in self-incompatible response and identified 127 times in 827 dominant SI phosphorylation residues. Functional annotation of the identified phosphoproteins revealed the major roles of these phosphoproteins in plant-pathogen interactions, cell wall modification, mRNA surveillance, RNA degradation, and plant hormone signal transduction. In particular, levels of homolog proteins ABF3, BKI1, BZR2/BSE1, and EIN2 were significantly increased in pistils pollinated with incompatible pollens. Abscisic acid and ethephon treatment partially inhibited seed set, while brassinolide promoted pollen germination and tube growth in SI response. Collectively, our results provided an overview of protein phosphorylation during compatible/incompatible pollination, which may be a potential component of B. napus SI responses.

In-depth Understanding of Camellia oleifera Self-incompatibility by Comparative Transcriptome, Proteome and Metabolome

Oil-tea tree (Camellia oleifera) is the most important edible oil tree species in China with late-acting self-incompatibility (LSI) properties. The mechanism of LSI is uncertain, which seriously hinders the research on its genetic characteristics, construction of genetic map, selection of cross breeding parents and cultivar arrangement. To gain insights into the LSI mechanism, we performed cytological, transcriptomic, proteomic and metabolomic studies on self- and cross-pollinated pistils. The studies identified 166,591 transcripts, 6851 proteins and 6455 metabolites. Transcriptomic analysis revealed 1197 differentially expressed transcripts between self- and cross-pollinated pistils and 47 programmed cell death (PCD)-control transcripts. Trend analysis by Pearson correlation categorized nine trend graphs linked to 226 differentially expressed proteins and 38 differentially expressed metabolites. Functional enrichment analysis revealed that the LSI was closely associated with PCD-related genes, mitogen-activated protein kinase (MAPK) signaling pathway, plant hormone signal transduction, ATP-binding cassette (ABC) transporters and ubiquitin-mediated proteolysis. These particular trends in transcripts, proteins and metabolites suggested the involvement of PCD in LSI. The results provide a solid genetic foundation for elucidating the regulatory network of PCD-mediated self-incompatibility in C. oleifera.

RNA-Seq analysis of compatible and incompatible styles of Pyrus species at the beginning of pollination

At the early stage of pollination, the difference in gene expression between compatibility and incompatibility is highly significant about the pollen-specific expression of the LRR gene, resistance, and defensin genes. In Rosaceae, incompatible pollen can penetrate into the style during the gametophytic self-incompatibility response. It is therefore considered a stylar event rather than a stigmatic event. In this study, we explored the differences in gene expression between compatibility and incompatibility in the early stage of pollination. The self-compatible pear variety "Jinzhuili" is a naturally occurring bud mutant from "Yali", a leading Chinese native cultivar exhibiting typical gametophytic self-incompatibility. We collected the styles of 'Yali' and 'Jinzhuili' at 0.5 and 2 h after self-pollination and then performed high-throughput sequencing. According to the KEGG analysis of the differentially expressed genes, several metabolic pathways, such as "Plant hormone signal transduction", "Plant-pathogen interaction", are the main pathways was the most represented pathway. Quantitative PCR was used to validate these differential genes. The expression levels of genes related to pollen growth and disease inhibition, such as LRR (Leucine-rich repeat extensin), resistance, defensin, and auxin, differed significantly between compatible and incompatible pollination. Interestingly, at 0.5 h, most of these genes were upregulated in the compatible pollination system compared with the incompatible pollination system. Calcium transport, which requires ATPase, also demonstrated upregulated expression. In summary, the self-incompatibility reaction was initiated when the pollen land on the stigma.

A gamma-thionin protein from apple, MdD1, is required for defence against S-RNase-induced inhibition of pollen tube prior to self/non-self recognition

Apple exhibits S-RNase-mediated self-incompatibility. Although the cytotoxic effect of S-RNase inside the self-pollen tube has been studied extensively, the underlying defence mechanism in pollen tube in Rosaceae remains unclear. On exposure to stylar S-RNase, plant defence responses are activated in the pollen tube; however, how these are regulated is currently poorly understood. Here, we show that entry of both self and non-self S-RNase into pollen tubes of apple (Malus domestica) stimulates jasmonic acid (JA) production, in turn inducing the accumulation of MdMYC2 transcripts, a transcription factor in the JA signalling pathway widely considered to be involved in plant defence processes. MdMYC2 acts as a positive regulator in the pollen tube activating expression of MdD1, a gene encoding a defence protein. Importantly, MdD1 was shown to bind to the RNase activity sites of S-RNase leading to inhibition of enzymatic activity. This work provides intriguing insights into an ancient defence mechanism present in apple pollen tubes where MdD1 likely acts as a primary line of defence to inhibit S-RNase cytotoxicity prior to self/non-self recognition.

Identifying Genome-Wide Sequence Variations and Candidate Genes Implicated in Self-Incompatibility by Resequencing Fragaria viridis

It is clear that the incompatibility system in Fragaria is gametophytic, however, the genetic mechanism behind this remains elusive. Eleven second-generation lines of Fragaria viridis with different compatibility were obtained by manual self-pollination, which can be displayed directly by the level of fruit-set rate. We sequenced two second-generation selfing lines with large differences in fruit-set rate: Ls-S₂-53 as a self-incompatible sequencing sample, and Ls-S₂-76 as a strong self-compatible sequencing sample. Fragaria vesca was used as a completely self-compatible reference sample, and the genome-wide variations were identified and subsequently annotated. The distribution of polymorphisms is similar on each chromosome between the two sequencing samples, however, the distribution regions and the number of homozygous variations are inconsistent. Expression pattern analysis showed that six candidate genes were significantly associated with self-incompatibility. Using F. vesca as a reference, we focused our attention on the gene FIP2-like (FH protein interacting protein), associated with actin cytoskeleton formation, as the resulting proteins in Ls-S₂-53 and Ls-S₂-76 have each lost a number of different amino acids. Suppression of FIP2-like to some extent inhibits germination of pollen grains and growth of pollen tubes by reducing F-actin of the pollen tube tips. Our results suggest that the differential distribution of homozygous variations affects F. viridis fruit-set rate and that the fully encoded FIP2-like can function normally to promote F-actin formation, while the new FIP2-like proteins with shortened amino acid sequences have influenced the (in)compatibility of two selfing lines of F. viridis.

Global Transcriptional Insights of Pollen-Pistil Interactions Commencing Self-Incompatibility and Fertilization in Tea [Camellia sinensis (L.) O. Kuntze]

This study explicates molecular insights commencing Self-Incompatibility (SI) and CC (cross-compatibility/fertilization) in self (SP) and cross (CP) pollinated pistils of tea. The fluorescence microscopy analysis revealed ceased/deviated pollen tubes in SP, while successful fertilization occurred in CP at 48 HAP. Global transcriptome sequencing of SP and CP pistils generated 109.7 million reads with overall 77.9% mapping rate to draft tea genome. Furthermore, concatenated de novo assembly resulted into 48,163 transcripts. Functional annotations and enrichment analysis (KEGG & GO) resulted into 3793 differentially expressed genes (DEGs). Among these, de novo and reference-based expression analysis identified 195 DEGs involved in pollen-pistil interaction. Interestingly, the presence of 182 genes [PT germination & elongation (67), S-locus (11), fertilization (43), disease resistance protein (30) and abscission (31)] in a major hub of the protein-protein interactome network suggests a complex signaling cascade commencing SI/CC. Furthermore, tissue-specific qRT-PCR analysis affirmed the localized expression of 42 DE putative key candidates in stigma-style and ovary, and suggested that LSI initiated in style and was sustained up to ovary with the active involvement of csRNS, SRKs & SKIPs during SP. Nonetheless, COBL10, RALF, FERONIA-rlk, LLG and MAPKs were possibly facilitating fertilization. The current study comprehensively unravels molecular insights of phase-specific pollen-pistil interaction during SI and fertilization, which can be utilized to enhance breeding efficiency and genetic improvement in tea.