Comparative transcriptome analysis of stylar canal cells identifies novel candidate genes implicated in the self-incompatibility response of Citrus clementina
© Caruso et al; licensee BioMed Central Ltd. 2012
Received: 7 December 2011
Accepted: 14 February 2012
Published: 14 February 2012
Reproductive biology in citrus is still poorly understood. Although in recent years several efforts have been made to study pollen-pistil interaction and self-incompatibility, little information is available about the molecular mechanisms regulating these processes. Here we report the identification of candidate genes involved in pollen-pistil interaction and self-incompatibility in clementine (Citrus clementina Hort. ex Tan.). These genes have been identified comparing the transcriptomes of laser-microdissected stylar canal cells (SCC) isolated from two genotypes differing for self-incompatibility response ('Comune', a self-incompatible cultivar and 'Monreal', a self- compatible mutation of 'Comune').
The transcriptome profiling of SCC indicated that the differential regulation of few specific, mostly uncharacterized transcripts is associated with the breakdown of self-incompatibility in 'Monreal'. Among them, a novel F-box gene showed a drastic up-regulation both in laser microdissected stylar canal cells and in self-pollinated whole styles with stigmas of 'Comune' in concomitance with the arrest of pollen tube growth. Moreover, we identify a non-characterized gene family as closely associated to the self-incompatibility genetic program activated in 'Comune'. Three different aspartic-acid rich (Asp-rich) protein genes, located in tandem in the clementine genome, were over-represented in the transcriptome of 'Comune'. These genes are tightly linked to a DELLA gene, previously found to be up-regulated in the self-incompatible genotype during pollen-pistil interaction.
The highly specific transcriptome survey of the stylar canal cells identified novel genes which have not been previously associated with self-pollen rejection in citrus and in other plant species. Bioinformatic and transcriptional analyses suggested that the mutation leading to self-compatibility in 'Monreal' affected the expression of non-homologous genes located in a restricted genome region. Also, we hypothesize that the Asp-rich protein genes may act as Ca2+ "entrapping" proteins, potentially regulating Ca2+ homeostasis during self-pollen recognition.
Among the citrus species, several pummelos (Citrus grandis L. Osbeck) and mandarin-like varieties are self-incompatible . Clementine mandarin (Citrus clementina Hort. ex Tan.), derived from an uncontrolled cross between a sweet orange and a mandarin, is probably the most widespread citrus species showing self-incompatibility (SI). It is characterized by gametophytic SI, with the pollen tubes stopping their growth in the upper or middle style [2, 3]. Moreover, in this species SI is coupled with a variable degree of parthenocarpy. SI and parthenocarpy in citrus are very important traits for fruit production because they result in seedless fruits, which have a higher value in the markets compared to the seeded ones. Therefore, understanding the molecular basis of SI would be important to plan marker-assisted breeding to obtain new seedless genotypes.
Despite the importance of this trait, the genetic basis are still poorly understood and the key genes of SI have not been identified yet. The study of populations segregating for SI might be definitely a powerful strategy to give new insights into its genetic basis. However, the difficulties to obtain and characterize appropriate populations, with a variable degree of parthenocarpy and female/male sterility observed in the progeny might limit this approach. Such kind of strategy for the identification of the S-locus was carried out analyzing several crosses among different citrus cultivar and accessions with Got-3 isozyme, which is thought to be linked with the S-locus , providing only a rough estimation of their possible S-genotype.
In recent times, different research groups attempted to better understand SI and pollen-pistil interaction in several citrus genotypes, mainly trying to characterize putative homologs of key genes and proteins of already characterized SI systems. Gentile and colleagues  reported the involvement of Ca2+-dependent transglutaminase (TGase) in the self-incompatible response in pummelo, as already reported for Rosaceae . Regarding the S-locus genes, a S-like RNase has been isolated from 'Zigui shatian' pummelo , however the authors suggested that this gene might play an important role during ovary senescence rather than in the incompatibility response. Another S-like RNase has been isolated from a mandarin variety and was partially characterized , but it's not still clear whether this gene is the key determinant for the self-incompatible response.
To overcome these limits, the transcriptome analysis of natural mutants displaying contrasting compatibility behaviour might be more effective to better understanding the molecular basis controlling the progamic phase in citrus. Over the last decade, the genome and/or transcriptome analysis of natural or induced citrus mutants have been a powerful strategy to study the molecular basis of agronomically important traits, such as ripening period, fruit pigmentation, seedlessness and other traits related to quality [9–12]. Concerning SI, a few citrus natural mutants displaying different sexual behaviour with respect to their original varieties have been identified and characterized [13–16]. Differences between the mutants and the original cultivars were related to differences in pollen , style  or ovary  functionality. In some cases, different behaviours during the progamic phase were associated to abnormal embryo development [13, 16]. Therefore, it seems clear the mechanisms preventing fertilization are different in these genotypes, so it's reasonable to hypothesize that the mutations affected different genes or pathways implicated in reproduction.
Recently we chose two clementine clones with contrasting behaviour relating to self-pollen recognition ('Comune', self-incompatible; and 'Monreal', self-compatible mutation of 'Comune' ) as a model to identify candidate genes implicated in pollen-pistil interaction. Histological assays and analysis of pollen tube kinetics were performed to study the pollen tube behaviour in the two genotypes, and to assess whether the breakdown of SI in 'Monreal' was caused by changes in pistil or pollen functionality. The analysis demonstrated that the 'Monreal' mutation affected pistil functions, since pollen tubes of the two varieties grew equally in the pistils of self-compatible mutant, while 'Comune' rejected the pollen of both varieties, recognizing the pollen of 'Monreal' as self-pollen . A first transcriptome comparison was conducted analyzing whole styles with stigmas, leading to the identification of a first set of genes differentially expressed in non-pollinated flowers and during pollen-tube elongation in self-pollination condition  including stress related genes, and transcripts related to Ca2+ and hormone signalling. Surprisingly, a relatively high number of gene tags of different classes of retrotransposons were isolated, indicating their possible activation in response to pollination. However, the cDNA-AFLP analysis covered only part of the transcriptome and, due to the presence of pollen tubes growing along the pistil, the isolated gene tags were presumably not pistil-specific.
Here we describe a complementary approach to identify another set of candidate genes in 'Comune' and 'Monreal', to provide a better view of the molecular aspects related to citrus reproduction. We used laser capture microdissection (LCM), which has been efficiently used to investigate several aspects of plant reproductive biology [18–21]. In our approach, LCM was coupled to microarray analysis to identify genes specifically expressed in the stylar canal cells (SCC). We demonstrate that our transcriptomic survey of SCC is an efficient strategy to discover candidate genes involved in the SI response in clementine.
Laser microdissection of SCC and microarray analysis
Biotin-labeled aRNA was hybridized on the Affymetrix Citrus GeneChip® which contains 30,171 probe sets representing up to 33,879 citrus transcripts from several citrus species and hybrids. The array is estimated to represent about 15,500 genes .
List of genes showing differential expression between 'Comune' (self-incompatible) and 'Monreal' (self-compatible) stylar canal cells
Over-represented in the self incompatible (Comune)
No homology with functionally annotated proteins
No homology with functionally annotated proteins
No homology with functionally annotated proteins
No homology with functionally annotated proteins
Oligopeptide transporter-ISP4-like protein
Thioredoxin (TRX)-like [2Fe-2S] Ferredoxin (Fd) family
PSBS; photosystem II 22 kDa protein
Over-represented in the self compatible (Monreal)
No homology with functionally annotated proteins
40S Ribosomal protein S24
Glycosyltransferase, CAZy family
GASA gibberellin regulated cysteine-rich protein
similar to uncharacterized Arabidopsis mitochondrial gene ATMG00030
An F-box is the most up-regulated gene in the SCC of the self-incompatible genotype
cit.7568 was the unigene which showed the highest fold-change value, with more than 37-fold over-representation in 'Comune' (Table 1). This unigene is represented in the GeneChip by a single clone of a C. sinensis callus library. Attempts to get a deduced amino acid sequence of the cit.7568 unigene led to no reliable predictions. Actually, partial cDNA amplification and sequencing of cit.7568 revealed that the transcript expressed in the SCC of 'Comune' was 209 bp longer, sharing higher similarity to a single EST from C. sinensis early-developing fruits [GenBank: EY696233]. Sequencing of the genomic region surrounding the Open Reading Frame (ORF) revealed allelic variation, and translation evidenced a non-functional allele, having 100% homology with the cit.7568 locus of the clementine haploid genome, and a functional allele, which encodes an intronless 268 amino acids F-box protein [GenBank: JN885720]. A 1.5 kb region spanning the putative F-box locus was sequenced but no polymorphisms were detected between 'Comune' and 'Monreal', indicating that the differential transcriptional regulation was not related to nucleotide polymorphisms. BlastX search showed that the predicted functional ORF shared homology to a poplar putative F-box protein [Uniprot: B9H9G5]. The putative Arabidospis ortholog (At5g04010) was annotated as a non-specific F-box. Interproscan identified a single F-box domain Skp2-like (IPR022364) and no additional protein motif structures. The F-box domain was similar to the one of SLEEPY (SLY1) F-box of Arabidopsis (At4g24210), that interacts with DELLA proteins and is involved in GA signalling . However, no evident homology with SLY1 was observed outside the F-box domain.
Three aspartic acid-rich (Asp-rich) protein encoding genes are preferentially expressed in self-incompatible condition
Putative aspartic acid-rich (Asp-rich) protein genes found in different species, with information on the number of associated ESTs supporting their expression in different tissues, and total number of ESTs for each species
EST total no.
Callus, seedling, leaf (healthy and subjected to biotic and abiotic stress), mixed tissues
Root challenged with nematode or abiotic stresses, seedling, callus, ovaries (abscission zone), cambium
Roots with iron deficiency, seedling, CTV infected leaves
Bark challenged with insects
Flowers, early seeds, late seeds and stems
Mixed library of flower, leaf and root
Mixed, root, cotyledon, leaf, seedlings treated with salicylic acid, whole plant
Mixed, mixed stressed tissues, root, pod, flower, tissue culture
Pod, wood, leaf, flower, cushion
Leaves sprayed with the defense elicitors
Mixed tissues from water stressed plant
Stem, mixed (meristematic region, very young fiber, roots, stem)
M. × domestica
Young leaves inoculated with Marssonina coronaria, flower, fruit, bud, leaf
Abiotic and biotic stress-treated leaves, callus tissue and root tissue
Flowers and flower buds
Leaves induced by Ectropis obliqua feeding
Flower organ (pappus bristles)
Tissue culture growing 2 days in hypertonic medium
Flowers after meiosis
P. × hybrida
Roots, developing ovaries
Carrying out similarity searches among the sequenced plant genomes we identified Asp-rich homologs in Arabidopsis thaliana, Arabidopsis lyrata, Populus trichocarpa, Medicago truncatula, Glycine max, Ricinus communis, Theobroma cacao and Manihot esculenta. As in citrus, the Asp-rich protein genes are usually in tandem, with the number of linked genes depending on each species and varying from 2 to 7. Exceptions are represented by poplar and cassava, where single copies are present. Moreover, soybean contains two putative Asp-rich proteins which are placed in different chromosomes, while in Medicago, in addition to 5 clustered Asp-rich located in chromosome 2, there are two unclustered identical Asp-rich in chromosome 7 (Medtr7g104220 and Medtr7g104490) which are elements of a duplicated region. M. truncatula and A. lyrata possess the highest number of Asp-rich paralogs (seven and six respectively), while A. thaliana displays three clustered Asp-rich proteins; two are predicted in the TAIR genome browser, release 9 (At1g47395 and At1g47400), while the third one has no Gene ID but was predicted using Fgenesh and is supported by 37 ESTs. In some cases, as in A. thaliana, gene predictions are supported by ESTs, while in others, such as Ricinus and Medicago, there are no EST corresponding to the prediction. The low EST support found in this species is probably due to the relatively low number of ESTs deposited in the databases, but also imply the highly specific role of this genes in specific metabolic processes.
Another search was performed against Uniprot and dbEST databases to identify other Asp-rich homologs and gene tags from additional species (Table 2). In the Asp-rich homologs, aspartate is still the prevalent residue, with an average percentage of 22% in the amino acids composition of the predicted proteins. Also, as the citrus Asp-rich proteins they are characterized by a short length (between 44 and 94 residues), by a variable number of poly-D/G domains (from one to four), and the majority presents the polyD/G + YDYAPAA at the C-terminus (Figure 3). As in the case of citrus, most of the identified ESTs were generated from stressed tissues, and in some cases from flowers, pointing out to a possible conserved role in different species. Interestingly, one EST belongs to a cDNA library from the jelly fish Clytia hemisphaerica suggesting that the role of this genes might not be plant-specific.
The genomic region surrounding the Asp-rich protein genes is conserved among genomes
Further structural analysis of the genomic regions surrounding the Asp-rich protein genes evidenced another remarkable feature. Specifically, we identified a DELLA protein gene about 40 kb downstream the Asp-rich protein genes. DELLA was found to be differentially expressed in a previous transcriptome profiling . In particular, its mRNA levels were higher in 'Comune' pollinated styles with stigmas compared to 'Monreal' ones.
Histological and transcriptional changes during pollen-pistil interaction
In this study, LCM was used to spefically isolate SCC of two clementine genotypes differing for the SI response. 'Comune' is a widespread self-incompatible variety, while 'Monreal' is a self-compatible variety originated from a spontaneous bud mutation . The microdissection of the SCC allowed to perform a highly specific study of the transcriptome of the cells implicated in the interaction between pistil and pollen tubes, with the main aim of identifying candidate genes involved in self-pollen rejection. The results of microarray analysis suggested that the differential regulation of few specific transcripts might have lead to the breakdown of SI in 'Monreal'. Based on functional information retrieved from the databases, these genes do not show any clear functional enrichment, and many of them have no annotation, making it difficult to compare our results with SI systems characterized on model plants. On the other hand, our experiment provided a new set of transcripts that are very likely to play a key role during the progamic phase in citrus. Searches in the databases revealed that most of these genes are not flower specific, suggesting that probably the mutation leading to the breakdown of SI did not affected the S-locus determinants. However, it is known that non S-locus genes are implicated in SI . Most of the preferentially expressed genes were low-represented in EST databases and showed a weak expression in whole styles with stigmas before the SI reaction occurred. This indicates the usefulness of the LCM in the identification of highly specific and/or low expressed genes.
Flowers for LCM were sampled one day after pollination, when pollen was just germinating. Self-pollinations were performed to assess whether the presence of pollen in the stigma might have elicited a differential response of the SCC in self-compatible vs. self-incompatible conditions. So, while pollen tubes where still in the stigma, we targeted the rectangular SCC located in the upper/middle style, which is the site where incompatibility occurs in 'Comune', without contamination of pollen tubes. Although we did not analyze mRNAs of microdissected SCC of non-pollinated flowers, qRT-PCR demonstrated that slight differences in gene expression were also present in whole styles with stigmas before pollination, indicating that, at least in the cases of cit.7568, cit.11563, cit.5456 and cit.5776, the different mRNA levels in the SCC between 'Comune' and 'Monreal' were not induced by the pollen germination in the stigma.
Among the identified genes, our analysis focused on four unigenes over-represented in 'Comune' (cit.7568, cit.11563, cit.5456 and cit.5776). A time course analysis coupled to histological observation of pollen tube elongation was performed to correlate pollen tube behaviour with the difference in gene expression between 'Comune' and 'Monreal'. Until the 4th DAP we could not observe any clear differences in the mRNA level of the 4 analyzed genes, and no clear differences were evident in pollen tube behaviour between the two varieties. Form the 5th DAP we observed an impressive up-regulation in the self-incompatible combination. Selected candidate genes showed a clear differential expression (even > 100 fold change) during pollen - style interaction (Figure 7), in agreement with histological observation (pollen tube growth in 'Monreal' and no growth in 'Comune'). The peak of up-regulation was evidenced in concomitance with the pollen tube arrest, suggesting that all the analyzed genes have a key role in the stop of pollen tube elongation. The up-regulation was clearly induced by pollen tubes in the styles, as confirmed with the comparison of the expression levels of non-pollinated and self-pollinated flowers at T6 (Figure 7). The candidate genes were weakly expressed in pollen tubes and no differences in the mRNA levels were observed between 'Comune' and 'Monreal' (Figure 8). As a result, it is unlikely that the transcriptome of pollen tubes influenced the mRNA levels detected in the time course analysis. In fact, our expression data demonstrate that the drastic changes in mRNA levels occurred in the stylar tissues as a response to self-incompatible pollen tubes. Also, the drastic up-regulation of the four genes should not represent a downstream response to SI, since these genes are already differentially regulated in SCC before SI occurs.
Bioinformatic analysis provided useful information to hypothesize the role of the four unigenes. cit.7568 probe set matched to a putative F-box protein gene sharing homology to an Arabidopsis F-box (At5g04010) annotated as non-specific . This protein belongs to the C2 group of the F-box superfamily, of which SLY1 belongs. However, the similarity with SLY1 regards the F-box domain, while the C-terminus, responsible for targeting the substrate for ubiquitination, does not share any homology with SLY1.
F-box protein superfamily is one of the largest in plants [31, 32]. These proteins are key regulators of proteolysis, conferring specificity to the SCF (Skp1, Cullin, F-box) E3 ubiquitin ligase complex, which is responsible for recognizing the substrate for ubiquitin-mediated protein degradation. Despite the fact that it is not possible to correlate the role of the novel clementine F-box gene with any characterized orthologs, database searches and transcriptional data revealed that the clementine F-box gene was not previously identified in clementine cDNA libraries and evidenced highly specific expression patterns, since it was almost 250 fold up-regulated in SCC (Figure 2) and about 50 to 300 fold up-regulated in whole styles with stigmas (Figure 7A) in concomitance with the arrest of pollen tube elongation. This drastic up-regulation, and the absence of identical ESTs in databases support the high specificity of this gene and point out to its possible involvement in highly specific proteolytic events occurring in the style during the self-incompatible response. Studies on the families where SI system have been already characterized (namely Brassicaceae, Solanaceae and Papaveraceae) showed that proteolysis has a key role during self-pollen rejection. Usually proteolytic events follow the initial SI signal perception leading to the eventual death of the male gametophyte . Functional analysis will be necessary to shed light on the role of the clementine F-box during pollen-pistil interaction.
Another striking finding was the identification of the three putative Asp-rich protein encoding genes, cit.11563, cit.5456 and cit.5776, up regulated in 'Comune'. In addition to the high similarity in their primary structure, they co-localize in the clementine genome (Figure 5). The clustered genes displayed similar expression patters (Figure 7B, C, D) which might indicate that they are commonly regulated as already reported for other clusters [34, 35]. Although their function has not been investigated, it is known that the expression of the Arabidopsis homologous At1g47400 is sharply up regulated in plants exposed to 3-(30,40-dihydroxyphenyl)-L-alanine (L-DOPA), a phytotoxic allelochemical . In other studies, this gene turns out to be strongly induced by iron deficiency  and by treatment with the polycyclic aromatic hydrocarbon phenanthrene . These results support the hypothesis that the genes encoding the Asp-rich proteins might be triggered in response to different types of stresses. Since other information is still lacking, we attempt to assign a putative function to the Citrus Asp-rich proteins and also propose a model which overall might explain the regulation of SI in the 'Comune' variety. Because of their richness of aspartic acid residues, the Asp-rich proteins are supposed to act as novel Ca2+ "entrapping" proteins. Although they do not show sequence similarity with other well known Ca2+ interacting proteins, such as calsequestrin, calreticulin and calmodulin, the Asp-rich proteins share with those proteins the aspartate residue abundance which has been related with the protein ability of Ca2+-binding . Therefore, assuming that Ca2+ levels play a decisive role during pollen-pistil interaction, as already reported for several plant species , the up-regulation of the Asp-rich encoding genes observed in 'Comune' might enhance the amount of proteins functioning as Ca2+-trap elements, and lead to an exceptional decrease in Ca2+ availability, thus contributing to switch off the signal cascade usually induced by the increase in cytosolic Ca2+ concentration or to the alteration of Ca2+ gradient needed for pollen tube elongation. Pollen tube growth could represent, among others, the downstream physiological process dramatically affected by this missing triggering event. To support this hypothesis, it is worthwhile to mention that the aforesaid At1g47400 gene is up regulated in a T-DNA insertional mutant of a P-type ATPase cation pump, the MALE GAMETOGENESIS IMPAIRED ANTHERS (MIA); the mutant shows reduced male fertility and imbalanced cation homeostasis . Moreover, a Ca2+ pump (auto-inhibited Ca2+-ATPase - ACA), which presumably pumps Ca2+ out the cytosol , is differentially regulated in the two genotypes during pollen-pistil interaction showing an up regulation in the 'Comune' genotype . Further mandatory work will be undertaken to validate both the supposed role and the functioning model of pollen-pistil interaction in Citrus genotypes and, if they are proved, this might represent a case in which specific regulatory mechanism involving different loci rather than the S-locus could be co-responsible of the SI determination.
The integration of the transcriptomic data with the synteny analysis denoted a specific genome region containing a cluster of genes activated during self-pollen rejection (Figure 5). Specifically, the Asp-rich protein genes are linked to a DELLA gene which was previously isolated in the self-pollinated 'Comune' styles with stigmas  and showed a preferential expression in the self-incompatible genotype. This raise intriguing questions about the possibility that the non-homologous genes located in this genome region might contribute to a common function related to self-pollen rejection. Examples of co-expressed and functionally related gene clusters in eukaryotes have emerged over the last decade [34, 43] and are likely to increase with the huge amount of data coming from the genome projects. The most investigated operon-like organizations in plants are secondary metabolic pathways , mostly implicated in plant defence response. Clustering appears to have occurred de novo through some form of convergent evolutionary process . In our case, it is unlikely that the genes are clustered by chance, since this cluster is conserved in at least two other plant species, cacao and castor bean (Figure 5). The collinearity of the citrus genome segment comprising the Asp-rich protein genes and DELLA with other segments of unrelated genomes support the hypothesis that selection might have favoured the linkage of these genes. The advantage of clustering is related with the fact that tightly linked genes might be co-regulated at the levels of nuclear organization and/or chromatin [43, 45]. The co-localization of the three up-regulated Asp-rich protein genes as well as DELLA in the scaffold 9 of the clementine genome v0.9 suggests that the mutation leading to self compatibility probably affected the functionality of tightly linked genes. Further transcriptional analyses will be carried out on the other genes surrounding the Asp-rich protein genes to study their possible role during different stages of the pollen-pistil interaction.
LCM coupled to microarray analysis was definitely a powerful tool to identify candidate genes involved in self pollen rejection which have not been previously associated to SI. Data at the transcriptome level strongly suggested that a restricted number of differentially regulated transcripts are associated with self-pollen recognition. Although functional information is missing, we hypothesize that proteolysis and Ca2+ homeostasis might be crucial for SI response in clementine, reflecting to some extent molecular events occurring in other SI systems. However, further work at the protein level will be necessary for understanding the role of the candidate genes during self-pollen recognition. Also, the unigenes represented in the GeneChip do not cover the whole transcriptome, consequently other strategies such as RNA-seq might help to improve the transcriptome coverage to identify additional genes implicated in self-pollen recognition.
Plant material was collected from adult trees of 'Comune' clementine and its self-compatible natural mutation 'Monreal', grown at the experimental station of Consiglio per la Ricerca e la sperimentazione in Agricoltura (CRA) placed in Palazzelli, Lentini (SR), Italy.
Laser microdissection of SCC
Styles with stigmas, collected 24 hours after self-pollination, were immediately snap-frozen in OCT embedding medium (Sakura Finetek, Zoeterwoude, Netherlands) in Peel-A-Way plastic embedding molds (Polysciences, Polysciences, Warrington, PA, USA). The embedded samples were stored at -80°C until used. Transversal sections 10 μm thick at the upper part of the style were cut with a Leica CM1900 cryostat (Leica Microsystems, Germany) at -20°C. Cryosections were mounted on PET-membrane-coated stainless steel slides (Leica Microsystems, Wetzlar, Germany) and processed as described previously . Each slide contained 15-20 style sections. A Leica AS Laser Microdissection system (Leica Microsystems) was used for the isolation of stylar canals from transversal sections. Canals from the stigma were discarded to avoid contamination with pollen or pollen tubes. Microdissection was performed using the X40 magnification lens with the following settings: aperture: 6; intensity: 46; speed: 2; offset: 40; bridge: medium. SCC were collected in the cap of a 0.5 ml microtube filled with RLT buffer from the RNeasy Plus Micro Kit (Qiagen, Hilden, Germany).
RNA isolation and amplification
Three biological replicates were prepared for each genotype. Each biological replicate consisted of bulks of about 200 microdissected areas (composed of an average of 50 cells) coming from two different molds. RNA isolation was performed from ~ 10,000 cells using the RNeasy Micro Kit (Qiagen) following the manufacturer's instructions. RNA samples were subjected to a two-round amplification TargetAmp™ 2-Round Biotin-aRNA Amplification Kit 3.0 (Epicentre Biotechnologies, Madison, WI, USA) following the manufacturer's instructions. Quality of the aRNAs was assayed by OD260/OD280 measurements, agarose gel electrophoresis and Agilent 2100 Bioanalyzer.
Microarray experiments were performed using the GeneChip® Citrus Genome Array (Affymetrix, Santa Clara, CA, USA). Ten micrograms of biotynilated aRNA were used for the chip hybridizations. Data obtained from the microarray experiment were processed using Robin , which consists in an easy to use graphical interface for microarray analysis functions from R/BioConductor. CEL files were normalized by Robust Multi-Array Analysis (RMA) method. All microarray data have been deposited in the Gene Expression Omnibus (GEO) database under accession GSE33014. The list of probe sets showing significant differential expression was calculated comparing the three replicates for each genotype using Rank Products method . Percentage of false-positives (PFP) cut-off, corresponding to the false discovery rate, was set at 0.06 and calculated based on 100 permutations. Probe set information was retrieved from Harvest: Citrus database http://harvest.ucr.edu, http://www.harvest-web.org. Differentially expressed unigenes were subjected to Basic Local Alignment Search Tool (BLAST) search against the GenBank NR database. Interproscan http://www.ebi.ac.uk/Tools/InterProScan/ was used to identify conserved domains in the deduced protein sequences.
Validation of differentially expressed genes
Two RNA isolations (one for each genotype) were performed from additional bulks of SCC for the validation of microarray results. Each RNA was extracted from 10,000 cells and amplified using TargetAmp™ 2-Round aRNA Amplification Kit 2.0 to generate unlabeled aRNA. Quality of the aRNAs was assayed by OD260/OD280 measurements and agarose gel electrophoresis. cDNA was synthesized from 500 ng aRNA using the Ready-to-go RT-PCR beads (GE Healthcare Technologies, Little Chalfont, UK) with primer dT. Before performing qRT-PCR, an RT-PCR was carried out to check the presence of possible unspecific products. The transcript levels were determined by real-time quantitative RT-PCR (qRT-PCR) using the ABI Prism 7000 Sequence Detection System (Applied Biosystems, Foster City, CA, USA) and the Power SYBR Green PCR Master Mix as previously described . Real time was performed using three cDNA synthesis for each genotype, and each amplification was repeated twice. Primers used are listed in the Additional file 3. Primer pairs that did not reveal single dissociation peak were discarded. mRNA levels were calculated by standard-curve with a 5-fold dilution series and normalized against citrus Ubiquitin unigene contained in the GeneChip (AFFX-Cit-ubq11-3_x_at).
Ten self-pollinated flowers for each genotype were collected every 24 hours from the first to the eighth day after pollination to monitor pollen germination and pollen tube kinetics. Squash preparations and microscopy observation were performed as described in Distefano et al. .
Time course analysis
Selfed 'Comune' and 'Monreal' styles with stigmas were collected from 0 (virgin flowers) to 8 DAP. To check whether the expression patterns were influenced by pollen-pistil interaction, control samples of non-pollinated flowers were collected at 1, 3 and 6 DAP. Bulks of about 25 styles with stigmas were immediately stored in liquid nitrogen after collection and were used for total RNA isolation using RNeasy Plant Mini Kit according to the product manual (Qiagen). Total RNA was treated with RNase-Free DNase set (Qiagen) according to the manufacturer's instructions. Total RNA was amplified by PCR to confirm the absence of genomic DNA in the samples. cDNA synthesis was carried out starting from 2 μg total RNA. Real time runs and the data analysis were performed as described above. Two biological replicates were used for each sampling date, and all amplifications were repeated twice. Time course was carried out to evaluate expression patterns of 4 unigenes preferentially expressed in 'Comune' (probe set names cit.7568, corresponding to the SLY-like F-box gene; and cit.11563, cit. 5456, cit.5776, corresponding to the putative Asp-rich protein genes) as well as of a DELLA gene previously isolated by cDNA-AFLP from 'Comune' self-pollinated styles with stigmas [GenBank: GH733267].
In vitro pollen tube germination and RNA isolation
About 30 flowers were collected just before anthesis and their petals and pistils were removed. The anthers were left to dehisce for 24 h at room temperature. The fresh pollen was cultured in Petri dishes in the liquid pollen germination and tube growth medium described by Mesejo et al. . After 24 h of culture, pollen tubes were pelleted by centrifugation and harvested for RNA isolation using RNeasy Plant Mini Kit (Qiagen).
Similarity searches and phylogenetic analysis of the Asp-rich protein genes
DNA sequences of citrus unigenes corresponding to the Asp-rich protein were blasted against the sequenced genomes available at the Phytozome website v 7.0 http://www.phytozome.net. Other sequences of Asp-rich homologs were retrieved from GenBank, Uniprot and the Cacao Genome Database http://www.cacaogenomedb.org.
ORF prediction was based on the genome browser data for each analyzed species as well as Fgenesh analysis http://www.softberry.com. All the predicted protein sequences (3 citrus Asp-rich plus 44 from other species) were aligned by ClustalX . Alignments were analyzed using Mega 5 package  and a phylogenetic tree based on p-distance and pairwise deletions of gaps was constructed using the neighbor joining method with 1,000 bootstrap replicates.
Identification of the cit.7568 F-box functional allele
A 1,448 bp fragment that span the cit.7568 locus was amplified in both genotypes using the primers provided in the Additional file 3. The region was divided in three PCR products, which were cloned into pGem-T Easy Vector (Promega, Madison, WI, USA) and sequenced, or directly sequenced, using an ABI310 genetic analyzer (Applied Biosystems). The sequence of the region spanning the cit.7568 locus was deposited in the NCBI GenBank database [GenBank: JN885720].
Day(s) after pollination
Expressed sequence tag
Laser capture microdissection
Optimal cutting temperature
Open reading frame
Real time quantitative RT-PCR
Stylar canal cells
We thank Dr. Giuseppina Las Casas, Sergio Currò, Nicoletta Zingale for technical assistance, and Dr. Giuseppe Russo (CRA-ACM) who gave us the possibility to work at the Palazzelli experimental station. This work was supported by: the Italian Ministry of the University - Project PRIN "The productive process in fruit tree species: molecular, physiological and agronomical aspects of floral incompatibility and strategies for its control"; Regione Siciliana - Project "Risorse Genetiche Vegetali Sicilia". Work at the Centro de Genómica (IVIA) was supported by the Spanish Ministerio de Ciencia e Innovación-FEDER grant PSE-060000-2009-008 and the INIA grant RTA 2008-00065-00-00.
- Ollitrault P, Froelicher Y, Dambier D, Luro F, Yamamoto M: Seedlessness and ploidy manipulation. Citrus Genetics, Breeding and Biotechnology. Edited by: Khan IA. Wallingford: CABI; 2007:197-218.View ArticleGoogle Scholar
- Soost RK: Incompatibility alleles in the genus Citrus. Proc Amer Soc Hort Sci. 1965, 87: 176-180.Google Scholar
- Distefano G, Las Casas G, La Malfa S, Gentile A, Tribulato E, Herrero M: Pollen tube behavior in different mandarin hybrids. J Amer Soc Hort Sci. 2009, 134: 583-588.Google Scholar
- Ngo BX, Kim JH, Wakana A, Isshiki S, Mori T: Estimation of self-incompatibility genotypes of citrus cultivars with Got-3 allozyme markers. J Japan Soc Hort Sci. 2011, 80: 284-294. 10.2503/jjshs1.80.284.View ArticleGoogle Scholar
- Gentile A, Antognoni F, Iorio RA, Distefano G, Las Casas G, La Malfa S, Serafini-Fracassini D, Del Duca S: Polyamines and transglutaminase activity are involved in compatible and self-incompatible pollination of Citrus grandis. Amino Acids. 2012, 42: 1025-1035.PubMedView ArticleGoogle Scholar
- Di Sandro A, Del Duca S, Verderio E, Hargreaves A, Scarpellini A, Cai G, Cresti M, Faleri C, Iorio RA, Hirose S, Furutani Y, Coutts IGC, Griffin M, Bonner PLR, Serafini-Fracassini D: An extracellular transglutaminase is required for apple pollen tube growth. Biochem J. 2010, 429: 261-271. 10.1042/BJ20100291.PubMedView ArticleGoogle Scholar
- Chai LJ, Ge XX, Xu Q, Deng XX: CgSL2, an S-like RNase gene in 'Zigui shatian' pummelo (Citrus grandisOsbeck), is involved in ovary senescence. Mol Biol Rep. 2011, 38: 1-8. 10.1007/s11033-010-0070-x.PubMedView ArticleGoogle Scholar
- Miao HX, Qin YH, Teixeira-da-Silva JA, Ye ZX, Hu GB: Cloning and expression analysis of S-RNase homologous gene in Citrus reticulata Blanco cv. Wuzishatangju. Plant Sci. 2011, 180: 358-367. 10.1016/j.plantsci.2010.10.012.PubMedView ArticleGoogle Scholar
- Ríos G, Naranjo MA, Iglesias DJ, Ruiz-Rivero O, Geraud M, Usach A, Talón M: Characterization of hemizygous deletions in Citrus using array-comparative genomic hybridization and microsynteny comparisons with the poplar genome. BMC Genomics. 2008, 9: 381-10.1186/1471-2164-9-381.PubMedPubMed CentralView ArticleGoogle Scholar
- Alos E, Roca M, Iglesias DJ, Minguez-Mosquera MI, Damasceno CMB, Thannhauser TW, Rose JKC, Talon M, Cercos M: An evaluation of the basis and consequences of a stay-green mutation in the navel negra citrus mutant using transcriptomic and proteomic profiling and metabolite analysis. Plant Physiol. 2008, 147: 1300-1315. 10.1104/pp.108.119917.PubMedPubMed CentralView ArticleGoogle Scholar
- Ríos G, Naranjo MA, Rodrigo MJ, Alós E, Zacarías L, Cercós M, Talón M: Identification of a GCC transcription factor responding to fruit colour change events in citrus through the transcriptomic analyses of two mutants. BMC Plant Biol. 2010, 10: 276-10.1186/1471-2229-10-276.PubMedPubMed CentralView ArticleGoogle Scholar
- Talon M, Cercos M, Iglesias DJ, Colmenero-Flores JM, Ibáñez V, Brumos J, Herrero-Ortega MA, Rios G, Terol J, Tadeo FR: Citrus genomics and breeding: identification of candidate genes by the use of mutants and microarrays. Acta Hort. 2011, 892: 19-25.View ArticleGoogle Scholar
- Honsho C, Kotsubo M, Fukuda Y, Hamabata Y, Kurogi Y, Nishiwaki A, Tetsumura T: Reproductive Characteristics for Self-compatibility and Seedlessness in 'Nishiuchi Konatsu', a Bud Mutation of Hyuganatsu (Citrus tamurana hort. ex Tanaka). HortSci. 2009, 44: 1547-1551.Google Scholar
- Ye WJ, Qin YH, Ye ZX, Teixeira-da-Silva JA, Zhang LX, Wu XY, Lin SQ, Hu GB: Seedless mechanism of a new mandarin cultivar 'Wuzishatangju'(Citrus reticulata Blanco). Plant Sci. 2009, 177: 19-27. 10.1016/j.plantsci.2009.03.005.View ArticleGoogle Scholar
- Distefano G, Caruso M, La Malfa S, Gentile A, Tribulato E: Histological and molecular analysis of pollen-pistil interaction in clementine. Plant Cell Rep. 2009, 28: 1439-1451. 10.1007/s00299-009-0744-9.PubMedView ArticleGoogle Scholar
- Chai LJ, Ge XX, Biswas MK, Xu Q, Deng XX: Self-sterility in the mutant 'Zigui shatian' pummelo (Citrus grandis Osbeck) is due to abnormal post-zygotic embryo development and not self-incompatibility. Plant Cell Tiss Org. 2011, 104: 1-11. 10.1007/s11240-010-9793-6.View ArticleGoogle Scholar
- Hodgson RW: Horticultural varieties of citrus. Citrus industry. Edited by: Reuther W, Webber HJ, Batchelor LD. Berkeley: University of California Press;1967:1: 431-591.Google Scholar
- Suwabe K, Suzuki G, Takahashi H, Shiono K, Endo M, Yano K, Fujita M, Masuko H, Saito H, Fujioka T, Kaneko F, Kazama T, Mizuta Y, Kawagishi-Kobayashi M, Tsutsumi N, Kurata N, Nakazono M, Watanabe M: Separated transcriptomes of male gametophyte and tapetum in rice: validity of a laser microdissection (LM) microarray. Plant Cell Physiol. 2008, 49: 1407-1416. 10.1093/pcp/pcn124.PubMedPubMed CentralView ArticleGoogle Scholar
- Hobo T, Suwabe K, Aya K, Suzuki G, Yano K, Ishimizu T, Fujita M, Kikuchi S, Hamada K, Miyano M, Fujioka T, Kaneko F, Kazama T, Mizuta Y, Takahashi H, Shiono K, Nakazono M, Tsutsumi N, Nagamura Y, Kurata N, Watanabe M, Matsuoka M: Various spatiotemporal expression profiles of anther-expressed genes in rice. Plant Cell Physiol. 2008, 49: 1417-1428. 10.1093/pcp/pcn128.PubMedPubMed CentralView ArticleGoogle Scholar
- Tang X, Zhang ZY, Zhang WJ, Zhao XM, Li X, Zhang D, Liu QQ, Tang WH: Global gene profiling of laser-captured pollen mother cells indicates molecular pathways and gene subfamilies involved in rice meiosis. Plant Physiol. 2010, 154: 1855-1870. 10.1104/pp.110.161661.PubMedPubMed CentralView ArticleGoogle Scholar
- Wuest SE, Vijverberg K, Schmidt A, Weiss M, Gheyselinck J, Lohr M, Wellmer F, Rahnenfuhrer J, von Mering C, Grossniklaus U: Arabidopsis female gametophyte gene expression map reveals similarities between plant and animal gametes. Curr Biol. 2010, 20: 506-512. 10.1016/j.cub.2010.01.051.PubMedView ArticleGoogle Scholar
- Ciampolini F, Cresti M, Sarfatti G, Tiezzi A: Ultrastructure of the stylar canal cells of Citrus limon (Rutaceae). Plant Syst Evol. 1981, 138: 263-274. 10.1007/BF00985190.View ArticleGoogle Scholar
- Distefano G, Gentile A, Herrero M: Pollen-pistil interactions and early fruiting in parthenocarpic citrus. Ann Bot. 2011, 108: 499-509. 10.1093/aob/mcr187.PubMedPubMed CentralView ArticleGoogle Scholar
- Aprile A, Federici C, Close TJ, De Bellis L, Cattivelli L, Roose ML: Expression of the H+-ATPase AHA10 proton pump is associated with citric acid accumulation in lemon juice sac cells. Funct Integr Genomics. 2011, 11: 551-563.PubMedView ArticleGoogle Scholar
- Breitling R, Armengaud P, Amtmann A, Herzyk P: Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004, 573: 83-92. 10.1016/j.febslet.2004.07.055.PubMedView ArticleGoogle Scholar
- Casson S, Spencer M, Walker K, Lindsey K: Laser capture microdissection for the analysis of gene expression during embryogenesis of Arabidopsis. Plant J. 2005, 42: 111-123. 10.1111/j.1365-313X.2005.02355.x.PubMedView ArticleGoogle Scholar
- Hong F, Breitling R: A comparison of meta-analysis methods for detecting differentially expressed genes in microarray experiments. Bioinformatics. 2008, 24: 374-382. 10.1093/bioinformatics/btm620.PubMedView ArticleGoogle Scholar
- Dill A, Thomas SG, Hu J, Steber CM, Sun TP: The Arabidopsis F-box protein SLEEPY1 targets gibberellin signaling repressors for gibberellin-induced degradation. Plant Cell. 2004, 16: 1392-1405. 10.1105/tpc.020958.PubMedPubMed CentralView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2. 0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
- Kondo K, Yamamoto M, Itahashi R, Sato T, Egashira H, Hattori T, Kowyama Y: Insights into the evolution of selfcompatibility in Lycopersicon from a study of stylar factors. Plant J. 2002, 30: 143-153. 10.1046/j.1365-313X.2002.01275.x.PubMedView ArticleGoogle Scholar
- Gagne JM, Downes BP, Shiu SH, Durski AM, Vierstra RD: The F-box subunit of the SCF E3 complex is encoded by a diverse superfamily of genes in Arabidopsis. Proc Natl Acad Sci USA. 2002, 99: 11519-11524. 10.1073/pnas.162339999.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuroda H, Takahashi N, Shimada H, Seki M, Shinozaki K, Matsui M: Classification and expression analysis of Arabidopsis F-box containing protein genes. Plant Cell Physiol. 2002, 43: 1073-1085. 10.1093/pcp/pcf151.PubMedView ArticleGoogle Scholar
- Zhang Y, Zhao Z, Xue Y: Roles of proteolysis in plant self-incompatibility. Ann Rev Plant Biol. 2009, 60: 21-42. 10.1146/annurev.arplant.043008.092108.View ArticleGoogle Scholar
- Hurst LD, Pal C, Lercher MJ: The evolutionary dynamics of eukaryotic gene order. Nat Rev Genet. 2004, 5: 299-310.PubMedView ArticleGoogle Scholar
- Osbourn AE, Field B: Operons. Cell Mol Life Sci. 2009, 66: 3755-3775. 10.1007/s00018-009-0114-3.PubMedPubMed CentralView ArticleGoogle Scholar
- Golitz A, Sugano M, Hiradate S, Fujii Y: Microarray analysis of Arabidopsis plants in response to allelochemical L-DOPA. Planta. 2011, 233: 231-240. 10.1007/s00425-010-1294-7.View ArticleGoogle Scholar
- Buckhout TJ, Yang TJ, Schmidt W: Early iron-deficiency-induced transcriptional changes in Arabidopsis roots as revealed by microarray analyses. BMC Genomics. 2009, 10: 147-10.1186/1471-2164-10-147.PubMedPubMed CentralView ArticleGoogle Scholar
- Weisman D, Alkio M, Colon-Carmona A: Transcriptional responses to polycyclic aromatic hydrocarbon-induced stress in Arabidopsis thaliana reveal the involvement of hormone and defense signaling pathways. BMC Plant Biol. 2010, 10: 59-10.1186/1471-2229-10-59.PubMedPubMed CentralView ArticleGoogle Scholar
- Shin DW, Ma J, Kim DH: The asp-rich region at the carboxyl terminus of calsequestrin binds to Ca2+ and interacts with triadin. FEBS Lett. 2000, 486: 178-182. 10.1016/S0014-5793(00)02246-8.PubMedView ArticleGoogle Scholar
- Ge LL, Tian HQ, Russell SD: Calcium function and distribution during fertilization in angiosperms. Am J Bot. 2007, 94: 1046-1060. 10.3732/ajb.94.6.1046.PubMedView ArticleGoogle Scholar
- Jakobsen MK, Poulsen LR, Schulz A, Fleurat-Lessard P, Moller A, Husted S, Schiott M, Amtmann A, Palmgren MG: Pollen development and fertilization in Arabidopsis is dependent on the male gametogenesis impaired anthers gene encoding a type V P-type ATPase. Genes Dev. 2005, 19: 2757-2769. 10.1101/gad.357305.PubMedPubMed CentralView ArticleGoogle Scholar
- Schiott M, Romanowsky SM, Baekgaard L, Jakobsen MK, Palmgren MG, Harper JF: A plant plasma membrane Ca2+ pump is required for normal pollen tube growth and fertilization. Proc Natl Acad Sci USA. 2004, 101: 9502-9507. 10.1073/pnas.0401542101.PubMedPubMed CentralView ArticleGoogle Scholar
- Osbourn A: Gene clusters for secondary metabolic pathways: an emerging theme in plant biology. Plant Physiol. 2010, 154: 531-535. 10.1104/pp.110.161315.PubMedPubMed CentralView ArticleGoogle Scholar
- Chu HY, Wegel E, Osbourn A: From hormones to secondary metabolism: the emergence of metabolic gene clusters in plants. Plant J. 2011, 66: 66-79. 10.1111/j.1365-313X.2011.04503.x.PubMedView ArticleGoogle Scholar
- Amoutzias G, Van de Peer Y: Together we stand: genes cluster to coordinate regulation. Dev Cell. 2008, 14: 640-642. 10.1016/j.devcel.2008.04.006.PubMedView ArticleGoogle Scholar
- Agusti J, Merelo P, Cercos M, Tadeo FR, Talon M: Comparative transcriptional survey between laser-microdissected cells from laminar abscission zone and petiolar cortical tissue during ethylene-promoted abscission in citrus leaves. BMC Plant Biol. 2009, 9: 127-10.1186/1471-2229-9-127.PubMedPubMed CentralView ArticleGoogle Scholar
- Lohse M, Nunes-Nesi A, Krüger P, Nagel A, Hannemann J, Giorgi FM, Childs L, Osorio S, Walther D, Selbig J, Sreenivasulu N, Stitt M, Fernie AR, Usadel B: Robin: an intuitive wizard application for R-based expression microarray quality assessment and analysis. Plant Physiol. 2010, 153: 642-651. 10.1104/pp.109.152553.PubMedPubMed CentralView ArticleGoogle Scholar
- Mesejo C, Martinez-Fuentes A, Reig C, Rivas F, Agusti M: The inhibitory effect of CuSO4 on pollen germination and pollen tube growth and its application for the production of seedless fruit. Plant Sci. 2006, 170: 37-43.View ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar