A nearest-neighboring-end algorithm for genetic mapping

MOTIVATION: High-throughput methods are beginning to make possible the genotyping of thousands of loci in thousands of individuals, which could be useful for tightly associating phenotypes to candidate loci. Current mapping algorithms cannot handle so many data without building hierarchies of framework maps. RESULTS: A version of Kruskal's minimum spanning tree algorithm can solve any genetic mapping problem that can be stated as marker deletion from a set of linkage groups. These include backcross, recombinant inbred, haploid and double-cross recombinational populations, in addition to conventional deletion and radiation hybrid populations. The algorithm progressively joins linkage groups at increasing recombination fractions between terminal markers, and attempts to recognize and correct erroneous joins at peaks in recombination fraction. The algorithm is O (mn3) for m individuals and n markers, but the mean run time scales close to mn2. It is amenable to parallel processing and has recovered true map order in simulations of large backcross, recombinant inbred and deletion populations with up to 37,005 markers. Simulations were used to investigate map accuracy in response to population size, allelic dominance, segregation distortion, missing data and random typing errors. It produced accurate maps when marker distribution was sufficiently uniform, although segregation distortion could induce translocated marker orders. The algorithm was also used to map 1003 loci in the F7 ITMI population of bread wheat, Triticum aestivum L. emend Thell., where it shortened an existing standard map by 16%, but it failed to associate blocks of markers properly across gaps within linkage groups. This was because it depends upon the rankings of recombination fractions at individual markers, and is susceptible to sampling error, typing error and joint selection involving the terminal markers of nearly finished linkage groups. Therefore, the current form of the algorithm is useful mainly to improve local marker ordering in linkage groups obtained in other ways. AVAILABILITY: The source code and supplemental data are http://www.iubio.bio.indiana.edu/soft/molbio/qtl/flipper/ CONTACT: ccrane@purdue.edu.     

Rat Chromosome 1: Regional localization of seven genes (Slc9a3, Srd5a1, Esr, Tcp1, Grik5, Tnnt3, Jak2) and anchoring of the genetic linkage map to the cytogenetic map

Seven genes were regionally localized on rat Chromosome (Chr) 1, from 1p11 to 1q42, and two of these genes were also included in a linkage map. This mapping work integrates the genetic linkage map and the cytogenetic map, and allows us to orient the linkage map with respect to the centromere, and to deduce the approximate position of the centromere in the linkage map. These mapping data also indicate that the Slc9a3 gene, encoding the Na⁺/H⁺ exchanger 3, is an unlikely candidate for the blood pressure loci assigned to rat Chr 1. These new localizations expand comparative mapping between rat Chr 1 and mouse or human chromosomes. Received: 21 March 1997 / Accepted: 3 May 1997     

Iterative sequence/secondary structure search for protein homologs: comparison with amino acid sequence alignments and application to fold recognition in genome databases

MOTIVATION: Sequence alignment techniques have been developed into extremely powerful tools for identifying the folding families and function of proteins in newly sequenced genomes. For a sufficiently low sequence identity it is necessary to incorporate additional structural information to positively detect homologous proteins. We have carried out an extensive analysis of the effectiveness of incorporating secondary structure information directly into the alignments for fold recognition and identification of distant protein homologs. A secondary structure similarity matrix based on a database of three-dimensionally aligned proteins was first constructed. An iterative application of dynamic programming was used which incorporates linear combinations of amino acid and secondary structure sequence similarity scores. Initially, only primary sequence information is used. Subsequently contributions from secondary structure are phased in and new homologous proteins are positively identified if their scores are consistent with the predetermined error rate. RESULTS: We used the SCOP40 database, where only PDB sequences that have 40% homology or less are included, to calibrate homology detection by the combined amino acid and secondary structure sequence alignments. Combining predicted secondary structure with sequence information results in a 8-15% increase in homology detection within SCOP40 relative to the pairwise alignments using only amino acid sequence data at an error rate of 0.01 errors per query; a 35% increase is observed when the actual secondary structure sequences are used. Incorporating predicted secondary structure information in the analysis of six small genomes yields an improvement in the homology detection of approximately 20% over SSEARCH pairwise alignments, but no improvement in the total number of homologs detected over PSI-BLAST, at an error rate of 0.01 errors per query. However, because the pairwise alignments based on combinations of amino acid and secondary structure similarity are different from those produced by PSI-BLAST and the error rates can be calibrated, it is possible to combine the results of both searches. An additional 25% relative improvement in the number of genes identified at an error rate of 0.01 is observed when the data is pooled in this way. Similarly for the SCOP40 dataset, PSI-BLAST detected 15% of all possible homologs, whereas the pooled results increased the total number of homologs detected to 19%. These results are compared with recent reports of homology detection using sequence profiling methods. AVAILABILITY: Secondary structure alignment homepage at http://lutece.rutgers.edu/ssas CONTACT: anders@rutchem.rutgers.edu; ronlevy@lutece.rutgers.edu Supplementary Information: Genome sequence/structure alignment results at http://lutece.rutgers.edu/ss_fold_predictions.     

3FunMap: full-sib family functional mapping of dynamic traits

MOTIVATION: Functional mapping that embeds the developmental mechanisms of complex traits shows great power to study the dynamic pattern of genetic effects triggered by individual quantitative trait loci (QTLs). A full-sib family, produced by crossing two heterozygous parents, is characteristic of uncertainties about cross-type at a locus and linkage phase between different loci. Integrating functional mapping into a full-sib family requires a model selection procedure capable of addressing these uncertainties. 3FunMap, written in VC++ 6.0, provides a flexible and extensible platform to perform full-sib functional mapping of dynamic traits. Functions in the package encompass linkage phase determination, marker map construction and the pattern identification of QTL segregation, dynamic tests of QTL effects, permutation tests and numerical simulation. We demonstrate the features of 3FunMap through real data analysis and computer simulation. AVAILABILITY: http://statgen.psu.edu/software.     

Meiotic linkage mapping of 52 genes onto the canine map does not identify significant levels of microrearrangement

In an effort to extend our understanding of the evolutionary relationship between the canine and human genomes, we have developed and positioned 52 new gene-associated polymorphic markers on the canine meiotic linkage map. Canine-specific PCR primers were developed from the consensus of published sequences of several mammalian genomes and were designed to span intronic regions, thus optimizing the probability that a polymorphic site was included. The resulting markers were analyzed on a panel of three-generation canine reference families and the data were incorporated into the current meiotic linkage map. The data were compared with those generated by three chromosome paint studies in an effort to understand the distribution and frequency of microrearrangements within the canine genome. Forty-eight of 52 genes map to a chromosomal region predicted to contain genes from the corresponding region of the human genome according to all published reciprocal chromosome paint studies. Meiotic linkage mapping data for three genes can be used to resolve discrepancies between the published reciprocal chromosome paint studies, and for an additional two genes, meiotic mapping data allow evolutionary breakpoints to be more precisely defined. We conclude that microrearrangements of evolutionarily conserved segments between the canine and human genomes are rare, occurring for less than 0.5% of gene data reported to date. In addition, we have found that the placement of genes on the meiotic linkage map is a useful mechanism for resolving discrepancies between existing data sets. Received: 7 February 2001 / Accepted: 9 May 2001     

Mendel, a database of nomenclature for sequenced plant genes

The Mendel database contains names for plant-wide families of sequenced plant genes. The names have either been approved by the Commission on Plant Gene Nomenclature (CPGN), an organization of the International Society for Plant Molecular Biology (ISPMB), or are identified as provisional or temporary names. Mendel also identifies the corresponding genes in individual species of plants. Mendel can be searched through the mirror sites at Cornell (http://genome. cornell.edu/cgi-bin/WebAce/webace?db=mendel) and Stanford (http://genome-www.stanford.edu/Mendel/). In addition, parts of Mendel can be downloaded from the CPGN Web site (http://mbclserver. rutgers.edu/CPGN/).     

FunMap: functional mapping of complex traits

SUMMARY: FunMap is a Web-based user interface designed to map quantitative trait loci (QTL) affecting function-valued traits or infinite-dimensional traits in well-structured pedigrees or natural populations. User input includes three files: longitudinal trait data, marker genotypes and/or a linkage map. This software allows for a systematic genome-wide scan and significance test of QTL throughout the map. The dynamic change of QTL effects during the time course of growth is automatically drawn, from which specific biological hypotheses regarding the genetic control mechanisms of growth and development can be tested. AVAILABILITY: http://web.biostat.ufl.edu/~cma/genetics/software.html     

A comparative map of the porcine and human genomes demonstrates ZOO-FISH and gene mapping-based chromosomal homologies

ZOO-FISH with chromosome-specific DNA libraries (CSLs) from individual flow-sorted human chromosomes was applied on porcine metaphase chromosomes to establish segment homology between the pig and human karyotypes. Forty-seven porcine chromosomal segments corresponding to all human chromosomes except the Y were delineated, resulting in a nearly complete coverage of the porcine karyotype. The syntenic segments detected were further confirmed by the gene mapping information available in the two species. A map demarcating physical boundaries of human homologies on individual pig chromosomes is complemented with a detail survey of the physical and genetic linkage mapping data in the two species. The resultant map, thus, provides a comprehensive and updated comparative status of the human and porcine genomes. Received: 9 September 1995 / Accepted: 4 December 1995     

Kernel-based association test

Association mapping (i.e., linkage disequilibrium mapping) is a powerful tool for positional cloning of disease genes. We propose a kernel-based association test (KBAT), which is a composite function of "P-values of single-locus association tests" and "kernel weights related to intermarker distances and/or linkage disequilibria." The KBAT is a general form of some current test statistics. This method can be applied to the study of candidate genes and can scan each chromosome using a moving average procedure. We evaluated the performance of the KBAT through simulation studies that considered evolutionary parameters, disease models, sample sizes, kernel functions, test statistics, window attributes, empirical P-value estimations, and genetic/physical maps. The results showed that the KBAT had a well-controlled false positive rate and high power compared to existing methods. In addition, the KBAT was also applied to analyze a genomewide data set from the Collaborative Study on the Genetics of Alcoholism. Important genes associated with alcoholism dependence were identified. In summary, the merits of the KBAT are multifold: the KBAT is robust against the inclusion of nuisance markers, is invariant to the map scale, and accommodates different types of genomic data, study designs, and study purposes. The proposed methods are packaged in the user-friendly software, KBAT, available at http://www.stat.sinica.edu.tw/hsinchou/genetics/association/KBAT.htm.     

A gene map of the rat derived from linkage analysis and related regions in the mouse and human genomes

We report the localization by linkage analysis in the rat genome of 148 new markers derived from 128 distinct known gene sequences, ESTs, and anonymous sequences selected in GenBank database on the basis of the presence of a repeated element. The composite linkage map of the rat contributed by our group integrates mapping information on a total of 370 different known genes, ESTs, and anonymous mouse or human sequences, and provides a valuable tool for comparative genome analysis. 206 and 254 homologous loci were identified in the mouse and human genomes respectively. Our linkage map, which combines both anonymous markers and gene markers, should facilitate the advancement of genetic studies for a wide variety of rat models characterized for complete phenotypes. The comparative genome mapping should define genetic regions in human likely to be homologous to susceptibility loci identified in rat and provide useful information for the identification of new potential candidates for genetic disorders. Received: 2 January 1999 / Accepted: 7 March 1999     

Haploview: analysis and visualization of LD and haplotype maps

SUMMARY: Research over the last few years has revealed significant haplotype structure in the human genome. The characterization of these patterns, particularly in the context of medical genetic association studies, is becoming a routine research activity. Haploview is a software package that provides computation of linkage disequilibrium statistics and population haplotype patterns from primary genotype data in a visually appealing and interactive interface. AVAILABILITY: http://www.broad.mit.edu/mpg/haploview/ CONTACT: jcbarret@broad.mit.edu     

The bovine B and C blood group systems are not likely to be the orthologues of human RH: an interesting twist in the comparative map

We tested the hypothesis that either the bovine B or C blood group system is the orthologue of human RH. A comparative linkage mapping strategy was applied, using blood typing and restriction fragment length polymorphism (RFLP) analysis of four loci linked to RH on HSA1; PGD, FGR, ALPL and FUCA1. Four sires with a total of 255 half-sib offspring were used for the linkage analysis. Strong support for linkage between ALPL, FUCA1 and FGR was obtained for all sire families (lod scores >11 for all pairwise comparisons). This new linkage group was assigned to bovine synteny group U17 based on previous somatic cell mapping of the FGR locus. The most favoured order is ALPL—FUCA1—FGR (2·18:1), with ALPL and FGR 5·4 cm and 6·3 cm , respectively, from FUCA1. The B and C blood group systems and PGD were genetically independent of each other and all other markers, indicating that neither B nor C is likely to be the bovine orthologue of human RH. However, given available comparative mapping data, there is some chance that the bovine orthologue of RH is on bovine synteny group U6. Although gene order appears to be conserved with humans, the differences in recombination rates between these three loci in cattle, humans and mice strongly suggest that it is not possible to use human map distances to predict map distances in cattle, making it imperative that bovine gene mappers continue to emphasize adding type I markers to the bovine linkage map.     

OmicsViz: Cytoscape plug-in for visualizing omics data across species

MOTIVATION: OmicsViz is a Cytoscape plug-in for mapping and visualizing large-scale omics datasets across species, including those with many-to-many mappings between homologs. This allows users to map their data onto pathways of related model organisms. Mapping schemas across species or different experimental protocols allow users to comparatively analyze the omics data. The data can also be visualized in parallel-coordinate plots. AVAILABILITY: The latest version of OmicsViz with documentation, tutorial and jar files can be downloaded from http://metnet.vrac.iastate.edu/MetNet_fcmodeler.htm     

A primary linkage map of the human chromosome 11q22-23 region

We have constructed a genetic map of the human chromosomal region 11q22-23 by multipoint linkage analysis of 13 DNA polymorphisms that we have condensed into eight loci. An analysis for linkage disequilibrium between tightly linked probe/enzyme systems allows us to make specific recommendations for future DNA typing at these loci. The resulting sex-averaged multipoint map spans approximately 80 cM and differs considerably from previously reported genetic maps of this region. Our mathematically derived "most likely order" of the markers is compatible with physical mapping data using somatic cell hybrids. The known localizations of at least 14 functional genes and several disease loci to 11q22-23, including ataxia telangiectasia, make the mapping of this region especially relevant to studies of disease pathogenesis.     

A radiation hybrid map of chicken Chromosome 4

The mapping resolution of the physical map for chicken Chromosome 4 (GGA4) was improved by a combination of radiation hybrid (RH) mapping and bacterial artificial chromosome (BAC) mapping. The ChickRH6 hybrid panel was used to construct an RH map of GGA4. Eleven microsatellites known to be located on GGA4 were included as anchors to the genetic linkage map for this chromosome. Based on the known conserved synteny between GGA4 and human Chromosomes 4 and X, sequences were identified for the orthologous chicken genes from these human chromosomes by BLAST analysis. These sequences were subsequently used for the development of STS markers to be typed on the RH panel. Using a logarithm of the odds (LOD) threshold of 5.0, nine linkage groups could be constructed which were aligned with the genetic linkage map of this chromosome. The resulting RH map consisted of the 11 microsatellite markers and 50 genes. To further increase the number of genes on the map and to provide additional anchor points for the physical BAC map of this chromosome, BAC clones were identified for 22 microsatellites and 99 genes. The combined RH and BAC mapping approach resulted in the mapping of 61 genes on GGA4 increasing the resolution of the chicken–human comparative map for this chromosome. This enhanced comparative mapping resolution enabled the identification of multiple rearrangements between GGA4 and human Chromosomes 4q and Xp.     

Physical and linkage mapping of mammary-derived expressed sequence tags in cattle

This study describes the physical and linkage mapping of 42 gene-associated markers developed from mammary gland-derived expressed sequence tags to the cattle genome. Of the markers, 25 were placed on the USDA reference linkage map and 37 were positioned on the Roslin 3000-rad radiation hybrid (RH) map, with 20 assignments shared between the maps. Although no novel regions of conserved synteny between the cattle and the human genomes were identified, the coverage was extended for bovine chromosomes 3, 7, 15, and 29 compared with previously published comparative maps between human and bovine genomes. Overall, these data improve the resolution of the human-bovine comparative maps and will assist future efforts to integrate bovine RH and linkage map data.     

Comparison off genetic distance and order off DNA markers in five populations of rice.

A group of about 300 evenly distributed DNA markers from a high density RFLP linkage map of rice constructed using an F2 population derived from a japonica variety, Nipponbare, and an indica variety, Kasalath, were used to evaluate gene order and genetic distance in four other rice mapping populations. The purpose of this study was to determine the degree to which information gained from the high density linkage map could be applied to other mapping populations, particularly with regard to its utility in bridging quantitative traits and molecular and physical mapping information. The mapping populations consisted of two F2 populations derived from Dao Ren Qiao/Fl-1084 and Kinandangputi/Fl-1007, recombinant inbred lines from Asominori/IR24, and a backcross population from Sasanishiki/Habataki//Sasanishiki. All DNA markers commonly mapped in the four populations showed the same linkage groups as in the Nipponbare/Kasalath linkage map with conserved linkage order. The genetic distance between markers among the different populations did not vary to a significant level in any of the 12 chromosomes. The differences in some markers could be attributed to the size of the population used in the construction of the linkage maps. Furthermore, the conservation of linkage order found in the distal region of chromosomes 11 and 12 was also confirmed in the RFLP maps based on the four populations of rice. These suggest that any major genetic information from the Nipponbare/Kasalath map can be expected to be approximately the same in other crosses or populations. This high density RFLP linkage map, which is being utilized in constructing a physical map of rice, can be very useful in interpreting genome structure with great accuracy in other populations. Key words : linkage map, japonica, indica, gene order, genetic distance.     

Flexible network reconstruction from relational databases with Cytoscape and CytoSQL

Background

Bisulfite sequencing using next generation sequencers yields genome-wide measurements of DNA methylation at single nucleotide resolution. Traditional aligners are not designed for mapping bisulfite-treated reads, where the unmethylated Cs are converted to Ts. We have developed BS Seeker, an approach that converts the genome to a three-letter alphabet and uses Bowtie to align bisulfite-treated reads to a reference genome. It uses sequence tags to reduce mapping ambiguity. Post-processing of the alignments removes non-unique and low-quality mappings.

Results

We tested our aligner on synthetic data, a bisulfite-converted Arabidopsis library, and human libraries generated from two different experimental protocols. We evaluated the performance of our approach and compared it to other bisulfite aligners. The results demonstrate that among the aligners tested, BS Seeker is more versatile and faster. When mapping to the human genome, BS Seeker generates alignments significantly faster than RMAP and BSMAP. Furthermore, BS Seeker is the only alignment tool that can explicitly account for tags which are generated by certain library construction protocols.

Conclusions

BS Seeker provides fast and accurate mapping of bisulfite-converted reads. It can work with BS reads generated from the two different experimental protocols, and is able to efficiently map reads to large mammalian genomes. The Python program is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker/BS_Seeker.html.