Comparison of the genome structure of the self-incompatibility (S) locus in interspecific pairs of S haplotypes

The determinants of recognition specificity of self-incompatibility in Brassica are SRK in the stigma and SP11/SCR in the pollen, both of which are encoded in the S locus. The nucleotide sequence analyses of many SRK and SP11/SCR alleles have identified several interspecific pairs of S haplotypes having highly similar sequences between B. oleracea and B. rapa. These interspecific pairs of S haplotypes are considered to be derived from common ancestors and to have maintained the same recognition specificity after speciation. In this study, the genome structures of three interspecific pairs of S haplotypes were compared by sequencing SRK, SP11/SCR, and their flanking regions. Regions between SRK and SP11/SCR in B. oleracea were demonstrated to be much longer than those of B. rapa and several retrotransposon-like sequences were identified in the S locus in B. oleracea. Among the seven retrotransposon-like sequences, six sequences were found to belong to the ty3 gypsy group. The gag sequences of the retrotransposon-like sequences were phylogenetically different from each other. In Southern blot analysis using retrotransposon-like sequences as probes, the B. oleracea genome showed more signals than the B. rapa genome did. These findings suggest a role for the S locus and genome evolution in self-incompatible plant species.     

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

Self-incompatibility in Brassica is controlled by a single multi-allelic locus (the S locus) which harbors at least two highly polymorphic genes, SLG and SRK. SRK is a putative transmembrane receptor kinase and its amino acid sequence of the extracellular domain of SRK (the S domain) exhibits high homology to that of SLG. The amino acid sequences of the SLGs of S8 and S46 haplotypes of B. rapa are very similar and those of S23 and S29 haplotypes of B. oleracea were also found to be almost identical. In both cases, SLG and the S domain of SRK of the same haplotype were less similar. This seems to contradict the idea that SLG and SRK of the same haplotype have the same self-recognition specificity. In the transmembrane-kinase domain, the SRK alleles of the S8 and S46 haplotypes had almost identical nucleotide sequences in spite of their lower homology in the S domain. Such a cluster of nucleotide substitutions is probably due to recombination or related events, although recombination in the S locus is thought to be suppressed. Based on our observations, the recognition mechanism and the evolution of self-incompatibility in Brassica are discussed.     

S locus genes and the evolution of self-fertility in Arabidopsis thaliana

Loss of self-incompatibility (SI) in Arabidopsis thaliana was accompanied by inactivation of genes required for SI, including S-LOCUS RECEPTOR KINASE (SRK) and S-LOCUS CYSTEINE-RICH PROTEIN (SCR), coadapted genes that constitute the SI specificity-determining S haplotype. Arabidopsis accessions are polymorphic for PsiSRK and PsiSCR, but it is unknown if the species harbors structurally different S haplotypes, either representing relics of ancestral functional and structurally heteromorphic S haplotypes or resulting from decay concomitant with or subsequent to the switch to self-fertility. We cloned and sequenced the S haplotype from C24, in which self-fertility is due solely to S locus inactivation, and show that this haplotype was produced by interhaplotypic recombination. The highly divergent organization and sequence of the C24 and Columbia-0 (Col-0) S haplotypes demonstrate that the A. thaliana S locus underwent extensive structural remodeling in conjunction with a relaxation of selective pressures that once preserved the integrity and linkage of coadapted SRK and SCR alleles. Additional evidence for this process was obtained by assaying 70 accessions for the presence of C24- or Col-0-specific sequences. Furthermore, analysis of SRK and SCR polymorphisms in these accessions argues against the occurrence of a selective sweep of a particular allele of SCR, as previously proposed.     

Interspecific pairs of class II S haplotypes having different recognition specificities between Brassica oleracea and Brassica rapa

There are several pairs of similar class I S haplotypes between Brassica oleracea and Brassica rapa. The similar S halotypes in these interspecific pairs have been reported to have the same recognition specificities. In the present study, three interspecific pairs showing a high sequence similarity were found in class II S haplotypes, i.e. between BoS-2b (B. oleracea S-2b) and BrS-44 (B. rapa S-44), between BoS-5 and BrS-40, and between BoS-15 and BrS-60. By pollination tests using interspecific hybrids between B. oleracea and B. rapa, BoS-5 and BoS-2b were revealed to have slightly and completely different recognition specificities from those of BrS-40 and BrS-44, respectively. The recognition reaction between SP11 and SRK of BoS-15 was suggested to be incomplete. The regions of class II SP11 and SRK important for self-recognition specificity and the diversification of class II S haplotypes are discussed herein.     

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.     

Intrahaplotype polymorphism at the Brassica S locus

The S locus receptor kinase and the S locus glycoproteins are encoded by genes located at the S locus, which controls the self-incompatibility response in Brassica. In class II self-incompatibility haplotypes, S locus glycoproteins can be encoded by two different genes, SLGA and SLGB. In this study, we analyzed the sequences of these genes in several independently isolated plants, all of which carry the same S haplotype (S(2)). Two groups of S(2) haplotypes could be distinguished depending on whether SRK was associated with SLGA or SLGB. Surprisingly, SRK alleles from the two groups could be distinguished at the sequence level, suggesting that recombination rarely occurs between haplotypes of the two groups. An analysis of the distribution of polymorphisms along the S domain of SRK showed that hypervariable domains I and II tend to be conserved within haplotypes but to be highly variable between haplotypes. This is consistent with these domains playing a role in the determination of haplotype specificity.     

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.     

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.     

Comparative analysis of a Brassica BAC clone containing several major aliphatic glucosinolate genes with its corresponding Arabidopsis sequence.

We compared the sequence of a 101-kb-long bacterial artificial chromosome (BAC) clone (B21H13) from Brassica oleracea with its homologous region in Arabidopsis thaliana. This clone contains a gene family involved in the synthesis of aliphatic glucosinolates. The A. thaliana homologs for this gene family are located on chromosome IV and correspond to three 2-oxoglutarate-dependent dioxygenase (AOP) genes. We found that B21H13 harbors 23 genes, whereas the equivalent region in Arabidopsis contains 37 genes. All 23 common genes have the same order and orientation in both Brassica and Arabidopsis. The 16 missing genes in the broccoli BAC clone were arranged in two major blocks of 5 and 7 contiguous genes, two singletons, and a twosome. The 118 exons comprising these 23 genes have high conservation between the two species. The arrangement of the AOP gene family in A. thaliana is as follows: AOP3 (GS-OHP) - AOP2 (GS-ALK) - pseudogene - AOP1. In contrast, in B. oleracea (broccoli and collard), two of the genes are duplicated and the third, AOP3, is missing. The remaining genes are arranged as follows: Bo-AOP2.1 (BoGSL-ALKa) - pseudogene - AOP2.2 (BoGSL-ALKb) - AOP1.1 - AOP1.2. When the survey was expanded to other Brassica accessions, we found variation in copy number and sequence for the Brassica AOP2 homologs. This study confirms that extensive rearrangements have taken place during the evolution of the Brassicacea at both gene and chromosomal levels.     

Origin and maintenance of a broad-spectrum disease resistance locus in Arabidopsis

The broad-spectrum mildew resistance genes RPW8.1 and RPW8.2 define a unique type of plant disease resistance (R) gene, and so far homologous sequences have been found in Arabidopsis thaliana only, which suggests a recent origin. In addition to RPW8.1 and RPW8.2, the RPW8 locus contains three homologs of RPW8, HR1, HR2, and HR3, which do not contribute to powdery mildew resistance. To investigate whether RPW8 has originated recently, and if so the processes involved, we have isolated and analyzed the syntenic RPW8 loci from Arabidopsis lyrata, and from Brassica rapa and B. oleracea. The A. lyrata locus contains four genes orthologous to HR1, HR2, HR3, and RPW8.2, respectively. Two syntenic loci have been characterized in Brassica; one locus contains three genes and is present in both B. oleracea and B. rapa, and the other locus contains a single gene and is detected in B. rapa only. The Brassica homologs have highest similarity to HR3. Sequence analyses suggested that the RPW8 gene family in Brassicaceae originated from an HR3-like ancestor gene through a series of duplications and that RPW8.1 and RPW8.2 evolved from functional diversification through positive selection several MYA. Examination of the sequence polymorphism of 32 A. thaliana accessions at the RPW8 locus and their disease reaction phenotypes revealed that the polymorphic RPW8 locus defines a major source of resistance to powdery mildew diseases. A possible evolutionary mechanism by which functional polymorphism at the AtRPW8 locus has been maintained in contemporary populations of A. thaliana is discussed.     

Diversification and alteration of recognition specificity of the pollen ligand SP11/SCR in self-incompatibility of Brassica and Raphanus

The recognition specificity of the pollen ligand of self-incompatibility (SP11/SCR) was investigated using Brassica rapa transgenic plants expressing SP11 transgenes, and SP11 of Raphanus sativus S-21 was found to have the same recognition specificity as that of B. rapa S-9. In a set of three S haplotypes, whose sequence identities of SP11 and SRK are fairly high, R. sativus S-6 showed the same recognition specificity as Brassica oleracea S-18 and a slightly different specificity from B. rapa S-52. B. oleracea S-18, however, showed a different specificity from B. rapa S-52. Using these similar S haplotypes, chimeric SP11 proteins were produced by domain swapping. Bioassay using the chimeric SP11 proteins revealed that the incompatibility response induction activity was altered by the replacement of Region III and Region V. Pollen grains of Brassica transgenic plants expressing chimeric SP11 of the B. oleracea SP11-18 sequence with Region III and Region V from B. rapa SP11-52 (chimeric BoSP11-18[52]) were partially incompatible with the B. rapa S-52 stigmas, and those expressing the R. sativus SP11-6 sequence with Region III and Region V from B. rapa SP11-52 (chimeric RsSP11-6[52]) were completely incompatible with the stigmas having B. rapa S-52. However, the transgenic plant expressing chimeric RsSP11-6(52) also showed incompatibility with B. oleracea S-18 stigmas. These results suggest that Regions III and Region V of SP11 are important for determining the recognition specificity, but not the sole determinant. A possible process of the generation of a new S haplotype is herein discussed.     

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.     

Genomic organization of the S core region and the S flanking regions of a class-II S haplotype in Brassica rapa

The nucleotide sequence of an 86.4-kb region that includes the SP11, SRK, and SLG genes of Brassica rapa S-60 (a class-II S haplotype) was determined. In the sequenced region, 13 putative genes were found besides SP11-60, SRK-60, and SLG-60. Five of these sequences were isolated as cDNAs, five were homologues of known genes, cDNAs, or ORFs, and three are hypothetical ORFs. Based on their nucleotide sequences, however, some of them are thought to be non-functional. Two regions of colinearity between the class-II S-60 and Brassica class-I S haplotypes were identified, i.e., S flanking region 1 which shows partial colinearity of non-genic sequences and S flanking region 2 which shows a high level of colinearity. The observed colinearity made it possible to compare the order of SP-11, SRK, and SLG genes in the S locus between the five sequenced S haplotypes. It emerged that the order of SRK and SLG in class-II S-60 is the reverse of that in the four class-I S haplotypes reported so far, and the order of SP11, SRK and SLG is the opposite of that in the class-I haplotype S-910. The possible gene designated as SAN1 (S locus Anther-expressed Non-coding RNA like-1), which is located in the region between SP11-60 and SRK-60, has features reminiscent of genes for non-coding RNAs (ncRNAs), but no homologous sequences were found in the databases. This sequence is transcribed in anthers but not in stigmas or leaves. These features of the genomic structure of S-60 are discussed with special reference to the characteristics of class-II S haplotypes.     

Regulation of sulfate uptake and expression of sulfate transporter genes in Brassica oleracea as affected by atmospheric H(2)S and pedospheric sulfate nutrition

Demand-driven signaling will contribute to regulation of sulfur acquisition and distribution within the plant. To investigate the regulatory mechanisms pedospheric sulfate and atmospheric H(2)S supply were manipulated in Brassica oleracea. Sulfate deprivation of B. oleracea seedlings induced a rapid increase of the sulfate uptake capacity by the roots, accompanied by an increased expression of genes encoding specific sulfate transporters in roots and other plant parts. More prolonged sulfate deprivation resulted in an altered shoot-root partitioning of biomass in favor of the root. B. oleracea was able to utilize atmospheric H(2)S as S-source; however, root proliferation and increased sulfate transporter expression occurred as in S-deficient plants. It was evident that in B. oleracea there was a poor shoot to root signaling for the regulation of sulfate uptake and expression of the sulfate transporters. cDNAs corresponding to 12 different sulfate transporter genes representing the complete gene family were isolated from Brassica napus and B. oleracea species. The sequence analysis classified the Brassica sulfate transporter genes into four different groups. The expression of the different sulfate transporters showed a complex pattern of tissue specificity and regulation by sulfur nutritional status. The sulfate transporter genes of Groups 1, 2, and 4 were induced or up-regulated under sulfate deprivation, although the expression of Group 3 sulfate transporters was not affected by the sulfate status. The significance of sulfate, thiols, and O-acetylserine as possible signal compounds in the regulation of the sulfate uptake and expression of the transporter genes is evaluated.     

The S15 self-incompatibility haplotype in Brassica oleracea includes three S gene family members expressed in stigmas.

Self-incompatibility in Brassica is controlled by a single, highly polymorphic locus that extends over several hundred kilobases and includes several expressed genes. Two stigma proteins, the S locus receptor kinase (SRK) and the S locus glycoprotein (SLG), are encoded by genes located at the S locus and are thought to be involved in the recognition of self-pollen by the stigma. We report here that two different SLG genes, SLGA and SLGB, are located at the S locus in the class II, pollen-recessive S15 haplotype. Both genes are interrupted by a single intron; however, SLGA encodes both soluble and membrane-anchored forms of SLG, whereas SLGB encodes only soluble SLG proteins. Thus, including SRK, the S locus in the S15 haplotype contains at least three members of the S gene family. The protein products of these three genes have been characterized, and each SLG glycoform was assigned to an SLG gene. Evidence is presented that the S2 and S5 haplotypes carry only one or the other of the SLG genes, indicating either that they are redundant or that they are not required for the self-incompatibility response.     

FISH-mapping of rDNAs and Arabidopsis BACs on pachytene complements of selected Brassicas.

To improve resolution of physical mapping on Brassica chromosomes, we have chosen the pachytene stage of meiosis where incompletely condensed bivalents are much longer than their counterparts at mitotic metaphase. Mapping with 5S and 45S rDNA sequences demonstrated the advantage of pachytene chromosomes in efficient physical mapping and confirmed the presence of a novel 5S rDNA locus in Brassica oleracea, initially identified by genetic mapping using restriction fragment length polymorphism (RFLP). Fluorescence in situ hybridization (FISH) analysis visualized the presence of the third 5S rDNA locus on the long arm of chromosome C2 and confirmed the earlier reports of two 45S rDNA loci in the B. oleracea genome. FISH mapping of low-copy sequences from the Arabidopsis thaliana bacterial artificial chromosome (BAC) clones on the B. oleracea chromosomes confirmed the expectation of efficient and precise physical mapping of meiotic bivalents based on data available from A. thaliana and indicated conserved organization of these two BAC sequences on two B. oleracea chromosomes. Based on the heterologous in situ hybridization with BACs and their mapping applied to long pachytene bivalents, a new approach in comparative analysis of Brassica and A. thaliana genomes is discussed.