Comparative genomic analysis of sequences sampled from a small region on soybean (Glycine max) molecular linkage group G.

Eight DNA markers spanning an interval of approximately 10 centimorgans (cM) on soybean (Glycine max) molecular linkage group G (MLG-G) were used to identify bacterial artificial chromosome (BAC) clones. Twenty-eight BAC clones in eight distinct contiguous groups (contigs) were isolated from this genome region, along with 59 BAC clones on 17 contigs homoeologous to those on MLG-G. BAC clones in four of the MLG-G contigs were also digested to produce subclones and detailed physical maps. All of the BAC-ends were sequenced, as were the subclones, to estimate proportions in different sequence categories, compare similarities among homoeologs, and explore microsynteny with Arabidopsis. Homoeologous BAC contigs were enriched in repetitive sequences compared with those on MLG-G or the soybean genome as a whole. Fingerprint and cross-hybridization comparisons between MLG-G and homoeologous contigs revealed cases of highly similar physical organization between soybean duplicates, as did DNA sequence comparisons. Twenty-seven out of 78 total sequences on soybean MLG-G showed significant similarity to Arabidopsis. The homologs mapped to six compact genome segments in Arabidopsis, with the longest containing seven homologs spanning two million base pairs. These results extend previous observations of large-scale duplication and selective gene loss in Arabidopsis, suggesting that networks of conserved synteny between Arabidopsis and other angiosperm families can stretch over long physical distances.     

BAC-end sequence-based SNPs and Bin mapping for rapid integration of physical and genetic maps in apple

A genome-wide BAC physical map of the apple, Malus × domestica Borkh., has been recently developed. Here, we report on integrating the physical and genetic maps of the apple using a SNP-based approach in conjunction with bin mapping. Briefly, BAC clones located at ends of BAC contigs were selected, and sequenced at both ends. The BAC end sequences (BESs) were used to identify candidate SNPs. Subsequently, these candidate SNPs were genetically mapped using a bin mapping strategy for the purpose of mapping the physical onto the genetic map. Using this approach, 52 (23%) out of 228 BESs tested were successfully exploited to develop SNPs. These SNPs anchored 51 contigs, spanning ～ 37 Mb in cumulative physical length, onto 14 linkage groups. The reliability of the integration of the physical and genetic maps using this SNP-based strategy is described, and the results confirm the feasibility of this approach to construct an integrated physical and genetic maps for apple.     

A plant-transformation-competent BIBAC/BAC-based map of rice for functional analysis and genetic engineering of its genomic sequence.

Sequencing of the rice genome has provided a platform for functional genomics research of rice and other cereal species. However, multiple approaches are needed to determine the functions of its genes and sequences and to use the genome sequencing results for genetic improvement of cereal crops. Here, we report a plant-transformation-competent, binary bacterial artificial chromosome (BIBAC) and bacterial artificial chromosome (BAC) based map of rice to facilitate these studies. The map was constructed from 20 835 BIBAC and BAC clones, and consisted of 579 overlapping BIBAC/BAC contigs. To facilitate functional analysis of chromosome 8 genomic sequence and cloning of the genes and QTLs mapped to the chromosome, we anchored the chromosomal contigs to the existing rice genetic maps. The chromosomal map consists of 11 contigs, 59 genetic markers, and 36 sequence tagged sites, spanning a total of ca. 38 Mb in physical length. Comparative analysis between the genetic and physical maps of chromosome 8 showed that there are 3 "hot" and 2 "cold" spots of genetic recombination along the chromosomal arms in addition to the "cold spot" in the centromeric region, suggesting that the sequence component contents of a chromosome may affect its local genetic recombination frequencies. Because of its plant transformability, the BIBAC/BAC map could provide a platform for functional analysis of the rice genome sequence and effective use of the sequencing results for gene and QTL cloning and molecular breeding.     

A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring Pinot Noir bacterial artificial chromosome contigs

The construction of a dense genetic map for Vitis vinifera and its anchoring to a BAC-based physical map is described: it includes 994 loci mapped onto 19 linkage groups, corresponding to the basic chromosome number of Vitis. Spanning 1245 cM with an average distance of 1.3 cM between adjacent markers, the map was generated from the segregation of 483 single-nucleotide polymorphism (SNP)-based genetic markers, 132 simple sequence repeats (SSRs), and 379 AFLP markers in a mapping population of 94 F(1) individuals derived from a V. vinifera cross of the cultivars Syrah and Pinot Noir. Of these markers, 623 were anchored to 367 contigs that are included in a physical map produced from the same clone of Pinot Noir and covering 352 Mbp. On the basis of contigs containing two or more genetically mapped markers, region-dependent estimations of physical and recombinational distances are presented. The markers used in this study include 118 SSRs common to an integrated map derived from five segregating populations of V. vinifera. The positions of these SSR markers in the two maps are conserved across all Vitis linkage groups. The addition of SNP-based markers introduces polymorphisms that are easy to database, are useful for evolutionary studies, and significantly increase the density of the map. The map provides the most comprehensive view of the Vitis genome reported to date and will be relevant for future studies on structural and functional genomics and genetic improvement.     

Genetic and physical mapping of new EST-derived SSRs on the A and B genome chromosomes of wheat

The availability of genetic maps and phenotypic data of segregating populations allows to localize and map agronomically important genes, and to identify closely associated molecular markers to be used in marker-assisted selection and positional cloning. The objective of the present work was to develop a durum wheat intervarietal genetic and physical map based on genomic microsatellite or genomic simple sequence repeats (gSSR) markers and expressed sequence tag (EST)-derived microsatellite (EST-SSR) markers. A set of 122 new EST-SSR loci amplified by 100 primer pairs was genetically mapped on the wheat A and B genome chromosomes. The whole map also comprises 149 gSSR markers amplified by 120 primer pairs used as anchor chromosome loci, two morphological markers (Black colour, Bla1, and spike glaucousness, Ws) and two seed storage protein loci (Gli-A2 and Gli-B2). The majority of SSR markers tested (182) was chromosome-specific. Out of 275 loci 241 loci assembled in 25 linkage groups assigned to the chromosomes of the A and B genome and 34 remained unlinked. A higher percentage of markers (54.4%), localized on the B genome chromosomes, in comparison to 45.6% distributed on the A genome. The whole map covered 1,605 cM. The B genome accounted for 852.2 cM of genetic distance; the A genome basic map spanned 753.1 cM with a minimum length of 46.6 cM for chromosome 5A and a maximum of 156.2 cM for chromosome 3A and an average value of 114.5 cM. The primer sets that amplified two or more loci mapped to homoeologous as well as to non-homoeologous sites. Out of 241 genetically mapped loci 213 (88.4%) were physically mapped by using the nulli-tetrasomic, ditelosomic and a stock of 58 deletion lines dividing the A and B genome chromosomes in 94 bins. No discrepancies concerning marker order were observed but the cytogenetic maps revealed in some cases small genetic distance covered large physical regions. Putative function for mapped SSRs were assigned by searching against GenBank nonredundant database using TBLASTX algorithms.     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

A BAC-based physical map of the apple genome

Genome-wide physical mapping is an essential step toward investigating the genetic basis of complex traits as well as pursuing genomics research of virtually all plant and animal species. We have constructed a physical map of the apple genome from a total of 74,281 BAC clones representing approximately 10.5x haploid genome equivalents. The physical map consists of 2702 contigs, and it is estimated to span approximately 927 Mb in physical length. The reliability of contig assembly was evaluated by several methods, including assembling contigs using variable stringencies, assembling contigs using fingerprints from individual libraries, checking consensus maps of contigs, and using DNA markers. Altogether, the results demonstrated that the contigs were properly assembled. The apple genome-wide BAC-based physical map represents the first draft genome sequence not only for any member of the large Rosaceae family, but also for all tree species. This map will play a critical role in advanced genomics research for apple and other tree species, including marker development in targeted chromosome regions, fine-mapping and isolation of genes/QTL, conducting comparative genomics analyses of plant chromosomes, and large-scale genomics sequencing.     

The Physical and Genetic Framework of the Maize B73 Genome

Maize is a major cereal crop and an important model system for basic biological research. Knowledge gained from maize research can also be used to genetically improve its grass relatives such as sorghum, wheat, and rice. The primary objective of the Maize Genome Sequencing Consortium (MGSC) was to generate a reference genome sequence that was integrated with both the physical and genetic maps. Using a previously published integrated genetic and physical map, combined with in-coming maize genomic sequence, new sequence-based genetic markers, and an optical map, we dynamically picked a minimum tiling path (MTP) of 16,910 bacterial artificial chromosome (BAC) and fosmid clones that were used by the MGSC to sequence the maize genome. The final MTP resulted in a significantly improved physical map that reduced the number of contigs from 721 to 435, incorporated a total of 8,315 mapped markers, and ordered and oriented the majority of FPC contigs. The new integrated physical and genetic map covered 2,120 Mb (93%) of the 2,300-Mb genome, of which 405 contigs were anchored to the genetic map, totaling 2,103.4 Mb (99.2% of the 2,120 Mb physical map). More importantly, 336 contigs, comprising 94.0% of the physical map (∼1,993 Mb), were ordered and oriented. Finally we used all available physical, sequence, genetic, and optical data to generate a golden path (AGP) of chromosome-based pseudomolecules, herein referred to as the B73 Reference Genome Sequence version 1 (B73 RefGen_v1).     

Construction and characterization of BAC libraries from major grapevine cultivars

Genome projects were initiated on grapevine (Vitis vinifera L., 2n=38, genome size 475 Mb) through the successful construction of four bacterial artificial chromosome (BAC) libraries from three major cultivars, Cabernet Sauvignon (Cabernet S), Syrah and two different clones of Pinot Noir (Pinot N). Depending on the library, the genome coverage represented 4.5–14.8 genome equivalents with clones having a mean insert size of 93–158 kb. BAC pools suitable for PCR screening were constructed for two of these BAC libraries [Cabernet S and Pinot N clone (cl) 115] and subsequently used to confirm the genome coverage of both libraries by PCR anchoring of 74 genetic markers sampled from the 19 linkage groups. For ten of these markers, two bands on separate BAC pools were differentiated that could correspond either to different alleles or to a duplication of the locus being studied. Finally, a preliminary assessment of the correspondence between genetic and physical distances was made through the anchoring of all the markers mapped along linkage group 1 of the V. vinifera genetic map. A pair of markers, 2.1 cM apart, anchored the same BAC clones, which allowed us to estimate that 1 cM corresponded in this particular region to a maximum length of 130 kb.     

An integrated map of Arabidopsis thaliana for functional analysis of its genome sequence.

The genome of the model plant species Arabidopsis thaliana has recently been sequenced. To accelerate its current genome research, we developed a whole-genome, BAC/BIBAC-based, integrated physical, genetic, and sequence map of the A. thaliana ecotype Columbia. This new map was constructed from the clones of a new plant-transformation-competent BIBAC library and is integrated with the existing sequence map. The clones were restriction fingerprinted by DNA sequencing gel-based electrophoresis, assembled into contigs, and anchored to an existing genetic map. The map consists of 194 BAC/BIBAC contigs, spanning 126 Mb of the 130-Mb Arabidopsis genome. A total of 120 contigs, spanning 114 Mb, were anchored to the chromosomes of Arabidopsis. Accuracy of the integrated map was verified using the existing physical and sequence maps and numerous DNA markers. Integration of the new map with the sequence map has enabled gap closure of the sequence map and will facilitate functional analysis of the genome sequence. The method used here has been demonstrated to be sufficient for whole-genome physical mapping from large-insert random bacterial clones and thus is applicable to rapid development of whole-genome physical maps for other species.     

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.     

Genomic analysis and marker development for the Tsn1 locus in wheat using bin-mapped ESTs and flanking BAC contigs

The wheat Tsn1 gene confers sensitivity to the host-selective toxin Ptr ToxA produced by the tan spot fungus (Pyrenophora tritici-repentis). The long-term goal of this research is to isolate Tsn1 using a positional cloning approach. Here, we evaluated 54 ESTs (expressed sequence tags) physically mapped to deletion bin 5BL 0.75–0.76, which is a gene-rich region containing Tsn1. Twenty-three EST loci were mapped as either PCR-based single-stranded conformational polymorphism or RFLP markers in a low-resolution wheat population. The genetic map corresponding to the 5BL 0.75–0.76 deletion bin spans 18.5 cM and contains 37 markers for a density of 2 markers/cM. The EST-based genetic map will be useful for tagging other genes, establishing colinearity with rice, and anchoring sequence ready BAC contigs of the 5BL 0.75–0.76 deletion bin. High-resolution mapping showed that EST-derived markers together with previously developed AFLP-derived markers delineated Tsn1 to a 0.8 cM interval. Flanking markers were used to screen the Langdon durum BAC library and contigs of 205 and 228 kb flanking Tsn1 were assembled, sequenced, and anchored to the genetic map. Recombination frequency averaged 760 kb/cM across the 228 kb contig, but no recombination was observed across the 205 kb contig resulting in an expected recombination frequency of more than 10 Mb/cM. Therefore, chromosome walking within the Tsn1 region may be difficult. However, the sequenced BACs allowed the identification of one microsatellite in each contig for which markers were developed and shown to be highly suitable for marker-assisted selection of Tsn1.     

Bacterial artificial chromosome-based physical map of Gibberella zeae (Fusarium graminearum).

Fusarium graminearum is the primary causal pathogen of Fusarium head blight of wheat and barley. To accelerate genomic analysis of F. graminearum, we developed a bacterial artificial chromosome (BAC)-based physical map and integrated it with the genome sequence and genetic map. One BAC library, developed in the HindIII restriction enzyme site, consists of 4608 clones with an insert size of approximately 107 kb and covers about 13.5 genome equivalents. The other library, developed in the BamHI restriction enzyme site, consists of 3072 clones with an insert size of approximately 95 kb and covers about 8.0 genome equivalents. We fingerprinted 2688 clones from the HindIII library and 1536 clones from the BamHI library and developed a physical map of F. graminearum consisting of 26 contigs covering 39.2 Mb. Comparison of our map with the F. graminearum genome sequence showed that the size of our physical map is equivalent to the 36.1 Mb of the genome sequence. We used 31 sequence-based genetic markers, randomly spaced throughout the genome, to integrate the physical map with the genetic map. We also end-sequenced 17 BamHI BAC clones and identified 4 clones that spanned gaps in the genome sequence. Our new integrated map is highly reliable and useful for a variety of genomics studies.     

Identification and High-Density Mapping of Gene-Rich Regions in Chromosome Group 5 of Wheat

The distribution of genes and recombination in the wheat genome was studied by comparing physical maps with the genetic linkage maps. The physical maps were generated by mapping 80 DNA and two phenotypic markers on an array of 65 deletion lines for homoeologous group 5 chromosomes. The genetic maps were constructed for chromosome 5B in wheat and 5D in Triticum tauschii. No marker mapped in the proximal 20% chromosome region surrounding the centromere. More than 60% of the long arm markers were present in three major clusters that physically encompassed <18% of the arm. Because 48% of the markers were cDNA clones and the distributions of the cDNA and genomic clones were similar, the marker distribution may represent the distribution of genes. The gene clusters were identified and allocated to very small chromosome regions because of a higher number of deletions in their surrounding regions. The recombination was suppressed in the centromeric regions and mainly occurred in the gene-rich regions. The bp/cM estimates varied from 118 kb for gene-rich regions to 22 Mb for gene-poor regions. The wheat genes present in these clusters are, therefore, amenable to molecular manipulations parallel to the plants with smaller genomes like rice.     

Integrated physical,genetic and genome map of chickpea (Cicer arietinum L.)

Physical map of chickpea was developed for the reference chickpea genotype (ICC 4958) using bacterial artificial chromosome (BAC) libraries targeting 71,094 clones (~12× coverage). High information content fingerprinting (HICF) of these clones gave high-quality fingerprinting data for 67,483 clones, and 1,174 contigs comprising 46,112 clones and 3,256 singletons were defined. In brief, 574 Mb genome size was assembled in 1,174 contigs with an average of 0.49 Mb per contig and 3,256 singletons represent 407 Mb genome. The physical map was linked with two genetic maps with the help of 245 BAC-end sequence (BES)-derived simple sequence repeat (SSR) markers. This allowed locating some of the BACs in the vicinity of some important quantitative trait loci (QTLs) for drought tolerance and reistance to Fusarium wilt and Ascochyta blight. In addition, fingerprinted contig (FPC) assembly was also integrated with the draft genome sequence of chickpea. As a result, ~965 BACs including 163 minimum tilling path (MTP) clones could be mapped on eight pseudo-molecules of chickpea forming 491 hypothetical contigs representing 54,013,992 bp (~54 Mb) of the draft genome. Comprehensive analysis of markers in abiotic and biotic stress tolerance QTL regions led to identification of 654, 306 and 23 genes in drought tolerance “QTL-hotspot” region, Ascochyta blight resistance QTL region and Fusarium wilt resistance QTL region, respectively. Integrated physical, genetic and genome map should provide a foundation for cloning and isolation of QTLs/genes for molecular dissection of traits as well as markers for molecular breeding for chickpea improvement.     

A genetic linkage map of an apple rootstock progeny anchored to the Malus genome sequence

An apple rootstock progeny raised from the cross between the very dwarfing ??M.27?? and the more vigorous ??M.116?? (??M.M.106???×???M.27??) was used for the construction of a linkage map comprising a total of 324 loci: 252 previously mapped SSRs, 71 newly characterised or previously unmapped SSR loci (including 36 amplified by 33 out of the 35 novel markers reported here), and the self-incompatibility locus. The map spanned the 17 linkage groups (LG) expected for apple covering a genetic distance of 1,229.5?cM, an estimated 91% of the Malus genome. Linkage groups were well populated and, although marker density ranged from 2.3 to 6.2?cM/SSR, just 15 gaps of more than 15?cM were observed. Moreover, only 17.5% of markers displayed segregation distortion and, unsurprisingly in a semi-compatible backcross, distortion was particularly pronounced surrounding the self-incompatibility locus (S) at the bottom of LG17. DNA sequences of 273 SSR markers and the S locus, representing a total of 314 loci in this investigation, were used to anchor to the ??Golden Delicious?? genome sequence. More than 260 of these loci were located on the expected pseudo-chromosome on the ??Golden Delicious?? genome or on its homeologous pseudo-chromosome. In total, 282.4?Mbp of sequence from 142 genome sequence scaffolds of the Malus genome were anchored to the ??M.27???×???M.116?? map, providing an interface between the marker data and the underlying genome sequence. This will be exploited for the identification of genes responsible for traits of agronomic importance such as dwarfing and water use efficiency.     

Physical mapping of the barley stem rust resistance gene rpg4

The barley stem rust resistance gene rpg4 was physically and genetically localized on two overlapping BAC clones covering an estimated 300-kb region of the long arm of barley chromosome 7(5H). Initially, our target was mapped within a 6.0-cM region between the previously described flanking markers MWG740 and ABG391. This region was then saturated by integrating new markers from several existing barley and rice maps and by using BAC libraries of barley cv. Morex and rice cv. Nipponbare. Physical/genetic distances in the vicinity of rpg4 were found to be 1.0 Mb/cM, which is lower than the average for barley (4 Mb/cM) and lower than that determined by translocation breakpoint mapping (1.8 Mb/cM). Synteny at high resolution levels has been established between the region of barley chromosome 7(5H) containing the rpg4 locus and the subtelomeric region of rice chromosome 3 between markers S16474 and E10757. This 1.7-cM segment of the rice genome was covered by two overlapping BAC clones, about 250 kb of total length. In barley the markers S16474 and E10757 genetically delimit rpg4, lying 0.6 cM distal and 0.4 cM proximal to the locus, respectively.     

Construction of a sequence-tagged high-density genetic map of papaya for comparative structural and evolutionary genomics in brassicales

A high-density genetic map of papaya (Carica papaya L.) was constructed using microsatellite markers derived from BAC end sequences and whole-genome shot gun sequences. Fifty-four F(2) plants derived from varieties AU9 and SunUp were used for linkage mapping. A total of 707 markers, including 706 microsatellite loci and the morphological marker fruit flesh color, were mapped into nine major and three minor linkage groups. The resulting map spanned 1069.9 cM with an average distance of 1.5 cM between adjacent markers. This sequence-based microsatellite map resolved the very large linkage group 2 (LG 2) of the previous high-density map using amplified fragment length polymorphism markers. The nine major LGs of our map represent papaya's haploid nine chromosomes with LG 1 of the sex chromosome being the largest. This map validates the suppression of recombination at the male-specific region of the Y chromosome (MSY) mapped on LG 1 and at potential centromeric regions of other LGs. Segregation distortion was detected in a large region on LG 1 surrounding the MSY region due to the abortion of the YY genotype and in a region of LG6 due to an unknown cause. This high-density sequence-tagged genetic map is being used to integrate genetic and physical maps and to assign genome sequence scaffolds to papaya chromosomes. It provides a framework for comparative structural and evolutional genomic research in the order Brassicales.     

Transcriptome-based single nucleotide polymorphism markers for genome mapping in Japanese pear (Pyrus pyrifolia Nakai)

The development of single nucleotide polymorphism (SNP) markers in Japanese pear (Pyrus pyrifolia Nakai) offers the opportunity to use DNA markers for marker-assisted selection in breeding programs because of their high abundance, codominant inheritance, and potential for automated high-throughput analysis. We developed a 1,536-SNP bead array without a reference genome sequence from more than 44,000 base changes on the basis of a large-scale expressed sequence tag (EST) analysis combined with 454 genome sequencing data of Japanese pear ‘Housui’. Among the 1,536 SNPs on the array, 756 SNPs were genotyped, and 609 SNP loci were mapped to linkage groups on a genetic linkage map of ‘Housui’, based on progeny of an interspecific cross between European pear (Pyrus communis L.) ‘Bartlett’ and ‘Housui’. The newly constructed genetic linkage map consists of 951 loci, comprising 609 new SNPs, 110 pear genomic simple sequence repeats (SSRs), 25 pear EST–SSRs, 127 apple SSRs, 61 pear SNPs identified by the “potential intron polymorphism” method, and 19 other loci. The map covers 22 linkage groups spanning 1341.9 cM with an average distance of 1.41 cM between markers and is anchored to reference genetic linkage maps of European pears and apples. A total of 514 contigs containing mapped SNP loci showed significant similarity to known proteins by functional annotation analysis.     

Aligned genetic linkage maps of apple rootstock cultivar ‘JM7’ and Malus sieboldii ‘Sanashi 63’ constructed with novel EST-SSRs

Identification of markers associated with genes of interest and quantitative trait loci (QTLs), combined with high-density genetic linkage maps, can help reduce labor and costs by enabling marker-assisted selection (MAS). In this study, a dwarfing apple rootstock cultivar ??JM7?? (Malus prunifolia × Malus pumila ??Malling 9??) and wild apple Malus sieboldii ??Sanashi 63?? (section Sorbomalus) were used for constructing genetic linkage maps. Here, a species from section Sorbomalus was used for the first time as a target species in a genome-wide mapping study. We also developed and mapped 137 novel-expressed sequence tag-simple sequence repeat (EST-SSR) markers. The genetic linkage maps of ??JM7?? and ??Sanashi 63?? consisted of 415 and 310 loci and spanned 998.0 and 981.8?cM, respectively, comparable to the reference map of Malus × domestica ??Discovery??. A BLASTN search revealed that all of the EST-SSR sequences used in this study exhibited very high homology to one or more previously characterized apple genome contigs. Although the most homologous contigs of 89 EST-SSRs were located within the same linkage groups (LGs) identified by mapping analysis, the other 48 EST-SSRs were aligned into contigs positioned in different LGs than those identified by mapping. When search criteria were expanded to include the five most homologous contigs of each EST-SSR, at least one of the top five contigs for 15 of these 48 EST-SSRs corresponded to the LG obtained by mapping. The maps of ??JM7?? and ??Sanashi 63?? may be useful for analyzing important rootstock characteristics and identifying markers for MAS.