MADS-box transcription factors (TFs) are important in floral organ specification as well as several other aspects of plant growth and development. Studies on stress resistance-related functions of MADS-box genes are very limited and no such functional studies in Brassica rapa have been reported. To gain insight into this gene family and to elucidate their roles in organ development and stress resistance, we performed genome-wide identification, characterization and expression analysis of MADS-box genes in B. rapa. Whole-genome survey of B. rapa revealed 167 MADS-box genes, which were categorized into type I (Mα, Mβ and Mγ) and type II (MIKCc and MIKC*) based on phylogeny, protein motif structure and exon-intron organization. Expression analysis of 89 MIKCc and 11 MIKC* genes was then carried out. In addition to those with floral and vegetative tissue expression, we identified MADS-box genes with constitutive expression patterns at different stages of flower development. More importantly, from a low temperature-treated whole-genome microarray data set, 19 BrMADS genes were found to show variable transcript abundance in two contrasting inbred lines of B. rapa. Among these, 13 BrMADS genes were further validated and their differential expression was monitored in response to cold stress in the same two lines via qPCR expression analysis. Additionally, the set of 19 BrMADS genes was analyzed under drought and salt stress, and 8 and 6 genes were found to be induced by drought and salt, respectively. The extensive annotation and transcriptome profiling reported in this study will be useful for understanding the involvement of MADS-box genes in stress resistance in addition to their growth and developmental functions, which ultimately provides the basis for functional characterization and exploitation of the candidate genes for genetic engineering of B. rapa.

Genome wide transcription analysis in response to stresses is essential to provide the basis of effective engineering strategies to improve stress tolerance in crop plants. In order to perform transcriptome analysis in Brassica rapa, we constructed a B. rapa oligo microarray, KBGP-24K, using sequence information from approximately 24,000 unigenes and analyzed cold (4 degrees C), salt (250 mM NaCl), and drought (air-dry) treated B. rapa plants. Among the B. rapa unigenes represented on the microarray, 417 (1.7%), 202 (0.8%), and 738 (3.1%) were identified as responsive genes that were differently expressed 5-fold or more at least once during a 48-h treatment with cold, salt, and drought, respectively. These results were confirmed by RT-PCR analysis. In the abiotic stress responsive genes identified, we found 56 transcription factor genes and 60 commonly responsive genes. It suggests that various transcriptional regulatory mechanisms and common signaling pathway are working together under the abiotic stresses in B. rapa. In conclusion, our new developed 24K oligo microarray will be a useful tool for transcriptome profiling and this work will provide valuable insight in the response to abiotic stress in B. rapa.

Simple Sequence Repeats (SSRs) represent short tandem duplications found within all eukaryotic organisms. To examine the distribution of SSRs in the genome of Brassica rapa ssp. pekinensis, SSRs from different genomic regions representing 17.7 Mb of genomic sequence were surveyed. SSRs appear more abundant in non-coding regions (86.6%) than in coding regions (13.4%). Comparison of SSR densities in different genomic regions demonstrated that SSR density was greatest within the 5'-flanking regions of the predicted genes. The proportion of different repeat motifs varied between genomic regions, with trinucleotide SSRs more prevalent in predicted coding regions, reflecting the codon structure in these regions. SSRs were also preferentially associated with gene-rich regions, with peri-centromeric heterochromatin SSRs mostly associated with retrotransposons. These results indicate that the distribution of SSRs in the genome is non-random. Comparison of SSR abundance between B. rapa and the closely related species Arabidopsis thaliana suggests a greater abundance of SSRs in B. rapa, which may be due to the proposed genome triplication. Our results provide a comprehensive view of SSR genomic distribution and evolution in Brassica for comparison with the sequenced genomes of A. thaliana and Oryza sativa.