A De Novo Transcriptome Assembly of Ceratopteris richardii Provides Insights into the Evolutionary Dynamics of Complex Gene Families in Land Plants

A conserved GRAS-domain transcriptional regulator links meristem indeterminacy to sex determination in Ceratopteris gametophytes

Most land plants alternate between generations of sexual gametophytes and asexual sporophytes. Unlike seed plants, fern gametophytes are free living and grow independently of their sporophytes. In homosporous ferns such as Ceratopteris, gametophytes derived from genetically identical spores exhibit sexual dimorphism, developing as either males or hermaphrodites. Males lack meristems and promote cell differentiation into sperm-producing antheridia. In contrast, hermaphrodites initiate multicellular meristems that stay undifferentiated, sustain cell division and prothallus expansion, and drive the formation of egg-producing archegonia. Once initiating the meristem, hermaphrodites secrete the pheromone antheridiogen, which triggers neighboring slower-growing gametophytes to develop as males, while the hermaphrodites themselves remain insensitive to antheridiogen. This strategy promotes outcrossing and prevents all individuals in the colony from becoming males. This study reveals that an evolutionarily conserved GRAS-domain transcriptional regulator (CrHAM), directly repressed by Ceratopteris microRNA171 (CrmiR171), promotes meristem development in Ceratopteris gametophytes and determines the male-to-hermaphrodite ratio in the colony. CrHAM preferentially accumulates within the meristems of hermaphrodites but is excluded from differentiated antheridia. CrHAM sustains meristem proliferation and cell division through conserved hormone pathways. In the meantime, CrHAM inhibits the antheridiogen-induced conversion of hermaphrodites to males by suppressing the male program expression and preventing meristem cells from differentiating into sperm-producing antheridia. This finding establishes a connection between meristem indeterminacy and sex determination in ferns, suggesting both conserved and diversified roles of meristem regulators in land plants.

Conserved autophagy and diverse cell wall composition: unifying features of vascular tissues in evolutionarily distinct plants

BACKGROUND AND AIMS

The formation of multifunctional vascular tissues represents a significant advancement in plant evolution. Differentiation of conductive cells is specific, involving two main pathways, namely protoplast clearance and cell wall modification. In xylogenesis, autophagy is a crucial process for complete protoplast elimination in tracheary elements, whose cell wall also undergoes strong changes. Knowledge pertaining to living sieve elements, which lose most of their protoplast during phloemogenesis, remains limited. We hypothesized that autophagy plays a crucial role, not only in complete cytoplasmic clearance in xylem but also in partial degradation in phloem. Cell wall elaborations of mature sieve elements are not so extensive. These analyses performed on evolutionarily diverse model species potentially make it possible to understand phloemogenesis to an equal extent to xylogenesis.

METHODS

We investigated the distribution of ATG8 protein, which is an autophagy marker, and cell wall components in the roots of ferns, gymnosperms and angiosperms (monocots, dicot herbaceous plants and trees). Furthermore, we conducted a bioinformatic analysis of complete data on ATG8 isoforms for Ceratopteris richardii.

KEY RESULTS

The presence of ATG8 protein was confirmed in both tracheary elements and sieve elements; however, the composition of cell wall components varied considerably among vascular tissues in the selected plants. Arabinogalactan proteins and β-1,4-galactan were detected in the roots of all studied species, suggesting their potential importance in phloem formation or function. In contrast, no evolutionary pattern was observed for xyloglucan, arabinan or homogalacturonan.

CONCLUSIONS

Our findings indicate that the involvement of autophagy in plants is universal during the development of tracheary elements that are dead at maturity and sieve elements that remain alive. Given the conserved nature of autophagy and its function in protoplast degradation for uninterrupted flow, autophagy might have played a vital role in the development of increasingly complex biological organizations, including the formation of vascular tissues. However, different cell wall compositions of xylem and phloem in different species might indicate diverse functionality and potential for substance transport, which is crucial in plant evolution.

Differential expression and co-localization of transcriptional factors during callus transition to differentiation for shoot organogenesis in the water fern Ceratopteris richardii

BACKGROUND AND AIMS

In flowering plants, regeneration can be achieved by a variety of approaches, and different sets of transcriptional factors are involved in these processes. However, regeneration in taxa other than flowering plants remains a mystery. Ceratopteris richardii is a representative fern capable of both direct and indirect organogenesis, and we aimed to investigate the genetics underlying the transition from callus proliferation to differentiation.

METHODS

Morphological and histological analyses were used to determine the type of regeneration involved. RNA sequencing and differential gene expression were used to investigate how the callus switches from proliferation to differentiation. Phylogenetic analysis and RNA in situ hybridization were used to understand whether transcriptional factors are involved in this transition.

KEY RESULTS

The callus formed on nascent leaves and subsequently developed the shoot pro-meristem and shoot meristem, thus completing indirect de novo shoot organogenesis in C. richardii. Genes were differentially expressed during the callus transition from proliferation to differentiation, indicating a role for photosynthesis, stimulus response and transmembrane signalling in this transition and the involvement of almost all cell layers that make up the callus. Transcriptional factors were either downregulated or upregulated, which were generally in many-to-many orthology with genes known to be involved in callus development in flowering plants, suggesting that the genetics of fern callus development are both conserved and divergent. Among them, an STM-like, a PLT-like and an ethylene- and salt-inducible ERF gene3-like gene were expressed simultaneously in the vasculature but not in the other parts of the callus, indicating that the vasculature played a role in the callus transition from proliferation to differentiation.

CONCLUSIONS

Indirect de novo shoot organogenesis could occur in ferns, and the callus transition from proliferation to differentiation required physiological changes, differential expression of transcriptional factors and involvement of the vasculature.

Whole Genome Duplication Events Likely Contributed to the Aquatic Adaptive Evolution of Parkerioideae

As the only aquatic lineage of Pteridaceae, Parkerioideae is distinct from many xeric-adapted species of the family and consists of the freshwater Ceratopteris species and the only mangrove ferns from the genus Acrostichum. Previous studies have shown that whole genome duplication (WGD) has occurred in Parkerioideae at least once and may have played a role in their adaptive evolution; however, more in-depth research regarding this is still required. In this study, comparative and evolutionary transcriptomics analyses were carried out to identify WGDs and explore their roles in the environmental adaptation of Parkerioideae. Three putative WGD events were identified within Parkerioideae, two of which were specific to Ceratopteris and Acrostichum, respectively. The functional enrichment analysis indicated that the lineage-specific WGD events have played a role in the adaptation of Parkerioideae to the low oxygen concentrations of aquatic habitats, as well as different aquatic environments of Ceratopteris and Acrostichum, such as the adaptation of Ceratopteris to reduced light levels and the adaptation of Acrostichum to high salinity. Positive selection analysis further provided evidence that the putative WGD events may have facilitated the adaptation of Parkerioideae to changes in habitat. Moreover, the gene family analysis indicated that the plasma membrane H+-ATPase (AHA), vacuolar H+-ATPase (VHA), and suppressor of K+ transport growth defect 1 (SKD1) may have been involved in the high salinity adaptation of Acrostichum. Our study provides new insights into the evolution and adaptations of Parkerioideae in different aquatic environments.

Comparison between the Gametophyte and the Sporophyte Transcriptomes of the Endangered Fern Vandenboschia speciosa

Genomic resources are essential to understanding the evolution and functional biology of organisms. Nevertheless, generating genomic resources from endangered species may be challenging due to the scarcity of available specimens and sampling difficulties. In this study, we compare the transcriptomes of the sporophyte and the gametophyte of the endangered fern Vandenboschia speciosa. After Illumina sequencing and de novo transcriptome assembly of the gametophyte, annotation proved the existence of cross-species contamination in the gametophyte sample. Thus, we developed an in silico decontamination step for the gametophyte sequences. Once the quality check of the decontaminated reads passed, we produced a de novo assembly with the decontaminated gametophyte reads (with 43,139 contigs) and another combining the sporophyte and in silico decontaminated gametophyte reads (with 42,918 contigs). A comparison of the enriched GO terms from the top 1000 most expressed transcripts from both tissues showed that the gametophyte GO term set was enriched in sequences involved in development, response to stress, and plastid organization, while the sporophyte GO term set had a larger representation of more general metabolic functions. This study complements the available genomic resources on the life cycle of the endangered fern Vandenboschia speciosa.

Cell Division and Meristem Dynamics in Fern Gametophytes

One of the most important questions in all multicellular organisms is how to define and maintain different cell fates during continuous cell division and proliferation. Plant meristems provide a unique research system to address this fundamental question because meristems dynamically maintain themselves and sustain organogenesis through balancing cell division and cell differentiation. Different from the gametophytes of seed plants that depend on their sporophytes and lack meristems, the gametophytes of seed-free ferns develop different types of meristems (including apical cell-based meristems and multicellular apical and marginal meristems) to promote independent growth and proliferation during the sexual gametophyte phase. Recent studies combining confocal time-lapse imaging and computational image analysis reveal the cellular basis of the initiation and proliferation of different types of meristems in fern gametophytes, providing new insights into the evolution of meristems in land plants. In this review, we summarize the recent progress in understanding the cell growth dynamics in fern gametophytes and discuss both conserved and diversified mechanisms underlying meristem cell proliferation in seed-free vascular plants.

Transcriptional analysis of Ceratopteris richardii young sporophyte reveals conservation of stem cell factors in the root apical meristem

Gene expression in roots has been assessed in different plant species in studies ranging from complete organs to specific cell layers, and more recently at the single cell level. While certain genes or functional categories are expressed in the root of all or most plant species, lineage-specific genes have also been discovered. An increasing amount of transcriptomic data is available for angiosperms, while a limited amount of data is available for ferns, and few studies have focused on fern roots. Here, we present a de novo transcriptome assembly from three different parts of the Ceratopteris richardii young sporophyte. Differential gene expression analysis of the root tip transcriptional program showed an enrichment of functional categories related to histogenesis and cell division, indicating an active apical meristem. Analysis of a diverse set of orthologous genes revealed conserved expression in the root meristem, suggesting a preserved role for different developmental roles in this tissue, including stem cell maintenance. The reconstruction of evolutionary trajectories for ground tissue specification genes suggests a high degree of conservation in vascular plants, but not for genes involved in root cap development, showing that certain genes are absent in Ceratopteris or have intricate evolutionary paths difficult to track. Overall, our results suggest different processes of conservation and divergence of genes involved in root development.

Positional cues and cell division dynamics drive meristem development and archegonium formation in Ceratopteris gametophytes

Fern gametophytes are autotrophic and independent of sporophytes, and they develop pluripotent meristems that drive prothallus development and sexual reproduction. To reveal cellular dynamics during meristem development in fern gametophytes, we performed long-term time-lapse imaging and determined the real-time lineage, identity and division activity of each single cell from meristem initiation to establishment in gametophytes of the fern Ceratopteris richardii. Our results demonstrate that in Ceratopteris gametophytes, only a few cell lineages originated from the marginal layer contribute to meristem initiation and proliferation, and the meristem lacks a distinguishable central zone or apical cell with low division activity. Within the meristem, cell division is independent of cell lineages and cells at the marginal layer are more actively dividing than inner cells. Furthermore, the meristem triggers differentiation of adjacent cells into egg-producing archegonia in a position-dependent manner. These findings advance the understanding of diversified meristem and gametophyte development in land plants.

Time-lapse imaging of the fern Ceratopteris richardii during meristem initiation and proliferation provides insights into the lineage, identity and division activity of each cell throughout the growth of gametophytes.

Cell growth dynamics in two types of apical meristems in fern gametophytes

In contrast to seed plants, the gametophytes of seed-free plants develop pluripotent meristems, which promote and sustain their independent growth and development. To date, the cellular basis of meristem development in gametophytes of seed-free ferns remains largely unknown. In this study, we used Woodsia obtusa, the blunt-lobe cliff fern, to quantitatively determine cell growth dynamics in two different types of apical meristems - the apical initial centered meristem and the multicellular apical meristem in gametophytes. Through confocal time-lapse live imaging and computational image analysis and quantification, we determined unique patterns of cell division and growth that sustain or terminate apical initials, dictate the transition from apical initials to multicellular apical meristems, and drive proliferation of apical meristems in ferns. Quantitative results showed that small cells correlated to active cell division in fern gametophytes. The marginal cells of multicellular apical meristems in fern gametophytes undergo division in both anticlinal and periclinal orientations, not only increasing cell numbers but also playing a dominant role in increasing cell layers during gametophyte development. All these findings provide insights into the function and regulation of meristems in gametophytes of seed-free vascular plants, suggesting both conserved and diversified mechanisms underlying meristem cell proliferation across land plants.

The Land-Sea Connection: Insights Into the Plant Lineage from a Green Algal Perspective

The colonization of land by plants generated opportunities for the rise of new heterotrophic life forms, including humankind. A unique event underpinned this massive change to earth ecosystems-the advent of eukaryotic green algae. Today, an abundant marine green algal group, the prasinophytes, alongside prasinodermophytes and nonmarine chlorophyte algae, is facilitating insights into plant developments. Genome-level data allow identification of conserved proteins and protein families with extensive modifications, losses, or gains and expansion patterns that connect to niche specialization and diversification. Here, we contextualize attributes according to Viridiplantae evolutionary relationships, starting with orthologous protein families, and then focusing on key elements with marked differentiation, resulting in patchy distributions across green algae and plants. We place attention on peptidoglycan biosynthesis, important for plastid division and walls; phytochrome photosensors that are master regulators in plants; and carbohydrate-active enzymes, essential to all manner of carbohydratebiotransformations. Together with advances in algal model systems, these areas are ripe for discovering molecular roles and innovations within and across plant and algal lineages.

The biology of C. richardii as a tool to understand plant evolution

The fern Ceratopteris richardii has been studied as a model organism for over 50 years because it is easy to grow and has a short life cycle. In particular, as the first homosporous vascular plant for which genomic resources were developed, C. richardii has been an important system for studying plant evolution. However, we know relatively little about the natural history of C. richardii. In this article, we summarize what is known about this aspect of C. richardii, and discuss how learning more about its natural history could greatly increase our understanding of the evolution of land plants.

Function and Regulation of microRNA171 in Plant Stem Cell Homeostasis and Developmental Programing

MicroRNA171 (miR171), a group of 21-nucleotide single-strand small RNAs, is one ancient and conserved microRNA family in land plants. This review focuses on the recent progress in understanding the role of miR171 in plant stem cell homeostasis and developmental patterning, and the regulation of miR171 by developmental cues and environmental signals. Specifically, miR171 regulates shoot meristem activity and phase transition through repressing the HAIRYMERISTEM (HAM) family genes. In the model species Arabidopsis, miR171 serves as a short-range mobile signal, which initiates in the epidermal layer of shoot meristems and moves downwards within a limited distance, to pattern the apical-basal polarity of gene expression and drive stem cell dynamics. miR171 levels are regulated by light and various abiotic stresses, suggesting miR171 may serve as a linkage between environmental factors and cell fate decisions. Furthermore, miR171 family members also demonstrate both conserved and lineage-specific functions in land plants, which are summarized and discussed here.

Quantitative live-imaging reveals the dynamics of apical cells during gametophyte development in ferns

Meristems in land plants share conserved functions but develop highly variable structures. Meristems in seed-free plants, including ferns, usually contain one or a few pyramid-/wedge-shaped apical cells (ACs) as initials, which are lacking in seed plants. It remained unclear how ACs promote cell proliferation in fern gametophytes and whether any persistent AC exists to sustain fern gametophyte development continuously. Here, we uncovered previously undefined ACs maintained even at late developmental stages in fern gametophytes. Through quantitative live-imaging, we determined division patterns and growth dynamics that maintain the persistent AC in Sphenomeris chinensis, a representative fern. The AC and its immediate progenies form a conserved cell packet, driving cell proliferation and prothallus expansion. At the apical centre of gametophytes, the AC and its adjacent progenies display small dimensions resulting from active cell division instead of reduced cell expansion. These findings provide insight into diversified meristem development in land plants.

Timing of meristem initiation and maintenance determines the morphology of fern gametophytes

The alternation of generations in land plants occurs between the sporophyte phase and the gametophyte phase. The sporophytes of seed plants develop self-maintained, multicellular meristems, and these meristems determine plant architecture. The gametophytes of seed plants lack meristems and are heterotrophic. In contrast, the gametophytes of seed-free vascular plants, including ferns, are autotrophic and free-living, developing meristems to sustain their independent growth and proliferation. Compared with meristems in the sporophytes of seed plants, the cellular mechanisms underlying meristem development in fern gametophytes remain largely unknown. Here, using confocal time-lapse live imaging and computational segmentation and quantification, we determined different patterns of cell divisions associated with the initiation and proliferation of two distinct types of meristems in gametophytes of two closely related Pteridaceae ferns, Pteris vittata and Ceratopteris richardii. Our results reveal how the simple timing of a switch between two meristems has considerable consequences for the divergent gametophyte morphologies of the two ferns. They further provide evolutionary insight into the function and regulation of gametophyte meristems in seed-free vascular plants.

Conditional stomatal closure in a fern shares molecular features with flowering plant active stomatal responses

Stomata evolved as plants transitioned from water to land, enabling carbon dioxide uptake and water loss to be controlled. In flowering plants, the most recently divergent land plant lineage, stomatal pores actively close in response to drought. In this response, the phytohormone abscisic acid (ABA) triggers signaling cascades that lead to ion and water loss in the guard cells of the stomatal complex, causing a reduction in turgor and pore closure. Whether this stimulus-response coupling pathway acts in other major land plant lineages is unclear, with some investigations reporting that stomatal closure involves ABA but others concluding that closure is passive. Here, we show that in the model fern Ceratopteris richardii active stomatal closure is conditional on sensitization by pre-exposure to either low humidity or exogenous ABA and is promoted by ABA. RNA-seq analysis and de novo transcriptome assembly reconstructed the protein-coding complement of the C. richardii genome, with coverage comparable to other plant models, enabling transcriptional signatures of stomatal sensitization and closure to be inferred. In both cases, changes in abundance of homologs of ABA, Ca²⁺, and ROS-related signaling components were observed, suggesting that the closure-response pathway is conserved in ferns and flowering plants. These signatures further suggested that sensitization is achieved by lowering the threshold required for a subsequent closure-inducing signal to trigger a response. We conclude that the canonical signaling network for active stomatal closure functioned in at least a rudimentary form in the stomata of the last common ancestor of ferns and flowering plants.

N-terminal region is required for functions of the HAM family member

Shoot meristems contain stem cells, and they sustain growth and development of the above-ground tissues in land plants. The HAIRY MERISTEM (HAM) family genes, encoding GRAS-domain transcriptional regulators, play essential roles in the control of shoot meristem development and stem cell homeostasis in several flowering plants. Similar to other GRAS proteins, the C-terminal regions of HAM family proteins across land plants are conserved, containing signature motifs that define the GRAS domain. In contrast, the N-terminal regions of HAM family proteins display substantial divergence in sequence and length. Whether the variable and divergent N-termini are required for the conserved functions of HAM proteins is unknown. Our recent work showed that CrHAM - the HAM homolog in the fern Ceratopteris richardii was able to replace the role of type-II HAM genes in Arabidopsis, maintaining established shoot apical meristems and promoting the initiation of new stem cell niches. Here, we provide additional information and show that CrHAM contains a much longer N-terminal region compared to Arabidopsis HAM proteins, which is conserved among different fern HAM homologs. The deletion of this region largely compromises the ability of CrHAM to replace the function of Arabidopsis HAM proteins in shoot meristems. These new data together with previous results suggest that, although lacking the sequence conservation among HAM homologs from different plant lineages, the N-termini are important for the conserved functions of HAM family proteins.

De novo transcriptome assembly, polymorphic SSR markers development and population genetics analyses for southern corn rust (Puccinia polysora)

Southern corn rust is a destructive maize disease caused by Puccinia polysora Underw that can lead to severe yield losses. However, genomic information and microsatellite markers are currently unavailable for this disease. In this study, we generated a total of 27,295,216 high-quality cDNA sequence reads using Illumina sequencing technology. These reads were assembled into 17,496 unigenes with an average length of 1015 bp. The functional annotation indicated that 8113 (46.37%), 1933 (11.04%) and 5516 (31.52%) unigenes showed significant similarity to known proteins in the NCBI Nr, Nt and Swiss-Prot databases, respectively. In addition, 2921 (16.70%) unigenes were assigned to KEGG database categories; 4218 (24.11%), to KOG database categories; and 6,603 (37.74%), to GO database categories. Furthermore, we identified 8,798 potential SSRs among 6653 unigenes. A total of 9 polymorphic SSR markers were developed to evaluate the genetic diversity and population structure of 96 isolates collected from Guangdong Province in China. Clonal reproduction of P. polysora in Guangdong was dominant. The YJ (Yangjiang) population had the highest genotypic diversity and the greatest number of the multilocus genotypes, followed by the HY (Heyuan), HZ (Huizhou) and XY (Xinyi) populations. These results provide valuable information for the molecular genetic analysis of P. polysora and related species.

De Novo Sporophyte Transcriptome Assembly and Functional Annotation in the Endangered Fern Species Vandenboschia speciosa (Willd.) G. Kunkel

We sequenced the sporophyte transcriptome of Killarney fern (Vandenboschia speciosa (Willd.) G. Kunkel). In addition to being a rare endangered Macaronesian-European endemism, this species has a huge genome (10.52 Gb) as well as particular biological features and extreme ecological requirements. These characteristics, together with the systematic position of ferns among vascular plants, make it of high interest for evolutionary, conservation and functional genomics studies. The transcriptome was constructed de novo and contained 36,430 transcripts, of which 17,706 had valid BLAST hits. A total of 19,539 transcripts showed at least one of the 7362 GO terms assigned to the transcriptome, whereas 6547 transcripts showed at least one of the 1359 KEGG assigned terms. A prospective analysis of functional annotation results provided relevant insights on genes involved in important functions such as growth and development as well as physiological adaptations. In this context, a catalogue of genes involved in the genetic control of plant development, during the vegetative to reproductive transition, in stress response as well as genes coding for transcription factors is given. Altogether, this study provides a first step towards understanding the gene expression of a significant fern species and the in silico functional and comparative analyses reported here provide important data and insights for further comparative evolutionary studies in ferns and land plants in general.

HAM Gene Family and Shoot Meristem Development

Land plants develop highly diversified shoot architectures, all of which are derived from the pluripotent stem cells in shoot apical meristems (SAMs). As sustainable resources for continuous organ formation in the aboveground tissues, SAMs play an important role in determining plant yield and biomass production. In this review, we summarize recent advances in understanding one group of key regulators - the HAIRY MERISTEM (HAM) family GRAS domain proteins - in shoot meristems. We highlight the functions of HAM family members in dictating shoot stem cell initiation and proliferation, the signaling cascade that shapes HAM expression domains in shoot meristems, and the conservation and diversification of HAM family members in land plants. We also discuss future directions that potentially lead to a more comprehensive view of the HAM gene family and stem cell homeostasis in land plants.

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.