The evolutionary history of plant T2/S-type ribonucleases

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

PREMISE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

The origin and the genetic regulation of the self-compatibility mechanism in clementine (Citrus clementina Hort. ex Tan.)

Self-incompatibility (SI) is a genetic mechanism common in flowering plants to prevent self-fertilization. Among citrus species, several pummelo, mandarin, and mandarin-like accessions show SI behavior. In these species, SI is coupled with a variable degree of parthenocarpy ensuring the production of seedless fruits, a trait that is highly appreciated by consumers. In Citrus, recent evidences have shown the presence of a gametophytic SI system based on S-ribonucleases (S-RNases) ability to impair self-pollen tube growth in the upper/middle part of the style. In the present study, we combined PCR analysis and next-generation sequencing technologies, to define the presence of S7- and S11-Rnases in the S-genotype of the Citrus clementina (Hort. ex Tan.), the self-incompatible 'Comune' clementine and its self-compatible natural mutant 'Monreal'. The reference genome of 'Monreal' clementine is presented for the first time, providing more robust results on the genetic sequence of the newly discovered S7-RNase. SNP discovery analysis coupled with the annotation of the variants detected enabled the identification of 7,781 SNPs effecting 5,661 genes in 'Monreal' compared to the reference genome of C. clementina. Transcriptome analysis of unpollinated pistils at the mature stage from both clementine genotypes revealed the lack of expression of S7-RNase in 'Monreal' suggesting its involvement in the loss of the SI response. RNA-seq analysis followed by gene ontology studies enabled the identification of 2,680 differentially expressed genes (DEGs), a significant number of those is involved in oxidoreductase and transmembrane transport activity. Merging of DNA sequencing and RNA data led to the identification of 164 DEGs characterized by the presence of at least one SNP predicted to induce mutations with a high effect on their amino acid sequence. Among them, four candidate genes referring to two Agamous-like MADS-box proteins, to MYB111 and to MLO-like protein 12 were validated. Moreover, the transcription factor MYB111 appeared to contain a binding site for the 2.0-kb upstream sequences of the S7- and S11-RNase genes. These results provide useful information about the genetic bases of SI indicating that SNPs present in their sequence could be responsible for the differential expression and the regulation of S7-RNase and consequently of the SI mechanism.

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.

Rapid evolution of T2/S-RNase genes in Fragaria linked to multiple transitions from self-incompatibility to self-compatibility

The T2/RNase gene family is widespread in eukaryotes, and particular members of this family play critical roles in the gametophytic self-incompatibility (GSI) system in plants. Wild diploid strawberry (Fragaria) species have diversified their sexual systems via self-incompatible and self-compatible traits, yet how these traits evolved in Fragaria remains elusive. By integrating the published and de novo assembled genomes and the newly generated RNA-seq data, members of the RNase T2 gene family were systematically identified in six Fragaria species, including three self-incompatible species (Fragaria nipponica, Fragaria nubicola, and Fragaria viridis) and three self-compatible species (Fragaria nilgerrensis, Fragaria vesca, and Fragaria iinumae). In total, 115 RNase T2 genes were identified in the six Fragaria genomes and can be classified into three classes (I-III) according to phylogenetic analysis. The identified RNase T2 genes could be divided into 22 homologous gene sets according to amino acid sequence similarity and phylogenetic and syntenic relationships. We found that extensive gene loss and pseudogenization coupled with small-scale duplications mainly accounted for variations in the RNase T2 gene numbers in Fragaria. Multiple copies of homologous genes were mainly generated from tandem and segmental duplication events. Furthermore, we newly identified five S-RNase genes in three self-incompatible Fragaria genomes, including two in F. nipponica, two in F. viridis, and one in F. nubicola, which fit for typical features of a pistil determinant, including highly pistil-specific expression, highly polymorphic proteins and alkaline isoelectric point (pI), while no S-RNase genes were found in all three self-compatible Fragaria species. Surprisingly, these T2/S-RNase genes contain at least one large intron (>10 kb). This study revealed that the rapid evolution of T2/S-RNase genes within the Fragaria genus could be associated with its sexual mode, and repeated evolution of the self-compatible traits in Fragaria was convergent via losses of S-RNase.

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.

PURC Provides Improved Sequence Inference for Polyploid Phylogenetics and Other Manifestations of the Multiple-Copy Problem

Inferring the true biological sequences from amplicon mixtures remains a difficult bioinformatics problem. The traditional approach is to cluster sequencing reads by similarity thresholds and treat the consensus sequence of each cluster as an "operational taxonomic unit" (OTU). Recently, this approach has been improved by model-based methods that correct PCR and sequencing errors in order to infer "amplicon sequence variants" (ASVs). To date, ASV approaches have been used primarily in metagenomics, but they are also useful for determining homeologs in polyploid organisms. To facilitate the usage of ASV methods among polyploidy researchers, we incorporated ASV inference alongside OTU clustering in PURC v2.0, a major update to PURC (Pipeline for Untangling Reticulate Complexes). In addition, PURC v2.0 features faster demultiplexing than the original version and has been updated to be compatible with Python 3. In this chapter we present results indicating that using the ASV approach is more likely to infer the correct biological sequences in comparison to the earlier OTU-based PURC and describe how to prepare sequencing data, run PURC v2.0 under several different modes, and interpret the output.

Evolutionary and structural aspects of Solanaceae RNases T2

Plant RNases T2 are involved in several physiological and developmental processes, including inorganic phosphate starvation, senescence, wounding, defense against pathogens, and the self-incompatibility system. Solanaceae RNases form three main clades, one composed exclusively of S-RNases and two that include S-like RNases. We identified several positively selected amino acids located in highly flexible regions of these molecules, mainly close to the B1 and B2 substrate-binding sites in S-like RNases and the hypervariable regions of S-RNases. These differences between S- and S-like RNases in the flexibility of amino acids in substrate-binding regions are essential to understand the RNA-binding process. For example, in the S-like RNase NT, two positively selected amino acid residues (Tyr156 and Asn134) are located at the most flexible sites on the molecular surface. RNase NT is induced in response to tobacco mosaic virus infection; these sites may thus be regions of interaction with pathogen proteins or viral RNA. Differential selective pressures acting on plant ribonucleases have increased amino acid variability and, consequently, structural differences within and among S-like RNases and S-RNases that seem to be essential for these proteins play different functions.

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.

Tissue-specific enhancement of OsRNS1 with root-preferred expression is required for the increase of crop yield

INTRODUCTION

Root development is a fundamental process that supports plant survival and crop productivity. One of the essential factors to consider when developing biotechnology crops is the selection of a promoter that can optimize the spatial-temporal expression of introduced genes. However, there are insufficient cases of suitable promoters in crop plants, including rice.

OBJECTIVES

This study aimed to verify the usefulness of a new rice root-preferred promoter to optimize the function of a target gene with root-preferred expression in rice.

METHODS

osrns1 mutant had defects in root development based on T-DNA insertional mutant screening and CRISPR technology. To optimize the function of OsRNS1, we generated OsRNS1-overexpression plants under two different promoters: a whole-plant expression promoter and a novel root-preferred expression promoter. Root growth, yield-related agronomic traits, RNA-seq, and reactive oxygen species (ROS) accumulation were analyzed for comparison.

RESULTS

OsRNS1 was found to be involved in root development through T-DNA insertional mutant analysis and gene editing mutant analysis. To understand the gain of function of OsRNS1, pUbi1::OsRNS1 was generated for the whole-plant expression, and both root growth defects and overall growth defects were found. To overcome this problem, a root-preferential overexpression line using Os1-CysPrxB promoter (Per) was generated and showed an increase in root length, plant height, and grain yield compared to wild-type (WT). RNA-seq analysis revealed that the response to oxidative stress-related genes was significantly up-regulated in both overexpression lines but was more obvious in pPer::OsRNS1. Furthermore, ROS levels in the roots were drastically decreased in pPer::OsRNS1 but were increased in the osrns1 mutants compared to WT.

CONCLUSION

The results demonstrated that the use of a root-preferred promoter effectively optimizes the function of OsRNS1 and is a useful strategy for improving root-related agronomic traits as well as ROS regulation.

Intraspecific breakdown of self-incompatibility in Physalis acutifolia (Solanaceae)

Variation in mating systems is prevalent throughout angiosperms, with many transitions between outcrossing and selfing above and below the species level. This study documents a new case of an intraspecific breakdown of self-incompatibility in a wild relative of tomatillo, Physalis acutifolia. We used controlled greenhouse crosses to identify self-incompatible (SI) and self-compatible (SC) individuals grown from seed sampled across seven sites across Arizona and New Mexico. We measured 14 flower and fruit traits to test for trait variation associated with mating system. We also quantified pollen tube growth in vivo and tested for the presence of the S-RNase proteins in SI and SC styles. We found that seed from six of the seven sites produced SI individuals that terminated self-pollen tubes in the style and showed detectable S-RNase expression. By contrast, seed from one Arizona site produced SC individuals with no S-RNase expression. These SC individuals displayed typical selfing-syndrome traits such as smaller corollas, reduced stigma-anther distances, and a smaller pollen-ovule ratio. We also found plasticity in self-incompatibility as most of the SI individuals became SC and lost S-RNase expression roughly after 6 months in the greenhouse. While fixed differences in mating systems are known among the SI wild species and the often SC domesticated tomatillos, our study is the first to demonstrate intraspecific variation in natural populations as well as variation in SI over an individual's lifespan.

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

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

S-RNase Alleles Associated With Self-Compatibility in the Tomato Clade: Structure, Origins, and Expression Plasticity

The self-incompatibility (SI) system in the Solanaceae is comprised of cytotoxic pistil S-RNases which are countered by S-locus F-box (SLF) resistance factors found in pollen. Under this barrier-resistance architecture, mating system transitions from SI to self-compatibility (SC) typically result from loss-of-function mutations in genes encoding pistil SI factors such as S-RNase. However, the nature of these mutations is often not well characterized. Here we use a combination of S-RNase sequence analysis, transcript profiling, protein expression and reproductive phenotyping to better understand different mechanisms that result in loss of S-RNase function. Our analysis focuses on 12 S-RNase alleles identified in SC species and populations across the tomato clade. In six cases, the reason for gene dysfunction due to mutations is evident. The six other alleles potentially encode functional S-RNase proteins but are typically transcriptionally silenced. We identified three S-RNase alleles which are transcriptionally silenced under some conditions but actively expressed in others. In one case, expression of the S-RNase is associated with SI. In another case, S-RNase expression does not lead to SI, but instead confers a reproductive barrier against pollen tubes from other tomato species. In the third case, expression of S-RNase does not affect self, interspecific or inter-population reproductive barriers. Our results indicate that S-RNase expression is more dynamic than previously thought, and that changes in expression can impact different reproductive barriers within or between natural populations.

Genome-wide identification and characterization of legume T2 Ribonuclease gene family and analysis of GmaRNS9, a soybean T2 Ribonuclease gene, function in nodulation

T2 ribonuclease family (RNaseT2) proteins are secretory and nonspecific endoribonucleases that have a large conserved biological role. Family members of RNaseT2 are found in every organism and carry out important biological functions. However, little is known about the functions of these proteins in legumes, including potential roles in symbiotic nodulation. This study aimed to characterize and perform bioinformatic analysis of RNaseT2 genes in four legume species that their genome was sequenced. In total, 60 RNaseT2 genes were identified and characterized. By analyzing their phylogeny, we divided these RNaseT2 into five clades. Expression analysis of RNaseT2 genes indicated that these genes are expressed in various tissues, and the most expression level was related to the pod, flower, and root. Moreover, GmaRNS9 expression analysis in soybean was consistent with in silico studies and demonstrated that this gene usually has high root tip expression. GmaRNS9 expression was reduced by Bradyrhizobium japonicum inoculation and nodule formation. Reduced expression of this gene was possibly controlled by the GmNARK gene either directly or pleiotropically through increased phosphorus requirements during increased nodulation. However, the nutrient stress (phosphate and nitrate starvation) led to an increase in the expression level of GmRNS9. In silico and quantitative gene expression analyses showed that RNaseT2 genes could play important roles in the growth and development of legumes as well as nodulation.

The maintenance of polymorphism in an ancient social supergene

Supergenes, regions of the genome with suppressed recombination between sets of functional mutations, contribute to the evolution of complex phenotypes in diverse systems. Excluding sex chromosomes, most supergenes discovered so far appear to be young, being found in one species or a few closely related species. Here, we investigate how a chromosome harbouring an ancient supergene has evolved over about 30 million years (Ma). The Formica supergene underlies variation in colony queen number in at least five species. We expand previous analyses of sequence divergence on this chromosome to encompass about 90 species spanning the Formica phylogeny. Within the nonrecombining region, the gene knockout contains 22 single nucleotide polymorphisms (SNPs) that are consistently differentiated between two alternative supergene haplotypes in divergent European Formica species, and we show that these same SNPs are present in most Formica clades. In these clades, including an early diverging Nearctic Formica clade, individuals with alternative genotypes at knockout also have higher differentiation in other portions of this chromosome. We identify hotspots of SNPs along this chromosome that are present in multiple Formica clades to detect genes that may have contributed to the emergence and maintenance of the genetic polymorphism. Finally, we infer three gene duplications on one haplotype, based on apparent heterozygosity within these genes in the genomes of haploid males. This study strengthens the evidence that this supergene originated early in the evolution of Formica and that just a few loci in this large region of suppressed recombination retain strongly differentiated alleles across contemporary Formica lineages.

RNase-based self-incompatibility in cacti

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

Primary restriction of S-RNase cytotoxicity by a stepwise ubiquitination and degradation pathway in Petunia hybrida

In self-incompatible Petunia species, the pistil S-RNase acts as cytotoxin to inhibit self-pollination but is polyubiquitinated by the pollen-specific nonself S-locus F-box (SLF) proteins and subsequently degraded by the ubiquitin-proteasome system (UPS), allowing cross-pollination. However, it remains unclear how S-RNase is restricted by the UPS. Using biochemical analyses, we first show that Petunia hybrida S₃ -RNase is largely ubiquitinated by K48-linked polyubiquitin chains at three regions, R I, R II and R III. R I is ubiquitinated in unpollinated, self-pollinated and cross-pollinated pistils, indicating its occurrence before PhS₃ -RNase uptake into pollen tubes, whereas R II and R III are exclusively ubiquitinated in cross-pollinated pistils. Transgenic analyses showed that removal of R II ubiquitination resulted in significantly reduced seed sets from cross-pollination and that of R I and R III to a lesser extent, indicating their increased cytotoxicity. Consistent with this, the mutated R II of PhS₃ -RNase resulted in a marked reduction of its degradation, whereas that of R I and R III resulted in less reduction. Taken together, we demonstrate that PhS₃ -RNase R II functions as a major ubiquitination region for its destruction and R I and R III as minor ones, revealing that its cytotoxicity is primarily restricted by a stepwise UPS mechanism for cross-pollination in P. hybrida.

Molecular characteristics of S-RNase alleles as the determinant of self-incompatibility in the style of Fragaria viridis

Strawberry (Fragaria spp.) is a member of the Rosoideae subfamily in the family Rosaceae. The self-incompatibility (SI) of some diploid species is a key agronomic trait that acts as a basic pollination barrier; however, the genetic mechanism underlying SI control in strawberry remains unclear. Two candidate S-RNases (S_a- and S_b-RNase) identified in the transcriptome of the styles of the self-incompatible Fragaria viridis 42 were confirmed to be SI determinants at the S locus following genotype identification and intraspecific hybridization using selfing progenies. Whole-genome collinearity and RNase T2 family analysis revealed that only an S locus exists in Fragaria; however, none of the compatible species contained S-RNase. Although the results of interspecific hybridization experiments showed that F. viridis (SI) styles could accept pollen from F. mandshurica (self-compatible), the reciprocal cross was incompatible. S_a and S_b-RNase contain large introns, and their noncoding sequences (promotors and introns) can be transcribed into long noncoding RNAs (lncRNAs). Overall, the genus Fragaria exhibits S-RNase-based gametophytic SI, and S-RNase loss occurs at the S locus of compatible germplasms. In addition, a type of SI-independent unilateral incompatibility exists between compatible and incompatible Fragaria species. Furthermore, the large introns and neighboring lncRNAs in S-RNase in Fragaria could offer clues about S-RNase expression strategies.

Spread of self-compatibility constrained by an intrapopulation crossing barrier

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

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Genome-wide identification and expression pattern analysis of the ribonuclease T2 family in Eucommia ulmoides

The 2′,3′-cycling ribonuclease (RNase) genes are catalysts of RNA cleavage and include the RNase T2 gene family. RNase T2 genes perform important roles in plants and have been conserved in the genome of eukaryotic organisms. In this study we identified 21 EURNS genes in Eucommia ulmoides Oliver (E. ulmoides) and analyzed their structure, chromosomal location, phylogenetic tree, gene duplication, stress-related cis-elements, and expression patterns in different tissues. The length of 21 predicted EURNS proteins ranged from 143 to 374 amino acids (aa), their molecular weight (MW) ranged from 16.21 to 42.38 kDa, and their isoelectric point (PI) value ranged from 5.08 to 9.09. Two classifications (class I and class III) were obtained from the conserved domains analysis and phylogenetic tree. EURNS proteins contained a total of 15 motifs. Motif 1, motif 2, motif 3, and motif 7 were distributed in multiple sequences and were similar to the conserved domain of RNase T2. EURNS genes with similar structure and the predicted EURNS proteins with conserved motif compositions are in the same group in the phylogenetic tree. The results of RT-PCR and transcription data showed that EURNS genes have tissue-specific expression and exhibited obvious trends in different developmental stages. Gene duplication analysis results indicated that segment duplication may be the dominant duplication mode in this gene family. This study provides a theoretical basis for research on the RNase T2 gene family and lays a foundation for the further study of EURNS genes.

Association of T2/S-RNase With Self-Incompatibility of Japanese Citrus Accessions Examined by Transcriptomic, Phylogenetic, and Genetic Approaches

Several citrus varieties show gametophytic self-incompatibility (GSI), which can contribute to seedless fruit production in several cultivars. This study investigated the genes regulating this trait through RNA-seq performed using styles collected from the flowers of Japanese citrus cultivars 'Hyuganatsu,' 'Tosabuntan,' 'Hassaku,' 'Banpeiyu,' and 'Sweet Spring'. We screened the transcripts of putative T2 RNases, i.e., the protein family including all S-RNases from S-RNase-based GSI plants, and constructed a phylogenetic tree using the screened T2 RNases and S-RNases retrieved from citrus genome databases and a public database. Three major clusters (class I-III) were formed, among which, the class III cluster contained family specific subclusters formed by S-RNase and a citrus-specific cluster monophyletic to the S-RNase clusters. From the citrus class III cluster, six transcripts were consistent with the S haplotypes previously determined in Japanese citrus accessions, sharing characteristics such as isoelectric point, extracellular localization, molecular weight, intron number and position, and tissue-specific expression with S-RNases. One T2 RNase gene in self-incompatible Hyuganatsu was significantly down-regulated in the styles of a self-compatible mutant of Hyuganatsu in RNA-seq and qPCR analyses. In addition, the inheritance pattern of some T2 RNase genes was consistent with the pattern of the S haplotype in the progeny population of Hyuganatsu and Tosabuntan. As all results supported citrus self-incompatibility being based on S-RNase, we believe that six T2 RNase genes were S-RNases. The homology comparison between the six T2 RNases and S-RNases recently reported in Chinese citrus revealed that three out of six T2 RNases were identical to S-RNases from Chinese citrus. Thus, the other three T2 RNases were finally concluded to be novel citrus S-RNases involved in self-incompatibility.

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

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

Phosphate-Starvation-Inducible S-Like RNase Genes in Rice Are Involved in Phosphate Source Recycling by RNA Decay

The fine-tuning of inorganic phosphate (Pi) for enhanced use efficiency has long been a challenging subject in agriculture, particularly in regard to rice as a major crop plant. Among ribonucleases (RNases), the RNase T2 family is broadly distributed across kingdoms, but little has been known on its substrate specificity compared to RNase A and RNase T1 families. Class I and class II of the RNase T2 family are defined as the S-like RNase (RNS) family and have showed the connection to Pi recycling in Arabidopsis. In this study, we first carried out a phylogenetic analysis of eight rice and five Arabidopsis RNS genes and identified mono-specific class I and dicot-specific class I RNS genes, suggesting the possibility of functional diversity between class I RNS family members in monocot and dicot species through evolution. We then compared the in silico expression patterns of all RNS genes in rice and Arabidopsis under normal and Pi-deficient conditions and further confirmed the expression patterns of rice RNS genes via qRT-PCR analysis. Subsequently, we found that most of the OsRNS genes were differentially regulated under Pi-deficient treatment. Association of Pi recycling by RNase activity in rice was confirmed by measuring total RNA concentration and ribonuclease activity of shoot and root samples under Pi-sufficient or Pi-deficient treatment during 21 days. The total RNA concentrations were decreased by < 60% in shoots and < 80% in roots under Pi starvation, respectively, while ribonuclease activity increased correspondingly. We further elucidate the signaling pathway of Pi starvation through upregulation of the OsRNS genes. The 2-kb promoter region of all OsRNS genes with inducible expression patterns under Pi deficiency contains a high frequency of P1BS cis-acting regulatory element (CRE) known as the OsPHR2 binding site, suggesting that the OsRNS family is likely to be controlled by OsPHR2. Finally, the dynamic transcriptional regulation of OsRNS genes by overexpression of OsPHR2, ospho2 mutant, and overexpression of OsPT1 lines involved in Pi signaling pathway suggests the molecular basis of OsRNS family in Pi recycling via RNA decay under Pi starvation.

Organellar and Secretory Ribonucleases: Major Players in Plant RNA Homeostasis

Organellar and secretory RNases, associated with different cellular compartments, are essential to maintain cellular homeostasis during development and in stress responses.

Tomato T2 ribonuclease LE is involved in the response to pathogens

T2 ribonucleases (RNases) are RNA-degrading enzymes that function in various cellular processes, mostly via RNA metabolism. T2 RNase-encoding genes have been identified in various organisms, from bacteria to mammals, and are most diverse in plants. The existence of T2 RNase genes in almost every organism suggests an important biological function that has been conserved through evolution. In plants, T2 RNases are suggested to be involved in phosphate scavenging and recycling, and are implicated in defence responses to pathogens. We investigated the function of the tomato T2 RNase LE, known to be induced by phosphate deficiency and wounding. The possible involvement of LE in pathogen responses was examined. Expression analysis showed LE induction during fungal infection and by stimuli known to be associated with pathogen inoculation, including oxalic acid and hydrogen peroxide. Analysis of LE-suppressed transgenic tomato lines revealed higher susceptibility to oxalic acid, a cell death-inducing factor, compared to the wild type. This elevated sensitivity of LE-suppressed lines was evidenced by visual signs of necrosis, and increased ion leakage and reactive oxygen species levels, indicating acceleration of cell death. Challenge of the LE-suppressed lines with the necrotrophic pathogen Botrytis cinerea resulted in accelerated development of disease symptoms compared to the wild type, associated with suppressed expression of pathogenesis-related marker genes. The results suggest a role for plant endogenous T2 RNases in antifungal activity.

Evolution of self-compatibility by a mutant Sm-RNase in citrus

Self-incompatibility (SI) is an important mechanism that prevents self-fertilization and inbreeding in flowering plants. The most widespread SI system utilizes S ribonucleases (S-RNases) and S-locus F-boxes (SLFs) as S determinants. In citrus, SI is ancestral, and Citrus maxima (pummelo) is self-incompatible, while Citrus reticulata (mandarin) and its hybrids are self-compatible (SC). Here, we identify nine highly polymorphic pistil-specific, developmentally expressed S-RNases from pummelo that segregate with S haplotypes in a gametophytic manner and cluster with authentic S-RNases. We provide evidence that these S-RNases function as the female S determinants in citrus. Moreover, we show that each S-RNase is linked to approximately nine SLFs. In an analysis of 117 citrus SLF and SFL-like (SLFL) genes, we reveal that they cluster into 12 types and that the S-RNases and intra-haplotypic SLF and SLFL genes co-evolved. Our data support the notion that citrus have a S locus comprising a S-RNase and several SLFs that fit the non-self-recognition model. We identify a predominant single nucleotide mutation, S_m-RNase, in SC citrus, which provides a 'natural' loss of function. We show that SI-SC transitions due to the S_m-RNase initially arose in mandarin, spreading to its hybrids and became fixed. Identification of an evolutionarily distant new genus utilizing the S-RNase-based SI system, >100 million years separated from the nearest S-RNase family, is a milestone for evolutionary comparative studies.

Nuclease and ribonuclease activities in response to salt stress: Identification of PvRNS3, a T2/S-like ribonuclease induced in common bean radicles by salt stress

The increase in soil salinization due to global climate change could cause large losses in crop productivity affecting, among other biological processes, to germination and seedling development. We have studied how salt stress affects nucleic acid degrading activities in radicles of common bean during seedling development. In radicles of common bean, a main nuclease of 37 kDa and two ribonucleases of 17 and 19 kDa were detected. Saline stress did not alter these three activities but induced a new ribonuclease of 16 kDa. All three ribonucleases are acidic enzymes that were inhibited by Zn. The 16 and 17 kDa ribonucleases are inhibited by guanilates. In the genome of common bean, we have identified 13 genes belonging to the T2 ribonuclease family and that are grouped in the 3 classes of T2 ribonucleases. The analysis of the expression of the 3 genes belonging to Class I (PvRNS1 to 3) and the unique gene from Class II (PvRNS4) in radicles showed that PvRNS3 is highly induced under salt stress.

Variation among S-locus haplotypes and among stylar RNases in almond

In many plant species, self-incompatibility systems limit self-pollination and mating among relatives. This helps maintain genetic diversity in natural populations but imposes constraints in agriculture and plant breeding. In almond [Prunus dulcis (Mill.) D.A. Webb], the specificity of self-incompatibility is mainly determined by stylar ribonuclease (S-RNase) and S-haplotype-specific F-box (SFB) proteins, both encoded within a complex locus, S. Prior to this research, a nearly complete sequence was available for one S-locus haplotype. Here, we report complete sequences for four haplotypes and partial sequences for 11 haplotypes. Haplotypes vary in sequences of genes (particularly S-RNase and SFB), distances between genes and numbers and positions of long terminal repeat transposons. Haplotype variation outside of the S-RNase and SFB genes may help maintain functionally important associations between S-RNase and SFB alleles. Fluorescence-based assays were developed to distinguish among some S-RNase alleles. With three-dimensional modelling of five S-RNase proteins, conserved active sites were identified and variation was observed in electrostatic potential and in the numbers, characteristics and positions of secondary structural elements, loop anchoring points and glycosylation sites. A hypervariable region on the protein surface and differences in the number, location and types of glycosylation sites may contribute to determining S-RNase specificity.

Interaction among ploidy, breeding system and lineage diversification

If particular traits consistently affect rates of speciation and extinction, broad macroevolutionary patterns can be interpreted as consequences of selection at high levels of the biological hierarchy. Identifying traits associated with diversification rates is difficult because of the wide variety of characters under consideration and the statistical challenges of testing for associations from comparative phylogenetic data. Ploidy (diploid vs polyploid states) and breeding system (self-incompatible vs self-compatible states) are both thought to be drivers of differential diversification in angiosperms. We fit 29 diversification models to extensive trait and phylogenetic data in Solanaceae and investigate how speciation and extinction rate differences are associated with ploidy, breeding system, and the interaction between these traits. We show that diversification patterns in Solanaceae are better explained by breeding system and an additional unobserved factor, rather than by ploidy. We also find that the most common evolutionary pathway to polyploidy in Solanaceae occurs via direct breakdown of self-incompatibility by whole genome duplication, rather than indirectly via breakdown followed by polyploidization. Comparing multiple stochastic diversification models that include complex trait interactions alongside hidden states enhances our understanding of the macroevolutionary patterns in plant phylogenies.

Factors to affect inbred beet plants while developing material for linear selection

Considering its capacities, the generative system of Beta vulgaris L. is regarded as highly productive. While inbreeding, the reproductive potential of cross-pollinated beet plants with gametophytic self-incompatibility (SI) changes significantly and is determined by a joint effect of multiple factors including the level of inbred depression. In the present study, original data have been obtained revealing relationships between inbred beet seed productivity, its self-incompatibility and microgametophyte parameters, which is crucial for developing and maintaining constant fertile beet lines. It has been discovered that inbred depression increases the number of sterile microgametes and anomalous pollen grains, reduces pollen fertility and the length of pollen tubes. As a result, the seed yield in inbred beet progeny, including SI ones, reduces significantly just after the third inbreeding. At the same time, highly productive inbred beet is characterized by a lower rate of pollen tube growth in vitro. In inbred plants, there is no close relationship between pollen viability and seed productivity, because the elimination of germinated male gametes and degeneration of seed embryos may go over the entire period of fertilization starting its progamic phase. The SI plants have more degenerating embryos than self-fertile ones, but seed vessel outgrowth in the seeds with abortive embryos makes them morphologically similar to fertile seeds. For that reason, when assessing inbred beet plants based on their self-incompatibility/self-fertility, one should consider the qualitative characteristics of the seeds. Using the method of recurrent selection based on such factors as seed productivity, pollen tube length and field germination rate increase the output of plant forms with a potentially high self-compatibility in their progeny. To support such genotypes in the progeny, one has to, starting from the third inbreeding, perform sib crossing to reduce the negative effect of inbred depression and self-incompatibility.

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

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

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

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

BDBM 1.0: A Desktop Application for Efficient Retrieval and Processing of High-Quality Sequence Data and Application to the Identification of the Putative Coffea S-Locus

Nowadays, bioinformatics is one of the most important areas in modern biology and the creation of high-quality scientific software supporting this recent research area is one of the core activities of many researchers. In this context, high-quality sequence datasets are needed to perform inferences on the evolution of species, genes, and gene families, or to get evidence for adaptive amino acid evolution, among others. Nevertheless, sequence data are very often spread over several databases, many useful genomes and transcriptomes are non-annotated, the available annotation is not for the desired coding sequence isoform, and/or is unlikely to be accurate. Moreover, although the FASTA text-based format is quite simple and usable by most software applications, there are a number of issues that may be critical depending on the software used to analyse such files. Therefore, researchers without training in informatics often use a fraction of all available data. The above issues can be addressed using already available software applications, but there is no easy-to-use single piece of software that allows performing all these tasks within the same graphical interface, such as the one here presented, named BDBM (Blast DataBase Manager). BDBM can be used to efficiently get gene sequences from annotated and non-annotated genomes and transcriptomes. Moreover, it can be used to look for alternatives to existing annotations and to easily create reliable custom databases. Such databases are essential to prepare high-quality datasets. The analyses that we have performed on the Coffea canephora genome using BDBM aimed at the identification of the S-locus region (that harbours the genes involved in gametophytic self-incompatibility) led to the conclusion that there are two likely regions, one on chromosome 2 (around region 6600000-6650000), and another on chromosome 5 (around 15830000-15930000). Such findings are discussed in the context of the Rubiaceae gametophytic self-incompatibility evolution.

Overcoming Self-Incompatibility in Diploid Potato Using CRISPR-Cas9

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

Immune Modulation by Human Secreted RNases at the Extracellular Space

The ribonuclease A superfamily is a vertebrate-specific family of proteins that encompasses eight functional members in humans. The proteins are secreted by diverse innate immune cells, from blood cells to epithelial cells and their levels in our body fluids correlate with infection and inflammation processes. Recent studies ascribe a prominent role to secretory RNases in the extracellular space. Extracellular RNases endowed with immuno-modulatory and antimicrobial properties can participate in a wide variety of host defense tasks, from performing cellular housekeeping to maintaining body fluid sterility. Their expression and secretion are induced in response to a variety of injury stimuli. The secreted proteins can target damaged cells and facilitate their removal from the focus of infection or inflammation. Following tissue damage, RNases can participate in clearing RNA from cellular debris or work as signaling molecules to regulate the host response and contribute to tissue remodeling and repair. We provide here an overall perspective on the current knowledge of human RNases’ biological properties and their role in health and disease. The review also includes a brief description of other vertebrate family members and unrelated extracellular RNases that share common mechanisms of action. A better knowledge of RNase mechanism of actions and an understanding of their physiological roles should facilitate the development of novel therapeutics.