Retention of triplicated phytoene synthase (PSY) genes in Brassica napus L. and its diploid progenitors during the evolution of the Brassiceae

The extent of genome redundancy exhibited by Brassica species provides a model to study the evolutionary fate of multi-copy genes and the effects of polyploidy in economically important crops. Phytoene synthase (PSY) catalyzes the first committed reaction of the carotenoid biosynthetic pathway, which has been shown to be rate-limiting in Brassica napus seeds. In Arabidopsis thaliana, a single PSY gene (AtPSY) regulates phytoene synthesis in all tissues. Considering that diploid Brassica genomes contain three Arabidopsis-like subgenomes, the objectives of the present work were to determine whether PSY gene families exist in B. napus (AACC) and its diploid progenitor species, Brassica rapa (AA) and Brassica oleracea (CC); to establish the level of retention of Brassica PSY genes; to map PSY gene family members in the A and C genomes and to compare Brassica PSY gene expression patterns. A total of 12 PSY homologues were identified, 6 in B. napus (BnaX.PSY.a-f) and 3 in B. rapa (BraA.PSY.a-c) and B. oleracea (BolC.PSY.a-c). Indeed, with six members, B. napus has the largest PSY gene family described to date. Sequence comparison between AtPSY and Brassica PSY genes revealed a highly conserved gene structure and identity percentages above 85% at the coding sequence (CDS) level. Altogether, our data indicate that PSY gene family expansion preceded the speciation of B. rapa and B. oleracea, dating back to the paralogous subgenome triplication event. In these three Brassica species, all PSY homologues are expressed, exhibiting overlapping redundancy and signs of subfunctionalization among photosynthetic and non-photosynthetic tissues. This evidence supports the hypothesis that functional divergence of PSY gene expression facilitates the accumulation of high levels of carotenoids in chromoplast-rich tissues. Thus, functional retention of triplicated Brassica PSY genes could be at least partially explained by the selective advantage provided by increased levels of gene product in floral organs. A better understanding of carotenogenesis in Brassica will aid in the future development of transgenic and conventional cultivars with carotenoid-enriched oil.     

Expression profiles and genomic organisation of group A protein phosphatase 2C genes in Brassica oleracea

We investigated the expression profiles and genomic organisation of the ABA‐responsive genes encoding protein phosphatases 2C (PP2C, group A members) in Brassica oleracea to better understand their functional and genetic relations. Gene expression profiling of drought responsive genes in B. oleracea and Arabidopsis thaliana revealed significant differences in the gene expression pattern of a key regulator of ABA signalling—ABI1 PP2C. This finding prompted us to study genetic relations within the PP2Cs group A in the Brassica species. Twenty homologous B. oleracea sequences were identified and characterised as putative PP2C group A members. Phylogenetic analysis revealed that the B. oleracea homologues were closely related to the particular members of the A. thaliana PP2C. The genetic analysis corroborated the presence of two to three gene copies in B. oleracea in comparison to the nine unique PP2C genes in the A. thaliana genome. Gene expression analyses showed significant differences in PP2C gene expression pattern in B. oleracea. Our results indicate that PP2C‐based drought stress signalling in B. oleracea has evolved distinctly. Different reactions of particular B. oleracea PP2C genes to drought stress and ABA treatment indicate low conservation of gene expression patterns and functional divergence between B. oleracea and A. thaliana homologous genes.     

Identification of anther-specific genes in a cruciferous model plant, Arabidopsis thaliana, by using a combination of Arabidopsis macroarray and mRNA derived from Brassica oleracea

Arabidopsis thaliana, a member of the Brassicaceae, is a model plant whose genome was the first higher plant genome to be sequenced. Because of the small size of the flowers, it is difficult to dissect and separate reproductive organs (anthers and pistils) at different developmental stages in A. thaliana. In order to perform genome-wide identification of anther-specific genes in A. thaliana, an Arabidopsis cDNA macroarray was hybridized to cDNA derived from anthers and pistils of another crucifer, Brassica oleracea. After scanning the signal intensity for each clone, and cluster analysis, 52 anther-specific genes were identified. These clones contained several anther-specific genes that have already been characterized, as well as novel anther-specific genes. In RT-PCR analysis with mRNA of A. thaliana and B. oleracea, the expression pattern of one-third of the clones was similar to that determined by cDNA macroarray. This system of heterologous hybridization analysis (Arabidopsis cDNA macroarray vs Brassica tissue-specific mRNA) should be applicable to other model species and their close relatives.     

Expression analysis of an Arabidopsis C2H2 zinc finger protein gene

C₂H₂ zinc finger protein genes encode nucleic acid-binding proteins involved in the regulation of gene activity. AtZFP1 (Arabidopsis thaliana zinc finger protein 1) is one member of a small family of C₂H₂ zinc finger-encoding sequences previously characterized from Arabidopsis. The genomic sequence corresponding to the AtZFP1 cDNA has been determined. Molecular analysis demonstrates that AtZFP1 is a unique, intronless gene which encodes a 1100 nucleotides mRNA highly expressed in roots and stems. A construct in which 2.5 kb of AtZFP1 upstream sequences is linked to the -glucuronidase gene was introduced into Arabidopsis by Agrobacterium-mediated transformation of roots. Histochemical analysis of transgenic Arabidopsis carrying the AtZFP1 promotor:-glucuronidase fusion shows good correlation with RNA blot hybridization analysis. This transgenic line will be a useful tool for analyzing the regulation of AtZFP1 to further our understanding of its function.     

Promiscuous and specific phospholipid binding by domains in ZAC,a membrane-associated Arabidopsis protein with an ARF GAP zinc finger and a C2 domain

Arabidopsis proteins were predicted which share an 80 residue zinc finger domain known from ADP-ribosylation factor GTPase-activating proteins (ARF GAPs). One of these is a 37 kDa protein, designated ZAC, which has a novel domain structure in which the N-terminal ARF GAP domain and a C-terminal C2 domain are separated by a region without homology to other known proteins. Zac promoter/-glucuronidase reporter assays revealed highest expression levels in flowering tissue, rosettes and roots. ZAC protein was immuno-detected mainly in association with membranes and fractionated with Golgi and plasma membrane marker proteins. ZAC membrane association was confirmed in assays by a fusion between ZAC and the green fluorescence protein and prompted an analysis of the in vitro phospholipid-binding ability of ZAC. Phospholipid dot-blot and liposome-binding assays indicated that fusion proteins containing the ZAC-C2 domain bind anionic phospholipids non-specifically, with some variance in Ca²⁺ and salt dependence. Similar assays demonstrated specific affinity of the ZAC N-terminal region (residues 1–174) for phosphatidylinositol 3-monophosphate (PI-3-P). Binding was dependent in part on an intact zinc finger motif, but proteins containing only the zinc finger domain (residues 1–105) did not bind PI-3-P. Recombinant ZAC possessed GTPase-activating activity on Arabidopsis ARF proteins. These data identify a novel PI-3-P-binding protein region and thereby provide evidence that this phosphoinositide is recognized as a signal in plants. A role for ZAC in the regulation of ARF-mediated vesicular transport in plants is discussed.     

Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3′ end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.     

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.     

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.     

A review of sources of resistance to turnip yellows virus (TuYV) in Brassica species

Turnip yellows virus (TuYV; previously known as beet western yellows virus) causes major diseases of Brassica species worldwide resulting in severe yield-losses in arable and vegetable crops. It has also been shown to reduce the quality of vegetables, particularly cabbage where it causes tip burn. Incidences of 100% have been recorded in commercial crops of winter oilseed rape (Brassica napus) and vegetable crops (particularly Brassica oleracea) in Europe. This review summarises the known sources of resistance to TuYV in B. napus (AACC genome), Brassica rapa (AA genome) and B. oleracea (CC genome). It also proposes names for the quantitative trait loci (QTLs) responsible for the resistances, Tu rnip Y ellows virus R esistance (TuYR), that have been mapped to at least the chromosome level in the different Brassica species. There is currently only one known source of resistance deployed commercially (TuYR1). This resistance is said to have originated in B. rapa and was introgressed into the A genome of oilseed rape via hybridisation with B. oleracea to produce allotetraploid (AACC) plants that were then backcrossed into oilseed rape. It has been utilised in the majority of known TuYV-resistant oilseed rape varieties. This has placed significant selection pressure for resistance-breaking mutations arising in TuYV. Further QTLs for resistance to TuYV (TuYR2-TuYR9) have been mapped in the genomes of B. napus, B. rapa and B. oleracea and are described here. QTLs from the latter two species have been introgressed into allotetraploid plants, providing for the first time, combined resistance from both the A and the C genomes for deployment in oilseed rape. Introgression of these new resistances into commercial oilseed rape and vegetable brassicas can be accelerated using the molecular markers that have been developed. The deployment of these resistances should lessen selection pressure for resistance-breaking isolates of TuYV and thereby prolong the effectiveness of each other and extant resistance.