Comparative analysis of S haplotypes with very similar SLG alleles in Brassica rapa and Brassica oleracea

Development of S Haplotype-Specific Markers to Identify Genotypes of Self-Incompatibility in Radish (Raphanus sativus L.)

Radish (Raphanus sativus L.), a root vegetable belonging to the Brassicaceae family, is considered one of the representative crops displaying sporophytic self-incompatibility (SSI). The utilization of a self-incompatibility system in F1 breeding can improve the efficiency of cross-combinations, leading to a reduction in breeding time and aiding in the development of novel F1 varieties. The successful implementation of this system necessitates the rapid and accurate identification of S haplotypes in parental lines. In this study, we identified a total of nine S haplotypes among 22 elite radish lines through Sanger sequencing. Subsequently, we obtained sequences for showing a 95% similarity to nine S haplotypes, along with sequences identified by other researchers using BLAST. Following this, multiple sequence alignment (MSA) was conducted to identify SRK and SLG sequence similarities, as well as polymorphisms within the class I and II groups. Subsequently, S haplotype-specific marker sets were developed, targeting polymorphic regions of SRK and SLG alleles. These markers successfully amplified each of the nine S haplotypes. These markers will play a crucial role in the rapid and precise identification of parental S haplotypes in the radish F1 breeding process, proving instrumental in the radish F1 purity test.

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.

S haplotype collection in Brassicaceae crops-an updated list of S haplotypes

Self-incompatibility is the system that inhibits pollen germination and pollen tube growth by self-pollen. This trait is important for the breeding of Brassica and Raphanus species. In these species, self-incompatibility is governed by the S locus, which contains three linked genes (a set called the S haplotype), i.e., S-locus receptor kinase, S-locus cysteine-rich protein/S-locus protein 11, and S-locus glycoprotein. A large number of S haplotypes have been identified in Brassica oleracea, B. rapa, and Raphanus sativus to date, and the nucleotide sequences of their many alleles have also been registered. In this state, it is important to avoid confusion between S haplotypes, i.e., an identical S haplotype with different names and a different S haplotype with an identical S haplotype number. To mitigate this issue, we herein constructed a list of S haplotypes that are easily accessible to the latest nucleotide sequences of S-haplotype genes, together with revisions to and an update of S haplotype information. Furthermore, the histories of the S-haplotype collection in the three species are reviewed, the importance of the collection of S haplotypes as a genetic resource is discussed, and the management of information on S haplotypes is proposed.

Analysis of S-locus and expression of S-alleles of self-compatible rapid-cycling Brassica oleracea 'TO1000DH3'

Brassica oleracea is a strictly self-incompatible (SI) plant, but rapid-cycling B. oleracea 'TO1000DH3' is self-compatible (SC). Self-incompatibility in Brassicaceae is controlled by multiple alleles of the S-locus. Three S-locus genes, S-locus glycoprotein (SLG), S-locus receptor kinase (SRK) and S-locus protein 11 or S-locus cysteine-rich (SP11/SCR), have been reported to date, all of which are classified into class I and II. In this study, we investigated the molecular mechanism behind alterations of SI to SC in rapid-cycling B. olerace 'TO1000DH3'. Class I SRK were identified by genomic DNA PCR and PCR-RFLP analysis using SRK specific markers and found to be homozygous. Cloning and sequencing of class I SRK revealed a normal kinase domain without any S-domain/transmembrane domain. Moreover, S-locus sequencing analysis revealed only an SLG sequence, but no SP11/SCR. Expression analysis showed no SRK expression in the stigma, although other genes involved in the SI recognition reaction (SLG, MLPK, ARC1, THL) were found to have normal expression in the stigma. Taken together, the above results suggest that structural aberrations such as deletion of the SI recognition genes may be responsible for the breakdown of SI in rapid-cycling B. oleracea 'TO1000DH3'.

Self-incompatibility in Brassicaceae crops: lessons for interspecific incompatibility

Most wild plants and some crops of the Brassicaceae express self-incompatibility, which is a mechanism that allows stigmas to recognize and discriminate against "self" pollen, thus preventing self-fertilization and inbreeding. Self-incompatibility in this family is controlled by a single S locus containing two multiallelic genes that encode the stigma-expressed S-locus receptor kinase and its pollen coat-localized ligand, the S-locus cysteine-rich protein. Physical interaction between receptor and ligand encoded in the same S locus activates the receptor and triggers a signaling cascade that results in inhibition of "self" pollen. Sequence information for many S-locus haplotypes in Brassica species has spurred studies of dominance relationships between S haplotypes and of S-locus structure, as well as the development of methods for S genotyping. Furthermore, molecular genetic studies have begun to identify genes that encode putative components of the self-incompatibility signaling pathway. In parallel, standard genetic analysis and QTL analysis of the poorly understood interspecific incompatibility phenomenon have been initiated to identify genes responsible for the inhibition of pollen from other species by the stigma. Herewith, we review recent studies of self-incompatibility and interspecific incompatibility, and we propose a model in which a universal pollen-inhibition pathway is shared by these two incompatibility systems.

Distribution of S haplotypes and its relationship with restorer-maintainers of self-incompatibility in cultivated Brassica napus

Brassica napus (AACC, 2n = 38) is a self-compatible amphidiploid plant that arose from the interspecies hybridization of two self-incompatible species, B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18). Self-incompatibility (S) haplotypes in one self-incompatible line and 124 cultivated B. napus lines were detected using S-locus-specific primers, and their relationships with restorer-maintainers were investigated. Two class I (S-I ( SLG ) a and S-I ( SLG ) b) and four class II (S-II ( SLG ) a, S-II ( SLG ) b, S-II ( SP11 ) a and S-II ( SP11 ) b) S haplotypes were observed, of which S-II ( SP11 ) b was newly identified. The nucleotide sequence of SP11 showed little similarity to the reported SP11 alleles. The lines were found to express a total of eleven S genotypes. The self-incompatible line had a specific genotype consisting of S-II ( SP11 ) a, similar to B. rapa S-60, and S-II ( SLG ) a, similar to B. oleracea S-15. Restorers expressed six genotypes: the most common genotype contained S-I ( SLG ) a, similar to B. rapa S-47, and S-II ( SLG ) b, similar to B. oleracea S-15. Maintainers expressed nine genotypes: the predominant genotype was homozygous for two S haplotypes, S-II ( SLG ) a and S-II ( SP11 ) b. One genotype was specific to restorers and four genotypes were specific to maintainers, whereas five genotypes were expressed in both restorers and maintainers. This suggests that there is no definitive correlation between the distribution of S genotypes and restorer-maintainers of self-incompatibility. The finding that restorers and maintainers express unique genotypes, and share some common genotypes, would be valuable for detecting the interaction of S haplotypes in inter- or intra-genomes as well as for developing markers-assisted selection in self-incompatibility hybrid breeding.

Effects of recombination on hitchhiking diversity in the Brassica self-incompatibility locus complex

In self-incompatibility, a number of S haplotypes are maintained by frequency-dependent selection, which results in trans-specific S haplotypes. The region of several kilobases (approximately 40-60 kb) from SP6 to SP2, including self-incompatibility-related genes and some adjacent genes in Brassica rapa, has high nucleotide diversity due to the hitchhiking effect, and therefore we call this region the "S-locus complex." Recombination in the S-locus complex is considered to be suppressed. We sequenced regions of >50 kb of the S-locus complex of three S haplotypes in B. rapa and found higher nucleotide diversity in intergenic regions than in coding regions. Two highly similar regions of >10 kb were found between BrS-8 and BrS-46. Phylogenetic analysis using trans-specific S haplotypes (called interspecific pairs) of B. rapa and B. oleracea suggested that recombination reduced the nucleotide diversity in these two regions and that the genes not involved in self-incompatibility in the S-locus complex and the kinase domain, but not the S domain, of SRK have also experienced recombination. Recombination may reduce hitchhiking diversity in the S-locus complex, whereas the region from the S domain to SP11 would disfavor recombination.

Evolution under strong balancing selection: how many codons determine specificity at the female self-incompatibility gene SRK in Brassicaceae?

Background

Molecular lock-and-key systems are common among reproductive proteins, yet their evolution remains a major puzzle in evolutionary biology. In the Brassicaceae, the genes encoding self-incompatibility have been identified, but technical challenges currently prevent detailed analyses of the molecular interaction between the male and female components. In the present study, we investigate sequence polymorphism in the female specificity determinant SRK of Arabidopsis halleri from throughout Europe. Using a comparative approach based on published SRK sequences in A. lyrata and Brassica, we track the signature of frequency-dependent selection acting on these genes at the codon level. Using simulations, we evaluate power and accuracy of our approach and estimate the proportion of codon sites involved in the molecular interaction.

Pollen recognition and rejection during the sporophytic self-incompatibility response: Brassica and beyond

Many hermaphrodite flowering plants avoid self-fertilization through genetic systems of self-incompatibility (SI). SI allows a plant to recognize and to reject self or self-related pollen, thereby preserving its ovules for outcrossing. Genes situated at the S-locus encode the 'male' (pollen) and 'female' (pistil) recognition determinants of SI. In sporophytic SI (SSI) the male determinant is expressed in the diploid anther, therefore haploid pollen grains behave with a diploid S phenotype. In Brassica, the male and the female determinants of SSI have been identified as a peptide ligand and its cognate receptor, respectively, and recent studies have identified downstream signalling molecules involved in pollen rejection. It now needs to be established whether the Brassica mechanism is universal in species with SSI, or unique to the Brassicaceae.

The S haplotypes lacking SLG in the genome of Brassica rapa

Self-incompatibility (SI) discriminating self- and non-self pollen is regulated by S-locus genes in Brassica. In most of the S haplotypes, a highly polymorphic S-locus glycoprotein ( SLG) gene is tightly linked to genes for the SI determinants, S-receptor kinase ( SRK) and SP11, although the precise function of SLG in SI has not been clarified. In the present study, we performed DNA gel blot analysis for S(32), S(33), and S(36) haplotypes of Brassica rapa showing normal SI phenotypes and concluded that there might be no SLG in their genome. RNA gel blot analysis of the SLG-less S haplotypes indicated the possible existence of eSRK transcripts in the stigma. These three S haplotypes are useful resources to discern the molecular mechanism of the SI reaction without SLG.

Coevolution of the S-locus genes SRK, SLG and SP11/SCR in Brassica oleracea and B. rapa

Brassica self-incompatibility (SI) is controlled by SLG and SRK expressed in the stigma and by SP11/SCR expressed in the anther. We determined the sequences of the S domains of 36 SRK alleles, 13 SLG alleles, and 14 SP11 alleles from Brassica oleracea and B. rapa. We found three S haplotypes lacking SLG genes in B. rapa, confirming that SLG is not essential for the SI recognition system. Together with reported sequences, the nucleotide diversities per synonymous and nonsynonymous site (pi(S) and pi(N)) at the SRK, SLG, and SP11 loci within B. oleracea were computed. The ratios of pi(N):pi(S) for SP11 and the hypervariable region of SRK were significantly >1, suggesting operation of diversifying selection to maintain the diversity of these regions. In the phylogenetic trees of 12 SP11 sequences and their linked SRK alleles, the tree topology was not significantly different between SP11 and SRK, suggesting a tight linkage of male and female SI determinants during the evolutionary course of these haplotypes. Genetic exchanges between SLG and SRK seem to be frequent; three such recent exchanges were detected. The evolution of S haplotypes and the effect of gene conversion on self-incompatibility are discussed.

Recognition specificity of self-incompatibility maintained after the divergence of Brassica oleracea and Brassica rapa

The determinants of recognition specificity of self-incompatibility in Brassica are SRK in the stigma and SP11/SCR in the pollen, respectively. In the pair of S haplotypes BrS46 (S46 in B. rapa) and BoS7 (S7 in B. oleracea), which have highly similar SRK alleles, the SP11 alleles were found to be similar, with 96.1% identity in the deduced amino acid sequence. Two other pairs of S haplotypes, BrS47 and BoS12, and BrS8 and BoS32, having highly similar SRK and SP11 alleles between the two species were also found. The haplotypes in each pair are considered to have been derived from a single S haplotype in the ancestral species. The allotetraploid produced by interspecific hybridization between homozygotes of BrS46 and BoS15 showed incompatibility with a BoS7 homozygote and compatibility with other B. oleracea S haplotypes in reciprocal crossings. This result indicates that BrS46 and BoS7 have maintained the same recognition specificity after the divergence of the two species and that amino acid substitutions found in such cases in both SRK alleles and SP11 alleles do not alter the recognition specificity. DNA blot analysis of SRK, SP11, SLG and other S-locus genes showed different DNA fragment sizes between the interspecific pairs of S haplotypes. A much lower level of sequence similarity was observed outside the genes of SRK and SP11 between BrS46 and BoS7. These results suggest that the DNA sequences of the regions intervening between the S-locus genes were diversified after or at the time of speciation. This is the first report demonstrating the presence of common S haplotypes in different plant species and presenting definite evidence of the trans-specific evolution of self-incompatibility genes.

Allele-specific receptor-ligand interactions in Brassica self-incompatibility

Genetic self-incompatibility in Brassica is determined by alleles of the transmembrane serine-threonine kinase SRK, which functions in the stigma epidermis, and of the cysteine-rich peptide SCR, which functions in pollen. Using tagged versions of SRK and SCR as well as endogenous stigma and pollen proteins, we show that SCR binds the SRK ectodomain and that this binding is allele specific. Thus, SRK and SCR function as a receptor-ligand pair in the recognition of self pollen. Specificity in the self-incompatibility response derives from allele-specific formation of SRK-SCR complexes at the pollen-stigma interface.

Molecular recognition and response in pollen and pistil interactions

Many bisexual flowering plants possess a reproductive strategy called self-incompatibility (SI) that enables the female tissue (the pistil) to reject self but accept non-self pollen for fertilization. Three different SI mechanisms are discussed, each controlled by two separate, highly polymorphic genes at the S-locus. For the Solanaceae and Papaveraceae types, the genes controlling female function in SI, the S-RNase gene and the S-gene, respectively, have been identified. For the Brassicaceae type, the gene controlling male function, SCR/SP11, and the gene controlling female function, SRK, have been identified. The S-RNase based mechanism involves degradation of RNA of self-pollen tubes; the S-protein based mechanism involves a signal transduction cascade in pollen, including a transient rise in [Ca(2+)]i and subsequent protein phosphorylation/dephosphorylation; and the SRK (a receptor kinase) based mechanism involves interaction of a pollen ligand, SCR/SP11, with SRK, followed by a signal transduction cascade in the stigmatic surface cell.

Cell-cell signaling in the self-incompatibility response

Significant progress towards understanding the molecular basis of self recognition in the self-incompatibility response of Brassica has been made during the past two years. The highly polymorphic molecules that determine the specificity of this interaction in the pollen and stigma have been identified. The structural features of these molecules suggest that a ligand-receptor-kinase interaction triggers the initiation of a signaling cascade within the stigma epidermis and the subsequent arrest of self-pollination.

Characterization of Brassica S-haplotypes lacking S-locus glycoprotein

Self-incompatibility (SI) in Brassica is regulated by a single multi-allelic locus, S, which contains highly polymorphic stigma-expressed genes, SLG and SRK. While SRK is shown to be the determinant of female SI specificity, SLG is thought to assist the function of SRK. Here we report that the SLG genes of self-incompatible S(18) and S(60) homozygotes of Brassica oleracea have an in-frame stop codon and a 23 bp deletion resulting in a frame-shift, respectively. The finding that these SLG genes do not encode functional SLG proteins suggests that SLG is not essential for SI. The possible role of SLG in SI was discussed.

Post-transcriptional maturation of the S receptor kinase of Brassica correlates with co-expression of the S-locus glycoprotein in the stigmas of two Brassica strains and in transgenic tobacco plants

The S-locus-encoded S receptor kinase (SRK) is an intrinsic plasma membrane protein that is viewed as the primary stigma determinant of specificity in the self-incompatibility response of Brassica spp. We analyzed two self-compatible mutant strains that express low levels of the S-locus glycoprotein (SLG), a cell wall-localized protein also encoded at the S locus that is coordinately expressed with SRK. We found that mutant stigmas synthesized wild-type levels of SRK transcripts but failed to produce SRK protein at any of the developmental stages analyzed. Furthermore, SRK was shown to form aberrant high-molecular mass aggregates when expressed alone in transgenic tobacco (Nicotiana tabacum) plants. This aggregation was prevented in tobacco plants that co-expressed SRK and SLG, but not in tobacco plants that co-expressed SRK and SLR1, an SLG-related secreted protein not encoded at the S locus. In analyses of protein extracts under reducing and non-reducing conditions, evidence of intermolecular association was obtained only for SLG, a fraction of which formed disulfide-linked oligomers and was membrane associated. The data indicate that, at least in plants carrying the S haplotypes we analyzed, SRK is an inherently unstable protein and that SLG facilitates its accumulation to physiologically relevant levels in Brassica stigmas.

Evolutionary dynamics of self-incompatibility alleles in Brassica

Self-incompatibility in Brassica entails the rejection of pollen grains that express specificities held in common with the seed parent. In Brassica, pollen specificity is encoded at the multipartite S-locus, a complex region comprising many expressed genes. A number of species within the Brassicaceae express sporophytic self-incompatibility, under which individual pollen grains bear specificities determined by one or both S-haplotypes of the pollen parent. Classical genetic and nucleotide-level analyses of the S-locus have revealed a dichotomy in sequence and function among S-haplotypes; in particular, all class I haplotypes show dominance over all class II haplotypes in determination of pollen specificity. Analysis of an evolutionary model that explicitly incorporates features of the Brassica system, including the class dichotomy, indicates that class II haplotypes may invade populations at lower rates and decline to extinction at higher rates than class I haplotypes. This analysis suggests convergence to an evolutionarily persistent state characterized by the maintenance in high frequency of a single class II haplotype together with many class I haplotypes, each in low frequency. This expectation appears to be consistent with empirical observations of high frequencies of relatively few distinct recessive haplotypes.

Sequence and structural diversity of the S locus genes from different lines with the same self-recognition specificities in Brassica oleracea

Self-incompatibility (SI) is a mechanism for preventing self-fertilization in flowering plants. In Brassica, it is controlled by a single multi-allelic locus, S, and it is believed that two highly polymorphic genes in the S locus, SLG and SRK, play central roles in self-recognition in stigmas. SRK is a putative receptor protein kinase, whose extracellular domain exhibits high similarity to SLG. We analyzed two pairs of lines showing cross-incompatibility (S(2) and S(2-b); S(13) and S(13-b)). In S(2) and S(2-b), SRKs were more highly conserved than SLGs. This was also the case with S(13) and S(13-b). This suggests that the SRKs of different lines must be conserved for the lines to have the same self-recognition specificity. In particular, SLG(2-b) showed only 88. 5% identity to SLG(2), which is comparable to that between the SLGs of different S haplotypes, while SRK(2-b) showed 97.3% identity to SRK(2) in the S domain. These findings suggest that the SLGs in these S haplotypes are not important for self-recognition in SI.