Origins of the amphiploid species <Emphasis Type="Italic">Brassica napus</Emphasis>L. investigated by chloroplast and nuclear molecular markers

Background  

The amphiploid species Brassica napus (oilseed rape, Canola) is a globally important oil crop yielding food, biofuels and industrial compounds such as lubricants and surfactants. Identification of the likely ancestors of each of the two genomes (designated A and C) found in B. napus would facilitate incorporation of novel alleles from the wider Brassica genepool in oilseed rape crop genetic improvement programmes. Knowledge of the closest extant relatives of the genotypes involved in the initial formation of B. napus would also allow further investigation of the genetic factors required for the formation of a stable amphiploid and permit the more efficient creation of fully fertile re-synthesised B. napus. We have used a combination of chloroplast and nuclear genetic markers to investigate the closest extant relatives of the original maternal progenitors of B. napus. This was based on a comprehensive sampling of the relevant genepools, including 83 accessions of A genome B. rapa L. (both wild and cultivated types), 94 accessions of B. napus and 181 accessions of C genome wild and cultivated B. oleracea L. and related species.     

Hybridisation and introgression between Brassica napus and B. rapa in the Netherlands

We used flow cytometry, chromosome counting and AFLP markers to investigate gene flow from the crop plant oilseed rape, Brassica napus (AACC) to wild B. rapa (AA) in the Netherlands. From 89 B. napus source populations investigated, all near cropping fields or at transhipment sites, only 19 contained a B. rapa population within a 2.5‐km radius. During our survey we found only three populations with F1 hybrids (AAC), as recognized by their nine extra chromosomes and by flow cytometry. These hybrids were all collected in mixed populations where the two species grew in close proximity. Populations with F1 hybrids were not close to crops, but instead were located on road verges with highly disturbed soils, in which both species were probably recruited from the soil seed bank. Many plants in the F2, BC1 or higher backcrosses are expected to carry one to eight C chromosomes. However, these plants were not observed among the hybrids. We further investigated introgression with molecular markers (AFLP) and compared sympatric B. rapa populations (near populations of B. napus) with control populations of B. rapa (no B. napus within at least 7 km). We found no difference between sympatric and control populations in the number of C markers in B. rapa, nor did we find that these sympatric populations closely resembled B. napus. Our data show that hybrids occur but also suggest no recent introgression of alleles from the crop plant B. napus into wild B. rapa in the Dutch populations studied.     

Genome-wide analysis and stress-responsive expression of CCCH zinc finger family genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Background

Ubiquitous CCCH nucleic acid-binding motif is found in a wide-variety of organisms. CCCH genes are involved in plant developmental processes and biotic and abiotic stress responses. Brassica rapa is a vital economic crop and classical model plant of polyploidy evolution, but the functions of CCCH genes in B. rapa are unclear.

Results

In this study, 103 CCCH genes in B. rapa were identified. A comparative analysis of the chromosomal position, gene structure, domain organization and duplication event between B. rapa and Arabidopsis thaliana were performed. Results showed that CCCH genes could be divided into 18 subfamilies, and segmental duplication might mainly contribute to this family expansion. C-X_7/8-C-X₅-C₃-H was the most commonly found motif, but some novel CCCH motifs were also found, along with some loses of typical CCCH motifs widespread in other plant species. The multifarious gene structures and domain organizations implicated functional diversity of CCCH genes in B. rapa. Evidence also suggested functional redundancy in at least one subfamily due to high conservation between members. Finally, the expression profiles of subfamily-IX genes indicated that they are likely involved in various stress responses.

Conclusion

This study provides the first genome-wide characterization of the CCCH genes in B. rapa. The results suggest that B. rapa CCCH genes are likely functionally divergent, but mostly involved in plant development and stress response. These results are expected to facilitate future functional characterization of this potential RNA-binding protein family in Brassica crops.

     

Sequence‐level comparative analysis of the Brassica napus genome around two stearoyl‐ACP desaturase loci

We conducted a sequence‐level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl‐ACP desaturase‐encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (～56%) is lower than has been reported previously between A. thaliana and Brassica (～66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 ± 3.9% and 85.8 ± 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 ± 3.1% between the B. napus A genome and B. rapa, and 93.1 ± 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.     

Gene transferability from transgenic Brassica napus L. to various subspecies and varieties of Brassica rapa

Gene transferability from transgenic rapeseed to various subspecies and varieties of Brassica rapa was assessed in this study. Artificial crossability was studied in 118 cultivars of 7 B. rapa subspecies and varieties with the transgenic rapeseed GT73 (Brassica napus) as the pollen donor. On average 5.7 seeds were obtained per pollination, with a range from 0.05 to 19.4. The heading type of B. rapa L. showed significantly higher crossability than non-heading types of B. rapa. The spontaneous outcrossing rate between B. rapa (female) and the transgenic rapeseed Ms8 × Rf3 (B. napus) (male) ranged from 0.039 to 0.406%, with an average of 0.19%. The fertilization process and the development of the hybrid seeds as shown by fluorescent staining techniques indicated that the number of adhered pollens on the stigma was reduced by 80%, the number of pollen tubes in the style was reduced by 2/3 and the fertilization time was delayed by over 20 h when pollinated with the transgenic rapeseed Ms8 × Rf3 in comparison with the bud self-pollination of B. rapa as control. About 10–70% of the interspecific hybrid embryos were aborted in the course of development. Some seeds looked cracked in mature pods, which showed germination abilities lower than 10%. The spontaneous outcrossing rates were much lower than the artificial crossability, and their survival fitness of the interspecific hybrid was very low, indicating that it should be possible to keep the adventitious presence of the off-plants under the allowed threshold, if proper measures are taken.     

Anthocyanin biosynthetic genes in Brassica rapa

Background

Anthocyanins are a group of flavonoid compounds. As a group of important secondary metabolites, they perform several key biological functions in plants. Anthocyanins also play beneficial health roles as potentially protective factors against cancer and heart disease. To elucidate the anthocyanin biosynthetic pathway in Brassica rapa, we conducted comparative genomic analyses between Arabidopsis thaliana and B. rapa on a genome-wide level.

Results

In total, we identified 73 genes in B. rapa as orthologs of 41 anthocyanin biosynthetic genes in A. thaliana. In B. rapa, the anthocyanin biosynthetic genes (ABGs) have expanded and most genes exist in more than one copy. The anthocyanin biosynthetic structural genes have expanded through whole genome and tandem duplication in B. rapa. More structural genes located upstream of the anthocyanin biosynthetic pathway have been retained than downstream. More negative regulatory genes are retained in the anthocyanin biosynthesis regulatory system of B. rapa.

Conclusions

These results will promote an understanding of the genetic mechanism of anthocyanin biosynthesis, as well as help the improvement of the nutritional quality of B. rapa through the breeding of high anthocyanin content varieties.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-426) contains supplementary material, which is available to authorized users.     

Development of public immortal mapping populations,molecular markers and linkage maps for rapid cycling <Emphasis Type="Italic">Brassica rapa</Emphasis> and <Emphasis Type="Italic">B. oleracea</Emphasis>

Publicly available genomic tools help researchers integrate information and make new discoveries. In this paper, we describe the development of immortal mapping populations of rapid cycling, self-compatible lines, molecular markers, and linkage maps for Brassica rapa and B. oleracea and make the data and germplasm available to the Brassica research community. The B. rapa population consists of 160 recombinant inbred (RI) lines derived from the cross of highly inbred lines of rapid cycling and yellow sarson B. rapa. The B. oleracea population consists of 155 double haploid (DH) lines derived from an F1 cross between two DH lines, rapid cycling and broccoli. A total of 120 RFLP probes, 146 SSR markers, and one phenotypic trait (flower color) were used to construct genetic linkage maps for both species. The B. rapa map consists of 224 molecular markers distributed along 10 linkage groups (A1–A10) with a total distance of 1125.3 cM and a marker density of 5.7 cM/marker. The B. oleracea genetic map consists of 279 molecular markers and one phenotypic marker distributed along nine linkage groups (C1–C9) with a total distance of 891.4 cM and a marker density of 3.2 cM/marker. A syntenic analysis with Arabidopsis thaliana identified collinear genomic blocks that are in agreement with previous studies, reinforcing the idea of conserved chromosomal regions across the Brassicaceae.     

Carotenoid biosynthetic genes in Brassica rapa: comparative genomic analysis,phylogenetic analysis,and expression profiling

Background

Carotenoids are isoprenoid compounds synthesized by all photosynthetic organisms. Despite much research on carotenoid biosynthesis in the model plant Arabidopsis thaliana, there is a lack of information on the carotenoid pathway in Brassica rapa. To better understand its carotenoid biosynthetic pathway, we performed a systematic analysis of carotenoid biosynthetic genes at the genome level in B. rapa.

Results

We identified 67 carotenoid biosynthetic genes in B. rapa, which were orthologs of the 47 carotenoid genes in A. thaliana. A high level of synteny was observed for carotenoid biosynthetic genes between A. thaliana and B. rapa. Out of 47 carotenoid biosynthetic genes in A. thaliana, 46 were successfully mapped to the 10 B. rapa chromosomes, and most of the genes retained more than one copy in B. rapa. The gene expansion was caused by the whole-genome triplication (WGT) event experienced by Brassica species. An expression analysis of the carotenoid biosynthetic genes suggested that their expression levels differed in root, stem, leaf, flower, callus, and silique tissues. Additionally, the paralogs of each carotenoid biosynthetic gene, which were generated from the WGT in B. rapa, showed significantly different expression levels among tissues, suggesting differentiated functions for these multi-copy genes in the carotenoid pathway.

Conclusions

This first systematic study of carotenoid biosynthetic genes in B. rapa provides insights into the carotenoid metabolic mechanisms of Brassica crops. In addition, a better understanding of carotenoid biosynthetic genes in B. rapa will contribute to the development of conventional and transgenic B. rapa cultivars with enriched carotenoid levels in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1655-5) contains supplementary material, which is available to authorized users.     

Beyond genomic variation - comparison and functional annotation of three Brassica rapa genomes: a turnip,a rapid cycling and a Chinese cabbage

Background

Brassica rapa is an economically important crop species. During its long breeding history, a large number of morphotypes have been generated, including leafy vegetables such as Chinese cabbage and pakchoi, turnip tuber crops and oil crops.

Results

To investigate the genetic variation underlying this morphological variation, we re-sequenced, assembled and annotated the genomes of two B. rapa subspecies, turnip crops (turnip) and a rapid cycling. We then analysed the two resulting genomes together with the Chinese cabbage Chiifu reference genome to obtain an impression of the B. rapa pan-genome. The number of genes with protein-coding changes between the three genotypes was lower than that among different accessions of Arabidopsis thaliana, which can be explained by the smaller effective population size of B. rapa due to its domestication. Based on orthology to a number of non-brassica species, we estimated the date of divergence among the three B. rapa morphotypes at approximately 250,000 YA, far predating Brassica domestication (5,000-10,000 YA).

Conclusions

By analysing genes unique to turnip we found evidence for copy number differences in peroxidases, pointing to a role for the phenylpropanoid biosynthesis pathway in the generation of morphological variation. The estimated date of divergence among three B. rapa morphotypes implies that prior to domestication there was already considerably divergence among B. rapa genotypes. Our study thus provides two new B. rapa reference genomes, delivers a set of computer tools to analyse the resulting pan-genome and uses these to shed light on genetic drivers behind the rich morphological variation found in B. rapa.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-250) contains supplementary material, which is available to authorized users.     

Characterization of directly transformed weedy <Emphasis Type="Italic">Brassica rapa</Emphasis> and introgressed <Emphasis Type="Italic">B. rapa</Emphasis> with Bt <Emphasis Type="Italic">cry1Ac</Emphasis> and <Emphasis Type="Italic">gfp</Emphasis> genes

Crop to weed transgene flow, which could result in more competitive weed populations, is an agricultural biosafety concern. Crop Brassica napus to weedy Brassica rapa hybridization has been extensively characterized to better understand the transgene flow and its consequences. In this study, weedy accessions of B. rapa were transformed with Bacillus thuringiensis (Bt) cry1Ac- and green fluorescence protein (gfp)-coding transgenes using Agrobacterium to assess ecological performance of the wild biotype relative to introgressed hybrids in which the transgenic parent was the crop. Regenerated transgenic B. rapa events were characterized by progeny analysis, Bt protein enzyme-linked immunosorbent assay (ELISA), Southern blot analysis, and GFP expression assay. GFP expression level and Bt protein concentration were significantly different between independent transgenic B. rapa events. Similar reproductive productivity was observed in comparison between transgenic B. rapa events and B. rapa × B. napus introgressed hybrids in greenhouse and field experiments. In the greenhouse, Bt transgenic plants experienced significantly less herbivory damage from the diamondback moth (Plutella xylostella). No differences were found in the field experiment under ambient, low, herbivore pressure. Directly transformed transgenic B. rapa plants should be a helpful experimental control to better understand crop genetic load in introgressed transgenic weeds.     

Multiple copies of eukaryotic translation initiation factors in Brassica rapa facilitate redundancy,enabling diversification through variation in splicing and broad‐spectrum virus resistance

Recessive strain‐specific resistance to a number of plant viruses in the Potyvirus genus has been found to be based on mutations in the eukaryotic translation initiation factor 4E (eIF4E) and its isoform, eIF(iso)4E. We identified three copies of eIF(iso)4E in a number of Brassica rapa lines. Here we report broad‐spectrum resistance to the potyvirus Turnip mosaic virus (TuMV) due to a natural mechanism based on the mis‐splicing of the eIF(iso)4E allele in some TuMV‐resistant B. rapa var. pekinensis lines. Of the splice variants, the most common results in a stop codon in intron 1 and a much truncated, non‐functional protein. The existence of multiple copies has enabled redundancy in the host plant's translational machinery, resulting in diversification and emergence of the resistance. Deployment of the resistance is complicated by the presence of multiple copies of the gene. Our data suggest that in the B. rapa subspecies trilocularis, TuMV appears to be able to use copies of eIF(iso)4E at two loci. Transformation of different copies of eIF(iso)4E from a resistant B. rapa line into an eIF(iso)4E knockout line of Arabidopsis thaliana proved misleading because it showed that, when expressed ectopically, TuMV could use multiple copies which was not the case in the resistant B. rapa line. The inability of TuMV to access multiple copies of eIF(iso)4E in B. rapa and the broad spectrum of the resistance suggest it may be durable.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Maternal Inheritance of Chloroplasts between <Emphasis Type="Italic">Brassica rapa</Emphasis> and F<Subscript>1</Subscript>-hybrids Demonstrated by cpDNA Markers Specific to Oilseed Rape and <Emphasis Type="Italic">B. rapa</Emphasis>

Controlled reciprocal crosses between B. rapa and F₁-hybrids (B. napus (♀) × B. rapa), giving 20 pair-crossings, were made to reveal possible irregularities in chloroplast inheritance during production of BC₁s. Despite the close relationship of chloroplasts in B. rapa and B. napus, development of PCR-based molecular markers specific to B. rapa chloroplasts and B. napus chloroplasts was successful. Offspring from each cross were investigated and among these, we found no irregular chloroplast inheritance, since their plastid genotypes in all cases were identical to that of their mother. With a certainty of 95% our data indicate that the probability that chloroplasts are being inherited paternally is less than 0.015. In oilseed rape, pollen-mediated transgene-dispersal poses a well-known risk. Our results support development of transplastomic oilseed rape as an approach to reduce transgene dispersal.     

High‐throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus,Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole‐genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C‐genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Identification of candidate genes of QTLs for seed weight in <Emphasis Type="Italic">Brassica napus</Emphasis> through comparative mapping among <Emphasis Type="Italic">Arabidopsis</Emphasis> and <Emphasis Type="Italic">Brassica</Emphasis>species

Background

Map-based cloning of quantitative trait loci (QTLs) in polyploidy crop species remains a challenge due to the complexity of their genome structures. QTLs for seed weight in B. napus have been identified, but information on candidate genes for identified QTLs of this important trait is still rare.

Results

In this study, a whole genome genetic linkage map for B. napus was constructed using simple sequence repeat (SSR) markers that covered a genetic distance of 2,126.4 cM with an average distance of 5.36 cM between markers. A procedure was developed to establish colinearity of SSR loci on B. napus with its two progenitor diploid species B. rapa and B. oleracea through extensive bioinformatics analysis. With the aid of B. rapa and B. oleracea genome sequences, the 421 homologous colinear loci deduced from the SSR loci of B. napus were shown to correspond to 398 homologous loci in Arabidopsis thaliana. Through comparative mapping of Arabidopsis and the three Brassica species, 227 homologous genes for seed size/weight were mapped on the B. napus genetic map, establishing the genetic bases for the important agronomic trait in this amphidiploid species. Furthermore, 12 candidate genes underlying 8 QTLs for seed weight were identified, and a gene-specific marker for BnAP2 was developed through molecular cloning using the seed weight/size gene distribution map in B. napus.

Conclusions

Our study showed that it is feasible to identify candidate genes of QTLs using a SSR-based B. napus genetic map through comparative mapping among Arabidopsis and B. napus and its two progenitor species B. rapa and B. oleracea. Identification of candidate genes for seed weight in amphidiploid B. napus will accelerate the process of isolating the mapped QTLs for this important trait, and this approach may be useful for QTL identification of other traits of agronomic significance.