S locus F-box brothers: multiple and pollen-specific F-box genes with S haplotype-specific polymorphisms in apple and Japanese pear

An S-locus F-box protein as pollen S determinant targets non-self S-RNase underlying self-incompatibility in Citrus

Self-incompatibility (SI) is a crucial mechanism that prevents self-fertilization and inbreeding in flowering plants. Citrus exhibits SI regulated by a polymorphic S-locus containing an S-RNase gene and multiple S-locus F-box (SLF) genes. It has been documented that S-RNase functions as the pistil S determinant, but there is no direct evidence that the SLF genes closely linked with S-RNase function as pollen S determinants in Citrus. This study assembled the genomes of two pummelo (Citrus grandis) plants, obtained three novel complete and well-annotated S-haplotypes, and isolated 36 SLF or SLF-like alleles on the S-loci. Phylogenetic analysis of 138 SLFs revealed that the SLF genes were classified into 12 types, including six types with divergent or missing alleles. Furthermore, transformation experiments verified that the conserved S6-SLF7a protein can lead to the transition of SI to self-compatibility by recognizing non-self S8-RNase in 'Mini-Citrus' plants (S7S8 and S8S29, Fortunella hindsii), a model plant for citrus gene function studies. In vitro assays demonstrated interactions between SLFs of different S haplotypes and the Skp1-Cullin1-F-box subunit CgSSK1 protein. This study provides direct evidence that SLF controls the pollen function in Citrus, demonstrating its role in the 'non-self recognition' SI system.

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

The influence of the pollination compatibility type on the pistil S-RNase expression in European pear (Pyrus communis)

The S-RNase gene plays an essential role in the gametophytic self-incompatibility (GSI) system of Pyrus. It codes for the stylar-expressed S-RNase protein which inhibits the growth of incompatible pollen tubes through cytotoxicity and the induction of programmed cell death in the pollen tube. While research on the Pyrus GSI system has primarily focused on the S-RNase gene, there is still a lack of insight into its spatiotemporal expression profile and the factors that regulate it. Previous studies have suggested that S-RNase expression in the style is influenced by pollination and is dependent on the compatibility type. We here continue on this basic hypothesis by analyzing the spatiotemporal expression of the S-RNase alleles in Pyrus communis "Conference" styles in response to different types of pollination; namely, upon full- and semi-compatible pollination and upon incompatible selfing. The results revealed that temporal dynamics of S-RNase expression are influenced by the pollen's compatibility type, indicating the presence of a signaling mechanism between pollen and style to control S-RNase production during pollen tube growth. In our experiment, S-RNase expression continuously decreased after cross-pollination and in the unpollinated control. However, after a fully incompatible pollination, S-RNase expression remained constant. Finally, semi-compatible pollination showed a initially constant S-RNase expression for both alleles followed by a strong decrease in expression. Based on these results and previous findings, we propose a regulatory mechanism to explain the effect of pollination and the associated compatibility type on S-RNase expression in the style. This proposed mechanism could be used as a starting point for future research.

Transcriptomic Analysis of Self-Incompatibility in Alfalfa

Alfalfa (Medicago sativa L.) is an important forage crop worldwide, but molecular genetics and breeding research in this species are hindered by its self-incompatibility (SI). Although the mechanisms underlying SI have been extensively studied in other plant families, SI in legumes, including alfalfa, remains poorly understood. Here, we determined that self-pollinated pollen tubes could germinate on the stigma of alfalfa, grow through the style, and reach the ovarian cavity, but the ovules collapsed ~48 h after self-pollination. A transcriptomic analysis of dissected pistils 24 h after self-pollination identified 941 differently expressed genes (DEGs), including 784 upregulated and 157 downregulated genes. A gene ontology (GO) analysis showed that the DEGs were highly enriched in functions associated with the regulation of pollen tube growth and pollen germination. A Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that pentose and glucuronate interconversion, plant hormone signal transduction, the spliceosome, and ribosomes might play important roles in SI. Our co-expression analysis showed that F-box proteins, serine/threonine protein kinases, calcium-dependent protein kinases (CDPKs), bHLHs, bZIPs, and MYB-related family proteins were likely involved in the SI response. Our study provides a catalog of candidate genes for further study to understand SI in alfalfa and related legumes.

Comparative physiochemical and transcriptomic analysis reveals the influences of cross-pollination on ovary and fruit development in pummelo (Citrus maxima)

'Shuijingmiyou' pummelo (SJ), one of the most popular fruits in Yunnan province of China, is of relatively low fruit shape (FS) quality. In this study, we compared the FS promoting effects of cross pollinations using pollens from seven pummelo varieties, and found that 'Guanximiyou' pummelo (GX) cross-pollination showed the best FS promoting effects on SJ fruits by shortening its fruit neck. To explore the underlying mechanism, physiochemical and transcriptomic differences between self- and cross-pollinated SJ ovaries (SJO and GXO) were investigated. Higher salicylic acid, gibberellin and indole acetic acid contents and superoxide dismutase, peroxidase and catalase activities, and lower polyphenol oxidase activity were determined in GXO compared with SJO. Enrichment analysis of the identified 578 differentially expressed genes (123 up-regulated and 455 down-regulated) in GXO showed that genes involved in solute transport, RNA biosynthesis, phytohormone action and cell wall organization were significantly enriched. The results obtained in this study will be helpful in understanding the influences of cross-pollination on pummelo ovary and fruit development, and can provide the basis for clarifying the underlying mechanism of cross-pollination improved fruit quality.

HT-B and S-RNase CRISPR-Cas9 double knockouts show enhanced self-fertility in diploid Solanum tuberosum

The Gametophytic Self-Incompatibility (GSI) system in diploid potato (Solanum tuberosum L.) poses a substantial barrier in diploid potato breeding by hindering the generation of inbred lines. One solution is gene editing to generate self-compatible diploid potatoes which will allow for the generation of elite inbred lines with fixed favorable alleles and heterotic potential. The S-RNase and HT genes have been shown previously to contribute to GSI in the Solanaceae family and self-compatible S. tuberosum lines have been generated by knocking out S-RNase gene with CRISPR-Cas9 gene editing. This study employed CRISPR-Cas9 to knockout HT-B either individually or in concert with S-RNase in the diploid self-incompatible S. tuberosum clone DRH-195. Using mature seed formation from self-pollinated fruit as the defining characteristic of self-compatibility, HT-B-only knockouts produced little or no seed. In contrast, double knockout lines of HT-B and S-RNase displayed levels of seed production that were up to three times higher than observed in the S-RNase-only knockout, indicating a synergistic effect between HT-B and S-RNase in self-compatibility in diploid potato. This contrasts with compatible cross-pollinations, where S-RNase and HT-B did not have a significant effect on seed set. Contradictory to the traditional GSI model, self-incompatible lines displayed pollen tube growth reaching the ovary, yet ovules failed to develop into seeds indicating a potential late-acting self-incompatibility in DRH-195. Germplasm generated from this study will serve as a valuable resource for diploid potato breeding.

Integration of semi-in vivo assays and multi-omics data reveals the effect of galloylated catechins on self-pollen tube inhibition in Camellia oleifera

Camellia oil extracted from the seeds of Camellia oleifera Abel. is a popular and high-quality edible oil, but its yield is limited by seed setting, which is mainly caused by self-incompatibility (SI). One of the obvious biological features of SI plants is the inhibition of self-pollen tubes; however, the underlying mechanism of this inhibition in C. oleifera is poorly understood. In this study, we constructed a semi-in vivo pollen tube growth test (SIV-PGT) system that can screen for substances that inhibit self-pollen tubes without interference from the genetic background. Combined with multi-omics analysis, the results revealed the important role of galloylated catechins in self-pollen tube inhibition, and a possible molecular regulatory network mediated by UDP-glycosyltransferase (UGT) and serine carboxypeptidase-like (SCPL) was proposed. In summary, galloylation of catechins and high levels of galloylated catechins are specifically involved in pollen tube inhibition under self-pollination rather than cross-pollination, which provides a new understanding of SI in C. oleifera. These results will contribute to sexual reproduction research on C. oleifera and provide theoretical support for improving Camellia oil yield in production.

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

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

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

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

A Dual-Successive-Screen Model at Pollen/Stigma and Pollen Tube/Ovary Explaining Paradoxical Self-Incompatibility Diagnosis in the Olive Tree—An Interpretative Update of the Literature

The ‘pollen test’ and ‘fruit set test’ following controlled crossing combinations of parents are the most commonly used methods for pollination incompatibility studies in Olea europaea L. Self-incompatibility (SI), with diagnoses based on the pollen test and pollen germination, indicating self-compatibility, is not always followed by fruit set in this species. To solve this dispute, we have reconciled all observations into a new model. Mismatches between field and laboratory data and between methods are resolved by the dual-successive-screen model (DSSM) supposing two different loci for the expression of the two SI mechanisms. Pollen/stigma is controlled by diallelic SI, or DSI, inferring two G1 and G2 compatibility/incompatibility (C/I) groups for varieties, whereas pollen tubes in ovaries are controlled by poly-allelic SI or PASI with twenty C/I groups. To explain the selfing of varieties, we have suggested that some determinants in the pollen tube and stigma are unstable and degrade (DS-D for degradation of S-determinant) after three to five days, enabling some pollen tubes to avoid being rejected, hence reaching ovules. DSI and PASI in the DSSM and DS-D mechanisms, plus the andromonoecy of the olive tree, complexify SI studies. Inferences from DSSM and DS-D mechanisms in olive orchard practice are detailed.

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.

Deciphering S-RNase Allele Patterns in Cultivated and Wild Accessions of Italian Pear Germplasm

The genus Pyrus is characterized by an S-RNase-based gametophytic self-incompatibility (GSI) system, a mechanism that promotes outbreeding and prevents self-fertilization. While the S-genotype of the most widely known pear cultivars was already described, little is known on the S-allele variability within local accessions. The study was conducted on 86 accessions encompassing most of the local Sicilian varieties selected for their traits of agronomic interest and complemented with some accessions of related wild species (P. pyrifolia Nakai, P. amygdaliformis Vill.) and some national and international cultivars used as references. The employment of consensus and specific primers enabled the detection of 24 S-alleles combined in 48 S-genotypes. Results shed light on the distribution of the S-alleles among accessions, with wild species and international cultivars characterized by a high diversity and local accessions showing a more heterogeneous distribution of the S-alleles, likely reflecting a more complex history of hybridization. The S-allele distribution was largely in agreement with the genetic structure of the studied collection. In particular, the “wild” genetic background was often characterized by the same S-alleles detected in P. pyrifolia and P. amygdaliformis. The analysis of the S-allele distribution provided novel insight into the contribution of the wild and international cultivars to the genetic background of the local Sicilian or national accessions. Furthermore, these results provide information that can be readily employed by breeders for the set-up of novel mating schemes.

Genome-Wide Identification and Transcriptional Expression of the METTL21C Gene Family in Chicken

The chicken is a common type of poultry that is economically important both for its medicinal and nutritional values. Previous studies have found that free-range chickens have more skeletal muscle mass. The methyltransferase-like 21C gene (METTL21C) plays an important role in muscle development; however, there have been few reports on the role of METTL21C in chickens. In this study, we performed a genome-wide identification of chicken METTL21C genes and analyzed their phylogeny, transcriptional expression profile, and real-time quantitative polymerase chain reaction (qPCR). We identified 10 GgMETTL21C genes from chickens, 11 from mice, and 32 from humans, and these genes were divided into six groups, which showed a large amount of variation among these three species. A total of 15 motifs were detected in METTL21C genes, and the intron phase of the gene structure showed that the METTL21C gene family was conservative in evolution. Further, both the transcript data and qPCR showed that a single gene’s (GgMETTL21C3) expression level increased with the muscle development of chickens, indicating that the METTL21C genes are involved in the development of chicken muscles. Our results provide some reference value for the subsequent study of the function of METTL21C.

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

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

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

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

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.

Platanus-allee is a de novo haplotype assembler enabling a comprehensive access to divergent heterozygous regions

The ultimate goal for diploid genome determination is to completely decode homologous chromosomes independently, and several phasing programs from consensus sequences have been developed. These methods work well for lowly heterozygous genomes, but the manifold species have high heterozygosity. Additionally, there are highly divergent regions (HDRs), where the haplotype sequences differ considerably. Because HDRs are likely to direct various interesting biological phenomena, many genomic analysis targets fall within these regions. However, they cannot be accessed by existing phasing methods, and we have to adopt costly traditional methods. Here, we develop a de novo haplotype assembler, Platanus-allee ( http://platanus.bio.titech.ac.jp/platanus2 ), which initially constructs each haplotype sequence and then untangles the assembly graphs utilizing sequence links and synteny information. A comprehensive benchmark analysis reveals that Platanus-allee exhibits high recall and precision, particularly for HDRs. Using this approach, previously unknown HDRs are detected in the human genome, which may uncover novel aspects of genome variability.

Genome-Wide Identification and Transcriptional Expression Profiles of the F-box Gene Family in Common Walnut (Juglans regia L.)

The common walnut (or Persian walnut), Juglans regia L., is an economically important temperate tree species valued for both its edible nut and high-quality wood. F-box gene family members are involved in plant development, which includes regulating plant development, reproduction, cellular protein degradation, response to biotic and abiotic stresses, and flowering. However, in common walnut (J. regia), there are no reports about the F-box gene family. Here, we report a genome-wide identification of J. regia F-box genes and analyze their phylogeny, duplication, microRNA, pathway, and transcriptional expression profile. In this study, 74 F-box genes were identified and clustered into three groups based on phylogenetic analysis and eight subfamilies based on special domains in common walnut. These common walnut F-box genes are distributed on 31 different pseudo-chromosomes. The gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and microRNA profiles showed that the F-box gene family might play a critical role in the flowering of common walnut. The expressions were significantly higher in female flowers and male flowers compared with leaf and hull tissues at a transcriptome level. The results revealed that the expressions of the F-box gene in female flowers were positively correlated with male flowers, but there was no correlation between any other tissue combinations in common walnut. Our results provided insight into the general characteristics of the F-box genes in common walnut.

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.

Overcoming Self-Incompatibility in Diploid Potato Using CRISPR-Cas9

Potato breeding can be redirected to a diploid inbred/F1 hybrid variety breeding strategy if self-compatibility can be introduced into diploid germplasm. However, the majority of diploid potato clones (Solanum spp.) possess gametophytic self-incompatibility that is primarily controlled by a single multiallelic locus called the S-locus which is composed of tightly linked genes, S-RNase (S-locus RNase) and multiple SLFs (S-locus F-box proteins), which are expressed in the style and pollen, respectively. Using S-RNase genes known to function in the Solanaceae gametophytic SI mechanism, we identified S-RNase alleles with flower-specific expression in two diploid self-incompatible potato lines using genome resequencing data. Consistent with the location of the S-locus in potato, we genetically mapped the S-RNase gene using a segregating population to a region of low recombination within the pericentromere of chromosome 1. To generate self-compatible diploid potato lines, a dual single-guide RNA (sgRNA) strategy was used to target conserved exonic regions of the S-RNase gene and generate targeted knockouts (KOs) using a Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9 (Cas9) approach. Self-compatibility was achieved in nine S-RNase KO T₀ lines which contained bi-allelic and homozygous deletions/insertions in both genotypes, transmitting self compatibility to T₁ progeny. This study demonstrates an efficient approach to achieve stable, consistent self-compatibility through S-RNase KO for use in diploid potato breeding approaches.

Characterization of 25 full-length S-RNase alleles, including flanking regions, from a pool of resequenced apple cultivars

Data obtained from Illumina resequencing of 63 apple cultivars were used to obtain full-length S-RNase sequences using a strategy based on both alignment and de novo assembly of reads. The reproductive biology of apple is regulated by the S-RNase-based gametophytic self-incompatibility system, that is genetically controlled by the single, multi-genic and multi-allelic S locus. Resequencing of apple cultivars provided a huge amount of genetic data, that can be aligned to the reference genome in order to characterize variation to a genome-wide level. However, this approach is not immediately adaptable to the S-locus, due to some peculiar features such as the high degree of polymorphism, lack of colinearity between haplotypes and extensive presence of repetitive elements. In this study we describe a dedicated procedure aimed at characterizing S-RNase alleles from resequenced cultivars. The S-genotype of 63 apple accessions is reported; the full length coding sequence was determined for the 25 S-RNase alleles present in the 63 resequenced cultivars; these included 10 previously incomplete sequences (S ₅ , S _6a , S _6b , S ₈ , S ₁₁ , S ₂₃ , S ₃₉ , S ₄₆ , S ₅₀ and S ₅₈ ). Moreover, sequence divergence clearly suggests that alleles S _6a and S _6b , proposed to be neutral variants of the same alleles, should be instead considered different specificities. The promoter sequences have also been analyzed, highlighting regions of homology conserved among all the alleles.

CRISPR/Cas9-mediated knockout of PiSSK1 reveals essential role of S-locus F-box protein-containing SCF complexes in recognition of non-self S-RNases during cross-compatible pollination in self-incompatible Petunia inflata

Function of Petunia PiSSK1. Self-incompatibility (SI), an inbreeding-preventing mechanism, is regulated in Petunia inflata by the polymorphic S-locus, which houses multiple pollen-specific S-locus F-box (SLF) genes and a single pistil-specific S-RNase gene. S ₂-haplotype and S ₃-haplotype possess the same 17 polymorphic SLF genes (named SLF1 to SLF17), and each SLF protein produced in pollen is assembled into an SCF (Skp1-Cullin1-F-box) E3 ubiquitin ligase complex. A complete suite of SLF proteins is thought to collectively interact with all non-self S-RNases to mediate their ubiquitination and degradation by the 26S proteasome, allowing cross-compatible pollination. For each SCF^SLF complex, the Cullin1 subunit (named PiCUL1-P) and Skp1 subunit (named PiSSK1), like the F-box protein subunits (SLFs), are pollen-specific, raising the possibility that they also evolved specifically to function in SI. Here we used CRISPR/Cas9-meditated genome editing to generate frame-shift indel mutations in PiSSK1 and examined the SI behavior of a T ₀ plant (S ₂ S ₃) with biallelic mutations in the pollen genome and two progeny plants (S ₂ S ₂) each homozygous for one of the indel alleles and not carrying the Cas9-containing T-DNA. Their pollen was completely incompatible with pistils of seven otherwise-compatible S-genotypes, but fully compatible with pistils of an S ₃ S ₃ transgenic plant in which production of S₃-RNase was completely suppressed by an antisense S ₃-RNase gene, and with pistils of immature flower buds, which produce little S-RNase. These results suggest that PiSSK1 specifically functions in SI and support the hypothesis that SLF-containing SCF complexes are essential for compatible pollination.

Reproduction in woody perennial Citrus: an update on nucellar embryony and self-incompatibility

Review on citrus reproduction. Citrus is one of the most important and widely grown fruit crops. It possesses several special reproductive characteristics, such as nucellar embryony and self-incompatibility. The special phenomenon of nucellar embryony in citrus, also known as the polyembryony, is a kind of sporophytic apomixis. During the past decade, the emergence of novel technologies and the construction of multiple citrus reference genomes have facilitated rapid advances to our understanding of nucellar embryony. Indeed, several research teams have preliminarily determined the genetic basis of citrus apomixis. On the other hand, the phenomenon of self-incompatibility that promotes genetic diversity by rejecting self-pollen and accepting non-self-pollen is difficult to study in citrus because the long juvenile period of citrus presents challenges to identifying candidate genes that control this phenomenon. In this review, we focus on advances to our understanding of reproduction in citrus from the last decade and discuss priorities for the coming decade.

SLFL Genes Participate in the Ubiquitination and Degradation Reaction of S-RNase in Self-compatible Peach

It has been proved that the gametophytic self-incompatibility (GSI), mainly exists in Rosaceae and Solanaceae, is controlled by S genes, which are two tightly linked genes located at highly polymorphic S-locus: the S-RNase for pistil specificity and the F-box gene (SFB/SLF) for pollen specificity, respectively. However, the roles of those genes in SI of peach are still a subject of extensive debate. In our study, we selected 37 representative varieties according to the evolution route of peach and identified their S genotypes. We cloned pollen determinant genes mutated PperSFB1m, PperSFB2m, PperSFB4m, and normal PperSFB2, and style determinant genes PperS1-RNase, PperS2-RNase, PperS2m-RNase, and PperS4-RNase. The mutated PperSFBs encode truncated SFB proteins due to a fragment insertion. The truncated PperSFBs and normal PperSFB2 interacted with PperS-RNases demonstrated by Y2H. Normal PperSFB2 was divided into four parts: box, box-V1, V1-V2, and HVa-HVb. The box domain of PperSFB2 did not interact with PperS-RNases, both of the box-V1 and V1-V2 had interactions with PperS-RNases, while the hypervariable region of PperSFB2 HVa-HVb only interacted with PperS2-RNase showed by Y2H and BiFC assay. Bioinformatics analysis of peach genome revealed that there were other F-box genes located at S-locus, and of which three F-box genes were specifically expressed in pollen, named as PperSLFL1, PperSLFL2, and PperSLFL3, respectively. In phylogenetic analysis PperSLFLs clustered with Maloideae SFBB genes, and PperSFB genes were clustered into the other group with other SFB genes of Prunus. Protein interaction analysis revealed that the three PperSLFLs interacted with PperSSK1 and PperS-RNases with no allelic specificity. In vitro ubiquitination assay showed that PperSLFLs could tag ubiquitin molecules onto PperS-RNases. The above results suggest that three PperSLFLs are the appropriate candidates for the "general inhibitor," which would inactivate the S-RNases in pollen tubes, involved in the self-incompatibility of peach.

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.

Electrostatic potentials of the S-locus F-box proteins contribute to the pollen S specificity in self-incompatibility in Petunia hybrida

Self-incompatibility (SI) is a self/non-self discrimination system found widely in angiosperms and, in many species, is controlled by a single polymorphic S-locus. In the Solanaceae, Rosaceae and Plantaginaceae, the S-locus encodes a single S-RNase and a cluster of S-locus F-box (SLF) proteins to control the pistil and pollen expression of SI, respectively. Previous studies have shown that their cytosolic interactions determine their recognition specificity, but the physical force between their interactions remains unclear. In this study, we show that the electrostatic potentials of SLF contribute to the pollen S specificity through a physical mechanism of 'like charges repel and unlike charges attract' between SLFs and S-RNases in Petunia hybrida. Strikingly, the alteration of a single C-terminal amino acid of SLF reversed its surface electrostatic potentials and subsequently the pollen S specificity. Collectively, our results reveal that the electrostatic potentials act as a major physical force between cytosolic SLFs and S-RNases, providing a mechanistic insight into the self/non-self discrimination between cytosolic proteins in angiosperms.

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

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

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.

A simple, high-throughput modeling approach reveals insights into the mechanism of gametophytic self-incompatibility

Specificity in the GSI response results from the S-haplotype-specific molecular interaction of S-locus F-box (SLF/SFB) and SRNase proteins in the self-incompatibility locus (S-locus). The answer to the question of how these two components of the S-locus (SRNase and SLF/SFB) interact has been gathered from several models. Since there is not enough evidence as to which one is the definitive model, none of them can be ruled out. Despite the identification of interacting protein elements, the mechanism by which SLF/SFB and SRNase interact to differently trigger the self-incompatibility among families and subfamilies remain uncertain. The high-throughput modeling approach demonstrates structural visions into the possible existence of a Collaborative Non-Self Recognition model in apple. These findings postulate several prospects for future investigation providing useful information to guide the implementation of breeding strategies.

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

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

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

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

Dissect style response to pollination using metabolite profiling in self-compatible and self-incompatible tomato species

Tomato style is the pathway for pollen germination and pollen tubes growth from the stigma to the ovules where fertilization occurs. It is essential to supplying the nutrients for pollen tube growth and guidance for the pollen tubes. To our knowledge, style also regulates gametophytic self-incompatibility (SI) in tomato species. This study identified the metabolites and monitored the metabolic changes of self-incompatible and self-compatible tomato with self-pollinated or unpollinated styles by gas chromatography-mass spectrometry (GC-MS). A total of 9 classes of compounds were identified in SI and self-compatibility (SC) self-pollinated and unpollinated styles which included amino acids, sugars, fatty acids/lipids, amines, organic acids, alcohols, nitriles, inorganic acids and other compounds. The contents of d-Mannose-6-phosphate, Cellobiose, Myristic acid, 2,4-Diaminobutyric acid, Inositol and Urea were significantly decreased and the rest did not significantly change in SI styles. But change of metabolites content significantly happened in SC styles. In addition, among the total 9 classes of compounds, the different metabolites accounted for a different proportion in amino acids, sugars, amines, organic acids and alcohols compared SC and SI. The result indicated that the physiological changes of styles existed differences in SC and SI after self pollination.

Molecular mechanism of the S-RNase-based gametophytic self-incompatibility in fruit trees of Rosaceae

Self-incompatibility (SI) is a major obstacle for stable fruit production in fruit trees of Rosaceae. SI of Rosaceae is controlled by the S locus on which at least two genes, pistil S and pollen S, are located. The product of the pistil S gene is a polymorphic and extracellular ribonuclease, called S-RNase, while that of the pollen S gene is a protein containing the F-box motif, SFB (S haplotype-specific F-box protein)/SFBB (S locus F-box brothers). Recent studies suggested that SI of Rosaceae includes two different systems, i.e., Prunus of tribe Amygdaleae exhibits a self-recognition system in which its SFB recognizes self-S-RNase, while tribe Pyreae (Pyrus and Malus) shows a non-self-recognition system in which many SFBB proteins are involved in SI, each recognizing subset of non-self-S-RNases. Further biochemical and biological characterization of the S locus genes, as well as other genes required for SI not located at the S locus, will help our understanding of the molecular mechanisms, origin, and evolution of SI of Rosaceae, and may provide the basis for breeding of self-compatible fruit tree cultivars.

Advances in Japanese pear breeding in Japan

The Japanese pear (Pyrus pyrifolia Nakai) is one of the most widely grown fruit trees in Japan, and it has been used throughout Japan's history. The commercial production of pears increased rapidly with the successive discoveries of the chance seedling cultivars 'Chojuro' and 'Nijisseiki' around 1890, and the development of new cultivars has continued since 1915. The late-maturing, leading cultivars 'Niitaka' and 'Shinko' were released during the initial breeding stage. Furthermore, systematic breeding by the Horticultural Research Station (currently, NARO Institute of Fruit Tree Science, National Agriculture and Food Research Organization (NIFTS)) began in 1935, which mainly aimed to improve fruit quality by focusing on flesh texture and black spot disease resistance. To date, 22 cultivars have been released, including 'Kosui', 'Hosui', and 'Akizuki', which are current leading cultivars from the breeding program. Four induced mutant cultivars induced by gamma irradiation, which exhibit some resistance to black spot disease, were released from the Institute of Radiation Breeding. Among these cultivars, 'Gold Nijisseiki' has become a leading cultivar. Moreover, 'Nansui' from the Nagano prefectural institute breeding program was released, and it has also become a leading cultivar. Current breeding objectives at NIFTS mainly combine superior fruit quality with traits related to labor and cost reduction, multiple disease resistance, or self-compatibility. Regarding future breeding, marker-assisted selection for each trait, QTL analyses, genome-wide association studies, and genomic selection analyses are currently in progress.

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.

S-genotype identification based on allele-specific PCR in Japanese pear

Gametophytic self-incompatibility in Japanese pear (Pyrus pyrifolia Nakai) is controlled by the single, multi-allelic S-locus. Information about the S-genotypes is important for breeding and the selection of pollen donors for fruit production. Rapid and reliable S-genotype identification system is necessary for efficient breeding of new cultivars in Japanese pear. We designed S allele-specific PCR primer pairs for ten previously reported S-RNase alleles (S (1)-S (9) and S (k)) as simple and reliable method. Specific nucleotide sequences were chosen to design the primers to amplify fragments of only the corresponding S alleles. The developed primer pairs were evaluated by using homozygous S-genotypes (S (1)/S (1)-S (9)/S (9) and S (4sm)/S (4sm)) and 14 major Japanese pear cultivars, and found that S allele-specific primer pairs can identify S-genotypes effectively. The S allele-specific primer pairs developed in this study will be useful for efficient S-genotyping and for marker-assisted selection in Japanese pear breeding programs.

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.

Gene duplication and genetic exchange drive the evolution of S-RNase-based self-incompatibility in Petunia

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16-20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16-20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations.

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.

Identification of SFBB-containing canonical and noncanonical SCF complexes in pollen of apple (Malus × domestica)

Gametophytic self-incompatibility (GSI) of Rosaceae, Solanaceae and Plantaginaceae is controlled by a single polymorphic S locus. The S locus contains at least two genes, S-RNase and F-box protein encoding gene SLF/SFB/SFBB that control pistil and pollen specificity, respectively. Generally, the F-box protein forms an E3 ligase complex, SCF complex with Skp1, Cullin1 (CUL1) and Rbx1, however, in Petunia inflata, SBP1 (S-RNase binding protein1) was reported to play the role of Skp1 and Rbx1, and form an SCFSLF-like complex for ubiquitination of non-self S-RNases. On the other hand, in Petunia hybrida and Petunia inflata of Solanaceae, Prunus avium and Pyrus bretschneideri of Rosaceae, SSK1 (SLF-interacting Skp1-like protein1) is considered to form the SCFSLF/SFB complex. Here, we isolated pollen-expressed apple homologs of SSK1 and CUL1, and named MdSSK1, MdCUL1A and MdCUL1B. MdSSK1 was preferentially expressed in pollen, but weakly in other organs analyzed, while, MdCUL1A and MdCUL1B were almost equally expressed in all the organs analyzed. MdSSK1 transcript abundance was significantly (>100 times) higher than that of MdSBP1. In vitro binding assays showed that MdSSK1 and MdSBP1 interacted with MdSFBB1-S9 and MdCUL1, and MdSFBB1-S9 interacted more strongly with MdSSK1 than with MdSBP1. The results suggest that both MdSSK1-containing SCFSFBB1 and MdSBP1-containing SCFSFBB1-like complexes function in pollen of apple, and the former plays a major role.

Sequence characterization and spatiotemporal expression patterns of PbS26-RNase gene in Chinese White Pear (Pyrus bretschneideri)

Many flowering plants exhibit an important intraspecific reproductive barrier phenomenon, that is, self-incompatibility (SI), in which S-RNase genes play a significant role. To clarify the specific function of S-RNase genes in Chinese pears, the full length cDNA of PbS₂₆-RNase was isolated by rapid amplification of cDNA ends (RACE) technology from Chinese white pear (Pyrus bretschneideri) cultivar “Hongpisu.” The cDNA sequence for PbS₂₆-RNase was deposited in GenBank under accession number EU081888. At the amino acid level, the PbS₂₆-RNase displayed the highest similarity (96.9%) with PcSa-RNase of P. communis, and only seven amino acid differences were present in the two S-RNases. Phylogenetic analysis of rosaceous S-RNases indicated that the PbS₂₆-RNase clustered with maloideous S-RNases, forming a subfamily-specific not a species-specific group. The PbS₂₆-RNase gene was specifically expressed in the style but not other tissues/organs. The expression level of the PbS₂₆-RNase gene rapidly increased at bell balloon stage (BBS), and then it dropped after pollination. However, the abundance of the PbS₂₆-RNase gene transcript in the style was greater after cross-pollination than after self-pollination. In addition, a method for rapidly detecting the PbS₂₆-RNase gene was developed via allele-specific primers design. The present study could provide a scientific basis for fully clarifying the mechanism of pear SI at the molecular level.

A segmental duplication encompassing S-haplotype triggers pollen-part self-compatibility in Japanese pear (Pyrus pyrifolia)

Self-compatible mutants of self-incompatible crops have been extensively studied for research and agricultural purposes. Until now, the only known pollen-part self-compatible mutants in Rosaceae subtribe Pyrinae, which contains many important fruit trees, were polyploid. This study revealed that the pollen-part self-compatibility of breeding selection 415-1, a recently discovered mutant of Japanese pear (Pyrus pyrifolia) derived from γ-irradiated pollen, is caused by a duplication of an S-haplotype. In the progeny of 415-1, some plants had three S-haplotypes, two of which were from the pollen parent. Thus, 415-1 was able to produce pollen with two S-haplotypes, even though it was found to be diploid: the relative nuclear DNA content measured by flow cytometry showed no significant difference from that of a diploid cultivar. Inheritance patterns of simple sequence repeat (SSR) alleles in the same linkage group as the S-locus (LG 17) showed that some SSRs closely linked to S-haplotypes were duplicated in progeny containing the duplicated S-haplotype. These results indicate that the pollen-part self-compatibility of 415-1 is not caused by a mutation of pollen S factors in either one of the S-haplotypes, but by a segmental duplication encompassing the S-haplotype. Consequently, 415-1 can produce S-heteroallelic pollen grains that are capable of breaking down self-incompatibility (SI) by competitive interaction between the two different S factors in the pollen grain. 415-1 is the first diploid pollen-part self-compatible mutant with a duplicated S-haplotype to be discovered in the Pyrinae. The fact that 415-1 is not polyploid makes it particularly valuable for further studies of SI mechanisms.

Genetic features of the spontaneous self-compatible mutant, 'Jin Zhui' (Pyrus bretschneideri Rehd.)

‘Jin Zhui’ is a spontaneous self-compatible mutant of ‘Ya Li’ (Pyrus bretschneideri Rehd. S₂₁S₃₄), the latter displaying a typical S-RNase-based gametophytic self-incompatibility (GSI). The pollen-part mutation (PPM) of ‘Jin Zhui’ might be due to a natural mutation in the pollen-S gene (S₃₄ haplotype). However, the molecular mechanisms behind these phenotypic changes are still unclear. In this study, we identified five SLF (S-Locus F-box) genes in ‘Ya Li’, while no nucleotide differences were found in the SLF genes of ‘Jin Zhui’. Further genetic analysis by S-RNase PCR-typing of selfed progeny of ‘Jin Zhui’ and ‘Ya Li’ × ‘Jin Zhui’ progeny showed three progeny classes (S₂₁S₂₁, S₂₁S₃₄ and S₃₄S₃₄) as opposed to the two classes reported previously (S₂₁S₃₄ and S₃₄S₃₄), indicating that the pollen gametes of ‘Jin Zhui’, bearing either the S₂₁- or S₃₄-haplotype, were able to overcome self-incompatibility (SI) barriers. Moreover, no evidence of pollen-S duplication was found. These findings support the hypothesis that loss of function of S-locus unlinked PPM expressed in pollen leads to SI breakdown in ‘Jin Zhui’, rather than natural mutation in the pollen-S gene (S₃₄ haplotype). Furthermore, abnormal meiosis was observed in a number of pollen mother cells (PMCs) in ‘Jin Zhui’, but not in ‘Ya Li’. These and other interesting findings are discussed.

Isolation and characterization of multiple F-box genes linked to the S9- and S10-RNase in apple (Malus × domestica Borkh.)

Using 11 consensus primer pairs designed from S-linked F-box genes of apple and Japanese pear, 10 new F-box genes (MdFBX21 to 30) were isolated from the apple cultivar 'Spartan' (S(9)S(10)). MdFBX21 to 23 and MdFBX24 to 30 were completely linked to the S(9) -RNase and S(10-)RNase, respectively, and showed pollen-specific expression and S-haplotype-specific polymorphisms. Therefore, these 10 F-box genes are good candidates for the pollen determinant of self-incompatibility in apple. Phylogenetic analysis and comparison of deduced amino acid sequences of MdFBX21 to 30 with those of 25 S-linked F-box genes previously isolated from apple showed that a deduced amino acid identity of greater than 88.0 % can be used as the tentative criterion to classify F-box genes into one type. Using this criterion, 31 of 35 F-box genes of apple were classified into 11 types (SFBB1-11). All types included F-box genes derived from S(3-) and S(9-)haplotypes, and seven types included F-box genes derived from S(3-), S(9-), and S(10-)haplotypes. Moreover, comparison of nucleotide sequences of S-RNases and multiple F-box genes among S(3-), S(9-), and S(10-)haplotypes suggested that F-box genes within each type showed high nucleotide identity regardless of the identity of the S-RNase. The large number of F-box genes as candidates for the pollen determinant and the high degree of conservation within each type are consistent with the collaborative non-self-recognition model reported for Petunia. These findings support that the collaborative non-self-recognition system also exists in apple.

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.

Identification of a canonical SCF(SLF) complex involved in S-RNase-based self-incompatibility of Pyrus (Rosaceae)

S-RNase-based self-incompatibility (SI) is an intraspecific reproductive barrier to prevent self-fertilization found in many species of the Solanaceae, Plantaginaceae and Rosaceae. In this system, S-RNase and SLF/SFB (S-locus F-box) genes have been shown to control the pistil and pollen SI specificity, respectively. Recent studies have shown that the SLF functions as a substrate receptor of a SCF (Skp1/Cullin1/F-box)-type E3 ubiquitin ligase complex to target S-RNases in Solanaceae and Plantaginaceae, but its role in Rosaceae remains largely undefined. Here we report the identification of two pollen-specific SLF-interacting Skp1-like (SSK) proteins, PbSSK1 and PbSSK2, in Pyrus bretschneideri from the tribe Pyreae of Rosaceae. Both yeast two-hybrid and pull-down assays demonstrated that they could connect PbSLFs to PbCUL1 to form a putative canonical SCF(SLF) (SSK/CUL1/SLF) complex in Pyrus. Furthermore, pull-down assays showed that the SSK proteins could bind SLF and CUL1 in a cross-species manner between Pyrus and Petunia. Additionally, phylogenetic analysis revealed that the SSK-like proteins from Solanaceae, Plantaginaceae and Rosaceae form a monoclade group, hinting their shared evolutionary origin. Taken together, with the recent identification of a canonical SCF(SFB) complex in Prunus of the tribe Amygdaleae of Rosaceae, our results show that a conserved canonical SCF(SLF/SFB) complex is present in Solanaceae, Plantaginaceae and Rosaceae, implying that S-RNase-based self-incompatibility shares a similar molecular and biochemical mechanism.

Identification of differentially expressed genes in pistils from self-incompatible Citrus reticulata by suppression subtractive hybridization

Self-incompatibility (SI) is one important factor that can result in Citrus seedlessness. However, the molecular mechanism of SI in Citrus is not clear yet. To isolate the pistil's SI-related genes, a suppression subtractive hybridization library was constructed using mature pistils of 'Wuzishatangju' mandarin (SI) as the tester and mature pistils of 'Shatangju' mandarin (self-compatibility, SC) as the driver. 229 differentially expressed cDNA clones from 967 positive clones were sequenced and identified. Differentially expressed ESTs are possibly involved in the SI reaction of 'Wuzishatangju' through a regulating signaling pathway, serine/threonine phosphatase activity, receptor kinase, embryonic development, gibberellin stimulus, or transcription. 11 out of 36 SI candidate genes displayed different expression patterns in various tissues and stages after self- and cross-pollination of 'Wuzishatangju'. The expression of CaBP (WY65), a senescence-protease (WY372), an unknown gene (WY283), and a WRKY (WY17) were up-regulated in the styles of 'Wuzishatangju' while higher expression of WY190 was observed in styles of 'Shatangju'. Highest expression levels of WY65, WY372, an annexin (WY598), the zinc-finger protein (WY376), a C2-protein (WY291), and an unknown gene (WY318) were detected in styles at 3 days after self-pollination of 'Wuzishatangju' while lowest levels were observed in styles at 3 days after cross-pollination of 'Wuzishatangju' × 'Shatangju'. The potential involvement of these genes in the SI reaction is discussed.

Interhaplotypic heterogeneity and heterochromatic features may contribute to recombination suppression at the S locus in apple (Malusxdomestica)

Gametophytic self-incompatibility (GSI) is controlled by a complex S locus containing the pistil determinant S-RNase and pollen determinant SFB/SLF. Tight linkage of the pistil and pollen determinants is necessary to guarantee the self-incompatibility (SI) function. However, multiple probable pollen determinants of apple and Japanese pear, SFBBs (S locus F-box brothers), exist in each S haplotype, and how these multiple genes maintain the SI function remains unclear. It is shown here by high-resolution fluorescence in situ hybridization (FISH) that SFBB genes of the apple S9 haplotype are physically linked to the S9-RNase gene, and the S locus is located in the subtelomeric region. FISH analyses also determined the relative order of SFBB genes and S-RNase in the S9 haplotype, and showed that gene order differs between the S9 and S3 haplotypes. Furthermore, it is shown that the apple S locus is located in a knob-like large heterochromatin block where DNA is highly methylated. It is proposed that interhaplotypic heterogeneity and the heterochromatic nature of the S locus help to suppress recombination at the S locus in apple.