Finding a Compatible Partner: Self-Incompatibility in European Pear (Pyrus communis); Molecular Control, Genetic Determination, and Impact on Fertilization and Fruit Set

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

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

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

Genome-wide identification and characterization of SRLK gene family reveal their roles in self-incompatibility of Erigeron breviscapus

Self-incompatibility (SI) is a reproductive protection mechanism that plants acquired during evolution to prevent self-recession. As the female determinant of SI specificity, SRK has been shown to be the only recognized gene on the stigma and plays important roles in SI response. Asteraceae is the largest family of dicotyledonous plants, many of which exhibit self-incompatibility. However, systematic studies on SRK gene family in Asteraceae are still limited due to lack of high-quality genomic data. In this study, we performed the first systematic genome-wide identification of S-locus receptor like kinases (SRLKs) in the self-incompatible Asteraceae species, Erigeron breviscapus, which is also a widely used perennial medicinal plant endemic to China.52 SRLK genes were identified in the E. breviscapus genome. Structural analysis revealed that the EbSRLK proteins in E. breviscapus are conserved. SRLK proteins from E. breviscapus and other SI plants are clustered into 7 clades, and the majority of the EbSRLK proteins are distributed in Clade I. Chromosomal and duplication analyses indicate that 65% of the EbSRLK genes belong to tandem repeats and could be divided into six tandem gene clusters. Gene expression patterns obtained in E. breviscapus multiple-tissue RNA-Seq data revealed differential temporal and spatial features of EbSRLK genes. Among these, two EbSRLK genes having high expression levels in tongue flowers were cloned. Subcellular localization assay demonstrated that both of their fused proteins are localized on the plasma membrane. All these results indicated that EbSRLK genes possibly involved in SI response in E. breviscapus. This comprehensive genome-wide study of the SRLK gene family in E. breviscapus provides valuable information for understanding the mechanism of SSI in Asteraceae.

MIK2 is a candidate gene of the S-locus for sporophytic self-incompatibility in chicory (Cichorium intybus, Asteraceae)

The Cichorium genus offers a unique opportunity to study the sporophytic self-incompatibility (SSI) system, being composed of species characterized by highly efficient self-incompatibility (e.g., C. intybus) and complete self-compatibility (e.g., C. endivia). To this end, the chicory genome was used to map seven previously identified SSI locus-associated markers. The region containing the S-locus was therefore restricted to an ~4 M bp window on chromosome 5. Among the genes predicted in this region, MDIS1 INTERACTING RECEPTOR LIKE KINASE 2 (ciMIK2) was particularly promising as a candidate for SSI. Its ortholog in Arabidopsis (atMIK2) is involved in pollen-stigma recognition reactions, and its protein structure is similar to that of S-receptor kinase (SRK), a key component of the SSI system in the Brassica genus. The amplification and sequencing of MIK2 in chicory and endive accessions revealed two contrasting scenarios. In C. endivia, MIK2 was fully conserved even when comparing different botanical varieties (i.e., smooth and curly endive). In C. intybus, 387 polymorphic positions and 3 INDELs were identified when comparing accessions of different biotypes all belonging to the same botanical variety (i.e., radicchio). The polymorphism distribution throughout the gene was uneven, with hypervariable domains preferentially localized in the LRR-rich extracellular region, putatively identified as the receptor domain. The gene was hypothesized to be under positive selection, as the nonsynonymous mutations were more than double the synonymous ones (dN/dS = 2.17). An analogous situation was observed when analyzing the first 500 bp of the MIK2 promoter: no SNPs were observed among the endive samples, whereas 44 SNPs and 6 INDELs were detected among the chicory samples. Further analyses are needed to confirm the role of MIK2 in SSI and to demonstrate whether the 23 species-specific nonsynonymous SNPs in the CDS and/or the species-specific 10 bp-INDEL found in a CCAAT box region of the promoter are responsible for the contrasting sexual behaviors of chicory and endive.

BAG6-A from Fragaria viridis pollen modulates gametophyte development in diploid strawberry

Male and female gametophyte development processes are essential steps in the life cycles of all land plants. Here, we characterized a gene, FviBAG6-A, screened from the Fragaria viridis (2 n = 2x=14) pollen cDNA library and physically interacted with S-RNase. Ubiquitinated of S_a-RNase might be determined by the interaction of FviBAG6-A in the ubiquitin-proteasome system during fertilization. We found that overexpression of FviBAG6-A in Arabidopsis caused shorter silique length, and decreased silique number. Moreover, overexpression of FviBAG6-A in Fragaria vesca (2 n = 2x=14) led to a greatly reduced seed number, with nearly 80% of the seeds aborted. Analyses of paraffin sections and reactive oxygen species (ROS) content revealed that the majority of severe pollen defects were likely due to the early degradation of the tapetum and middle layer as a result of ROS accumulation and abnormal development of the uninucleate megaspore mother. Moreover, the FviBAG6-A interact with the E3 ligase SIZ1 and contribute to the SUMOylation of FviBAG6-A , which may be induced by the high level of ROS content, further promoting gametophyte abortion in strawberry transgenic lines. This study characterized the FviBAG6-A and reveals its novel function in gametophyte development.

Genome-wide identification of members of the Skp1 family in almond (Prunus dulcis), cloning and expression characterization of PsdSSK1

Skp1 (S-phase kinase-associated protein 1) is the core gene of SCF ubiquitin ligase, which mediates protein degradation, thereby regulating biological processes such as cell cycle progression, transcriptional regulation, and signal transduction. A variety of plant Skp1 gene family studies have been reported. However, the almond Skp1 gene family has not yet been studied. In this study, we identified 18 members of the Prunus dulcis PdSkp1 family that were unevenly distributed across six chromosomes of the almond genome. Phylogenetic tree analysis revealed that the PdSkp1 members can be divided into three groups: I, II, and III. PdSkp1 members in each subfamily have relatively conserved motif types and exon/intron numbers. There were three pairs of fragment duplication genes and one pair of tandem repeat genes, and their functions were highly evolutionarily conserved. Transcriptome data showed that PdSkp1 is expressed in almond flower tissues, and that its expression shows significant change during cross-pollination. Fluorescence quantitative results showed that eight PdSkp1 genes had different expression levels in five tissues of almond, i.e., branches, leaves, flower buds, flesh, and cores. In addition, we cloned a PsdSSK1 gene based on PdSkp1. The cloned PsdSSK1 showed the same protein sequence as PdSkp1-12. Results of qPCR and western blot analysis showed high expression of PsdSSK1 in almond pollen. In conclusion, we report the first clone of the key gene SSK1 that controls self-incompatibility in almonds. Our research lays a foundation for future functional research on PdSkp1 members, especially for exploring the mechanism of almond self-incompatibility.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01278-9.

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

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

Programmed Cell Death May Be Involved in the Seedless Phenotype Formation of Oil Palm

Oil palm (Elaeis guineensis Jacq.) is a well-known vegetable oil-yielding crop. Seedlessness is one of the most prominent traits in oil palm due to its low processing costs and high oil content. Nevertheless, an extensive study on molecular mechanisms regulating seedless phenotype formation in oil palm is very limited so far. In this study, stigma, style, and ovary from seedless and seeded (Tenera and Pisifera) oil palm trees were used to investigate the possible mechanism. Results showed that non-pollination resulted in no fruits, and self- and cross-pollinations resulted in seedless fruits, while boron treatment had no effect on seedless phenotype formation, implying that seedless trees have incomplete self and outcrossing incompatibility. Furthermore, the transcriptome data analysis highlighted eight programmed cell death (PCD) genes and three groups of PCD-related genes: 4-coumarate-CoA ligase (4CL), S-RNase, and MADS-box. The majority of these genes were significantly up-regulated in the stigma and style of Seedless palm trees compared to Tenera and Pisifera. In addition, the co-expression network analysis confirmed the significant correlation among these genes. Moreover, two simple sequence repeats (SSR) markers (S41 and S44) were developed to identify the seedless phenotype. The up-regulation of 4CL and MADS-box TFs activated the expression of PCD genes; on the other hand, S-RNase resulted in pollen tube RNA degradation and triggered PCD. While the link between PCD and seedless phenotype formation in oil palm has not been extensively studied to date, these findings suggest a role of PCD in pollen tube lethality, leading to double fertilization failure and the seedless phenotype.

The Multipartite Mitochondrial Genome of Marama (Tylosema esculentum)

Tylosema esculentum (marama bean), a wild legume from tropical Africa, has long been considered as a potential crop for local farmers due to its rich nutritional value. Genomics research of marama is indispensable for the domestication and varietal improvement of the bean. The chloroplast genome of marama has been sequenced and assembled previously using a hybrid approach based on both Illumina and PacBio data. In this study, a similar method was used to assemble the mitochondrial genome of marama. The mitochondrial genome of the experimental individual has been confirmed to have two large circles OK638188 and OK638189, which do not recombine according to the data. However, they may be able to restructure into five smaller circles through recombination on the 4 pairs of long repeats (>1 kb). The total length of marama mitogenome is 399,572 bp. A 9,798 bp DNA fragment has been found that is homologous to the chloroplast genome of marama, accounting for 2.5% of the mitogenome. In the Fabaceae family, the mitogenome of Millettia pinnata is highly similar to marama, including for both the genes present and the total size. Some genes including cox2, rpl10, rps1, and sdh4 have been lost during the evolution of angiosperms and are absent in the mitogenomes of some legumes. However, these remain intact and functional in marama. Another set of genes, rpl2, rps2, rps7, rps11, rps13, and rps19 are either absent, or present as pseudogenes, in the mitogenome of marama.

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.

Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (F_ROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that F_ROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.

Phenotypic assessment of pear varieties in mutual pollination

The self-fertility and cross-fertility of pear varieties have been studied in the conditions of the of the Krasnodar region. It was established that most of the not varieties are self-fertile. The varieties are distinguished: Ladozhskaya, Malyshka and Talgar beauty, with partial self-fertility, in which a low percentage of fruit set and its fluctuation over the years are noted. Based on the pollination performed and the analysis carried out, good and acceptable pollinators for promising pear varieties were selected. It has been established that in order to form full-fledged seeds in the fruit of the pear, a full cross-pollination with good pollinators is necessary. The possibilities of individual varieties (a high percentage of fruit tying) with mutual pollination are presented. The analysis of the negative influence of technogenic factors during flowering on the inter-pollination of varieties has been carried out. Cross-pollinated varieties have a narrow selective ability to pollinate varieties, which varies over the years and depends on the growing conditions. Proved the need, when laying gardens correctly choose varieties-pollinators in order to increase the yield of varieties. According to the research results, a self-fertile variety Iyulskaya early was identified. The best pollinators were selected for promising pear varieties for the conditions of the Kuban.

RNA-Seq analysis of compatible and incompatible styles of Pyrus species at the beginning of pollination

At the early stage of pollination, the difference in gene expression between compatibility and incompatibility is highly significant about the pollen-specific expression of the LRR gene, resistance, and defensin genes. In Rosaceae, incompatible pollen can penetrate into the style during the gametophytic self-incompatibility response. It is therefore considered a stylar event rather than a stigmatic event. In this study, we explored the differences in gene expression between compatibility and incompatibility in the early stage of pollination. The self-compatible pear variety "Jinzhuili" is a naturally occurring bud mutant from "Yali", a leading Chinese native cultivar exhibiting typical gametophytic self-incompatibility. We collected the styles of 'Yali' and 'Jinzhuili' at 0.5 and 2 h after self-pollination and then performed high-throughput sequencing. According to the KEGG analysis of the differentially expressed genes, several metabolic pathways, such as "Plant hormone signal transduction", "Plant-pathogen interaction", are the main pathways was the most represented pathway. Quantitative PCR was used to validate these differential genes. The expression levels of genes related to pollen growth and disease inhibition, such as LRR (Leucine-rich repeat extensin), resistance, defensin, and auxin, differed significantly between compatible and incompatible pollination. Interestingly, at 0.5 h, most of these genes were upregulated in the compatible pollination system compared with the incompatible pollination system. Calcium transport, which requires ATPase, also demonstrated upregulated expression. In summary, the self-incompatibility reaction was initiated when the pollen land on the stigma.

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

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

Genetic Diversity and Structure of the Portuguese Pear (Pyrus communis L.) Germplasm

A rich heritage of traditional pear varieties is kept in national Portuguese collections. Out of these varieties, “Rocha” dominates national pear production. Although a noticeable phenotypic variation among clones of this variety has been reported, little is known about its genetic variability, as to date molecular studies have been performed on a single “Rocha” clone. Eleven Simple Sequence Repeats (SSR) markers were used to assess the genetic diversity of 130 local cultivars, 80 of them being “Rocha” clones. The results allowed the differentiation of 75 genotypes of which 29 are “Rocha”. Three synonyms groups and four homonymous groups of other local varieties were confirmed. A Bayesian model-based clustering approach identified two distinct clusters. Using flow cytometry, six cultivars were found to be triploids. These results show high genetic variability among “Rocha” clones. In conclusion, there is a need for different “Rocha” clones to be preserved to enable the correct selection of the multiplication material.