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

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.

A pollen selection system links self and interspecific incompatibility in the Brassicaceae

Self-incompatibility and recurrent transitions to self-compatibility have shaped the extant mating systems underlying the nonrandom mating critical for speciation in angiosperms. Linkage between self-incompatibility and speciation is illustrated by the shared pollen rejection pathway between self-incompatibility and interspecific unilateral incompatibility (UI) in the Brassicaceae. However, the pollen discrimination system that activates this shared pathway for heterospecific pollen rejection remains unknown. Here we show that Stigma UI3.1, the genetically identified stigma determinant of UI in Arabidopsis lyrata × Arabidopsis arenosa crosses, encodes the S-locus-related glycoprotein 1 (SLR1). Heterologous expression of A. lyrata or Capsella grandiflora SLR1 confers on some Arabidopsis thaliana accessions the ability to discriminate against heterospecific pollen. Acquisition of this ability also requires a functional S-locus receptor kinase (SRK), whose ligand-induced dimerization activates the self-pollen rejection pathway in the stigma. SLR1 interacts with SRK and interferes with SRK homomer formation. We propose a pollen discrimination system based on competition between basal or ligand-induced SLR1-SRK and SRK-SRK complex formation. The resulting SRK homomer levels would be sensed by the common pollen rejection pathway, allowing discrimination among conspecific self- and cross-pollen as well as heterospecific pollen. Our results establish a mechanistic link at the pollen recognition phase between self-incompatibility and interspecific incompatibility.

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.

Phase separation of S-RNase promotes self-incompatibility in Petunia hybrida

Self-incompatibility (SI) is an intraspecific reproductive barrier widely present in angiosperms. The SI system with the broadest occurrence in angiosperms is based on an S-RNase linked to a cluster of multiple S-locus F-box (SLF) genes found in the Solanaceae, Plantaginaceae, Rosaceae, and Rutaceae. Recent studies reveal that non-self S-RNase is degraded by the Skip Cullin F-box (SCF)SLF-mediated ubiquitin-proteasome system in a collaborative manner in Petunia, but how self-RNase functions largely remains mysterious. Here, we show that S-RNases form S-RNase condensates (SRCs) in the self-pollen tube cytoplasm through phase separation and the disruption of SRC formation breaks SI in self-incompatible Petunia hybrida. We further find that the pistil SI factors of a small asparagine-rich protein HT-B and thioredoxin h together with a reduced state of the pollen tube all promote the expansion of SRCs, which then sequester several actin-binding proteins, including the actin polymerization factor PhABRACL, the actin polymerization activity of which is reduced by S-RNase in vitro. Meanwhile, we find that S-RNase variants lacking condensation ability fail to recruit PhABRACL and are unable to induce actin foci formation required for pollen tube growth inhibition. Taken together, our results demonstrate that phase separation of S-RNase promotes SI response in P. hybrida, revealing a new mode of S-RNase action.

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 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.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Overexpression of the Ginkgo biloba dihydroflavonol 4-reductase gene GbDFR6 results in the self-incompatibility-like phenotypes in transgenic tobacco

Although flavonoids play multiple roles in plant growth and development, the involvement in plant self-incompatibility (SI) have not been reported. In this research, the fertility of transgenic tobacco plants overexpressing the Ginkgo biloba dihydroflavonol 4-reductase gene, GbDFR6, were investigated. To explore the possible physiological defects leading to the failure of embryo development in transgenic tobacco plants, functions of pistils and pollen grains were examined. Transgenic pistils pollinated with pollen grains from another tobacco plants (either transgenic or wild-type), developed full of well-developed seeds. In contrast, in self-pollinated transgenic tobacco plants, pollen-tube growth was arrested in the upper part of the style, and small abnormal seeds developed without fertilization. Although the mechanism remains unclear, our research may provide a valuable method to create SI tobacco plants for breeding.

Intra- and inter-specific reproductive barriers in the tomato clade

Tomato (Solanum lycopersicum L.) domestication and later introduction into Europe resulted in a genetic bottleneck that reduced genetic variation. Crosses with other wild tomato species from the Lycopersicon clade can be used to increase genetic diversity and improve important agronomic traits such as stress tolerance. However, many species in the Lycopersicon clade have intraspecific and interspecific incompatibility, such as gametophytic self-incompatibility and unilateral incompatibility. In this review, we provide an overview of the known incompatibility barriers in Lycopersicon. We begin by addressing the general mechanisms self-incompatibility, as well as more specific mechanisms in the Rosaceae, Papaveraceae, and Solanaceae. Incompatibility in the Lycopersicon clade is discussed, including loss of self-incompatibility, species exhibiting only self-incompatibility and species presenting both self-compatibility and self-incompatibility. We summarize unilateral incompatibility in general and specifically in Lycopersicon, with details on the 'self-compatible x self-incompatible' rule, implications of self-incompatibility in unilateral incompatibility and self-incompatibility-independent pathways of unilateral incompatibility. Finally, we discuss advances in the understanding of compatibility barriers and their implications for tomato breeding.

Acetylation of inorganic pyrophosphatase by S-RNase signaling induces pollen tube tip swelling by repressing pectin methylesterase

Self-incompatibility (SI) is a widespread genetically determined system in flowering plants that prevents self-fertilization to promote gene flow and limit inbreeding. S-RNase-based SI is characterized by the arrest of pollen tube growth through the pistil. Arrested pollen tubes show disrupted polarized growth and swollen tips, but the underlying molecular mechanism is largely unknown. Here, we demonstrate that the swelling at the tips of incompatible pollen tubes in pear (Pyrus bretschneideri [Pbr]) is mediated by the SI-induced acetylation of the soluble inorganic pyrophosphatase (PPA) PbrPPA5. Acetylation at Lys-42 of PbrPPA5 by the acetyltransferase GCN5-related N-acetyltransferase 1 (GNAT1) drives accumulation of PbrPPA5 in the nucleus, where it binds to the transcription factor PbrbZIP77, forming a transcriptional repression complex that inhibits the expression of the pectin methylesterase (PME) gene PbrPME44. The function of PbrPPA5 as a transcriptional repressor does not require its PPA activity. Downregulating PbrPME44 resulted in increased levels of methyl-esterified pectins in growing pollen tubes, leading to swelling at their tips. These observations suggest a mechanism for PbrPPA5-driven swelling at the tips of pollen tubes during the SI response. The targets of PbrPPA5 include genes encoding cell wall-modifying enzymes, which are essential for building a continuous sustainable mechanical structure for pollen tube growth.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

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.

Haplotype-resolved genomes of two buckwheat crops provide insights into their contrasted rutin concentrations and reproductive systems

BACKGROUND

Two widely cultivated annual buckwheat crops, Fagopyrum esculentum and F. tataricum, differ from each other in both rutin concentration and reproductive system. However, the underlying genetic mechanisms remain poorly elucidated.

RESULTS

Here, we report the first haplotype-resolved chromosome-level genome assemblies of the two species. Two haplotype genomes of F. esculentum were assembled as 1.23 and 1.19 Gb with N50 = 9.8 and 12.4 Mb, respectively; the two haplotype genomes of F. tataricum were 453.7 and 446.2 Mb with N50 = 50 and 30 Mb, respectively. We further annotated protein-coding genes of each haplotype genome based on available gene sets and 48 newly sequenced transcriptomes. We found that more repetitive sequences, especially expansion of long terminal repeat retrotransposons (LTR-RTs), contributed to the large genome size of F. esculentum. Based on the well-annotated sequences, gene expressions, and luciferase experiments, we identified the sequence mutations of the promoter regions of two key genes that are likely to have greatly contributed to the high rutin concentration and selfing reproduction in F. tartaricum.

CONCLUSIONS

Our results highlight the importance of high-quality genomes to identify genetic mutations underlying phenotypic differences between closely related species. F. tataricum may have been experienced stronger selection than F. esculentum through choosing these two non-coding alleles for the desired cultivation traits. These findings further suggest that genetic manipulation of the non-coding promoter regions could be widely employed for breeding buckwheat and other crops.

Genome-wide screening of the RNase T2 gene family and functional analyses in jujube (Ziziphus jujuba Mill.)

BACKGROUND

Ribonuclease (RNase T2) plays crucial roles in plant evolution and breeding. However, there have been few studies on the RNase T2 gene family in Ziziphus jujuba Mill., one of important dried fruit tree species. Recently, the released sequences of the reference genome of jujube provide a good chance to perform genome-wide identification and characterization of ZjRNase gene family in the jujube.

RESULTS

In this study, we identified four members of RNase T2 in jujube distributed on three chromosomes and unassembled chromosomes. They all contained two conserved sites (CASI and CASII). Analysis of the phylogenetic relationships revealed that the RNase T2 genes in jujube could be divided into two groups: ZjRNase1 and ZjRNase2 belonged to class I, while ZjRNase3 and ZjRNase4 belonged to class II. Only ZjRNase1 and ZjRNase2 expression were shown by the jujube fruit transcriptome analysis. So ZjRNase1 and ZjRNase2 were selected functional verification by overexpression transformation of Arabidopsis. The overexpression of these two genes led to an approximately 50% reduction in seed number, which deserve further attention. Moreover, the leaves of the ZjRNase1 overexpression transgenic lines were curled and twisted. Overexpression of ZjRNase2 resulted in shortened and crisp siliques and the production of trichomes, and no seeds were produced.

CONCLUSION

In summary, these findings will provide new insights into the molecular mechanisms of low number of hybrid seeds in jujube and a reference for the future molecular breeding of jujube.

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.

Gene coexpression analysis reveals key pathways and hub genes related to late-acting self-incompatibility in Camellia oleifera

INTRODUCTION

Self-incompatibility (SI) is an important strategy for plants to maintain abundant variation to enhance their adaptability to the environment. Camellia oleifera is one of the most important woody oil plants and is widely cultivated in China. Late acting self-incompatibility (LSI) in C. oleifera results in a relatively poor fruit yield in the natural state, and understanding of the LSI mechanism remains limited.

METHODS

To better understand the molecular expression and gene coexpression network in the LSI reaction in C. oleifera, we conducted self- and cross-pollination experiments at two different flower bud developmental stages (3-4 d before flowering and 1 d before flowering), and cytological observation, fruit setting rate (FSR) investigation and RNA-Seq analysis were performed to investigate the mechanism of the male -female interaction and identify hub genes responsible for the LSI in C. oleifera.

RESULTS

Based on the 21 ovary transcriptomes, a total of 7669 DEGs were identified after filtering out low-expression genes. Weighted gene coexpression network analysis (WGCNA) divided the DEGs into 15 modules. Genes in the blue module (1163 genes) were positively correlated with FSR, and genes in the pink module (339 genes) were negatively correlated with FSR. KEGG analysis indicated that flavonoid biosynthesis, plant MAPK signaling pathways, ubiquitin-mediated proteolysis, and plant-pathogen interaction were the crucial pathways for the LSI reaction. Fifty four transcription factors (TFs) were obtained in the two key modules, and WRKY and MYB were potentially involved in the LSI reaction in C. oleifera. Network establishment indicated that genes encoding G-type lectin S-receptor-like serine (lecRLK), isoflavone 3'-hydroxylase-like (CYP81Q32), cytochrome P450 87A3-like (CYP87A3), and probable calcium-binding protein (CML41) were the hub genes that positively responded to the LSI reaction. The other DEGs inside the two modules, including protein RALF-like 10 (RALF), F-box and pectin acetylesterase (MTERF5), might also play vital roles in the LSI reaction in C. oleifera.

DISCUSSION

Overall, our study provides a meaningful resource for gene network studies of the LSI reaction process and subsequent analyses of pollen-pistil interactions and TF roles in the LSI reaction, and it also provides new insights for exploring the mechanisms of the LSI response.

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.

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.