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.     

BAC-HAPPY Mapping (BAP Mapping): A New and Efficient Protocol for Physical Mapping

Physical and linkage mapping underpin efforts to sequence and characterize the genomes of eukaryotic organisms by providing a skeleton framework for whole genome assembly. Hitherto, linkage and physical “contig” maps were generated independently prior to merging. Here, we develop a new and easy method, BAC HAPPY MAPPING (BAP mapping), that utilizes BAC library pools as a HAPPY mapping panel together with an Mbp-sized DNA panel to integrate the linkage and physical mapping efforts into one pipeline. Using Arabidopsis thaliana as an exemplar, a set of 40 Sequence Tagged Site (STS) markers spanning ∼10% of chromosome 4 were simultaneously assembled onto a BAP map compiled using both a series of BAC pools each comprising 0.7x genome coverage and dilute (0.7x genome) samples of sheared genomic DNA. The resultant BAP map overcomes the need for polymorphic loci to separate genetic loci by recombination and allows physical mapping in segments of suppressed recombination that are difficult to analyze using traditional mapping techniques. Even virtual “BAC-HAPPY-mapping” to convert BAC landing data into BAC linkage contigs is possible.     

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.     

Application of alternative models to identify QTL for growth traits in an F2 Duroc x Pietrain pig resource population

Background

Several lines of evidence including allozyme analysis, restriction digest patterns and sequencing of mtDNA as well as mini- and micro-satellite allele frequencies indicate that Atlantic salmon (Salmo salar) from North America and Europe are genetically distinct. These observations are supported by karyotype analysis, which revealed that North American Atlantic salmon have 27 pairs of chromosomes whereas European salmon have 29 pairs. We set out to construct a linkage map for a North American Atlantic salmon family and to compare this map with the well developed map for European Atlantic salmon.

Results

We used microsatellite markers, which had previously been mapped in the two Atlantic salmon SALMAP mapping families from the River Tay, Scotland, to carry out linkage analysis in an Atlantic salmon family (NB1) whose parents were derived from the Saint John River stock in New Brunswick, Canada. As large differences in recombination rates between female and male Atlantic salmon have been noted, separate genetic maps were constructed for each sex. The female linkage map comprises 218 markers in 37 linkage groups while the male map has 226 markers in 28 linkage groups. We combined 280 markers from the female and male maps into 27 composite linkage groups, which correspond to the haploid number of chromosomes in Atlantic salmon from the Western Atlantic.

Conclusions

A comparison of the composite NB1 and SALMAP linkage maps revealed the reason for the difference in the chromosome numbers between European and North American Atlantic salmon: Linkage groups AS-4 and AS-32 in the Scottish salmon, which correspond to chromosomes Ssa-6 and Ssa-22, are combined into a single NB1 linkage group as are linkage groups AS-21 and AS-33 (corresponding to chromosomes Ssa-26 and Ssa-28). The comparison of the linkage maps also suggested some additional chromosomal rearrangements, but it will require finer mapping, potentially using SNPs, to test these predictions. Our results provide the first comparison of the genomic architecture of Atlantic salmon from North America and Europe with respect to chromosome organization.     

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.     

An integrated map of Oryza sativa L. chromosome 5

The developments of molecular marker-based genetic linkage maps are now routine. Physical maps based on contigs of large insert genomic clones have been established in several plant species. However, integration of genetic, physical, and cytological maps is still a challenge for most plant species. Here we present an integrated map of rice (Oryza sativa L.) chromosome 5, developed by fluorescence in situ hybridization mapping of 18 bacterial artificial chromosome (BAC) clones or PI-derived artificial chromosome (PAC) clones on meiotic pachytene chromosomes. Each BAC/PAC clone was anchored by a restriction fragment length polymorphism marker mapped to the rice genetic linkage map. This molecular cytogenetic map shows the genetic recombination and sequence information of a physical map, correlated to the cytological features of rice chromosome 5. Detailed comparisons of the distances between markers on genetic, cytological, and physical maps, revealed the distributions of recombination events and molecular organization of the chromosomal features of rice chromosome 5 at the pachytene stage. Discordance of distances between the markers was found among the different maps. Our results revealed that neither the recombination events nor the degree of chromatin condensation were evenly distributed along the entire length of chromosome 5. Detailed comparisons of the correlative positions of markers on the genetic, cytological, and physical maps of rice chromosome 5 provide insight into the molecular architecture of rice chromosome 5, in relation to its cytological features and recombination events on the genetic map. The prospective applications of such an integrated cytogenetic map are discussed.     

A linkage map of Atlantic salmon (Salmo salar) reveals an uncommonly large difference in recombination rate between the sexes

A genetic linkage map of the Atlantic salmon (Salmo salar) was constructed, using 54 microsatellites and 473 amplified fragment length polymorphism (AFLP) markers. The mapping population consisted of two full-sib families within one paternal half-sib family from the Norwegian breeding population. A mapping strategy was developed that facilitated the construction of separate male and female maps, while retaining all the information contributed by the dominant AFLP markers. By using this strategy, we were able to map a significant number of the AFLP markers for which all informative offspring had two heterozygous parents; these markers then served as bridges between the male and female maps. The female map spanned 901 cM and had 33 linkage groups, while the male spanned 103 cM and had 31 linkage groups. Twenty-five linkage groups were common between the two maps. The construction of the genetic map revealed a large difference in recombination rate between females and males. The ratio of female recombination rate vs. male recombination rate was 8.26, the highest ratio reported for any vertebrate. This map constitutes the first linkage map of Atlantic salmon, one of the most important aquaculture species worldwide.     

Identification of the sex-determining locus of Atlantic salmon (Salmo salar) on chromosome 2

We have integrated data from linkage mapping, physical mapping and karyotyping to gain a better understanding of the sex-determining locus, SEX, in Atlantic salmon (Salmo salar). SEX has been mapped to Atlantic salmon linkage group 1 (ASL1) and is associated with several microsatellite markers. We have used probes designed from the flanking regions of these sex-linked microsatellite markers to screen a bacterial artificial chromosome (BAC) library, representing an 11.7x coverage of the Atlantic salmon genome, which has been HindIII fingerprinted and assembled into contigs. BACs containing sex-linked microsatellites and their related contigs have been identified and representative BACs have been placed on the Atlantic salmon chromosomes by fluorescent in situ hybridization (FISH). This identified chromosome 2, a large metacentric, as the sex chromosome. By positioning several BACs on this chromosome by FISH, it was possible to orient ASL1 with respect to chromosome 2. The region containing SEX appears to lie on the long arm between marker Ssa202DU and a region of heterochromatin identified by DAPI staining. BAC end-sequencing of clones within sex-linked contigs revealed five hitherto unmapped genes along the sex chromosome. We are using an in silico approach coupled with physical probing of the BAC library to extend the BAC contigs to provide a physical map of ASL1, with a view to sequencing chromosome 2 and, in the process, identifying the sex-determining gene.     

Construction of a high-density,high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.     

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.     

Integration of physical and genetic maps in apple confirms whole-genome and segmental duplications in the apple genome

A total of 355 simple sequence repeat (SSR) markers were developed, based on expressed sequence tag (EST) and bacterial artificial chromosome (BAC)-end sequence databases, and successfully used to construct an SSR-based genetic linkage map of the apple. The consensus linkage map spanned 1143 cM, with an average density of 2.5 cM per marker. Newly developed SSR markers along with 279 SSR markers previously published by the HiDRAS project were further used to integrate physical and genetic maps of the apple using a PCR-based BAC library screening approach. A total of 470 contigs were unambiguously anchored onto all 17 linkage groups of the apple genome, and 158 contigs contained two or more molecular markers. The genetically mapped contigs spanned ～421 Mb in cumulative physical length, representing 60.0% of the genome. The sizes of anchored contigs ranged from 97 kb to 4.0 Mb, with an average of 995 kb. The average physical length of anchored contigs on each linkage group was ～24.8 Mb, ranging from 17.0 Mb to 37.73 Mb. Using BAC DNA as templates, PCR screening of the BAC library amplified fragments of highly homologous sequences from homoeologous chromosomes. Upon integrating physical and genetic maps of the apple, the presence of not only homoeologous chromosome pairs, but also of multiple locus markers mapped to adjacent sites on the same chromosome was detected. These findings demonstrated the presence of both genome-wide and segmental duplications in the apple genome and provided further insights into the complex polyploid ancestral origin of the apple.     

Chromosomal differences between European and North American Atlantic salmon discovered by linkage mapping and supported by fluorescence in situ hybridization analysis

ABSTRACT: BACKGROUND: Geographical isolation has generated a distinct difference between Atlantic salmon of European and North American Atlantic origin. The European Atlantic salmon generally has 29 pairs of chromosomes and 74 chromosome arms whereas it has been reported that the North American Atlantic salmon has 27 chromosome pairs and an NF of 72. In order to predict the major chromosomal rearrangements causing these differences, we constructed a dense linkage map for Atlantic salmon of North American origin and compared it with the well-developed map for European Atlantic salmon. RESULTS: The presented male and female genetic maps for the North American subspecies of Atlantic salmon, contains 3,662 SNPs located on 27 linkage groups. The total lengths of the female and male linkage maps were 2,153 cM and 968 cM respectively, with males characteristically showing recombination only at the telomeres. We compared these maps with recently published SNP maps from European Atlantic salmon, and predicted three chromosomal reorganization events that we then tested using fluorescence in situ hybridization (FISH) analysis. The proposed rearrangements, which define the differences in the karyotypes of the North American Atlantic salmon relative to the European Atlantic salmon, include the translocation of the p arm of ssa01 to ssa23 and polymorphic fusions: ssa26 with ssa28, and ssa08 with ssa29. CONCLUSIONS: This study identified major chromosomal differences between European and North American Atlantic salmon. However, while gross structural differences were significant, the order of genetic markers at the fine-resolution scale was remarkably conserved. This is a good indication that information from the International Cooperation to Sequence the Atlantic salmon Genome, which is sequencing a European Atlantic salmon, can be transferred to Atlantic salmon from North America.     

Palaeohexaploid ancestry for Caryophyllales inferred from extensive gene-based physical and genetic mapping of the sugar beet genome (Beta vulgaris)

Sugar beet (Beta vulgaris) is an important crop plant that accounts for 30% of the world's sugar production annually. The genus Beta is a distant relative of currently sequenced taxa within the core eudicotyledons; the genomic characterization of sugar beet is essential to make its genome accessible to molecular dissection. Here, we present comprehensive genomic information in genetic and physical maps that cover all nine chromosomes. Based on this information we identified the proposed ancestral linkage groups of rosids and asterids within the sugar beet genome. We generated an extended genetic map that comprises 1127 single nucleotide polymorphism markers prepared from expressed sequence tags and bacterial artificial chromosome (BAC) end sequences. To construct a genome-wide physical map, we hybridized gene-derived oligomer probes against two BAC libraries with 9.5-fold cumulative coverage of the 758 Mbp genome. More than 2500 probes and clones were integrated both in genetic maps and the physical data. The final physical map encompasses 535 chromosomally anchored contigs that contains 8361 probes and 22 815 BAC clones. By using the gene order established with the physical map, we detected regions of synteny between sugar beet (order Caryophyllales) and rosid species that involves 1400-2700 genes in the sequenced genomes of Arabidopsis, poplar, grapevine, and cacao. The data suggest that Caryophyllales share the palaeohexaploid ancestor proposed for rosids and asterids. Taken together, we here provide extensive molecular resources for sugar beet and enable future high-resolution trait mapping, gene identification, and cross-referencing to regions sequenced in other plant species.     

A comparative genetic map of the turkey genome

Genetic markers (microsatellites and SNPs) were used to create and compare maps of the turkey and chicken genomes. A physical map of the chicken genome was built by comparing sequences of turkey markers with the chicken whole-genome sequence by BLAST analysis. A genetic linkage map of the turkey genome (Meleagris gallopavo) was developed by segregation analysis of genetic markers within the University of Minnesota/Nicholas Turkey Breeding Farms (UMN/NTBF) resource population. This linkage map of the turkey genome includes 314 loci arranged into 29 linkage groups. An additional 40 markers are tentatively placed within linkage groups based on two-point LOD scores and 16 markers remain unlinked. Total map distance contained within linkage groups is 2,011 cM with the longest linkage group (47 loci) measuring 413.3 cM. Average marker interval over the 29 linkage groups was 6.4 cM. All but one turkey linkage group could be aligned with the physical map of the chicken genome. The present genetic map of the turkey provides a comparative framework for future genomic studies.     

Comparative Genomics Identifies Candidate Genes for Infectious Salmon Anemia (ISA) Resistance in Atlantic Salmon (<Emphasis Type="Italic">Salmo salar</Emphasis>)

Infectious salmon anemia (ISA) has been described as the hoof and mouth disease of salmon farming. ISA is caused by a lethal and highly communicable virus, which can have a major impact on salmon aquaculture, as demonstrated by an outbreak in Chile in 2007. A quantitative trait locus (QTL) for ISA resistance has been mapped to three microsatellite markers on linkage group (LG) 8 (Chr 15) on the Atlantic salmon genetic map. We identified bacterial artificial chromosome (BAC) clones and three fingerprint contigs from the Atlantic salmon physical map that contains these markers. We made use of the extensive BAC end sequence database to extend these contigs by chromosome walking and identified additional two markers in this region. The BAC end sequences were used to search for conserved synteny between this segment of LG8 and the fish genomes that have been sequenced. An examination of the genes in the syntenic segments of the tetraodon and medaka genomes identified candidates for association with ISA resistance in Atlantic salmon based on differential expression profiles from ISA challenges or on the putative biological functions of the proteins they encode. One gene in particular, HIV-EP2/MBP-2, caught our attention as it may influence the expression of several genes that have been implicated in the response to infection by infectious salmon anemia virus (ISAV). Therefore, we suggest that HIV-EP2/MBP-2 is a very strong candidate for the gene associated with the ISAV resistance QTL in Atlantic salmon and is worthy of further study.     

Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait

Black poplar (Populus nigra L.) is a tree of ecological and economic interest. A better knowledge of P. nigra genome is needed for an effective protection and use of its genetic resources. The main objective of this study is the construction of a highly informative genetic map of P. nigra species including genes of adaptive and economic interest. Two genotypes originated from contrasted natural Italian populations were crossed to generate a F₁ mapping pedigree of 165 individuals. Amplification fragment length polymorphism (AFLP), simple sequence repeat (SSR), and single nucleotide polymorphism (SNP) markers were used to genotype 92 F₁ individuals, and the pseudo-test-cross strategy was applied for linkage analysis. The female parent map included 368 markers (274 AFLPs, 91 SSRs, and 3 SNPs) and spanned 2,104 cM with 20 linkage groups, and the male parent map, including 317 markers (205 AFLPs, 106 SSRs, 5 SNPs, and sex trait), spanned 2,453 cM with 23 main linkage groups. The sex, as morphological trait, was mapped on the linkage group XIX of the male parent map. The generated maps are among the most informative in SSRs when compared to the Populus maps published so far and allow a complete alignment with the 19 haploid chromosomes of Populus sequence genome. These genetic maps provide informative tools for a better understanding of P. nigra genome structure and genetic improvement of this ecologically and economically important European tree species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Construction and characterization of a BAC library from the Coffea arabica genotype Timor Hybrid CIFC 832/2

Most of the world’s coffee production originates from Coffea arabica, an allotetraploid species with low genetic diversity and for which few genomic resources are available. Genomic libraries with large DNA fragment inserts are useful tools for the study of plant genomes, including the production of physical maps, integration studies of physical and genetic maps, genome structure analysis and gene isolation by positional cloning. Here, we report the construction and characterization of a Bacterial Artificial Chromosome (BAC) library from C. arabica Timor Hybrid CIFC 832/2, a parental genotype for several modern coffee cultivars. The BAC library consists of 56,832 clones with an average insert size of 118 kb, which represents a dihaploid genome coverage of five to sixfold. The content of organellar DNA was estimated at 1.04 and 0.5 % for chloroplast and mitochondrial DNA, respectively. The BAC library was screened for the NADPH-dependent mannose-6-phosphate reductase gene (CaM6PR) with markers positioned on four linkage groups of a partial C. arabica genetic map. A mixed approach using PCR and membrane hybridization of BAC pools allowed for the discovery of nine BAC clones with the CaM6PR gene and 53 BAC clones that were anchored to the genetic map with simple sequence repeat markers. This library will be a useful tool for future studies on comparative genomics and the identification of genes and regulatory elements controlling major traits in this economically important crop species.     

Development of Chromosome-specific Cytogenetic Markers and Merging of Linkage Fragments in Papaya

Carica papaya L. is a tropical and sub-tropical fruit-tree crop with a small genome and nine pairs of chromosomes. The transgenic cultivar ‘SunUp’ has been sequenced and three high-density genetic maps are available for mapping agronomically and economically-important traits. However, the small size and similar morphology of papaya chromosomes hinder their identification and few cytological resources are available for integration of genetic and cytogenetic information. Fluorescence in situ hybridization (FISH) was performed on mitotic metaphase chromosomes using BAC clones harboring mapped simple sequence repeat (SSR) markers as probes. A total of 104 BAC clones covering all 12 linkage groups (LGs) were tested and 12 of them, that gave a single specific signal, were chosen as representative of the 12 LGs of the SSR genetic map. This set of chromosome-specific DNA markers acted as a foundation for papaya chromosome karyotyping and re-assigning orientation of LGs. Chromosome-specific markers allowed us to assign the minor LGs 10, 11, and 12 to major LGs 8, 9, and 7, respectively. We thus reduced the number of LGs in the genetic map to nine, corresponding to the haploid number of papaya chromosomes. We also tested the relative order of DNA markers on minor LGs 10 and 11 to place them on top of LGs 8 and 9 in the correct orientation. Ribosomal DNAs (rDNAs), a set of major cytogenetic markers, were positioned on specific papaya chromosomes. The 25S rDNA showed strong signals at the constriction site of a single pair of chromosomes identified as LG 2 by LG 2-specific BAC clone. The 5S rDNA showed strong signals on two pairs of chromosomes that are syntenic with LG 4- and LG 5-specific BAC clones. This integrated map will facilitate genome assembly, quantitative trait locus (QTL) mapping, and the study of cytological, physical and genetic distance relationships between papaya chromosomes.