S-RNases and sexual incompatibility: structure, functions, and evolutionary perspectives

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

Evolutionary and structural aspects of Solanaceae RNases T2

Plant RNases T2 are involved in several physiological and developmental processes, including inorganic phosphate starvation, senescence, wounding, defense against pathogens, and the self-incompatibility system. Solanaceae RNases form three main clades, one composed exclusively of S-RNases and two that include S-like RNases. We identified several positively selected amino acids located in highly flexible regions of these molecules, mainly close to the B1 and B2 substrate-binding sites in S-like RNases and the hypervariable regions of S-RNases. These differences between S- and S-like RNases in the flexibility of amino acids in substrate-binding regions are essential to understand the RNA-binding process. For example, in the S-like RNase NT, two positively selected amino acid residues (Tyr156 and Asn134) are located at the most flexible sites on the molecular surface. RNase NT is induced in response to tobacco mosaic virus infection; these sites may thus be regions of interaction with pathogen proteins or viral RNA. Differential selective pressures acting on plant ribonucleases have increased amino acid variability and, consequently, structural differences within and among S-like RNases and S-RNases that seem to be essential for these proteins play different functions.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

Fine-tuning the transcriptional regulatory model of adaptation response to phosphate stress in maize (Zea mays L.)

The post green revolution agriculture is based on generous application of fertilizers and high-yielding genotypes that are suited for such high input regimes. Cereals, like maize (Zea mays L.) are capable of utilizing less than 20% of the applied inorganic phosphate (Pi) - a non-renewable fertilizer resource. A greater understanding of the molecular mechanisms underlying the acquisition, transportation and utilization of Pi may lead to engineering genotypes with high phosphorus use efficiency. In this study, we carried out functional domain similarity analysis, promoter analysis and comparative transcriptional expression profiling of 12 selected Pi responsive genes in the Pi stress tolerant maize inbred line HKI-163 under sufficient and deficient Pi conditions. Pi starvation led to significant increase in root length; marked proliferation of root hairs and lesser number of crown roots. Eleven genes were significantly up or down regulated in Pi deficient condition. The putative acid phosphatase, ZmACP5 expression was up regulated by 162.81 and 74.40 fold in root and leaf tissues, respectively. The RNase, ZmRNS1 showed 115 fold up regulation in roots under Pi deprivation. Among the two putative high affinity Pi transporters ZmPht1;4 was found specific to root, whereas ZmPht2 was found to be up regulated in both root and leaf tissues. The genes involved in Pi homeostasis pathway (ZmSIZ1, SPX1 and Pho2) were up regulated in root and leaf. In light of the expression profiling of selected regulatory genes, an updated model of transcriptional regulation under Pi starvation in maize has been presented.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01155-x.

Genome-wide identification and characterization of legume T2 Ribonuclease gene family and analysis of GmaRNS9, a soybean T2 Ribonuclease gene, function in nodulation

T2 ribonuclease family (RNaseT2) proteins are secretory and nonspecific endoribonucleases that have a large conserved biological role. Family members of RNaseT2 are found in every organism and carry out important biological functions. However, little is known about the functions of these proteins in legumes, including potential roles in symbiotic nodulation. This study aimed to characterize and perform bioinformatic analysis of RNaseT2 genes in four legume species that their genome was sequenced. In total, 60 RNaseT2 genes were identified and characterized. By analyzing their phylogeny, we divided these RNaseT2 into five clades. Expression analysis of RNaseT2 genes indicated that these genes are expressed in various tissues, and the most expression level was related to the pod, flower, and root. Moreover, GmaRNS9 expression analysis in soybean was consistent with in silico studies and demonstrated that this gene usually has high root tip expression. GmaRNS9 expression was reduced by Bradyrhizobium japonicum inoculation and nodule formation. Reduced expression of this gene was possibly controlled by the GmNARK gene either directly or pleiotropically through increased phosphorus requirements during increased nodulation. However, the nutrient stress (phosphate and nitrate starvation) led to an increase in the expression level of GmRNS9. In silico and quantitative gene expression analyses showed that RNaseT2 genes could play important roles in the growth and development of legumes as well as nodulation.

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

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

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.

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

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

The apple MdPTI1L kinase is phosphorylated by MdOXI1 during S-RNase-induced reactive oxygen species signaling in pollen tubes

Apple (Malus domestica) exhibits classic S-RNase-mediated gametophytic self-incompatibility. Previous studies have shown that the S-RNase secreted from style cells could trigger signal transduction and defense responses mediated by Ca²⁺ and reactive oxygen species (ROS) after entering into the pollen tube. In this study, we investigated the downstream genes activated by ROS during S-RNase-mediated gametophytic self-incompatibility in pollen tubes. A substantial increase in ROS, as well as up-regulated expression of a serine-threonine protein kinase gene, OXIDATIVE SIGNAL-INDUCIBLE1 (MdOXI1), was detected in apple pollen tubes treated with self-S-RNase. A kinase assay-linked phosphoproteomics (KALIP) analysis suggested that MdOXI1 could bind and phosphorylate the downstream protein kinase Pto-interacting protein 1-like (MdPTI1L). The phosphorylation level of MdPTI1L was significantly reduced after silencing MdOXI1 with antisense oligonucleotides in the pollen tube. Silencing of either MdOXI1 or MdPTI1L alleviated the inhibitory effect of self-S-RNase on pollen tube growth. Our results thus indicate that MdPTI1L is phosphorylated by MdOXI1 in the pollen tube and participates in the ROS signaling pathway triggered by S-RNase.

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.

Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (F_ROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that F_ROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.

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

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

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.

Finding a Compatible Partner: Self-Incompatibility in European Pear (Pyrus communis); Molecular Control, Genetic Determination, and Impact on Fertilization and Fruit Set

Pyrus species display a gametophytic self-incompatibility (GSI) system that actively prevents fertilization by self-pollen. The GSI mechanism in Pyrus is genetically controlled by a single locus, i.e., the S-locus, which includes at least two polymorphic and strongly linked S-determinant genes: a pistil-expressed S-RNase gene and a number of pollen-expressed SFBB genes (S-locus F-Box Brothers). Both the molecular basis of the SI mechanism and its functional expression have been widely studied in many Rosaceae fruit tree species with a particular focus on the characterization of the elusive SFBB genes and S-RNase alleles of economically important cultivars. Here, we discuss recent advances in the understanding of GSI in Pyrus and provide new insights into the mechanisms of GSI breakdown leading to self-fertilization and fruit set. Molecular analysis of S-genes in several self-compatible Pyrus cultivars has revealed mutations in both pistil- or pollen-specific parts that cause breakdown of self-incompatibility. This has significantly contributed to our understanding of the molecular and genetic mechanisms that underpin self-incompatibility. Moreover, the existence and development of self-compatible mutants open new perspectives for pear production and breeding. In this framework, possible consequences of self-fertilization on fruit set, development, and quality in pear are also reviewed.

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.

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.

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

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

Allelic diversity of S-RNase alleles in diploid potato species

The S-ribonuclease sequences of 16 S-alleles derived from diploid types of Solanum are presented. A phylogenetic analysis and partial phenotypic analysis support the conclusion that these are functional S-alleles. S-Ribonucleases (S-RNases) control the pistil specificity of the self-incompatibility (SI) response in the genus Solanum and several other members of the Solanaceae. The nucleotide sequences of S-RNases corresponding to a large number of S-alleles or S-haplotypes have been characterised. However, surprisingly, few S-RNase sequences are available for potato species. The identification of new S-alleles in diploid potato species is desirable as these stocks are important sources of traits such as biotic and abiotic resistance. S-RNase sequences are reported here from three distinct diploid types of potato: cultivated Solanum tuberosum Group Phureja, S. tuberosum Group Stenotomum, and the wild species Solanum okadae. Partial S-RNase sequences were obtained from pistil RNA by RT-PCR or 3'RACE (Rapid Amplification of cDNA Ends) using a degenerate primer. Full-length sequences were obtained for two alleles by 5'RACE. Database searches with these sequences identified 16 S-RNases in total, all of which are novel. The sequence analysis revealed all the expected features of functional S-RNases. Phylogenetic analysis with selected published S-RNase and S-like-RNase sequences from the Solanaceae revealed extensive trans-generic evolution of the S-RNases and a clear distinction from S-like-RNases. Pollination tests were used to confirm the self-incompatibility status and cross-compatibility relationships of the S. okadae accessions. All the S. okadae accessions were found to be self-incompatible as expected with crosses amongst them exhibiting both cross-compatibility and semi-compatibility consistent with the S-genotypes determined from the S-RNase sequence data. The progeny analysis of four semi-compatible crosses examined by allele-specific PCR provided further confirmation that these are functional S-RNases.

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.

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.

SCF(SLF)-mediated cytosolic degradation of S-RNase is required for cross-pollen compatibility in S-RNase-based self-incompatibility in Petunia hybrida

Many flowering plants adopt self-incompatibility (SI) to maintain their genetic diversity. In species of Solanaceae, Plantaginaceae, and Rosaceae, SI is genetically controlled by a single S-locus with multiple haplotypes. The S-locus has been shown to encode S-RNases expressed in pistil and multiple SLF (S-locus F-box) proteins in pollen controlling the female and male specificity of SI, respectively. S-RNases appear to function as a cytotoxin to reject self-pollen. In addition, SLFs have been shown to form SCF (SKP1/Cullin1/F-box) complexes to serve as putative E3 ubiquitin ligase to interact with S-RNases. Previously, two different mechanisms, the S-RNase degradation and the S-RNase compartmentalization, have been proposed as the restriction mechanisms of S-RNase cytotoxicity allowing compatible pollination. In this study, we have provided several lines of evidence in support of the S-RNase degradation mechanism by a combination of cellular, biochemical and molecular biology approaches. First, both immunogold labeling and subcellular fractionation assays showed that two key pollen SI factors, PhS3L-SLF1 and PhSSK1 (SLF-interacting SKP1-like1) from Petunia hybrida, a Solanaceous species, are co-localized in cytosols of both pollen grains and tubes. Second, PhS3L-RNases are mainly detected in the cytosols of both self and non-self-pollen tubes after pollination. Third, we found that PhS-RNases selectively interact with PhSLFs by yeast two-hybrid and co-immunoprecipitation assays. Fourth, S-RNases are specifically degraded in compatible pollen tubes by non-self SLF action. Taken together, our results demonstrate that SCF(SLF-mediated) non-self S-RNase degradation occurs in the cytosol of pollen tube through the ubiquitin/26S proteasome system serving as the major mechanism to neutralize S-RNase cytotoxicity during compatible pollination in P. hybrida.

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.

S-RNase-based self-incompatibility in Petunia inflata

BACKGROUND

For the Solanaceae-type self-incompatibility, also possessed by Rosaceae and Plantaginaceae, the specificity of self/non-self interactions between pollen and pistil is controlled by two polymorphic genes at the S-locus: the S-locus F-box gene (SLF or SFB) controls pollen specificity and the S-RNase gene controls pistil specificity.

SCOPE

This review focuses on the work from the authors' laboratory using Petunia inflata (Solanaceae) as a model. Here, recent results on the identification and functional studies of S-RNase and SLF are summarized and a protein-degradation model is proposed to explain the biochemical mechanism for specific rejection of self-pollen tubes by the pistil.

CONCLUSIONS

The protein-degradation model invokes specific degradation of non-self S-RNases in the pollen tube mediated by an SLF, and can explain compatible versus incompatible pollination and the phenomenon of competitive interaction, where SI breaks down in pollen carrying two different S-alleles. In Solanaceae, Plantaginaceae and subfamily Maloideae of Rosaceae, there also exist multiple S-locus-linked SLF/SFB-like genes that potentially function as the pollen S-gene. To date, only three such genes, all in P. inflata, have been examined, and they do not function as the pollen S-gene in the S-genotype backgrounds tested. Interestingly, subfamily Prunoideae of Rosaceae appears to possess only a single SLF/SFB gene, and competitive interaction, observed in Solanaceae, Plantaginaceae and subfamily Maloideae, has not been observed. Thus, although the cytotoxic function of S-RNase is an integral part of SI in Solanaceae, Plantaginaceae and Rosaceae, the function of SLF/SFB may have diverged. This highlights the complexity of the S-RNase-based SI mechanism. The review concludes by discussing some key experiments that will further advance our understanding of this self/non-self discrimination mechanism.

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.

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.

Zebrafish RNase T2 genes and the evolution of secretory ribonucleases in animals

Background

Members of the Ribonuclease (RNase) T2 family are common models for enzymological studies, and their evolution has been well characterized in plants. This family of acidic RNases is widespread, with members in almost all organisms including plants, animals, fungi, bacteria and even some viruses. While several biological functions have been proposed for these enzymes in plants, their role in animals is unknown. Interestingly, in vertebrates most of the biological roles of plant RNase T2 proteins are carried out by members of a different family, RNase A. Still, RNase T2 proteins are conserved in these animals

Results

As a first step to shed light on the role of animal RNase T2 enzymes, and to understand the evolution of these proteins while co-existing with the RNase A family, we characterized RNase Dre1 and RNase Dre2, the two RNase T2 genes present in the zebrafish (Danio rerio) genome. These genes are expressed in most tissues examined, including high expression in all stages of embryonic development, and their expression corresponds well with the presence of acidic RNase activities in every tissue analyzed. Embryo expression seems to be a conserved characteristic of members of this family, as other plant and animal RNase T2 genes show similar high expression during embryo development. While plant RNase T2 proteins and the vertebrate RNase A family show evidences of radiation and gene sorting, vertebrate RNase T2 proteins form a monophyletic group, but there is also another monophyletic group defining a fish-specific RNase T2 clade.

Conclusion

Based on gene expression and phylogenetic analyses we propose that RNase T2 enzymes carry out a housekeeping function. This conserved biological role probably kept RNase T2 enzymes in animal genomes in spite of the presence of RNases A. A hypothetical role during embryo development is also discussed.

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.

Compatible pollinations in Solanum chacoense decrease both S-RNase and S-RNase mRNA

Gametophytic self-incompatibility (GSI) allows plants to block fertilization by haploid pollen whose S-allele constitution matches one of the two S-alleles in the diploid styles. GSI in Solanum chacoense requires a stylar S-RNase, first secreted from cells of the transmitting tract then imported into incompatible (self) pollen tubes. However, the molecular mechanisms allowing compatible pollen to evade S-RNase attack are less clear, as compatible pollen tubes also import S-RNase. Using styles of the same age and size in order to lower the degree of inter-style variability in S-RNase levels, we observe reduction of up to 30% of the total non-self stylar S-RNase in vivo during compatible crosses, whereas no degradation of self S-RNases is detected. This marked difference in stylar S-RNase levels dovetails with measurements of pollen-specific Lat52 mRNA, which decreases four-fold in incompatible compared to compatible crosses. Unexpectedly, we also find evidence for a reciprocal signaling mechanism from compatible pollen to the cells of the transmitting tract that results in a roughly three-fold decrease in S-RNase transcript levels. These findings reveal a previously unsuspected feedback loop that may help reinforce the compatible reaction.

Basic RNases of wild almond (Prunus webbii): cloning and characterization of six new S-RNase and one "non-S RNase" genes

In order to investigate the S-RNase allele structure of a Prunus webbii population from the Montenegrin region of the Balkans, we analyzed 10 Prunus webbii accessions. We detected 10 different S-RNase allelic variants and obtained the nucleotide sequences for six S-RNases. The BLAST analysis showed that these six sequences were new Prunus webbii S-RNase alleles. It also revealed that one of sequenced alleles, S(9)-RNase, coded for an amino acid sequence identical to that for Prunus dulcis S(14)-RNase, except for a single conservative amino acid replacement in the signal peptide region. Another, S(3)-RNase, was shown to differ by only three amino acid residues from Prunus salicina Se-RNase. The allele S(7)-RNase was found to be inactive by stylar protein isoelectric focusing followed by RNase-specific staining, but the reason for the inactivity was not at the coding sequence level. Further, in five of the 10 analyzed accessions, we detected the presence of one active basic RNase (marked PW(1)) that did not amplify with S-RNase-specific DNA primers. However, it was amplified with primers designed from the PA1 RNase nucleotide sequence (basic "non-S RNase" of Prunus avium) and the obtained sequence showed high homology (80%) with the PA1 allele. Although homologs of PA1 "non-S RNases" have been reported in four other Prunus species, this is the first recorded homolog in Prunus webbii. The evolutionary implications of the data are discussed.

The Prunus self-incompatibility locus (S locus) is seldom rearranged

Self-incompatibility enables flowering plants to discriminate between self- and non-selfpollen. In Prunus, the 2 genes determining specificity are the S-RNase (the female determinant that is a glycoprotein with ribonuclease activity) and the SFB (the male determinant, a protein with an F-box motif). In all Prunus S haplotypes characterized so far, with the exception of Prunus armeniaca S(2) haplotype, the 2 genes have opposite transcription orientations. Nevertheless, the relative transcription orientation observed in P. armeniaca S(2) haplotype has been postulated to be the one present in all S haplotypes from this species. We show that this is not the case by demonstrating that that the relative transcription orientation of the pollen and pistil genes of the P. armeniaca S(17) haplotype is that which is commonly found in Prunus. Using a phylogenetic approach, we show that the relative transcription orientation of the S-RNase and SFB genes is seldom changed (less than once every 380 million years). This contrasts with the Brassica sporophytic S locus where chromosomal rearrangements are often observed in the region between the pollen and pistil genes.

Identification of major lysine residues of S(3)-RNase of Petunia inflata involved in ubiquitin-26S proteasome-mediated degradation in vitro

S-RNase-based self-incompatibility has been identified in three flowering plant families, including the Solanaceae, and this self/non-self recognition mechanism between pollen and pistil is controlled by two polymorphic genes at the S-locus, S-RNase and S-locus F-box (SLF). S-RNase is produced in the pistil and taken up by pollen tubes in a non-S-haplotype-specific manner. How an allelic product of SLF interacts with self and non-self S-RNases to result in growth inhibition of self pollen tubes is not completely understood. One model predicts that SLF targets non-self S-RNases for ubiquitin/26S proteasome-mediated degradation, thereby only allowing self S-RNase to exert cytotoxic activity inside a pollen tube. To test this model, we studied whether any of the 20 lysine residues in S(3)-RNase of Petunia inflata might be targets for ubiquitination. We identified six lysines near the C-terminus for which mutation to arginine significantly reduced ubiquitination and degradation of the mutant S(3)-RNase, GST:S(3)-RNase (K141-164R) in pollen tube extracts. We further showed that GST:S(3)-RNase (K141-164R) and GST:S(3)-RNase had similar RNase activity, suggesting that their degradation was probably not caused by an ER-associated protein degradation pathway that removes mis-folded proteins. Finally, we showed that PiSBP1 (P. inflata S-RNase binding protein 1), a potential RING-HC subunit of the PiSLF (P. inflata SLF)-containing E3-like complex, could target S-RNase for ubiquitination in vitro. All these results suggest that ubiquitin/26S proteasome-dependent degradation of S-RNase may be an integral part of the S-RNase-based self-incompatibility mechanism.

Trans-specific S-RNase and SFB alleles in Prunus self-incompatibility haplotypes

Self-incompatibility in the genus Prunus is controlled by two genes at the S-locus, S-RNase and SFB. Both genes exhibit the high polymorphism and high sequence diversity characteristic of plant self-incompatibility systems. Deduced polypeptide sequences of three myrobalan and three domestic plum S-RNases showed over 97% identity with S-RNases from other Prunus species, including almond, sweet cherry, Japanese apricot and Japanese plum. The second intron, which is generally highly polymorphic between alleles was also remarkably well conserved within these S-allele pairs. Degenerate consensus primers were developed and used to amplify and sequence the co-adapted polymorphic SFB alleles. Sequence comparisons also indicated high degrees of polypeptide sequence identity between three myrobalan and the three domestic plum SFB alleles and the corresponding Prunus SFB alleles. We discuss these trans-specific allele identities in terms of S-allele function, evolution of new allele specificities and Prunus taxonomy and speciation.

Different positively selected sites at the gametophytic self-incompatibility pistil S-RNase gene in the Solanaceae and Rosaceae (Prunus, Pyrus, and Malus)

In this work we perform a comparative study on the location of positively selected sites (those likely responsible for defining specificity differences) at the S-RNase gene, the pistil component of the gametophytic self-incompatibility system. For Plantaginaceae and Rosaceae (Prunus and Pyrus/Malus) this is the first study of this kind. A clear sign of positive selection was observed for 13, 17, and 27 amino acid sites in Solanaceae, Prunus, and Pyrus/Malus, respectively, using two different methodologies. In Plantaginaceae no clear positively selected sites were identified. Possible reasons for this result are discussed. Indirect experimental evidence suggests that the identified positively selected amino acid sites play a role in specificity determination. The percentage of positively selected sites is similar in Solanaceae and Rosaceae but the location of those sites is different.

Self-incompatibility of Prunus tenella and evidence that reproductively isolated species of Prunus have different SFB alleles coupled with an identical S-RNase allele

Many species of Prunus display an S-RNase-based gametophytic self-incompatibility (SI), controlled by a single highly polymorphic multigene complex termed the S-locus. This comprises tightly linked stylar- and pollen-expressed genes that determine the specificity of the SI response. We investigated SI of Prunus tenella, a wild species found in small, isolated populations on the Balkan peninsula, initially by pollination experiments and identifying stylar-expressed RNase alleles. Nine P. tenella S-RNase alleles (S(1)-S(9)) were cloned; their sequence analysis showed a very high ratio of non-synonymous to synonymous nucleotide substitutions (K(a)/K(s)) and revealed that S-RNase alleles of P. tenella, unlike those of Prunus dulcis, show positive selection in all regions except the conserved regions and that between C2 and RHV. Remarkably, S(8)-RNase, was found to be identical to S(1)-RNase from Prunus avium, a species that does not interbreed with P. tenella and, except for just one amino acid, to S(11) of P. dulcis. However, the corresponding introns and S-RNase-SFB intergenic regions showed considerable differences. Moreover, protein sequences of the pollen-expressed SFB alleles were not identical, harbouring 12 amino-acid replacements between those of P. tenella SFB(8) and P. avium SFB(1). Implications of this finding for hypotheses about the evolution of new S-specificities are discussed.

A nonself recognition gene complex in Neurospora crassa

Nonself recognition is exemplified in the fungal kingdom by the regulation of cell fusion events between genetically different individuals (heterokaryosis). The het-6 locus is one of approximately 10 loci that control heterokaryon incompatibility during vegetative growth of N. crassa. Previously, it was found that het-6-associated incompatibility in Oak Ridge (OR) strains involves two contiguous genes, het-6 and un-24. The OR allele of either gene causes "strong" incompatibility (cell death) when transformed into Panama (PA)-background strains. Several remarkable features of the locus include the nature of these incompatibility genes (het-6 is a member of a repetitive gene family and un-24 also encodes the large subunit of ribonucleotide reductase) and the observation that un-24 and het-6 are in severe linkage disequilibrium. Here, we identify "weak" (slow, aberrant growth) incompatibility activities by un-24PA and het-6PA when transformed separately into OR strains, whereas together they exhibit an additive, strong effect. We synthesized strains with the new allelic combinations un-24PA het-6OR and un-24OR het-6PA, which are not found in nature. These strains grow normally and have distinct nonself recognition capabilities but may have reduced fitness. Comparing the Oak Ridge and Panama het-6 regions revealed a paracentric inversion, the architecture of which provides insights into the evolution of the un-24-het-6 gene complex.

Characterization of RNASET2, the first human member of the Rh/T2/S family of glycoproteins

Ribonucleases are ubiquitous enzymes involved in RNA metabolism and are classified in several families on the basis of their structural, catalytic, and biological properties. Here, we describe characterization of the only human member of the Rh/T2/S family of acid hydrolases so far described, named RNASET2. This protein was previously reported to have an interesting biological function in the control of tumourigenesis and metastatization. We show that RNASET2 is present in multiple forms in human cell lines and mouse tissues, one of which represents the full length, glycosylated and secreted form, while the others are proteolytic products. RNASET2 is endowed with catalytic activity as demonstrated with purified recombinant protein expressed in the Baculovirus Expression Vector System and in a human cell line ectopically expressing various types of constructs. Furthermore, we document for this protein a lysosomal localization as described for other members of the Rh/T2/S family of ribonucleases. The results presented herein represent a further advancement toward the molecular understanding of the tumour suppressive properties of the human RNASET2 protein.

Evolution under tight linkage to mating type

Recent large-scale sequencing studies of mating type loci in a number of organisms offer insight into the origin and evolution of these genomic regions. Extensive tracts containing genes with a wide diversity of functions typically cosegregate with mating type. Cases in which mating type determination entails complementarity between distinct transcription units may descend from systems in which close physical linkage facilitated the coordinated expression and cosegregation of the interacting genes. In response to the particular selection pressures associated with the maintenance of more than one mating type, this nucleus of low recombination may expand over evolutionary time, engulfing neighboring tracts bearing genes with no direct role in reproduction. This scenario is consistent with the present-day structure of some mating type loci, including regulators of homomorphic self-incompatibility in angiosperms (S-loci). Recombination suppression and enforced S-locus heterozygosity accelerate the accumulation of genetic load and promote genetic associations between S-alleles and degenerating genes in cosegregating tracts. This S-allele-specific load may influence the evolution of self-incompatibility systems.

Analysis of the S-locus structure in Prunus armeniaca L. Identification of S-haplotype specific S-RNase and F-box genes

The gametophytic self-incompatibility (GSI) system in Rosaceae has been proposed to be controlled by two genes located in the S -locusan S-RNase and a recently described pollen expressed S -haplotype specific F-box gene (SFB). However, in apricot (Prunus armeniaca L.) these genes had not been identified yet. We have sequenced 21 kb in total of the S -locus region in 3 different apricot S -haplotypes. These fragments contain genes homologous to the S-RNase and F-box genes found in other Prunus species, preserving their basic gene structure features and defined amino acid domains. The physical distance between the F-box and the S-RNase genes was determined exactly in the S2-haplotype (2.9 kb) and inferred approximately in the S 1-haplotype (< 49 kb) confirming that these genes are linked. Sequence analysis of the 5' flanking regions indicates the presence of a conserved region upstream of the putative TATA box in the S-RNase gene. The three identified S-RNase alleles (S1, S2 and S4) had a high allelic sequence diversity (75.3 amino acid identity), and the apricot F-box allelic variants (SFB1, SFB2 and SFB4) were also highly haplotype-specific (79.4 amino acid identity). Organ specific-expression was also studied, revealing that S1- and S2-RNases are expressed in style tissues, but not in pollen or leaves. In contrast, SFB1 and SFB2 are only expressed in pollen, but not in styles or leaves. Taken together, these results support these genes as candidates for the pistil and pollen S-determinants of GSI in apricot.