Integrating physical and genetic maps: from genomes to interaction networks

Physical and genetic mapping data have become as important to network biology as they once were to the Human Genome Project. Integrating physical and genetic networks currently faces several challenges: increasing the coverage of each type of network; establishing methods to assemble individual interaction measurements into contiguous pathway models; and annotating these pathways with detailed functional information. A particular challenge involves reconciling the wide variety of interaction types that are currently available. For this purpose, recent studies have sought to classify genetic and physical interactions along several complementary dimensions, such as ordered versus unordered, alleviating versus aggravating, and first versus second degree.     

Integration of the genetic and physical maps of the chicken macrochromosomes

A large amount of genetic mapping information has been obtained in the chicken from the East Lansing, Compton and Wageningen reference populations. Physical mapping information has however, been more limited. We have mapped 14 new clones, both genetically and physically, and all 14 have been assigned to macrochromosomes. The orientation of linkage groups E01C01C11W01 (Chr 1), E06C02W02 (Chr 2), E02C03W03 (Chr 3), E05C04W04 (Chr 4), E07E34C05W05 (Chr 5), E11C10W06 (Chr 6), E45C07W07 (Chr 7) and E43C12W11 (Chr 8) has been established. Here we present integrated maps of the eight macrochromosomes and the Z chromosome of the chicken and correlate genetic with physical distances for chromosomes 1-3 and the Z sex chromosome.     

Alignment of genetic and physical maps of Gibberella zeae

We previously published a genetic map of Gibberella zeae (Fusarium graminearum sensu lato) based on a cross between Kansas strain Z-3639 (lineage 7) and Japanese strain R-5470 (lineage 6). In this study, that genetic map was aligned with the third assembly of the genomic sequence of G. zeae strain PH-1 (lineage 7) using seven structural genes and 108 sequenced amplified fragment length polymorphism markers. Several linkage groups were combined based on the alignments, the nine original linkage groups were reduced to six groups, and the total size of the genetic map was reduced from 1,286 to 1,140 centimorgans. Nine supercontigs, comprising 99.2% of the genomic sequence assembly, were anchored to the genetic map. Eight markers (four markers from each parent) were not found in the genome assembly, and four of these markers were closely linked, suggesting that >150 kb of DNA sequence is missing from the PH-1 genome assembly. The alignments of the linkage groups and supercontigs yielded four independent sets, which is consistent with the four chromosomes reported for this fungus. Two proposed heterozygous inversions were confirmed by the alignments; otherwise, the colinearity of the genetic and physical maps was high. Two of four regions with segregation distortion were explained by the two selectable markers employed in making the cross. The average recombination rates for each chromosome were similar to those previously reported for G. zeae. Despite an inferred history of genetic isolation of lineage 6 and lineage 7, the chromosomes of these lineages remain homologous and are capable of recombination along their entire lengths, even within the inversions. This genetic map can now be used in conjunction with the physical sequence to study phenotypes (e.g., fertility and fitness) and genetic features (e.g., centromeres and recombination frequency) that do not have a known molecular signature in the genome.     

A graphical representation of genetic and physical maps: the Marey map.

A novel, simultaneous, visual representation of sex-specific genetic maps and physical maps is introduced. Such maps, called Marey maps, provide direct comparisons of multiple genetic maps and elucidate the relationship of recombination frequency to physical distance.     

Comparative visualization of genetic and physical maps with Strudel

Data visualization can play a key role in comparative genomics, for example, underpinning the investigation of conserved synteny patterns. Strudel is a desktop application that allows users to easily compare both genetic and physical maps interactively and efficiently. It can handle large datasets from several genomes simultaneously, and allows all-by-all comparisons between these. Availability and implementation: Installers for Strudel are available for Windows, Linux, Solaris and Mac OS X at http://bioinf.scri.ac.uk/strudel/.     

ODS2: a multiplatform software application for creating integrated physical and genetic maps

A contig map is a physical map that shows the native order of a library of overlapping genomic clones. One common method for creating such maps involves using hybridization to detect clone overlaps. False- positive and false-negative hybridization errors, the presence of chimeric clones, and gaps in library coverage lead to ambiguity and error in the clone order. Genomes with good genetic maps, such as Neurospora crassa, provide a means for reducing ambiguities and errors when constructing contig maps if clones can be anchored with genetic markers to the genetic map. A software application called ODS2 for creating contig maps based on clone-clone hybridization data is presented. This application is also designed to exploit partial ordering information provided by anchorage of clones to a genetic map. This information, along with clone-clone hybridization data, is used by a clone ordering algorithm and is represented graphically, allowing users to interactively align physical and genetic maps. ODS2 has a graphical user interface and is implemented entirely in Java, so it runs on multiple platforms. Other features include the flexibility of storing data in a local file or relational database and the ability to create full or minimum tiling contig maps.     

Integration of the Aedes aegypti mosquito genetic linkage and physical maps

Two approaches were used to correlate the Aedes aegypti genetic linkage map to the physical map. STS markers were developed for previously mapped RFLP-based genetic markers so that large genomic clones from cosmid libraries could be found and placed to the metaphase chromosome physical maps using standard FISH methods. Eight cosmids were identified that contained eight RFLP marker sequences, and these cosmids were located on the metaphase chromosomes. Twenty-one cDNAs were mapped directly to metaphase chromosomes using a FISH amplification procedure. The chromosome numbering schemes of the genetic linkage and physical maps corresponded directly and the orientations of the genetic linkage maps for chromosomes 2 and 3 were inverted relative to the physical maps. While the chromosome 2 linkage map represented essentially 100% of chromosome 2, approximately 65% of the chromosome 1 linkage map mapped to only 36% of the short p-arm and 83% of the chromosome 3 physical map contained the complete genetic linkage map. Since the genetic linkage map is a RFLP cDNA-based map, these data also provide a minimal estimate for the size of the euchromatic regions. The implications of these findings on positional cloning in A. aegypti are discussed.     

Spline methods for the comparison of physical and genetic maps.

The first genetic maps were constructed by linkage analysis. Physical mapping techniques, such as radiation hybrids and complete sequencing, produce a different picture. For the purposes of population genetics, clinical genetics, and genetic epidemiology, it is important to harmonize and amalgamate existing genetic and physical maps. Among other things, comparisons of the two kinds of maps promotes better understanding of the wide variation in local recombination rates per unit physical length of DNA. The current paper presents methods for estimating recombination intensity as a function of physical distance along a chromosome. Genetic map distance is the integral of intensity. We derive fast reliable estimation algorithms based on a Poisson process model, penalized likelihoods, and cubic spline interpolation. Our methods provide a rigorous and statistically sound foundation for comparing physical and genetic maps. To illustrate the possibilities, we apply the methods to published recombination data on CEPH families and the complete sequences of chromosomes 21 and 22. Our results are in good agreement with previous studies and the biological data.     

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.     

Map error reduction: using genetic and sequence-based physical maps to order closely linked markers

The Marshfield comprehensive genetic maps are frequently used for linkage and association studies, however, for some regions of these maps the marker order has low level of likelihood ratio support. In order to investigate the level of statistical support and the accuracy of the genetic maps compared to sequence-based physical maps, two approximately 30 cM autosomal regions were selected. The first region was selected from chromosome 3 and consisted predominately of draft sequence. The second region was selected from chromosome 21 and consisted of finished sequence data. The physical order of these markers was based upon their position on Celera (CEL) and Human Genome Project-Santa Cruz (HGP-sc) sequence-based physical maps. The chromosome 3 and 21 regions contained 100 and 61 markers, respectively, on the Marshfield genetic map. The genetic and physical map order was consistent for 88.9 and 89.2% of the markers in the region on chromosome 3 and 21, respectively. Using a novel scoring criterion to assess inconsistent marker order between genetic and physical maps, it was determined that the physical order was likely the correct order for 3.3 and 7.1% of the markers in the chromosome 3 and 21 regions, respectively. To increase the accuracy of the order of markers selected for fine mapping a method is presented which combines information from genetic and sequence-based physical maps.     

A graph-theoretic approach to comparing and integrating genetic, physical and sequence-based maps

For many species, multiple maps are available, often constructed independently by different research groups using different sets of markers and different source material. Integration of these maps provides a higher density of markers and greater genome coverage than is possible using a single study. In this article, we describe a novel approach to comparing and integrating maps by using abstract graphs. A map is modeled as a directed graph in which nodes represent mapped markers and edges define the order of adjacent markers. Independently constructed graphs representing corresponding maps from different studies are merged on the basis of their common loci. Absence of a path between two nodes indicates that their order is undetermined. A cycle indicates inconsistency among the mapping studies with regard to the order of the loci involved. The integrated graph thus produced represents a complete picture of all of the mapping studies that comprise it, including all of the ambiguities and inconsistencies among them. The objective of this representation is to guide additional research aimed at interpreting these ambiguities and inconsistencies in locus order rather than presenting a "consensus order" that ignores these problems.     

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.     

IGG: A tool to integrate GeneChips for genetic studies

To facilitate genetic studies using high-throughput genotyping technologies, we have developed an open source tool to integrate genotype data across the Affymetrix and Illumina platforms. It can efficiently integrate a large amount of data from various GeneChips, add genotypes of the HapMap Project into a specific project, flexibly trim and export the integrated data with different formats of popular genetic analysis tools, and highly control the quality of genotype data. Furthermore, this tool has sufficiently simplified its usage through its user-friendly graphic interface and is independent of third-party databases. IGG has successfully been applied to a genome-wide linkage scan in a Charcot-Marie-Tooth disease pedigree by integrating three types of GeneChips and HapMap project genotypes.     

On the firing maps of a general class of forced integrate and fire neurons

Integrate and fire processes are fundamental mechanisms causing excitable and oscillatory behavior. Van der Pol [Philos. Mag. (7) 2 (11) (1926) 978] studied oscillations caused by these processes, which he called 'relaxation oscillations' and pointed out their relevance, not only to engineering, but also to the understanding of biological phenomena [Acta Med. Scand. Suppl. CVIII (108) (1940) 76], like cardiac rhythms and arrhythmias. The complex behavior of externally stimulated integrate and fire oscillators has motivated the study of simplified models whose dynamics are determined by iterations of 'firing circle maps' that can be studied in terms of Poincaré's rotation theory [Chaos 1 (1991) 20; Chaos 1 (1991) 13; SIAM J. Appl. Math. 41 (3) (1981) 503]. In order to apply this theory to understand the responses and bifurcation patterns of forced systems, it is fundamental to determine the regions in parameter space where the different regularity properties (e.g., continuity and injectivity) of the firing maps are satisfied. Methods for carrying out this regularity analysis for linear systems, have been devised and the response of integrate and fire neurons (with linear accumulation) to a cyclic input has been analyzed [SIAM J. Appl. Math. 41 (3) (1981) 503]. In this paper we are concerned with the most general class of forced integrate and fire systems, modelled by one first-order differential equation. Using qualitative analysis we prove theorems on which we base a new method of regularity analysis of the firing map, that, contrasting with methods previously reported in the literature, does not requires analytic knowledge of the solutions of the differential equation and therefore it is also applicable to non-linear integrate and fire systems. To illustrate this new methodology, we apply it to determine the regularity regions of a non-linear example whose firing maps undergo bifurcations that were unknown for the previously studied linear systems.     

Dynamic genetic features of chromosomes revealed by comparison of soybean genetic and sequence-based physical maps

Despite the intensive soybean [Glycine max (L.) Merrill] genome studies, the high chromosome number (20) of the soybean plant relative to many other major crops has hindered the development of a high-resolution genomewide genetic map derived from a single population. Here, we report such a map, which was constructed in an F₁₅ population derived from a cross between G. max and G. soja lines using indel polymorphisms detected via a G. soja genome resequencing. By targeting novel indel markers to marker-poor regions, all marker intervals were reduced to under 6 cM on a genome scale. Comparison of the Williams 82 soybean reference genome sequence and our genetic map indicated that marker orders of 26 regions were discrepant with each other. In addition, our comparison showed seven misplaced and two absent markers in the current Williams 82 assembly and six markers placed on the scaffolds that were not incorporated into the pseudomolecules. Then, we showed that, by determining the missing sequences located at the presumed beginning points of the five major discordant segments, these observed discordant regions are mostly errors in the Williams 82 assembly. Distributions of the recombination rates along the chromosomes were similar to those of other organisms. Genotyping of indel markers and genome resequencing of the two parental lines suggested that some marker-poor chromosomal regions may represent introgression regions, which appear to be prevalent in soybean. Given the even and dense distribution of markers, our genetic map can serve as a bridge between genomics research and breeding programs.     

From E‐MAPs to module maps: dissecting quantitative genetic interactions using physical interactions

Recent technological breakthroughs allow the quantification of hundreds of thousands of genetic interactions (GIs) in Saccharomyces cerevisiae. The interpretation of these data is often difficult, but it can be improved by the joint analysis of GIs along with complementary data types. Here, we describe a novel methodology that integrates genetic and physical interaction data. We use our method to identify a collection of functional modules related to chromosomal biology and to investigate the relations among them. We show how the resulting map of modules provides clues for the elucidation of function both at the level of individual genes and at the level of functional modules.     

Exploratory integration of peanut genetic and physical maps and possible contributions from Arabidopsis

Arachis hypogaea is a widely cultivated crop both as an oilseed and protein source. The genomic analysis of Arachis species hitherto has been limited to the construction of genetic maps; the most comprehensive one contains 370 loci over 2,210 cM in length. However, no attempt has been made to analyze the physical structure of the peanut genome. To investigate the practicality of physical mapping in peanut, we applied a total of 117 oligonucleotide-based probes (overgos) derived from genetically mapped RFLP probes onto peanut BAC filters containing 182,784 peanut large-insert DNA clones in a multiplex experimental design; 91.5% of the overgos identified at least one BAC clone. In order to gain insights into the potential value of Arabidopsis genome sequence for studies in divergent species with complex genomes such as peanut, we employed 576 Arabidopsis-derived overgos selected on the basis of maximum homology to orthologous sequences in other plant taxa to screen the peanut BAC library. A total of 353 (61.3%) overgos detected at least one peanut BAC clone. This experiment represents the first steps toward the creation of a physical map in peanut and illustrates the potential value of leveraging information from distantly related species such as Arabidopsis for both practical applications such as comparative map-based cloning and shedding light on evolutionary relationships. We also evaluated the possible correlation between functional categories of Arabidopsis overgos and their success rates in hybridization to the peanut BAC library.Electronic Supplementary Material Supplementary material is available for this article at