Haplotype-resolved chromosomal-level assembly of wasabi (Eutrema japonicum) genome.

In Japan, wasabi (Eutrema japonicum) is an important traditional condiment, and is recognized as an endemic species. In the present study, we generated a chromosome-level and haplotype-resolved reference genome for E. japonicum using PacBio CLR (continuous long reads), Illumina, and Hi-C sequencing data. The genome consists of 28 chromosomes that contain 1,512.1 Mb of sequence data, with a scaffold N50 length of 55.67 Mb. We also reported the subgenome and haplotype assignment of the 28 chromosomes by read-mapping and phylogenic analysis. Three validation methods (Benchmarking Universal Single-Copy Orthologs, Merqury, and Inspector) indicated that our obtained genome sequences were a high-quality and high-completeness genome assembly. Comparison of genome assemblies from previously published genomes showed that our obtained genome was of higher quality. Therefore, our genome will serve as a valuable genetic resource for both chemical ecology and evolution research of the genera Eutrema and Brassicaceae, as well as for wasabi breeding.

Practical synthesis of 1,3-benzoazaphosphole analogues.

Despite the promising pharmacological activity and material properties of five-membered heterocyclic compounds containing phosphorus and nitrogen, synthetic examples of them have been rather limited due to the instability of phosphorus toward air and water. In this study, 1,3-benzoazaphosphol analogues were selected as target molecules, and various synthetic routes were examined to establish a fundamental technology for the introduction of phosphorus groups into aromatic rings and formation of five-membered rings containing phosphorus and nitrogen by cyclization. As a result, we found that 2-aminophenyl(phenyl)phosphine is an extremely promising synthetic intermediate with high stability and easy handling. Furthermore, 2-methyl-3-phenyl-2,3-dihydro-1H-benzo[d][1,3]azaphosphole and 3-phenyl-2,3-dihydro-1H-benzo[d][1,3]azaphosphole-2-thione as synthetically useful 1,3-benzoazaphosphol analogues were successfully synthesized by using 2-aminophenyl(phenyl)phosphine as a key intermediate.

P-Glycoprotein-Mediated Pharmacokinetic Interactions Increase Pimozide hERG Channel Inhibition.

Pimozide, an antipsychotic drug, is a potent inhibitor of the hERG channel. A case of death due to cardiac arrest has been reported in a boy who received pimozide together with sertraline and aripiprazole. In this study, we focused on drug-drug interactions and investigated the relationships between transporter-mediated intracellular accumulation and the hERG inhibitory effect of pimozide. The accumulation of pimozide in cardiomyocyte-derived AC16 cells was significantly increased by sertraline and aripiprazole, which are thought to have a P-glycoprotein (P-gp) inhibitory effect, and under P-gp siRNA conditions. These results suggest P-gp inhibition increases pimozide accumulation in AC16 cells. We introduced the hERG plasmid into AC16 cells and investigated the concentration-dependent hERG inhibitory effect of pimozide from within AC16 cells. Addition of 10 nM or more pimozide significantly inhibited the hERG current with concentration dependence. These results indicate P-gp-mediated pharmacokinetic interaction increases pimozide accumulation in AC16 cells, and the subsequent elevated pimozide levels within the cells may result in an increased risk of hERG channel inhibition. Our present study calls attention to the risks associated with the combined use of cardiotoxic P-gp substrate(s) and P-gp inhibitory medicines.

Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis.

Land plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is a continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning spermatogenesis remain to be elucidated. We successfully identified candidate genes involved in spermatogenesis, deeply divergent BLD10s, by computational analyses combining multiple methods and omics data. We then examined the functions of BLD10s in the liverwort Marchantia polymorpha and the moss Physcomitrium patens. MpBLD10 and PpBLD10 are required for normal basal body (BB) and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, and thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. We identified orthologs of MpBLD10 and PpBLD10 in diverse Streptophyta and found that MpBLD10 and PpBLD10 are orthologous to BLD10/CEP135 family proteins, which function in BB assembly. However, BLD10s evolved especially quickly in land plants and MpBLD10 might have acquired additional functions in spermatozoid formation through rapid molecular evolution.

DOMINANT AWN INHIBITOR Encodes the ALOG Protein Originating from Gene Duplication and Inhibits AWN Elongation by Suppressing Cell Proliferation and Elongation in Sorghum.

The awn, a needle-like structure extending from the tip of the lemma in grass species, plays a role in environmental adaptation and fitness. In some crops, awns appear to have been eliminated during domestication. Although numerous genes involved in awn development have been identified, several dominant genes that eliminate awns are also known to exist. For example, in sorghum (Sorghum bicolor), the dominant awn-inhibiting gene has been known since 1921; however, its molecular features remain uncharacterized. In this study, we conducted quantitative trait locus analysis and a genome-wide association study of awn-related traits in sorghum and identified DOMINANT AWN INHIBITOR (DAI), which encodes the ALOG family protein on chromosome 3. DAI appeared to be present in most awnless sorghum cultivars, likely because of its effectiveness. Detailed analysis of the ALOG protein family in cereals revealed that DAI originated from a duplication of its twin paralog (DAIori) on chromosome 10. Observations of immature awns in near-isogenic lines revealed that DAI inhibits awn elongation by suppressing both cell proliferation and elongation. We also found that only DAI gained a novel function to inhibit awn elongation through an awn-specific expression pattern distinct from that of DAIori. Interestingly, heterologous expression of DAI with its own promoter in rice inhibited awn elongation in the awned cultivar Kasalath. We found that DAI originated from gene duplication, providing an interesting example of gain-of-function that occurs only in sorghum but shares its functionality with rice and sorghum.

Expression and localization of tight junction-related proteins in adult rat pituitary stem/progenitor cell niches.

Pituitary endocrine cells are supplied by Sox2-expressing stem/progenitor cells in the anterior lobe of the adult pituitary gland. These SOX2-positive cells are maintained in two types of microenvironments (niches): the marginal cell layer (MCL)-niche and the parenchymal-niche. Recently, we isolated dense SOX2-positive cell clusters from the parenchymal-niche by taking advantage of their resistance to protease treatment as parenchymal stem/progenitor cell (PS)-clusters. In the present study, by analyzing these isolated PS-clusters, we attempted to identify novel structural characteristics of pituitary stem/progenitor cell niches. Quantitative real-time PCR showed that tight junction-related genes were distinctly expressed in the isolated PS-clusters. Immunocytostaining showed that the tight junction molecules, ZO-1 and occludin, were localized in the apical membrane facing the pseudo-follicle-like structure of the isolated PS-clusters regardless of the expression of S100ß, which distinguishes the sub-population of SOX2-positive cells. Furthermore, immunohistochemistry of the pituitary glands of adult rats clearly demonstrated that ZO-1 and occludin were densely present in the parenchymal-niche encircling the pseudo-follicle, while they were observed in the apical membrane in the MCL-niche facing the residual lumen. Collectively, these tight junction-related proteins might be involved in the architecture and maintenance of the plasticity of pituitary stem/progenitor cell niches.

Assessment of Greenhouse Tomato Anthesis Rate Through Metabolomics Using LASSO Regularized Linear Regression Model.

While the high year-round production of tomatoes has been facilitated by solar greenhouse cultivation, these yields readily fluctuate in response to changing environmental conditions. Mathematic modeling has been applied to forecast phenotypes of tomatoes using environmental measurements (e.g., temperature) as indirect parameters. In this study, metabolome data, as direct parameters reflecting plant internal status, were used to construct a predictive model of the anthesis rate of greenhouse tomatoes. Metabolome data were obtained from tomato leaves and used as variables for linear regression with the least absolute shrinkage and selection operator (LASSO) for prediction. The constructed model accurately predicted the anthesis rate, with an R2 value of 0.85. Twenty-nine of the 161 metabolites were selected as candidate markers. The selected metabolites were further validated for their association with anthesis rates using the different metabolome datasets. To assess the importance of the selected metabolites in cultivation, the relationships between the metabolites and cultivation conditions were analyzed via correspondence analysis. Trigonelline, whose content did not exhibit a diurnal rhythm, displayed major contributions to the cultivation, and is thus a potential metabolic marker for predicting the anthesis rate. This study demonstrates that machine learning can be applied to metabolome data to identify metabolites indicative of agricultural traits.

Temporal changes in transcripts of miniature inverted-repeat transposable elements during rice endosperm development.

The repression of transcription from transposable elements (TEs) by DNA methylation is necessary to maintain genome integrity and prevent harmful mutations. However, under certain circumstances, TEs may escape from the host defense system and reactivate their transcription. In Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), DNA demethylases target the sequences derived from TEs in the central cell, the progenitor cell for the endosperm in the female gametophyte. Genome-wide DNA demethylation is also observed in the endosperm after fertilization. In the present study, we used a custom microarray to survey the transcripts generated from TEs during rice endosperm development and at selected time points in the embryo as a control. The expression patterns of TE transcripts are dynamically up- and downregulated during endosperm development, especially those of miniature inverted-repeat TEs (MITEs). Some TE transcripts were directionally controlled, whereas the other DNA transposons and retrotransposons were not. We also discovered the NUCLEAR FACTOR Y binding motif, CCAAT, in the region near the 5' terminal inverted repeat of Youren, one of the transcribed MITEs in the endosperm. Our results uncover dynamic changes in TE activity during endosperm development in rice.

De novo assembly and annotation of the mangrove cricket genome.

OBJECTIVES: The mangrove cricket, Apteronemobius asahinai, shows endogenous activity rhythms that synchronize with the tidal cycle (i.e., a free-running rhythm with a period of ~ 12.4 h [the circatidal rhythm]). Little is known about the molecular mechanisms underlying the circatidal rhythm. We present the draft genome of the mangrove cricket to facilitate future molecular studies of the molecular mechanisms behind this rhythm. DATA DESCRIPTION: The draft genome contains 151,060 scaffolds with a total length of 1.68 Gb (N50: 27 kb) and 92% BUSCO completeness. We obtained 28,831 predicted genes, of which 19,896 (69%) were successfully annotated using at least one of two databases (UniProtKB/SwissProt database and Pfam database).

Laser Vaporization of Intracoronary Thrombus and Identifying Plaque Morphology in ST-Segment Elevation Myocardial Infarction as Assessed by Optical Coherence Tomography.

OBJECTIVES: We evaluated the thrombus-vaporizing effect of excimer laser coronary angioplasty (ELCA) in patients with ST-segment elevation myocardial infarction (STEMI) by optical coherence tomography (OCT). BACKGROUND: Larger intracoronary thrombus elevates the risk of interventional treatment and mortality in patients with STEMI. METHODS: A total of 92 patients with STEMI who presented within 24 hours from the onset and underwent ELCA following manual aspiration thrombectomy (MT) were analyzed. RESULTS: The mean baseline thrombolysis in myocardial infarction flow grade was 0.4 ± 0.6, which subsequently improved to 2.3 ± 0.7 after MT (p < 0.0001) and 2.7 ± 0.5 after ELCA (p=0.0001). The median residual thrombus volume after MT was 65.7 mm3, which significantly reduced to 47.5 mm3 after ELCA (p < 0.0001). Plaque rupture was identified by OCT in only 22 cases (23.9%) after MT, but was distinguishable in 36 additional cases after ELCA (total: 58 cases; 63.0%). Ruptured lesions contained a higher proportion of red thrombus than nonruptured lesions (75.9% vs. 43.3%, p=0.001). Significantly larger thrombus burden after MT (69.6 mm3 vs. 56.3 mm3, p < 0.05) and greater thrombus reduction by ELCA (21.2 mm3 vs. 11.8 mm3, p < 0.01) were observed in ruptured lesions than nonruptured lesions. CONCLUSIONS: ELCA effectively vaporized intracoronary thrombus in patients with STEMI even after MT. Lesions with plaque rupture contained larger thrombus burden that was frequently characterized by red thrombus and more effectively reduced by ELCA.

Drug resistance via radixin-mediated increase of P-glycoprotein membrane expression during SNAI1-induced epithelial-mesenchymal transition in HepG2 cells.

OBJECTIVES: Epithelial-mesenchymal transition (EMT) plays a role in cancer metastasis as well as in drug resistance through various mechanisms, including increased drug efflux mediated by P-glycoprotein (P-gp). In this study, we investigated the activation mechanism of P-gp, including its regulatory factors, during EMT in hepatoblastoma-derived HepG2 cells. METHODS: HepG2 cells were transfected with SNAI1 using human adenovirus serotype 5 vector. We quantified mRNA and protein expression levels using qRT-PCR and western blot analysis, respectively. P-gp activity was evaluated by uptake assay, and cell viability was assessed by an MTT assay. KEY FINDINGS: P-gp protein expression on plasma membrane was higher in SNAI1-transfected cells than in Mock cells, although there was no difference in P-gp protein level in whole cells. Among the scaffold proteins such as ezrin, radixin and moesin (ERM), only radixin was increased in SNAI1-transfected cells. Uptake of both Rho123 and paclitaxel was decreased in SNAI1-transfected cells, and this decrease was blocked by verapamil, a P-gp inhibitor. The reduced susceptibility of SNAI1-transfected cells to paclitaxel was reversed by elacridar, another P-gp inhibitor. CONCLUSIONS: Increased expression of radixin during SNAI1-induced EMT leads to increased P-gp membrane expression in HepG2 cells, enhancing P-gp function and thereby increasing drug resistance.

Combination of genetic analysis and ancient literature survey reveals the divergence of traditional Brassica rapa varieties from Kyoto, Japan.

Since ancient times, humans have bred several plants that we rely on today. However, little is known about the divergence of most of these plants. In the present study, we investigated the divergence of Mibuna (Brassica rapa L. subsp. nipposinica L. H. Bailey), a traditional leafy vegetable in Kyoto (Japan), by combining genetic analysis and a survey of ancient literature. Mibuna is considered to have been bred 200 years ago from Mizuna, another traditional leafy vegetable in Kyoto. Mibuna has simple spatulate leaves, whereas Mizuna has characteristic serrated leaves. The quantitative trait loci (QTL) and gene expression analyses suggested that the downregulation of BrTCP15 expression contributed to the change in the leaf shape from serrated to simple spatulate. Interestingly, the SNP analysis indicated that the genomic region containing the BrTCP15 locus was transferred to Mibuna by introgression. Furthermore, we conducted a survey of ancient literature to reveal the divergence of Mibuna and found that hybridization between Mizuna and a simple-leaved turnip might have occurred in the past. Indeed, the genomic analysis of multiple turnip cultivars showed that one of the cultivars, Murasakihime, has almost the same sequence in the BrTCP15 region as Mibuna. These results suggest that the hybridization between Mizuna and turnip has resulted in the establishment of Mibuna.

Rapid Increase of Gastrointestinal P-Glycoprotein Functional Activity in Response to Etoposide Stimulation.

We previously reported that exposure of human colon adenocarcinoma (Caco-2) cells to the bitter substance phenylthiocarbamide (PTC) rapidly enhanced the transport function of P-glycoprotein (P-gp). In this study, we investigated the short-term effect of etoposide, another bitter-tasting P-gp substrate, on P-gp transport function in the same cell line. We found that etoposide exposure significantly increased both the P-gp protein level in the plasma membrane fraction and the efflux rate of rhodamine123 (Rho123) in Caco-2 cells within 10 min. The efflux ratio (ratio of the apparent permeability coefficient in the basal-to-apical direction to that in the apical-to-basal direction) of Rho123 in etoposide-treated cells was also significantly increased compared with the control. These results indicated that etoposide rapidly enhances P-gp function in Caco-2 cells. In contrast, P-gp expression in whole cells at both the mRNA and protein level was unchanged by etoposide exposure, compared with the levels in non-treated cells. Furthermore, etoposide increased the level of phosphorylated ezrin, radixin and moesin (P-ERM) proteins in the plasma membrane fraction of Caco-2 cells within 10 min. P-gp functional changes were blocked by YM022, an inhibitor of cholecystokinin (CCK) receptor. These results suggest that etoposide induces release of CCK, causing activation of the CCK receptor followed by phosphorylation of ERM proteins, which recruit intracellular P-gp for trafficking to the gastrointestinal membrane, thereby increasing the functional activity of P-gp.

Effect of Paternal Genome Excess on the Developmental and Gene Expression Profiles of Polyspermic Zygotes in Rice.

Polyploid zygotes with a paternal gamete/genome excess exhibit arrested development, whereas polyploid zygotes with a maternal excess develop normally. These observations indicate that paternal and maternal genomes synergistically influence zygote development via distinct functions. In this study, to clarify how paternal genome excess affects zygotic development, the developmental and gene expression profiles of polyspermic rice zygotes were analyzed. The results indicated that polyspermic zygotes were mostly arrested at the one-cell stage after karyogamy had completed. Through comparison of transcriptomes between polyspermic zygotes and diploid zygotes, 36 and 43 genes with up-regulated and down-regulated expression levels, respectively, were identified in the polyspermic zygotes relative to the corresponding expression in the diploid zygotes. Notably, OsASGR-BBML1, which encodes an AP2 transcription factor possibly involved in initiating rice zygote development, was expressed at a much lower level in the polyspermic zygotes than in the diploid zygotes.

Functional Alterations of Multidrug Resistance-Associated Proteins 2 and 5, and Breast Cancer Resistance Protein upon Snail-Induced Epithelial-Mesenchymal Transition in HCC827 Cells.

Our previous report indicated that Snail-induced epithelial-mesenchymal transition (EMT) enhanced P-glycoprotein (P-gp) function and drug resistance to P-gp substrate anticancer drug in a human non-small cell lung cancer (NSCLC) cell line, HCC827. Our objective is to evaluate the changes in the mRNA and protein expression levels and the functions of multidrug resistance-associated protein (MRP) 2, MRP5 and breast cancer resistance protein (BCRP). Snail-expressing HCC827 cells showed increased mRNA levels of Snail and a mesenchymal marker vimentin, and decreased mRNA levels of an epithelial marker E-cadherin after transduction, indicating that Snail had induced EMT consistent with our previous reports. The mRNA level of MRP2 was significantly decreased, while that of MRP5 remained unchanged, in Snail-expressing cells. The expression levels of MRP2 and MRP5 proteins in whole-cell homogenate were unchanged in Snail-expressing cells, but MRP5 protein showed significantly increased membrane localization. Snail-transduction increased the efflux transport of 5-(and-6)-carboxy-2',7'-dichlorofluorescein (CDCF), a substrate of MRP2, 3 and 5. This increase was blocked by MK571, which inhibits MRP1, 2, and 5. Toxicity of cisplatin, a substrate of MRP2 and 5, was significantly decreased in Snail-expressing cells. BCRP mRNA and protein levels were both decreased in Snail-expressing cells, which showed an increase in the intracellular accumulation of 7-ethyl-10-hydroxycamptothecin (SN-38), a BCRP substrate, resulting in reduced viability. These results suggested that MRP5 function appears to be increased via an increase in membrane localization, whereas the BCRP function is decreased via a decrease in the expression level in HCC827 cells with Snail-induced EMT.

Physiological Roles of ERM Proteins and Transcriptional Regulators in Supporting Membrane Expression of Efflux Transporters as Factors of Drug Resistance in Cancer.

One factor contributing to the malignancy of cancer cells is the acquisition of drug resistance during chemotherapy via increased expression of efflux transporters, such as P-glycoprotein (P-gp), multidrug resistance-associated proteins (MRPs), and breast cancer resistance protein (BCRP). These transporters operate at the cell membrane, and are anchored in place by the scaffold proteins ezrin (Ezr), radixin (Rdx), and moesin (Msn) (ERM proteins), which regulate their functional activity. The identity of the regulatory scaffold protein(s) differs depending upon the transporter, and also upon the tissue in which it is expressed, even for the same transporter. Another factor contributing to malignancy is metastatic ability. Epithelial-mesenchymal transition (EMT) is the first step in the conversion of primary epithelial cells into mesenchymal cells that can be transported to other organs via the blood. The SNAI family of transcriptional regulators triggers EMT, and SNAI expression is used is an indicator of malignancy. Furthermore, EMT has been suggested to be involved in drug resistance, since drug excretion from cancer cells is promoted during EMT. We showed recently that ERM proteins are induced by a member of the SNAI family, Snail. Here, we first review recent progress in research on the relationship between efflux transporters and scaffold proteins, including the question of tissue specificity. In the second part, we review the relationship between ERM scaffold proteins and the transcriptional regulatory factors that induce their expression.

Genome editing in plants using CRISPR type I-D nuclease.

Genome editing in plants has advanced greatly by applying the clustered regularly interspaced short palindromic repeats (CRISPRs)-Cas system, especially CRISPR-Cas9. However, CRISPR type I-the most abundant CRISPR system in bacteria-has not been exploited for plant genome modification. In type I CRISPR-Cas systems, e.g., type I-E, Cas3 nucleases degrade the target DNA in mammals. Here, we present a type I-D (TiD) CRISPR-Cas genome editing system in plants. TiD lacks the Cas3 nuclease domain; instead, Cas10d is the functional nuclease in vivo. TiD was active in targeted mutagenesis of tomato genomic DNA. The mutations generated by TiD differed from those of CRISPR/Cas9; both bi-directional long-range deletions and short indels mutations were detected in tomato cells. Furthermore, TiD can be used to efficiently generate bi-allelic mutant plants in the first generation. These findings indicate that TiD is a unique CRISPR system that can be used for genome engineering in plants.

Gastrointestinal absorption of pimozide is enhanced by inhibition of P-glycoprotein.

Drug-drug interaction was suggested to have played a role in the recent death due to cardiac arrest of a patient taking pimozide, sertraline and aripiprazole antipsychotic/antidepressant combination therapy. Here, we investigated the possible involvement of P-glycoprotein (P-gp)-mediated interaction among these drugs, using in vitro methods. ATPase assay confirmed that pimozide is a P-gp substrate, and might act as a P-gp inhibitor at higher concentrations. The maximum transport rate (Jmax) and half-saturation concentration (Kt) for the carrier-mediated transport estimated by means of pimozide efflux assay using P-gp-overexpressing LLC-GA5-CoL150 cells were 84.9 ± 8.9 pmol/min/mg protein, and 10.6 ± 4.7 µM, respectively. These results indicate that pimozide is a good P-gp substrate, and it appears to have the potential to cause drug-drug interactions in the digestive tract at clinically relevant gastrointestinal concentrations. Moreover, sertraline or aripiprazole significantly decreased the efflux ratio of pimozide in LLC-GA5-CoL150 cells. Transport studies using Caco-2 cell monolayers were consistent with the results in LLC-GA5-CoL150 cells, and indicate that P-gp-mediated drug-drug interaction may occur in the gastrointestinal tract. Thus, P-gp inhibition by sertraline and/or aripiprazole may increase the gastrointestinal permeability of co-administered pimozide, resulting in an increased blood concentration of pimozide, which is known to be associated with an increased risk of QT prolongation, a life-threatening side effect.

Double-Locking Mechanism of Self-Compatibility in Arabidopsis thaliana: The Synergistic Effect of Transcriptional Depression and Disruption of Coding Region in the Male Specificity Gene.

Self-compatibility in Arabidopsis thaliana represents the relatively recent disruption of ancestral obligate cross pollination, recognized as one of the prevalent evolutionary pathways in flowering plants, as noted by Darwin. Our previous study found that inversion of the male specificity gene (SP11/SCR) disrupted self-incompatibility, which was restored by overexpressing the SCR with the reversed inversion. However, SCR in A. thaliana has other mutations aside from the pivotal inversion, in both promoter and coding regions, with probable effects on transcriptional regulation. To examine the functional consequences of these mutations, we conducted reciprocal introductions of native promoters and downstream sequences from orthologous loci of self-compatible A. thaliana and self-incompatible A. halleri. Use of this inter-species pair enabled us to expand the scope of the analysis to transcriptional regulation and deletion in the intron, in addition to inversion in the native genomic background. Initial analysis revealed that A. thaliana has a significantly lower basal expression level of SCR transcripts in the critical reproductive stage compared to that of A. halleri, suggesting that the promoter was attenuated in inducing transcription in A. thaliana. However, in reciprocal transgenic experiments, this A. thaliana promoter was able to restore partial function if coupled with the functional A. halleri coding sequence, despite extensive alterations due to the self-compatible mode of reproduction in A. thaliana. This represents a synergistic effect of the promoter and the inversion resulting in fixation of self-compatibility, primarily enforced by disruption of SCR. Our findings elucidate the functional and evolutionary context of the historical transition in A. thaliana thus contributing to the understanding of the molecular events leading to development of self-compatibility.