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

Transposable elements cause the loss of self-incompatibility in citrus

Self-incompatibility (SI) is a widespread prezygotic mechanism for flowering plants to avoid inbreeding depression and promote genetic diversity. Citrus has an S-RNase-based SI system, which was frequently lost during evolution. We previously identified a single nucleotide mutation in Sm-RNase, which is responsible for the loss of SI in mandarin and its hybrids. However, little is known about other mechanisms responsible for conversion of SI to self-compatibility (SC) and we identify a completely different mechanism widely utilized by citrus. Here, we found a 786-bp miniature inverted-repeat transposable element (MITE) insertion in the promoter region of the FhiS2-RNase in Fortunella hindsii Swingle (a model plant for citrus gene function), which does not contain the Sm-RNase allele but are still SC. We demonstrate that this MITE plays a pivotal role in the loss of SI in citrus, providing evidence that this MITE insertion prevents expression of the S-RNase; moreover, transgenic experiments show that deletion of this 786-bp MITE insertion recovers the expression of FhiS2-RNase and restores SI. This study identifies the first evidence for a role for MITEs at the S-locus affecting the SI phenotype. A family-wide survey of the S-locus revealed that MITE insertions occur frequently adjacent to S-RNase alleles in different citrus genera, but only certain MITEs appear to be responsible for the loss of SI. Our study provides evidence that insertion of MITEs into a promoter region can alter a breeding strategy and suggests that this phenomenon may be broadly responsible for SC in species with the S-RNase system.

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.

Whole-Genome DNA Methylation Analysis in Brassica rapa subsp. perviridis in Response to Albugo candida Infection

DNA methylation is an epigenetic mark associated with several mechanisms in plants including immunity mechanisms. However, little is known about the regulatory role of DNA methylation in the resistance response of Brassica species against fungal diseases. White rust, caused by the fungus Albugo candida, is one of the most widespread and destructive diseases of all the cultivated Brassica species, particularly Brassica rapa L. and Brassica juncea (L.) Czern and Coss. Here, we investigate whole-genome DNA methylation modifications of B. rapa subsp. perviridis in response to white rust. As a result, 233 and 275 differentially methylated regions (DMRs) in the susceptible cultivar "Misugi" and the resistant cultivar "Nanane" were identified, respectively. In both cultivars, more than half of the DMRs were associated with genes (DMR-genes). Gene expression analysis showed that 13 of these genes were also differentially expressed between control and infected samples. Gene ontology enrichment analysis of DMR genes revealed their involvement in various biological processes including defense mechanisms. DMRs were unevenly distributed around genes in susceptible and resistant cultivars. In "Misugi," DMRs tended to be located within genes, while in "Nanane," DMRs tended to be located up and downstream of the genes. However, CG DMRs were predominantly located within genes in both cultivars. Transposable elements also showed association with all three sequence contexts of DMRs but predominantly with CHG and CHH DMRs in both cultivars. Our findings indicate the occurrence of DNA methylation modifications in B. rapa in response to white rust infection and suggest a potential regulatory role of DNA methylation modification in defense mechanisms which could be exploited to improve disease resistance.

Genome Triplication Leads to Transcriptional Divergence of FLOWERING LOCUS C Genes During Vernalization in the Genus Brassica

The genus Brassica includes oil crops, vegetables, condiments, fodder crops, and ornamental plants. Brassica species underwent a whole genome triplication event after speciation between ancestral species of Brassica and closely related genera including Arabidopsis thaliana. Diploid species such as Brassica rapa and Brassica oleracea have three copies of genes orthologous to each A. thaliana gene, although deletion in one or two of the three homologs has occurred in some genes. The floral transition is one of the crucial events in a plant's life history, and time of flowering is an important agricultural trait. There is a variation in flowering time within species of the genus Brassica, and this variation is largely dependent on a difference in vernalization requirements. In Brassica, like in A. thaliana, the key gene of vernalization is FLOWERING LOCUS C (FLC). In Brassica species, the vernalization response including the repression of FLC expression by cold treatment and the enrichment of the repressive histone modification tri-methylated histone H3 lysine 27 (H3K27me3) at the FLC locus is similar to A. thaliana. B. rapa and B. oleracea each have four paralogs of FLC, and the allotetraploid species, Brassica napus, has nine paralogs. The increased number of paralogs makes the role of FLC in vernalization more complicated; in a single plant, paralogs vary in the expression level of FLC before and after vernalization. There is also variation in FLC expression levels between accessions. In this review, we focus on the regulatory circuits of the vernalization response of FLC expression in the genus Brassica.

Plant mating systems: self-incompatibility and evolutionary transitions to self-fertility in the mustard family

Flowering plants have evolved diverse mechanisms that promote outcrossing. The most widespread of these outbreeding devices are self-incompatibility systems, the highly selective prefertilization mating barriers that prevent self-fertilization by disrupting pollen-pistil interactions. Despite the advantages of outcrossing, loss of self-incompatibility has occurred repeatedly in many plant families. In the mustard family, the highly polymorphic receptors and ligands that mediate the recognition and inhibition of self-pollen in self-incompatibility have been characterized and the 3D structure of the receptor-ligand complex has been solved. Sequence analyses and empirical studies in self-incompatible and self-compatible species are elucidating the genetic basis of switches from the outcrossing to selfing modes of mating and beginning to provide clues to the diversification of the self recognition repertoire.

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.

Trait transitions in explicit ecological and genomic contexts: plant mating systems as case studies

Plants are astonishingly diverse in how they reproduce sexually, and the study of plant mating systems provides some of the most compelling cases of parallel and independent evolutionary transitions. In this chapter, we review how the massive amount of genomic data being produced is allowing long-standing predictions from ecological and evolutionary theory to be put to test. After a review of theoretical predictions about the importance of considering the genomic architecture of the mating system, we focus on a set of recent discoveries on how the mating system is controlled in a variety of model and non-model species. In parallel, genomic approaches have revealed the complex interaction between the evolution of genes controlling mating systems and genome evolution, both genome-wide and in the mating system control region. In several cases, major transitions in the mating system can be clearly associated with important ecological changes, hence illuminating an important interplay between ecological and genomic approaches. We also list a number of major unsolved questions that remain for the field, and highlight foreseeable conceptual developments that are likely to play a major role in our understanding of how plant mating systems evolve in Nature.

Self-(in)compatibility inheritance and allele-specific marker development in yellow mustard (Sinapis alba)

Yellow mustard (Sinapis alba) has a sporophytic self-incompatibility reproduction system. Genetically stable self-incompatible (SI) and self-compatible (SC) inbred lines have recently been developed in this crop. Understanding the S haplotype of different inbred lines and the inheritance of the self-(in)compatibility (SI/SC) trait is very important for breeding purposes. In this study, we used the S-locus gene-specific primers in Brassica rapa and Brassica oleracea to clone yellow mustard S-locus genes of SI lines Y514 and Y1130 and SC lines Y1499 and Y1501. The PCR amplification results and DNA sequences of the S-locus genes revealed that Y514 carried the class I S haplotype, while Y1130, Y1499, and Y1501 had the class II S haplotype. The results of our genetic studies indicated that self-incompatibility was dominant over self-compatibility and controlled by a one-gene locus in the two crosses of Y514 × Y1499 and Y1130 × Y1501. Of the five S-locus gene polymorphic primer pairs, Sal-SLGI and Sal-SRKI each generated one dominant marker for the SI phenotype of Y514; Sal-SLGII and Sal-SRKII produced dominant marker(s) for the SC phenotype of Y1501 and Y1499; Sal-SP11II generated one dominant marker for Y1130. These markers co-segregated with the SI/SC phenotype in the F₂ populations of the two crosses. In addition, co-dominant markers were developed by mixing the two polymorphic primer pairs specific for each parent in the multiplex PCR, which allowed zygosity to be determined in the F₂ populations. The SI/SC allele-specific markers have proven to be very useful for the selection of the desirable SC genotypes in our yellow mustard breeding program.

The genome of the pear (Pyrus bretschneideri Rehd.)

The draft genome of the pear (Pyrus bretschneideri) using a combination of BAC-by-BAC and next-generation sequencing is reported. A 512.0-Mb sequence corresponding to 97.1% of the estimated genome size of this highly heterozygous species is assembled with 194× coverage. High-density genetic maps comprising 2005 SNP markers anchored 75.5% of the sequence to all 17 chromosomes. The pear genome encodes 42,812 protein-coding genes, and of these, ∼28.5% encode multiple isoforms. Repetitive sequences of 271.9 Mb in length, accounting for 53.1% of the pear genome, are identified. Simulation of eudicots to the ancestor of Rosaceae has reconstructed nine ancestral chromosomes. Pear and apple diverged from each other ∼5.4–21.5 million years ago, and a recent whole-genome duplication (WGD) event must have occurred 30–45 MYA prior to their divergence, but following divergence from strawberry. When compared with the apple genome sequence, size differences between the apple and pear genomes are confirmed mainly due to the presence of repetitive sequences predominantly contributed by transposable elements (TEs), while genic regions are similar in both species. Genes critical for self-incompatibility, lignified stone cells (a unique feature of pear fruit), sorbitol metabolism, and volatile compounds of fruit have also been identified. Multiple candidate SFB genes appear as tandem repeats in the S-locus region of pear; while lignin synthesis-related gene family expansion and highly expressed gene families of HCT, C3′H, and CCOMT contribute to high accumulation of both G-lignin and S-lignin. Moreover, alpha-linolenic acid metabolism is a key pathway for aroma in pear fruit.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this genomic region.

The self-compatibility mechanism in Brassica napus L. is applicable to F1 hybrid breeding

Brassica napus, an allopolyploid species having the A genome of B. rapa and the C genome of B. oleracea, is self-compatible, although both B. rapa and B. oleracea are self-incompatible. We have previously reported that SP11/SCR alleles are not expressed in anthers, while SRK alleles are functional in the stigma in B. napus cv. 'Westar', which has BnS-1 similar to B. rapa S-47 and BnS-6 similar to B. oleracea S-15. This genotype is the most frequent S genotype in B. napus, and we hypothesized that the loss of the function of SP11 is the primary cause of the self-compatibility of 'Westar'. To verify this hypothesis, we transformed 'Westar' plants with the SP11 allele of B. rapa S-47. All the transgenic plants and their progeny were completely self-incompatible, demonstrating self-compatibility to be due to the S haplotype having the non-functional SP11 allele in the A genome, which suppresses a functional recessive SP11 allele in the C genome. An artificially synthesized B. napus line having two recessive SP11 alleles was developed by interspecific hybridization between B. rapa and B. oleracea. This line was self-incompatible, but F(1) hybrids between this line and 'Westar' were self-compatible. These results suggest that the self-compatibility mechanism of 'Westar' is applicable to F(1) seed production in B. napus.

Analysis of target sequences of DDM1s in Brassica rapa by MSAP

DNA methylation is an important epigenetic modification regulating gene expression and transposon silencing. Although epigenetic regulation is involved in some agricultural traits, there has been relatively little research on epigenetic modifications of genes in Brassica rapa, which includes many important vegetables. In B. rapa, orthologs of DDM1, a chromatin remodeling factor required for maintenance of DNA methylation, have been characterized and DNA hypomethylated knock-down plants by RNAi (ddm1-RNAi plants) have been generated. In this study, we investigated differences of DNA methylation status at the genome-wide level between a wild-type (WT) plant and a ddm1-RNAi plant by methylation-sensitive amplification polymorphism (MSAP) analysis. MSAP analysis detected changes of DNA methylation of many repetitive sequences in the ddm1-RNAi plant. Search for body methylated regions in the WT plant revealed no difference in gene body methylation levels between the WT plant and the ddm1-RNAi plant. These results indicate that repetitive sequences are preferentially methylated by DDM1 genes in B. rapa.

Genomic microstructure and differential expression of the genes encoding UDP-glucose:sinapate glucosyltransferase (UGT84A9) in oilseed rape (Brassica napus)

In oilseed rape (Brassica napus), the glucosyltransferase UGT84A9 catalyzes the formation of 1-O-sinapoyl-beta-glucose, which feeds as acyl donor into a broad range of accumulating sinapate esters, including the major antinutritive seed component sinapoylcholine (sinapine). Since down-regulation of UGT84A9 was highly efficient in decreasing the sinapate ester content, the genes encoding this enzyme were considered as potential targets for molecular breeding of low sinapine oilseed rape. B. napus harbors two distinguishable sequence types of the UGT84A9 gene designated as UGT84A9-1 and UGT84A9-2. UGT84A9-1 is the predominantly expressed variant, which is significantly up-regulated during the seed filling phase, when sinapate ester biosynthesis exhibits strongest activity. In the allotetraploid genome of B. napus, UGT84A9-1 is represented by two loci, one derived from the Brassica C-genome (UGT84A9a) and one from the Brassica A-genome (UGT84A9b). Likewise, for UGT84A9-2 two loci were identified in B. napus originating from both diploid ancestor genomes (UGT84A9c, Brassica C-genome; UGT84A9d, Brassica A-genome). The distinct UGT84A9 loci were genetically mapped to linkage groups N15 (UGT84A9a), N05 (UGT84A9b), N11 (UGT84A9c) and N01 (UGT84A9d). All four UGT84A9 genomic loci from B. napus display a remarkably low micro-collinearity with the homologous genomic region of Arabidopsis thaliana chromosome III, but exhibit a high density of transposon-derived sequence elements. Expression patterns indicate that the orthologous genes UGT84A9a and UGT84A9b should be considered for mutagenesis inactivation to introduce the low sinapine trait into oilseed rape.

Comparative analyses reveal distinct sets of lineage-specific genes within Arabidopsis thaliana

Background

The availability of genome and transcriptome sequences for a number of species permits the identification and characterization of conserved as well as divergent genes such as lineage-specific genes which have no detectable sequence similarity to genes from other lineages. While genes conserved among taxa provide insight into the core processes among species, lineage-specific genes provide insights into evolutionary processes and biological functions that are likely clade or species specific.

Results

Comparative analyses using the Arabidopsis thaliana genome and sequences from 178 other species within the Plant Kingdom enabled the identification of 24,624 A. thaliana genes (91.7%) that were termed Evolutionary Conserved (EC) as defined by sequence similarity to a database entry as well as two sets of lineage-specific genes within A. thaliana. One of the A. thaliana lineage-specific gene sets share sequence similarity only to sequences from species within the Brassicaceae family and are termed Conserved Brassicaceae-Specific Genes (914, 3.4%, CBSG). The other set of A. thaliana lineage-specific genes, the Arabidopsis Lineage-Specific Genes (1,324, 4.9%, ALSG), lack sequence similarity to any sequence outside A. thaliana. While many CBSGs (76.7%) and ALSGs (52.9%) are transcribed, the majority of the CBSGs (76.1%) and ALSGs (94.4%) have no annotated function. Co-expression analysis indicated significant enrichment of the CBSGs and ALSGs in multiple functional categories suggesting their involvement in a wide range of biological functions. Subcellular localization prediction revealed that the CBSGs were significantly enriched in proteins targeted to the secretory pathway (412, 45.1%). Among the 107 putatively secreted CBSGs with known functions, 67 encode a putative pollen coat protein or cysteine-rich protein with sequence similarity to the S-locus cysteine-rich protein that is the pollen determinant controlling allele specific pollen rejection in self-incompatible Brassicaceae species. Overall, the ALSGs and CBSGs were more highly methylated in floral tissue compared to the ECs. Single Nucleotide Polymorphism (SNP) analysis showed an elevated ratio of non-synonymous to synonymous SNPs within the ALSGs (1.99) and CBSGs (1.65) relative to the EC set (0.92), mainly caused by an elevated number of non-synonymous SNPs, indicating that they are fast-evolving at the protein sequence level.

Conclusions

Our analyses suggest that while a significant fraction of the A. thaliana proteome is conserved within the Plant Kingdom, evolutionarily distinct sets of genes that may function in defining biological processes unique to these lineages have arisen within the Brassicaceae and A. thaliana.

Independent S-locus mutations caused self-fertility in Arabidopsis thaliana

A common yet poorly understood evolutionary transition among flowering plants is a switch from outbreeding to an inbreeding mode of mating. The model plant Arabidopsis thaliana evolved to an inbreeding state through the loss of self-incompatibility, a pollen-rejection system in which pollen recognition by the stigma is determined by tightly linked and co-evolving alleles of the S-locus receptor kinase (SRK) and its S-locus cysteine-rich ligand (SCR). Transformation of A. thaliana, with a functional AlSRKb-SCRb gene pair from its outcrossing relative A. lyrata, demonstrated that A. thaliana accessions harbor different sets of cryptic self-fertility–promoting mutations, not only in S-locus genes, but also in other loci required for self-incompatibility. However, it is still not known how many times and in what manner the switch to self-fertility occurred in the A. thaliana lineage. Here, we report on our identification of four accessions that are reverted to full self-incompatibility by transformation with AlSRKb-SCRb, bringing to five the number of accessions in which self-fertility is due to, and was likely caused by, S-locus inactivation. Analysis of S-haplotype organization reveals that inter-haplotypic recombination events, rearrangements, and deletions have restructured the S locus and its genes in these accessions. We also perform a Quantitative Trait Loci (QTL) analysis to identify modifier loci associated with self-fertility in the Col-0 reference accession, which cannot be reverted to full self-incompatibility. Our results indicate that the transition to inbreeding occurred by at least two, and possibly more, independent S-locus mutations, and identify a novel unstable modifier locus that contributes to self-fertility in Col-0.

The mating system adopted by a species has a profound influence on extent of polymorphism, population structure, and evolutionary potential. In flowering plants, the switch from outbreeding to inbreeding has occurred repeatedly, yet little is known about the underlying genetic events. This is true even for the model species A. thaliana, a highly self-fertile member of the crucifer family. In this family, outbreeding is enforced by a self-incompatibility system controlled by the S locus, which involves the recognition of pollen by the stigma to prevent self-fertilization and familial inbreeding. We recently demonstrated that A. thaliana accessions may be reverted to full or partial self-incompatibility by transformation with S-locus genes isolated from its close self-incompatible relative A. lyrata. Despite much recent debate, however, we still do not know how A. thaliana became self-fertile. Here, we use our recently established A. thaliana transgenic self-incompatible experimental model to address these issues. Analysis of the S locus in accessions that can be reverted to full self-incompatibility demonstrates that self-fertility in A. thaliana arose by at least two independent S-locus mutations. Furthermore, analysis of an accession that expresses only partial self-incompatibility shows that self-fertility is associated with an unstable allele at a locus unlinked to the S locus.

Molecular population genetics of the SRK and SCR self-incompatibility genes in the wild plant species Brassica cretica (Brassicaceae)

Self-incompatibility (SI) in plants is a classic example of a trait evolving under strong frequency-dependent selection. As a consequence, population genetic theory predicts that the S locus, which controls SI, should maintain numerous alleles, display a high level of nucleotide diversity, and, in structured populations, show a lower level of among-population differentiation compared to neutral loci. Population-level investigations of DNA sequence variation at the S locus have recently been carried out in the genus Arabidopsis, largely confirming results from theoretical models of S-locus evolutionary dynamics, but no comparable studies have been done in wild Brassica species. In this study, we sequenced parts of the S-locus genes SRK and SCR, two tightly linked genes that are directly involved in the determination of SI specificity in samples from four natural populations of the wild species Brassica cretica. The amount and distribution of nucleotide diversity, as well as the frequency spectrum of putative functional haplotypes, observed at the S locus in B. cretica fit very well with expectations from theoretical models, providing strong evidence for frequency-dependent selection acting on the S locus in a wild Brassica species.

The pattern of amplification and differentiation of Ty1-copia and Ty3-gypsy retrotransposons in Brassicaceae species

One of the causes of genome size expansion is considered to be amplification of retrotransposons. We determined nucleotide sequences of 24 PCR products for each of six retrotransposons in Brassica rapa and Brassica oleracea. Phylogenetic trees of these sequences showed species-specific clades. We also sequenced STF7a homologs and Tto1 homologs, 24 PCR products each, in nine diploids and three allopolyploids, and constructed phylogenetic trees. In these phylogenetic trees, species-specific clades of diploid species were also formed, but retrotransposons of allopolyploids were clustered into the clades of their original genomes, indicating that these two retrotransposons amplified after speciation of the nine diploids. Genetic variation in these retrotransposons may have arisen before emergence of allopolyploid species. There was a positive correlation between the genome size and the average number of substitutions of STF7a and Tto1 homologs in at least seven diploids. The implications of these results in the genome evolution of Brassicaceae are herein discussed.

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

Hypomethylation and transcriptional reactivation of retrotransposon-like sequences in ddm1 transgenic plants of Brassica rapa

DNA methylation and histone modification play important roles in regulating gene expression. The DDM1 gene in Arabidopsis thaliana (AtDDM1) is required for the maintenance of DNA methylation level and histone H3 methylation pattern. We isolated DDM1 homologs of Brassica rapa, BrDDM1a and BrDDM1b, which have 84.4% and 84.1% deduced amino acid sequence identities with AtDDM1, respectively. Both the BrDDM1a and BrDDM1b genes were found to be expressed in vegetative and reproductive tissues. B. rapa ddm1-RNAi transgenic plants with reduced levels of BrDDM1a/BrDDM1b expression showed genome-wide and non-tissue-specific demethylation. These results suggest that the BrDDM1a and BrDDM1b genes are orthologs of AtDDM1 and are required for the maintenance of DNA methylation as is AtDDM1. Despite genome-wide demethylation, developmental abnormalities were not found in the ddm1-RNAi transgenic plants. Dominance relationships of SP11/SCR alleles, the determinant of pollen recognition specificity in Brassica self-incompatibility, in S heterozygotes in B. rapa were not influenced by the low level of the BrDDM1 expression. Transcriptional reactivation of retrotransposon-like sequences observed in the ddm1-RNAi transgenic plants indicates that BrDDM1a and BrDDM1b participate in silencing of retrotransposons. Hypomethylation states of the ddm1-RNAi transgenic plants were inherited by plants of the next generation even by plants which had lost the RNAi construct by segregation. Remethylation was observed in a few progenies. Efficiencies of remethylation in the progenies without the RNAi construct were different between 18S rDNA, BoSTF12a/15a, and BrTto1 sequences.

How far are we from unravelling self-incompatibility in grasses?

The genetic and physiological mechanisms involved in limiting self-fertilization in angiosperms, referred to as self-incompatibility (SI), have significant effects on population structure and have potential diversification and evolutionary consequences. Up to now, details of the underlying genetic control and physiological basis of SI have been elucidated in two different gametophytic SI (GSI) systems, the S-RNase SI and the Papaver SI systems, and the sporophytic SI (SSI) system (Brassica). In the grass family (Poaceae), which contains all the cereal and major forage crops, SI has been known for half a century to be controlled gametophytically by two multiallelic and independent loci, S and Z. But still none of the gene products for S and Z is known and only limited information on related biochemical responses is available. Here we compare current knowledge of grass SI with that of other well-characterized SI systems and speculate about the relationship between SSI and grass SI. Additionally, we discuss comparative mapping as a tool for the further investigation of grass SI.

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.

Self-compatibility in Brassica napus is caused by independent mutations in S-locus genes

Brassica napus is an amphidiploid species with the A genome from Brassica rapa and the C genome from Brassica oleracea. Although B. rapa, B. oleracea and artificially synthesized amphidiploids with the AC genome are self-incompatible, B. napus is self-compatible. Six S genotypes were identified in B. napus, five of which had class I S haplotypes from one species and a class II S haplotype from the other species, and mutations causing self-compatibility were identified in three of these S genotypes. The most predominant S genotype (BnS-1;BnS-6), which is that of cv. 'Westar', had a class I S haplotype similar to B. rapa S-47 (BrS-47) and a class II S haplotype similar to B. oleracea S-15 (BoS-15). The stigmas of 'Westar' rejected the pollen grains of both BrS-47 and BoS-15, while reciprocal crossings were compatible. Insertion of a DNA fragment of about 3.6 kb was found in the promoter region of the SP11/SCR allele of BnS-1, and transcripts of SP11/SCR were not detected in 'Westar'. The nucleotide sequence of the SP11 genomic DNA of BnS-6 was 100% identical to that of SP11 of BoS-15. Class I SP11 alleles from one species showed dominance over class II SP11 alleles from the other species in artificially synthesized B. napus lines, suggesting that the non-functional dominant SP11 allele suppressed the expression of the recessive SP11 allele in 'Westar'. Two other S genotypes in B. napus also had non-functional class I S haplotypes together with recessive BnS-6. These observations suggest independent origins of self-compatibility in B. napus.

S locus F-box brothers: multiple and pollen-specific F-box genes with S haplotype-specific polymorphisms in apple and Japanese pear

Although recent findings suggest that the F-box genes SFB/SLF control pollen-part S specificity in the S-RNase-based gametophytic self-incompatibility (GSI) system, how these genes operate in the system is unknown, and functional variation of pollen S genes in different species has been reported. Here, we analyzed the S locus of two species of Maloideae: apple (Malus domestica) and Japanese pear (Pyrus pyrifolia). The sequencing of a 317-kb region of the apple S9 haplotype revealed two similar F-box genes. Homologous sequences were isolated from different haplotypes of apple and Japanese pear, and they were found to be polymorphic genes derived from the S locus. Since each S haplotype contains two or three related genes, the genes were named SFBB for S locus F-box brothers. The SFBB genes are specifically expressed in pollen, and variable regions of the SFBB genes are under positive selection. In a style-specific mutant S haplotype of Japanese pear, the SFBB genes are retained. Apart from their multiplicity, SFBB genes meet the expected characteristics of pollen S. The unique multiplicity of SFBB genes as the pollen S candidate is discussed in the context of mechanistic variation in the S-RNase-based GSI system.

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.