Integrating sequence with FPC fingerprint maps

Recent advances in both clone fingerprinting and draft sequencing technology have made it increasingly common for species to have a bacterial artificial clone (BAC) fingerprint map, BAC end sequences (BESs) and draft genomic sequence. The FPC (fingerprinted contigs) software package contains three modules that maximize the value of these resources. The BSS (blast some sequence) module provides a way to easily view the results of aligning draft sequence to the BESs, and integrates the results with the following two modules. The MTP (minimal tiling path) module uses sequence and fingerprints to determine a minimal tiling path of clones. The DSI (draft sequence integration) module aligns draft sequences to FPC contigs, displays them alongside the contigs and identifies potential discrepancies; the alignment can be based on either individual BES alignments to the draft, or on the locations of BESs that have been assembled into the draft. FPC also supports high-throughput fingerprint map generation as its time-intensive functions have been parallelized for Unix-based desktops or servers with multiple CPUs. Simulation results are provided for the MTP, DSI and parallelization. These features are in the FPC V9.3 software package, which is freely available.     

LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes

Background  

Physical maps are the substrate of genome sequencing and map-based cloning and their construction relies on the accurate assembly of BAC clones into large contigs that are then anchored to genetic maps with molecular markers. High Information Content Fingerprinting has become the method of choice for large and repetitive genomes such as those of maize, barley, and wheat. However, the high level of repeated DNA present in these genomes requires the application of very stringent criteria to ensure a reliable assembly with the FingerPrinted Contig (FPC) software, which often results in short contig lengths (of 3-5 clones before merging) as well as an unreliable assembly in some difficult regions. Difficulties can originate from a non-linear topological structure of clone overlaps, low power of clone ordering algorithms, and the absence of tools to identify sources of gaps in Minimal Tiling Paths (MTPs).     

SynBrowse: a synteny browser for comparative sequence analysis

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu     

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).     

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.     

Microsatellite discovery from BAC end sequences and genetic mapping to anchor the soybean physical and genetic maps.

Whole-genome sequencing of the soybean (Glycine max (L.) Merr. 'Williams 82') has made it important to integrate its physical and genetic maps. To facilitate this integration of maps, we screened 3290 microsatellites (SSRs) identified from BAC end sequences of clones comprising the 'Williams 82' physical map. SSRs were screened against 3 mapping populations. We found the AAT and ACT motifs produced the greatest frequency of length polymorphisms, ranging from 17.2% to 32.3% and from 11.8% to 33.3%, respectively. Other useful motifs include the dinucleotide repeats AG, AT, and AG, with frequency of length polymorphisms ranging from 11.2% to 18.4% (AT), 12.4% to 20.6% (AG), and 11.3% to 16.4% (GT). Repeat lengths less than 16 bp were generally less useful than repeat lengths of 40-60 bp. Two hundred and sixty-five SSRs were genetically mapped in at least one population. Of the 265 mapped SSRs, 60 came from BAC singletons not yet placed into contigs of the physical map. One hundred and ten originated in BACs located in contigs for which no genetic map location was previously known. Ninety-five SSRs came from BACs within contigs for which one or more other BACs had already been mapped. For these fingerprinted contigs (FPC) a high percentage of the mapped markers showed inconsistent map locations. A strategy is introduced by which physical and genetic map inconsistencies can be resolved using the preliminary 4x assembly of the whole genome sequence of soybean.     

High-resolution comprehensive radiation hybrid maps of the porcine chromosomes 2p and 9p compared with the human chromosome 11

We are constructing high-resolution, chromosomal 'test' maps for the entire pig genome using a 12,000-rad WG-RH panel (IMNpRH2(12,000-rad))to provide a scaffold for the rapid assembly of the porcine genome sequence. Here we present an initial, comparative map of human chromosome (HSA) 11 with pig chromosomes (SSC) 2p and 9p. Two sets of RH mapping vectors were used to construct the RH framework (FW) maps for SSC2p and SSC9p. One set of 590 markers, including 131 microsatellites (MSs), 364 genes/ESTs, and 95 BAC end sequences (BESs) was typed on the IMNpRH2(12,000-rad) panel. A second set of 271 markers (28 MSs, 138 genes/ESTs, and 105 BESs) was typed on the IMpRH(7,000-rad) panel. The two data sets were merged into a single data-set of 655 markers of which 206 markers were typed on both panels. Two large linkage groups of 72 and 194 markers were assigned to SSC2p, and two linkage groups of 84 and 168 markers to SSC9p at a two-point LOD score of 10. A total of 126 and 114 FW markers were ordered with a likelihood ratio of 1000:1 to the SSC2p and SSC9p RH(12,000-rad) FW maps, respectively, with an accumulated map distance of 4046.5 cR(12,000 )and 1355.2 cR(7,000 )for SSC2p, and 4244.1 cR(12,000) and 1802.5 cR(7,000) for SSC9p. The kb/cR ratio in the IMNpRH2(12,000-rad) FW maps was 15.8 for SSC2p, and 15.4 for SSC9p, while the ratio in the IMpRH(7,000-rad) FW maps was 47.1 and 36.3, respectively, or an approximately 3.0-fold increase in map resolution in the IMNpRH(12,000-rad) panel over the IMpRH(7,000-rad) panel. The integrated IMNpRH(12,000-rad) andIMpRH(7,000-rad) maps as well as the genetic and BAC FPC maps provide an inclusive comparative map between SSC2p, SSC9p and HSA11 to close potential gaps between contigs prior to sequencing, and to identify regions where potential problems may arise in sequence assembly.     

A Single Molecule Scaffold for the Maize Genome

About 85% of the maize genome consists of highly repetitive sequences that are interspersed by low-copy, gene-coding sequences. The maize community has dealt with this genomic complexity by the construction of an integrated genetic and physical map (iMap), but this resource alone was not sufficient for ensuring the quality of the current sequence build. For this purpose, we constructed a genome-wide, high-resolution optical map of the maize inbred line B73 genome containing >91,000 restriction sites (averaging 1 site/∼23 kb) accrued from mapping genomic DNA molecules. Our optical map comprises 66 contigs, averaging 31.88 Mb in size and spanning 91.5% (2,103.93 Mb/∼2,300 Mb) of the maize genome. A new algorithm was created that considered both optical map and unfinished BAC sequence data for placing 60/66 (2,032.42 Mb) optical map contigs onto the maize iMap. The alignment of optical maps against numerous data sources yielded comprehensive results that proved revealing and productive. For example, gaps were uncovered and characterized within the iMap, the FPC (fingerprinted contigs) map, and the chromosome-wide pseudomolecules. Such alignments also suggested amended placements of FPC contigs on the maize genetic map and proactively guided the assembly of chromosome-wide pseudomolecules, especially within complex genomic regions. Lastly, we think that the full integration of B73 optical maps with the maize iMap would greatly facilitate maize sequence finishing efforts that would make it a valuable reference for comparative studies among cereals, or other maize inbred lines and cultivars.     

Integrated cytogenetic map of mitotic metaphase chromosome 9 of maize: resolution,sensitivity, and banding paint development

To study the correlation of the sequence positions on the physical DNA finger print contig (FPC) map and cytogenetic maps of pachytene and somatic maize chromosomes, sequences located along the chromosome 9 FPC map approximately every 10 Mb were selected to place on maize chromosomes using fluorescent in situ hybridization (FISH). The probes were produced as pooled polymerase chain reaction products based on sequences of genetic markers or repeat-free portions of mapped bacterial artificial chromosome (BAC) clones. Fifteen probes were visualized on chromosome 9. The cytological positions of most sequences correspond on the pachytene, somatic, and FPC maps except some probes at the pericentromeric regions. Because of unequal condensation of mitotic metaphase chromosomes, being lower at pericentromeric regions and higher in the arms, probe positions are displaced to the distal ends of both arms. The axial resolution of FISH on somatic chromosome 9 varied from 3.3 to 8.2 Mb, which is 12-30 times lower than on pachytene chromosomes. The probe collection can be used as chromosomal landmarks or as a "banding paint" for the physical mapping of sequences including transgenes and BAC clones and for studying chromosomal rearrangements.     

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.     

iMap: a database-driven utility to integrate and access the genetic and physical maps of maize

MOTIVATION: Because of the unique biological features, a bioinformatic platform for the integrated genetic and physical map of maize is required for storing, integrating, accessing and visualizing the underlying data. RESULTS: The goal of the Maize Mapping Project is to develop a fully integrated genetic and physical map for maize. To display this integrated map, we have developed iMap. iMap has three main components: a relational database (iMapDB), a map graphic browser (iMap Viewer) and a search utility (iMap Search). iMapDB is populated with current genetic and physical map data, describing relationships among genetic loci, molecular markers and bacterial artificial chromosome (BAC) contigs. The database also contains integrated information produced by applying a set of anchoring rules to assign BAC contigs to specific locations on the genetic map. The iMap Viewer and iMap Search functions are combined in the user interface to allow viewing and retrieving many types of genetic and physical map data. The iMap Viewer features side-by-side chromosome-based displays of the genetic map and associated BAC contigs. For each genetic locus, information about marker type or contig can be viewed via pop-up windows that feature links to external data resources. Searches can be conducted for genetic locus, probe or sequence accession number; search results include relevant map positions, anchored BAC contigs and links to the graphical display of relevant chromosomes. iMap can be accessed at http://www.maizemap.org AVAILABILITY: The iMap utility package is available for non-commercial use upon request from the authors.     

The Mouse Genome Project and human genetics. A report from the 5th International Mouse Genome Mapping Workshop, Lunteren, Holland.

Genome-wide mapping efforts are moving toward the establishment of a 1-cM genetic map of the entire mouse genome. The bulk of linkage groups conserved between the mouse and the human genomes has been identified. Microsatellite mapping has had a major impact on the development of genome-wide genetic maps and, in particular, on genome-wide searches for polygenic disease loci. Some substantial regions of the mouse genome have a marker density of 1 cM or less and many of these regions are now physically mapped. Embryonic YAC contigs have been established in some physically mapped regions. A unitary, global mouse mapping database--the Mouse Genome Database--is under development along with associated software tools. Chromosome committees are having a major impact on the establishment and verification of chromosome maps through the preparation of published annual reports.     

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.     

Genome Annotation Resource Fields--GARFIELD: a genome browser for Felis catus

Annotation features from the 1.9-fold whole-genome shotgun (WGS) sequences of domestic cat have been organized into an interactive web application, Genome Annotation Resource Fields (GARFIELD) (http://lgd.abcc.ncifcrf.gov) at the Laboratory of Genomic Diversity and Advanced Biomedical Computing Center (ABCC) at The National Cancer Institute (NCI). The GARFIELD browser allows the user to view annotations on a per chromosome basis with unplaced contigs provided on placeholder chromosomes. Various tracks on the browser allow display of annotations. A Genes track on the browser includes 20 285 regions that align to genes annotated in other mammalian genomes: Homo sapiens, Pan troglodytes, Mus musculus, Rattus norvegicus, Bos taurus, and Canis familiaris. Also available are tracks that display the contigs that make up the chromosomes and representations of their GC content and repetitive elements as detected using the RepeatMasker (http://www.repeatmasker.org). Data from the browser can be downloaded in FASTA and GFF format, and users can upload their own data to the display. The Felis catus sequences and their chromosome assignments and additional annotations incorporate data analyzed and produced by a multicenter collaboration between NCI, ABCC, Agencourt Biosciences Corporation, Broad Institute of Harvard and Massachusetts Institute of Technology, National Human Genome Research Institute, National Center for Biotechnology and Information, and Texas A&M.     

Expression analysis and physical mapping of low-molecular-weight glutenin loci in hexaploid wheat (Triticum aestivum L.).

Gliadins and glutenins are storage proteins important in determining the bread-, noodle-, and pasta-making quality of wheat. Glutenins consist of HMW and LMW subunits. The Glu-A3, Glu-B3, and Glu-D3 loci on the short arms of chromosomes 1A, 1B, and 1D, respectively, are the major loci for LMW glutenins. To construct physical maps of the Glu-3 loci, a set of 24 high-density filters representing a 3.1x genome coverage hexaploid wheat BAC library was screened by hybridization using a probe made of 3 LMW glutenin sequences. After 2 rounds of hybridization, a subset of 536 BAC clones were selected and fingerprinted. Three developing seed cDNA libraries were also constructed. A total of 5000-6000 ESTs were generated from each library, assembled into contigs and searched by homology for LMW glutenin sequences. In total, 90 full-length LMW glutenin sequences were found to cluster into 8 distinct groups representing at least 21 different LMW glutenin subunits. A set of 24 pairs of PCR primers was designed from these groups and used as markers on the BAC clones. The combined fingerprinting and marker data were used to build the physical maps using FPC software. A total of 91 contigs comprising 254 clones were obtained and 282 clones remained singletons.     

Extended cytogenetic maps of sheep chromosome 1 and their cattle and river buffalo homoeologues: comparison with the OAR1 RH map and human chromosomes 2, 3, 21 and 1q

Cytogenetic maps are useful tools for several applications, such as the physical anchoring of linkage and RH maps or genome sequence contigs to specific chromosome regions or the analysis of chromosome rearrangements. Recently, a detailed RH map was reported in OAR1. In the present study, we selected 38 markers equally distributed in this RH map for identification of ovine genomic DNA clones within the ovine BAC library CHORI-243 using the virtual sheep genome browser and performed FISH mapping for both comparison of OAR1 and homoeologous chromosomes BBU1q-BBU6 and BTA1-BTA3 and considerably extending the cytogenetic maps of the involved species-specific chromosomes. Comparison of the resulting maps with human-identified homology with HSA2q, HSA3, HSA21 and HSA1q reveals complex chromosome rearrangements differentiating human and bovid chromosomes. In addition, we identified 2 new small human segments from HSA2q and HSA3q conserved in the telomeric regions of OAR1p and homoeologous chromosome regions of BTA3 and BBU6, and OAR1q, respectively. Evaluation of the present OAR1 cytogenetic map and the OAR1 RH map supports previous RH assignments with 2 main exceptions. The 2 loci BMS4011 and CL638002 occupy inverted positions in these 2 maps.