Brassinosteroid signals cooperate with katanin-mediated microtubule severing to control stamen filament elongation.

Proper stamen filament elongation is essential for pollination and plant reproduction. Plant hormones are extensively involved in every stage of stamen development; however, the cellular mechanisms by which phytohormone signals couple with microtubule dynamics to control filament elongation remain unclear. Here, we screened a series of Arabidopsis thaliana mutants showing different microtubule defects and revealed that only those unable to sever microtubules, lue1 and ktn80.1234, displayed differential floral organ elongation with less elongated stamen filaments. Prompted by short stamen filaments and severe decrease in KTN1 and KTN80s expression in qui-2 lacking five BZR1-family transcription factors (BFTFs), we investigated the crosstalk between microtubule severing and brassinosteroid (BR) signaling. The BFTFs transcriptionally activate katanin-encoding genes, and the microtubule-severing frequency was severely reduced in qui-2. Taken together, our findings reveal how BRs can regulate cytoskeletal dynamics to coordinate the proper development of reproductive organs.

A Quantitative Computational Framework for Allopolyploid Single-Cell Data Integration and Core Gene Ranking in Development.

Polyploidization drives regulatory and phenotypic innovation. How the merger of different genomes contributes to polyploid development is a fundamental issue in evolutionary developmental biology and breeding research. Clarifying this issue is challenging because of genome complexity and the difficulty in tracking stochastic subgenome divergence during development. Recent single-cell sequencing techniques enabled probing subgenome-divergent regulation in the context of cellular differentiation. However, analyzing single-cell data suffers from high error rates due to high dimensionality, noise, and sparsity, and the errors stack up in polyploid analysis due to the increased dimensionality of comparisons between subgenomes of each cell, hindering deeper mechanistic understandings. In this study, we develop a quantitative computational framework, called "pseudo-genome divergence quantification" (pgDQ), for quantifying and tracking subgenome divergence directly at the cellular level. Further comparing with cellular differentiation trajectories derived from single-cell RNA sequencing data allows for an examination of the relationship between subgenome divergence and the progression of development. pgDQ produces robust results and is insensitive to data dropout and noise, avoiding high error rates due to multiple comparisons of genes, cells, and subgenomes. A statistical diagnostic approach is proposed to identify genes that are central to subgenome divergence during development, which facilitates the integration of different data modalities, enabling the identification of factors and pathways that mediate subgenome-divergent activity during development. Case studies have demonstrated that applying pgDQ to single-cell and bulk tissue transcriptomic data promotes a systematic and deeper understanding of how dynamic subgenome divergence contributes to developmental trajectories in polyploid evolution.

McAN: a novel computational algorithm and platform for constructing and visualizing haplotype networks.

Haplotype networks are graphs used to represent evolutionary relationships between a set of taxa and are characterized by intuitiveness in analyzing genealogical relationships of closely related genomes. We here propose a novel algorithm termed McAN that considers mutation spectrum history (mutations in ancestry haplotype should be contained in descendant haplotype), node size (corresponding to sample count for a given node) and sampling time when constructing haplotype network. We show that McAN is two orders of magnitude faster than state-of-the-art algorithms without losing accuracy, making it suitable for analysis of a large number of sequences. Based on our algorithm, we developed an online web server and offline tool for haplotype network construction, community lineage determination, and interactive network visualization. We demonstrate that McAN is highly suitable for analyzing and visualizing massive genomic data and is helpful to enhance the understanding of genome evolution. Availability: Source code is written in C/C++ and available at https://github.com/Theory-Lun/McAN and https://ngdc.cncb.ac.cn/biocode/tools/BT007301 under the MIT license. Web server is available at https://ngdc.cncb.ac.cn/bit/hapnet/. SARS-CoV-2 dataset are available at https://ngdc.cncb.ac.cn/ncov/. Contact: songshh@big.ac.cn (Song S), zhaowm@big.ac.cn (Zhao W), baoym@big.ac.cn (Bao Y), zhangzhang@big.ac.cn (Zhang Z), ybxue@big.ac.cn (Xue Y).

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.

The Snapdragon Genomes Reveal the Evolutionary Dynamics of the S-Locus Supergene.

The genus Antirrhinum has been used as a model to study self-incompatibility extensively. The multi-allelic S-locus, carrying a pistil S-RNase and dozens of S-locus F-box (SLF) genes, underlies the genetic control of self-incompatibility (SI) in Antirrhinum hispanicum. However, there have been limited studies on the genomic organization of the S-locus supergene due to a lack of high-quality genomic data. Here, we present the chromosome-level reference and haplotype-resolved genome assemblies of a self-incompatible A. hispanicum line, AhS7S8. For the first time, 2 complete A. hispanicum S-haplotypes spanning ∼1.2 Mb and containing a total of 32 SLFs were reconstructed, whereas most of the SLFs derived from retroelement-mediated proximal or tandem duplication ∼122 Mya. Back then, the S-RNase gene and incipient SLFs came into linkage to form the pro-type of type-1 S-locus in the common ancestor of eudicots. Furthermore, we detected a pleiotropic cis-transcription factor (TF) associated with regulating the expression of SLFs, and two miRNAs may control the expression of this TF. Interspecific S-locus and intraspecific S-haplotype comparisons revealed the dynamic nature and polymorphism of the S-locus supergene mediated by continuous gene duplication, segmental translocation or loss, and TE-mediated transposition events. Our data provide an excellent resource for future research on the evolutionary studies of the S-RNase-based self-incompatibility system.

Evolutionary rewiring of the wheat transcriptional regulatory network by lineage-specific transposable elements.

More than 80% of the wheat genome consists of transposable elements (TEs), which act as major drivers of wheat genome evolution. However, their contributions to the regulatory evolution of wheat adaptations remain largely unclear. Here, we created genome-binding maps for 53 transcription factors (TFs) underlying environmental responses by leveraging DAP-seq in Triticum urartu, together with epigenomic profiles. Most TF binding sites (TFBSs) located distally from genes are embedded in TEs, whose functional relevance is supported by purifying selection and active epigenomic features. About 24% of the non-TE TFBSs share significantly high sequence similarity with TE-embedded TFBSs. These non-TE TFBSs have almost no homologous sequences in non-Triticeae species and are potentially derived from Triticeae-specific TEs. The expansion of TE-derived TFBS linked to wheat-specific gene responses, suggesting TEs are an important driving force for regulatory innovations. Altogether, TEs have been significantly and continuously shaping regulatory networks related to wheat genome evolution and adaptation.

N-terminal acetylation orchestrates glycolate-mediated ROS homeostasis to promote rice thermoresponsive growth.

Climate warming poses a significant threat to global crop production and food security. However, our understanding of the molecular mechanisms governing thermoresponsive development in crops remains limited. Here we report that the auxiliary subunit of N-terminal acetyltransferase A (NatA) in rice OsNAA15 is a prerequisite for rice thermoresponsive growth. OsNAA15 produces two isoforms OsNAA15.1 and OsNAA15.2, via temperature-dependent alternative splicing. Among the two, OsNAA15.1 is more likely to form a stable and functional NatA complex with the potential catalytic subunit OsNAA10, leading to a thermoresponsive N-terminal acetylome. Intriguingly, while OsNAA15.1 promotes plant growth under elevated temperatures, OsNAA15.2 exhibits an inhibitory effect. We identified two glycolate oxidases (GLO1/5) as major substrates from the thermoresponsive acetylome. These enzymes are involved in hydrogen peroxide (H2O2) biosynthesis via glycolate oxidation. N-terminally acetylated GLO1/5 undergo their degradation through the ubiquitin-proteasome system. This leads to reduced reactive oxygen species (ROS) production, thereby promoting plant growth, particularly under high ambient temperatures. Conclusively, our findings highlight the pivotal role of N-terminal acetylation in orchestrating the glycolate-mediated ROS homeostasis to facilitate thermoresponsive growth in rice.

An atlas of wheat epigenetic regulatory elements reveals subgenome divergence in the regulation of development and stress responses.

Wheat (Triticum aestivum) has a large allohexaploid genome. Subgenome-divergent regulation contributed to genome plasticity and the domestication of polyploid wheat. However, the specificity encoded in the wheat genome determining subgenome-divergent spatio-temporal regulation has been largely unexplored. The considerable size and complexity of the genome are major obstacles to dissecting the regulatory specificity. Here, we compared the epigenomes and transcriptomes from a large set of samples under diverse developmental and environmental conditions. Thousands of distal epigenetic regulatory elements (distal-epiREs) were specifically linked to their target promoters with coordinated epigenomic changes. We revealed that subgenome-divergent activity of homologous regulatory elements is affected by specific epigenetic signatures. Subgenome-divergent epiRE regulation of tissue specificity is associated with dynamic modulation of H3K27me3 mediated by Polycomb complex and demethylases. Furthermore, quantitative epigenomic approaches detected key stress responsive cis- and trans-acting factors validated by DNA Affinity Purification and sequencing, and demonstrated the coordinated interplay between epiRE sequence contexts, epigenetic factors, and transcription factors in regulating subgenome divergent transcriptional responses to external changes. Together, this study provides a wealth of resources for elucidating the epiRE regulomics and subgenome-divergent regulation in hexaploid wheat, and gives new clues for interpreting genetic and epigenetic interplay in regulating the benefits of polyploid wheat.

Molecular and evolutionary mechanisms of self-incompatibility in angiosperms.

As an intraspecific outcrossing mechanism, self-incompatibility (SI) widely adopted by hermaphroditic plants is usually controlled by a polymorphic multi-allelic S locus. Typically, six molecular types of SI have been found, including type-I controlled by the pistil S S-RNase and pollen S SLFs commonly spread in Plantaginaceae, Solanaceae, Rosaceae and Rutaceae, type-II by SRK and SCR in Brassicaceae, type-III by PrsS and PrpS in Papaveraceae, type-IV by CYP-GLO2-KFB-CCM-PUM in Primulaceae, type-V by TsSPH1-TsYUC6-TsBAHD in Turneraceae and type-VI by HPS10-S and DUF247I-S in Poaceae, with type-I characterized as a non-self recognition system but types-II, -III and -VI self ones. Furthermore, remarkable progresses have been made in their origin and evolutionary mechanisms recently. Among them, type-I SI possessed a single origin in the most recent common ancestor of eudicots and types II-V dynamically evolved following its losses, while type-VI SI exclusively existed in monocot Poaceae may be regained after the loss of the ancient type-I. Here, we mainly review the molecular and evolutionary mechanisms of angiosperm SI systems, thus providing a helpful reference for their theoretical research and breeding application.

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

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

Villin controls the formation and enlargement of punctate actin foci in pollen tubes.

Self-incompatibility (SI) in the poppy Papaver rhoeas triggers dramatic alterations in actin within pollen tubes. However, how these actin alterations are mechanistically achieved remains largely unexplored. Here, we used treatment with the Ca2+ ionophore A23187 to mimic the SI-induced elevation in cytosolic Ca2+ and trigger formation of the distinctive F-actin foci. Live-cell imaging revealed that this remodeling involves F-actin fragmentation and depolymerization, accompanied by the rapid formation of punctate actin foci and subsequent increase in their size. We established that actin foci are generated and enlarged from crosslinking of fragmented actin filament structures. Moreover, we show that villins associate with actin structures and are involved in this actin reorganization process. Notably, we demonstrate that Arabidopsis VILLIN5 promotes actin depolymerization and formation of actin foci by fragmenting actin filaments, and controlling the enlargement of actin foci via bundling of actin filaments. Our study thus uncovers important novel insights about the molecular players and mechanisms involved in forming the distinctive actin foci in pollen tubes.

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

Selection and gene flow shape genomic islands that control floral guides.

Genomes of closely-related species or populations often display localized regions of enhanced relative sequence divergence, termed genomic islands. It has been proposed that these islands arise through selective sweeps and/or barriers to gene flow. Here, we genetically dissect a genomic island that controls flower color pattern differences between two subspecies of Antirrhinum majus, A.m.striatum and A.m.pseudomajus, and relate it to clinal variation across a natural hybrid zone. We show that selective sweeps likely raised relative divergence at two tightly-linked MYB-like transcription factors, leading to distinct flower patterns in the two subspecies. The two patterns provide alternate floral guides and create a strong barrier to gene flow where populations come into contact. This barrier affects the selected flower color genes and tightly-linked loci, but does not extend outside of this domain, allowing gene flow to lower relative divergence for the rest of the chromosome. Thus, both selective sweeps and barriers to gene flow play a role in shaping genomic islands: sweeps cause elevation in relative divergence, while heterogeneous gene flow flattens the surrounding "sea," making the island of divergence stand out. By showing how selective sweeps establish alternative adaptive phenotypes that lead to barriers to gene flow, our study sheds light on possible mechanisms leading to reproductive isolation and speciation.

RETINOBLASTOMA-RELATED Genes Specifically Control Inner Floral Organ Morphogenesis and Pollen Development in Rice.

RETINOBLASTOMA-RELATED (RBR) is an essential gene in plants, but its molecular function outside of its role in cell cycle entry remains poorly understood. We characterized the functions of OsRBR1 and OsRBR2 in plant growth and development in rice using both forward- and reverse-genetics methods. The two genes were coexpressed and performed redundant roles in vegetative organs but exhibited separate functions in flowers. OsRBR1 was highly expressed in the floral meristem and regulated the expression of floral homeotic genes to ensure floral organ formation. Mutation of OsRBR1 caused loss of floral meristem identity, resulting in the replacement of lodicules, stamens, and the pistil with either a panicle-like structure or whorls of lemma-like organs. OsRBR2 was preferentially expressed in stamens and promoted pollen formation. Mutation of OsRBR2 led to deformed anthers without pollen. Similar to the protein interaction between AtRBR and AtMSI1 that is essential for floral development in Arabidopsis, OsMSI1 was identified as an interaction partner of OsRBR1 and OsRBR2. OsMSI1 was ubiquitously expressed and appears to be essential for development in rice (Oryza sativa), as the mutation of OsMSI1 was lethal. These results suggest that OsRBR1 and OsRBR2 function with OsMSI1 in reproductive development in rice. This work characterizes further functions of RBRs and improves current understanding of specific regulatory pathways of floral specification and pollen formation in rice.

Ectopic expression of nucleolar DEAD-Box RNA helicase OsTOGR1 confers improved heat stress tolerance in transgenic Chinese cabbage.

KEY MESSAGE: The DEAD-Box RNA helicase OsTOGR1 positively regulates heat stress tolerance in Chinese cabbage. Non-heading Chinese cabbage (Brassica rapa L. ssp. chinensis) is primarily cultivated vegetable crop in Asian countries. Heat stress is one of the major threats for its growth and yield. Numerous regulatory genes in various crops have shown to contribute thermotolerance. Among them, Thermotolerant growth required 1 (TOGR1) is an important DEAD-box RNA helicase. To examine whether its role is conserved in other crops, we constructed pCAMBIA1300-pHSP:OsTOGR1 expression vector driven by the rice small heat shock protein promoter (pHSP17.9) and successfully produced transgenic non-heading Chinese cabbage plants expressing OsTOGR1 gene via Agrobacterium-mediated vacuum infiltration transformation. In total, we generated three independent transgenic cabbage lines expressing TOGR1 gene. Expression and integration of TOGR1 was confirmed by PCR, RT-PCR and qPCR in T1 and T2 generations. The relative leaf electrical conductivity of transgenic seedlings was reduced subjected to high temperature (38 °C) compared to heat shock treatment (46 °C). In addition, hypocotyl length of transgenic seedlings increased compared to wild-type plants under high temperature and heat shock treatment. Furthermore, the transgenic plants exhibited higher chlorophyll content than wild-type plants under high temperature and heat shock treatment. The transgenic seeds displayed better germination under heat shock treatment. Tested heat stress-responsive genes were also up-regulated in the transgenic plants subjected to high temperature or heat shock treatment. To the best of our knowledge, this is the first report on describing the role of DAED-Box RNA helicases in improving heat stress tolerance of transgenic plants.

The 2019 novel coronavirus resource.

An ongoing outbreak of a novel coronavirus infection in Wuhan, China since December 2019 has led to 31,516 infected persons and 638 deaths across 25 countries (till 16:00 on February 7, 2020). The virus causing this pneumonia was then named as the 2019 novel coronavirus (2019-nCoV) by the World Health Organization. To promote the data sharing and make all relevant information of 2019-nCoV publicly available, we construct the 2019 Novel Coronavirus Resource (2019nCoVR, https://bigd.big.ac.cn/ncov). 2019nCoVR features comprehensive integration of genomic and proteomic sequences as well as their metadata information from the Global Initiative on Sharing All Influenza Data, National Center for Biotechnology Information, China National GeneBank, National Microbiology Data Center and China National Center for Bioinformation (CNCB)/National Genomics Data Center (NGDC). It also incorporates a wide range of relevant information including scientific literatures, news, and popular articles for science dissemination, and provides visualization functionalities for genome variation analysis results based on all collected 2019-nCoV strains. Moreover, by linking seamlessly with related databases in CNCB/NGDC, 2019nCoVR offers virus data submission and sharing services for raw sequence reads and assembled sequences. In this report, we provide comprehensive descriptions on data deposition, management, release and utility in 2019nCoVR, laying important foundations in aid of studies on virus classification and origin, genome variation and evolution, fast detection, drug development and pneumonia precision prevention and therapy.

Evolutionary history of the human multigene families reveals widespread gene duplications throughout the history of animals.

BACKGROUND: The hypothesis that vertebrates have experienced two ancient, whole genome duplications (WGDs) is of central interest to evolutionary biology and has been implicated in evolution of developmental complexity. Three-way and Four-way paralogy regions in human and other vertebrate genomes are considered as vital evidence to support this hypothesis. Alternatively, it has been proposed that such paralogy regions are created by small-scale duplications that occurred at different intervals over the evolution of life. RESULTS: To address this debate, the present study investigates the evolutionary history of multigene families with at least three-fold representation on human chromosomes 1, 2, 8 and 20. Phylogenetic analysis and the tree topology comparisons classified the members of 36 multigene families into four distinct co-duplicated groups. Gene families falling within the same co-duplicated group might have duplicated together, whereas genes belong to different co-duplicated groups might have distinct evolutionary origins. CONCLUSION: Taken together with previous investigations, the current study yielded no proof in favor of WGDs hypothesis. Rather, it appears that the vertebrate genome evolved as a result of small-scale duplication events, that cover the entire span of the animals' history.

Nucleolar DEAD-Box RNA Helicase TOGR1 Regulates Thermotolerant Growth as a Pre-rRNA Chaperone in Rice.

Plants have evolved a considerable number of intrinsic tolerance strategies to acclimate to ambient temperature increase. However, their molecular mechanisms remain largely obscure. Here we report a DEAD-box RNA helicase, TOGR1 (Thermotolerant Growth Required1), prerequisite for rice growth themotolerance. Regulated by both temperature and the circadian clock, its expression is tightly coupled to daily temperature fluctuations and its helicase activities directly promoted by temperature increase. Located in the nucleolus and associated with the small subunit (SSU) pre-rRNA processome, TOGR1 maintains a normal rRNA homeostasis at high temperature. Natural variation in its transcript level is positively correlated with plant height and its overexpression significantly improves rice growth under hot conditions. Our findings reveal a novel molecular mechanism of RNA helicase as a key chaperone for rRNA homeostasis required for rice thermotolerant growth and provide a potential strategy to breed heat-tolerant crops by modulating the expression of TOGR1 and its orthologs.

A High Temperature-Dependent Mitochondrial Lipase EXTRA GLUME1 Promotes Floral Phenotypic Robustness against Temperature Fluctuation in Rice (Oryza sativa L.).

The sessile plants have evolved diverse intrinsic mechanisms to control their proper development under variable environments. In contrast to plastic vegetative development, reproductive traits like floral identity often show phenotypic robustness against environmental variations. However, it remains obscure about the molecular basis of this phenotypic robustness. In this study, we found that eg1 (extra glume1) mutants of rice (Oryza savita L.) showed floral phenotypic variations in different growth locations resulting in a breakdown of floral identity robustness. Physiological and biochemical analyses showed that EG1 encodes a predominantly mitochondria-localized functional lipase and functions in a high temperature-dependent manner. Furthermore, we found that numerous environmentally responsive genes including many floral identity genes are transcriptionally repressed in eg1 mutants and OsMADS1, OsMADS6 and OsG1 genetically act downstream of EG1 to maintain floral robustness. Collectively, our results demonstrate that EG1 promotes floral robustness against temperature fluctuation by safeguarding the expression of floral identify genes through a high temperature-dependent mitochondrial lipid pathway and uncovers a novel mechanistic insight into floral developmental control.

Electrostatic potentials of the S-locus F-box proteins contribute to the pollen S specificity in self-incompatibility in Petunia hybrida.

Self-incompatibility (SI) is a self/non-self discrimination system found widely in angiosperms and, in many species, is controlled by a single polymorphic S-locus. In the Solanaceae, Rosaceae and Plantaginaceae, the S-locus encodes a single S-RNase and a cluster of S-locus F-box (SLF) proteins to control the pistil and pollen expression of SI, respectively. Previous studies have shown that their cytosolic interactions determine their recognition specificity, but the physical force between their interactions remains unclear. In this study, we show that the electrostatic potentials of SLF contribute to the pollen S specificity through a physical mechanism of 'like charges repel and unlike charges attract' between SLFs and S-RNases in Petunia hybrida. Strikingly, the alteration of a single C-terminal amino acid of SLF reversed its surface electrostatic potentials and subsequently the pollen S specificity. Collectively, our results reveal that the electrostatic potentials act as a major physical force between cytosolic SLFs and S-RNases, providing a mechanistic insight into the self/non-self discrimination between cytosolic proteins in angiosperms.