Comparative physical mapping of segments of the genome of Brassica oleracea var. alboglabra that are homoeologous to sequenced regions of chromosomes 4 and 5 of Arabidopsis thaliana

Due to their relatedness to Arabidopsis thaliana (Arabidopsis), the cultivated Brassica species represent the first group of crops with which to evaluate comparative genomics approaches to understanding biological processes and manipulating traits. We have constructed a high-quality binary BAC library (JBo) from genomic DNA of Brassica oleracea var. alboglabra, in order to underpin such investigations. Using the Arabidopsis genome sequence and clones from the JBo library, we have analysed aspects of gene conservation and microsynteny between a 222 kb region of the genome of Arabidopsis and homoeologous segments of the genome of B. oleracea. All 19 predicted genes tested were found to hybridize to clones in the JBo library, indicating a high level of gene conservation. Further analyses and physical mapping with the BAC clones identified allowed us to construct clone contig maps and analyse in detail the gene content and organization in the set of paralogous segments identified in the genome of B. oleracea. Extensive divergence of gene content was observed, both between the B. oleracea paralogous segments and between them and their homoeologous segment within the genome of Arabidopsis. However, the genes present show highly conserved collinearity with their orthologues in the genome of Arabidopsis. We have identified one example of a Brassica gene in a non-collinear position and one rearrangement. Some of the genes not present in the discernible homoeologous regions appear to be located elsewhere in the B. oleracea genome. The implications of our findings for comparative map-based cloning of genes from crop species are discussed.     

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.     

Conservation of the microstructure of genome segments in Brassica napus and its diploid relatives

The cultivated Brassica species are the group of crops most closely related to Arabidopsis thaliana (Arabidopsis). They represent models for the application in crops of genomic information gained in Arabidopsis and provide an opportunity for the investigation of polyploid genome formation and evolution. The scientific literature contains contradictory evidence for the dynamics of the evolution of polyploid genomes. We aimed at overcoming the inherent complexity of Brassica genomes and clarify the effects of polyploidy on the evolution of genome microstructure in specific segments of the genome. To do this, we have constructed bacterial artificial chromosome (BAC) libraries from genomic DNA of B. rapa subspecies trilocularis (JBr) and B. napus var Tapidor (JBnB) to supplement an existing BAC library from B. oleracea. These allowed us to analyse both recent polyploidization (under 10,000 years in B. napus) and more ancient polyploidization events (ca. 20 Myr for B. rapa and B. oleracea relative to Arabidopsis), with an analysis of the events occurring on an intermediate time scale (over the ca. 4 Myr since the divergence of the B. rapa and B. oleracea lineages). Using the Arabidopsis genome sequence and clones from the JBr library, we have analysed aspects of gene conservation and microsynteny between six regions of the genome of B. rapa with the homoeologous regions of the genomes of B. oleracea and Arabidopsis. Extensive divergence of gene content was observed between the B. rapa paralogous segments and their homoeologous segments within the genome of Arabidopsis. A pattern of interspersed gene loss was identified that is similar, but not identical, to that observed in B. oleracea. The conserved genes show highly conserved collinearity with their orthologues across genomes, but a small number of species-specific rearrangements were identified. Thus the evolution of genome microstructure is an ongoing process. Brassica napus is a recently formed polyploid resulting from the hybridization of B. rapa (containing the Brassica A genome) and B. oleracea (containing the Brassica C genome). Using clones from the JBnB library, we have analysed the microstructure of the corresponding segments of the B. napus genome. The results show that there has been little or no change to the microstructure of the analysed segments of the Brassica A and C genomes as a consequence of the hybridization event forming natural B. napus. The observations indicate that, upon polyploid formation, these segments of the genome did not undergo a burst of evolution discernible at the scale of microstructure.     

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.     

Three homologous genes encoding sn-glycerol-3-phosphate acyltransferase 4 exhibit different expression patterns and functional divergence in Brassica napus

Brassica napus is an allotetraploid (AACC) formed from the fusion of two diploid progenitors, Brassica rapa (AA) and Brassica oleracea (CC). Polyploidy and genome-wide rearrangement during the evolution process have resulted in genes that are present as multiple homologs in the B. napus genome. In this study, three B. napus homologous genes encoding endoplasmic reticulum-bound sn-glycerol-3-phosphate acyltransferase 4 (GPAT4) were identified and characterized. Although the three GPAT4 homologs share a high sequence similarity, they exhibit different expression patterns and altered epigenetic features. Heterologous expression in yeast further revealed that the three BnGPAT4 homologs encoded functional GPAT enzymes but with different levels of polypeptide accumulation. Complementation of the Arabidopsis (Arabidopsis thaliana) gpat4 gpat8 double mutant line with individual BnGPAT4 homologs suggested their physiological roles in cuticle formation. Analysis of gpat4 RNA interference lines of B. napus revealed that the BnGPAT4 deficiency resulted in reduced cutin content and altered stomatal structures in leaves. Our results revealed that the BnGPAT4 homologs have evolved into functionally divergent forms and play important roles in cutin synthesis and stomatal development.     

Variation and selection at the CAULIFLOWER floral homeotic gene accompanying the evolution of domesticated Brassica oleracea

The evolution of plant morphologies during domestication events provides clues to the origin of crop species and the evolutionary genetics of structural diversification. The CAULIFLOWER gene, a floral regulatory locus, has been implicated in the cauliflower phenotype in both Arabidopsis thaliana and Brassica oleracea. Molecular population genetic analysis indicates that alleles carrying a nonsense mutation in exon 5 of the B. oleracea CAULIFLOWER (BoCAL) gene are segregating in both wild and domesticated B. oleracea subspecies. Alleles carrying this nonsense mutation are nearly fixed in B. oleracea ssp. botrytis (domestic cauliflower) and B. oleracea ssp. italica (broccoli), both of which show evolutionary modifications of inflorescence structures. Tests for selection indicate that the pattern of variation at this locus is consistent with positive selection at BoCAL in these two subspecies. This nonsense polymorphism, however, is also present in both B. oleracea ssp. acephala (kale) and B. oleracea ssp. oleracea (wild cabbage). These results indicate that specific alleles of BoCAL were selected by early farmers during the domestication of modified inflorescence structures in B. oleracea.     

Comparative analysis of a transposon-rich Brassica oleracea BAC clone with its corresponding sequence in A. thaliana

We compared the sequence of a 96.7 Kb-long BAC clone (B19N3) from Brassica oleracea (broccoli) with its corresponding regions in Arabidopsis thaliana. B19N3 contains eight genes and 15 transposable elements (TEs). The first two genes in this clone, Bo1 and Bo2, have its corresponding region at the end of chromosome V of Arabidopsis (24 Mb). The third gene, Bo3, corresponds to an ortholog at the opposite end (2.6 Mb) of the same chromosome. The other five genes, Bo4 to Bo8 also have a corresponding region on the same chromosome but at 7.7 Mb . These five genes are colinear with those found in the corresponding region of Arabidopsis, which contains, however, 15 genes. Therefore, a cluster of 10 genes is missing in B. oleracea clone (B19N3). All five genes in common have the same order and orientation in the genomes of both species. Their 36 exons constituting the eight homologous genes have high conservation in size and sequence identity in both species. Among these, there is a major gene involved in aliphatic glucosinolate biosynthesis, BoGSL-ELONG (Bo4). Similar to A. thaliana, this gene, has a tandem duplicate, Bo5. A contig for this region was constructed by primer walking and BAC-end-sequencing, revealing general gene colinearity between both species. During the 20 million years separating A. thaliana from B. oleracea from a common ancestor both genomes have diverged by chromosomal rearrangements and differential TE activity. These events, in addition to changes in chromosome number are responsible for the evolution of the genomes of both species. In spite of these changes, both species conserve general colinearity for their corresponding genes.     

Molecular and functional characterization of broccoli EMBRYONIC FLOWER 2 genes

Polycomb group (PcG) proteins regulate major developmental processes in Arabidopsis. EMBRYONIC FLOWER 2 (EMF2), the VEFS domain-containing PcG gene, regulates diverse genetic pathways and is required for vegetative development and plant survival. Despite widespread EMF2-like sequences in plants, little is known about their function other than in Arabidopsis and rice. To study the role of EMF2 in broccoli (Brassica oleracea var. italica cv. Elegance) development, we identified two broccoli EMF2 (BoEMF2) genes with sequence homology to and a similar gene expression pattern to that in Arabidopsis (AtEMF2). Reducing their expression in broccoli resulted in aberrant phenotypes and gene expression patterns. BoEMF2 regulates genes involved in diverse developmental and stress programs similar to AtEMF2 in Arabidopsis. However, BoEMF2 differs from AtEMF2 in the regulation of flower organ identity, cell proliferation and elongation, and death-related genes, which may explain the distinct phenotypes. The expression of BoEMF2.1 in the Arabidopsis emf2 mutant (Rescued emf2) partially rescued the mutant phenotype and restored the gene expression pattern to that of the wild type. Many EMF2-mediated molecular and developmental functions are conserved in broccoli and Arabidopsis. Furthermore, the restored gene expression pattern in Rescued emf2 provides insights into the molecular basis of PcG-mediated growth and development.     

Genome evolution in Arabidopsis/Brassica: conservation and divergence of ancient rearranged segments and their breakpoints

Since the tetraploidization of the Arabidopsis thaliana ancestor 30-35 million years ago (Mya), a wave of chromosomal rearrangements have modified its genome architecture. The dynamics of this process is unknown, as it has so far been impossible to date individual rearrangement events. In this paper, we present evidence demonstrating that the majority of rearrangements occurred before the Arabidopsis-Brassica split 20-24 Mya, and that the segmental architecture of the A. thaliana genome is predominantly conserved in Brassica. This finding is based on the conservation of four rearrangement breakpoints analysed by fluorescence in situ hybridization (FISH) and RFLP mapping of three A. thaliana chromosomal regions. For this purpose, 95 Arabidopsis bacterial artificial chromosomes (BACs) spanning a total of 8.25 Mb and 81 genetic loci for 36 marker genes were studied in the Brassica oleracea genome. All the regions under study were triplicated in the B. oleracea genome, confirming the hypothesis of Brassica ancestral genome triplication. However, whilst one of the breakpoints was conserved at one locus, it was not at the two others. Further comparison of their organization may indicate that the evolution of the hexaploid Brassica progenitor proceeded by several events, separated in time. Genetic mapping and reprobing with rDNA allowed assignment of the regions to particular Brassica chromosomes. Based on this study of regional organization and evolution, a new insight into polyploidization/diploidization cycles is proposed.     

Contrasting genome organisation: two regions of the Brassica oleracea genome compared with collinear regions of the Arabidopsis thaliana genome.

Brassica crop species are of worldwide importance and are closely related to the model plant Arabidopsis thaliana for which the complete genome sequence has recently been established. We investigated collinearity of marker order by comparing two contrasting regions of the Brassica oleracea genome with homologous regions of A. thaliana. Although there is widespread replication of marker loci in both A. thaliana and B. oleracea, we found that a combination of genetic markers mapped in B. oleracea, including RFLPs, CAPS, and SSRs allowed comparison and interpretation of medium-scale chromosomal organisation and rearrangements. The interpretation of data was facilitated by hybridising probes onto the whole A. thaliana genome, as represented by BAC contigs. Twenty marker loci were sampled from the whole length of the shortest B. oleracea linkage group, 06, and 21 from a 30.4-cM section of the longest linkage group, 03. There is evidence of locus duplication on linkage group 06. Locus order is well conserved between a putative duplicated region of 10.5 cM and a discrete region comprising 25 cM of A. thaliana chromosome I. This was supported by evidence from seven paralogous loci, three of which were duplicated in a 30.6-cM region of linkage group 06. The pattern of locus order for the remainder of linkage group 06 and the sampled section of linkage group 03 was more complex when compared with the A. thaliana genome. Although there was some conservation of locus order between markers on linkage group 03 and approximately 9 cM of A. thaliana chromosome I, this was superimposed upon a complex pattern of additional loci that were replicated in both A. thaliana and B. oleracea. The results are discussed in the context of the ability to use collinear information to assist map-based cloning.     

Multiple flowering time QTLs within several Brassica species could be the result of duplicated copies of one ancestral gene.

Quantitative trait locus (QTL) analysis was used to study the evolution of genes controlling the timing of flowering in four Brassica genomes that are all extensively replicated. Comparative mapping showed that a chromosomal region from the top of Arabidopsis thaliana chromosome 5 corresponded to three homoeologous copies in each of the diploid species Brassica nigra, B. oleracea, and B. rapa and six copies in the amphidiploid B. juncea. QTLs were detected in two of the three replicated segments in each diploid genome and in three of the six replicated segments in B. juncea. These results indicate that, for the studied trait, multiple QTLs resulting from genome duplication is the rule rather than the exception. Brassica homologues to two candidate genes (CO and FLC) identified from the corresponding A. thaliana region were mapped. CO homologues mapped close to the QTL peaks in eight of nine QTLs, while FLC homologues mapped farther away in those cases where the mapping resolution allowed a comparison. Thus, our data are consistent with the hypothesis that all the major QTLs we detected in the different species of Brassica could be the result of duplicated copies of the same ancestral gene, possibly the ancestor of CO.     

Expression, mapping, and genetic variability of Brassica napus disease resistance gene analogues.

Numerous sequences analogous to resistance (R) genes exist in plant genomes and could be involved in resistance traits. The aim of this study was to identify a large number of Brassica napus sequences related to R genes and also to test the adequacy of specific PCR-based tools for studying them. Different consensus primers were compared for their efficiency in amplifying resistance-gene analogues (RGAs) related to the nucleotide-binding-site subgroup of R genes. Specific primers were subsequently designed to fine-study the different RGAs and we tested their efficiency in three species related to B. napus: Brassica oleracea, Brassica rapa, and Arabidopsis thaliana. Forty-four B. napus RGAs were identified. Among 29 examined, at least one-third were expressed. Eighteen RGAs were mapped on 10 of the 19 B. napus linkage groups. The high variability within these sequences permitted discrimination of each genotype within a B. napus collection. The RGA-specific primers amplified RGAs in the B. oleracea and B. rapa genomes, but the sequences appear to be poorly conserved in A. thaliana. Specific RGA primers are a precise tool for studying known-sequence RGAs. These sequences represent interesting markers that could be correlated with resistance traits in B. napus or related Brassica genomes.     

Chromosomal mapping of <Emphasis Type="Italic"> Brassica oleracea </Emphasis>based on ESTs from <Emphasis Type="Italic"> Arabidopsis thaliana</Emphasis>: complexity of the comparative map

Expressed sequence tags (ESTs) from the Arabidopsis thaliana sequencing project were used to construct a genetic RFLP map for Brassica oleracea. Of the 110 A. thaliana ESTs tested, 95 were found to be informative RFLP probes in map construction. In total, 212 new loci corresponding to the 95 ESTs were added to the existing genetic map of B. oleracea. The enriched map covers all nine basic linkage groups and confirms that the chromosomes of B. oleracea and A. thaliana are similar in linear organization. However, varying levels of sequence conservation between the chromosomes of B. oleracea and A. thaliana were detected in different regions of the genomes. Long conserved regions encompassing entire chromosome arms in both genomes were identified; these are probably shared by descent. On the other hand, extensive rearrangements were observed in numerous chromosome regions, producing a mosaic of A. thaliana-like segments in the genome of Brassica. The presence of extensive chromosome duplication in A. thaliana was taken into consideration in the construction of the comparative maps of B. oleracea and A. thaliana.