Molecular mechanisms of self-incompatibility in Brassicaceae and Solanaceae.

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

Homology-Based Interactions between Small RNAs and Their Targets Control Dominance Hierarchy of Male Determinant Alleles of Self-Incompatibility in Arabidopsis lyrata.

Self-incompatibility (SI) is conserved among members of the Brassicaceae plant family. This trait is controlled epigenetically by the dominance hierarchy of the male determinant alleles. We previously demonstrated that a single small RNA (sRNA) gene is sufficient to control the linear dominance hierarchy in Brassica rapa and proposed a model in which a homology-based interaction between sRNAs and target sites controls the complicated dominance hierarchy of male SI determinants. In Arabidopsis halleri, male dominance hierarchy is reported to have arisen from multiple networks of sRNA target gains and losses. Despite these findings, it remains unknown whether the molecular mechanism underlying the dominance hierarchy is conserved among Brassicaceae. Here, we identified sRNAs and their target sites that can explain the linear dominance hierarchy of Arabidopsis lyrata, a species closely related to A. halleri. We tested the model that we established in Brassica to explain the linear dominance hierarchy in A. lyrata. Our results suggest that the dominance hierarchy of A. lyrata is also controlled by a homology-based interaction between sRNAs and their targets.

MYB transcription factor gene involved in sex determination in Asparagus officinalis.

Dioecy is a plant mating system in which individuals of a species are either male or female. Although many flowering plants evolved independently from hermaphroditism to dioecy, the molecular mechanism underlying this transition remains largely unknown. Sex determination in the dioecious plant Asparagus officinalis is controlled by X and Y chromosomes; the male and female karyotypes are XY and XX, respectively. Transcriptome analysis of A. officinalis buds showed that a MYB-like gene, Male Specific Expression 1 (MSE1), is specifically expressed in males. MSE1 exhibits tight linkage with the Y chromosome, specific expression in early anther development and loss of function on the X chromosome. Knockout of the MSE1 orthologue in Arabidopsis induces male sterility. Thus, MSE1 acts in sex determination in A. officinalis.

Ca2+-activated reactive oxygen species production by Arabidopsis RbohH and RbohJ is essential for proper pollen tube tip growth.

In flowering plants, pollen germinates on the stigma and pollen tubes grow through the style to fertilize the ovules. Enzymatic production of reactive oxygen species (ROS) has been suggested to be involved in pollen tube tip growth. Here, we characterized the function and regulation of the NADPH oxidases RbohH and RbohJ (Respiratory burst oxidase homolog H and J) in pollen tubes in Arabidopsis thaliana. In the rbohH and rbohJ single mutants, pollen tube tip growth was comparable to that of the wild type; however, tip growth was severely impaired in the double mutant. In vivo imaging showed that ROS accumulation in the pollen tube was impaired in the double mutant. Both RbohH and RbohJ, which contain Ca(2+) binding EF-hand motifs, possessed Ca(2+)-induced ROS-producing activity and localized at the plasma membrane of the pollen tube tip. Point mutations in the EF-hand motifs impaired Ca(2+)-induced ROS production and complementation of the double mutant phenotype. We also showed that a protein phosphatase inhibitor enhanced the Ca(2+)-induced ROS-producing activity of RbohH and RbohJ, suggesting their synergistic activation by protein phosphorylation and Ca(2+). Our results suggest that ROS production by RbohH and RbohJ is essential for proper pollen tube tip growth, and furthermore, that Ca(2+)-induced ROS positive feedback regulation is conserved in the polarized cell growth to shape the long tubular cell.

A pollen coat-inducible autoinhibited Ca2+-ATPase expressed in stigmatic papilla cells is required for compatible pollination in the Brassicaceae.

In the Brassicaceae, intraspecific non-self pollen (compatible pollen) can germinate and grow into stigmatic papilla cells, while self-pollen or interspecific pollen is rejected at this stage. However, the mechanisms underlying this selective acceptance of compatible pollen remain unclear. Here, using a cell-impermeant calcium indicator, we showed that the compatible pollen coat contains signaling molecules that stimulate Ca(2+) export from the papilla cells. Transcriptome analyses of stigmas suggested that autoinhibited Ca(2+)-ATPase13 (ACA13) was induced after both compatible pollination and compatible pollen coat treatment. A complementation test using a yeast Saccharomyces cerevisiae strain lacking major Ca(2+) transport systems suggested that ACA13 indeed functions as an autoinhibited Ca(2+) transporter. ACA13 transcription increased in papilla cells and in transmitting tracts after pollination. ACA13 protein localized to the plasma membrane and to vesicles near the Golgi body and accumulated at the pollen tube penetration site after pollination. The stigma of a T-DNA insertion line of ACA13 exhibited reduced Ca(2+) export, as well as defects in compatible pollen germination and seed production. These findings suggest that stigmatic ACA13 functions in the export of Ca(2+) to the compatible pollen tube, which promotes successful fertilization.

Efficient expression of SRK intracellular domain by a modeling-based protein engineering.

S-locus protein kinase (SRK) is a receptor kinase that plays a critical role in self-recognition in the Brassicaceae self-incompatibility (SI) response. SRK is activated by binding of its ligand S-locus protein 11 (SP11) and subsequently induced phosphorylation of the intracellular kinase domain. However, a detailed activation mechanism of SRK is still largely unknown because of the difficulty in stably expressing SRK recombinant proteins. Here, we performed modeling-based protein engineering of the SRK kinase domain for stable expression in Escherichia coli. The engineered SRK intracellular domain was expressed about 54-fold higher production than wild type SRK, without loss of the kinase activity, suggesting it could be useful for further biochemical and structural studies.

Cullin1-P is an Essential Component of Non-Self Recognition System in Self-Incompatibility in Petunia.

Self-incompatibility (SI) in flowering plants is a genetic reproductive barrier to distinguish self- and non-self pollen to promote outbreeding. In Solanaceae, self-pollen is rejected by the ribonucleases expressed in the styles (S-RNases), via its cytotoxic function. On the other side, the male-determinant is the S-locus F-box proteins (SLFs) expressed in pollen. Multiple SLFs collaboratively detoxify non-self S-RNases, therefore, non-self recognition is the mode of self-/non-self discrimination in Solanaceae. It is considered that SLFs function as a substrate-recognition module of the Skp1-Cullin1-F-box (SCF) complex that inactivates non-self S-RNases via their polyubiquitination, which leads to degradation by 26S proteasome. In fact, PhSSK1 (Petunia hybrida SLF-interacting Skp1-like1) was identified as a specific component of SCFSLF and was shown to be essential for detoxification of S-RNase in Petunia However, different molecules are proposed as the candidate Cullin1, another component of SCFSLF, and there is as yet no definite conclusion. Here, we identified five Cullin1s from the expressed sequence tags (ESTs) derived from the male reproductive organ in Petunia Among them, only PhCUL1-P was co-immunoprecipitated with S7-SLF2. In vitro protein-binding assay suggested that PhSSK1 specifically forms a complex with PhCUL1-P in an SLF-dependent manner. Knockdown of PhCUL1-P suppressed fertility of transgenic pollen in cross-compatible pollination in the functional S-RNase-dependent manner. These results suggested that SCFSLF selectively uses PhCUL1-P. Phylogeny of Cullin1s indicates that CUL1-P is recruited into the SI machinery during the evolution of Solanaceae, suggesting that the SI components have evolved differently among species in Solanaceae and Rosaceae, despite both families sharing the S-RNase-based SI.

Trans-acting small RNA determines dominance relationships in Brassica self-incompatibility.

A diploid organism has two copies of each gene, one inherited from each parent. The expression of two inherited alleles is sometimes biased by the effects known as dominant/recessive relationships, which determine the final phenotype of the organism. To explore the mechanisms underlying these relationships, we have examined the monoallelic expression of S-locus protein 11 genes (SP11), which encode the male determinants of self-incompatibility in Brassica. We previously reported that SP11 expression was monoallelic in some S heterozygotes, and that the promoter regions of recessive SP11 alleles were specifically methylated in the anther tapetum. Here we show that this methylation is controlled by trans-acting small non-coding RNA (sRNA). We identified inverted genomic sequences that were similar to the recessive SP11 promoters in the flanking regions of dominant SP11 alleles. These sequences were specifically expressed in the anther tapetum and processed into 24-nucleotide sRNA, named SP11 methylation inducer (Smi). Introduction of the Smi genomic region into the recessive S homozygotes triggered the methylation of the promoter of recessive SP11 alleles and repressed their transcription. This is an example showing sRNA encoded in the flanking region of a dominant allele acts in trans to induce transcriptional silencing of the recessive allele. Our finding may provide new insights into the widespread monoallelic gene expression systems.

Evolution of self-compatibility in Arabidopsis by a mutation in the male specificity gene.

Ever since Darwin's pioneering research, the evolution of self-fertilisation (selfing) has been regarded as one of the most prevalent evolutionary transitions in flowering plants. A major mechanism to prevent selfing is the self-incompatibility (SI) recognition system, which consists of male and female specificity genes at the S-locus and SI modifier genes. Under conditions that favour selfing, mutations disabling the male recognition component are predicted to enjoy a relative advantage over those disabling the female component, because male mutations would increase through both pollen and seeds whereas female mutations would increase only through seeds. Despite many studies on the genetic basis of loss of SI in the predominantly selfing plant Arabidopsis thaliana, it remains unknown whether selfing arose through mutations in the female specificity gene (S-receptor kinase, SRK), male specificity gene (S-locus cysteine-rich protein, SCR; also known as S-locus protein 11, SP11) or modifier genes, and whether any of them rose to high frequency across large geographic regions. Here we report that a disruptive 213-base-pair (bp) inversion in the SCR gene (or its derivative haplotypes with deletions encompassing the entire SCR-A and a large portion of SRK-A) is found in 95% of European accessions, which contrasts with the genome-wide pattern of polymorphism in European A. thaliana. Importantly, interspecific crossings using Arabidopsis halleri as a pollen donor reveal that some A. thaliana accessions, including Wei-1, retain the female SI reaction, suggesting that all female components including SRK are still functional. Moreover, when the 213-bp inversion in SCR was inverted and expressed in transgenic Wei-1 plants, the functional SCR restored the SI reaction. The inversion within SCR is the first mutation disrupting SI shown to be nearly fixed in geographically wide samples, and its prevalence is consistent with theoretical predictions regarding the evolutionary advantage of mutations in male components.

Ubiquitin-proteasome-mediated degradation of S-RNase in a solanaceous cross-compatibility reaction.

Many plants have a self-incompatibility (SI) system in which the rejection of self-pollen is determined by multiple haplotypes at a single locus, termed S. In the Solanaceae, each haplotype encodes a single ribonuclease (S-RNase) and multiple S-locus F-box proteins (SLFs), which function as the pistil and pollen SI determinants, respectively. S-RNase is cytotoxic to self-pollen, whereas SLFs are thought to collaboratively recognize non-self S-RNases in cross-pollen and detoxify them via the ubiquitination pathway. However, the actual mechanism of detoxification remains unknown. Here we isolate the components of a SCF(SLF) (SCF = SKP1-CUL1-F-box-RBX1) from Petunia pollen. The SCF(SLF) polyubiquitinates a subset of non-self S-RNases in vitro. The polyubiquitinated S-RNases are degraded in the pollen extract, which is attenuated by a proteasome inhibitor. Our findings suggest that multiple SCF(SLF) complexes in cross-pollen polyubiquitinate non-self S-RNases, resulting in their degradation by the proteasome.

Transcriptional characteristics and differences in Arabidopsis stigmatic papilla cells pre- and post-pollination.

Pollination is an important early step in sexual plant reproduction. In Arabidopsis thaliana, sequential pollination events, from pollen adhesion onto the stigma surface to pollen tube germination and elongation, occur on the stigmatic papilla cells. Following successful completion of these events, the pollen tube penetrates the stigma and finally fertilizes a female gametophyte. The pollination events are thought to be initiated and regulated by interactions between papilla cells and pollen. Here, we report the characterization of gene expression profiles of unpollinated (UP), compatible pollinated (CP) and incompatible pollinated (IP) papilla cells in A. thaliana. Based on cell type-specific transcriptome analysis from a combination of laser microdissection and RNA sequencing, 15,475, 17,360 and 16,918 genes were identified as expressed in UP, CP and IP papilla cells, respectively, and, of these, 14,392 genes were present in all three data sets. Differentially expressed gene (DEG) analyses identified 147 and 71 genes up-regulated in CP and IP papilla cells, respectively, and 115 and 46 genes down-regulated. Gene Ontology and metabolic pathway analyses revealed that papilla cells play an active role as the female reproductive component in pollination, particularly in information exchange, signal transduction, internal physiological changes and external morphological modification. This study provides fundamental information on the molecular mechanisms involved in pollination in papilla cells, furthering our understanding of the reproductive role of papilla cells.

Cytoplasmic Ca2+ changes dynamically during the interaction of the pollen tube with synergid cells.

The directional growth of the pollen tube from the stigma to the embryo sac in the ovules is regulated by pollen-pistil interactions based on intercellular communication. Although pollen tube growth is regulated by the cytoplasmic Ca(2+) concentration ([Ca(2+)](cyt)), it is not known whether [Ca(2+)](cyt) is involved in pollen tube guidance and reception. Using Arabidopsis expressing the GFP-based Ca(2+)-sensor yellow cameleon 3.60 (YC3.60) in pollen tubes and synergid cells, we monitored Ca(2+) dynamics in these cells during pollen tube guidance and reception under semi-in vivo fertilization conditions. In the pollen tube growing towards the micropyle, pollen tubes initiated turning within 150 µm of the micropylar opening; the [Ca(2+)](cyt) in these pollen tube tips was higher than in those not growing towards an ovule in assays with myb98 mutant ovules, in which pollen tube guidance is disrupted. These results suggest that attractants secreted from the ovules affect Ca(2+) dynamics in the pollen tube. [Ca(2+)](cyt) in synergid cells did not change when the pollen tube grew towards the micropyle or entered the ovule. Upon pollen tube arrival at the synergid cell, however, [Ca(2+)](cyt) oscillation began at the micropylar pole of the synergid, spreading towards the chalazal pole. Finally, [Ca(2+)](cyt) in the synergid cell reached a maximum at pollen tube rupture. These results suggest that signals from the pollen tube induce Ca(2+) oscillations in synergid cells, and that this Ca(2+) oscillation is involved in the interaction between the pollen tube and synergid cell.

Dominance relationships between self-incompatibility alleles controlled by DNA methylation.

In crucifers, the pollen S-determinant gene, SP11, is sporophytically expressed in the anther tapetum, and the pollen self-incompatibility phenotype is determined by the dominance relationships between the two S-haplotypes it carries. We report here that 5' promoter sequences of recessive SP11 alleles are specifically methylated in the tapetum before the initiation of SP11 transcription. These results suggest that tissue-specific monoallelic de novo DNA methylation is involved in determining the dominance interactions that determine the cruciferous self-incompatibility phenotype.

Physiological and genetic analysis of CO2-induced breakdown of self-incompatibility in Brassica rapa.

Self-incompatibility (SI) of the Brassicaceae family can be overcome by CO2 gas treatment. This method has been used for decades as an effective means to obtain a large amount of inbred seeds which can then be used for F1 hybrid seed production; however, the molecular mechanism by which CO2 alters the SI pathway has not been elucidated. In this study, to obtain new insights into the mechanism of CO2-induced SI breakdown, the focus was on two inbred lines of Brassica rapa (syn. campestris) with different CO2 sensitivity. Physiological examination using X-ray microanalysis suggested that SI breakdown in the CO2-sensitive line was accompanied by a significant accumulation of calcium at the pollen-stigma interface. Pre-treatment of pollen or pistil with CO2 gas before pollination showed no effect on the SI reaction, suggesting that some physiological process after pollination is necessary for SI to be overcome. Genetic analyses using F1 progeny of a CO2-sensitive × CO2-insensitive cross suggested that CO2 sensitivity is a semi-dominant trait in these lines. Analysis of F2 progeny suggested that CO2 sensitivity could be a quantitative trait, which is controlled by more than one gene. Quantitative trait locus (QTL) analyses identified two major loci, BrSIO1 and BrSIO2, which work additively in overcoming SI during CO2 treatment. No QTL was detected at the loci previously shown to affect SI stability, suggesting that CO2 sensitivity is determined by novel genes. The QTL data presented here should be useful for determining the responsible genes, and for the marker-assisted selection of desirable parental lines with stable but CO2-sensitive SI in F1 hybrid breeding.

Cell type-specific transcriptome of Brassicaceae stigmatic papilla cells from a combination of laser microdissection and RNA sequencing.

Pollination is an early and critical step in plant reproduction, leading to successful fertilization. It consists of many sequential processes, including adhesion of pollen grains onto the surface of stigmatic papilla cells, foot formation to strengthen pollen-stigma interaction, pollen hydration and germination, and pollen tube elongation and penetration. We have focused on an examination of the expressed genes in papilla cells, to increase understanding of the molecular systems of pollination. From three representative species of Brassicaceae (Arabidopsis thaliana, A. halleri and Brassica rapa), stigmatic papilla cells were isolated precisely by laser microdissection, and cell type-specific gene expression in papilla cells was determined by RNA sequencing. As a result, 17,240, 19,260 and 21,026 unigenes were defined in papilla cells of A. thaliana, A. halleri and B. rapa, respectively, and, among these, 12,311 genes were common to all three species. Among the17,240 genes predicted in A. thaliana, one-third were papilla specific while approximately half of the genes were detected in all tissues examined. Bioinformatics analysis revealed that genes related to a wide range of reproduction and development functions are expressed in papilla cells, particularly metabolism, transcription and membrane-mediated information exchange. These results reflect the conserved features of general cellular function and also the specific reproductive role of papilla cells, highlighting a complex cellular system regulated by a diverse range of molecules in these cells. This study provides fundamental biological knowledge to dissect the molecular mechanisms of pollination in papilla cells and will shed light on our understanding of plant reproduction mechanisms.

The transcription factor OsNAC4 is a key positive regulator of plant hypersensitive cell death.

The hypersensitive response (HR) is a common feature of plant immune responses and a type of programmed cell death. However, little is known about the induction mechanism of HR cell death. We report that overexpression of OsNAC4, which encodes a plant-specific transcription factor, leads to HR cell death accompanied by the loss of plasma membrane integrity, nuclear DNA fragmentation and typical morphological changes. In OsNAC4 knock-down lines, HR cell death is markedly decreased in response to avirulent bacterial strains. After induction by an avirulent pathogen recognition signal, OsNAC4 is translocated into the nucleus in a phosphorylation-dependent manner. A microarray analysis showed that the expression of 139 genes including OsHSP90 and IREN, encoding a Ca(2+)-dependent nuclease, were different between the OsNAC4 knock-down line and control line during HR cell death. During the induction of HR cell death, OsHSP90 is involved in the loss of plasma membrane integrity, whereas IREN causes nuclear DNA fragmentation. Overall, our results indicate that two important events occurring during HR cell death are regulated by independent pathways.

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

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

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

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

SHI family transcription factors regulate an interspecific barrier.

Pre-zygotic interspecies incompatibility in angiosperms is an important mechanism to prevent unfavourable hybrids between species. Here we report our identification of STIGMATIC PRIVACY 2 (SPRI2), a transcription factor that has a zinc-finger domain and regulates interspecies barriers in Arabidopsis thaliana, via genome-wide association study. Knockout analysis of SPRI2/SRS7 and its paralogue SPRI2-like/SRS5 demonstrated their necessity in rejecting male pollen from other species within female pistils. Additionally, they govern mRNA transcription of xylan O-acetyltransferases (TBL45 and TBL40) related to cell wall modification, alongside SPRI1, a pivotal transmembrane protein for interspecific pollen rejection. SPRI2/SRS7 is localized as condensed structures in the nucleus formed via liquid-liquid phase separation (LLPS), and a prion-like sequence in its amino-terminal region was found to be responsible for the formation of the condensates. The LLPS-regulated SPRI2/SRS7 discovered in this study may contribute to the establishment of interspecific reproductive barriers through the transcriptional regulation of cell wall modification genes and SPRI1.

Dominance in self-compatibility between subgenomes of allopolyploid Arabidopsis kamchatica shown by transgenic restoration of self-incompatibility.

The evolutionary transition to self-compatibility facilitates polyploid speciation. In Arabidopsis relatives, the self-incompatibility system is characterized by epigenetic dominance modifiers, among which small RNAs suppress the expression of a recessive SCR/SP11 haplogroup. Although the contribution of dominance to polyploid self-compatibility is speculated, little functional evidence has been reported. Here we employ transgenic techniques to the allotetraploid plant A. kamchatica. We find that when the dominant SCR-B is repaired by removing a transposable element insertion, self-incompatibility is restored. This suggests that SCR was responsible for the evolution of self-compatibility. By contrast, the reconstruction of recessive SCR-D cannot restore self-incompatibility. These data indicate that the insertion in SCR-B conferred dominant self-compatibility to A. kamchatica. Dominant self-compatibility supports the prediction that dominant mutations increasing selfing rate can pass through Haldane's sieve against recessive mutations. The dominance regulation between subgenomes inherited from progenitors contrasts with previous studies on novel epigenetic mutations at polyploidization termed genome shock.