The genome of the mesopolyploid crop species Brassica rapaThe Brassica rapa Genome Sequencing Project Consortium 1 Abstract:The Brassicaceae family which includes Arabidopsis thaliana, is a natural priority for reaching beyond botanical models to more deeply sample angiosperm genomic and functional diversity. Here we report the draft genome sequence and its annoation of Brassica rapa, one of the two ancestral species of oilseed rape. We modeled 41,174 protein-coding genes in the B. rapa genome. B. rapa has experienced only the second genome triplication reported to date, with its close relationship to A. thaliana providing a useful outgroup for investigating many consequences of triplication for its structural and functional evolution. The extent of gene loss (fractionation) among triplicated genome segments varies, with one copy containing a greater proportion of genes expected to have been present in its ancestor (70%) than the remaining two (46% and 36%). Both a generally rapid evolutionary rate, and specific copy number amplifications of particular gene families, may contribute to the remarkable propensity of Brassica species for the development of new morphological variants. The B. rapa genome provides a new resource for comparative and evolutionary analysis of the Brassicaceae genomes and also a platform for genetic improvement of Brassica oil and vegetable crops.2

No abstract

Genic microsatellite markers, also known as functional markers, are preferred over anonymous markers as they reveal the variation in transcribed genes among individuals. In this study, we developed a total of 707 expressed sequence tag-derived simple sequence repeat markers (EST-SSRs) and used for development of a high-density integrated map using four individual mapping populations of B. rapa. This map contains a total of 1426 markers, consisting of 306 EST-SSRs, 153 intron polymorphic markers, 395 bacterial artificial chromosome-derived SSRs (BAC-SSRs), and 572 public SSRs and other markers covering a total distance of 1245.9 cM of the B. rapa genome. Analysis of allelic diversity in 24 B. rapa germplasm using 234 mapped EST-SSR markers showed amplification of 2 alleles by majority of EST-SSRs, although amplification of alleles ranging from 2 to 8 was found. Transferability analysis of 167 EST-SSRs in 35 species belonging to cultivated and wild brassica relatives showed 42.51% (Sysimprium leteum) to 100% (B. carinata, B. juncea, and B. napus) amplification. Our newly developed EST-SSRs and high-density linkage map based on highly transferable genic markers would facilitate the molecular mapping of quantitative trait loci and the positional cloning of specific genes, in addition to marker-assisted selection and comparative genomic studies of B. rapa with other related species.

Fine mapping of six seed glucosinolate QTL (J2Gsl1, J3Gsl2, J9Gsl3, J16Gsl4, J17Gsl5 and J3Gsl6) (Ramchiary et al. in Theor Appl Genet 116:77-85, 2007a) was undertaken by the candidate gene approach. Based on the DNA sequences from Arabidopsis and Brassica oleracea for the different genes involved in the aliphatic glucosinolate biosynthesis, candidate genes were amplified and sequenced from high to low glucosinolate Brassica juncea lines Varuna and Heera, respectively. Of the 20 paralogues identified, 17 paralogues belonging to six gene families were mapped to 12 of the 18 linkage groups of B. juncea genome. Co-mapping of candidate genes with glucosinolate QTL revealed that the candidate gene BjuA.GSL-ELONG.a mapped to the QTL interval of J2Gsl1, BjuA.GSL-ELONG.c, BjuA.GSL-ELONG.d and BjuA.Myb28.a mapped to the QTL interval of J3Gsl2, BjuA.GSL-ALK.a mapped to the QTL interval of J3Gsl6 and BjuB.Myb28.a mapped to the QTL interval of J17Gsl5. The QTL J9Gsl3 and J16Gsl4 did not correspond to any of the mapped candidate genes. The functionality and contribution of different candidate genes/QTL was assessed by allelic variation study using phenotypic data of 785 BC(4)DH lines. It was observed that BjuA.Myb28.a and J9Gsl3 contributed significantly to the base level glucosinolate production while J16Gsl4, probably GSL-PRO, BjuA.GSL-ELONG.a and BjuA.GSL-ELONG.c contributed to the C3, C4 and C5 elongation pathways, respectively. Three A genome QTL: J2Gsl1harbouring BjuA.GSL-ELONG.a, J3Gsl2 harbouring both BjuA.GSL-ELONG.c and BjuA.Myb28.a and J9Gsl3, possibly the 'Bronowski genes', were identified as most important loci for breeding low glucosinolate B. juncea. We observed two-step genetic control of seed glucosinolate in B. juncea mainly effected by these three A genome QTL. This study, therefore, provides clues to the genetic mechanism of 'Bronowski genes' controlling the glucosinolate trait and also provides efficient markers for marker-assisted introgression of low glucosinolate trait in B. juncea.

Quantitative trait loci (QTL) analysis of yield influencing traits was carried out in Brassica juncea (AABB) using a doubled haploid (DH) mapping population of 123 lines derived from a cross between Varuna (a line representing the Indian gene pool) and Heera (representing the east European gene pool) to identify potentially useful alleles from both the parents. The existing AFLP based map of B. juncea was further saturated with RFLP and SSR markers which led to the identification of the linkage groups belonging to the A (B. rapa) and B (B. nigra) genome components of B. juncea. For QTL dissection, the DH lines were evaluated at three different environments and phenotyped for 12 quantitative traits. A total of 65 QTL spread over 13 linkage groups (LG) were identified from the three environments. QTL analysis showed that the A genome has contributed more than the B genome to productivity (68% of the total QTL detected) suggesting a more prominent role of the A genome towards domestication of this crop. The east European line, Heera, carried favorable alleles for 42% of the detected QTL and the remaining 58% were in the Indian gene pool line, Varuna. We observed clustering of major QTL in a few linkage groups, particularly in J7 and J10 of the A genome, with QTL of different traits having agronomically antagonistic allelic effects co-mapping to the same genetic interval. QTL analysis also identified some well-separated QTL which could be readily transferred between the two pools. Based on the QTL analysis, we propose that improvement in yield could be achieved more readily by heterosis breeding rather than by pure line breeding.

Foxtail millet (Setaria italica), the second largest cultivated millet crop after pearl millet, is utilized for food and forage globally. Further, it is also considered as a model crop for studying agronomic, nutritional and biofuel traits. In the present study, a genome-wide association study (GWAS) was performed for ten important agronomic traits in 142 foxtail millet core eco-geographically diverse genotypes using 10 K SNPs developed through GBS-ddRAD approach. Number of SNPs on individual chromosome ranged from 844 (chromosome 5) to 2153 (chromosome 8) with an average SNP frequency of 25.9 per Mb. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations was found to decay rapidly with the genetic distance of 177 Kb. However, for individual chromosome, LD decay distance ranged from 76 Kb (chromosome 6) to 357 Kb (chromosome 4). GWAS identified 81 MTAs (marker-trait associations) for ten traits across the genome. High confidence MTAs for three important agronomic traits including FLW (flag leaf width), GY (grain yield) and TGW (thousand-grain weight) were identified. Significant pyramiding effect of identified MTAs further supplemented its importance in breeding programs. Desirable alleles and superior genotypes identified in the present study may prove valuable for foxtail millet improvement through marker-assisted selection.

Seed glucosinolate content in Brassica juncea is a complex quantitative trait. A recurrent selection backcross (RSB) method with a doubled haploid (DH) generation interspersing backcross generations was used for the introgression of low glucosinolate alleles from an east European gene pool B. juncea line, Heera into an Indian gene pool variety, Varuna. Phenotypic comparisons among the DH populations derived from early to advanced backcrosses revealed a shift in the mean values for various glucosinolates with the advancement of backcrossing, indicating a change in the selective values of the alleles with change in the genetic background due to the existence of epistasis and context dependencies. QTL mapping for various seed glucosinolates from early (F(1)DH) and advanced generation (BC(4)DH) populations confirmed the presence of epistasis and context dependency. The common QTL detected in both F(1)DH and BC(4)DH changed their R (2) values from the former to the later generation. Some of the QTL detected in the F(1)DH became irrelevant in the BC(4)DH population. Further, new QTL were detected in the BC(4)DH population for various glucosinolates. A validation study on a population of low glucosinolate DH lines derived from all the backcross generations of the RSB breeding programme revealed that the QTL detected in BC(4)DH were the 'true' QTL. Using glucosinolate as an example, the study provides strong evidence for the importance of the RSB method for the identification of the 'true' QTL which would be significant for marker assisted introgression of a complex quantitative trait whose expression is influenced by epistatic interactions.

BackgroundThe 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.ResultsWe 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.ConclusionsWe 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.