Identification of a non-S RNase, a possible ancestral form of S-RNases, in Prunus

Auto-phylo v2 and auto-phylo-pipeliner: building advanced, flexible, and reusable pipelines for phylogenetic inferences, estimation of variability levels and identification of positively selected amino acid sites

The vast amount of genome sequence data that is available, and that is predicted to drastically increase in the near future, can only be efficiently dealt with by building automated pipelines. Indeed, the Earth Biogenome Project will produce high-quality reference genome sequences for all 1.8 million named living eukaryote species, providing unprecedented insight into the evolution of genes and gene families, and thus on biological issues. Here, new modules for gene annotation, further BLAST search algorithms, further multiple sequence alignment methods, the adding of reference sequences, further tree rooting methods, the estimation of rates of synonymous and nonsynonymous substitutions, and the identification of positively selected amino acid sites, have been added to auto-phylo (version 2), a recently developed software to address biological problems using phylogenetic inferences. Additionally, we present auto-phylo-pipeliner, a graphical user interface application that further facilitates the creation and running of auto-phylo pipelines. Inferences on S-RNase specificity, are critical for both cross-based breeding and for the establishment of pollination requirements. Therefore, as a test case, we develop an auto-phylo pipeline to identify amino acid sites under positive selection, that are, in principle, those determining S-RNase specificity, starting from both non-annotated Prunus genomes and sequences available in public databases.

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.

Molecular and genetic characterization of a self-compatible apple cultivar, 'CAU-1'

In this study, we characterized a naturally occurring self-compatible apple cultivar, 'CAU-1' (S₁S₉), and studied the underlying mechanism that causes its compatibility. Analyses of both fruit set rate and seed number after self-pollination or cross-pollination with 'Fuji' (S₁S₉), and of pollen tube growth, demonstrated that 'CAU-1' is self-compatible. Genetic analysis by S-RNase PCR-typing of selfed progeny of 'CAU-1' revealed the presence of all progeny classes (S₁S₁, S₁S₉, and S₉S₉). Moreover, no evidence of S-allele duplication was found. These findings support the hypothesis that loss of function of an S-locus unlinked pollen-part mutation (PPM) expressed in pollen, rather than a natural mutation in the pollen-S gene (S₁- and S₉- haplotype), leads to SI breakdown in 'CAU-1'. In addition, there were no significant differences in pollen morphology or fertility between 'Fuji' and 'CAU-1'. However, we found that the effect of S₁- and S₉-RNase on the SI behavior of pollen could not be addressed better in 'CAU-1' than in 'Fuji'. Furthermore, we found that a pollen-expressed hexose transporter, MdHT1, interacted with S-RNases and showed significantly less expression in 'CAU-1' than in 'Fuji' pollen tubes. These findings support the hypothesis that MdHT1 may participate in S-RNase internalization during the SI process, and decrease of MdHT1 expression in 'CAU-1' hindered the release of self S-RNase into the cytoplasm of pollen tubes, thereby protecting pollen from the cytotoxicity of S-RNase, finally probably resulting in self-compatibility. Together, these findings indicate that S-locus external factors are required for gametophytic SI in the Rosaceae subtribe Pyrinae.

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.

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.

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.

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

It has been argued that the common ancestor of about 75% of all dicots possessed an S-RNase-based gametophytic self-incompatibility (GSI) system. S-RNase genes should thus be found in most plant families showing GSI. The S-RNase gene (or a duplicate) may also acquire a new function and thus genes belonging to the S-RNase lineage may also persist in plant families without GSI. Nevertheless, sequences that belong to the S-RNase lineage have been found in the Solanaceae, Scrophulariaceae, Rosaceae, Cucurbitaceae, and Fabaceae plant families only. Here we search for new sequences that may belong to the S-RNase lineage, using both a phylogenetic and a much faster and simpler amino acid pattern-based approach. We show that the two methods have an apparently similar false-negative rate of discovery (approximately 10%). The amino acid pattern-based approach produces about 15% false positives. Genes belonging to the S-RNase lineage are found in three new plant families, namely, the Rubiaceae, Euphorbiaceae, and Malvaceae. Acquisition of a new function by genes belonging to the S-RNase lineage is shown to be a frequent event. A putative S-RNase sequence is identified in Lotus, a plant genus for which molecular studies on GSI are lacking. The hypothesis of a single origin for S-RNase-based GSI (before the split of the Asteridae and Rosidae) is further supported by the finding of genes belonging to the S-RNase lineage in some of the oldest lineages of the Asteridae and Rosidae, and by Baysean constrained tree analyses.

Purification and characterization of a non-S-RNase and S-RNases from styles of Japanese pear (Pyrus pyrifolia)

In this study we biochemically characterized stylar ribonucleases (RNases) of Japanese pear (Pyrus pyrifolia), which exhibits S-RNase-based gametophytic self-incompatibility. We separated the RNase fractions NS-1, NS-2, and NS-3 from stylar extracts of the cultivar Nijisseiki (S(2)S(4)). The RNase in each fraction was purified to homogeneity through a series of chromatographic steps. Chemical analysis of the proteins revealed that the basic RNases in the NS-2 and NS-3 fractions were the S(4)- and S(2)-RNases, respectively. Five additional S-RNases were purified from other cultivars. An acidic RNase in the NS-1 fraction was also purified from other cultivars, and identified as a non-S-allele-associated RNase (non-S-RNase). The non-S-RNase is composed of 203 amino acids, is non-glycosylated and is a N-terminal-pyroglutamylated enzyme of the RNase T(2) family. The substrate specificities and optimum pH levels of the non-S-RNase and S-RNases were similar. Interestingly, the specific activity of the non-S-RNase was 7.5-221-fold higher than those of the S-RNases when tolura yeast RNA was used as the substrate. The specific activity of the S(2)-RNase was 8.8-28.6-fold lower than those of the other S-RNases. These differences in specific activities among the stylar RNases are discussed.

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.

Origin and dissemination of the pollen-part mutated SC haplotype which confers self-compatibility in apricot (Prunus armeniaca)

In China, its centre of origin, apricot (Prunus armeniaca) is self-incompatible. However, most European cultivars are self-compatible. In most cases, self-compatibility is a result of a loss-of-function mutation within the pollen gene (SFB) in the SC haplotype. Controlled pollinations performed in this work revealed that the cross 'Ceglédi óriás' (S8S9)x'Ceglédi arany' (SCS9) set well, as expected, but the reciprocal cross did not. Apricot S8, S9 and SC haplotypes were analysed using a multilevel approach including fruit set evaluation, pollen tube growth analysis, RNase activity assays, polymerase chain reaction (PCR) analysis and DNA sequencing of the S-RNase and SFB alleles. SFB8 was revealed to be the first known progenitor allele of a naturally occurring self-compatibility allele in Prunus, and consequently SC=The first intron of SC-RNase is a phase one intron, indicating its more recent evolutionary origin compared with the second intron. Sequence analysis of different cultivars revealed that more single nucleotide polymorphisms accumulated in SC-RNase than in SFBC. New methods were designed to allow high-throughput analysis of S genotypes of apricot cultivars and selections. S-RNase sequence data from various sources helped to elucidate the putative origin and dissemination of self-compatibility in apricot conferred by the SC haplotype.

Self-compatible peach (Prunus persica) has mutant versions of the S haplotypes found in self-incompatible Prunus species

This study demonstrates that self-compatible (SC) peach has mutant versions of S haplotypes that are present in self-incompatible (SI) Prunus species. All three peach S haplotypes, S (1), S (2), and S (2m), found in this study encode mutated pollen determinants, SFB, while only S (2m) has a mutation that affects the function of the pistil determinant S-RNase. A cysteine residue in the C5 domain of the S (2m)-RNase is substituted by a tyrosine residue, thereby reducing RNase stability. The peach SFB mutations are similar to the SFB mutations found in SC haplotypes of sweet cherry (P. avium) and Japanese apricot (P. mume). SFB (1) of the S (1) haplotype, a mutant version of almond (P. dulcis) S (k) haplotype, encodes truncated SFB due to a 155 bp insertion. SFB (2) of the S (2) and S (2m) haplotypes, both of which are mutant versions of the S (a) haplotype in Japanese plum (P. salicina), encodes a truncated SFB due to a 5 bp insertion. Thus, regardless of the functionality of the pistil determinant, all three peach S haplotypes are SC haplotypes. Our finding that peach has mutant versions of S haplotypes that function in almond and Japanese plum, which are phylogenetically close and remote species, respectively, to peach in the subfamily Prunoideae of the Roasaceae, provides insight into the SC/SI evolution in Prunus. We discuss the significance of SC pollen part mutation in peach with special reference to possible differences in the SI mechanisms between Prunus and Solanaceae.

Inheritance and interactions of incompatibility alleles in the tetraploid sour cherry

Three progenies of sour cherry (Prunus cerasus) were analysed to correlate self-(in)compatibility status with S-RNase phenotype in this allotetraploid hybrid of sweet and ground cherry. Self-(in)compatibility was assessed in the field and by monitoring pollen tube growth after selfing. The S-RNase phenotypes were determined by isoelectric focusing of stylar proteins and staining for RNase activity and, for the parents, confirmed by PCR. Seedling phenotypes were generally consistent with disomic segregation of S-RNase alleles. The genetic arrangements of the parents were deduced to be 'Köröser' (self-incompatible) S1S4.S(B) S(D), 'Schattenmorelle' (self-compatible) S6S13.S(B)S(B), and clone 43.87 (self-compatible) S4S13.S(B)S(B), where "." separates the two homologous genomes. The presence of S4 and S6 alleles at the same locus led to self-incompatibility, whereas S13 and S(B) at homologous loci led to self-compatibility. The failure of certain heteroallelic genotypes in the three crosses or in the self-incompatible seedlings indicates that S4 and S6 are dominant to S(B). However, the success of S13S(B) pollen on styles expressing corresponding S-RNases indicates competitive interaction or lack of pollen-S components. In general, the universal compatibility of S13S(B) pollen may explain the frequent occurrence of S13 and S(B) together in sour cherry cultivars. Alleles S(B) and S(D), that are presumed to derive from ground cherry, and S13, presumably from sweet cherry, were sequenced. Our findings contribute to an understanding of inheritance of self-(in)compatibility, facilitate screening of progenies for self-compatibility and provide a basis for studying molecular interactions in heteroallelic pollen.