- Research article
- Open Access
No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes
BMC Plant Biology volume 15, Article number: 129 (2015)
Fabaceae species are important in agronomy and livestock nourishment. They have a long breeding history, and most cultivars have lost self-incompatibility (SI), a genetic barrier to self-fertilization. Nevertheless, to improve legume crop breeding, crosses with wild SI relatives of the cultivated varieties are often performed. Therefore, it is fundamental to characterize Fabaceae SI system(s). We address the hypothesis of Fabaceae gametophytic (G)SI being RNase based, by recruiting the same S-RNase lineage gene of Rosaceae, Solanaceae or Plantaginaceae SI species.
We first identify SSK1 like genes (described only in species having RNase based GSI), in the Trifolium pratense, Medicago truncatula, Cicer arietinum, Glycine max, and Lupinus angustifolius genomes. Then, we characterize the S-lineage T2-RNase genes in these genomes. In T. pratense, M. truncatula, and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. In M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses also cluster with S-pollen genes. In T. pratense the S-RNase lineage genes show, however, expression in tissues not involved in GSI. Moreover, levels of diversity are lower than those observed for other S-RNase genes. The M. truncatula and C. arietinum S-RNase and S-pollen like genes phylogenetically related to Pyrinae S-genes, are also expressed in tissues other than those involved in GSI. To address if other T2-RNases could be determining Fabaceae GSI, here we obtained a style with stigma transcriptome of Cytisus striatus, a species that shows significant difference on the percentage of pollen growth in self and cross-pollinations. Expression and polymorphism analyses of the C. striatus S-RNase like genes revealed that none of these genes, is the S-pistil gene.
We find no evidence for Fabaceae GSI being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. There is no evidence that T2-RNase lineage genes could be determining GSI in C. striatus. Therefore, to characterize the Fabaceae S-pistil gene(s), expression analyses, levels of diversity, and segregation analyses in controlled crosses are needed for those genes showing high expression levels in the tissues where GSI occurs.
Useful agronomic traits can be found in wild populations of crop species. Nevertheless, a large fraction of species with hermaphroditic flowers have developed genetic mechanisms that allow the pistil to recognize and reject pollen from genetically related individuals (self-incompatibility; ), and this may affect the efficient incorporation of such traits into crop varieties. Self-incompatibility is, in general, evolutionarily advantageous, because it promotes cross-fertilization, and thus inbreeding depression avoidance.
Fabaceae is an economically important plant family with a large number of self-incompatible species (62.3% in Caesalpinioideae, 66.7% in Mimosoideae, and 22.1% in Papilionoideae sub families; ), that have been reported often as showing self-incompatibility of the gametophytic type (GSI; [1-9]). In GSI, if the specificity of the haploid pollen grain matches either one of the diploid pistil, an incompatible reaction occurs, leading to the degradation of the pollen tube within the pistil . It should be noted, however, that in all Fabaceae species where pollen tube growth was assessed in controlled crosses, only in species of the genus Trifolium the GSI reaction seems to be complete and takes place in the stlyle [3,11] as observed in Rosaceae (Rosidae; for a review see [12,13]), Solanaceae (Asteridae; ) and Plantaginaceae (Asteridae; [15,16]) SI species. In other species such as Vicia faba , Lotus corniculatus , Cytisus striatus , Coronilla emerus and Colutea arborescens  there is, however a significant difference on the percentage of pollen growth in self and cross-pollinations. In C. striatus, one of the species here studied, the percentage of ovules that are penetrated by pollen tubes is 72% in hand self-pollinated flowers compared with the 90.6% when hand cross-pollinations are performed . These authors have shown that an important fraction of self pollen grains collapse along the style, as observed in Rosaceae, Solanaceae and Plantaginaceae SI species.
Although the molecular characterization of the Fabaceae S-locus has never been performed, some authors have suggested that in Fabaceae GSI is RNase based [1,2,4-9]. Nevertheless, there are other GSI systems, such as that present in Papaveraceae [for a review see ]. Moreover, late-acting SI (LSI), so called because rejection of self-pollen takes place either in the ovary prior to fertilization, or in the first divisions of the zygote , has been described in Fabaceae [18,22-24]. It should be noted that, LSI can also be of the gametophytic type . In Fabaceae, however, the genetic basis of the different mechanisms that control LSI are mostly unknown, and thus, in this work we only address the possibility that Fabaceae GSI is determined by a S-RNase gene that clusters with those of the well characterized Rosaceae [12,13], Solanaceae  and Plantaginaceae [15,16] species. The most common ancestor of Fabaceae (Rosidae) and Rosaceae species lived about 89–91 million years ago (MYA; ). Since, according to phylogenetic analyses of the T2-RNases, RNase based GSI has evolved only once, before the split of the Asteridae and Rosidae, about 120 MYA [26-28], at least some Fabaceae SI species are expected to have this system. Therefore, in principle, a homology based approach could be used to identify the putative pistil S-gene in Fabaceae species.
Three amino acid patterns (amino acid patterns 1 and 2 that are exclusively found in proteins encoded by S-RNase lineage genes, and amino acid pattern 4 that is not found in any of the proteins encoded by S-RNase lineage genes), allow the distinction of S-RNase lineage genes from other T2 -RNase genes [28,29]. These patterns can be used to easily identify putative S-lineage genes using blast searches. The results can be further refined by selecting only those genes that encode basic proteins (isoelectric point higher than 7.5) since S-RNases have an isoelectric point between 8 and 10 . Furthermore, the number of introns can also be used to select S-lineage genes since S-RNases have one or two introns only (Figure one in ). Phylogenetic analyses where a set of reference genes are used, can then be performed to show that such genes belong, indeed, to the S-lineage. Nevertheless, in order to show that the identified genes are the pistil S-gene, it is necessary to show that they are highly expressed in pistils, although they can show lower expression in stigma and styles (see references in ). In Malus fusca where a large number of transcriptomes (flowers, pedicel, petal, stigma, style, ovary, stamen, filaments, anthers pollen, fruit, embryo and seed) have been analysed the same pattern is observed (CP Vieira, personal communication). Moreover, it is necessary to show that they have high polymorphism levels, that there is evidence for positive selection, and that in controlled crosses they co-segregate with S-locus alleles (see references in ).
The pollen component(s), always an F-box protein, has been identified as one gene in Prunus (Rosaceae; the gene is called SFB [32-37]), but multiple genes in Pyrinae (Rosaceae; the genes are called SFBBs [38-45]) and Solanaceae (called SLFs; [46-48]). F-box genes belong to a large gene family, and so far, no typical amino acid patterns have been reported for S-locus F-box protein sequences. Therefore, in non-characterized species, it is difficult to identify the pollen S-gene(s) using sequence data alone. In contrast to the S-RNase gene, Pyrinae SFBB genes show low polymorphism and high divergence [41-45]. Pollen S-gene(s) is (are), however, expected to be mainly expressed in the pollen [32,33,40,46,47].
Although the mechanism of self pollen tubes recognition is different when one or multiple S-pollen genes are involved [35,49], SSK1 (SKP1 like) proteins are involved in the self-incompatibility reaction in Rosaceae, Solanaceae and Plantaginaceae species, where GSI systems are well characterized. SKP1 like proteins are adapters that connect diverse F-box proteins to the SCF complex, and that are necessary in a wide range of cellular processes involving proteosome degradation (see references in ). SSK1 proteins have been described only in species having RNase based GSI [50-53], and thus, their presence has been suggested as a marker for RNase based GSI . These proteins are highly conserved and have a unique C-terminus, composed of a 5–9 amino acid residues following the conventional “WAFE” motif that is found in most plant SKP1 proteins . Therefore, the genes encoding such proteins can be easily retrieved using blast searches. In Solanaceae, Plantaginaceae, and Pyrinae, SSK1 proteins are expressed in pollen only [50-53], but in Prunus they are also expressed in styles .
To identify T2-RNases that could be S-locus candidate genes in Fabaceae subfamily Papilionoideae, in this work, we characterized the S-lineage T2-RNase genes in five genomes of species belonging to three major subclades: Trifolium pratense, Medicago truncatula, and Cicer arietinum from the inverted-repeat-lacking clade (IRLC), Glycine max from the millettioid clade, and Lupinus angustifolius from the genistoid clade. Trifolium and Medicago are the most closely related genera, and they share the most recent common ancestor, about 24 MYA . Cicer is diverging from these two genera for about 27 MY. Glycine is diverging from species of the IRLC clade for about 54 MY, and Lupinus is diverging from these for about 56 MY . Except for T. pratense, all these species are self-compatible. Nevertheless, the S-locus region could, in principle, be present, although the S-locus genes are expected to be non-functional . Compatible with this view, sequences closely related to the SSK1 genes are here identified in T. pratense, M. truncatula, C. arietinum, and G. max genomes. In T. pratense, M. truncatula and C. arietinum we identify S-RNase lineage genes that in phylogenetic analyses cluster with Pyrinae S-RNases. Furthermore, in M. truncatula and C. arietinum genomes, where large scaffolds are available, these sequences are surrounded by F-box genes that in phylogenetic analyses cluster with S-pollen genes. Nevertheless, none of these genes show expression only in tissues related with GSI. Moreover, T. pratense genes present levels of diversity lower than those of the characterized S-RNase genes. We also obtained a style with stigma transcriptome for Cytisus striatus, a species where self-pollen grains have been reported to collapse along the style, although partially . Once again, we found two genes that encode proteins showing the typical features of SSK1 genes and three T2-RNase like sequences, but none of these genes shows expression and variability levels compatible with being the S-RNase gene. Thus, we find no evidence for RNase based GSI in C. striatus. The data here presented supports the hypothesis that Fabaceae GSI is not determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. Alternative hypotheses are here discussed regarding the presence of SSK1 genes and Fabaceae GSI system.
SSK1 like genes in Fabaceae
SSK1 genes(s) are restricted to species having RNase based GSI [50-53]. The presence/absence of this gene(s) has been reported as a diagnosis marker for the presence/absence of RNase based GSI [50-53]. The protein encoded by SSK1 has an unique C-terminus, composed of 5–9 amino acid residues, following the conventional “WAFE” motif . In Rosaceae, this amino acid tail shows the conserved sequence “GVDED” (Additional file 5 in ). In Solanaceae and Plantaginaceae this motif is not so well conserved but a D residue is always found at the last position of the motif. It should be noted that most of the Fabaceae genomes that are available are from self-compatible species, and thus, SSK1 genes may be non-functional, or not involved in SI pathway. Therefore, when retrieving the sequences we allowed for some variability regarding these motifs (see Methods).
When using these features and the NCBI flowering plant species database, we retrieved 21 sequences from Solanaceae (three), Plantaginaceae (one), Rosaceae (eight), Fabaceae (five), Malvaceae (one), Rutaceae (one), Euphorbiaceae (one) and Salicaceae (one) species. Two other sequences, cy54873-cy21397 (this gene is the result of merging two sequences - cy54873g1 and cy21397g1 that overlap in a 22 bp region at the end of one and beginning of the other; PRJNA279853; http://evolution.ibmc.up.pt/node/77; http://dx.doi.org/10.5061/dryad.71rn0) and cy41479g1 (PRJNA279853; http://evolution.ibmc.up.pt/node/77; http://dx.doi.org/10.5061/dryad.71rn0) were identified in the C. striatus style with stigma transcriptome. These C. striatus sequences are incomplete at the 5′ region, since using blastx, the first 77 amino acids of SSK1 proteins are not present in these sequences. On the other hand, these sequences are complete at the 3′ region since their putative amino acid sequence presents the Rosaceae GVDED motif after the WAFE motif.
The phylogenetic relationship of the 23 SSK1 sequences, as well as the C-terminus sequence motif of the proteins they encode is presented in Figure 1 (see also Additional file 1). Fabaceae SSK1 like genes are more closely related to Rosaceae SSK1 sequences than to those from Solanaceae and Plantaginaceae (Figure 1), according to the known relationship of the plant families. It should be noted that only the two C. striatus deduced proteins present the Rosaceae GVDED motif after the WAFE motif. The T. pratense ASHM01022027.1, and G. max XM_003545885 genes encode proteins that present the WAFExxxxD motif, described for Solanaceae and Plantaginaceae SSK1. The presence of SSK1 genes in Fabaceae is, thus, consistent with the claims of RNase based GSI in Fabaceae.
SSK1 proteins showing the Rosaceae motif are also found in Hevea brasiliensis (Euphorbiaceae) and Populus trigonocarpa (Salicaceae). None of these species, or species of these families, has been described as having GSI. Furthermore, in Citrus clementina SSK1 like proteins present a proline instead of a glutamic acid in the Rosaceae WAFEGVDED motif. Citrus species present GSI and cytological analysis showed that growth of pollen tubes is arrested in different regions depending on the species analysed . In C. clementina pollen tubes are arrested in the upper styles . RNase activity has been identified in stigmas and pistils of C. reticulata [59,60] and also in ovaries of C. grandis , but the genetic mechanism is not clear yet . Indeed, in the comparative transcriptome analyses of stylar cells of a self-incompatible and a self-compatible cultivar of C. clementina, no T2-RNases where identified , rising doubts if GSI is RNase based in C. clementina. In T. cacao (Malvaceae) a SSK1 like protein with the same pattern as in C. clementina has also been identified. In this species self-pollen tubes grow to the ovary without inhibition, and self-incompatibility occurs at the embryo sac , and not in the style. Nevertheless, other Malvaceae species such as diploid species of the Tarasa genera present GSI (Table 1 in ), although the genetic mechanism is unknown.
T. pratense, M. truncatula, C. arietinum, G. max and L. angustifólio T2-RNase S-lineage genes
Given the evidence for the presence of RNase based GSI in Fabaceae (see above), we attempted to identify the S-RNase gene in Fabaceae species. Three main criteria were used to first identify putative S-RNase lineage genes in the T. pratense, M. truncatula, C. arietinum, G. max and L. angustifolius genomes, namely: 1) similarity at the amino acid level with S-RNases from Malus and/or Prunus (Methods); 2) the gene must encode a protein where amino acid pattern 4 is absent, once this pattern is found in proteins encoded by non-S-RNase lineage genes only [28,29]; and 3) the gene must encode a protein with an isoelectric point higher than 7.5, since S-RNases are always basic proteins [26,30]. Except for T. pratense, the genomes here analyzed are from self-compatible species. Nevertheless, the S-locus region could also be present, although the S-genes could show mutations that disrupt the coding region. For instance, in Rosaceae, mutated versions of the S-RNase and/or SFB genes have been described in self-compatible species . Table 1 summarizes the features of all gene sequences longer than 500 bp showing similarity at the amino acid level with S-RNases from Malus and/or Prunus. Although intron number was not used as a criterion for the selection of the genes, all these genes have one or two introns in the same location as those of the S-RNases . Three T. pratense (TP1, Tp5, and TP15, Table 1), two M. truncatula (Mt8 and Mt23, Table 1), five C. arietinum (Ca3, Ca6, Ca7, Ca12, Ca13, Table 1), and one G. max (Gm2, Table 1) genes are likely non-functional, since they present stop codons in their putative coding region. The number of putative S-lineage genes in T. pratense, M. truncatula, and C. arietinum (species from the IRLC clade) is about three times larger than in G. max (millettioid clade ) or L. angustifolius (from the genistoid clade). Although in C. arietinum the large number of T2-RNase lineage genes can be attributed to recent gene duplications, most of the T. pratense, and M. truncatula gene duplications are old (Figure 2, and Additional file 2). Three Lotus corniculatus, two L. japonicus, one Pisum sativum, one Cajanus cajan, one Lens culinaris, and one Cyamopsis tetragonoloba T2-RNase sequences that code for putative proteins without amino acid pattern 4, and that code for basic proteins were also included in the phylogenetic analyses (Additional file 3).
According to the phylogenetic analyses, the Fabaceae sequences that show amino acid patterns 1 and 2 (T. pratense Tp5, Tp8, Tp10, Tp11, Tp12, and Tp14, M. truncatula Mt12 and Mt13, C. arietinum Ca1, Ca3, Ca4, Ca10, Ca15, Ca17, and Ca18, L. corniculatus Lc3, and L. japonicus Lj4; Table 1 and Additional file 3), that are present in Rosaceae, Solanaceae, Plantaginaceae and Rubiaceae S-RNases [28,29], do not cluster toghether (Figure 2, and Additional file 2). Furthermore, Fabaceae genes - Tp6, Tp3, Ca4, Mt3, Mt17 and Mt18, in two of the alignment methods used (Figure 2, and Additional file 2B), cluster with Pyrinae S-RNases. Mt17 and Mt18 are neighbour genes (they are 3805 bp apart; Table 1). Mt17 is 56164 bp apart from Mt3 (Table 1). These genes could also represent the Fabaceae S-RNase. Although, the phylogenetic relationship of M. truncatula Mt20 gene and Plantaginaceae S-RNases depends on the alignment method used, we also included this gene in the following analyses.
Expression patterns of T. pratense Tp3, and Tp6, C. arietinum Ca4 and M. truncatula Mt3, Mt17, Mt18, and Mt20 genes
S-RNase expression is highest in pistils, although it can show lower expression in stigma and styles (CP Vieira, personal communication; see above; and [29-31,67]). For T. pratense we address the expression of genes Tp3, and Tp6 using cDNA of styles with stigmas, ovaries, and leaves. T3 gene shows expression in styles with stigmas, ovaries, and leaves (Figure 3A). For T6 gene, expression is observed in the styles with stigmas, and in leaves (Figure 3B). Since T. pratense is a SI species, these genes are thus, likely not S-RNases. Accordingly, levels of silent site (synonymous sites and non-coding positions) diversity for Tp3 and Tp6 genes are 0.008 and 0.011, respectively (based on five individuals and a genomic region of 447 bp and 414 bp, respectively). S-RNases show levels of silent variability higher than 0.23 .
Genes similar to the S-RNase but that are not involved in GSI may, in principle, show expression in other tissues. Indeed, S-RNase lineage 1 genes in Malus (Rosaceae) are expressed in embryo and seeds (Vieira CP, unpublished). This is in contrast to the S-RNase gene expression that is restricted to the stigma, styles and pistils of flowers at anthesis [29,30,67]. Therefore, genes showing expression in tissues other than the stigma, styles and pistils of flowers at anthesis are unlikely to be S-RNases. For C. arietinum Ca4 gene, blast searches against NCBI EST database shows that this gene is expressed in etiolated seedlings [GenBank:XM_004486248]). Thus, this gene is likely a gene not involved in GSI.
According to M. truncatula Gene Expression Atlas (Material and Methods) Mt20 ([GenBank:Mtr.49135.1.S1_at]) also shows expression in leaf and root tissues, among other tissues analysed. Since Mt3, Mt17 and Mt18 genes are not represented in the Affymetrix GeneChip, used in M. truncatula Gene Expression Atlas (Material and Methods), we addressed their expression using blastn and the SRA experiment sets for M. truncatula (99 RNA-Seq data sets from SRP033257 study from a mixed sample of M. truncatula root knot galls infected with Meloidogyne hapla (a nematode)). We find evidence for expression of the three genes in this large RNA-seq data set (Additional file 4). Therefore, according to gene expression, none of these genes seems to be determining pistil GSI specificity.
F-box genes in the vicinity of the C. arietinum Ca4 and M. truncatula Mt3, Mt17, Mt1, and Mt20 genes
At the S-locus region, the S-RNase gene is always surrounded by the S-pollen gene(s), that can be one gene as in Prunus (called SFB; [32-37], or multiple genes as in Pyrinae (called SFBBs; [38-41,45,47], and Solanaceae (called SLFs [14,46,47]). It should be noted that in Prunus, other F-box genes called SLFLs, not involved in GSI specificity determination  are also found surrounding the S-RNase gene [32,33]. Therefore, as an attempt to identify the S-locus in Fabaceae species, we identified all SFBBs/ SLFs, SLFLs, and SFB like genes in the vicinity (1 Mb) of the C. arietinum Ca4, and M. truncatula Mt3, Mt17, Mt18, and Mt20 genes (Figure 4, see Methods). For those gene sequences larger than 500 bp, phylogenetic inferences using reference genes (see Methods) show that C. arietinum Ca1_5 and M. truncatula Mt2_10, Mt2_11, and Mt7_7 are F-box genes that belong to the Malus, Solanaceae, and Plantaginaceae S-pollen and Prunus S- like pollen genes clade (Figure 5, and Additional file 5).
Expression pattern of the C. arietinum Ca1_5 and M. truncatula Mt2_10, Mt2_11, and Mt7_7 genes
Prunus SFB, Petunia and Antirrhinum SLFs, and Malus SFBB (S-pollen genes determining GSI specificity) genes have expression restricted to pollen and anthers [39-41,46,47,70]. Genes showing similarity to SLFs but that are not involved in GSI specificity determination (called SLFL) have also been described, but they have a broader pattern of expresion. For instance, in Prunus, SLFL genes are expressed in pollen and anthers but also in the style [32,33]. Furthermore, in Malus, SLFL genes are expressed in pollen, and anthers, but also in pistils, leaves, and seeds (Vieira CP, unpublished). Therefore, we addressed the expression pattern of C. arietinum Ca1_5 and M. truncatula Mt2_10, Mt2_11, and Mt7_7 genes.
C. arietinum Ca1_5 gene is expressed in etiolated seedlings ([GenBank:NW_004515210]), as the S-RNase like sequence located in its vicinity. Although we do not know if this gene is also expressed in pollen and anthers, because of its expression in seeds it is likely not involved in GSI. M. truncatula Mt7_7, and Mt2_11 genes, according to Gene Expression Atlas (Material and Methods), are expressed in leafs, petiole, stems, flowers, and roots, among other tissues analyzed (Mt7_7 -Mtr.14778.1.S1 _at, and Mt2_11 - Mtr.2939.1.S1_at). For Mt2_11 gene an EST ([GenBank:CA990259.1]) also supports expression of this gene in immature seeds 11 to 19 days after pollination. Mt2_10 gene is not represented in the Affymetrix GeneChip, and there is no EST data for this gene. Therefore, we addressed their expression using blastn and the SRA SRP033257 experiment data sets for M. truncatula (a mixed sample of M. truncatula root knot galls infected with M. hapla). We find evidence for expression of this gene in this large RNA-seq data set (Additional file 4). Therefore according to gene expression, none of these genes seems to be determining S-pollen GSI specificity.
T2-RNases from the C. striatus style with stigma transcriptome
Since we found no evidence in the available Fabaceae genomes for S-RNase like genes that could be involved in GSI specificity, we performed a transcriptome analysis of C. striatus styles with stigmas. This species has been described as having partial GSI . Five C. striatus sequences obtained from the style with stigma transcriptome show similarity with S-RNases (Table 2; PRJNA279853; http://evolution.ibmc.up.pt/node/77; http://dx.doi.org/10.5061/dryad.71rn0). CsRNase4, and CsRNase5 genes encode proteins with amino acid pattern 4, that is absent from all known S-RNases [28,29]. These genes encode putative acidic proteins (with an isoelectric point of 4.63 and 4.92, respectively), in contrast with S-RNases that are always basic proteins [26,30]. Furthermore, they share at least 85% amino acid similarity with other Fabaceae proteins that are expressed in tissues other than pistils (G. max [GenBank:XP_003518732.1], and [GenBank:XP_001235183.1], respectively). Moreover, these genes have three introns, and known S-RNases have only one or two introns . Therefore CsRNase4, and CsRNase5 genes are not S-RNases.
CsRNase1, and CsRNase2 genes code for proteins that do not present amino acid pattern 4, like the S-RNase gene (Table 2). Because the CsRNase3 coding sequence is incomplete, it is not possible to ascertain whether the protein encoded by this gene shows the amino acid pattern 4. Phylogenetic analyses of CsRNase1, and CsRNase2 genes, together with the sequences of other Fabaceae S-lineage genes, Rosaceae, Solanaceae, and Plantaginaceae S-RNases show that none of these genes belong to the known S-RNase gene lineages (Figure 6A, Additional file 6). CsRNase3 gene, however, clusters with Pyrinae S-RNases, and thus could represent a putative S-RNase gene (Figure 6B). For CsRNase3 gene, in the 266 bp region available, there are no introns. Accordingly, in the corresponding region there are no introns at the S-RNase gene. Nevertheless, unlike the S-RNases, CsRNase3 gene is expressed in ovaries, petals, leaves and fruits (Figure 7A). Moreover, levels of silent site diversity for this gene are moderate (π = 0.0233; based on a genomic region of 133 bp and five individuals of C. striatus from the Marecos population), but lower than that of the S-RNase gene (higher than 0.23; ). Thus, CsRNase3 gene does not present the expected features of a S-RNase gene.
Since we could not find any S-RNase candidate belonging to the Rosaceae, Solanaceae and Plantaginaceae S-RNase lineage genes, we characterized the CsRNase1 and CsRNase2 genes, that do not belong to any of the known S-RNases lineages. CsRNase1 gene is one of the most expressed genes (see Fragments Per Kilobase of target transcript length per Million reads mapped (FPKM) at position 24 in Additional file 7), but their genomic sequence revealed three introns (Additional file 8A). Moreover, CsRNase1 gene is expressed in ovaries, petals, pistils, leaves and fruits (Figure 7B), in contrast with the S-RNases that are expressed mainly in pistils [29,30,67]. Furthermore, levels of silent site (synonymous sites and non-coding positions) variability for this gene are low (π = 0.0006; based on a genomic region of 1020 bp and five individuals of C. striatus from the Marecos population) which is in sharp contrast with the expectation of high levels of variability at the S-RNase gene . Therefore, the overall evidence is that the CsRNase1 gene is not a S-pistil gene. For CsRNase2 gene the genomic sequence revealed five introns (Additional file 8B), it shows expression in ovaries, petals, pistils, leaves and fruits (Figure 7C), and low levels of silent site (synonymous sites and non-coding positions) variability (π = 0.0157; based on a genomic region of 1147 bp and five individuals of C. striatus from the Marecos population). Therefore, CsRNase2 gene is also not a S-pistil gene.
Phylogenetic analysis of T2-RNase genes from five Fabaceae genomes and one pistil transcriptome revealed more than six S-lineage genes. The two T. pratense genes that are phylogenetically related with Pyrinae S-RNases show, however, expression and polymorphism levels incompatible with being involved in GSI. Although the breeding system of the T. pratense individuals used in the polymorphism analyses was not characterized, in the literature all individuals analysed are SI [11,71,72]. Furthermore, red clover is described as being difficult to self, because of low seed set after selfing . Furthermore, the sequences obtained for the Portuguese population for the two T. pratense genes phylogenetically related with Pyrinae S-RNases, are very similar to those of the individual used for the T. pratense genome. Furthermore, none of the Fabaceae T2-RNase genes phylogenetically related with known S-RNases, revealed expression patterns compatible with a candidate Fabaceae S-pistil gene. It could be argued that only T. pratense is a self-incompatible species [71,73], and that the S-locus region may not be present in the other available genomes. Nevertheless, the presence of the same gene lineages in the T. pratense, M. truncatula and/or C. arietinum suggests that this is not the case. In Rosaceae, SC species still present the S-locus region, but S-RNase and SFB genes are non-functional . Nevertheless, mutations at loci involved in GSI but that are unlinked to the S-locus are also observed [74,75]. A similar pattern is also described in other SI systems such as that present in Brassicaceae family. For instance, the S-locus is present in the genome of the SC Arabidopsis thaliana, but the genes determining S-specificity are non-functional [76,77]. It should be noted, however, that the SI loss in M. truncatula is at least twice as old as that of A. thaliana. Therefore, genomes of SC species can also help in the identification of the putative S-locus genes.
The presence of Fabaceae sequences that cluster with Pyrinae S-RNases and S-pollen genes supports the hypothesis that we have identified the orthologous Pyrinae S-locus region in Fabaceae. These genes in Fabaceae seem to be performing functions other than GSI. Nevertheless, to exclude these genes as being the ones determining GSI, segregation analyses from controlled crosses are needed to show that these genes do not segregate as S-locus genes.
The fact that in Fabaceae, the Rosaceae, Solanaceae, and Plantaginaceae S-RNase gene lineages seem not to be involved in GSI, raises the hypothesis that in Trifolium GSI could be not RNase based. This hypothesis has been suggested before, based on the observation that on M. truncatula chromosome 1, that is largely syntenic to linkage group HG1 of T. pratense, where the S-locus has been mapped, there are no T2-RNases exhibiting significant similarity to Solanaceae, Rosaceae and Plantaginaceae S-RNases. The same observation has been reported for the numerous T2-RNase like sequences in the M. truncatula genome, even for those located near F-box genes, like the S-RNases . Nevertheless, under the current hypothesis, RNase based GSI evolved only once [26-28]. It is, however, conceivable that the ancestral S-locus has been duplicated during evolution. The presence of Fabaceae sequences presenting motifs 1 and 2 along the phylogeny support this hypothesis. In C. striatus, however, none of the T2-RNase genes expressed in pistils is determining GSI. Thus, there is no evidence to suggest that other T2-RNase lineage genes could be determining Fabaceae GSI. If this is the case, Fabaceae GSI has evolved the novo from T2-RNase unrelated genes, and thus, the information on Solanaceae, Rosaceae, Plantaginaceae and Rubiaceae S-RNases is not useful for the identification of the Fabaceae S-locus. It is expected that the S-pistil gene is highly expressed in the tissue where GSI occurs, and transcriptome analyses of this tissue can produce a list of genes showing high expression levels, such as that we present for C. striatus (see Additional file 7). Nevertheless, expression analyses, levels of diversity, and segregation analyses in controlled crosses will be needed to identify which gene(s) is(are) involved in S-pistil specificity.
It should be noted that in several Fabaceae species, we find SSK1 like genes with the typical features of those found in S-RNase based SI species from other plant families. It is conceivable that SSK1 like genes will be present in species where T2-RNase genes belonging to the S-lineage are present, even though such genes may not be involved in RNase based GSI. This must be the case for C. striatus. Moreover, the presence of a SSK1 like gene in C. clementina where no T2-RNases were identified from the transcriptome analyses of stylar cells of a self-incompatible and a self-compatible cultivar  offers support to this hypothesis.
The possibility that the frequency of self-incompatible species is overestimated in Fabaceae should not be also, ruled out. Indeed, the presence of binucleate pollen (typically associated with GSI), as well as fruit and seed production, are frequently used to assess the breeding system of a species. Nevertheless, other processes known to occur in Fabaceae species can affect fruit and seed production. For instance, Papilionoideae species have a membrane at the stigmatic surface that needs to be disrupted for pollen grain germination. In species of this subfamily the flower’s own pollen can cover the stigma at the bud stage [7,78-81], but it does not germinate while the stigmatic surface is intact [5,7,82]. With flowering maturation this stigmatic membrane in SI species must be scratched by a pollinator that visits the flower [7,19,83]. Moreover, late-acting self-incompatibility (LSI) has been described in many Fabaceae species such as Medicago sativa , Vicia faba , Pisum sativum , and Colutea arborescens  from the IRLC clade; Lotus corniculatus  and Coronilla emerus , both from the robinoid clade; Phaseolus vulgaris  from the millettioid clade; Dalbergia miscolobium  and D. retusa  from the dalbergioid clade; as well as in Genista hirsuta, Adenocarpus complicatus, Retama sphaerocarpa, Cytisus striatus, C. grandiflorus [7,83], and C. multiflorus [83,87] from the genistoid clade. In Fabaceae, LSI is due to multiple causes such as disharmony in endosperm/embryo development , differential growth rate of the pollen tubes within the ovaries [18,24], embryonic abortion [22,23] and inbreeding depression . Although the genetics and physiology of LSI is still poorly understood, it is clear that it can be genetically determined , and that LSI and GSI can co-occur, as it happens in C. striatus . Indeed, LSI implies similar growth of pollen tubes in the style following self- and cross-pollination (see for instance [88,89]), and in this species there is a significant difference in the percentage of pollen growth in self and cross-pollinations. Therefore, besides LSI, an additional partial GSI system has been inferred in C. striatus . Similar inferences have been made for V. faba , L. corniculatus , C. emerus and Colutea arborescens .
There is no evidence for Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes determining GSI in Fabaceae species. LSI is frequent in this family and may co-occur with GSI. Nevertheless, so far, in Fabaceae, only Trifolium species have been described as presenting GSI only. Thus, LSI or LSI in combination with GSI, will be likely the major hurdle when attempting to efficiently incorporating traits of agronomical interest from wild populations into crop varieties.
SSK1 like genes
To identify SSK1 like sequences in flowering plants we have used NCBI’s Pattern hit initiated blastp using as query A. hispanicum SSK1 ([GenBank:ABC84197.1]) and the pattern WAFExxxxD, as well as Pyrus x bretschneideri SSK1 like ([GenBank:CCH26218.1]), and Prunus avium SSK1 like ([GenBank:AFJ21661.1]) proteins and the pattern GVDED. For the non-annotated T. pratense genome ([GenBank:PRJNA200547]; ) we have obtained all putative open reading frames longer than 100 bp (getorf; http://emboss.sourceforge.net; ). Then we used local tblastn , with an Expect value of (e) < 0.05, and as query the above Rosaceae SSK1 like proteins.
T. pratense, M. truncatula, C. arietinum, G. max and L. angustifolius S-RNase lineage genes
Since four out of the five genomes here studied are from self-compatible species, S-pistil genes may be present as non-annotated pseudogenes. Therefore, putative open reading frames longer than 100 bp (getorf; http://emboss.sourceforge.net; ) were obtained for T. pratense ([GenBank:PRJNA200547]; ), M. truncatula ([GenBank:PRJNA30099], [GenBank:PRJNA10791], ; http://www.medicagohapmap.org), C. arietinum ([GenBank:PRJNA190909], [GenBank: PRJNA175619], ; http://cicar.comparative-legumes.org), G. max ([GenBank: PRJNA483899], [GenBank:PRJNA19861], ; http://www.Soybase.org) and L. angustifolius ([GenBank:PRJNA179231]; ; http://lupinus.comparative-legumes.org) genomes. Then, T2-RNase lineage sequences (including putative pseudogenes) of these species were identified and annotated by homology using local tblastn , using an Expect value of (e) < 0.05, and as query the M. domestica S2-RNase ([GenBank:AAA79841.1]), and P. persica S1-RNase ([GenBank:BAF42768.1]) proteins. If the inferred genes have been annotated before, the original name and accession number is indicated for that gene. Only sequences larger than 500 bp, and not presenting pattern 4 (absent in all S-RNases; ), were considered. In some cases, sequences were curated by introduction of sequence gaps to extend recognizable homology with the query sequence. Other Fabaceae T2-RNase sequences from M. sativa, Pisum sativum, Lens culinaris, (also belonging to IRLC), Lotus corniculatus, L. japonicus (from the robinoid clade), Cajanus cajan, (from the millettiod clade), Cyamopsis tetragonoloba (from the indigoferoid clade), and Arachis hypogaea (from the dalbergioid clade) were obtained from GenBank, using tblastn, an Expect value (e) < 0.05, and the above M. domestica, and P. persica sequences as query (Additional file 3). For all peptides, isoelectric points were calculated using ExPASy . Given the large number of genes analysed, for the sake of simplicity, in this work, we use short gene codes rather than the long mostly non-informative gene names. The correspondences between gene codes and gene names are given in Table 1, and Additional file 3.
F-box SFBB - and SFB - like genes in the vicinity of C. arietinum and M. truncatula S -RNase like genes
Putative open reading frames longer than 100 bp (getorf; http://emboss.sourceforge.net; ) were obtained for the 500 Kb of the C. arietinum and M. truncatula regions surrounding putative S-RNase lineage genes. F-box genes were identified and annotated by homology using local tblastn , an Expect value of (e) < 0.05), and the M. domestica SFBB3-beta ([GenBank:AB270796.1]), P. avium SFB3 ([GenBank:AY571665.1]), and P. axillaris S19-SLF ([GenBank:AY766154.1]) proteins. The correspondences between gene codes and gene names are given in Additional file 9.
Five data sets were used: 1- SSK1 like genes from flowering plants (that includes as reference sequences from Solanaceae, Plantaginaceae and Rosaceae SSK1 like genes), 2- Fabaceae S-RNase like genes that encode proteins with an isoelectric point higher than 7.5 (S-RNases are always basic proteins; ), with the exception of the Mt5, Mt26, Ca14 and Ca19 sequences that result in the introduction of many alignment gaps in the resulting alignment. Reference sequences are Solanaceae, Plantaginaceae and Rosaceae S-RNase genes, 3- C. arietinum and M. truncatula F-box SFBB - and SFB - like genes in the vicinity of S-RNase lineage genes. Reference sequences are Solanaceae and Plantaginaceae SLFs, Malus SFBBs and Prunus SFB, and Rosaceae S-pollen like genes (genes similar to S-pollen genes but that are not involved in GSI specificity), 4- C. striatus CsRNase1, and CsRNase2 genes. Reference sequences are Fabaceae S-RNase like genes that encode proteins with an isoelectric point higher than 7.5, Solanaceae, Plantaginaceae and Rosaceae S-RNase genes, and 5- C. striatus CsRNase3 gene. Reference sequences are Solanaceae, Plantaginaceae and Rosaceae S-RNase genes. With the exception of data set 5 (because of the size (264 bp) of C. striatus CsRNase3 sequence), sequences in the data sets were aligned with the ClustalW2, Muscle and T-coffee alignment algorithms as implemented in ADOPS . Only codons with a support value above two are used for phylogenetic reconstruction. Bayesian trees were obtained using MrBayes 3.1.2 , as implemented in the ADOPS pipeline, using the GTR model of sequence evolution, allowing for among-site rate variation and a proportion of invariable sites. Third codon positions were allowed to have a gamma distribution shape parameter different from that of first and second codon positions. Two independent runs of 2,000,000 generations with four chains each (one cold and three heated chains) were set up. The average standard deviation of split frequencies was always about 0.01 and the potential scale reduction factor for every parameter about 1.00 showing that convergence has been achieved. Trees were sampled every 100th generation and the first 5000 samples were discarded (burn-in). The remaining trees were used to compute the Bayesian posterior probabilities of each clade of the consensus tree.
In the phylogenetic analyses that include C. striatus CsRNase3 gene we used the MEGA 5 software . The alignment was performed using ClustalW, and for the phylogenetic reconstruction we used pairwise deletion and minimum evolution method. We run 10000 bootstrap replications, using maximum composite likelihood method, and including transitions + transversions substitutions, and all codons.
Expression of T. pratense Tp3 and Tp6 genes in styles with stigmas, ovaries, petals and leaves
To collect enough material for the cDNA synthesis of style with stigma (since in T. pratense each individual has less than three inflorescences with less than 50 flowers at anthesis), we have mixed the plant material obtained from two different individuals. These individuals present an amplification product of the expected size, obtained from genomic DNA (extracted from leaves, using the method of Ingram et al. ), using specific primers for Tp3 and Tp6 genes (Additional file 10), and standard amplification conditions of 35 cycles of denaturation at 94°C for 30 s, primer annealing temperature according to Additional file 10 for 30 s, and primer extension at 72°C for 2 min. More than 500 styles with stigmas were collected from these two individuals, that were frozen in liquid nitrogen and stored at −80°. For one of these individuals we also collected ovaries, and leaves. Total RNA was extracted using TRIzol® (Invitrogen, Spain) according to the manufacturer’s instructions and treated with DNase I (Turbo RNase-Free) (Ambion, Portugal). RNA quantity was assessed by NanoDrop v.1.0 (Thermo Scientific). cDNA was synthesized with SuperScript® III First-Strand Synthesis System for RT-PCR from Invitrogen. Elongation factor 1-α (Elf1-α) was used as positive control for cDNA synthesis. Standard amplification conditions as described above were used.
Levels of diversity at T. pratense Tp3 and Tp6 genes
To determine levels of diversity for Tp3 and Tp6 genes, genomic DNA from leaves of five T. pratense individuals of a Porto population (assigned as TpPorto1to TpPorto5) was extracted using the method of Ingram et al. . For each individual, genomic DNA was used in PCR reactions using primers 1821 F + 1821R, and 5450 F + 5450R, to amplify Tp3 and Tp6 genes, respectively (Additional file 10). Standard amplification conditions were 35 cycles of denaturation at 94°C for 30 seconds, primer annealing according to supplementary Additional file 10 for 30 s, and primer extension at 72°C for 3 min. The amplification products were cloned, using the TA cloning kit (Invitrogen, Carlsbad, CA). For each amplification product, the insert of an average of 10 colonies was cut separately with RsaI, and Sau3AI restriction enzymes. For each restriction pattern three colonies were sequenced in order to obtain a consensus sequence. The ABI PRISM BigDye cycle-sequencing kit (Perkin Elmer, Foster City, CA), and specific primers, or the primers for the M13 forward and reverse priming sites of the pCR2.1 vector, were used to prepare the sequencing reactions. Sequencing runs were performed by STABVIDA (Lisboa, Portugal). DNA sequences were deposited in GenBank (accession numbers KR054719 - KR054728). Nucleotide sequences were aligned using ClustalW algorithm as implemented in MEGA 5 . Analyses of DNA polymorphism were performed using DnaSP (version 4.1) .
Expression of M. truncatula Mt3, Mt17, Mt18, Mt20, Mt7_7, Mt2_10, and Mt2_11 genes
For the genes of interest, using blast at M. truncatula gene expression atlas (http://mtgea.noble.org/v3/; Affymetrix GeneChip Medicago Genome Array; ) we identify Probeset ID and the expression pattern associated with that probe. For the genes not represented in the M. truncatula gene expression atlas, we used blastn and the SRA SRP033257 experiment sets for M. truncatula (99 RNA-Seq data sets from a mixed sample of M. truncatula root knot galls infected with Meloidogyne hapla (a plant-nematode)).
Cytisus striatus style with stigma transcriptome
C. striatus has been described as having partial GSI, since a fraction (about 27%) of self-pollen tubes after hand self-pollination, stop growing along the style and the ovary . For one C. striatus individual (assigned as Cs1), from a population at Marecos (Valongo, Portugal), 400 flower buds ranging from 1.8 to 2 cm (the size of pre-anthesis stages;  were dissected to collect the styles with stigmas, that were frozen in liquid nitrogen and stored at −80°. Total RNA was extracted as described above. RNA quantity was assessed by NanoDrop v.1.0 (Thermo Scientific) and RNA quality by BioRad’s Experion System. A total RNA sample of approximately 2.691 μg ,with RQI of 7.1, and a 260/280 nm absorption ratio 2.08 was obtained. Total RNA was processed for Illumina RNA-Seq, at BGI (Hong Kong, China).
Only high quality reads were provided by BGI. Before assembly, adaptor sequences were removed from raw reads. FASTQC reports were then generated and based on this information the resulting reads were trimmed at both ends. Nucleotide positions with a score lower than 20 were masked (replaced by an N). These analyses were performed using the FASTQ tools implemented in the Galaxy platform [104-106]. The resulting high-quality reads were used in the subsequent transcriptome assembly using Trinity with default parameters . The Transcriptome project has been deposited at GenBank PRJNA279853, and the assembled transcriptome at http://evolution.ibmc.up.pt/node/77, or http://dx.doi.org/10.5061/dryad.71rn0. All contigs were used as queries for tblastn searches using local blast , and the SSK1 and S-RNase query sequences reported above. Fragments Per Kilobase of target transcript length per Million reads mapped (FPKM) values were estimated using the eXpress software  as implemented in Trinity. BLAST2Go  was used to determine PFAM (protein families) codes for the 100 most expressed genes.
The genomic sequence of the C. striatus S-lineage T2-RNases
To determine intron number for C. striatus CsRNase1, CsRNase2, and CsRNase3, primers were designed (Additional file 10) based on the sequences obtained from the transcriptome. Genomic DNA was extracted from leaves of the Cs1 individual, as described above, and used as template in PCR reactions. Standard amplification conditions were 35 cycles of denaturation at 94°C for 30 seconds, primer annealing according to Additional file 10 for 30 s, and primer extension at 72°C for 3 min. The amplification products were cloned, and sequenced as described above. The genomic sequences for C. striatus CsRNase1 and CsRNase2 genes of individual Cs1 were deposited at GenBank (accession numbers KR054703, and KR054709).
Expression of the C. striatus S-lineage T2-RNase genes in pollen, ovaries, petals, pistils, leaves and fruits
Pollen, ovaries, petals, pistils, leaves and fruits from individual Cs1 were collected and immediately frozen in liquid nitrogen and stored at −80°. Total RNA and cDNA synthesis was performed as described above. Elongation factor 1-α (Elf1-α) was used as positive control for cDNA synthesis. Primers CytSRN-62 F + CytisusRNase531R, CytR2-cons142F + CytR2-445R, and Cy10F + Cy10R were used for the amplification of the CsRNase1, CsRNase2, and CsRNase3 genes, respectively (Additional file 10). Standard amplification conditions were 35 cycles of denaturation at 94°C for 30 s, primer annealing temperature according to Additional file 10 for 30 s, and primer extension at 72°C for 2 min.
Nucleotide diversity at C. striatus S-lineage genes
To determine levels of diversity for CsRNase1, CsRNase2, and CsRNase3 genes, genomic DNA from leaves of four C. striatus individuals of the Marecos population (assigned as Cs2 to Cs5) was extracted as described above. For each individual, genomic DNA was used in PCR reactions using the same primers and conditions described above. The amplification products were cloned, as described above. For each amplification product, the insert of an average of 10 colonies was cut separately with RsaI, and Sau3AI restriction enzymes. For each restriction pattern three colonies were sequenced in order to obtain a consensus sequence. Sequencing has been performed as described above. DNA sequences were deposited in GenBank (accession numbers KR054704 - KR054707, KR054710 - KR054713, and KR054714 - KR054718, respectively). Nucleotide sequences were aligned using ClustalW algorithm as implemented in MEGA 5 . Analyses of DNA polymorphism were performed using DnaSP (version 4.1) .
Availability of supporting data
The C. striatus assembled transcriptome, supporting the results of this article is available in the [Cytisus striatus style with stigma transcriptome] repository [http://evolution.ibmc.up.pt/node/77], and at Dryad [http://dx.doi.org/10.5061/dryad.71rn0].
The data used to perform the phylogenetic analyses is available at Dryad [http://dx.doi.org/10.5061/dryad.71rn0].
Million years ago
Fragments per kilobase of target transcript length per million reads mapped
Isoelectric point (IP)
De Nettancourt D. Incompatibility in angiosperms. Berlin: Springer-Verlag; 1977.
Arroyo MTK, Armesto JJ, Villagran C. Plant phenological patterns in the high andean cordillera of central chile. J Ecol. 1981;69(1):205–23.
Atwood SS. Genetics of cross-incompatibility among self-incompatible plants of Trifolium Repens. J Am Soc Agron. 1940;32(12):955–68.
Heslop-Harrison J, Heslop-Harrison Y. Pollen-stigma interaction in the leguminosae: constituents of the stylar fluid and stigma secretion of Trifolium pratense L. Ann Bot. 1982;49(6):729–35.
Shivanna K, Owens S: Pollen-pistil interactions (Papilionoideae). Advances in Legume Biology: Monograph of Systematic Botany Missouri: Missouri Botanical Garden 1989;(29):157–182.
Weller S, Donoghue M, Charlesworth D. The evolution of self-incompatibility in flowering plants: a phylogenetic approach. In: Experimental and molecular approaches to plant biosystematics. 53rd ed. St Louis, Mo: Missouri Botanical Garden ((Monographs in Systematic Botany; 1995. p. 355–82.
Rodríguez-Riaño T, Ortega-Olivencia A, Devesa JA. Reproductive biology in two Genisteae (Papilionoideae) endemic of the western Mediterranean region: Cytisus striatus and Retama sphaerocarpa. Can J Bot. 1999;77(6):809–20.
Igic B, Lande R, Kohn JR. Loss of self‐incompatibility and its evolutionary consequences. Int J Plant Sci. 2008;169(1):93–104.
Casey NM, Milbourne D, Barth S, Febrer M, Jenkins G, Abberton MT, et al. The genetic location of the self-incompatibility locus in white clover (Trifolium repens L.). Theor Appl Genet. 2010;121(3):567–76.
Kao TH, Tsukamoto T. The molecular and genetic bases of S-RNase-based self-incompatibility. Plant Cell. 2004;16:S72–83.
Leduc N, Douglas G, Monnier M, Connolly V. Pollination in vitro: effects on the growth of pollen tubes, seed set and gametophytic self-incompatibility in Trifolium pratense L. and T. repens L. Theor Appl Genet. 1990;80(5):657–64.
Tao R, Iezzoni AF. The S-RNase-based gametophytic self-incompatibility system in Prunus exhibits distinct genetic and molecular features. Sci Hortic (Amsterdam). 2010;124(4):423–33.
De Franceschi P, Dondini L, Sanzol J. Molecular bases and evolutionary dynamics of self-incompatibility in the Pyrinae (Rosaceae). J Exp Bot. 2012;63(11):4015–32.
Wang N, Kao TH. Self-incompatibility in Petunia: a self/nonself-recognition mechanism employing S-locus F-box proteins and S-RNase to prevent inbreeding. Wiley Interdiscip Rev Dev Biol. 2011;1(2):267–75.
Xue Y, Carpenter R, Dickinson HG, Coen ES. Origin of allelic diversity in Antirrhinum S locus RNases. Plant Cell. 1996;8(5):805–14.
Vieira C, Charlesworth D. Molecular variation at the self-incompatibility locus in natural populations of the genera Antirrhinum and Misopates. Heredity. 2002;88(3):172–81.
Rowlands D. The nature of the breeding system in the field bean (V. faba L) and its relationship to breeding for yield. Heredity. 1958;12:113–26.
Bubar JS. Differences between self-incompatibility and self-sterility. Nature. 1959;183(4658):411–2.
Galloni M, Podda L, Vivarelli D, Cristofolini G. Pollen presentation, pollen-ovule ratios, and other reproductive traits in Mediterranean Legumes (Fam. Fabaceae-Subfam. Faboideae). Plant Syst Evol. 2007;266(3–4):147–64.
Wilkins KA, Poulter NS, Franklin-Tong VE. Taking one for the team: self-recognition and cell suicide in pollen. J Exp Bot. 2014;65(5):1331–42.
Allen AM, Hiscock SJ. Evolution and phylogeny of self-incompatibility systems in Angiosperms. In: Self-incompatibility in flowering plants. Berlin Heidelberg: Springer; 2008. p. 73–101.
Briggs C, Westoby M, Selkirk P, Oldfield R. Embryology of early abortion due to limited maternal resources in Pisum sativum L. Ann Bot. 1987;59(6):611–9.
Sage TL, Webster BD. Seed abortion in Phaseolus vulgaris L. Bot Gaz. 1990;151:167–75.
Brink R, Cooper D. Partial self-incompatibility in Medicago sativa. Proc Natl Acad Sci U S A. 1938;24(11):497–9.
Wikstrom N, Savolainen V, Chase MW. Evolution of the angiosperms: calibrating the family tree. Proc R Soc B - Biol Sci. 2001;268(1482):2211–20.
Igic B, Kohn JR. Evolutionary relationships among self-incompatibility RNases. Proc Natl Acad Sci U S A. 2001;98(23):13167–71.
Steinbachs JE, Holsinger KE. S-RNase-mediated gametophytic self-incompatibility is ancestral in eudicots. Mol Biol Evol. 2002;19(6):825–9.
Vieira J, Fonseca NA, Vieira CP. An S-RNase-based gametophytic self-incompatibility system evolved only once in eudicots. J Mol Evol. 2008;67(2):179–90.
Nowak MD, Davis AP, Anthony F, Yoder AD. Expression and trans-specific polymorphism of self-incompatibility RNases in Coffea (Rubiaceae). PLoS One. 2011;6(6):e21019.
Roalson EH, McCubbin AG. S-RNases and sexual incompatibility: structure, functions, and evolutionary perspectives. Mol Phylogenet Evol. 2003;29(3):490–506.
Vieira J, Ferreira PG, Aguiar B, Fonseca NA, Vieira CP. Evolutionary patterns at the RNase based gametophytic self - incompatibility system in two divergent Rosaceae groups (Maloideae and Prunus). BMC Evol Biol. 2010;10:200.
Ushijima K, Sassa H, Dandekar AM, Gradziel TM, Tao R, Hirano H. Structural and transcriptional analysis of the self-incompatibility locus of almond: identification of a pollen-expressed F-box gene with haplotype-specific polymorphism. Plant Cell. 2003;15(3):771–81.
Entani T, Iwano M, Shiba H, Che FS, Isogai A, Takayama S. Comparative analysis of the self-incompatibility (S-) locus region of Prunus mume: identification of a pollen-expressed F-box gene with allelic diversity. Genes Cells. 2003;8(3):203–13.
Ikeda K, Igic B, Ushijima K, Yamane H, Hauck N, Nakano R, et al. Primary structural features of the S haplotype-specific F-box protein, SFB, in Prunus. Sex Plant Reprod. 2004;16(5):235–43.
Sonneveld T, Tobutt KR, Vaughan SP, Robbins TP. Loss of pollen-S function in two self-compatible selections of Prunus avium is associated with deletion/mutation of an S haplotype-specific F-box gene. Plant Cell. 2005;17(1):37–51.
Nunes MD, Santos RA, Ferreira SM, Vieira J, Vieira CP. Variability patterns and positively selected sites at the gametophytic self-incompatibility pollen SFB gene in a wild self-incompatible Prunus spinosa (Rosaceae) population. New Phytol. 2006;172(3):577–87.
Vieira J, Santos RA, Ferreira SM, Vieira CP. Inferences on the number and frequency of S-pollen gene (SFB) specificities in the polyploid Prunus spinosa. Heredity. 2008;101(4):351–8.
Cheng J, Han Z, Xu X, Li T. Isolation and identification of the pollen-expressed polymorphic F-box genes linked to the S-locus in apple (Malus × domestica). Sex Plant Reprod. 2006;19(4):175–83.
Kakui H, Tsuzuki T, Koba T, Sassa H. Polymorphism of SFBB-gamma and its use for S genotyping in Japanese pear (Pyrus pyrifolia). Plant Cell Rep. 2007;26(9):1619–25.
Sassa H, Kakui H, Miyamoto M, Suzuki Y, Hanada T, Ushijima K, et al. S locus F-Box brothers: multiple and pollen-specific F-box genes with S haplotype-specific polymorphisms in apple and Japanese pear. Genetics. 2007;175(4):1869–81.
Minamikawa M, Kakui H, Wang S, Kotoda N, Kikuchi S, Koba T, et al. Apple S locus region represents a large cluster of related, polymorphic and pollen-specific F-box genes. Plant Mol Biol. 2010;74(1–2):143–54.
De Franceschi P, Pierantoni L, Dondini L, Grandi M, Sansavini S, Sanzol J. Evaluation of candidate F-box genes for the pollen S of gametophytic self-incompatibility in the Pyrinae (Rosaceae) on the basis of their phylogenomic context. Tree Genet Genomes. 2011;7(4):663–83.
Kakui H, Kato M, Ushijima K, Kitaguchi M, Kato S, Sassa H. Sequence divergence and loss-of-function phenotypes of S locus F-box brothers genes are consistent with non-self recognition by multiple pollen determinants in self-incompatibility of Japanese pear (Pyrus pyrifolia). Plant J. 2011;68(6):1028–38.
Okada K, Tonaka N, Taguchi T, Ichikawa T, Sawamura Y, Nakanishi T, et al. Related polymorphic F-box protein genes between haplotypes clustering in the BAC contig sequences around the S-RNase of Japanese pear. J Exp Bot. 2011;62(6):1887–902.
Aguiar B, Vieira J, Cunha AE, Fonseca NA, Reboiro-Jato D, Reboiro-Jato M, et al. Patterns of evolution at the gametophytic self-incompatibility Sorbus aucuparia (Pyrinae) S pollen genes support the non-self recognition by multiple factors model. J Exp Bot. 2013;64(8):2423–34.
Wheeler D, Newbigin E. Expression of 10 S-class SLF-like genes in Nicotiana alata pollen and its implications for understanding the pollen factor of the S locus. Genetics. 2007;177(4):2171–80.
Kubo K, Entani T, Takara A, Wang N, Fields AM, Hua Z, et al. Collaborative non-self recognition system in S-RNase-based self-incompatibility. Science. 2010;330(6005):796–9.
Williams JS, Der JP, Kao T-h. Transcriptome analysis reveals the same 17 S-Locus F-Box genes in two haplotypes of the self-incompatibility locus of Petunia inflata. The Plant Cell Online. 2014;26(7):2873–88.
Luu D-T, Qin X, Laublin G, Yang Q, Morse D, Cappadocia M. Rejection of S-heteroallelic pollen by a dual-specific S-RNase in Solanum chacoense predicts a multimeric SI pollen xomponent. Genetics. 2001;159(1):329–35.
Huang J, Zhao L, Yang Q, Xue Y. AhSSK1, a novel SKP1‐like protein that interacts with the S‐locus F‐box protein SLF. The Plant J. 2006;46(5):780–93.
Hua Z, Kao TH. Identification and characterization of components of a putative Petunia S-locus F-box-containing E3 ligase complex involved in S-RNase-based self-incompatibility. Plant Cell. 2006;18(10):2531–53.
Zhao L, Huang J, Zhao Z, Li Q, Sims TL, Xue Y. The Skp1‐like protein SSK1 is required for cross‐pollen compatibility in S‐RNase‐based self‐incompatibility. The Plant J. 2010;62(1):52–63.
Xu C, Li M, Wu J, Guo H, Li Q, Zhang Y, et al. Identification of a canonical SCFSLF complex involved in S-RNase-based self-incompatibility of Pyrus (Rosaceae). Plant Mol Biol. 2013;81(3):245–57.
Matsumoto D, Tao R. Yeast Two-Hybrid screening for the general inhibitor detoxifying S-RNase in Prunus. Acta Hortic. 2012;967:167–70.
Lavin M, Herendeen PS, Wojciechowski MF. Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the tertiary. Syst Biol. 2005;54(4):575–94.
Tsuchimatsu T, Suwabe K, Shimizu-Inatsugi R, Isokawa S, Pavlidis P, Städler T, et al. Evolution of self-compatibility in Arabidopsis by a mutation in the male specificity gene. Nature. 2010;464(7293):1342–6.
Ngo BX, Wakana A, Kim JH, Mori T, Sakai K. Estimation of self-incompatibility S genotypes of Citrus cultivars and plants based on controlled pollination with restricted number of pollen grains. J Fac Agric Kyushu Univ. 2010;55(1):67–72.
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(9):1439–51.
Roiz L, Goren R, Shoseyov O. Stigmatic RNase in calamondin (Citrus reticulata var. austera x Fortunella sp.). Physiol Plantarum. 1995;94(4):585–90.
H-x M, Y-h Q, da Silva JA T, Ye Z-x, Hu G-b. Cloning and expression analysis of S-RNase homologous gene in Citrus reticulata Blanco cv. Wuzishatangju Plant Sci. 2011;180(2):358–67.
Chai L, Ge X, Xu Q, Deng X. CgSL2, an S-like RNase gene in ‘Zigui shatian’pummelo (Citrus grandis Osbeck), is involved in ovary senescence. Mol Biol Rep. 2011;38(1):1–8.
Miao H-X, Qin Y-H, Ye Z-X, Hu G-B. Molecular characterization and expression analysis of ubiquitin-activating enzyme E1 gene in Citrus reticulata. Gene. 2013;513(2):249–59.
Caruso M, Merelo P, Distefano G, La Malfa S, Piero ARL, Tadeo FR, et al. Comparative transcriptome analysis of stylar canal cells identifies novel candidate genes implicated in the self-incompatibility response of Citrus clementina. BMC Plant Biol. 2012;12(1):20.
Ford CS, Wilkinson MJ. Confocal observations of late-acting self-incompatibility in Theobroma cacao L. Sex Plant Reprod. 2012;25(3):169–83.
Tate JA, Simpson BB. Breeding system evolution in Tarasa (Malvaceae) and selection for reduced pollen grain size in the polyploid species. Am J Bot. 2004;91(2):207–13.
Tao R, Watari A, Hanada T, Habu T, Yaegaki H, Yamaguchi M, et al. Self-compatible peach (Prunus persica) has mutant versions of the S haplotypes found in self-incompatible Prunus species. Plant Mol Biol. 2007;63(1):109–23.
Broothaerts W, Janssens GA, Proost P, Broekaert WF. cDNA cloning and molecular analysis of two self-incompatibility alleles from apple. Plant Mol Biol. 1995;27(3):499–511.
Vieira J, Morales-Hojas R, Santos RA, Vieira CP. Different positively selected sites at the gametophytic self-incompatibility pistil S-RNase gene in the Solanaceae and Rosaceae (Prunus, Pyrus, and Malus). J Mol Evol. 2007;65(2):175–85.
Matsumoto D, Yamane H, Tao R. Characterization of SLFL1, a pollen-expressed F-box gene located in the Prunus S locus. Sex Plant Reprod. 2008;21(2):113–21.
Sassa H, Kakui H, Minamikawa M. Pollen-expressed F-box gene family and mechanism of S-RNase-based gametophytic self-incompatibility (GSI) in Rosaceae. Sex Plant Reprod. 2010;23(1):39–43.
Dhar R, Sharma N, Sharma B. Ovule abortion in relation to breeding system in four Trifolium species. Curr Sci. 2006;91(4):482–5.
Ištvánek J, Jaroš M, Křenek A, Řepková J. Genome assembly and annotation for red clover (Trifolium pratense; Fabaceae). Am J Bot. 2014;101(2):327–37.
Lawrence M. Number of incompatibility alleles in clover and other species. Heredity. 1996;76(6):610–5.
Vilanova S, Badenes ML, Burgos L, Martínez-Calvo J, Llácer G, Romero C. Self-compatibility of two apricot selections is associated with two pollen-part mutations of different nature. Plant Physiol. 2006;142(2):629–41.
Zuriaga E, Muñoz-Sanz JV, Molina L, Gisbert AD, Badenes ML, Romero C. An S-Locus independent pollen factor confers self-compatibility in ‘Katy’ Apricot. PLoS One. 2013;8(1):e53947.
Bechsgaard JS, Castric V, Charlesworth D, Vekemans X, Schierup MH. The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr. Mol Biol Evol. 2006;23(9):1741–50.
Boggs NA, Nasrallah JB, Nasrallah ME. Independent S-locus mutations caused self-fertility in Arabidopsis thaliana. PLoS Genet. 2009;5(3):e1000426.
Asmussen C. Pollination biology of the sea pea, Lathyrus japonicus: floral characters and activity and flight patterns of bumblebees. Flora. 1993;188(2):227–37.
López J, Rodríguez-Riaño T, Ortega-Olivencia A, Devesa JA, Ruiz T. Pollination mechanisms and pollen-ovule ratios in some Genisteae (Fabaceae) from Southwestern Europe. Plant Syst Evol. 1999;216(1–2):23–47.
Rodríguez-Riaño T. Biología floral y reproductiva en Fabaceae de Extremadura. Badajoz, Spain: Universidad de Extremadura; 1997.
Rodet G, Vaissière BE, Brévault T, Grossa J-PT. Status of self-pollen in bee pollination efficiency of white clover (Trifolium repens L.). Oecologia. 1998;114(1):93–9.
Gibbs P, Sassaki R. Reproductive biology of Dalbergia miscolobium Benth. (Leguminosae-Papilionoideae) in SE Brazil: the effects of pistillate sorting on fruit-set. Ann Bot. 1998;81(6):735–40.
Rodríguez-Riaño T, Ortega-Olivencia A, Devesa JA. Reproductive biology in Cytisus multiflorus (Fabaceae). In: Annales Botanici Fennici: 2004. Helsinki: Societas Biologica Fennica Vanamo; 1964. p. 179–88.
Cooper D, Brink R. Somatoplastic sterility as a cause of seed failure after interspecific hybridization. Genetics. 1940;25(6):593–617.
Miri R, Bubar J. Self-incompatibility as an outcrossing mechanism in birdsfoot trefoil (Lotus corniculatus). Can J Plant Sci. 1966;46(4):411–8.
Seavey SR, Bawa KS. Late-acting self-incompatibility in Angiosperms. Bot Rev. 1986;52(2):195–219.
Valtueña FJ, Rodríguez-Riaño T, Espinosa F, Ortega-Olivencia A. Self-sterility in two Cytisus species (Leguminosae, Papilionoideae) due to early-acting inbreeding depression. Am J Bot. 2010;97(1):123–35.
Waser NM, Price MV. Reproductive costs of self-pollination in Ipomopsis aggregata (Polemoniaceae) - are ovules usurped? Am J Bot. 1991;78(8):1036–43.
Gibbs PE, Bianchi M. Postpollination events in species of Chorisia (Bombacaceae) and Tabebuia (Bignoniaceae) with late-acting self-incompatibility. Bot Acta. 1993;106(1):64–71.
Rice P, Longden I, Bleasby A. EMBOSS: the European molecular biology open software suite. Trends Genet. 2000;16(6):276–7.
Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36 suppl 2:W5–9.
Young ND, Debellé F, Oldroyd GE, Geurts R, Cannon SB, Udvardi MK, et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011;480(7378):520–4.
Varshney RK, Song C, Saxena RK, Azam S, Yu S, Sharpe AG, et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat Biotech. 2013;31(3):240–6.
Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, et al. Genome sequence of the palaeopolyploid soybean. Nature. 2010;463(7278):178–83.
Yang H, Tao Y, Zheng Z, Shao D, Li Z, Sweetingham MW, et al. Rapid development of molecular markers by next-generation sequencing linked to a gene conferring phomopsis stem blight disease resistance for marker-assisted selection in lupin (Lupinus angustifolius L.) breeding. Theor Appl Genet. 2013;126(2):511–22.
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, de Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40(W1):W597–603.
Reboiro-Jato D, Reboiro-Jato M, Fdez-Riverola F, Vieira CP, Fonseca NA, Vieira J. ADOPS - Automatic Detection Of Positively Selected Sites. J Integr Bioinform. 2012;9(3):200.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.
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(10):2731–9.
Ingram GC, Doyle S, Carpenter R, Schultz EA, Simon R, Coen ES. Dual role for fimbriata in regulating floral homeotic genes and cell division in Antirrhinum. Embo J. 1997;16(21):6521–34.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19(18):2496–7.
Benedito VA, Torres‐Jerez I, Murray JD, Andriankaja A, Allen S, Kakar K, et al. A gene expression atlas of the model legume Medicago truncatula. The Plant J. 2008;55(3):504–13.
Rodríguez-Riaño T, Valtueña FJ, Ortega-Olivencia A. Megasporogenesis, megagametogenesis and ontogeny of the aril in Cytisus striatus and C. multiflorus (Leguminosae: Papilionoideae). Ann Bot. 2006;98(4):777–91.
Blankenberg D, Kuster GV, Coraor N, Ananda G, Lazarus R, Mangan M, et al. Galaxy: a web‐based genome analysis tool for experimentalists. Current protocols in molecular biology. 2010;19:10 11–21.
Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, et al. Galaxy: a platform for interactive large-scale genome analysis. Genome Res. 2005;15(10):1451–5.
Goecks J, Nekrutenko A, Taylor J. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11(8):R86.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.
Roberts A, Pachter L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods. 2013;10(1):71–3.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.
This work was funded by FEDER funds through the Operational Competitiveness Programme – COMPETE and by National Funds through FCT – Fundação para a Ciência e a Tecnologia under the project FCOMP-01-0124-FEDER-028299 (PTDC/BIA-BIC/3830/2012). Bruno Aguiar is the recipient of a PhD grant (SFRH/BD/69207/2010) from FCT.
The authors declare that they have no competing interests.
BA, JV, and CPV performed the bioinformatics and phylogenetic analyses, and drafted the manuscript. BA and AEC performed RNA extractions, cDNA synthesis, and PCR amplifications. JV and CPV conceived, coordinated, and supervised the experiments. All authors have read and approved the final manuscript.
Bruno Aguiar and Jorge Vieira contributed equally to this work.
Bayesian phylogenetic trees showing the relationship of SSK1 like genes in flowering plants. Sequences were aligned using ClustalW2 (A), and T-coffee (B) algorithms. The tree was rooted using O. sativa ([GenBank:AP003824]) and C. maxima ([GenBank:FJ851401]) genes. Numbers below the branches represent posterior credibility values above 60.
Bayesian phylogenetic trees showing the relationship of Fabaceae S-RNase lineage genes and Prunus , Pyrinae, Solanaceae and Plantaginaceae S-RNases. Sequences were aligned using ClustalW2 (A), and T-coffee (B) algorithms. The reference sequences are shaded.
Fabaceae T2-RNases available in GenBank not presenting amino acid pattern 4.
Reads from the SRP033257 experiment of M. truncatula (RNA-Seq data sets from a mixed sample of M. truncatula root knot galls infected with Meloidogyne hapla (a plant-nematode)) supporting the expression of the Mt3 , Mt17 , Mt18 , and Mt_10 genes.
Bayesian phylogenetic trees showing the relationship of the F-box SFB - SFBB - and SLFL - like genes surrounding the C. arietinum Ca4 , M. truncatula Mt3 , Mt17 , Mt18 , and Mt20 genes, and S- pollen genes from Prunus , Malus , Solanaceae and Plantaginaceae, and Prunus S -like genes (shaded sequences). Sequences were aligned using ClustalW2 (A), and T-coffee (B) algorithms. The tree was rooted using A. thaliana F-box/kelch-repeat ([GenBank:NM111499]) gene. Numbers below the branches represent posterior credibility values above 60.
Bayesian phylogenetic trees, showing the relationship of the C. striatus CsRNase1 and CsRNase2 genes with Fabaceae S-RNase lineage genes and Prunus , Pyrinae, Solanaceae and Plantaginaceae S-RNases (shaded sequences). Sequences were aligned using ClustalW2 (A), and T-coffee (B) algorithms. Numbers below the branches represent posterior credibility values above 60.
The 100 most expressed genes of the C. striatus stigma with style transcriptome.
Representation of the genomic region of C. striatus CsRNase 1 (A) and CsRNase 2 (B) genes. Lines represent introns, grey boxes represent exons, and arrows indicate the most external primers used.
Correspondences between gene codes and gene names for F-box SFBB- and SFB - like genes in the vicinity of C. arietinum and M. truncatula S-RNase like genes.
Primers used in this work.
About this article
Cite this article
Aguiar, B., Vieira, J., Cunha, A.E. et al. No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes. BMC Plant Biol 15, 129 (2015). https://doi.org/10.1186/s12870-015-0497-2
- Gametophytic self-incompatibility
- Molecular evolution
- S-RNase like genes
- Trifolium pratense
- Medicago truncatula
- Cicer arietinum
- Cytisus striatus