Integration of linkage maps for the Amphidiploid Brassica napus and comparative mapping with Arabidopsis and Brassica rapa

A 1,681-locus consensus genetic map of cultivated cucumber including 67 NB-LRR resistance gene homolog and ten gene loci

BACKGROUND

Cucumber is an important vegetable crop that is susceptible to many pathogens, but no disease resistance (R) genes have been cloned. The availability of whole genome sequences provides an excellent opportunity for systematic identification and characterization of the nucleotide binding and leucine-rich repeat (NB-LRR) type R gene homolog (RGH) sequences in the genome. Cucumber has a very narrow genetic base making it difficult to construct high-density genetic maps. Development of a consensus map by synthesizing information from multiple segregating populations is a method of choice to increase marker density. As such, the objectives of the present study were to identify and characterize NB-LRR type RGHs, and to develop a high-density, integrated cucumber genetic-physical map anchored with RGH loci.

RESULTS

From the Gy14 draft genome, 70 NB-containing RGHs were identified and characterized. Most RGHs were in clusters with uneven distribution across seven chromosomes. In silico analysis indicated that all 70 RGHs had EST support for gene expression. Phylogenetic analysis classified 58 RGHs into two clades: CNL and TNL. Comparative analysis revealed high-degree sequence homology and synteny in chromosomal locations of these RGH members between the cucumber and melon genomes. Fifty-four molecular markers were developed to delimit 67 of the 70 RGHs, which were integrated into a genetic map through linkage analysis. A 1,681-locus cucumber consensus map including 10 gene loci and spanning 730.0 cM in seven linkage groups was developed by integrating three component maps with a bin-mapping strategy. Physically, 308 scaffolds with 193.2 Mbp total DNA sequences were anchored onto this consensus map that covered 52.6% of the 367 Mbp cucumber genome.

CONCLUSIONS

Cucumber contains relatively few NB-LRR RGHs that are clustered and unevenly distributed in the genome. All RGHs seem to be transcribed and shared significant sequence homology and synteny with the melon genome suggesting conservation of these RGHs in the Cucumis lineage. The 1,681-locus consensus genetic-physical map developed and the RGHs identified and characterized herein are valuable genomics resources that may have many applications such as quantitative trait loci identification, map-based gene cloning, association mapping, marker-assisted selection, as well as assembly of a more complete cucumber genome.

Evidence from genome-wide simple sequence repeat markers for a polyphyletic origin and secondary centers of genetic diversity of Brassica juncea in China and India

The oilseed Brassica juncea is an important crop with a long history of cultivation in India and China. Previous studies have suggested a polyphyletic origin of B. juncea and more than one migration from the primary to secondary centers of diversity. We investigated molecular genetic diversity based on 99 simple sequence repeat markers in 119 oilseed B. juncea varieties from China, India, Europe, and Australia to test whether molecular differentiation follows Vavilov's proposal of secondary centers of diversity in India and China. Two distinct groups were identified by markers in the A genome, and the same two groups were confirmed by markers in the B genome. Group 1 included accessions from central and western India, in addition to those from eastern China. Group 2 included accessions from central and western China, as well as those from northern and eastern India. European and Australian accessions were found only in Group 2. Chinese accessions had higher allelic diversity per accession (Group 1) and more private alleles per accession (Groups 1 and 2) than those from India. The marker data and geographic distribution of Groups 1 and 2 were consistent with two independent migrations of B. juncea from its center of origin in the Middle East and neighboring regions along trade routes to western China and northern India, followed by regional adaptation. Group 1 migrated further south and west in India, and further east in China, than Group 2. Group 2 showed diverse agroecological adaptation, with yellow-seeded spring-sown types in central and western China and brown-seeded autumn-sown types in India.

High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L

Background

High density genetic maps built with SNP markers that are polymorphic in various genetic backgrounds are very useful for studying the genetics of agronomical traits as well as genome organization and evolution. Simultaneous dense SNP genotyping of segregating populations and variety collections was applied to oilseed rape (Brassica napus L.) to obtain a high density genetic map for this species and to study the linkage disequilibrium pattern.

Results

We developed an integrated genetic map for oilseed rape by high throughput SNP genotyping of four segregating doubled haploid populations. A very high level of collinearity was observed between the four individual maps and a large number of markers (>59%) was common to more than two maps. The precise integrated map comprises 5764 SNP and 1603 PCR markers. With a total genetic length of 2250 cM, the integrated map contains a density of 3.27 markers (2.56 SNP) per cM. Genotyping of these mapped SNP markers in oilseed rape collections allowed polymorphism level and linkage disequilibrium (LD) to be studied across the different collections (winter vs spring, different seed quality types) and along the linkage groups. Overall, polymorphism level was higher and LD decayed faster in spring than in “00” winter oilseed rape types but this was shown to vary greatly along the linkage groups.

Conclusions

Our study provides a valuable resource for further genetic studies using linkage or association mapping, for marker assisted breeding and for Brassica napus sequence assembly and genome organization analyses.

KNOX1 genes regulate lignin deposition and composition in monocots and dicots

Plant secondary cell walls are deposited mostly in vascular tissues such as xylem vessels, tracheids, and fibers. These cell walls are composed of a complex matrix of compounds including cellulose, hemicellulose, and lignin. Lignin functions primarily to maintain the structural and mechanical integrity of both the transport vessel and the entire plant itself. Since lignin has been identified as a major source of biomass for biofuels, regulation of secondary cell wall biosynthesis has been a topic of much recent investigation. Biosynthesis and patterning of lignin involves many developmental and environmental cues including evolutionarily conserved transcriptional regulatory modules and hormonal signals. Here, we investigate the role of the class I Knotted1-like-homeobox (KNOX) genes and gibberellic acid in the lignin biosynthetic pathway in a representative monocot and a representative eudicot. Knotted1 overexpressing mutant plants showed a reduction in lignin content in both maize and tobacco. Expression of four key lignin biosynthesis genes was analyzed and revealed that KNOX1 genes regulate at least two steps in the lignin biosynthesis pathway. The negative regulation of lignin both in a monocot and a eudicot by the maize Kn1 gene suggests that lignin biosynthesis may be preserved across large phylogenetic distances. The evolutionary implications of regulation of lignification across divergent species are discussed.

Accessing complex crop genomes with next-generation sequencing

Many important crop species have genomes originating from ancestral or recent polyploidisation events. Multiple homoeologous gene copies, chromosomal rearrangements and amplification of repetitive DNA within large and complex crop genomes can considerably complicate genome analysis and gene discovery by conventional, forward genetics approaches. On the other hand, ongoing technological advances in molecular genetics and genomics today offer unprecedented opportunities to analyse and access even more recalcitrant genomes. In this review, we describe next-generation sequencing and data analysis techniques that vastly improve our ability to dissect and mine genomes for causal genes underlying key traits and allelic variation of interest to breeders. We focus primarily on wheat and oilseed rape, two leading examples of major polyploid crop genomes whose size or complexity present different, significant challenges. In both cases, the latest DNA sequencing technologies, applied using quite different approaches, have enabled considerable progress towards unravelling the respective genomes. Our ability to discover the extent and distribution of genetic diversity in crop gene pools, and its relationship to yield and quality-related traits, is swiftly gathering momentum as DNA sequencing and the bioinformatic tools to deal with growing quantities of genomic data continue to develop. In the coming decade, genomic and transcriptomic sequencing, discovery and high-throughput screening of single nucleotide polymorphisms, presence-absence variations and other structural chromosomal variants in diverse germplasm collections will give detailed insight into the origins, domestication and available trait-relevant variation of polyploid crops, in the process facilitating novel approaches and possibilities for genomics-assisted breeding.

Identification of candidate genes of QTLs for seed weight in Brassica napus through comparative mapping among Arabidopsis and Brassica 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.

Integrated consensus genetic and physical maps of flax (Linum usitatissimum L.)

Three linkage maps of flax (Linum usitatissimum L.) were constructed from populations CDC Bethune/Macbeth, E1747/Viking and SP2047/UGG5-5 containing between 385 and 469 mapped markers each. The first consensus map of flax was constructed incorporating 770 markers based on 371 shared markers including 114 that were shared by all three populations and 257 shared between any two populations. The 15 linkage group map corresponds to the haploid number of chromosomes of this species. The marker order of the consensus map was largely collinear in all three individual maps but a few local inversions and marker rearrangements spanning short intervals were observed. Segregation distortion was present in all linkage groups which contained 1-52 markers displaying non-Mendelian segregation. The total length of the consensus genetic map is 1,551 cM with a mean marker density of 2.0 cM. A total of 670 markers were anchored to 204 of the 416 fingerprinted contigs of the physical map corresponding to ~274 Mb or 74 % of the estimated flax genome size of 370 Mb. This high resolution consensus map will be a resource for comparative genomics, genome organization, evolution studies and anchoring of the whole genome shotgun sequence.

A multi-population consensus genetic map reveals inconsistent marker order among maps likely attributed to structural variations in the apple genome

Genetic maps serve as frameworks for determining the genetic architecture of quantitative traits, assessing structure of a genome, as well as aid in pursuing association mapping and comparative genetic studies. In this study, a dense genetic map was constructed using a high-throughput 1,536 EST-derived SNP GoldenGate genotyping platform and a global consensus map established by combining the new genetic map with four existing reliable genetic maps of apple. The consensus map identified markers with both major and minor conflicts in positioning across all five maps. These major inconsistencies among marker positions were attributed either to structural variations within the apple genome, or among mapping populations, or genotyping technical errors. These also highlighted problems in assembly and anchorage of the reference draft apple genome sequence in regions with known segmental duplications. Markers common across all five apple genetic maps resulted in successful positioning of 2875 markers, consisting of 2033 SNPs and 843 SSRs as well as other specific markers, on the global consensus map. These markers were distributed across all 17 linkage groups, with an average of 169±33 marker per linkage group and with an average distance of 0.70±0.14 cM between markers. The total length of the consensus map was 1991.38 cM with an average length of 117.14±24.43 cM per linkage group. A total of 569 SNPs were mapped onto the genetic map, consisting of 140 recombinant individuals, from our recently developed apple Oligonucleotide pool assays (OPA). The new functional SNPs, along with the dense consensus genetic map, will be useful for high resolution QTL mapping of important traits in apple and for pursuing comparative genetic studies in Rosaceae.

Three-point appraisal of genetic linkage maps

This paper develops a simple diagnostic for the investigation of uncertainty within genetic linkage maps using a Bayesian procedure. The method requires only the genotyping data and the proposed genetic map, and calculates the posterior probability for the possible orders of any set of three markers, accounting for the presence of genotyping error (mistyping) and for missing genotype data. The method uses a Bayesian approach to give insight into conflicts between the order in the proposed map and the genotype scores. The method can also be used to assess the accuracy of a genetic map at different genomic scales and to assess alternative potential marker orders. Simulation and two case studies were used to illustrate the method. In the first case study, the diagnostic revealed conflicts in map ordering for short inter-marker distances that were resolved at a distance of 8-12 cM, except for a set of markers at the end of the linkage group. In the second case study, the ordering did not resolve as distances increase, which could be attributed to regions of the map where many individuals were untyped.

Promoter variation and transcript divergence in Brassicaceae lineages of FLOWERING LOCUS T

Brassica napus (AACC, 2n = 38), an oil crop of world-wide importance, originated from interspecific hybridization of B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18), and has six FLOWERING LOCUS T (FT) paralogues. Two located on the homeologous chromosomes A2 and C2 arose from a lineage distinct from four located on A7 and C6. A set of three conserved blocks A, B and C, which were found to be essential for FT activation by CONSTANS (CO) in Arabidopsis, was identified within the FT upstream region in B. napus and its progenitor diploids. However, on chromosome C2, insertion of a DNA transposable element (TE) and a retro-element in FT upstream blocks A and B contributed to significant structural divergence between the A and C genome orthologues. Phylogenetic analysis of upstream block A indicated the conserved evolutionary relationships of distinct FT genes within Brassicaceae. We conclude that the ancient At-α whole genome duplication contributed to distinct ancestral lineages for this key adaptive gene, which co-exist within the same genus. FT-A2 was found to be transcribed in all leaf samples from different developmental stages in both B. rapa and B. napus, whereas FT-C2 was not transcribed in either B. napus or B. oleracea. Silencing of FT-C2 appeared to result from TE insertion and consequent high levels of cytosine methylation in TE sequences within upstream block A. Interestingly, FT-A7/C6 paralogues were specifically silenced in winter type B. napus but abundantly expressed in spring type cultivars under vernalization-free conditions. Motif prediction indicated the presence of two CO protein binding sites within all Brassica block A and additional sites for FT activation in block C. We propose that the ancestral whole genome duplications have contributed to more complex mechanisms of floral regulation and niche adaptation in Brassica compared to Arabidopsis.

A genetic linkage map of Brassica carinata constructed with a doubled haploid population

Brassica carinata is an important oilseed crop with unique favourable traits that are desirable for other Brassica crops. However, given the limited research into genetic resources in B. carinata, knowledge of the genetic structure of this species is relatively poor. Nine homozygous, genetically distinct accessions of B. carinata were obtained via microspore culture, from which two divergent doubled haploid (DH) lines were used to develop a DH mapping population that consisted of 183 lines. The mapping population showed segregation of multiple traits of interest. A genetic map was constructed with PCR-based markers, and a total of 212 loci, which covered 1,703 cM, were assigned to eight linkage groups in the B genome and nine linkage groups in the C genome, which allowed comparison with genetic maps of other important Brassica species that contain the B/C genome(s). Loci for two Mendelian-inherited traits related to pigmentation (petal and anther tip colour) and one quantitative trait (seed coat colour) were identified using the linkage map. The significance of the mapping population in the context of genetic improvement of Brassica crops is discussed.

A high density consensus genetic map of tetraploid cotton that integrates multiple component maps through molecular marker redundancy check

A consensus genetic map of tetraploid cotton was constructed using six high-density maps and after the integration of a sequence-based marker redundancy check. Public cotton SSR libraries (17,343 markers) were curated for sequence redundancy using 90% as a similarity cutoff. As a result, 20% of the markers (3,410) could be considered as redundant with some other markers. The marker redundancy information had been a crucial part of the map integration process, in which the six most informative interspecific Gossypium hirsutum×G. barbadense genetic maps were used for assembling a high density consensus (HDC) map for tetraploid cotton. With redundant markers being removed, the HDC map could be constructed thanks to the sufficient number of collinear non-redundant markers in common between the component maps. The HDC map consists of 8,254 loci, originating from 6,669 markers, and spans 4,070 cM, with an average of 2 loci per cM. The HDC map presents a high rate of locus duplications, as 1,292 markers among the 6,669 were mapped in more than one locus. Two thirds of the duplications are bridging homoeologous A(T) and D(T) chromosomes constitutive of allopolyploid cotton genome, with an average of 64 duplications per A(T)/D(T) chromosome pair. Sequences of 4,744 mapped markers were used for a mutual blast alignment (BBMH) with the 13 major scaffolds of the recently released Gossypium raimondii genome indicating high level of homology between the diploid D genome and the tetraploid cotton genetic map, with only a few minor possible structural rearrangements. Overall, the HDC map will serve as a valuable resource for trait QTL comparative mapping, map-based cloning of important genes, and better understanding of the genome structure and evolution of tetraploid cotton.

Identification and mapping of a major dominant quantitative trait locus controlling seeds per silique as a single Mendelian factor in Brassica napus L

One putative quantitative trait locus (QTL) for seeds per silique (SS), cqSS.A8, was identified using a double haploid (DH) population in Brassica napus, and near-isogenic lines (NILs; BC(3)F(1)) for cqSS.A8 were developed. However, the flanking markers from cqSS.A8 showed no significant difference using single-marker analysis, even though the frequency distribution of SS in the BC(3)F(1) was bimodal, suggesting that one novel locus existed. In this study, we characterized the effects of this locus in the NILs and used a published linkage map to determine its location. A three-step approach was designed for mapping the locus in the NILs (BC(3)F(2)): (1) determining the individual BC(3)F(2) genotype at the locus using a progeny test; (2) identifying amplified fragment length polymorphism (AFLP) markers linked to the locus using a combination of AFLP and bulked segregant analysis; and (3) determining the location and effects of this locus. QTL analysis in the BC(3)F(2) revealed that this locus explained 85.8 and 55.7 % of phenotypic variance for SS and SL, respectively. Its additive and dominant effects on SS were 6.1 and 5.7, respectively. The locus was validated using a DH population by composite interval mapping and located to linkage group C9 (designated as qSS.C9). Mapping qSS.C9 was undertaken using 230 extremely low-SS plants of a BC(4)F(1) population containing 807 plants. We found that qSS.C9 delimited a 1.005-Mb interval including 218 predicted genes in the reference Brassica rapa (Chiifu-401). These results will greatly facilitate map-based cloning of qSS.C9 and seed yield improvement in rapeseed.

Assigning Brassica microsatellite markers to the nine C-genome chromosomes using Brassica rapa var. trilocularis-B. oleracea var. alboglabra monosomic alien addition lines

Brassica rapa var. trilocularis-B. oleracea var. alboglabra monosomic alien addition lines (MAALs) were used to assign simple sequence repeat (SSR) markers to the nine C-genome chromosomes. A total of 64 SSR markers specific to single C-chromosomes were identified. The number of specific markers for each chromosome varied from two (C3) to ten (C4, C7 and C9), where the designation of the chromosomes was according to Cheng et al. (Genome 38:313-319, 1995). Seventeen additional SSRs, which were duplicated on 2-5 C-chromosomes, were also identified. Using the SSR markers assigned to the previously developed eight MAALs and recently obtained aneuploid plants, a new Brassica rapa-B. oleracea var. alboglabra MAAL carrying the alien chromosome C7 was identified and developed. The application of reported genetically mapped SSR markers on the nine MAALs contributed to the determination of the correspondence between numerical C-genome cytological (Cheng et al. in Genome 38:313-319, 1995) and linkage group designations. This correspondence facilitates the integration of C-genome genetic information that has been generated based on the two designation systems and accordingly increases our knowledge about each chromosome. The present study is a significant contribution to genetic linkage analysis of SSR markers and important agronomic traits in B. oleracea and to the potential use of the MAALs in plant breeding.

Map-based cloning of a recessive genic male sterility locus in Brassica napus L. and development of its functional marker

We previously mapped one male-sterile gene (Bnms3) from an extensively used recessive genic male sterility line (9012AB) in Brassica napus to a 0.14-cM genomic region. In this study, two highly homologous BAC contigs possibly containing the candidate BnMs3 gene were identified using a map-based cloning strategy. A BnMs3-linked SCAR marker (DM1) capable of differentiating the subgenomes between B. rapa and the B. oleracea aided mapping of BnMs3 on the contig derived from the B. napus chromosome C9. One representative BAC clone was sequenced from each of the two contigs and resulted in a larger number of markers according to the sequence difference between the two clones. To isolate BnMs3, these markers were then analyzed in another two BC(1) populations with different genetic backgrounds. This assay allowed for a delimitation of the mutated functional region of BnMs3 to a 9.3-kb DNA fragment. Gene prediction suggested that one complete open reading frame (ORF, ORF2) and partial CDS fragments of ORF1 and ORF3 reside in this fragment. Sequence comparison and genetic transformation eventually indicated that ORF1 (designated as BnaC9.Tic40), an analogue of the Arabidopsis gene AT5G16620 which encodes a translocon of the inner envelope of chloroplasts 40 (Tic40), is the only candidate gene of BnMs3. Furthermore, two distinct mutation types in ORF1 both causing the male-sterile phenotype were individually revealed from 9012A and the temporary maintainer line T45. The molecular mechanism of this male sterility as well as the application of BnMs3-associated functional and cosegregated markers in true breeding programs was also discussed.

Exploiting comparative mapping among Brassica species to accelerate the physical delimitation of a genic male-sterile locus (BnRf) in Brassica napus

The recessive genic male sterility (RGMS) line 9012AB has been used as an important pollination control system for rapeseed hybrid production in China. Here, we report our study on physical mapping of one male-sterile locus (BnRf) in 9012AB by exploiting the comparative genomics among Brassica species. The genetic maps around BnRf from previous reports were integrated and enriched with markers from the Brassica A7 chromosome. Subsequent collinearity analysis of these markers contributed to the identification of a novel ancestral karyotype block F that possibly encompasses BnRf. Fourteen insertion/deletion markers were further developed from this conserved block and genotyped in three large backcross populations, leading to the construction of high-resolution local genetic maps where the BnRf locus was restricted to a less than 0.1-cM region. Moreover, it was observed that the target region in Brassica napus shares a high collinearity relationship with a region from the Brassica rapa A7 chromosome. A BnRf-cosegregated marker (AT3G23870) was then used to screen a B. napus bacterial artificial chromosome (BAC) library. From the resulting 16 positive BAC clones, one (JBnB089D05) was identified to most possibly contain the BnRf (c) allele. With the assistance of the genome sequence from the Brassica rapa homolog, the 13.8-kb DNA fragment covering both closest flanking markers from the BAC clone was isolated. Gene annotation based on the comparison of microcollinear regions among Brassica napus, B. rapa and Arabidopsis showed that five potential open reading frames reside in this fragment. These results provide a foundation for the characterization of the BnRf locus and allow a better understanding of the chromosome evolution around BnRf.

Seed colour loci, homoeology and linkage groups of the C genome chromosomes revealed in Brassica rapa-B. oleracea monosomic alien addition lines

BACKGROUND AND AIMS

Brassica rapa and B. oleracea are the progenitors of oilseed rape B. napus. The addition of each chromosome of B. oleracea to the chromosome complement of B. rapa results in a series of monosomic alien addition lines (MAALs). Analysis of MAALs determines which B. oleracea chromosomes carry genes controlling specific phenotypic traits, such as seed colour. Yellow-seeded oilseed rape is a desirable breeding goal both for food and livestock feed end-uses that relate to oil, protein and fibre contents. The aims of this study included developing a missing MAAL to complement an available series, for studies on seed colour control, chromosome homoeology and assignment of linkage groups to B. oleracea chromosomes.

METHODS

A new batch of B. rapa-B. oleracea aneuploids was produced to generate the missing MAAL. Seed colour and other plant morphological features relevant to differentiation of MAALs were recorded. For chromosome characterization, Snow's carmine, fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) were used.

KEY RESULTS

The final MAAL was developed. Morphological traits that differentiated the MAALs comprised cotyledon number, leaf morphology, flower colour and seed colour. Seed colour was controlled by major genes on two B. oleracea chromosomes and minor genes on five other chromosomes of this species. Homoeologous pairing was largely between chromosomes with similar centromeric positions. FISH, GISH and a parallel microsatellite marker analysis defined the chromosomes in terms of their linkage groups. Conclusions A complete set of MAALs is now available for genetic, genomic, evolutionary and breeding perspectives. Defining chromosomes that carry specific genes, physical localization of DNA markers and access to established genetic linkage maps contribute to the integration of these approaches, manifested in the confirmed correspondence of linkage groups with specific chromosomes. Applications include marker-assisted selection and breeding for yellow seeds.

InterStoreDB: a generic integration resource for genetic and genomic data

Associating phenotypic traits and quantitative trait loci (QTL) to causative regions of the underlying genome is a key goal in agricultural research. InterStoreDB is a suite of integrated databases designed to assist in this process. The individual databases are species independent and generic in design, providing access to curated datasets relating to plant populations, phenotypic traits, genetic maps, marker loci and QTL, with links to functional gene annotation and genomic sequence data. Each component database provides access to associated metadata, including data provenance and parameters used in analyses, thus providing users with information to evaluate the relative worth of any associations identified. The databases include CropStoreDB, for management of population, genetic map, QTL and trait measurement data, SeqStoreDB for sequence-related data and AlignStoreDB, which stores sequence alignment information, and allows navigation between genetic and genomic datasets. Genetic maps are visualized and compared using the CMAP tool, and functional annotation from sequenced genomes is provided via an EnsEMBL-based genome browser. This framework facilitates navigation of the multiple biological domains involved in genetics and genomics research in a transparent manner within a single portal. We demonstrate the value of InterStoreDB as a tool for Brassica research. InterStoreDB is available from: http://www.interstoredb.org.

Non-random distribution of extensive chromosome rearrangements in Brassica napus depends on genome organization

Chromosome rearrangements are common, but their dynamics over time, mechanisms of occurrence and the genomic features that shape their distribution and rate are still poorly understood. We used allohaploid Brassica napus (AC, n = 19) as a model to analyze the effect of genomic features on the formation and diversity of meiotically driven chromosome rearrangements. We showed that allohaploid B. napus meiosis leads to extensive new structural diversity. Almost every allohaploid offspring carried a unique combination of multiple rearrangements throughout the genome, and was thus structurally differentiated from both its haploid parent and its sister plants. This large amount of genome reshuffling was remarkably well-tolerated in the heterozygous state, as neither male nor female fertility were strongly reduced, and meiosis behavior was normal in most cases. We also used a quantitative statistical model, which accounted for 75% of the observed variation in rearrangement rates, to show that the distribution of meiotically driven chromosome rearrangements was not random but was shaped by three principal genomic features. In descending order of importance, the rate of marker loss increased strongly with genetic distance from the centromere, the degree of collinearity between chromosomes, and the genome of origin (A < C). Overall, our results demonstrate that B. napus accumulates a large number of genetic changes, but these rearrangements are not randomly distributed in the genome. The structural genetic diversity produced by the allohaploid pathway and its role in the evolution of polyploid species compared to diploid meiosis are discussed.

Molecular mapping of Arabidopsis thaliana lipid-related orthologous genes in Brassica napus

Quantitative Trait Loci (QTL) for oil content has been previously analyzed in a SG-DH population from a cross between a Chinese cultivar and a European cultivar of Brassica napus. Eight QTL with additive and epistatic effects, and with environmental interactions were evaluated. Here we present an integrated linkage map of this population predominantly based on informative markers derived from Brassica sequences, including 249 orthologous A. thaliana genes, where nearly half (112) are acyl lipid metabolism related genes. Comparative genomic analysis between B. napus and A. thaliana revealed 33 colinearity regions. Each of the conserved A. thaliana segments is present two to six times in the B. napus genome. Approximately half of the mapped lipid-related orthologous gene loci (76/137) were assigned in these conserved colinearity regions. QTL analysis for seed oil content was performed using the new map and phenotypic data from 11 different field trials. Nine significant QTL were identified on linkage groups A1, A5, A7, A9, C2, C3, C6 and C8, together explaining 57.79% of the total phenotypic variation. A total of 14 lipid related candidate gene loci were located in the confidence intervals of six of these QTL, of which ten were assigned in the conserved colinearity regions and felled in the most frequently overlapped QTL intervals. The information obtained from this study demonstrates the potential role of the suggested candidate genes in rapeseed kernel oil accumulation.

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.

Saturation of an intra-gene pool linkage map: towards a unified consensus linkage map for fine mapping and synteny analysis in common bean

Map-based cloning and fine mapping to find genes of interest and marker assisted selection (MAS) requires good genetic maps with reproducible markers. In this study, we saturated the linkage map of the intra-gene pool population of common bean DOR364 × BAT477 (DB) by evaluating 2,706 molecular markers including SSR, SNP, and gene-based markers. On average the polymorphism rate was 7.7% due to the narrow genetic base between the parents. The DB linkage map consisted of 291 markers with a total map length of 1,788 cM. A consensus map was built using the core mapping populations derived from inter-gene pool crosses: DOR364 × G19833 (DG) and BAT93 × JALO EEP558 (BJ). The consensus map consisted of a total of 1,010 markers mapped, with a total map length of 2,041 cM across 11 linkage groups. On average, each linkage group on the consensus map contained 91 markers of which 83% were single copy markers. Finally, a synteny analysis was carried out using our highly saturated consensus maps compared with the soybean pseudo-chromosome assembly. A total of 772 marker sequences were compared with the soybean genome. A total of 44 syntenic blocks were identified. The linkage group Pv6 presented the most diverse pattern of synteny with seven syntenic blocks, and Pv9 showed the most consistent relations with soybean with just two syntenic blocks. Additionally, a co-linear analysis using common bean transcript map information against soybean coding sequences (CDS) revealed the relationship with 787 soybean genes. The common bean consensus map has allowed us to map a larger number of markers, to obtain a more complete coverage of the common bean genome. Our results, combined with synteny relationships provide tools to increase marker density in selected genomic regions to identify closely linked polymorphic markers for indirect selection, fine mapping or for positional cloning.

Physical Mapping in a Triplicated Genome: Mapping the Downy Mildew Resistance Locus Pp523 in Brassica oleracea L

We describe the construction of a BAC contig and identification of a minimal tiling path that encompass the dominant and monogenically inherited downy mildew resistance locus Pp523 of Brassica oleracea L. The selection of BAC clones for construction of the physical map was carried out by screening gridded BAC libraries with DNA overgo probes derived from both genetically mapped DNA markers flanking the locus of interest and BAC-end sequences that align to Arabidopsis thaliana sequences within the previously identified syntenic region. The selected BAC clones consistently mapped to three different genomic regions of B. oleracea. Although 83 BAC clones were accurately mapped within a ∼4.6 cM region surrounding the downy mildew resistance locus Pp523, a subset of 33 BAC clones mapped to another region on chromosome C8 that was ∼60 cM away from the resistance gene, and a subset of 63 BAC clones mapped to chromosome C5. These results reflect the triplication of the Brassica genomes since their divergence from a common ancestor shared with A. thaliana, and they are consonant with recent analyses of the C genome of Brassica napus. The assembly of a minimal tiling path constituted by 13 (BoT01) BAC clones that span the Pp523 locus sets the stage for map-based cloning of this resistance gene.

Syntenic relationships between cucumber (Cucumis sativus L.) and melon (C. melo L.) chromosomes as revealed by comparative genetic mapping

BACKGROUND

Cucumber, Cucumis sativus L. (2n = 2 × = 14) and melon, C. melo L. (2n = 2 × = 24) are two important vegetable species in the genus Cucumis (family Cucurbitaceae). Both species have an Asian origin that diverged approximately nine million years ago. Cucumber is believed to have evolved from melon through chromosome fusion, but the details of this process are largely unknown. In this study, comparative genetic mapping between cucumber and melon was conducted to examine syntenic relationships of their chromosomes.

RESULTS

Using two melon mapping populations, 154 and 127 cucumber SSR markers were added onto previously reported F(2)- and RIL-based genetic maps, respectively. A consensus melon linkage map was developed through map integration, which contained 401 co-dominant markers in 12 linkage groups including 199 markers derived from the cucumber genome. Syntenic relationships between melon and cucumber chromosomes were inferred based on associations between markers on the consensus melon map and cucumber draft genome scaffolds. It was determined that cucumber Chromosome 7 was syntenic to melon Chromosome I. Cucumber Chromosomes 2 and 6 each contained genomic regions that were syntenic with melon chromosomes III+V+XI and III+VIII+XI, respectively. Likewise, cucumber Chromosomes 1, 3, 4, and 5 each was syntenic with genomic regions of two melon chromosomes previously designated as II+XII, IV+VI, VII+VIII, and IX+X, respectively. However, the marker orders in several syntenic blocks on these consensus linkage maps were not co-linear suggesting that more complicated structural changes beyond simple chromosome fusion events have occurred during the evolution of cucumber.

CONCLUSIONS

Comparative mapping conducted herein supported the hypothesis that cucumber chromosomes may be the result of chromosome fusion from a 24-chromosome progenitor species. Except for a possible inversion, cucumber Chromosome 7 has largely remained intact in the past nine million years since its divergence from melon. Meanwhile, many structural changes may have occurred during the evolution of the remaining six cucumber chromosomes. Further characterization of the genomic nature of Cucumis species closely related to cucumber and melon might provide a better understanding of the evolutionary history leading to modern cucumber.