Transcriptome analysis reveals the same 17 S-locus F-box genes in two haplotypes of the self-incompatibility locus of Petunia inflata

Long-read genome sequencing reveals the sequence characteristics of pear self-incompatibility locus

The S-RNase-based self-incompatibility locus (S-locus) in Petunia species contains 16-20 F-box genes, which collaboratively function in the recognition and subsequent degradation of non-self S-RNases, while distinguishing them from self S-RNase. However, the number of S-locus F-box genes (SFBBs) physically interacted with non-self S-RNases remains uncertain in Pyrus species. Utilizing Pacbio long-read sequencing, we successfully assembled the genome of pear cultivar 'Yali' (Pyrus bretschneideri), and identified 19 SFBBs from the Pyrus S17-locus spanning approximately 1.78 Mb. Additionally, we identified 17-21 SFBBs from other Pyrus and Malus S-loci spanning a range of 1.35 to 2.64 Mb. Based on the phylogenetic analysis, it was determined that Pyrus and Malus SFBBs could be classified into 22 groups, denoted as I to XXII. At amino acid level, SFBBs within a given group exhibited average identities ranged from 88.9% to 97.9%. Notably, all 19 SFBBs from the S17-locus co-segregated with S17-RNase, with 18 of them being specifically expressed in pollen. Consequently, these 18 pollen-specifically expressed SFBBs are considered potential candidates for the pollen-S determinant. Intriguingly, out of the 18 pollen-specifically expressed SFBBs, eight demonstrated interactions with at least one non-self S-RNase, while the remaining SFBBs failed to recognize any S-RNase. These findings provide compelling evidence supporting the existence of a collaborative non-self-recognition system governing self-incompatibility in pear species.

Rapid detection of RNase-based self-incompatibility in Lysimachia monelli (Primulaceae)

PREMISE

Primroses famously employ a system that simultaneously expresses distyly and filters out self-pollen. Other species in the Primulaceae family, including Lysimachia monelli (blue pimpernel), also express self-incompatibility (SI), but involving a system with distinct features and an unknown molecular genetic basis.

METHODS

We utilize a candidate-based transcriptome sequencing (RNA-seq) approach, relying on candidate T2/S-RNase Class III and S-linked F-box-motif-containing genes and harnessing the unusual evolutionary and genetic features of SI, to examine whether an RNase-based mechanism underlies SI in L. monelli. We term this approach "SI detection with RNA-seq" (SIDR).

RESULTS

The results of sequencing, crossing, population genetics, and molecular evolutionary features each support a causal association linking the recovered genotypes with SI phenotypes. The finding of RNase-based SI in Primulaceae (Ericales) all but cements the long-held view that this mechanism was present in the ancestral pentapetal eudicot, whose descendants now comprise two-thirds of angiosperms. It also significantly narrows the plausible maximum age for the heterostyly evolution within the family.

CONCLUSIONS

SIDR is powerful, flexible, inexpensive, and most critically enables work in often-neglected species. It may be used with or without candidate genes to close enormous gaps in understanding the genetic basis of SI and the history of breeding system evolution.

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.

The role of promiscuous molecular recognition in the evolution of RNase-based self-incompatibility in plants

How do biological networks evolve and expand? We study these questions in the context of the plant collaborative-non-self recognition self-incompatibility system. Self-incompatibility evolved to avoid self-fertilization among hermaphroditic plants. It relies on specific molecular recognition between highly diverse proteins of two families: female and male determinants, such that the combination of genes an individual possesses determines its mating partners. Though highly polymorphic, previous models struggled to pinpoint the evolutionary trajectories by which new specificities evolved. Here, we construct a novel theoretical framework, that crucially affords interaction promiscuity and multiple distinct partners per protein, as is seen in empirical findings disregarded by previous models. We demonstrate spontaneous self-organization of the population into distinct "classes" with full between-class compatibility and a dynamic long-term balance between class emergence and decay. Our work highlights the importance of molecular recognition promiscuity to network evolvability. Promiscuity was found in additional systems suggesting that our framework could be more broadly applicable.

Mapping quantitative trait loci associated with self-(in)compatibility in goji berries (Lycium barbarum)

BACKGROUND

Goji (Lycium barbarum L.) is a perennial deciduous shrub widely distributed in arid and semiarid regions of Northwest China. It is highly valued for its medicinal and functional properties. Most goji varieties are naturally self-incompatible, posing challenges in breeding and cultivation. Self-incompatibility is a complex genetic trait, with ongoing debates regarding the number of self-incompatible loci. To date, no genetic mappings has been conducted for S loci or other loci related to self-incompatibility in goji.

RESULTS

We used genome resequencing to create a high-resolution map for detecting de novo single-nucleotide polymorphisms (SNP) in goji. We focused on 229 F1 individuals from self-compatible '13-19' and self-incompatible 'new 9' varieties. Subsequently, we conducted a quantitative trait locus (QTL) analysis on traits associated with self-compatibility in goji berries. The genetic map consisted of 249,327 SNPs distributed across 12 linkage groups (LGs), spanning a total distance of 1243.74 cM, with an average interval of 0.002 cM. Phenotypic data related to self-incompatibility, such as average fruit weight, fruit rate, compatibility index, and comparable compatibility index after self-pollination and geitonogamy, were collected for the years 2021-2022, as well as for an extra year representing the mean data from 2021 to 2022 (2021/22). A total of 43 significant QTL, corresponding to multiple traits were identified, accounting for more than 11% of the observed phenotypic variation. Notably, a specific QTL on chromosome 2 consistently appeared across different years, irrespective of the relationship between self-pollination and geitonogamy. Within the localization interval, 1180 genes were annotated, including Lba02g01102 (annotated as an S-RNase gene), which showed pistil-specific expression. Cloning of S-RNase genes revealed that the parents had two different S-RNase alleles, namely S1S11 and S2S8. S-genotype identification of the F1 population indicated segregation of the four S-alleles from the parents in the offspring, with the type of S-RNase gene significantly associated with self-compatibility.

CONCLUSIONS

In summary, our study provides valuable insights into the genetic mechanism underlying self-compatibility in goji berries. This highlights the importance of further positional cloning investigations and emphasizes the importance of integration of marker-assisted selection in goji breeding programs.

The Identification and Analysis of the Self-Incompatibility Pollen Determinant Factor SLF in Lycium barbarum

Self-incompatibility is a widespread genetic mechanism found in flowering plants. It plays a crucial role in preventing inbreeding and promoting outcrossing. The genes that control self-incompatibility in plants are typically determined by the S-locus female determinant factor and the S-locus male determinant factor. In the Solanaceae family, the male determinant factor is often the SLF gene. In this research, we cloned and analyzed 13 S2-LbSLF genes from the L. barbarum genome, which are located on chromosome 2 and close to the physical location of the S-locus female determinant factor S-RNase, covering a region of approximately 90.4 Mb. The amino acid sequence identity of the 13 S2-LbSLFs is 58.46%, and they all possess relatively conserved motifs and typical F-box domains, without introns. A co-linearity analysis revealed that there are no tandemly repeated genes in the S2-LbSLF genes, and that there are two pairs of co-linear genes between S2-LbSLF and the tomato, which also belongs to the Solanaceae family. A phylogenetic analysis indicates that the S2-LbSLF members can be divided into six groups, and it was found that the 13 S2-LbSLFs are clustered with the SLF genes of tobacco and Petunia inflata to varying degrees, potentially serving as pollen determinant factors regulating self-incompatibility in L. barbarum. The results for the gene expression patterns suggest that S2-LbSLF is only expressed in pollen tissue. The results of the yeast two-hybrid assay showed that the C-terminal region of S2-LbSLFs lacking the F-box domain can interact with S-RNase. This study provides theoretical data for further investigation into the functions of S2-LbSLF members, particularly for the identification of pollen determinant factors regulating self-incompatibility in L. barbarum.

Hormonal Signaling during dPCD: Cytokinin as the Determinant of RNase-Based Self-Incompatibility in Solanaceae

Research into molecular mechanisms of self-incompatibility (SI) in plants can be observed in representatives of various families, including Solanaceae. Earlier studies of the mechanisms of S-RNase-based SI in petunia (Petunia hybrida E. Vilm.) demonstrate that programmed cell death (PCD) is an SI factor. These studies suggest that the phytohormon cytokinin (CK) is putative activator of caspase-like proteases (CLPs). In this work, data confirming this hypothesis were obtained in two model objects-petunia and tomato (six Solanaceae representatives). The exogenous zeatin treatment of tomato and petunia stigmas before a compatible pollination activates CLPs in the pollen tubes in vivo, as shown via the intravital imaging of CLP activities. CK at any concentration slows down the germination and growth of petunia and tomato male gametophytes both in vitro and in vivo; shifts the pH of the cytoplasm (PHc) to the acid region, thereby creating the optimal conditions for CLP to function and inhibiting the F-actin formation and/or destructing the cytoskeleton in pollen tubes to point foci during SI-induced PCD; and accumulates in style tissues during SI response. The activity of the ISOPENTENYLTRANSFERASE 5 (IPT5) gene at this moment exceeds its activity in a cross-compatible pollination, and the levels of expression of the CKX1 and CKX2 genes (CK OXIDASE/DEHYDROGENASE) are significantly lower in self-incompatible pollination. All this suggests that CK plays a decisive role in the mechanism underlying SI-induced PCD.

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

PlantTribes2: Tools for comparative gene family analysis in plant genomics

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

Fine-Mapping and Comparative Genomic Analysis Reveal the Gene Composition at the S and Z Self-incompatibility Loci in Grasses

Self-incompatibility (SI) is a genetic mechanism of hermaphroditic plants to prevent inbreeding after self-pollination. Allogamous Poaceae species exhibit a unique gametophytic SI system controlled by two multi-allelic and independent loci, S and Z. Despite intense research efforts in the last decades, the genes that determine the initial recognition mechanism are yet to be identified. Here, we report the fine-mapping of the Z-locus in perennial ryegrass (Lolium perenne L.) and provide evidence that the pollen and stigma components are determined by two genes encoding DUF247 domain proteins (ZDUF247-I and ZDUF247-II) and the gene sZ, respectively. The pollen and stigma determinants are located side-by-side and were genetically linked in 10,245 individuals of two independent mapping populations segregating for Z. Moreover, they exhibited high allelic diversity as well as tissue-specific gene expression, matching the expected characteristics of SI determinants known from other systems. Revisiting the S-locus using the latest high-quality whole-genome assemblies revealed a similar gene composition and structure as found for Z, supporting the hypothesis of a duplicated origin of the two-locus SI system of grasses. Ultimately, comparative genomic analyses across a wide range of self-compatible and self-incompatible Poaceae species revealed that the absence of a functional copy of at least one of the six putative SI determinants is accompanied by a self-compatible phenotype. Our study provides new insights into the origin and evolution of the unique gametophytic SI system in one of the largest and economically most important plant families.

Global transcriptome dissection of pollen-pistil interactions induced self-incompatibility in dragon fruit (Selenicereus spp.)

Self-incompatibility (SI) is a major issue in dragon fruit (Selenicereus spp.) breeding and production. Therefore, a better understanding of the dragon fruit SI mechanism is needed to improve breeding efficiency and ultimate production costs. To reveal the underlying mechanisms of SI in dragon fruit, plant anatomy, de novo RNA sequencing-based transcriptomic analysis, and multiple bioinformatic approaches were used to analyze gene expression in the pistils of the self-pollinated and cross-pollinated dragon fruit flowers at different intervals of time after pollination. Using fluorescence microscopy, we observed that the pollen of 'Hongshuijing', a self-incompatible dragon fruit variety (S. monacanthus), germinated on its own stigma. However, the pollen tube elongation has ceased at 1/2 of the style, confirming that dragon fruit experiences gametophyte self-incompatibility (GSI). We found that the pollen tube elongation in vitro was inhibited by self-style glycoproteins in the SI variety, indicating that glycoproteins were involved in SI. That is to say the female S factor should be homologous of S-RNase or PrsS (P. rhoeas stigma S factor), both of which are glycoproteins and are the female S factors of the two known GSI mechanism respectively. Bioinformatics analyses indicated that among the 43,954 assembled unigenes from pistil, there were six S-RNase genes, while 158 F-box genes were identified from a pollen transcriptomic dataset. There were no P. rhoeas type S genes discovered. Thus, the identified S-RNase and F-box represent the candidate female and male S genes, respectively. Analysis of differentially expressed genes (DEGs) between the self and cross-pollinated pistils at different time intervals led to the identification of 6,353 genes. We then used a weighted gene co-expression network analysis (WGCNA) to find some non-S locus genes in SI responses in dragon fruit. Additionally, 13 transcription factors (TFs) (YABBY4, ANL2, ERF43, ARF2, BLH7, KNAT6, PIF3, two OBF1, two HY5 and two LHY/CCA) were identified to be involved in dragon fruit GSI. Thus, we uncovered candidate S and non-S genes and predicted more SI-related genes for a more detailed investigation of the molecular mechanism of dragon fruit SI. Our findings suggest that dragon fruit possesses a GSI system and involves some unique regulators. This study lays the groundwork for future research into SI mechanisms in dragon fruit and other plant species.

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

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

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

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

Development of first linkage map for Silphium integrifolium (Asteraceae) enables identification of sporophytic self-incompatibility locus

Silphium integrifolium (Asteraceae) has been identified as a candidate for domestication as a perennial oilseed crop and is assumed to have sporophytic self-incompatibility system-the genetic basis of which is not well understood in the Asteraceae. To address this gap, we sought to map the genomic location of the self-recognition locus (S-locus) in this species. We used a biparental population and genotyping-by-sequencing to create the first genetic linkage map for this species, which contained 198 SNP markers and resolved into the correct number of linkage groups. Then we developed a novel crossing scheme and set of analysis methods in order to infer S-locus genotypes for a subset of these individuals, allowing us to map the trait. Finally, we evaluated potential genes of interest using synteny analysis with the annual sunflower (Helianthus annuus) and lettuce (Lactuca sativa) genomes. Our results confirm that S. integrifolium does indeed have a sporophytic self-incompatibility system. Our method is effective and efficient, allowed us to map the S. integrifolium S-locus using fewer resources than existing methods, and could be readily applied to other species.

The proteome reveals the involvement of serine/threonine kinase in the recognition of self- incompatibility in almond

The self-incompatibility recognition mechanism determines whether the gametophyte is successfully fertilized between pollen tube SCF (SKP1-CUL1-F-box-RBX1) protein and pistil S-RNase protein during fertilization is unclear. In this study, the pistils of two almond cultivars 'Wanfeng' and 'Nonpareil' were used as the experimental materials after self- and nonself/cross-pollination, and pistils from the stamen-removed flowers were used as controls. We used fluorescence microscopy to observe the development of pollen tubes after pollination and 4D-LFQ to detect the protein expression profiles of 'Wanfeng' and 'Nonpareil' pistils and in controls. The results showed that it took 24-36 h for the development of the pollen tube to 1/3 of the pistil, and a total of 7684 differentially accumulated proteins (DAPs) were identified in the pistil after pollinating for 36 h, of which 7022 were quantifiable. Bioinformatics analysis based on the function of DAPs, identified RNA polymerases (4 DAPs), autophagy (3 DAPs), oxidative phosphorylation (3 DAPs), and homologous recombination (2 DAPs) pathways associated with the self-incompatibility process. These results were confirmed by parallel reaction monitoring (PRM), protein interaction and bioinformatics analysis. Taken together, these results provide the involvement of serine/threonine kinase protein in the reaction of pollen tube recognition the nonself- and the self-S-RNase protein. SIGNIFICANCE: Gametophytic self-incompatibility (GSI) is controlled by the highly polymorphic S locus or S haplotype, with two linked self-incompatibility genes, one encoding the S-RNase protein of the pistil S-determinant and the other encoding the F-box/SLF/SFB (S haplotype-specific F-box protein) protein of the pollen S-determinant. The recognition mechanism between pollen tube SCF protein and pistil S-RNase protein is divided into nonself- and self-recognition hypothesis mechanisms. At present, two hypothetical mechanisms cannot explain the recognition between pollen and pistil well, so the mechanism of gametophytic self-incompatibility recognition is still not fully revealed. In this experiment, we investigated the molecular mechanism of pollen-pistil recognition in self-incompatibility using self- and nonself-pollinated pistils of almond cultivars 'Wanfeng' and 'Nonpareil'. Based on our results, we proposed a potential involvement of the MARK2 (serine/threonine kinase) protein in the reaction of pollen tube recognition of the nonself- and the self-S-RNase protein. It provides a new way to reveal how almond pollen tubes recognize the self and nonself S-RNase enzyme protein.

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

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

RNase-based self-incompatibility in cacti

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

A nonS-locus F-box gene breaks self-incompatibility in diploid potatoes

Potato is the third most important staple food crop. To address challenges associated with global food security, a hybrid potato breeding system, aimed at converting potato from a tuber-propagated tetraploid crop into a seed-propagated diploid crop through crossing inbred lines, is under development. However, given that most diploid potatoes are self-incompatible, this represents a major obstacle which needs to be addressed in order to develop inbred lines. Here, we report on a self-compatible diploid potato, RH89-039-16 (RH), which can efficiently induce a mating transition from self-incompatibility to self-compatibility, when crossed to self-incompatible lines. We identify the S-locusinhibitor (Sli) gene in RH, capable of interacting with multiple allelic variants of the pistil-specific S-ribonucleases (S-RNases). Further, Sli gene functions like a general S-RNase inhibitor, to impart SC to RH and other self-incompatible potatoes. Discovery of Sli now offers a path forward for the diploid hybrid breeding program.

Diploid potatoes are typically self-incompatible, complicating efforts to breed diploid cultivars. Here the authors report map-based cloning of the S-locus inhibitor (Sli) gene in potato which encodes a non S-locus F-box protein that is expressed in pollen and can functions like a general S-RNase inhibitor to overcome self-incompatibility.

Spread of self-compatibility constrained by an intrapopulation crossing barrier

Mating system transitions from self-incompatibility (SI) to self-compatibility (SC) are common in plants. In the absence of high levels of inbreeding depression, SC alleles are predicted to spread due to transmission advantage and reproductive assurance. We characterized mating system and pistil-expressed SI factors in 20 populations of the wild tomato species Solanum habrochaites from the southern half of the species range. We found that a single SI to SC transition is fixed in populations south of the Rio Chillon valley in central Peru. In these populations, SC correlated with the presence of the hab-6 S-haplotype that encodes a low activity S-RNase protein. We identified a single population segregating for SI/SC and hab-6. Intrapopulation crosses showed that hab-6 typically acts in the expected codominant fashion to confer SC. However, we found one specific S-haplotype (hab-10) that consistently rejects pollen of the hab-6 haplotype, and results in SI hab-6/hab-10 heterozygotes. We suggest that the hab-10 haplotype could act as a genetic mechanism to stabilize mixed mating in this population by presenting a disadvantage for the hab-6 haplotype. This barrier may represent a mechanism allowing for the persistence of SI when an SC haplotype appears in or invades a population.

Ornithine decarboxylase genes contribute to S-RNase-independent pollen rejection

Unilateral incompatibility (UI) manifests as pollen rejection in the pistil, typically when self-incompatible (SI) species are pollinated by self-compatible (SC) relatives. In the Solanaceae, UI occurs when pollen lack resistance to stylar S-RNases, but other, S-RNase-independent mechanisms exist. Pistils of the wild tomato Solanum pennellii LA0716 (SC) lack S-RNase yet reject cultivated tomato (Solanum lycopersicum, SC) pollen. In this cross, UI results from low pollen expression of a farnesyl pyrophosphate synthase gene (FPS2) in S. lycopersicum. Using pollen from fps2-/- loss-of-function mutants in S. pennellii, we identified a pistil factor locus, ui3.1, required for FPS2-based pollen rejection. We mapped ui3.1 to an interval containing 108 genes situated on the IL 3-3 introgression. This region includes a cluster of ornithine decarboxylase (ODC2) genes, with four copies in S. pennellii, versus one in S. lycopersicum. Expression of ODC2 transcript was 1,034-fold higher in S. pennellii than in S. lycopersicum styles. Pistils of odc2-/- knockout mutants in IL 3-3 or S. pennellii fail to reject fps2 pollen and abolish transmission ratio distortion (TRD) associated with FPS2. Pollen of S. lycopersicum express low levels of FPS2 and are compatible on IL 3-3 pistils, but incompatible on IL 12-3 × IL 3-3 hybrids, which express both ODC2 and ui12.1, a locus thought to encode the SI proteins HT-A and HT-B. TRD observed in F2 IL 12-3 × IL 3-3 points to additional ODC2-interacting pollen factors on both chromosomes. Thus, ODC2 genes contribute to S-RNase independent UI and interact genetically with ui12.1 to strengthen pollen rejection.

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

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

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

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

Ubiquitination of S4-RNase by S-LOCUS F-BOX LIKE2 Contributes to Self-Compatibility of Sweet Cherry 'Lapins'

Recent studies have shown that loss of pollen-S function in S₄' pollen from sweet cherry (Prunus avium) is associated with a mutation in an S haplotype-specific F-box4 (SFB4) gene. However, how this mutation leads to self-compatibility is unclear. Here, we examined this mechanism by analyzing several self-compatible sweet cherry varieties. We determined that mutated SFB4 (SFB4') in S4' pollen (pollen harboring the SFB4' gene) is approximately 6 kD shorter than wild-type SFB4 due to a premature termination caused by a four-nucleotide deletion. SFB4' did not interact with S-RNase. However, a protein in S4' pollen ubiquitinated S-RNase, resulting in its degradation via the 26S proteasome pathway, indicating that factors in S4' pollen other than SFB4 participate in S-RNase recognition and degradation. To identify these factors, we used S₄-RNase as a bait to screen S4' pollen proteins. Our screen identified the protein encoded by S ₄ -SLFL2, a low-polymorphic gene that is closely linked to the S-locus. Further investigations indicate that SLFL2 ubiquitinates S-RNase, leading to its degradation. Subcellular localization analysis showed that SFB4 is primarily localized to the pollen tube tip, whereas SLFL2 is not. When S ₄ -SLFL2 expression was suppressed by antisense oligonucleotide treatment in wild-type pollen tubes, pollen still had the capacity to ubiquitinate S-RNase; however, this ubiquitin-labeled S-RNase was not degraded via the 26S proteasome pathway, suggesting that SFB4 does not participate in the degradation of S-RNase. When SFB4 loses its function, S₄-SLFL2 might mediate the ubiquitination and degradation of S-RNase, which is consistent with the self-compatibility of S4' pollen.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1 Mb of the S₂ -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S_6a -haplotype of P. inflata and the S_N -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

NaTrxh is an essential protein for pollen rejection in Nicotiana by increasing S-RNase activity

In self-incompatible Solanaceae, the pistil protein S-RNase contributes to S-specific pollen rejection in conspecific crosses, as well as to rejecting pollen from foreign species or whole clades. However, S-RNase alone is not sufficient for either type of pollen rejection. We describe a thioredoxin (Trx) type h from Nicotiana alata, NaTrxh, which interacts with and reduces S-RNase in vitro. Here, we show that expressing a redox-inactive mutant, NaTrxhSS , suppresses both S-specific pollen rejection and rejection of pollen from Nicotiana plumbaginifolia. Biochemical experiments provide evidence that NaTrxh specifically reduces the Cys155 -Cys185 disulphide bond of SC10 -Rnase, resulting in a significant increase of its ribonuclease activity. This reduction and increase in S-RNase activity by NaTrxh helps to explain why S-RNase alone could be insufficient for pollen rejection.

Intraspecific Genetic Variation Underlying Postmating Reproductive Barriers between Species in the Wild Tomato Clade (Solanum sect. Lycopersicon)

A goal of speciation genetics is to understand how the genetic components underlying interspecific reproductive barriers originate within species. Unilateral incompatibility (UI) is a postmating prezygotic barrier in which pollen rejection in the female reproductive tract (style) occurs in only one direction of an interspecific cross. Natural variation in the strength of UI has been observed among populations within species in the wild tomato clade. In some cases, molecular loci underlying self-incompatibility (SI) are associated with this variation in UI, but the mechanistic connection between these intra- and inter-specific pollen rejection behaviors is poorly understood in most instances. We generated an F2 population between SI and SC genotypes of a single species, Solanum pennellii, to examine the genetic basis of intraspecific variation in UI against other species, and to determine whether loci underlying SI are genetically associated with this variation. We found that F2 individuals vary in the rate at which UI rejection occurs. One large effect QTL detected for this trait co-localized with the SI-determining S-locus. Moreover, individuals that expressed S-RNase-the S-locus protein involved in SI pollen rejection-in their styles had much more rapid UI responses compared with those without S-RNase protein. Our analysis shows that intraspecific variation at mate choice loci-in this case at loci that prevent self-fertilization-can contribute to variation in the expression of interspecific isolation, including postmating prezygotic barriers. Understanding the nature of such intraspecific variation can provide insight into the accumulation of these barriers between diverging lineages.

Intraspecific Genetic Variation Underlying Postmating Reproductive Barriers between Species in the Wild Tomato Clade (Solanum sect. Lycopersicon)

A goal of speciation genetics is to understand how the genetic components underlying interspecific reproductive barriers originate within species. Unilateral incompatibility (UI) is a postmating prezygotic barrier in which pollen rejection in the female reproductive tract (style) occurs in only one direction of an interspecific cross. Natural variation in the strength of UI has been observed among populations within species in the wild tomato clade. In some cases, molecular loci underlying self-incompatibility (SI) are associated with this variation in UI, but the mechanistic connection between these intra- and inter-specific pollen rejection behaviors is poorly understood in most instances. We generated an F2 population between SI and SC genotypes of a single species, Solanum pennellii, to examine the genetic basis of intraspecific variation in UI against other species, and to determine whether loci underlying SI are genetically associated with this variation. We found that F2 individuals vary in the rate at which UI rejection occurs. One large effect QTL detected for this trait co-localized with the SI-determining S-locus. Moreover, individuals that expressed S-RNase-the S-locus protein involved in SI pollen rejection-in their styles had much more rapid UI responses compared with those without S-RNase protein. Our analysis shows that intraspecific variation at mate choice loci-in this case at loci that prevent self-fertilization-can contribute to variation in the expression of interspecific isolation, including postmating prezygotic barriers. Understanding the nature of such intraspecific variation can provide insight into the accumulation of these barriers between diverging lineages.

Self-(In)compatibility Systems: Target Traits for Crop-Production, Plant Breeding, and Biotechnology

Self-incompatibility (SI) mechanisms prevent self-fertilization in flowering plants based on specific discrimination between self- and non-self pollen. Since this trait promotes outcrossing and avoids inbreeding it is a widespread mechanism of controlling sexual plant reproduction. Growers and breeders have effectively exploited SI as a tool for manipulating domesticated crops for thousands of years. However, only within the past thirty years have studies begun to elucidate the underlying molecular features of SI. The specific S-determinants and some modifier factors controlling SI have been identified in the sporophytic system exhibited by Brassica species and in the two very distinct gametophytic systems present in Papaveraceae on one side and in Solanaceae, Rosaceae, and Plantaginaceae on the other. Molecular level studies have enabled SI to SC transitions (and vice versa) to be intentionally manipulated using marker assisted breeding and targeted approaches based on transgene integration, silencing, and more recently CRISPR knock-out of SI-related factors. These scientific advances have, in turn, provided a solid basis to implement new crop production and plant breeding practices. Applications of self-(in)compatibility include widely differing objectives such as crop yield and quality improvement, marker-assisted breeding through SI genotyping, and development of hybrids for overcoming intra- and interspecific reproductive barriers. Here, we review scientific progress as well as patented applications of SI, and also highlight future prospects including further elucidation of SI systems, deepening our understanding of SI-environment relationships, and new perspectives on plant self/non-self recognition.

Genome-wide identification and characterization of the ALOG gene family in Petunia

BACKGROUND

The ALOG (Arabidopsis LSH1 and Oryza G1) family of proteins, namely DUF640 (domain of unknown function 640) domain proteins, were found in land plants. Functional characterization of a few ALOG members in model plants such as Arabidopsis and rice suggested they play important regulatory roles in plant development. The information about its evolution, however, is largely limited, and there was no any report on the ALOG genes in Petunia, an important ornamental species.

RESULTS

The ALOG genes were identified in four species of Petunia including P. axillaris, P. inflata, P. integrifolia, and P. exserta based on the genome and/or transcriptome databases, which were further confirmed by cloning from P. hybrida 'W115' (Mitchel diploid), a popular laboratorial petunia line susceptible to genetic transformation. Phylogenetic analysis indicated that Petunia ALOG genes (named as LSHs according to their closest Arabidopsis homologs) were grouped into four clades, which can be further divided into eight groups, and similar exon-intron structure and motifs are reflected in the same group. The PhLSH genes of hybrid petunia 'W115' were mainly derived from P. axillaris. The qPCR analysis revealed distinct spatial expression patterns among them suggesting potentially functional diversification. Moreover, over-expressing PhLSH7a and PhLSH7b in Arabidopsis uncovered their functions in the development of both vegetative and reproductive organs.

CONCLUSIONS

Petunia genome includes 11 ALOG genes that can be divided into eight distinct groups, and they also show different expression patterns. Among these genes, PhLSH7b and PhLSH7a play significant roles in plant growth and development, especially in fruit development. Our results provide new insight into the evolution of ALOG gene family and have laid a good foundation for the study of petunia LSH gene in the future.

Identification and Characterization of the FLOWERING LOCUS T/TERMINAL FLOWER 1 Gene Family in Petunia

The phosphatidylethanolamine-binding protein (PEBP) gene family exists in all eukaryote kingdoms, with three subfamilies: FT (FLOWERING LOCUS T)-like, TFL1 (TERMINAL FLOWER 1)-like, and MFT (MOTHER OF FT AND TFL1)-like. FT genes promote flowering, TFL1 genes act as a repressor of the floral transition, and MFT genes have functions in flowering promotion and regulating seed germination. We identified and characterized orthologs of the Arabidopsis FT/TFL1 gene family in petunia to elucidate their expression patterns and evolution. Thirteen FT/TFL1-like genes were isolated from petunia, with the five FT-like genes mainly expressed in leaves. The circadian rhythms of five FT-like genes and PhCO (petunia CONSTANS ortholog) were figured out. The expression of PhFT1 was contrary to that of PhFT2, PhFT3, PhFT4, and PhFT5. PhCO had a circadian clock different from Arabidopsis CO, but coincided with PhFT1; it decreased in daytime and accumulated at night. Two of the FT-like genes with differential circadian rhythm and higher expression levels, PhFT1 and PhFT4, were used to transform Arabidopsis. Eventually, overexpressing PhFT1 strongly delayed flowering, whereas overexpression of PhFT4 produced extremely early-flowering phenotype. Different from previous reports, PhTFL1a, PhTFL1b, and PhTFL1c were relatively highly expressed in roots. Taken together, this study demonstrates that petunia FT-like genes, like FT, are able to respond to photoperiod. The expression pattern of FT/TFL1 gene family in petunia contributes to a new insight into the functional evolution of this gene family.

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.

Mate Selection in Self-Compatible Wild Tobacco Results from Coordinated Variation in Homologous Self-Incompatibility Genes

In flowering plants, intraspecific mate preference is frequently related to mating systems: the rejection of self pollen in self-incompatible (SI) plants that prevents inbreeding is one of the best described examples. However, in other mating systems, more nuanced patterns of pollen rejection occur. In the self-compatible (SC) Nicotiana attenuata, in which SI is not found and all crosses are compatible, certain pollen genotypes are consistently selected in mixed pollinations. However, the molecular mechanisms of this polyandrous mate selection remain unknown. Style-expressed NaS-like-RNases and pollen-expressed NaSLF-like genes, homologous to SI factors in Solanaceae, were identified and examined for a role in N. attenuata's mate selection. A comparison of two NaS-like-RNases and six NaSLF-like genes among 26 natural accessions revealed specific combinations of co-expression and direct protein-protein interactions. To evaluate their role in mate selection, we silenced the expression of specific NaS-like-RNases and NaSLF-like proteins and conducted diagnostic binary mixed pollinations and mixed pollinations with 14 different non-self pollen donors. Styles expressing particular combinations of NaS-like-RNases selected mates from plants with corresponding NaS-like-RNase expression patterns, while styles lacking NaS-like-RNase expression were non-selective in their fertilizations, which reflected the genotype ratios of pollen mixtures deposited on the stigmas. DNA methylation could account for some of the observed variation in stylar NaS-like-RNase patterns. We conclude that the S-RNase-SLF recognition mechanism plays a central role in polyandrous mate selection in this self-compatible species. These results suggest that after the SI-SC transition, natural variation of SI homologous genes was repurposed to mediate intraspecific mate selection.

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.

S-Locus F-Box Proteins Are Solely Responsible for S-RNase-Based Self-Incompatibility of Petunia Pollen

Self-incompatibility (SI) in Petunia is regulated by a polymorphic S-locus. For each S-haplotype, the S-locus contains a pistil-specific S-RNase gene and multiple pollen-specific S-locus F-box (SLF) genes. Both gain-of-function and loss-of-function experiments have shown that S-RNase alone regulates pistil specificity in SI. Gain-of-function experiments on SLF genes suggest that the entire suite of encoded proteins constitute the pollen specificity determinant. However, clear-cut loss-of-function experiments must be performed to determine if SLF proteins are essential for SI of pollen. Here, we used CRISPR/Cas9 to generate two frame-shift indel alleles of S₂ -SLF1 (SLF1 of S₂ -haplotype) in S₂S₃ plants of P. inflata and examined the effect on the SI behavior of S₂ pollen. In the absence of a functional S₂-SLF1, S₂ pollen was either rejected by or remained compatible with pistils carrying one of eight normally compatible S-haplotypes. All results are consistent with interaction relationships between the 17 SLF proteins of S₂ -haplotype and these eight S-RNases that had been determined by gain-of-function experiments performed previously or in this work. Our loss-of-function results provide definitive evidence that SLF proteins are solely responsible for SI of pollen, and they reveal their diverse and complex interaction relationships with S-RNases to maintain SI while ensuring cross-compatibility.

Multiple model species selection for transcriptomics analysis of non-model organisms

Background

Transcriptomic sequencing (RNA-seq) related applications allow for rapid explorations due to their high-throughput and relatively fast experimental capabilities, providing unprecedented progress in gene functional annotation, gene regulation analysis, and environmental factor verification. However, with increasing amounts of sequenced reads and reference model species, the selection of appropriate reference species for gene annotation has become a new challenge.

Methods

We proposed a novel approach for finding the most effective reference model species through taxonomic associations and ultra-conserved orthologous (UCO) gene comparisons among species. An online system for multiple species selection (MSS) for RNA-seq differential expression analysis was developed, and comprehensive genomic annotations from 291 reference model eukaryotic species were retrieved from the RefSeq, KEGG, and UniProt databases.

Results

Using the proposed MSS pipeline, gene ontology and biological pathway enrichment analysis can be efficiently achieved, especially in the case of transcriptomic analysis of non-model organisms. The results showed that the proposed method solved problems related to limitations in annotation information and provided a roughly twenty-fold reduction in computational time, resulting in more accurate results than those of traditional approaches of using a single model reference species or the large non-redundant reference database.

Conclusions

Selection of appropriate reference model species helps to reduce missing annotation information, allowing for more comprehensive results than those obtained with a single model reference species. In addition, adequate model species selection reduces the computational time significantly while retaining the same order of accuracy. The proposed system indeed provides superior performance by selecting appropriate multiple species for transcriptomic analysis compared to traditional approaches.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2278-z) contains supplementary material, which is available to authorized users.

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.

Genome-wide identification and characterization of the SBP-box gene family in Petunia

Background

SQUAMOSA PROMOTER BINDING PROTEIN (SBP)-box genes encode a family of plant-specific transcription factors (TFs) that play important roles in many growth and development processes including phase transition, leaf initiation, shoot and inflorescence branching, fruit development and ripening etc. The SBP-box gene family has been identified and characterized in many species, but has not been well studied in Petunia, an important ornamental genus.

Results

We identified 21 putative SPL genes of Petunia axillaris and P. inflata from the reference genome of P. axillaris N and P. inflata S6, respectively, which were supported by the transcriptome data. For further confirmation, all the 21 genes were also cloned from P. hybrida line W115 (Mitchel diploid). Phylogenetic analysis based on the highly conserved SBP domains arranged PhSPLs in eight groups, analogous to those from Arabidopsis and tomato. Furthermore, the Petunia SPL genes had similar exon-intron structure and the deduced proteins contained very similar conserved motifs within the same subgroup. Out of 21 PhSPL genes, fourteen were predicted to be potential targets of PhmiR156/157, and the putative miR156/157 response elements (MREs) were located in the coding region of group IV, V, VII and VIII genes, but in the 3’-UTR regions of group VI genes. SPL genes were also identified from another two wild Petunia species, P. integrifolia and P. exserta, based on their transcriptome databases to investigate the origin of PhSPLs. Phylogenetic analysis and multiple alignments of the coding sequences of PhSPLs and their orthologs from wild species indicated that PhSPLs were originated mainly from P. axillaris. qRT-PCR analysis demonstrated differential spatiotemperal expression patterns of PhSPL genes in petunia and many were expressed predominantly in the axillary buds and/or inflorescences. In addition, overexpression of PhSPL9a and PhSPL9b in Arabidopsis suggested that these genes play a conserved role in promoting the vegetative-to-reproductive phase transition.

Conclusion

Petunia genome contains at least 21 SPL genes, and most of the genes are expressed in different tissues. The PhSPL genes may play conserved and diverse roles in plant growth and development, including flowering regulation, leaf initiation, axillary bud and inflorescence development. This work provides a comprehensive understanding of the SBP-box gene family in Petunia and lays a significant foundation for future studies on the function and evolution of SPL genes in petunia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4537-9) contains supplementary material, which is available to authorized users.

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.

A farnesyl pyrophosphate synthase gene expressed in pollen functions in S-RNase-independent unilateral incompatibility

Multiple independent and overlapping pollen rejection pathways contribute to unilateral interspecific incompatibility (UI). In crosses between tomato species, pollen rejection usually occurs when the female parent is self-incompatible (SI) and the male parent self-compatible (SC) (the 'SI × SC rule'). Additional, as yet unknown, UI mechanisms are independent of self-incompatibility and contribute to UI between SC species or populations. We identified a major quantitative trait locus on chromosome 10 (ui10.1) which affects pollen-side UI responses in crosses between cultivated tomato, Solanum lycopersicum, and Solanum pennelliiLA0716, both of which are SC and lack S-RNase, the pistil determinant of S-specificity in Solanaceae. Here we show that ui10.1 is a farnesyl pyrophosphate synthase gene (FPS2) expressed in pollen. Expression is about 18-fold higher in pollen of S. pennellii than in S. lycopersicum. Pollen with the hypomorphic S. lycopersicum allele is selectively eliminated on pistils of the F₁ hybrid, leading to transmission ratio distortion in the F₂ progeny. CRISPR/Cas9-generated knockout mutants (fps2) in S. pennelliiLA0716 are self-sterile due to pollen rejection, but mutant pollen is fully functional on pistils of S. lycopersicum. F₂ progeny of S. lycopersicum × S. pennellii (fps2) show reversed transmission ratio distortion due to selective elimination of pollen bearing the knockout allele. Overexpression of FPS2 in S. lycopersicum pollen rescues the pollen elimination phenotype. FPS2-based pollen selectivity does not involve S-RNase and has not been previously linked to UI. Our results point to an entirely new mechanism of interspecific pollen rejection in plants.

Use of Domain-Swapping to Identify Candidate Amino Acids Involved in Differential Interactions between Two Allelic Variants of Type-1 S-Locus F-Box Protein and S3-RNase in Petunia inflata

Petunia inflata possesses a self-incompatibility (SI) mechanism, which involves S-RNase and multiple S-locus F-box (SLF) genes at the polymorphic S-locus. For a given S-haplotype, each SLF is thought to interact with some of its non-self S-RNases, but not with its self S-RNase. In this work, we studied an allelic pair of SLF1, S2-SLF1 and S3-SLF1, which differ in 44 amino acids and show differential interactions with S3-RNase. We first used an in vivo transgenic assay to determine whether four chimeric proteins of S2-SLF1 and S3-SLF1, each with one of the three functional domains swapped, interact with S3-RNase. The results narrowed the candidate amino acids for specific interaction of S2-SLF1 with S3-RNase to the 16 in domain FD3. We then examined seven additional chimeric proteins by dividing FD3 into two subdomains and four mini-domains (A, B, C and D). The results further narrowed the candidate amino acids to four in mini-domain A and four in mini-domain D. Molecular modeling of interactions between S3-RNase and S2-SLF1 revealed that three of these eight are at the interaction surface, and all three are conserved in S1-SLF1 and S6a-SLF1, both of which interact with S3-RNase based on the in vivo transgenic assay. Three of the chimeric proteins were used for the in vivo transgenic assay to determine whether FD3 alone contains the amino acids required for S2-SLF1 to interact with S7-RNase and S13-RNase. The results revealed the diversity and complexity of interactions between SLF proteins and S-RNases.

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

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

The evolutionary history of plant T2/S-type ribonucleases

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

Duplicated pollen-pistil recognition loci control intraspecific unilateral incompatibility in Brassica rapa

In plants, cell-cell recognition is a crucial step in the selection of optimal pairs of gametes to achieve successful propagation of progeny. Flowering plants have evolved various genetic mechanisms, mediated by cell-cell recognition, to enable their pistils to reject self-pollen, thus preventing inbreeding and the consequent reduced fitness of progeny (self-incompatibility, SI), and to reject foreign pollen from other species, thus maintaining species identity (interspecific incompatibility)¹. In the genus Brassica, the SI system is regulated by an S-haplotype-specific interaction between a stigma-expressed female receptor (S receptor kinase, SRK) and a tapetum cell-expressed male ligand (S locus protein 11, SP11), encoded by their respective polymorphic genes at the S locus^2-6. However, the molecular mechanism for recognition of foreign pollen, leading to reproductive incompatibility, has not yet been identified. Here, we show that recognition between a novel pair of proteins, a pistil receptor SUI1 (STIGMATIC UNILATERAL INCOMPATIBILITY 1) and a pollen ligand PUI1 (POLLEN UNILATERAL INCOMPATIBILITY 1), triggers unilateral reproductive incompatibility between plants of two geographically distant self-incompatible Brassica rapa lines, even though crosses would be predicted to be compatible based on the S haplotypes of pollen and stigma. Interestingly, SUI1 and PUI1 are similar to the SI genes, SRK and SP11, respectively, and are maintained as cryptic incompatibility genes in these two populations. The duplication of the SRK and SP11 followed by reciprocal loss in different populations would provide a molecular mechanism of the emergence of a reproductive barrier in allopatry.

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

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

The genetic architecture of shoot-root covariation during seedling emergence of a desert tree, Populus euphratica

The coordination of shoots and roots is critical for plants to adapt to changing environments by fine-tuning energy production in leaves and the availability of water and nutrients from roots. To understand the genetic architecture of how these two organs covary during developmental ontogeny, we conducted a mapping experiment using Euphrates poplar (Populus euphratica), a so-called hero tree able to grow in the desert. We geminated intraspecific F₁ seeds of Euphrates Poplar individually in a tube to obtain a total of 370 seedlings, whose shoot and taproot lengths were measured repeatedly during the early stage of growth. By fitting a growth equation, we estimated asymptotic growth, relative growth rate, the timing of inflection point and duration of linear growth for both shoot and taproot growth. Treating these heterochronic parameters as phenotypes, a univariate mapping model detected 19 heterochronic quantitative trait loci (hQTLs), of which 15 mediate the forms of shoot growth and four mediate taproot growth. A bivariate mapping model identified 11 pleiotropic hQTLs that determine the covariation of shoot and taproot growth. Most QTLs detected reside within the region of candidate genes with various functions, thus confirming their roles in the biochemical processes underlying plant growth.

Mating system transitions in Solanum habrochaites impact interactions between populations and species

In plants, transitions in mating system from outcrossing to self-fertilization are common; however, the impact of these transitions on interspecific and interpopulation reproductive barriers is not fully understood. We examined the consequences of mating system transition for reproductive barriers in 19 populations of the wild tomato species Solanum habrochaites. We identified S. habrochaites populations with self-incompatible (SI), self-compatible (SC) and mixed population (MP) mating systems, and characterized pollen-pistil interactions among S. habrochaites populations and between S. habrochaites and other tomato species. We examined the relationship between mating system, floral morphology, interspecific and interpopulation compatibility and pistil SI factors. We documented five distinct phenotypic groups by combining reproductive behavior with molecular data. Transitions from SI to MP were not associated with weakened interspecific reproductive barriers or loss of known pistil SI factors. However, transitions to SC at the northern range margin were accompanied by loss of S-RNase, smaller flowers, and weakened (or absent) interspecific pollen-pistil barriers. Finally, we identified a subset of SC populations that exhibited a partial interpopulation reproductive barrier with central SI populations. Our results support the hypothesis that shifts in mating system, followed by additional loss-of-function mutations, impact reproductive barriers within and between species.

Pollen-Pistil Interactions and Their Role in Mate Selection

Pollen-pistil interactions contribute to mate selection at the postmating, prezygotic level.

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

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

All 17 S-locus F-box proteins of the S2 - and S3 -haplotypes of Petunia inflata are assembled into similar SCF complexes with a specific function in self-incompatibility

The collaborative non-self-recognition model for S-RNase-based self-incompatibility predicts that multiple S-locus F-box proteins (SLFs) produced by pollen of a given S-haplotype collectively mediate ubiquitination and degradation of all non-self S-RNases, but not self S-RNases, in the pollen tube, thereby resulting in cross-compatible pollination but self-incompatible pollination. We had previously used pollen extracts containing GFP-fused S2 -SLF1 (SLF1 with an S2 -haplotype) of Petunia inflata for co-immunoprecipitation (Co-IP) and mass spectrometry (MS), and identified PiCUL1-P (a pollen-specific Cullin1), PiSSK1 (a pollen-specific Skp1-like protein) and PiRBX1 (a conventional Rbx1) as components of the SCF(S) (2-) (SLF) (1) complex. Using pollen extracts containing PiSSK1:FLAG:GFP for Co-IP/MS, we identified two additional SLFs (SLF4 and SLF13) that were assembled into SCF(SLF) complexes. As 17 SLF genes (SLF1 to SLF17) have been identified in S2 and S3 pollen, here we examined whether all 17 SLFs are assembled into similar complexes and, if so, whether these complexes are unique to SLFs. We modified the previous Co-IP/MS procedure, including the addition of style extracts from four different S-genotypes to pollen extracts containing PiSSK1:FLAG:GFP, to perform four separate experiments. The results taken together show that all 17 SLFs and an SLF-like protein, SLFLike1 (encoded by an S-locus-linked gene), co-immunoprecipitated with PiSSK1:FLAG:GFP. Moreover, of the 179 other F-box proteins predicted by S2 and S3 pollen transcriptomes, only a pair with 94.9% identity and another pair with 99.7% identity co-immunoprecipitated with PiSSK1:FLAG:GFP. These results suggest that SCF(SLF) complexes have evolved specifically to function in self-incompatibility.