Characterization of the 'Xiangshui' lemon transcriptome by de novo assembly to discover genes associated with self-incompatibility

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

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 origin and the genetic regulation of the self-compatibility mechanism in clementine (Citrus clementina Hort. ex Tan.)

Self-incompatibility (SI) is a genetic mechanism common in flowering plants to prevent self-fertilization. Among citrus species, several pummelo, mandarin, and mandarin-like accessions show SI behavior. In these species, SI is coupled with a variable degree of parthenocarpy ensuring the production of seedless fruits, a trait that is highly appreciated by consumers. In Citrus, recent evidences have shown the presence of a gametophytic SI system based on S-ribonucleases (S-RNases) ability to impair self-pollen tube growth in the upper/middle part of the style. In the present study, we combined PCR analysis and next-generation sequencing technologies, to define the presence of S7- and S11-Rnases in the S-genotype of the Citrus clementina (Hort. ex Tan.), the self-incompatible 'Comune' clementine and its self-compatible natural mutant 'Monreal'. The reference genome of 'Monreal' clementine is presented for the first time, providing more robust results on the genetic sequence of the newly discovered S7-RNase. SNP discovery analysis coupled with the annotation of the variants detected enabled the identification of 7,781 SNPs effecting 5,661 genes in 'Monreal' compared to the reference genome of C. clementina. Transcriptome analysis of unpollinated pistils at the mature stage from both clementine genotypes revealed the lack of expression of S7-RNase in 'Monreal' suggesting its involvement in the loss of the SI response. RNA-seq analysis followed by gene ontology studies enabled the identification of 2,680 differentially expressed genes (DEGs), a significant number of those is involved in oxidoreductase and transmembrane transport activity. Merging of DNA sequencing and RNA data led to the identification of 164 DEGs characterized by the presence of at least one SNP predicted to induce mutations with a high effect on their amino acid sequence. Among them, four candidate genes referring to two Agamous-like MADS-box proteins, to MYB111 and to MLO-like protein 12 were validated. Moreover, the transcription factor MYB111 appeared to contain a binding site for the 2.0-kb upstream sequences of the S7- and S11-RNase genes. These results provide useful information about the genetic bases of SI indicating that SNPs present in their sequence could be responsible for the differential expression and the regulation of S7-RNase and consequently of the SI mechanism.

The Pitaya Flower Tissue's Gene Differential Expression Analysis between Self-Incompatible and Self-Compatible Varieties for the Identification of Genes Involved in Self-Incompatibility Regulation

Self-incompatible pitaya varieties have low fruit-setting rates under natural conditions, leading to higher production costs and hindering industrial prosperity. Through transcriptome sequencing, we obtained the 36,900 longest transcripts (including 9167 new transcripts) from 60 samples of flowers. Samples were collected pre- and post-pollination (at 0 h, 0.5 h, 2 h, 4 h, and 12 h) from two varieties of pitaya (self-compatible Jindu No. 1 and self-incompatible Cu Sha). Using the RNA-Seq data and comparison of reference genomes, we annotated 28,817 genes in various databases, and 1740 genes were optimized in their structure for annotation. There were significant differences in the expression of differentially expressed genes (DEGs) in the pitaya stigmas under different pollination types, especially at the late post-pollination stage, where the expression of protease genes increasedal significantly under cross-pollination. We identified DEGs involved in the ribosomal, ubiquitination-mediated, and phyto-signaling pathways that may be involved in pitaya SI regulation. Based on the available transcriptome data and bioinformatics analysis, we tentatively identified HuS-RNase2 as a candidate gynogenetic S gene in the pitaya GSI system.

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.

Genome-Wide Analysis of the RNase T2 Family and Identification of Interacting Proteins of Four ClS-RNase Genes in ‘XiangShui’ Lemon

S-RNase plays vital roles in the process of self-incompatibility (SI) in Rutaceae plants. Data have shown that the rejection phenomenon during self-pollination is due to the degradation of pollen tube RNA by S-RNase. The cytoskeleton microfilaments of pollen tubes are destroyed, and other components cannot extend downwards from the stigma and, ultimately, cannot reach the ovary to complete fertilisation. In this study, four S-RNase gene sequences were identified from the ‘XiangShui’ lemon genome and ubiquitome. Sequence analysis revealed that the conserved RNase T2 domains within S-RNases in ‘XiangShui’ lemon are the same as those within other species. Expression pattern analysis revealed that S₃-RNase and S₄-RNase are specifically expressed in the pistils, and spatiotemporal expression analysis showed that the S₃-RNase expression levels in the stigmas, styles and ovaries were significantly higher after self-pollination than after cross-pollination. Subcellular localisation analysis showed that the S₁-RNase, S₂-RNase, S₃-RNase and S₄-RNase were found to be expressed in the nucleus according to laser confocal microscopy. In addition, yeast two-hybrid (Y2H) assays showed that S₃-RNase interacted with F-box, Bifunctional fucokinase/fucose pyrophosphorylase (FKGP), aspartic proteinase A1, RRP46, pectinesterase/pectinesterase inhibitor 51 (PME51), phospholipid:diacylglycerol acyltransferase 1 (PDAT1), gibberellin receptor GID1B, GDT1-like protein 4, putative invertase inhibitor, tRNA ligase, PAP15, PAE8, TIM14-2, PGIP1 and p24beta2. Moreover, S₃-RNase interacted with TOPP4. Therefore, S₃-RNase may play an important role in the SI of ‘XiangShui’ lemon.

Investigating the Mechanism of Unilateral Cross Incompatibility in Longan (Dimocarpus longan Lour.) Cultivars (Yiduo × Shixia)

Longan (Dimocarpus longan Lour.) is an important subtropical fruit tree in China. Nearly 90% of longan fruit imports from Thailand are from the cultivar Yiduo. However, we have observed that there exists a unilateral cross incompatibility (UCI) when Yiduo is used as a female parent and Shixia (a famous Chinese cultivar) as a male parent. Here, we performed a comparative transcriptome analysis coupled with microscopy of pistils from two reciprocal pollination combinations [Shixia♂ × Yiduo♀(SY) and Yiduo♀ × Shixia♂(YS)] 4, 8, 12, and 24 h after pollination. We also explored endogenous jasmonic acid (JA) and jasmonyl isoleucine (JA-Ile) levels in pistils of the crosses. The microscopic observations showed that the UCI was sporophytic. The endogenous JA and JA-Ile levels were higher in YS than in SY at the studied time points. We found 7,251 differentially expressed genes from the transcriptome analysis. Our results highlighted that genes associated with JA biosynthesis and signaling, pollen tube growth, cell wall modification, starch and sucrose biosynthesis, and protein processing in endoplasmic reticulum pathways were differentially regulated between SY and YS. We discussed transcriptomic changes in the above-mentioned pathways regarding the observed microscopic and/or endogenous hormone levels. This is the first report on the elaboration of transcriptomic changes in longan reciprocal pollination combination showing UCI. The results presented here will enable the longan breeding community to better understand the mechanisms of UCI.

Downregulated expression of S2-RNase attenuates self-incompatibility in "Guiyou No. 1" pummelo

Self-incompatibility (SI) substantially restricts the yield and quality of citrus. Therefore, breeding and analyzing self-compatible germplasm is of great theoretical and practical significance for citrus. Here, we focus on the mechanism of a self-compatibility mutation in 'Guiyou No. 1' pummelo (Citrus maxima), which is a spontaneous mutant of 'Shatian' pummelo (Citrus maxima, self-incompatibility). The rate of fruit set and the growth of pollen tubes in the pistil confirmed that a spontaneous mutation in the pistil is responsible for the self-compatibility of 'Guiyou No. 1'. Segregation ratios of the S genotype in F₁ progeny, expression analysis, and western blotting validated that the reduced levels of S₂-RNase mRNA contribute to the loss of SI in 'Guiyou No. 1'. Furthermore, we report a phased assembly of the 'Guiyou No. 1' pummelo genome and obtained two complete and well-annotated S haplotypes. Coupled with an analysis of SV variations, methylation levels, and gene expression, we identified a candidate gene (CgHB40), that may influence the regulation of the S₂-RNase promoter. Our data provide evidence that a mutation that affects the pistil led to the loss of SI in 'Guiyou No. 1' by influencing a poorly understood mechanism that affects transcriptional regulation. This work significantly advances our understanding of the genetic basis of the SI system in citrus and provides information on the regulation of S-RNase genes.

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.

Genomic insights into citrus domestication and its important agronomic traits

Citrus originated in Southeast Asia, and it has become one of the most important fruit crops worldwide. Citrus has a long and obscure domestication history due to its clonal propagation, long life cycle, wide sexual compatibility, and complex genetic background. As the genomic information of both wild and cultivated citrus becomes available, their domestication history and underlying traits or genes are becoming clear. This review outlines the genomic features of wild and cultivated species. We propose that the reduction of citric acid is a critical trait for citrus domestication. The genetic model representing the change during domestication may be associated with a regulatory complex known as WD-repeat-MYB-bHLH-WRKY (WMBW), which is involved in acidification and anthocyanin accumulation. The reduction in or loss of anthocyanins may be due to a hitchhiking effect of fruit acidity selection, in which mutation occurs in the common regulator of these two pathways in some domesticated types. Moreover, we have summarized the domestication traits and candidate genes for breeding purposes. This review represents a comprehensive summary of the genes controlling key traits of interest, such as acidity, metabolism, and disease resistance. It also sheds light on recent advances in early flowering from transgenic studies and provides a new perspective for fast breeding of citrus. Our review lays a foundation for future research on fruit acidity, flavor, and disease resistance in citrus.

TMT-Based Quantitative Proteomic Analysis Reveals the Crucial Biological Pathways Involved in Self-Incompatibility Responses in Camellia oleifera

Camellia oleifera is a valuable woody oil plant belonging to the Theaceae, Camellia oil extracted from the seed is an excellent edible oil source. Self-incompatibility (SI) in C. oleifera results in low fruit set, and our knowledge about the mechanism remains limited. In the present study, the Tandem mass tag (TMT) based quantitative proteomics was employed to analyze the dynamic change of proteins response to self- and cross-pollinated in C. oleifera. A total of 6,616 quantified proteins were detected, and differentially abundant proteins (DAPs) analysis identified a large number of proteins. Combined analysis of differentially expressed genes (DEGs) and DAPs of self- and cross-pollinated pistils based on transcriptome and proteome data revealed that several candidate genes or proteins involved in SI of C. oleifera, including polygalacturonase inhibitor, UDP-glycosyltransferase 92A1-like, beta-D-galactosidase, S-adenosylmethionine synthetase, xyloglucan endotransglucosylase/hydrolase, ABC transporter G family member 36-like, and flavonol synthase. Venn diagram analysis identified 11 proteins that may participate in pollen tube growth in C. oleifera. Our data also revealed that the abundance of proteins related to peroxisome was altered in responses to SI in C. oleifera. Moreover, the pathway of lipid metabolism-related, flavonoid biosynthesis and splicesome were reduced in self-pollinated pistils by the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. In summary, the results of the present study lay the foundation for learning the regulatory mechanism underlying SI responses as well as provides valuable protein resources for the construction of self-compatibility C. oleifera through genetic engineering in the future.

Evolution of self-compatibility by a mutant Sm-RNase in citrus

Self-incompatibility (SI) is an important mechanism that prevents self-fertilization and inbreeding in flowering plants. The most widespread SI system utilizes S ribonucleases (S-RNases) and S-locus F-boxes (SLFs) as S determinants. In citrus, SI is ancestral, and Citrus maxima (pummelo) is self-incompatible, while Citrus reticulata (mandarin) and its hybrids are self-compatible (SC). Here, we identify nine highly polymorphic pistil-specific, developmentally expressed S-RNases from pummelo that segregate with S haplotypes in a gametophytic manner and cluster with authentic S-RNases. We provide evidence that these S-RNases function as the female S determinants in citrus. Moreover, we show that each S-RNase is linked to approximately nine SLFs. In an analysis of 117 citrus SLF and SFL-like (SLFL) genes, we reveal that they cluster into 12 types and that the S-RNases and intra-haplotypic SLF and SLFL genes co-evolved. Our data support the notion that citrus have a S locus comprising a S-RNase and several SLFs that fit the non-self-recognition model. We identify a predominant single nucleotide mutation, S_m-RNase, in SC citrus, which provides a 'natural' loss of function. We show that SI-SC transitions due to the S_m-RNase initially arose in mandarin, spreading to its hybrids and became fixed. Identification of an evolutionarily distant new genus utilizing the S-RNase-based SI system, >100 million years separated from the nearest S-RNase family, is a milestone for evolutionary comparative studies.

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.

Comparative transcriptional survey between self-incompatibility and self-compatibility in Citrus reticulata Blanco

Seedlessness is an excellent economical trait, and self-incompatibility (SI) is one of important factors resulting in seedless fruit in Citrus. However, SI molecular mechanism in Citrus is still unclear. In this study, RNA-Seq technology was used to identify differentially expressed genes related to SI reaction of 'Wuzishatangju' (Citrus reticulata Blanco). A total of 35.67GB raw RNA-Seq data was generated and was de novo assembled into 50,364 unigenes with an average length of 897bp and N50 value of 1549. Twenty-three candidate unigenes related to SI were analyzed using qPCR at different tissues and stages after self- and cross-pollination. Seven pollen S genes (Unigene0050323, Unigene0001060, Unigene0004230, Unigene0004222, Unigene0012037, Unigene0048889 and Unigene0004272), three pistil S genes (Unigene0019191, Unigene0040115, Unigene0036542) and three genes (Unigene0038751, Unigene0031435 and Unigene0029897) associated with the pathway of ubiquitin-mediated proteolysis were identified. Unigene0031435, Unigene0038751 and Unigene0029897 are probably involved in SI reaction of 'Wuzishatangju' based on expression analyses. The present study provides a new insight into the molecular mechanism of SI in Citrus at the transcriptional level.

Genome-wide identification and functional analysis of S-RNase involved in the self-incompatibility of citrus

S-RNase-based self-incompatibility is found in Solanaceae, Rosaceae, and Scrophulariaceae, and is the most widespread mechanism that prevents self-fertilization in plants. Although 'Shatian' pummelo (Citrus grandis), a traditional cultivated variety, possesses the self-incompatible trait, the role of S-RNases in the self-incompatibility of 'Shatian' pummelo is poorly understood. To identify genes associated with self-incompatibility in citrus, we identified 16 genes encoding homologs of ribonucleases in the genomes of sweet orange (Citrus sinensis) and clementine mandarin (Citrus clementine). We preliminarily distinguished S-RNases from S-like RNases with a phylogenetic analysis that classified these homologs into three groups, which is consistent with the previous reports. Expression analysis provided evidence that CsRNS1 and CsRNS6 are S-like RNase genes. The expression level of CsRNS1 was increased during fruit development. The expression of CsRNS6 was increased during the formation of embryogenic callus. In contrast, we found that CsRNS3 possessed several common characteristics of the pistil determinant of self-incompatibility: it has an alkaline isoelectric point (pI), harbors only one intron, and is specifically expressed in style. We obtained a cDNA encoding CgRNS3 from 'Shatian' pummelo and found that it is homolog to CsRNS3 and that CgRNS3 exhibited the same expression pattern as CsRNS3. In an in vitro culture system, the CgRNS3 protein significantly inhibited the growth of self-pollen tubes from 'Shatian' pummelo, but after a heat treatment, this protein did not significantly inhibit the elongation of self- or non-self-pollen tubes. In conclusion, an S-RNase gene, CgRNS3, was obtained by searching the genomes of sweet orange and clementine for genes exhibiting sequence similarity to ribonucleases followed by expression analyses. Using this approach, we identified a protein that significantly inhibited the growth of self-pollen tubes, which is the defining property of an S-RNase.

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.

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.

Transcriptomic comparison of the self-pollinated and cross-pollinated flowers of Erigeron breviscapus to analyze candidate self-incompatibility-associated genes

BACKGROUND

Self-incompatibility (SI) is a widespread and important mating system that promotes outcrossing in plants. Erigeron breviscapus, a medicinal herb used widely in traditional Chinese medicine, is a self-incompatible species of Asteraceae. However, the genetic characteristics of SI responses in E. breviscapus remain largely unknown. To understand the possible mechanisms of E. breviscapus in response to SI, we performed a comparative transcriptomic analysis with capitulum of E. breviscapus after self- and cross-pollination, which may provide valuable information for analyzing the candidate SI-associated genes of E. breviscapus.

METHODS

Using a high-throughput next-generation sequencing (Illumina) approach, the transcriptionexpression profiling of the different genes of E. breviscapus were obtained, some results were verified by quantitative real time PCR (qRT-PCR).

RESULTS

After assembly, 63,485 gene models were obtained (average gene size 882 bp; N50 = 1485 bp), among which 38,540 unigenes (60.70% of total gene models) were annotated by comparisons with four public databases (Nr, Swiss-Prot, KEGG and COG): 38,338 unigenes (60.38% of total gene models) showed high homology with sequences in the Nr database. Differentially expressed genes were identified among the three cDNA libraries (non-, self- and cross-pollinated capitulum of E. breviscapus), and approximately 230 genes might be associated with SI responses. Several these genes were upregulated in self-pollinated capitulum but downregulated in cross-pollinated capitulum, such as SRLK (SRK-like) and its downstream signal factor, MLPK. qRT-PCR confirmed that the expression patterns of EbSRLK1 and EbSRLK3 genes were not closely related to SI of E. breviscapus.

CONCLUSIONS

This work represents the first large-scale analysis of gene expression in the self-pollinated and cross-pollinated flowers of E. breviscapus. A larger number of notable genes potentially involved in SI responses showed differential expression, including genes playing crucial roles in cell-cell communication, signal transduction and the pollination process. We thus hypothesized that those genes showing differential expression and encoding critical regulators of SI responses, such as MLPK, ARC1, CaM, Exo70A1, MAP, SF21 and Nod, might affect SI responses in E. breviscapus. Taken together, our study provides a pool of SI-related genes in E. breviscapus and offers a valuable resource for elucidating the mechanisms of SI in Asteraceae.

CrWSKP1, an SKP1-like Gene, Is Involved in the Self-Incompatibility Reaction of "Wuzishatangju" (Citrus reticulata Blanco)

Plant S-phase kinase-associated protein 1 (SKP1) genes play crucial roles in plant development and differentiation. However, the role of SKP1 in citrus is unclear. Herein, we described a novel SKP1-like gene, designated as CrWSKP1, from "Wuzishatangju" (Citrus reticulata Blanco). The cDNA sequence of CrWSKP1 is 779 base pairs (bp) and contains an open reading frame (ORF) of 477 bp. The genomic sequence of the CrWSKP1 gene is 1296 bp with two exons and one intron. CrWSKP1 has high identity with SKP1-like genes from other plant species within two conserved regions. Approximately 85% of pollen tubes of self-pollinated CrWSKP1 transgenic tobaccos became twisted at four days after self-pollination. Pollen tube numbers of self-pollinated CrWSKP1 transformants entering into ovules were significantly fewer than that of the control. Seed number of self-pollinated CrWSKP1 transformants was significantly reduced. These results suggested that the CrWSKP1 is involved in the self-incompatibility (SI) reaction of "Wuzishatangju".

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.