An S-RNase-based gametophytic self-incompatibility system evolved only once in eudicots

Drivers of diversification in Linum (Linaceae) by means of chromosome evolution: correlations with biogeography, breeding system and habit

BACKGROUND AND AIMS

Chromosome evolution leads to hybrid dysfunction and recombination patterns and has thus been proposed as a major driver of diversification in all branches of the tree of life, including flowering plants. In this study we used the genus Linum (flax species) to evaluate the effects of chromosomal evolution on diversification rates and on traits that are important for sexual reproduction. Linum is a useful study group because it has considerable reproductive polymorphism (heterostyly) and chromosomal variation (n = 6-36) and a complex pattern of biogeographical distribution.

METHODS

We tested several traditional hypotheses of chromosomal evolution. We analysed changes in chromosome number across the phylogenetic tree (ChromEvol model) in combination with diversification rates (ChromoSSE model), biogeographical distribution, heterostyly and habit (ChromePlus model).

KEY RESULTS

Chromosome number evolved across the Linum phylogeny from an estimated ancestral chromosome number of n = 9. While there were few apparent incidences of cladogenesis through chromosome evolution, we inferred up to five chromosomal speciation events. Chromosome evolution was not related to heterostyly but did show significant relationships with habit and geographical range. Polyploidy was negatively correlated with perennial habit, as expected from the relative commonness of perennial woodiness and absence of perennial clonality in the genus. The colonization of new areas was linked to genome rearrangements (polyploidy and dysploidy), which could be associated with speciation events during the colonization process.

CONCLUSIONS

Chromosome evolution is a key trait in some clades of the Linum phylogeny. Chromosome evolution directly impacts speciation and indirectly influences biogeographical processes and important plant traits.

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

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

The origin and evolution of RNase T2 family and gametophytic self-incompatibility system in plants

Ribonuclease (RNase) T2 genes are found widely in both eukaryotes and prokaryotes, and genes from this family have been revealed to have various functions in plants. In particular, S-RNase is known to be the female determinant in the S-RNase-based gametophytic self-incompatibility (GSI) system. However, the origin and evolution of the RNase T2 gene family and GSI system are not well understood. In this study, 785 RNase T2 genes were identified in 81 sequenced plant genomes representing broad-scale diversity and divided into three subgroups (Class I, II, and III) based on phylogenetic and synteny network analysis. Class I was found to be of ancient origin and to emerge in green algae, Class II was shown to originate with the appearance of angiosperms, while Class III was discovered to be eudicot-specific. Each of the three major classes could be further classified into several subclasses of which some subclasses were found to be lineage-specific. Furthermore, duplication, deletion, or inactivation of the S/S-like-locus was revealed to be linked to repeated loss and gain of self-incompatibility in different species from distantly related plant families with GSI. Finally, the origin and evolutionary history of S-locus in Rosaceae species was unraveled with independent loss and gain of S-RNase occurred in different subfamilies of Rosaceae. Our findings provide insights into the origin and evolution of the RNase T2 family and the GSI system in plants.

Origin, loss, and regain of self-incompatibility in angiosperms

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

RNase-based self-incompatibility in cacti

Approximately one-half of all flowering plants express genetically based physiological mechanisms that prevent self-fertilisation. One such mechanism, termed RNase-based self-incompatibility, employs ribonucleases as the pistil component. Although it is widespread, it has only been characterised in a handful of distantly related families, partly due to the difficulties presented by life history traits of many plants, which complicate genetic research. Many species in the cactus family are known to express self-incompatibility but the underlying mechanisms remain unknown. We demonstrate the utility of a candidate-based RNA-seq approach, combined with some unusual features of self-incompatibility-causing genes, which we use to uncover the genetic basis of the underlying mechanisms. Specifically, we assembled transcriptomes from Schlumbergera truncata (crab cactus or false Christmas cactus), and interrogated them for tissue-specific expression of candidate genes, structural characteristics, correlation with expressed phenotype(s), and phylogenetic placement. The results were consistent with operation of the RNase-based self-incompatibility mechanism in Cactaceae. The finding yields additional evidence that the ancestor of nearly all eudicots possessed RNase-based self-incompatibility, as well as a clear path to better conservation practices for one of the most charismatic plant families.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Genome-wide identification and expression pattern analysis of the ribonuclease T2 family in Eucommia ulmoides

The 2′,3′-cycling ribonuclease (RNase) genes are catalysts of RNA cleavage and include the RNase T2 gene family. RNase T2 genes perform important roles in plants and have been conserved in the genome of eukaryotic organisms. In this study we identified 21 EURNS genes in Eucommia ulmoides Oliver (E. ulmoides) and analyzed their structure, chromosomal location, phylogenetic tree, gene duplication, stress-related cis-elements, and expression patterns in different tissues. The length of 21 predicted EURNS proteins ranged from 143 to 374 amino acids (aa), their molecular weight (MW) ranged from 16.21 to 42.38 kDa, and their isoelectric point (PI) value ranged from 5.08 to 9.09. Two classifications (class I and class III) were obtained from the conserved domains analysis and phylogenetic tree. EURNS proteins contained a total of 15 motifs. Motif 1, motif 2, motif 3, and motif 7 were distributed in multiple sequences and were similar to the conserved domain of RNase T2. EURNS genes with similar structure and the predicted EURNS proteins with conserved motif compositions are in the same group in the phylogenetic tree. The results of RT-PCR and transcription data showed that EURNS genes have tissue-specific expression and exhibited obvious trends in different developmental stages. Gene duplication analysis results indicated that segment duplication may be the dominant duplication mode in this gene family. This study provides a theoretical basis for research on the RNase T2 gene family and lays a foundation for the further study of EURNS genes.

Association of T2/S-RNase With Self-Incompatibility of Japanese Citrus Accessions Examined by Transcriptomic, Phylogenetic, and Genetic Approaches

Several citrus varieties show gametophytic self-incompatibility (GSI), which can contribute to seedless fruit production in several cultivars. This study investigated the genes regulating this trait through RNA-seq performed using styles collected from the flowers of Japanese citrus cultivars 'Hyuganatsu,' 'Tosabuntan,' 'Hassaku,' 'Banpeiyu,' and 'Sweet Spring'. We screened the transcripts of putative T2 RNases, i.e., the protein family including all S-RNases from S-RNase-based GSI plants, and constructed a phylogenetic tree using the screened T2 RNases and S-RNases retrieved from citrus genome databases and a public database. Three major clusters (class I-III) were formed, among which, the class III cluster contained family specific subclusters formed by S-RNase and a citrus-specific cluster monophyletic to the S-RNase clusters. From the citrus class III cluster, six transcripts were consistent with the S haplotypes previously determined in Japanese citrus accessions, sharing characteristics such as isoelectric point, extracellular localization, molecular weight, intron number and position, and tissue-specific expression with S-RNases. One T2 RNase gene in self-incompatible Hyuganatsu was significantly down-regulated in the styles of a self-compatible mutant of Hyuganatsu in RNA-seq and qPCR analyses. In addition, the inheritance pattern of some T2 RNase genes was consistent with the pattern of the S haplotype in the progeny population of Hyuganatsu and Tosabuntan. As all results supported citrus self-incompatibility being based on S-RNase, we believe that six T2 RNase genes were S-RNases. The homology comparison between the six T2 RNases and S-RNases recently reported in Chinese citrus revealed that three out of six T2 RNases were identical to S-RNases from Chinese citrus. Thus, the other three T2 RNases were finally concluded to be novel citrus S-RNases involved in self-incompatibility.

The identification of the Rosa S-locus and implications on the evolution of the Rosaceae gametophytic self-incompatibility systems

In Rosaceae species, two gametophytic self-incompatibility (GSI) mechanisms are described, the Prunus self-recognition system and the Maleae (Malus/Pyrus/Sorbus) non-self- recognition system. In both systems the pistil component is a S-RNase gene, but from two distinct phylogenetic lineages. The pollen component, always a F-box gene(s), in the case of Prunus is a single gene, and in Maleae there are multiple genes. Previously, the Rosa S-locus was mapped on chromosome 3, and three putative S-RNase genes were identified in the R. chinensis ‘Old Blush’ genome. Here, we show that these genes do not belong to the S-locus region. Using R. chinensis and R. multiflora genomes and a phylogenetic approach, we identified the S-RNase gene, that belongs to the Prunus S-lineage. Expression patterns support this gene as being the S-pistil. This gene is here also identified in R. moschata, R. arvensis, and R. minutifolia low coverage genomes, allowing the identification of positively selected amino acid sites, and thus, further supporting this gene as the S-RNase. Furthermore, genotype–phenotype association experiments also support this gene as the S-RNase. For the S-pollen GSI component we find evidence for multiple F-box genes, that show the expected expression pattern, and evidence for diversifying selection at the F-box genes within an S-haplotype. Thus, Rosa has a non-self-recognition system, like in Maleae species, despite the S-pistil gene belonging to the Prunus S-RNase lineage. These findings are discussed in the context of the Rosaceae GSI evolution. Knowledge on the Rosa S-locus has practical implications since genes controlling floral and other ornamental traits are in linkage disequilibrium with the S-locus.

SKP1-like protein, CrSKP1-e, interacts with pollen-specific F-box proteins and assembles into SCF-type E3 complex in 'Wuzishatangju' (Citrus reticulata Blanco) pollen

S-ribonuclease (S-RNase)-based self-incompatibility (SI) mechanisms have been extensively studied in Solanaceae, Rosaceae and Plantaginaceae. S-RNase-based SI is controlled by two closely related genes, S-RNase and S-locus F-box (SLF), located at a polymorphic S-locus. In the SI system, the SCF-type (SKP1-CUL1-F-box-RBX1) complex functions as an E3 ubiquitin ligase complex for ubiquitination of non-self S-RNase. Pummelo (Citrus grandis) and several mandarin cultivars are suggested to utilize an S-RNase-based SI system. However, the molecular mechanism of the non-S-factors involved in the SI reaction is not straightforward in Citrus. To investigate the SCF-type E3 complex responsible for the SI reaction in mandarin, SLF, SKP1-like and CUL1 candidates potentially involved in the SI reaction of ‘Wuzishatangju’ (Citrus reticulata Blanco) were identified based on the genome-wide identification and expression analyses. Sixteen pollen-specific F-box genes (CrFBX1-CrFBX16), one pollen-specific SKP1-like gene (CrSKP1-e) and two CUL1 genes (CrCUL1A and CrCUL1B) were identified and cloned from ‘Wuzishatangju’. Yeast two-hybrid (Y2H) and in vitro binding assays showed that five CrFBX proteins could bind to CrSKP1-e, which is an ortholog of SSK1 (SLF-interacting-SKP1-like), a non-S-factor responsible for the SI reaction. Luciferase complementation imaging (LCI) and in vitro binding assays also showed that CrSKP1-e interacts with the N-terminal region of both CrCUL1A and CrCUL1B. These results indicate that CrSKP1-e may serve as a functional member of the SCF-type E3 ubiquitin ligase complex in ‘Wuzishatangju’.

Identification of a Novel SBP1-Containing SCFSFB Complex in Wild Dwarf Almond (Prunus tenella)

S-RNase-based gametophytic self-incompatibility (SI), in which specificities of pistil and pollen are determined by S-RNase and the S locus F-box protein, respectively, has been discovered in the Solanaceae, Plantaginaceae, and Rosaceae families, but some underlying molecular mechanisms remain elusive and controversial. Previous studies discovered SI in wild dwarf almond (Prunus tenella), and pistil S (S-RNase) and pollen S (SFB) determinant genes have been investigated. However, the SCF (SKP1–Cullin1–F-box-Rbx1) complex, which serves as an E3 ubiquitin ligase on non-self S-RNase, has not been investigated. In the current study, PetSSK1 (SLF-interacting-SKP1-like1), SBP1 (S-RNase binding protein 1), CUL1, and SFB genes (S-haplotype-specific F-box) were identified in an accession (ZB1) of P. tenella. Yeast two-hybrid assays revealed interactions between PetSBP1 and PetCUL1 and between PetSBP1 and PetSFBs (SFB16 and SFB17), and subsequent pull-down assays confirmed these interactions, suggesting a novel SBP1-containing SCF^SFB complex in wild dwarf almond. Moreover, despite a putative interaction between PetSSK1 and PetCUL1, we revealed that PetSSK1 does not interact with PetSFB16 or PetSFB17, and thus the canonical SSK1-containing SCF^SFB complex could not be identified. This suggests a novel molecular mechanism of gametophytic SI in Prunus species.

Predicting Specificities Under the Non-self Gametophytic Self-Incompatibility Recognition Model

Non-self gametophytic self-incompatibility (GSI) recognition system is characterized by the presence of multiple F-box genes tandemly located in the S-locus, that regulate pollen specificity. This reproductive barrier is present in Solanaceae, Plantaginacea and Maleae (Rosaceae), but only in Petunia functional assays have been performed to get insight on how this recognition mechanism works. In this system, each of the encoded S-pollen proteins (called SLFs in Solanaceae and Plantaginaceae /SFBBs in Maleae) recognizes and interacts with a sub-set of non-self S-pistil proteins, called S-RNases, mediating their ubiquitination and degradation. In Petunia there are 17 SLF genes per S-haplotype, making impossible to determine experimentally each SLF specificity. Moreover, domain -swapping experiments are unlikely to be performed in large scale to determine S-pollen and S-pistil specificities. Phylogenetic analyses of the Petunia SLFs and those from two Solanum genomes, suggest that diversification of SLFs predate the two genera separation. Here we first identify putative SLF genes from nine Solanum and 10 Nicotiana genomes to determine how many gene lineages are present in the three genera, and the rate of origin of new SLF gene lineages. The use of multiple genomes per genera precludes the effect of incompleteness of the genome at the S-locus. The similar number of gene lineages in the three genera implies a comparable effective population size for these species, and number of specificities. The rate of origin of new specificities is one per 10 million years. Moreover, here we determine the amino acids positions under positive selection, those involved in SLF specificity recognition, using 10 Petunia S-haplotypes with more than 11 SLF genes. These 16 amino acid positions account for the differences of self-incompatible (SI) behavior described in the literature. When SLF and S-RNase proteins are divided according to the SI behavior, and the positively selected amino acids classified according to hydrophobicity, charge, polarity and size, we identified fixed differences between SI groups. According to the in silico 3D structure of the two proteins these amino acid positions interact. Therefore, this methodology can be used to infer SLF/S-RNase specificity recognition.

Recognition of S-RNases by an S locus F-box like protein and an S haplotype-specific F-box like protein in the Prunus-specific self-incompatibility system

S-RNase was demonstrated to be predominantly recognized by an S locus F-box-like protein and an S haplotype-specific F-box-like protein in compatible pollen tubes of sweet cherry. Self-incompatibility (SI) is a reproductive barrier that rejects self-pollen and inhibits self-fertilization to promote outcrossing. In Solanaceae and Rosaceae, S-RNase-based gametophytic SI (GSI) comprises S-RNase and F-box protein(s) as the pistil and pollen S determinants, respectively. Compatible pollen tubes are assumed to detoxify the internalized cytotoxic S-RNases to maintain growth. S-RNase detoxification is conducted by the Skp1-cullin1-F-box protein complex (SCF) formed by pollen S determinants, S locus F-box proteins (SLFs), in Solanaceae. In Prunus, the general inhibitor (GI), but not pollen S determinant S haplotype-specific F-box protein (SFB), is hypothesized to detoxify S-RNases. Recently, SLF-like proteins 1-3 (SLFL1-3) were suggested as GI candidates, although it is still possible that other proteins function predominantly in GI. To identify the other GI candidates, we isolated four other pollen-expressed SLFL and SFB-like (SFBL) proteins PavSLFL6, PavSLFL7A, PavSFBL1, and PavSFBL2 in sweet cherry. Binding assays with four PavS-RNases indicated that PavSFBL2 bound to PavS^{1, 6}-RNase while the others bound to nothing. PavSFBL2 was confirmed to form an SCF complex in vitro. A co-immunoprecipitation assay using the recombinant PavS⁶-RNase as bait against pollen extracts and a mass spectrometry analysis identified the SCF complex components of PavSLFLs and PavSFBL2, M-locus-encoded glutathione S-transferase (MGST), DnaJ-like protein, and other minor proteins. These results suggest that SLFLs and SFBLs could act as predominant GIs in Prunus-specific S-RNase-based GSI.

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Inferences on specificity recognition at the Malus×domestica gametophytic self-incompatibility system

In Malus × domestica (Rosaceae) the product of each SFBB gene (the pollen component of the gametophytic self-incompatibility (GSI) system) of a S-haplotype (the combination of pistil and pollen genes that are linked) interacts with a sub-set of non-self S-RNases (the pistil component), but not with the self S-RNase. To understand how the Malus GSI system works, we identified 24 SFBB genes expressed in anthers, and determined their gene sequence in nine M. domestica cultivars. Expression of these SFBBs was not detected in the petal, sepal, filament, receptacle, style, stigma, ovary or young leaf. For all SFBBs (except SFBB15), identical sequences were obtained only in cultivars having the same S-RNase. Linkage with a particular S-RNase was further established using the progeny of three crosses. Such data is needed to understand how other genes not involved in GSI are affected by the S-locus region. To classify SFBBs specificity, the amino acids under positive selection obtained when performing intra-haplotypic analyses were used. Using this information and the previously identified S-RNase positively selected amino acid sites, inferences are made on the S-RNase amino acid properties (hydrophobicity, aromatic, aliphatic, polarity, and size), at these positions, that are critical features for GSI specificity determination.

The evolutionary history of plant T2/S-type ribonucleases

A growing number of T2/S-RNases are being discovered in plant genomes. Members of this protein family have a variety of known functions, but the vast majority are still uncharacterized. We present data and analyses of phylogenetic relationships among T2/S-RNases, and pay special attention to the group that contains the female component of the most widespread system of self-incompatibility in flowering plants. The returned emphasis on the initially identified component of this mechanism yields important conjectures about its evolutionary context. First, we find that the clade involved in self-rejection (class III) is found exclusively in core eudicots, while the remaining clades contain members from other vascular plants. Second, certain features, such as intron patterns, isoelectric point, and conserved amino acid regions, help differentiate S-RNases, which are necessary for expression of self-incompatibility, from other T2/S-RNase family members. Third, we devise and present a set of approaches to clarify new S-RNase candidates from existing genome assemblies. We use genomic features to identify putative functional and relictual S-loci in genomes of plants with unknown mechanisms of self-incompatibility. The widespread occurrence of possible relicts suggests that the loss of functional self-incompatibility may leave traces long after the fact, and that this manner of molecular fossil-like data could be an important source of information about the history and distribution of both RNase-based and other mechanisms of self-incompatibility. Finally, we release a public resource intended to aid the search for S-locus RNases, and help provide increasingly detailed information about their taxonomic distribution.

Lack of S-RNase-Based Gametophytic Self-Incompatibility in Orchids Suggests That This System Evolved after the Monocot-Eudicot Split

Self-incompatibility (SI) is found in approximately 40% of flowering plant species and at least 100 families. Although orchids belong to the largest angiosperm family, only 10% of orchid species present SI and have gametophytic SI (GSI). Furthermore, a majority (72%) of Dendrobium species, which constitute one of the largest Orchidaceae genera, show SI and have GSI. However, nothing is known about the molecular mechanism of GSI. The S-determinants of GSI have been well characterized at the molecular level in Solanaceae, Rosaceae, and Plantaginaceae, which use an S-ribonuclease (S-RNase)-based system. Here, we investigate the hypothesis that Orchidaceae uses a similar S-RNase to those described in Rosaceae, Solanaceae, and Plantaginaceae SI species. In this study, two SI species (Dendrobium longicornu and D. chrysanthum) were identified using fluorescence microscopy. Then, the S-RNase- and SLF-interacting SKP1-like1 (SSK1)-like genes present in their transcriptomes and the genomes of Phalaenopsis equestris, D. catenatum, Vanilla shenzhenica, and Apostasia shenzhenica were investigated. Sequence, phylogenetic, and tissue-specific expression analyses revealed that none of the genes identified was an S-determinant, suggesting that Orchidaceae might have a novel SI mechanism. The results also suggested that RNase-based GSI might have evolved after the split of monocotyledons (monocots) and dicotyledons (dicots) but before the split of Asteridae and Rosidae. This is also the first study to investigate S-RNase-based GSI in monocots. However, studies on gene identification, differential expression, and segregation analyses in controlled crosses are needed to further evaluate the genes with high expression levels in GSI tissues.

The Jujube Genome Provides Insights into Genome Evolution and the Domestication of Sweetness/Acidity Taste in Fruit Trees

Jujube (Ziziphus jujuba Mill.) belongs to the Rhamnaceae family and is a popular fruit tree species with immense economic and nutritional value. Here, we report a draft genome of the dry jujube cultivar 'Junzao' and the genome resequencing of 31 geographically diverse accessions of cultivated and wild jujubes (Ziziphus jujuba var. spinosa). Comparative analysis revealed that the genome of 'Dongzao', a fresh jujube, was ~86.5 Mb larger than that of the 'Junzao', partially due to the recent insertions of transposable elements in the 'Dongzao' genome. We constructed eight proto-chromosomes of the common ancestor of Rhamnaceae and Rosaceae, two sister families in the order Rosales, and elucidated the evolutionary processes that have shaped the genome structures of modern jujubes. Population structure analysis revealed the complex genetic background of jujubes resulting from extensive hybridizations between jujube and its wild relatives. Notably, several key genes that control fruit organic acid metabolism and sugar content were identified in the selective sweep regions. We also identified S-locus genes controlling gametophytic self-incompatibility and investigated haplotype patterns of the S locus in the jujube genomes, which would provide a guideline for parent selection for jujube crossbreeding. This study provides valuable genomic resources for jujube improvement, and offers insights into jujube genome evolution and its population structure and domestication.

Cullin1-P is an Essential Component of Non-Self Recognition System in Self-Incompatibility in Petunia

Self-incompatibility (SI) in flowering plants is a genetic reproductive barrier to distinguish self- and non-self pollen to promote outbreeding. In Solanaceae, self-pollen is rejected by the ribonucleases expressed in the styles (S-RNases), via its cytotoxic function. On the other side, the male-determinant is the S-locus F-box proteins (SLFs) expressed in pollen. Multiple SLFs collaboratively detoxify non-self S-RNases, therefore, non-self recognition is the mode of self-/non-self discrimination in Solanaceae. It is considered that SLFs function as a substrate-recognition module of the Skp1-Cullin1-F-box (SCF) complex that inactivates non-self S-RNases via their polyubiquitination, which leads to degradation by 26S proteasome. In fact, PhSSK1 (Petunia hybrida SLF-interacting Skp1-like1) was identified as a specific component of SCF^SLF and was shown to be essential for detoxification of S-RNase in Petunia However, different molecules are proposed as the candidate Cullin1, another component of SCF^SLF, and there is as yet no definite conclusion. Here, we identified five Cullin1s from the expressed sequence tags (ESTs) derived from the male reproductive organ in Petunia Among them, only PhCUL1-P was co-immunoprecipitated with S₇-SLF2. In vitro protein-binding assay suggested that PhSSK1 specifically forms a complex with PhCUL1-P in an SLF-dependent manner. Knockdown of PhCUL1-P suppressed fertility of transgenic pollen in cross-compatible pollination in the functional S-RNase-dependent manner. These results suggested that SCF^SLF selectively uses PhCUL1-P. Phylogeny of Cullin1s indicates that CUL1-P is recruited into the SI machinery during the evolution of Solanaceae, suggesting that the SI components have evolved differently among species in Solanaceae and Rosaceae, despite both families sharing the S-RNase-based SI.

Non-self- and self-recognition models in plant self-incompatibility

The mechanisms by which flowering plants choose their mating partners have interested researchers for a long time. Recent findings on the molecular mechanisms of non-self-recognition in some plant species have provided new insights into self-incompatibility (SI), the trait used by a wide range of plant species to avoid self-fertilization and promote outcrossing. In this Review, we compare the known SI systems, which can be largely classified into non-self- or self-recognition systems with respect to their molecular mechanisms, their evolutionary histories and their modes of evolution. We review previous controversies on haplotype evolution in the gametophytic SI system of Solanaceae species in light of a recently elucidated non-self-recognition model. In non-self-recognition SI systems, the transition from self-compatibility (SC) to SI may be more common than previously thought. Reversible transition between SI and SC in plants may have contributed to their adaptation to diverse and fluctuating environments.

Recognition of a wide-range of S-RNases by S locus F-box like 2, a general-inhibitor candidate in the Prunus-specific S-RNase-based self-incompatibility system

Many species in the Rosaceae, the Solanaceae, and the Plantaginaceae exhibit S-RNase-based gametophytic self-incompatibility (GSI). This system comprises S-ribonucleases (S-RNases) as the pistil S determinant and a single or multiple F-box proteins as the pollen S determinants. In Prunus, pollen specificity is determined by a single S haplotype-specific F-box protein (SFB). The results of several studies suggested that SFB exerts cognate S-RNase cytotoxicity, and a hypothetical general inhibitor (GI) is assumed to detoxify S-RNases in non-specific manner unless it is affected by SFB. Although the identity of the GI is unknown, phylogenetic and evolutionary analyses have indicated that S locus F-box like 1-3 (or S locus F-box with low allelic sequence polymorphism 1-3; SLFL1-3), which are encoded by a region of the Prunus genome linked to the S locus, are good GI candidates. Here, we examined the biochemical characteristics of SLFL1-3 to determine whether they have appropriate GI characteristics. Pull-down assays and quantitative expression analyses indicated that Prunus avium SLFL1-3 mainly formed a canonical SCF complex with PavSSK1 and PavCul1A. Binding assays with PavS(1,3,4,6)-RNases showed that PavSLFL1, PavSLFL2, and PavSLFL3 bound to PavS(3)-RNase, all PavS-RNases tested, and none of the PavS-RNases tested, respectively. Together, these results suggested that SLFL2 has the appropriate characteristics to be the GI in sweet cherry pollen, while SLFL1 may redundantly work with SLFL2 to detoxify all S-RNases. We discuss the possible roles of SLFL1-3 as the GI in the Prunus-specific S-RNase-based GSI mechanism.

Detection of Self Incompatibility Genotypes in Prunus africana: Characterization, Evolution and Spatial Analysis

In flowering plants, self-incompatibility is an effective genetic mechanism that prevents self-fertilization. Most Prunus tree species exhibit a homomorphic gametophytic self-incompatibility (GSI) system, in which the pollen phenotype is encoded by its own haploid genome. To date, no identification of S-alleles had been done in Prunus africana, the only member of the genus in Africa. To identify S-RNase alleles and hence determine S-genotypes in African cherry (Prunus africana) from Mabira Forest Reserve, Uganda, primers flanking the first and second intron were designed and these amplified two bands in most individuals. PCR bands on agarose indicated 26 and 8 different S-alleles for second and first intron respectively. Partial or full sequences were obtained for all these fragments. Comparison with published S-RNase data indicated that the amplified products were S-RNase alleles with very high interspecies homology despite the high intraspecific variation. Against expectations for a locus under balancing selection, frequency and spatial distribution of the alleles in a study plot was not random. Implications of the results to breeding efforts in the species are discussed, and mating experiments are strongly suggested to finally prove the functionality of SI in P. africana.

Insights into the Prunus-Specific S-RNase-Based Self-Incompatibility System from a Genome-Wide Analysis of the Evolutionary Radiation of S Locus-Related F-box Genes

Self-incompatibility (SI) is an important plant reproduction mechanism that facilitates the maintenance of genetic diversity within species. Three plant families, the Solanaceae, Rosaceae and Plantaginaceae, share an S-RNase-based gametophytic SI (GSI) system that involves a single S-RNase as the pistil S determinant and several F-box genes as pollen S determinants that act via non-self-recognition. Previous evidence has suggested a specific self-recognition mechanism in Prunus (Rosaceae), raising questions about the generality of the S-RNase-based GSI system. We investigated the evolution of the pollen S determinant by comparing the sequences of the Prunus S haplotype-specific F-box gene (SFB) with those of its orthologs in other angiosperm genomes. Our results indicate that the Prunus SFB does not cluster with the pollen S of other plants and diverged early after the establishment of the Eudicots. Our results further indicate multiple F-box gene duplication events, specifically in the Rosaceae family, and suggest that the Prunus SFB gene originated in a recent Prunus-specific gene duplication event. Transcriptomic and evolutionary analyses of the Prunus S paralogs are consistent with the establishment of a Prunus-specific SI system, and the possibility of subfunctionalization differentiating the newly generated SFB from the original pollen S determinant.

A Comprehensive Study of Molecular Evolution at the Self-Incompatibility Locus of Rosaceae

The family Rosaceae includes a range of important fruit trees, most of which have the S-RNase-based self-incompatibility (SI). Several models have been developed to explain how pollen (SLF) and pistil (S-RNase) components of the S-locus interact. It was discovered in 2010 that additional SLF proteins are involved in pollen specificity, and a Collaborative Non-Self Recognition model has been proposed for SI in Solanaceae; however, the validity of such model remains to be elucidated for other species. The results of this study support the divergent evolution of the S-locus genes from two Rosaceae subfamilies, Prunoideae/Amygdaloideae and Maloideae, The difference identified in the selective pressures between the two lineages provides evidence for positive selection at specific sites in both the S-RNase and the SLF proteins. The evolutionary findings of this study support the role of multiple SLF proteins leading to a Collaborative Non-Self Recognition model for SI in the Maloideae. Furthermore, the identification of the sites responsible for SI specificity determination and the mapping of these sites onto the modelled tertiary structure of ancestor proteins provide useful information for rational functional redesign and protein engineering for the future engineering of new functional alleles providing increased diversity in the SI system in the Maloideae.

NnSR1, a class III non-S-RNase specifically induced in Nicotiana alata under phosphate deficiency, is localized in endoplasmic reticulum compartments

A combined strategy of phosphate (Pi) remobilization from internal and external RNA sources seems to be conserved in plants exposed to Pi starvation. Thus far, the only ribonucleases (RNases) reported to be induced in Nicotiana alata undergoing Pi deprivation are extracellular S-like RNase NE and NnSR1. NnSR1 is a class III non S-RNase of unknown subcellular location. Here, we examine the hypothesis that NnSR1 is an intracellular RNase derived from the self-incompatibility system with specific expression in self-incompatible Nicotiana alata. NnSR1 was not induced in self-compatible Nicotiana species exposed to Pi deprivation. NnSR1 conjugated with a fluorescent protein and transiently expressed in Arabidopsis protoplasts and Nicotiana leaves showed that the fusion protein co-localized with an endoplasmic reticulum (ER) marker. Subcellular fractionation by ultracentrifugation of roots exposed to Pi deprivation revealed that the native NnSR1 migrated in parallel with the BiP protein, a typical ER marker. To our knowledge, NnSR1 is the first class III RNase reported to be localized in ER compartments. The induction of NnSR1 was detected earlier than the extracellular RNase NE, suggesting that intracellular RNA may be the first source of Pi used by the cell under Pi stress.

No evidence for Fabaceae Gametophytic self-incompatibility being determined by Rosaceae, Solanaceae, and Plantaginaceae S-RNase lineage genes

BACKGROUND

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.

RESULTS

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.

CONCLUSION

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.

Convergent evolution at the gametophytic self-incompatibility system in Malus and Prunus

S-RNase-based gametophytic self-incompatibility (GSI) has evolved once before the split of the Asteridae and Rosidae. This conclusion is based on the phylogenetic history of the S-RNase that determines pistil specificity. In Rosaceae, molecular characterizations of Prunus species, and species from the tribe Pyreae (i.e., Malus, Pyrus, Sorbus) revealed different numbers of genes determining S-pollen specificity. In Prunus only one pistil and pollen gene determine GSI, while in Pyreae there is one pistil but multiple pollen genes, implying different specificity recognition mechanisms. It is thus conceivable that within Rosaceae the genes involved in GSI in the two lineages are not orthologous but possibly paralogous. To address this hypothesis we characterised the S-RNase lineage and S-pollen lineage genes present in the genomes of five Rosaceae species from three genera: M. × domestica (apple, self-incompatible (SI); tribe Pyreae), P. persica (peach, self-compatible (SC); Amygdaleae), P. mume (mei, SI; Amygdaleae), Fragaria vesca (strawberry, SC; Potentilleae), and F. nipponica (mori-ichigo, SI; Potentilleae). Phylogenetic analyses revealed that the Malus and Prunus S-RNase and S-pollen genes belong to distinct gene lineages, and that only Prunus S-RNase and SFB-lineage genes are present in Fragaria. Thus, S-RNase based GSI system of Malus evolved independently from the ancestral system of Rosaceae. Using expression patterns based on RNA-seq data, the ancestral S-RNase lineage gene is inferred to be expressed in pistils only, while the ancestral S-pollen lineage gene is inferred to be expressed in tissues other than pollen.

NnSR1, a class III non-S-RNase constitutively expressed in styles, is induced in roots and stems under phosphate deficiency in Nicotiana alata

BACKGROUND AND AIMS

Non-S-ribonucleases (non-S-RNases) are class III T2 RNases constitutively expressed in styles of species with S-RNase-based self-incompatibility. So far, no function has been attributed to these RNases. The aim of this work is to examine if NnSR1, a non-S-RNase from Nicotiana alata, is induced under conditions of phosphate (Pi) deprivation. The hypothesis is that under Pi-limited conditions, non-S-RNase functions may resemble the role of S-like RNases. To date, the only RNases reported to be induced by Pi deficiency are class I and class II S-like RNases, which are phylogenetically different from the class III clade of RNases.

METHODS

Gene and protein expression of NnSR1 were assayed in plants grown hydroponically with and without Pi, by combining RT-PCR, immunoblot and enzymatic activity approaches.

KEY RESULTS

NnSR1 transcripts were detected in roots 7 d after Pi deprivation and remained stable for several days. Transcript expression was correlated based on Pi availability in the culture medium. Antiserum against a peptide based on a hypervariable domain of NnSR1 recognized NnSR1 in roots and stems but not leaves exposed to Pi shortage. NnSR1 was not detected in culture medium and was pelleted with the microsomal fraction, suggesting that it was membrane-associated or included in large compartments. The anti-NnSR1 inhibited selectively the enzymatic activity of a 31-kDa RNase indicating that NnSR1 was induced in an enzymatically active form.

CONCLUSIONS

The induction of NnSR1 indicates that there is a general recruitment of all classes of T2 RNases in response to Pi shortage. NnSR1 appears to have regained ancestral functions of class III RNases related to strategies to cope with Pi limitation and also possibly with other environmental challenges. This constitutes the first report for a specific function of class III RNases other than S-RNases.

Patterns of evolution at the gametophytic self-incompatibility Sorbus aucuparia (Pyrinae) S pollen genes support the non-self recognition by multiple factors model

S-RNase-based gametophytic self-incompatibility evolved once before the split of the Asteridae and Rosidae. In Prunus (tribe Amygdaloideae of Rosaceae), the self-incompatibility S-pollen is a single F-box gene that presents the expected evolutionary signatures. In Malus and Pyrus (subtribe Pyrinae of Rosaceae), however, clusters of F-box genes (called SFBBs) have been described that are expressed in pollen only and are linked to the S-RNase gene. Although polymorphic, SFBB genes present levels of diversity lower than those of the S-RNase gene. They have been suggested as putative S-pollen genes, in a system of non-self recognition by multiple factors. Subsets of allelic products of the different SFBB genes interact with non-self S-RNases, marking them for degradation, and allowing compatible pollinations. This study performed a detailed characterization of SFBB genes in Sorbus aucuparia (Pyrinae) to address three predictions of the non-self recognition by multiple factors model. As predicted, the number of SFBB genes was large to account for the many S-RNase specificities. Secondly, like the S-RNase gene, the SFBB genes were old. Thirdly, amino acids under positive selection-those that could be involved in specificity determination-were identified when intra-haplotype SFBB genes were analysed using codon models. Overall, the findings reported here support the non-self recognition by multiple factors model.

ADOPS - Automatic Detection Of Positively Selected Sites

Maximum-likelihood methods based on models of codon substitution have been widely used to infer positively selected amino acid sites that are responsible for adaptive changes. Nevertheless, in order to use such an approach, software applications are required to align protein and DNA sequences, infer a phylogenetic tree and run the maximum-likelihood models. Therefore, a significant effort is made in order to prepare input files for the different software applications and in the analysis of the output of every analysis. In this paper we present the ADOPS (Automatic Detection Of Positively Selected Sites) software. It was developed with the goal of providing an automatic and flexible tool for detecting positively selected sites given a set of unaligned nucleotide sequence data. An example of the usefulness of such a pipeline is given by showing, under different conditions, positively selected amino acid sites in a set of 54 Coffea putative S-RNase sequences. ADOPS software is freely available and can be downloaded from http://sing.ei.uvigo.es/ADOPS.

Self-sterility in flowering plants: preventing self-fertilization increases family diversification rates

BACKGROUND AND SCOPE

New data are presented on the distribution and frequency of self-sterility (SS) - predominantly pre-zygotic self-incompatibility (SI) systems - in flowering plants and the hypothesis is tested that families with self-sterile taxa have higher net diversification rates (DRs) than those with exclusively self-compatible taxa using both absolute and relative rate tests.

KEY RESULTS

Three major forms of SI systems (where pollen is rejected at the stigmatic, stylar or ovarian interface) are found to occur in the oldest families of flowering plants, with times of divergence >100 million years before the present (mybp), while post-fertilization SS and heterostyly appear in families with crown ages of 81 and 87 mybp, respectively. It is also founnd that many (22) angiosperm families exhibit >1 SI phenotype and that the distribution of different types of SS does not show strong phylogenetic clustering, collectively suggesting that SS and SI systems have evolved repeatedly de novo in angiosperm history. Families bearing self-sterile taxa have higher absolute DRs using all available calibrations of the angiosperm tree, and this affect is caused mostly by the high DR of families with homomorphic SI systems (in particular stigmatic SI) or those in which multiple SS/SI phenotypes have been observed (polymorphic). Lastly, using sister comparisons, it is further demonstrated that in 29 of 38 sister pairs (including 95 families), the self-sterile sister group had higher species richness and DR than its self-compatible sister based on either the total number of taxa in the clade with SS or only the estimated fraction to harbour SS based on literature surveys.

CONCLUSIONS

Collectively, these analyses point to the importance of SS, particularly pre-zygotic SI in the evolution of flowering plants.

Identification of a Skp1-like protein interacting with SFB, the pollen S determinant of the gametophytic self-incompatibility in Prunus

Many species in Rosaceae, Solanaceae, and Plantaginaceae exhibit S-RNase-based self-incompatibility (SI). In this system, the pistil and pollen specificities are determined by S-RNase and the S locus F-box protein, respectively. The pollen S determinant F-box protein in Prunus (Rosaceae) is referred to by two different terms, SFB (for S-haplotype-specific F-box protein) and SLF (for S locus F box), whereas it is called SLF in Solanaceae and Plantaginaceae. Prunus SFB is thought to be a molecule indispensable for its cognate S-RNase to exert cytotoxicity and to arrest pollen tube growth in incompatible reactions. Although recent studies have demonstrated the molecular function of SCF(SLF) in the SI reaction of Solanaceae and Plantaginaceae, how SFB participates in the Prunus SI mechanism remains to be elucidated. Here we report the identification of sweet cherry (Prunus avium) SFB (PavSFB)-interacting Skp1-like1 (PavSSK1) using a yeast (Saccharomyces cerevisiae) two-hybrid screening against the pollen cDNA library. Phylogenetic analysis showed that PavSSK1 belongs to the same clade as Antirrhinum hispanicum SLF-interacting Skp1-like1 and Petunia hybrida SLF-interacting Skp1-like1 (PhSSK1). In yeast, PavSSK1 interacted not only with PavSFBs from different S haplotypes and Cullin1-likes (PavCul1s), but also with S-locus F-box-likes. A pull-down assay confirmed the interactions between PavSSK1 and PavSFB and between PavSSK1 and PavCul1s. These results collectively indicate that PavSSK1 could be a functional component of the SCF complex and that PavSFB may function as a component of the SCF complex. We discuss the molecular function of PavSFB in self-/nonself-recognition in the gametophytic SI of Prunus.

Molecular bases and evolutionary dynamics of self-incompatibility in the Pyrinae (Rosaceae)

The molecular bases of the gametophytic self-incompatibility (GSI) system of species of the subtribe Pyrinae (Rosaceae), such as apple and pear, have been widely studied in the last two decades. The characterization of S-locus genes and of the mechanisms underlying pollen acceptance or rejection have been topics of major interest. Besides the single pistil-side S determinant, the S-RNase, multiple related S-locus F-box genes seem to be involved in the determination of pollen S specificity. Here, we collect and review the state of the art of GSI in the Pyrinae. We emphasize recent genomic data that have contributed to unveiling the S-locus structure of the Pyrinae, and discuss their consistency with the models of self-recognition that have been proposed for Prunus and the Solanaceae. Experimental data suggest that the mechanism controlling pollen-pistil recognition specificity of the Pyrinae might fit well with the collaborative 'non-self' recognition system proposed for Petunia (Solanaceae), whereas it presents relevant differences with the mechanism exhibited by the species of the closely related genus Prunus, which uses a single evolutionarily divergent F-box gene as the pollen S determinant. The possible involvement of multiple pollen S genes in the GSI system of Pyrinae, still awaiting experimental confirmation, opens up new perspectives to our understanding of the evolution of S haplotypes, and of the evolution of S-RNase-based GSI within the Rosaceae family. Whereas S-locus genes encode the players determining self-recognition, pollen rejection in the Pyrinae seems to involve a complex cascade of downstream cellular events with significant similarities to programmed cell death.

Assembly, gene annotation and marker development using 454 floral transcriptome sequences in Ziziphus celata (Rhamnaceae), a highly endangered, Florida endemic plant

Large-scale DNA sequence data may enable development of genetic resources in endangered species, thereby facilitating conservation efforts. Ziziphus celata, a federally endangered, self-incompatible plant species occurring in Florida, USA, is one species for which genetic resources are necessary to facilitate new introductions and augmentations essential for recovery of the species. We used 454 pyrosequencing of a Z. celata normalized floral cDNA library to create a genomic resource for gene and marker discovery. A half-plate GS-FLX Titanium run yielded 655 337 reads averaging 250 bp. A total of 474 025 reads were assembled de novo into 84 645 contigs averaging 408 bp, while 181 312 reads remained unassembled. Forty-seven and 43% of contig consensus sequences had BLAST matches to known proteins in the Uniref50 and TAIR9 annotated protein databases, respectively; many contigs fully represented orthologous proteins in TAIR9. A total of 22 707 unique genes were sequenced, indicating substantial coverage of the Z. celata transcriptome. We detected single-nucleotide polymorphisms and simple sequence repeats (SSRs) and developed thousands of SSR primers for use in future genetic studies. As a first step towards understanding self-incompatibility in Z. celata, we identified sequences belonging to the gene family encoding self-incompatibility. This study demonstrates the efficacy of 454 transcriptome sequencing for rapid gene and marker discovery in an endangered plant.

Expression and trans-specific polymorphism of self-incompatibility RNases in coffea (Rubiaceae)

Self-incompatibility (SI) is widespread in the angiosperms, but identifying the biochemical components of SI mechanisms has proven to be difficult in most lineages. Coffea (coffee; Rubiaceae) is a genus of old-world tropical understory trees in which the vast majority of diploid species utilize a mechanism of gametophytic self-incompatibility (GSI). The S-RNase GSI system was one of the first SI mechanisms to be biochemically characterized, and likely represents the ancestral Eudicot condition as evidenced by its functional characterization in both asterid (Solanaceae, Plantaginaceae) and rosid (Rosaceae) lineages. The S-RNase GSI mechanism employs the activity of class III RNase T2 proteins to terminate the growth of "self" pollen tubes. Here, we investigate the mechanism of Coffea GSI and specifically examine the potential for homology to S-RNase GSI by sequencing class III RNase T2 genes in populations of 14 African and Madagascan Coffea species and the closely related self-compatible species Psilanthus ebracteolatus. Phylogenetic analyses of these sequences aligned to a diverse sample of plant RNase T2 genes show that the Coffea genome contains at least three class III RNase T2 genes. Patterns of tissue-specific gene expression identify one of these RNase T2 genes as the putative Coffea S-RNase gene. We show that populations of SI Coffea are remarkably polymorphic for putative S-RNase alleles, and exhibit a persistent pattern of trans-specific polymorphism characteristic of all S-RNase genes previously isolated from GSI Eudicot lineages. We thus conclude that Coffea GSI is most likely homologous to the classic Eudicot S-RNase system, which was retained since the divergence of the Rubiaceae lineage from an ancient SI Eudicot ancestor, nearly 90 million years ago.

Functional gametophytic self-incompatibility in a peripheral population of Solanum peruvianum (Solanaceae)

The transition from self-incompatibility to self-compatibility is a common transition in angiosperms often reported in populations at the edge of species range limits. Geographically distinct populations of wild tomato species (Solanum section Lycopersicon (Solanaceae)) have been described as polymorphic for mating system with both self-incompatible and self-compatible populations. Using controlled pollinations and sequencing of the S-RNase mating system gene, we test the compatibility status of a population of S. peruvianum located near its southern range limit. Pollinations among plants of known genotypes revealed strong self-incompatibility; fruit set following compatible pollinations was significantly higher than following incompatible pollinations for all tested individuals. Sequencing of the S-RNase gene in parents and progeny arrays was also as predicted under self-incompatibility. Molecular variation at the S-RNase locus revealed a diverse set of alleles, and heterozygosity in over 500 genotyped individuals. We used controlled crosses to test the specificity of sequences recovered in this study; in all cases, results were consistent with a unique allelic specificity for each tested sequence, including two alleles sharing 92% amino-acid similarity. Site-specific patterns of selection at the S-RNase gene indicate positive selection in regions of the gene associated with allelic specificity determination and purifying selection in previously characterized conserved regions. Further, there is broad convergence between the present and previous studies in specific amino-acid positions inferred to be evolving under positive selection.

The genetic location of the self-incompatibility locus in white clover (Trifolium repens L.)

White clover (Trifolium repens L.) is a forage legume of considerable economic importance in temperate agricultural systems. It has a strong self-incompatibility system. The molecular basis of self-incompatibility in T. repens is unknown, but it is under the control of a single locus, which is expressed gametophytically. To locate the self-incompatibility locus (S locus) in T. repens, we carried out cross-pollination experiments in an F(1) mapping population and constructed a genetic linkage map using amplified fragment length polymorphism and simple sequence repeat markers. As the first step in a map-based cloning strategy, we locate for the first time the S locus in T. repens on a genetic linkage map, on the homoeologous linkage group pair 1 (E), which is broadly syntenic to Medicago truncatula L. chromosome 1. On the basis of this syntenic relationship, the possibility that the S locus may or may not possess an S-RNase gene is discussed.

Evolutionary patterns at the RNase based gametophytic self - incompatibility system in two divergent Rosaceae groups (Maloideae and Prunus)

Background

Within Rosaceae, the RNase based gametophytic self-incompatibility (GSI) system has been studied at the molecular level in Maloideae and Prunus species that have been diverging for, at least, 32 million years. In order to understand RNase based GSI evolution within this family, comparative studies must be performed, using similar methodologies.

Result

It is here shown that many features are shared between the two species groups such as levels of recombination at the S-RNase (the S-pistil component) gene, and the rate at which new specificities arise. Nevertheless, important differences are found regarding the number of ancestral lineages and the degree of specificity sharing between closely related species. In Maloideae, about 17% of the amino acid positions at the S-RNase protein are found to be positively selected, and they occupy about 30% of the exposed protein surface. Positively selected amino acid sites are shown to be located on either side of the active site cleft, an observation that is compatible with current models of specificity determination. At positively selected amino acid sites, non-conservative changes are almost as frequent as conservative changes. There is no evidence that at these sites the most drastic amino acid changes may be more strongly selected.

Conclusions

Many similarities are found between the GSI system of Prunus and Maloideae that are compatible with the single origin hypothesis for RNase based GSI. The presence of common features such as the location of positively selected amino acid sites and lysine residues that may be important for ubiquitylation, raise a number of issues that, in principle, can be experimentally addressed in Maloideae. Nevertheless, there are also many important differences between the two Rosaceae GSI systems. How such features changed during evolution remains a puzzling issue.

Petunia nectar proteins have ribonuclease activity

Plants requiring an insect pollinator often produce nectar as a reward for the pollinator's visitations. This rich secretion needs mechanisms to inhibit microbial growth. In Nicotiana spp. nectar, anti-microbial activity is due to the production of hydrogen peroxide. In a close relative, Petunia hybrida, limited production of hydrogen peroxide was found; yet petunia nectar still has anti-bacterial properties, suggesting that a different mechanism may exist for this inhibition. The nectar proteins of petunia plants were compared with those of ornamental tobacco and significant differences were found in protein profiles and function between these two closely related species. Among those proteins, RNase activities unique to petunia nectar were identified. The genes corresponding to four RNase T2 proteins from Petunia hybrida that show unique expression patterns in different plant tissues were cloned. Two of these enzymes, RNase Phy3 and RNase Phy4 are unique among the T2 family and contain characteristics similar to both S- and S-like RNases. Analysis of amino acid patterns suggest that these proteins are an intermediate between S- and S-like RNases, and support the hypothesis that S-RNases evolved from defence RNases expressed in floral parts. This is the first report of RNase activities in nectar.

RNase T2 genes from rice and the evolution of secretory ribonucleases in plants

The plant RNase T2 family is divided into two different subfamilies. S-RNases are involved in rejection of self-pollen during the establishment of self-incompatibility in three plant families. S-like RNases, on the other hand, are not involved in self-incompatibility, and although gene expression studies point to a role in plant defense and phosphate recycling, their biological roles are less well understood. Although S-RNases have been subjects of many phylogenetic studies, few have included an extensive analysis of S-like RNases, and genome-wide analyses to determine the number of S-like RNases in fully sequenced plant genomes are missing. We characterized the eight RNase T2 genes present in the Oryza sativa genome; and we also identified the full complement of RNase T2 genes present in other fully sequenced plant genomes. Phylogenetics and gene expression analyses identified two classes among the S-like RNase subfamily. Class I genes show tissue specificity and stress regulation. Inactivation of RNase activity has occurred repeatedly throughout evolution. On the other hand, Class II seems to have conserved more ancestral characteristics; and, unlike other S-like RNases, genes in this class are conserved in all plant species analyzed and most are constitutively expressed. Our results suggest that gene duplication resulted in high diversification of Class I genes. Many of these genes are differentially expressed in response to stress, and we propose that protein characteristics, such as the increase in basic residues can have a defense role independent of RNase activity. On the other hand, constitutive expression and phylogenetic conservation suggest that Class II S-like RNases may have a housekeeping role.

RNase-based gametophytic self-incompatibility evolution: Questioning the hypothesis of multiple independent recruitments of the S-pollen gene

Multiple independent recruitments of the S-pollen component (always an F-box gene) during RNase-based gametophytic self-incompatibility evolution have recently been suggested. Therefore, different mechanisms could be used to achieve the rejection of incompatible pollen in different plant families. This hypothesis is, however, mainly based on the interpretation of phylogenetic analyses, using a small number of divergent nucleotide sequences. In this work we show, based on a large collection of F-box S-like sequences, that the inferred relationship of F-box S-pollen and F-box S-like sequences is dependent on the sequence alignment software and phylogenetic method used. Thus, at present, it is not possible to address the phylogenetic relationship of F-box S-pollen and S-like sequences from different plant families. In Petunia and Malus/Pyrus the putative S-pollen gene(s) show(s) variability patterns different than expected for an S-pollen gene, raising the question of false identification. Here we show that in Petunia, the unexpected features of the putative S-pollen gene are not incompatible with this gene's being the S-pollen gene. On the other hand, it is very unlikely that the Pyrus SFBB-gamma gene is involved in specificity determination.

Gene classification based on amino acid motifs and residues: the DLX (distal-less) test case

Background

Comparative studies using hundreds of sequences can give a detailed picture of the evolution of a given gene family. Nevertheless, retrieving only the sequences of interest from public databases can be difficult, in particular, when working with highly divergent sequences. The difficulty increases substantially when one wants to include in the study sequences from many (or less well studied) species whose genomes are non-annotated or incompletely annotated.

Methodology/Principal Findings

In this work we evaluate the usefulness of different approaches of gene retrieval and classification, using the distal-less (DLX) gene family as a test case. Furthermore, we evaluate whether the use of a large number of gene sequences from a wide range of animal species, the use of multiple alternative alignments, and the use of amino acids aligned with high confidence only, is enough to recover the accepted DLX evolutionary history.

Conclusions/Significance

The canonical DLX homeobox gene sequence here derived, together with the characteristic amino acid variants here identified in the DLX homeodomain region, can be used to retrieve and classify DLX genes in a simple and efficient way. A program is made available that allows the easy retrieval of synteny information that can be used to classify gene sequences. Maximum likelihood trees using hundreds of sequences can be used for gene identification. Nevertheless, for the DLX case, the proposed DLX evolutionary is not recovered even when multiple alignment algorithms are used.