Mechanism of self/nonself-discrimination in Brassica self-incompatibility

Fine mapping of a major co-localized QTL associated with self-incompatibility identified in two F2 populations (broccoli × cauliflower and cauliflower × Chinese kale)

KEY MESSAGE

A major QTL responsible for self-incompatibility was stably identified in two F2 populations. Through fine mapping and qRT-PCR analysis, ARK3 emerged as the most promising candidate gene, playing a pivotal role in regulating self-incompatibility in Brassica oleracea. Self-incompatibility (SI) is a common phenomenon in Brassica oleracea species, which can maintain genetic diversity but will also limit seed production. Although the S locus has been extensively studied in Arabidopsis and some Brassicaceae crops, map-based cloning of self-incompatibility genes has not been conducted in Brassica oleracea, such as cauliflower and broccoli. In the present study, we identified a major co-localized QTL on chromosome C6 that control SI in two F2 populations derived from intervarietal crosses: broccoli × cauliflower (CL_F2) and cauliflower × Chinese kale (CJ_F2). Subsequently, this QTL was narrowed down to 168.5 Kb through fine mapping using 3,429 F2:3 progenies and 12 available KASP markers. Within this 168.5 Kb region, BolC6t39084H, a homologue of Arabidopsis ARK3, could be a candidate gene that plays a key role in regulating SI in B. oleracea species. This finding can pave the way for an in-depth understanding of the molecular mechanisms underlying SI, and will contribute to the seed production of B. oleracea vegetables.

Analysis of randomly mutated AlSRKb genes reveals that most loss-of-function mutations cause defects in plasma membrane localization

Only very limited information is available on why some nonsynonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on SRK which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 300 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed using error-prone PCR and the genotype and self-incompatibility phenotype of each transformant were determined. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. Our findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of AlSRKb function. In addition, we examined whether the RandomForest and Extreme Gradient Boosting methods could predict the self-incompatibility phenotypes of SRK mutants.

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.

LORE receptor homomerization is required for 3-hydroxydecanoic acid-induced immune signaling and determines the natural variation of immunosensitivity within the Arabidopsis genus

The S-domain-type receptor-like kinase (SD-RLK) LIPOOLIGOSACCHARIDE-SPECIFIC REDUCED ELICITATION (LORE) from Arabidopsis thaliana is a pattern recognition receptor that senses medium-chain 3-hydroxy fatty acids, such as 3-hydroxydecanoic acid (3-OH-C10:0), to activate pattern-triggered immunity. Here, we show that LORE homomerization is required to activate 3-OH-C10:0-induced immune signaling. Fluorescence lifetime imaging in Nicotiana benthamiana demonstrates that AtLORE homomerizes via the extracellular and transmembrane domains. Co-expression of AtLORE truncations lacking the intracellular domain exerts a dominant negative effect on AtLORE signaling in both N. benthamiana and A. thaliana, highlighting that homomerization is essential for signaling. Screening for 3-OH-C10:0-induced reactive oxygen species production revealed natural variation within the Arabidopsis genus. Arabidopsis lyrata and Arabidopsis halleri do not respond to 3-OH-C10:0, although both possess a putative LORE ortholog. Both LORE orthologs have defective extracellular domains that bind 3-OH-C10:0 to a similar level as AtLORE, but lack the ability to homomerize. Thus, ligand binding is independent of LORE homomerization. Analysis of AtLORE and AlyrLORE chimera suggests that the loss of AlyrLORE homomerization is caused by several amino acid polymorphisms across the extracellular domain. Our findings shed light on the activation mechanism of LORE and the loss of 3-OH-C10:0 perception within the Arabidopsis genus.

Self-incompatibility phenotypes of SRK mutants can be predicted with high accuracy

Only very limited information is available on why some non-synonymous variants severely alter gene function while others have no effect. To identify the characteristic features of mutations that strongly influence gene function, this study focused on S-locus receptor kinase, SRK, which encodes a highly polymorphic receptor kinase expressed in stigma papillary cells that underlies a female determinant of self-incompatibility in Brassicaceae. A set of 299 Arabidopsis thaliana transformants expressing mutated SRKb from A. lyrata was constructed and analyzed to determine the genotype and self-incompatibility phenotype of each transformant. Almost all the transformants showing the self-incompatibility defect contained mutations in AlSRKb that altered localization to the plasma membrane. The observed mutations occurred in amino acid residues that were highly conserved across S haplotypes and whose predicted locations were in the interior of the protein. These mutations were likely to underlie the self-incompatibility defect as they caused significant changes to amino acid properties. Such findings suggested that mutations causing the self-incompatibility defect were more likely to result from changes to AlSRKb biosynthesis than from loss of function. In addition, this study showed the RandomForest and Extreme Gradient Boosting methods could predict self-incompatibility phenotypes of SRK mutants with high accuracy.

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.

The self-compatibility is acquired after polyploidization: a case study of Brassica napus self-incompatible trilinear hybrid breeding system

Self-incompatibility plays a vital role in angiosperms, by preventing inbreeding depression and maintaining genetic diversity within populations. Following polyploidization, many angiosperm species transition from self-incompatibility to self-compatibility. Here, we investigated the S-locus in Brassicaceae and identified distinct origins for the sRNA loci, SMI and SMI2 (SCR Methylation Inducer 1 and 2), within the S-locus. The SMI loci were found to be widespread in Cruciferae, whereas the SMI2 loci were exclusive to Brassica species. Additionally, we discovered four major S-haplotypes (BnS-1, BnS-6, BnS-7, and BnS-1300) in rapeseed. Overexpression of BnSMI-1 in self-incompatible Brassica napus ('S-70S1300S6 ') resulted in a significant increase in DNA methylation in the promoter regions of BnSCR-6 and BnSCR-1300, leading to self-compatibility. Conversely, by overexpressing a point mutation of BnSmi-1 in the 'S-70S1300S6 ' line, we observed lower levels of DNA methylation in BnSCR-6 and BnSCR-1300 promoters. Furthermore, the overexpression of BnSMI2-1300 in the 'SI-326S7S6 ' line inhibited the expression of BnSCR-7 through transcriptional repression of the Smi2 sRNA from the BnS-1300 haplotype. Our study demonstrates that the self-compatibility of rapeseed is determined by S-locus sRNA-mediated silencing of SCR after polyploidization, which helps to further breed self-incompatible or self-compatible rapeseed lines, thereby facilitating the utilization of heterosis.

Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. SI is controlled by a single S-locus with multiple haplotypes (S-haplotypes). When the pistil and pollen share the same S-haplotype, the pollen is recognized as self and rejected by the pistil. This review introduces our research on Brassicaceae and Solanaceae SI systems to identify the S-determinants encoded at the S-locus and uncover the mechanisms of self/nonself-discrimination and pollen rejection. The recognition mechanisms of SI systems differ between these families. A self-recognition system is adopted by Brassicaceae, whereas a collaborative nonself-recognition system is used by Solanaceae. Work by our group and subsequent studies indicate that plants have evolved diverse SI systems.

Comprehensive computational analysis of the SRK-SP11 molecular interaction underlying self-incompatibility in Brassicaceae using improved structure prediction for cysteine-rich proteins

Plants employ self-incompatibility (SI) to promote cross-fertilization. In Brassicaceae, this process is regulated by the formation of a complex between the pistil determinant S receptor kinase (SRK) and the pollen determinant S-locus protein 11 (SP11, also known as S-locus cysteine-rich protein, SCR). In our previous study, we used the crystal structures of two eSRK-SP11 complexes in Brassica rapa S8 and S9 haplotypes and nine computationally predicted complex models to demonstrate that only the SRK ectodomain (eSRK) and SP11 pairs derived from the same S haplotype exhibit high binding free energy. However, predicting the eSRK-SP11 complex structures for the other 100 + S haplotypes and genera remains difficult because of SP11 polymorphism in sequence and structure. Although protein structure prediction using AlphaFold2 exhibits considerably high accuracy for most protein monomers and complexes, 46% of the predicted SP11 structures that we tested showed < 75 mean per-residue confidence score (pLDDT). Here, we demonstrate that the use of curated multiple sequence alignment (MSA) for cysteine-rich proteins significantly improved model accuracy for SP11 and eSRK-SP11 complexes. Additionally, we calculated the binding free energies of the predicted eSRK-SP11 complexes using molecular dynamics (MD) simulations and observed that some Arabidopsis haplotypes formed a binding mode that was critically different from that of B. rapa S8 and S9. Thus, our computational results provide insights into the haplotype-specific eSRK-SP11 binding modes in Brassicaceae at the residue level. The predicted models are freely available at Zenodo, https://doi.org/10.5281/zenodo.8047768.

Mechanisms of prezygotic post-pollination reproductive barriers in plants

Hybridisation between individuals of different species can lead to maladapted or inviable progeny due to genetic incompatibilities between diverging species. On the other hand, mating with close relatives, or self-fertilisation may lead to inbreeding depression. Thus, both too much or too little divergence may lead to problems and the organisms have to carefully choose mating partners to avoid both of these pitfalls. In plants this choice occurs at many stages during reproduction, but pollen-pistil interactions play a particularly important role in avoiding inbreeding and hybridisation with other species. Interestingly, the mechanisms involved in avoidance of selfing and interspecific hybridisation may work via shared molecular pathways, as self-incompatible species tend to be more 'choosy' with heterospecific pollen compared to self-compatible ones. This review discusses various prezygotic post-pollination barriers to interspecific hybridisation, with a focus on the mechanisms of pollen-pistil interactions and their role in the maintenance of species integrity.

Paradigms of receptor kinase signaling in plants

Plant receptor kinases (RKs) function as key plasma-membrane localized receptors in the perception of molecular ligands regulating development and environmental response. Through the perception of diverse ligands, RKs regulate various aspects throughout the plant life cycle from fertilization to seed set. Thirty years of research on plant RKs has generated a wealth of knowledge on how RKs perceive ligands and activate downstream signaling. In the present review, we synthesize this body of knowledge into five central paradigms of plant RK signaling: (1) RKs are encoded by expanded gene families, largely conserved throughout land plant evolution; (2) RKs perceive many different kinds of ligands through a range of ectodomain architectures; (3) RK complexes are typically activated by co-receptor recruitment; (4) post-translational modifications fulfill central roles in both the activation and attenuation of RK-mediated signaling; and, (5) RKs activate a common set of downstream signaling processes through receptor-like cytoplasmic kinases (RLCKs). For each of these paradigms, we discuss key illustrative examples and also highlight known exceptions. We conclude by presenting five critical gaps in our understanding of RK function.

Contrasting self-recognition rejection systems for self-incompatibility in Brassica and Papaver

Self-incompatibility (SI) plays a pivotal role in whether self-pollen is accepted or rejected. Most SI systems employ two tightly linked loci encoding highly polymorphic pollen (male) and pistil (female) S-determinants that control whether self-pollination is successful or not. In recent years our knowledge of the signalling networks and cellular mechanisms involved has improved considerably, providing an important contribution to our understanding of the diverse mechanisms used by plant cells to recognise each other and elicit responses. Here, we compare and contrast two important SI systems employed in the Brassicaceae and Papaveraceae. Both use 'self-recognition' systems, but their genetic control and S-determinants are quite different. We describe the current knowledge about the receptors and ligands, and the downstream signals and responses utilized to prevent self-seed set. What emerges is a common theme involving the initiation of destructive pathways that block the key processes that are required for compatible pollen-pistil interactions.

Plant reproduction: Stigma receptors regulate reactive oxygen species to establish pollination barriers

Exciting new research highlights how stigmatic receptors purposed for recognizing self-incompatible pollen interact with the FERONIA pathway to regulate stigmatic reactive oxygen species production to enforce a barrier against self-, intra- and interspecific pollen.

MLPK function is not required for self-incompatibility in the S29 haplotype of Brassica rapa L

S²⁹ haplotype does not require the MLPK function for self-incompatibility in Brassica rapa. Self-incompatibility (SI) in Brassicaceae is regulated by the self-recognition mechanism, which is based on the S-haplotype-specific direct interaction of the pollen-derived ligand, SP11/SCR, and the stigma-side receptor, SRK. M locus protein kinase (MLPK) is known to be one of the positive effectors of the SI response. MLPK directly interacts with SRK, and is phosphorylated by SRK in Brassica rapa. In Brassicaceae, MLPK was demonstrated to be essential for SI in B. rapa and Brassica napus, whereas it is not essential for SI in Arabidopsis thaliana (with introduced SRK and SP11/SCR from related SI species). Little is known about what determines the need for MLPK in SI of Brassicaceae. In this study, we investigated the relationship between S-haplotype diversity and MLPK function by analyzing the SI phenotypes of different S haplotypes in a mlpk/mlpk mutant background. The results have clarified that in B. rapa, all the S haplotypes except the S²⁹ we tested need the MLPK function, but the S²⁹ haplotype does not require MLPK for the SI. Comparative analysis of MLPK-dependent and MLPK-independent S haplotype might provide new insight into the evolution of S-haplotype diversity and the molecular mechanism of SI in Brassicaceae.

Generation of Arabidopsis thaliana transformants showing the self-recognition activity of Brassica rapa

Self-incompatibility in the Brassicaceae family is governed by SRK and SCR, which are two highly polymorphic genes located at the S-locus. Previously, the Arabidopsis lyrata SRK and SCR genes were introduced into Arabidopsis thaliana to generate self-incompatible lines. However, there are no reports showing that Brassica SRK and SCR genes confer self-incompatibility in A. thaliana. Doing so would further advance the mechanistic understanding of self-incompatibility in Brassicaceae. Therefore, we attempted to generate A. thaliana transformants showing the self-recognition activity of Brassica rapa by introducing BrSCR along with a chimeric BrSRK (BrSRK chimera, in which the kinase domain of BrSRK was replaced with that of AlSKR-b). We found that the BrSRK chimera and BrSCR of B. rapa S-9 and S-46 haplotypes, but not those of S-29, S-44, and S-60 haplotypes, conferred self-recognition activity in A. thaliana. Analyses of A. thaliana transformants expressing mutant variants of the BrSRK-9 chimera and BrSCR-9 revealed that mutations at the amino acid residues involved in BrSRK9-BrSCR9 interaction caused defects in the self-incompatibility response. The method developed in this study for generating self-incompatible A. thaliana transformants showing B. rapa self-recognition activity will be useful for analysis of self-recognition mechanisms in Brassicaceae.

Cell-cell signaling during the Brassicaceae self-incompatibility response

Self-incompatibility (SI) is a mechanism that many plant families employ to prevent self-fertilization. In the Brassicaceae, the S-haplotype-specific interaction of the pollen-borne ligand, and a stigma-specific receptor protein kinase triggers a signaling cascade that culminates in the rejection of self-pollen. While the upstream molecular components at the receptor level of the signaling pathway have been extensively studied, the intracellular responses beyond receptor activation were not as well understood. Recent research has uncovered several key molecules and signaling events that operate in concert for the manifestation of the self-incompatible responses in Brassicaceae stigmas. Here, we review the recent discoveries in both the compatible and self-incompatible pathways and provide new perspectives on the early stages of Brassicaceae pollen-pistil interactions.

Genetic and Molecular Characterization of a Self-Compatible Brassica rapa Line Possessing a New Class II S Haplotype

Most flowering plants have evolved a self-incompatibility (SI) system to maintain genetic diversity by preventing self-pollination. The Brassica species possesses sporophytic self-incompatibility (SSI), which is controlled by the pollen- and stigma-determinant factors SP11/SCR and SRK. However, the mysterious molecular mechanism of SI remains largely unknown. Here, a new class II S haplotype, named BrS-325, was identified in a pak choi line '325', which was responsible for the completely self-compatible phenotype. To obtain the entire S locus sequences, a complete pak choi genome was gained through Nanopore sequencing and de novo assembly, which provided a good reference genome for breeding and molecular research in B. rapa. S locus comparative analysis showed that the closest relatives to BrS-325 was BrS-60, and high sequence polymorphism existed in the S locus. Meanwhile, two duplicated SRKs (BrSRK-325a and BrSRK-325b) were distributed in the BrS-325 locus with opposite transcription directions. BrSRK-325b and BrSCR-325 were expressed normally at the transcriptional level. The multiple sequence alignment of SCRs and SRKs in class II S haplotypes showed that a number of amino acid variations were present in the contact regions (CR II and CR III) of BrSCR-325 and the hypervariable regions (HV I and HV II) of BrSRK-325s, which may influence the binding and interaction between the ligand and the receptor. Thus, these results suggested that amino acid variations in contact sites may lead to the SI destruction of a new class II S haplotype BrS-325 in B. rapa. The complete SC phenotype of '325' showed the potential for practical breeding application value in B. rapa.

Genetic Diversity of Genes Controlling Unilateral Incompatibility in Japanese Cultivars of Chinese Cabbage

In recent years, unilateral incompatibility (UI), which is an incompatibility system for recognizing and rejecting foreign pollen that operates in one direction, has been shown to be closely related to self-incompatibility (SI) in Brassica rapa. The stigma- and pollen-side recognition factors (SUI1 and PUI1, respectively) of this UI are similar to those of SI (stigma-side SRK and pollen-side SP11), indicating that SUI1 and PUI1 interact with each other and cause pollen-pistil incompatibility only when a specific genotype is pollinated. To clarify the genetic diversity of SUI1 and PUI1 in Japanese B. rapa, here we investigated the UI phenotype and the SUI1/PUI1 sequences in Japanese commercial varieties of Chinese cabbage. The present study showed that multiple copies of nonfunctional PUI1 were located within and in the vicinity of the UI locus region, and that the functional SUI1 was highly conserved in Chinese cabbage. In addition, we found a novel nonfunctional SUI1 allele with a dominant negative effect on the functional SUI1 allele in the heterozygote.

Molecular mechanisms of plant peptide binding to receptors

Plants have evolved diverse peptide hormones and cognate receptors to orchestrate plant growth and development. Secreted peptide ligands are mainly sensed by membrane receptor kinases that mediate cell-cell communication. The secreted peptides are categorized into two groups: small linear post-translationally modified peptides and cysteine-rich peptides. The small linear peptides are recognized by the corresponding receptors and co-receptors in a conserved manner. By contrast, the cysteine-rich peptides are perceived by various types of receptor proteins using diverse binding modes. Recent studies have revealed the molecular and mechanistic origins of peptide recognition and receptor activation. This review summarizes plant-peptide binding modes and receptor-activation mechanisms that have been structurally characterized in recent studies.

The role of diverse LURE-type cysteine-rich peptides as signaling molecules in plant reproduction

In angiosperm sexual reproduction, the male pollen tube undergoes a series of interactions with female tissues. For efficient growth and precise guidance, the pollen tube perceives extracellular ligands. In recent decades, various types of secreted cysteine-rich peptides (CRPs) have been identified as peptide ligands that regulate diverse angiosperm reproduction processes, including pollen tube germination, growth, guidance, and rupture. Notably, in two distant core eudicot plants, multiple LURE-type CRPs were found to be secreted from egg-accompanying synergid cells, and these CRPs act as a cocktail of pollen tube attractants for the final step of pollen tube guidance. LURE-type CRPs have species-preferential activity, even among close relatives, and exhibit remarkably divergent molecular evolution with conserved cysteine frameworks, demonstrating that they play a key role in species recognition in pollen tube guidance. In this review, I focus on "reproductive CRPs," particularly LURE-type CRPs, which underlie common but species-specific mechanisms in angiosperm sexual reproduction, and discuss their action, functional regulation, receptors, and evolution.

Spatiogenetic characterization of S receptor kinase (SRK) alleles in naturalized populations of Raphanus sativus L. var. raphanistroides on Yakushima island

In various coastal areas of Japan, naturalized radish populations are observed. Radish is a cruciferous plant and exhibits self-incompatibility, involving a system controlled by a single locus with multiple S alleles. Although the S allele diversity of radish cultivars and wild radishes has been characterized, the S allele distribution in naturalized populations has not yet been analyzed in relation to the positions of the plants in situ. Here, we show the S allele distribution in naturalized radish populations of Yakushima, a small island in the East China Sea, with positions of the plants. Radish plants were sampled in coastal areas in Yakushima, and their S alleles were detected and characterized. Most of the S alleles had been previously identified in radish cultivars. However, four novel S alleles, which may be unique to Yakushima, were also found. Moreover, seeds in siliques from plants growing in the study areas were sampled, and S allele determination in DNA extracted from these seeds suggested that the plants had exchanged their pollen among their close neighbors. There was also a problem in that the PCR amplification of some SRK alleles was difficult because of their sequence diversity in the naturalized populations, as occurs in cultivars. Our results suggest that the exchange of S alleles between cultivars and naturalized populations occurs and that S alleles in naturalized populations are highly diverse. The methodology established in our study should be applicable to other self-incompatible species to dissect the diversity of S allele distribution in naturalized populations.

The RGF/GLV/CLEL Family of Short Peptides Evolved Through Lineage-Specific Losses and Diversification and Yet Conserves Its Signaling Role Between Vascular Plants and Bryophytes

Short secreted plant peptides act as key signaling molecules and control a plethora of developmental and physiological processes. The ROOT GROWTH FACTOR (RGF)/GOLVEN (GLV)/CLE-Like (CLEL) family of peptides was discovered to be involved in root development in Arabidopsis thaliana. In contrast to active research efforts, which have been revealing receptors and downstream signaling components, little attention has been paid to evolutionary processes that shaped the RGF signaling system as we know it in angiosperms today. As a first step toward understanding how RGF signaling emerged and evolved, this study aimed to elucidate the phylogenetic distribution and functional conservation of RGF-like sequences. Using publicly available, genome and transcriptome data, RGF-like sequences were searched in 27 liverworts, 22 mosses, 8 hornworts, 23 lycophytes, 23 ferns, 38 gymnosperms, and 8 angiosperms. This led to the identification of more than four hundreds of RGF-like sequences in all major extant land plant lineages except for hornworts. Sequence comparisons within and between taxonomic groups identified lineage-specific characters. Notably, one of the two major RGF subgroups, represented by A. thaliana RGF6/GLV1/CLEL6, was found only in vascular plants. This subgroup, therefore, likely emerged in a common ancestor of vascular plants after its divergence from bryophytes. In bryophytes, our results infer independent losses of RGF-like sequences in mosses and hornworts. On the other hand, a single, highly similar RGF-like sequence is conserved in liverworts, including Marchantia polymorpha, a genetically tractable model species. When constitutively expressed, the M. polymorpha RGF-like sequence (MpRGF) affected plant development and growth both in A. thaliana and M. polymorpha. This suggests that MpRGF can exert known RGF-like effects and that MpRGF is under transcriptional control so that its potent activities are precisely controlled. These data suggest that RGFs are conserved as signaling molecules in both vascular plants and bryophytes and that lineage-specific diversification has increased sequence variations of RGFs. All together, our findings form a basis for further studies into RGF peptides and their receptors, which will contribute to our understandings of how peptide signaling pathways evolve.