The role of FRIGIDA and FLOWERING LOCUS C genes in flowering time of Brassica rapa leafy vegetables.

There is a wide variation of flowering time among lines of Brassica rapa L. Most B. rapa leafy (Chinese cabbage etc.) or root (turnip) vegetables require prolonged cold exposure for flowering, known as vernalization. Premature bolting caused by low temperature leads to a reduction in the yield/quality of these B. rapa vegetables. Therefore, high bolting resistance is an important breeding trait, and understanding the molecular mechanism of vernalization is necessary to achieve this goal. In this study, we demonstrated that BrFRIb functions as an activator of BrFLC in B. rapa. We showed a positive correlation between the steady state expression levels of the sum of the BrFLC paralogs and the days to flowering after four weeks of cold treatment, suggesting that this is an indicator of the vernalization requirement. We indicate that BrFLCs are repressed by the accumulation of H3K27me3 and that the spreading of H3K27me3 promotes stable FLC repression. However, there was no clear relationship between the level of H3K27me3 in the BrFLC and the vernalization requirement. We also showed that if there was a high vernalization requirement, the rate of repression of BrFLC1 expression following prolonged cold treatments was lower.

Next-generation sequencing-based bulked segregant analysis for QTL mapping in the heterozygous species Brassica rapa.

KEY MESSAGE: An improved protocol of QTL-seq, an NGS-based method for bulked segregant analysis we previously developed in rice, allowed successful mapping of QTLs of interest in the highly heterozygous genome of B. rapa, demonstrating the power of this elegant method for genetic analyses in heterozygous species of economic importance. Recent advances in next-generation sequencing (NGS) and the various NGS-based methods developed for rapidly identifying candidate genes of interest have accelerated genetic analysis mainly in the model plants rice and Arabidopsis. Brassica rapa includes several economically important crops such as Chinese cabbage, turnip and various leafy vegetables. The application of NGS-based approaches for the analysis of B. rapa has been limited mainly due to its highly heterozygous genome and poor quality of the reference genome sequence currently available for this species. In this study, we have improved QTL-seq, a method for NGS-based bulked segregant analysis we previously developed in rice, extending its applicability for accelerating the genetic analysis and molecular breeding of B. rapa. Addition of new filters to the original QTL-seq pipeline allowed removal of spurious single-nucleotide polymorphisms caused by alignment/sequencing errors and variability between parents, significantly improving accuracy of the analysis. As proof of principle, we successfully applied the new approach to identify candidate genomic regions controlling flowering and trichome formation using segregating F2 progeny obtained from crosses made between cultivars of B. rapa showing contrasting phenotypes for these traits. We strongly believe that the improved QTL-seq method reported here will extend the applicability of NGS-based genetic analysis not only to B. rapa but also to other plant species of economic importance with heterozygous genomes.

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa.

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2 mm; FB2), early large floral buds (2-4 mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Whole-Genome Resequencing of a Worldwide Collection of Rapeseed Accessions Reveals the Genetic Basis of Ecotype Divergence.

Rapeseed (Brassica napus), an important oilseed crop, has adapted to diverse climate zones and latitudes by forming three main ecotype groups, namely winter, semi-winter, and spring types. However, genetic variations underlying the divergence of these ecotypes are largely unknown. Here, we report the global pattern of genetic polymorphisms in rapeseed determined by resequencing a worldwide collection of 991 germplasm accessions. A total of 5.56 and 5.53 million single-nucleotide polymorphisms (SNPs) as well as 1.86 and 1.92 million InDels were identified by mapping reads to the reference genomes of "Darmor-bzh" and "Tapidor," respectively. We generated a map of allelic drift paths that shows splits and mixtures of the main populations, and revealed an asymmetric evolution of the two subgenomes of B. napus by calculating the genetic diversity and linkage disequilibrium parameters. Selective-sweep analysis revealed genetic changes in genes orthologous to those regulating various aspects of plant development and response to stresses. A genome-wide association study identified SNPs in the promoter regions of FLOWERING LOCUS T and FLOWERING LOCUS C orthologs that corresponded to the different rapeseed ecotype groups. Our study provides important insights into the genomic footprints of rapeseed evolution and flowering-time divergence among three ecotype groups, and will facilitate screening of molecular markers for accelerating rapeseed breeding.

Genome-Wide Identification, Evolutionary and Expression Analyses of the GALACTINOL SYNTHASE Gene Family in Rapeseed and Tobacco.

Galactinol synthase (GolS) is a key enzyme in raffinose family oligosaccharide (RFO) biosynthesis. The finding that GolS accumulates in plants exposed to abiotic stresses indicates RFOs function in environmental adaptation. However, the evolutionary relationships and biological functions of GolS family in rapeseed (Brassica napus) and tobacco (Nicotiana tabacum) remain unclear. In this study, we identified 20 BnGolS and 9 NtGolS genes. Subcellular localization predictions showed that most of the proteins are localized to the cytoplasm. Phylogenetic analysis identified a lost event of an ancient GolS copy in the Solanaceae and an ancient duplication event leading to evolution of GolS4/7 in the Brassicaceae. The three-dimensional structures of two GolS proteins were conserved, with an important DxD motif for binding to UDP-galactose (uridine diphosphate-galactose) and inositol. Expression profile analysis indicated that BnGolS and NtGolS genes were expressed in most tissues and highly expressed in one or two specific tissues. Hormone treatments strongly induced the expression of most BnGolS genes and homologous genes in the same subfamilies exhibited divergent-induced expression. Our study provides a comprehensive evolutionary analysis of GolS genes among the Brassicaceae and Solanaceae as well as an insight into the biological function of GolS genes in hormone response in plants.

A GLABRA1 ortholog on LG A9 controls trichome number in the Japanese leafy vegetables Mizuna and Mibuna (Brassica rapa L. subsp. nipposinica L. H. Bailey): evidence from QTL analysis.

Brassica rapa show a wide range of morphological variations. In particular, the leaf morphologies of the Japanese traditional leafy vegetables Mizuna and Mibuna (Brassica rapa L. subsp. nipposinica L. H. Bailey) are distinctly different, even though they are closely related cultivars that are easy to cross. In addition to the differences in the gross morphology of leaves, some cultivars of Mibuna (Kyo-nishiki) have many trichomes on its leaves, whereas Mizuna (Kyo-mizore) does not. To identify the genes responsible for the different number of trichomes, we performed a quantitative trait loci (QTL) analysis of Mizuna and Mibuna. To construct linkage maps for these cultivars, we used RNA-seq data to develop cleaved amplified polymorphic sequence (CAPS) markers. We also performed a restriction site-associated DNA sequencing (RAD-seq) analysis to detect single-nucleotide polymorphisms (SNPs). Two QTL analyses were performed in different years, and both analyses indicated that the largest effect was found on LG A9. Expression analyses showed that a gene homologous to GLABRA1 (GL1), a transcription factor implicated in trichome development in Arabidopsis thaliana, and the sequences 3'-flanking (downstream) of BrGL1, differed considerably between Mizuna (Kyo-mizore) and Mibuna (Kyo-nishiki). These results indicate that BrGL1 on LG A9 is one of the candidate genes responsible for the difference in trichome number between Mizuna and Mibuna. Detecting genes that are responsible for morphological variations allows us to better understand the breeding history of Mizuna and Mibuna.

Natural variation in CIRCADIAN CLOCK ASSOCIATED 1 is associated with flowering time in Brassica rapa.

Flowering time is a very important agronomic trait and the development of molecular markers associated with this trait can facilitate crop breeding. CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), a core oscillator component of circadian rhythms that affect metabolic pathways in plants, has been implicated in flowering time control in species of Brassica. CCA1 gene sequences from three Brassica rapa inbred lines, showing either early flowering or late flowering phenotypes, were analyzed and a high level of sequence variation was identified, especially within the fourth intron. Using this information, three PCR primer sets were designed and tested using various inbred lines of B. rapa. The usage of InDel markers was further validated by evaluation of flowering time and high resolution melting (HRM) analysis. Both methods, PCR and HRM, validated the use of newly developed markers. Additional sequence analyses of Brassica plants with diploid (AA, BB, or CC) and allotetraploid genomes further confirmed a large number of sequence polymorphisms in the CCA1 gene, including insertions/deletions in the fourth intron. Our results demonstrated that sequence variations in CCA1 can be used to develop valuable trait-related molecular markers for Brassica crop breeding.

Metabolic Profiling in Chinese Cabbage (Brassica rapa L. subsp. pekinensis) Cultivars Reveals that Glucosinolate Content Is Correlated with Carotenoid Content.

A total of 38 bioactive compounds, including glucosinolates, carotenoids, tocopherols, sterols, and policosanols, were characterized from nine varieties of Chinese cabbage (Brassica rapa L. subsp. pekinensis) to determine their phytochemical diversity and analyze their abundance relationships. The metabolite profiles were evaluated with principal component analysis (PCA), Pearson correlation analysis, and hierarchical clustering analysis (HCA). PCA and HCA identified two distinct varieties of Chinese cabbage (Cheonsangcheonha and Waldongcheonha) with higher levels of glucosinolates and carotenoids. Pairwise comparisons of the 38 metabolites were calculated using Pearson correlation coefficients. The HCA, which used the correlation coefficients, clustered metabolites that are derived from closely related biochemical pathways. Significant correlations were discovered between chlorophyll and carotenoids. Additionally, aliphatic glucosinolate and carotenoid levels were positively correlated. The Cheonsangcheonha and Waldongcheonha varieties appear to be good candidates for breeding because they have high glucosinolate and carotenoid levels.

Dissecting the complex molecular evolution and expression of polygalacturonase gene family in Brassica rapa ssp. chinensis.

Polygalacturonases (PGs) participate in pectin disassembly of cell wall and belong to one of the largest hydrolase families in plants. In this study, we identified 99 PG genes in Brassica rapa. Comprehensive analysis of phylogeny, gene structures, physico-chemical properties and coding sequence evolution demonstrated that plant PGs should be classified into seven divergent clades and each clade's members had specific sequence and structure characteristics, and/or were under specific selection pressures. Genomic distribution and retention rate analysis implied duplication events and biased retention contributed to PG family's expansion. Promoter divergence analysis using "shared motif method" revealed a significant correlation between regulatory and coding sequence evolution of PGs, and proved Clades A and E were of ancient origin. Quantitative real-time PCR analysis showed that expression patterns of PGs displayed group specificities in B. rapa. Particularly, nearly half of PG family members, especially those of Clades C, D and F, closely relates to reproductive development. Most duplicates showed similar expression profiles, suggesting dosage constraints accounted for preservation after duplication. Promoter-GUS assay further indicated PGs' extensive roles and possible redundancy during reproductive development. This work can provide a scientific classification of plant PGs, dissect the internal relationships between their evolution and expressions, and promote functional researches.

Comprehensive analysis of the flowering genes in Chinese cabbage and examination of evolutionary pattern of CO-like genes in plant kingdom.

In plants, flowering is the most important transition from vegetative to reproductive growth. The flowering patterns of monocots and eudicots are distinctly different, but few studies have described the evolutionary patterns of the flowering genes in them. In this study, we analysed the evolutionary pattern, duplication and expression level of these genes. The main results were as follows: (i) characterization of flowering genes in monocots and eudicots, including the identification of family-specific, orthologous and collinear genes; (ii) full characterization of CONSTANS-like genes in Brassica rapa (BraCOL genes), the key flowering genes; (iii) exploration of the evolution of COL genes in plant kingdom and construction of the evolutionary pattern of COL genes; (iv) comparative analysis of CO and FT genes between Brassicaceae and Grass, which identified several family-specific amino acids, and revealed that CO and FT protein structures were similar in B. rapa and Arabidopsis but different in rice; and (v) expression analysis of photoperiod pathway-related genes in B. rapa under different photoperiod treatments by RT-qPCR. This analysis will provide resources for understanding the flowering mechanisms and evolutionary pattern of COL genes. In addition, this genome-wide comparative study of COL genes may also provide clues for evolution of other flowering genes.

Comparison of a silver nanoparticle-based method and the modified spectrophotometric methods for assessing antioxidant capacity of rapeseed varieties.

The antioxidant capacity of 15 rapeseed varieties was determined by the proposed silver nanoparticle-based (AgNP) method and three modified assays: ferric reducing antioxidant power (FRAP), 2,2'-diphenyl-1-picrylhydrazyl (DPPH) and Folin-Ciocalteu reducing capacity (FC). The average antioxidant capacities of the studied rapeseed cultivars ranged between 5261-9462, 3708-7112, 18864-31245 and 5816-9937 µmol sinapic acid (SA)/100g for AgNP, FRAP, DPPH and FC methods, respectively. There are significant, positive correlations between antioxidant capacities of the studied rapeseed cultivars determined by four analytical methods (r=0.5971-0.9149, p<0.05). The comparable precision for the proposed AgNP method (RSD=1.4-4.4%) and the modified FRAP, DPPH and FC methods (RSD=1.0-4.4%, 0.7-2.1% and 0.8-3.6%, respectively), demonstrate the benefit of the AgNP method in the routine analysis of antioxidant capacity of rapeseed cultivars. The principal component analysis (PCA) and hierarchical cluster analysis (HCA) were used for discrimination the quality of the studied rapeseed varieties based on their antioxidant potential determined by different analytical methods. Three main groups were identified by HCA, while the classification and characterisation of rapeseed varieties within each of these groups were obtained from PCA. The chemometric analyses demonstrated that, rapeseed variety S13 had the highest antioxidant capacity, thus this cultivar should be considered as the richest source of natural antioxidants.

RNA-seq based SNPs in some agronomically important oleiferous lines of Brassica rapa and their use for genome-wide linkage mapping and specific-region fine mapping.

BACKGROUND: Brassica rapa (AA) contains very diverse forms which include oleiferous types and many vegetable types. Genome sequence of B. rapa line Chiifu (ssp. pekinensis), a leafy vegetable type, was published in 2011. Using this knowledge, it is important to develop genomic resources for the oleiferous types of B. rapa. This will allow more involved molecular mapping, in-depth study of molecular mechanisms underlying important agronomic traits and introgression of traits from B. rapa to major oilseed crops - B. juncea (AABB) and B. napus (AACC). The study explores the availability of SNPs in RNA-seq generated contigs of three oleiferous lines of B. rapa - Candle (ssp. oleifera, turnip rape), YSPB-24 and Tetra (ssp. trilocularis, Yellow sarson) and their use in genome-wide linkage mapping and specific-region fine mapping using a RIL population between Chiifu and Tetra. RESULTS: RNA-seq was carried out on the RNA isolated from young inflorescences containing unopened floral buds, floral axis and small leaves, using Illumina paired-end sequencing technology. Sequence assembly was carried out using the Velvet de-novo programme and the assembled contigs were organised against Chiifu gene models, available in the BRAD-CDS database. RNA-seq confirmed the presence of more than 17,000 single-copy gene models described in the BRAD database. The assembled contigs and the BRAD gene models were analyzed for the presence of SSRs and SNPs. While the number of SSRs was limited, more than 0.2 million SNPs were observed between Chiifu and the three oleiferous lines. Assays for SNPs were designed using KASPar technology and tested on a F7-RIL population derived from a Chiifu x Tetra cross. The design of the SNP assays were based on three considerations - the 50 bp flanking region of the SNPs should be strictly similar, the SNP should have a read-depth of ≥7 and no exon/intron junction should be present within the 101 bp target region. Using these criteria, a total of 640 markers (580 for genome-wide mapping and 60 for specific-region mapping) marking as many genes were tested for mapping. Out of 640 markers that were tested, 594 markers could be mapped unambiguously which included 542 markers for genome-wide mapping and 42 markers for fine mapping of the tet-o locus that is involved with the trait tetralocular ovary in the line Tetra. CONCLUSION: A large number of SNPs and PSVs are present in the transcriptome of B. rapa lines for genome-wide linkage mapping and specific-region fine mapping. Criteria used for SNP identification delivered markers, more than 93% of which could be successfully mapped to the F7-RIL population of Chiifu x Tetra cross.

Deciphering the diploid ancestral genome of the Mesohexaploid Brassica rapa.

The genus Brassica includes several important agricultural and horticultural crops. Their current genome structures were shaped by whole-genome triplication followed by extensive diploidization. The availability of several crucifer genome sequences, especially that of Chinese cabbage (Brassica rapa), enables study of the evolution of the mesohexaploid Brassica genomes from their diploid progenitors. We reconstructed three ancestral subgenomes of B. rapa (n = 10) by comparing its whole-genome sequence to ancestral and extant Brassicaceae genomes. All three B. rapa paleogenomes apparently consisted of seven chromosomes, similar to the ancestral translocation Proto-Calepineae Karyotype (tPCK; n = 7), which is the evolutionarily younger variant of the Proto-Calepineae Karyotype (n = 7). Based on comparative analysis of genome sequences or linkage maps of Brassica oleracea, Brassica nigra, radish (Raphanus sativus), and other closely related species, we propose a two-step merging of three tPCK-like genomes to form the hexaploid ancestor of the tribe Brassiceae with 42 chromosomes. Subsequent diversification of the Brassiceae was marked by extensive genome reshuffling and chromosome number reduction mediated by translocation events and followed by loss and/or inactivation of centromeres. Furthermore, via interspecies genome comparison, we refined intervals for seven of the genomic blocks of the Ancestral Crucifer Karyotype (n = 8), thus revising the key reference genome for evolutionary genomics of crucifers.

Field evaluation of in situ remediation of Cd-contaminated soil using four additives, two foliar fertilisers and two varieties of pakchoi.

This study was conducted to determine the optimal planting mode for pakchoi (Brassica rapa chinensis) in Cd-contaminated soil to reduce the accumulation of Cd in the edible parts while maintaining yields. Four additives (red mud (RM), silicon calcium fertiliser (SC), spodium (SP) and calcium magnesium phosphate (CMP)), two foliar fertilisers (Ca and Zn) and two varieties of pakchoi (Aijiaohuang (AJ) and Baixuegongzhu (BX)) were used in this study. The results show that the addition of SC and RM had an effect, but the other additives did not appear to increase the biomasses of AJ and BX. In some cases, the growth responses of AJ and BX to the same treatment were different. Extra additions of Ca or Zn to additive-treated pakchoi did not help the additives stimulate the growth of AJ and BX, except for SC-treated AJ and BX and SP-treated AJ. The SC and CMP additives significantly reduced the available Cd concentration in both the AJ soil and the BX soil; however, they did not significantly decrease the Cd concentration in the aboveground parts of AJ and BX. The RM treatments (for both levels) and some treatments containing SP reduced the available Cd concentration in the soils and reduced the accumulation of Cd in the two pakchoi varieties. Additions of Ca or Zn fertiliser significantly reduced the Cd concentration in the aboveground parts of AJ and BX. However, when Ca or Zn was sprayed on the additive-treated AJ and BX, they did not help the additives reduce the Cd accumulation in the aboveground parts of AJ and BX, except for the additive CMP. This study shows that RM may be an optimal amendment to reduce the accumulation of Cd in the edible part of pakchoi while simultaneously maintaining yields. The utilisation of Ca or Zn as a foliar fertiliser to additive-treated pakchoi showed positive effects only under some conditions.

Effect of modified atmosphere packaging on the quality and bioactive compounds of Chinese cabbage (Brasicca rapa L. ssp. chinensis).

BACKGROUND: The perishability of Brassica chinensis poses a major challenge to distribution and marketing. The aim of this work was to select a suitable modified atmosphere packaging to retain the overall quality and bioactive compounds during storage. RESULTS: Four types of biorientated polypropylene packaging (BOPP)--BOPP03, BOPP04, BOPP05 and BOPP06--with different perforations were evaluated regarding the maintenance of quality parameters (weight loss, leaf yellowing, colour L*, C*, h°), decay, chlorophyll a, chlorophyll b, bioactive compounds (carotenoids, ascorbic acid, total phenolic compounds), antioxidant scavenging activity, overall appearance and odour evaluation, at 10°C at 2, 4, 6, 8 and 10 days. Leaves were packed in BOPP (two 2-mm holes) and packed and unpacked leaves were included for comparison. The modified atmosphere created (2% O2 and 7% CO2) inside the BOPP05 reduced leaf yellowing (higher h°), improved the overall appearance with acceptable odour, moderately maintained chlorophyll a and b, bioactive compounds and antioxidant scavenging activity, and remained marketable for up to 10 days at 10°C. Gas composition within the packages influenced the retention of bioactive compounds and overall quality. CONCLUSION: Application of BOPP05 is a promising method for extending the shelf life of B. chinensis leaves in order to promote its utilisation and commercialisation via urban fresh-produce markets.

[Identification of rapeseed varieties by infrared spectroscopy combined with spectral retrieval].

The infrared spectra of 75 samples of 15 rapeseed varieties were obtained by means of Fourier transform infrared (FT-IR) spectroscopy technique. The results show that the infrared spectra of the samples are very similar. By the spectral professional software Omnic8.0, two spectral databases were constructed. Lib01 includes the average spectra of 9 samples from Yunnan, while Lib02 is constructed from the average spectra of all 15 varieties. The correlation search and the squared differential difference retrieval of the spectra of Yunnan samples were performed with the spectral database Lib01 in full spectrum range, respectively. The obtained matching values reflect that the samples have the closest correlation with each other only within the same species, with the correct rate of 77.8% for the correlation search, and that of 82.2% for the latter. The squared differential difference retrieval for the 15 sample spectra was carried out with Lib02 in full spectral range and in the specified range of 1 700-950 cm(-1), yielding correct rate of 82.7% and 90.6% respectively. These results show that FTIR combined with the spectral retrieval method can identify the different varieties of rapeseed, and the squared differential difference retrieval in the specified spectral range is a more effective and simple way for distinguishing different varieties of rapeseed.

Identification of Chinese cabbage sentrin as a suppressor of Bax-induced cell death in yeast.

Studies into the cell death program termed apoptosis have resulted in new information regarding how cells control and execute their own demise, including insights into the mechanism by which death-preventing factors can inhibit Bax-induced caspase activation. We investigated high temperature stress-induced cell death in Brassica rapa. Using a yeast functional screening from a Brassica rapa cDNA library, the BH5-127 EST clone encoding an apoptotic suppressor peptide was identified. However, a phylogenic tree showed that BH5-127 clusters within a clade containing SUMO-1 (Small Ubiquitin-like Modifier- 1). BH5-127 was confirmed similar to have function to SUMO-1 as Fas suppression. Expression of BH5-127 showed that substantial suppression of cell death survived on SD-galactose-Leu--Ura- medium. The results suggest that BrSE (Brassica rapa Sentrin EST, BH5-127) is one of the important regulatory proteins in programming cell death, especially in the seedling stage of Chinese cabbage.

Altered patterns of fractionation and exon deletions in Brassica rapa support a two-step model of paleohexaploidy.

The genome sequence of the paleohexaploid Brassica rapa shows that fractionation is biased among the three subgenomes and that the least fractionated subgenome has approximately twice as many orthologs as its close (and relatively unduplicated) relative Arabidopsis than had either of the other two subgenomes. One evolutionary scenario is that the two subgenomes with heavy gene losses (I and II) were in the same nucleus for a longer period of time than the third subgenome (III) with the fewest gene losses. This "two-step" hypothesis is essentially the same as that proposed previously for the eudicot paleohexaploidy; however, the more recent nature of the B. rapa paleohexaploidy makes this model more testable. We found that subgenome II suffered recent small deletions within exons more frequently than subgenome I, as would be expected if the genes in subgenome I had already been near maximally fractionated before subgenome III was introduced. We observed that some sequences, before these deletions, were flanked by short direct repeats, a unique signature of intrachromosomal illegitimate recombination. We also found, through simulations, that short--single or two-gene--deletions appear to dominate the fractionation patterns in B. rapa. We conclude that the observed patterns of the triplicated regions in the Brassica genome are best explained by a two-step fractionation model. The triplication and subsequent mode of fractionation could influence the potential to generate morphological diversity--a hallmark of the Brassica genus.

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.).

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines 'Lynx-037DH' and 'Monty-028DH'. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.

Development of full-length cDNAs from Chinese cabbage (Brassica rapa Subsp. pekinensis) and identification of marker genes for defence response.

Arabidopsis belongs to the Brassicaceae family and plays an important role as a model plant for which researchers have developed fine-tuned genome resources. Genome sequencing projects have been initiated for other members of the Brassicaceae family. Among these projects, research on Chinese cabbage (Brassica rapa subsp. pekinensis) started early because of strong interest in this species. Here, we report the development of a library of Chinese cabbage full-length cDNA clones, the RIKEN BRC B. rapa full-length cDNA (BBRAF) resource, to accelerate research on Brassica species. We sequenced 10 000 BBRAF clones and confirmed 5476 independent clones. Most of these cDNAs showed high homology to Arabidopsis genes, but we also obtained more than 200 cDNA clones that lacked any sequence homology to Arabidopsis genes. We also successfully identified several possible candidate marker genes for plant defence responses from our analysis of the expression of the Brassica counterparts of Arabidopsis marker genes in response to salicylic acid and jasmonic acid. We compared gene expression of these markers in several Chinese cabbage cultivars. Our BBRAF cDNA resource will be publicly available from the RIKEN Bioresource Center and will help researchers to transfer Arabidopsis-related knowledge to Brassica crops.