Experimental design and error detection for polyploid radiation hybrid mapping

Evaluation of genotyping concordance for commercial bovine SNP arrays using quality-assurance samples

SNP arrays are widely used in genetic research and agricultural genomics applications, and the quality of SNP genotyping data is of paramount importance. In the present study, SNP genotyping concordance and discordance were evaluated for commercial bovine SNP arrays based on two types of quality assurance (QA) samples provided by Neogen GeneSeek. The genotyping discordance rates (GDRs) between chips were on average between 0.06% and 0.37% based on the QA type I data and between 0.05% and 0.15% based on the QA type II data. The average genotyping error rate (GER) pertaining to single SNP chips, based on the QA type II data, varied between 0.02% and 0.08% per SNP and between 0.01% and 0.06% per sample. These results indicate that genotyping concordance rate was high (i.e. from 99.63% to 99.99%). Nevertheless, mitochondrial and Y chromosome SNPs had considerably elevated GDRs and GERs compared to the SNPs on the 29 autosomes and X chromosome. The majority of genotyping errors resulted from single allotyping errors, which also included the opposite instances for allele 'dropout' (i.e. from AB to AA or BB). Simultaneous allotyping errors on both alleles (e.g. mistaking AA for BB or vice versa) were relatively rare. Finally, a list of SNPs with a GER greater than 1% is provided. Interpretation of association effects of these SNPs, for example in genome-wide association studies, needs to be taken with caution. The genotyping concordance information needs to be considered in the optimal design of future bovine SNP arrays.

An update of the goat genome assembly using dense radiation hybrid maps allows detailed analysis of evolutionary rearrangements in Bovidae

Background

The domestic goat (Capra hircus), an important livestock species, belongs to a clade of Ruminantia, Bovidae, together with cattle, buffalo and sheep. The history of genome evolution and chromosomal rearrangements on a small scale in ruminants remain speculative. Recently completed goat genome sequence was released but is still in a draft stage. The draft sequence used a variety of assembly packages, as well as a radiation hybrid (RH) map of chromosome 1 as part of its validation.

Results

Using an improved RH mapping pipeline, whole-genome dense maps of 45,953 SNP markers were constructed with statistical confidence measures and the saturated maps provided a fine map resolution of approximate 65 kb. Linking RH maps to the goat sequences showed that the assemblies of scaffolds/super-scaffolds were globally accurate. However, we observed certain flaws linked to the process of anchoring chromosome using conserved synteny with cattle. Chromosome assignments, long-range order, and orientation of the scaffolds were reassessed in an updated genome sequence version. We also present new results exploiting the updated goat genome sequence to understand genomic rearrangements and chromosome evolution between mammals during species radiations. The sequence architecture of rearrangement sites between the goat and cattle genomes presented abundant segmental duplication on regions of goat chromosome 9 and 14, as well as new insertions in homologous cattle genome regions. This complex interplay between duplicated sequences and Robertsonian translocations highlights the rearrangement mechanism of centromeric nonallelic homologous recombination (NAHR) in mammals. We observed that species-specific shifts in ANKRD26 gene duplication are coincident with breakpoint reuse in divergent lineages and this gene family may play a role in chromosome stabilization in chromosome evolution.

Conclusions

We generated dense maps of the complete whole goat genome. The chromosomal maps allowed us to anchor and orientate assembled genome scaffolds along the chromosomes, annotate chromosome rearrangements and thereby get a better understanding of the genome evolution of ruminants and other mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-625) contains supplementary material, which is available to authorized users.

Statistical confidence measures for genome maps: application to the validation of genome assemblies

MOTIVATION

Genome maps are imperative to address the genetic basis of the biology of an organism. While a growing number of genomes are being sequenced providing the ultimate genome maps-this being done at an even faster pace now using new generation sequencers-the process of constructing intermediate maps to build and validate a genome assembly remains an important component for producing complete genome sequences. However, current mapping approach lack statistical confidence measures necessary to identify precisely relevant inconsistencies between a genome map and an assembly.

RESULTS

We propose new methods to derive statistical measures of confidence on genome maps using a comparative model for radiation hybrid data. We describe algorithms allowing to (i) sample from a distribution of maps and (ii) exploit this distribution to construct robust maps. We provide an example of application of these methods on a dog dataset that demonstrates the interest of our approach.

AVAILABILITY

Methods are implemented in two freely available softwares: Carthagene (http://www.inra.fr/mia/T/CarthaGene/) and a companion software (metamap, available at: http://snp.toulouse.inra.fr/~servin/index.cgi/Metamap).

Enhancing radiation hybrid mapping through whole genome amplification

Radiation hybrid (RH) mapping is limited by the inherent genomic instability of RH clones entailing both, limited DNA sample amounts and genomic heterogeneity of the clones. Here the instability of RH clones is quantified and the suitability of the multiple strand displacement whole genome amplification method (WGA) for radiation hybrid mapping is assessed. To quantify the instability of RH clones, eleven clones of a 10,000(Rad) rhesus macaque radiation hybrid panel were passaged ten times and analyzed by interspersed repeat sequence specific quantitative PCR and by genotyping of 46 macaque chromosome 5 STS markers. The quantitative PCR data indicate an average loss of 55% of the donor DNA over 10 passages. Over the same period, a dropout of 46.2% of the STS markers was observed. These data indicate a genome wide half-life of the donor DNA of 8.7 passages and of 10.6 passages for the chromosome 5 markers. The genotyping data of the genomic RH DNA were compared to three sets of WGA experiments: 1) single wgaDNA amplifications, 2) six WGA replicates, and 3) re-amplification of wga DNA. The assays demonstrated concordance rates of 97.6%, 98% and 99.3%, respectively, and indicated the marker specificity of some repeated WGA dropouts. The study confirms that WGA is suitable for RH mapping studies should enable the accurate analysis of almost an infinite numbers of markers. WGA will allow the analysis of earliest RH clone passages, thus limiting their heterogeneity and RH mapping artifacts.

A high-resolution radiation hybrid map of rhesus macaque chromosome 5 identifies rearrangements in the genome assembly

A 10,000-rad radiation hybrid (RH) cell panel of the rhesus macaque was generated to construct a comprehensive RH map of chromosome 5. The map represents 218 markers typed in 185 RH clones. The 4846-cR map has an average marker spacing of 798 kb. Alignments of the RH map to macaque and human genome sequences confirm a large inversion and reveal a previously unreported telomeric inversion. The macaque genome sequence indicates small translocations from the ancestral homolog of macaque chromosome 5 to macaque chromosomes 1 and 6. The RH map suggests that these are probably assembly artifacts. Unlike the genome sequence, the RH mapping data indicate the conservation of synteny between macaque chromosome 5 and human chromosome 4. This study shows that the 10,000-rad panel is appropriate for the generation of a high-resolution whole-genome RH map suitable for the verification of the rhesus genome assembly.

High resolution radiation hybrid maps of bovine chromosomes 19 and 29: comparison with the bovine genome sequence assembly

Background

High resolution radiation hybrid (RH) maps can facilitate genome sequence assembly by correctly ordering genes and genetic markers along chromosomes. The objective of the present study was to generate high resolution RH maps of bovine chromosomes 19 (BTA19) and 29 (BTA29), and compare them with the current 7.1X bovine genome sequence assembly (bovine build 3.1). We have chosen BTA19 and 29 as candidate chromosomes for mapping, since many Quantitative Trait Loci (QTL) for the traits of carcass merit and residual feed intake have been identified on these chromosomes.

Results

We have constructed high resolution maps of BTA19 and BTA29 consisting of 555 and 253 Single Nucleotide Polymorphism (SNP) markers respectively using a 12,000 rad whole genome RH panel. With these markers, the RH map of BTA19 and BTA29 extended to 4591.4 cR and 2884.1 cR in length respectively. When aligned with the current bovine build 3.1, the order of markers on the RH map for BTA19 and 29 showed inconsistencies with respect to the genome assembly. Maps of both the chromosomes show that there is a significant internal rearrangement of the markers involving displacement, inversion and flips within the scaffolds with some scaffolds being misplaced in the genome assembly. We also constructed cattle-human comparative maps of these chromosomes which showed an overall agreement with the comparative maps published previously. However, minor discrepancies in the orientation of few homologous synteny blocks were observed.

Conclusion

The high resolution maps of BTA19 (average 1 locus/139 kb) and BTA29 (average 1 locus/208 kb) presented in this study suggest that by the incorporation of RH mapping information, the current bovine genome sequence assembly can be significantly improved. Furthermore, these maps can serve as a potential resource for fine mapping QTL and identification of causative mutations underlying QTL for economically important traits.

Development of a deer mouse whole-genome radiation hybrid panel and comparative mapping of Mus chromosome 11 loci

A 5000-rad whole-genome radiation hybrid cell panel (BW5000) was developed for mapping the deer mouse (Peromyscus maniculatus bairdii) genome. The panel consists of 103 cell lines and has an estimated marker retention frequency of 63.9% (range, 28%-88%) based on PCR typing of 30 Type I (coding gene) and 25 Type II (microsatellite) markers. Using the composite Mus map, Type I markers were selected from six Mus chromosomes, 22 of which are on Mus Chr 11. Fifteen of the Mus Chr 11 markers were simultaneously mapped on an interspecific (P. maniculatus x P. polionotus) backcross panel to test the utility of the radiation hybrid panel, create a framework map, and help establish gene order. The radiation hybrids have effectively detected linkage in the deer mouse genome between markers as far apart as 6.7 cM and resolved markers that are, in the Mus genome, as close as 0.2 Mb. Combined results from both panels have indicated a high degree of gene order conservation of the telomeric 64 cM of Mus Chr 11 in the deer mouse genome. The remaining centromeric portion also shows gene order conservation with the deer mouse but as a separate linkage group. This indicates a translocation of that portion of Mus Chr 11 in P. maniculatus and is consistent with rearrangement breakpoints observed between Mus and other mammalian genomes, including rat and human. Furthermore, this separate linkage group is likely to reside in a chromosomal region of inversion polymorphism between P. maniculatus and P. polionotus.

Haplotype sharing transmission/disequilibrium tests that allow for genotyping errors

The present study introduces new Haplotype Sharing Transmission/Disequilibrium Tests (HS-TDTs) that allow for random genotyping errors. We evaluate the type I error rate and power of the new proposed tests under a variety of scenarios and perform a power comparison among the proposed tests, the HS-TDT and the single-marker TDT. The results indicate that the HS-TDT shows a significant increase in type I error when applied to data in which either Mendelian inconsistent trios are removed or Mendelian inconsistent markers are treated as missing genotypes, and the magnitude of the type I error increases both with an increase in sample size and with an increase in genotyping error rate. The results also show that a simple strategy, that is, merging each rare haplotype to a most similar common haplotype, can control the type I error inflation for a wide range of genotyping error rates, and after merging rare haplotypes, the power of the test is very similar to that without merging the rare haplotypes. Therefore, we conclude that a simple strategy may make the HS-TDT robust to genotyping errors. Our simulation results also show that this strategy may also be applicable to other haplotype-based TDTs.

Expression of G-protein alpha subunit genes in the vestibular periphery of Rattus norvegicus and their chromosomal mapping

OBJECTIVE

Heterotrimeric G-proteins play an important role in mediating signals transduced across the cell membrane by membrane-bound receptors. The precise role of G-proteins and their coupled receptors in the physiology of the vestibular neuroepithelium is not well understood. The purpose of this study was to better define the role of these proteins by examining their expression in the rat vestibular periphery and characterizing their chromosomal location.

MATERIAL AND METHODS

To characterize G-protein alpha subunit gene expression in the target tissue of interest, we performed polymerase chain reaction (PCR) using degenerate G-protein primers corresponding to conserved regions in the G-protein alpha subunit coding sequence on a normalized rat vestibular cDNA library. PCR amplicons were cloned and 50 clones were randomly selected and sequenced. Radiation hybrid (RH) mapping was used to determine the chromosomal location of G alpha(olf) and two previously identified G-protein alpha subunits--G alpha(i2) and G alpha(i2(vest))--in the rat genome.

RESULTS

The following G-protein alpha subunits were identified in the normalized cDNA library: G alpha(olf), G alpha(s), G alpha(o) and G alpha(s2). G alpha(olf) maps to chromosome 18 between markers D18Mit17b and D18Mgh2. G alpha(i2) maps to chromosome 8 between markers D8Rat65 and D8Mgh2. G alpha(i2(vest)) maps to chromosome 1 between markers D1Rat132 and D1Rat202. These chromosomal locations in the rat genome are syntenic to chromosomal regions in which the homologous G-protein alpha subunit genes have been localized in the human and mouse genomes, further validating RH mapping as an effective and accurate tool. We were unable to RH map the location of G alpha(o) due to its extensive homology with the hamster gene.

CONCLUSION

The characterization of G-protein alpha subunit gene expression in the vestibular periphery and the chromosomal localization of these genes in the rat revealed that a diverse group of these second messengers are expressed.

The first-generation whole-genome radiation hybrid map in the horse identifies conserved segments in human and mouse genomes

A first-generation radiation hybrid (RH) map of the equine (Equus caballus) genome was assembled using 92 horse x hamster hybrid cell lines and 730 equine markers. The map is the first comprehensive framework map of the horse that (1) incorporates type I as well as type II markers, (2) integrates synteny, cytogenetic, and meiotic maps into a consensus map, and (3) provides the most detailed genome-wide information to date on the organization and comparative status of the equine genome. The 730 loci (258 type I and 472 type II) included in the final map are clustered in 101 RH groups distributed over all equine autosomes and the X chromosome. The overall marker retention frequency in the panel is approximately 21%, and the possibility of adding any new marker to the map is approximately 90%. On average, the mapped markers are distributed every 19 cR (4 Mb) of the equine genome--a significant improvement in resolution over previous maps. With 69 new FISH assignments, a total of 253 cytogenetically mapped loci physically anchor the RH map to various chromosomal segments. Synteny assignments of 39 gene loci complemented the RH mapping of 27 genes. The results added 12 new loci to the horse gene map. Lastly, comparison of the assembly of 447 equine genes (256 linearly ordered RH-mapped and additional 191 FISH-mapped) with the location of draft sequences of their human and mouse orthologs provides the most extensive horse-human and horse-mouse comparative map to date. We expect that the foundation established through this map will significantly facilitate rapid targeted expansion of the horse gene map and consequently, mapping and positional cloning of genes governing traits significant to the equine industry.

Radiation hybrid mapping of five muscarinic acetylcholine receptor subtype genes in Rattus norvegicus

Acetylcholine is the main neurotransmitter of the vestibular efferent system and a wide variety of muscarinic and nicotinic acetylcholine receptors are expressed in the vestibular periphery. The role of these receptors and in particular the role of muscarinic acetylcholine receptors in the physiology of the vestibular neuroepithelium is not understood. Congenic and consomic rats are a convenient way to investigate the involvement of candidate genes in the manifestation of defined traits. To use congenic or consomic rats to elucidate the roles of these receptors in vestibular physiology or pathology the chromosomal location of the genes encoding these receptors has to be determined. Using radiation hybrid (RH) mapping and a rat RH map server (www.rgd.mcw.edu/RHMAP SERVER/), we determined the chromosomal locations of the muscarinic acetylcholine receptor genes in the rat (Rattus norvegicus). The m1-m5 muscarinic subtypes mapped to the following chromosomes: Chrm1, chromosome 1; Chrm2, chromosome 4; Chrm3, chromosome 17; Chrm4, chromosome 3; and Chrm5, chromosome 3. With the chromosomal location for each of these muscarinic subtypes known, it is now possible to develop congenic and consomic strains of rats that can be used to study the functions of each of these subtypes.

Detection and integration of genotyping errors in statistical genetics

Detection of genotyping errors and integration of such errors in statistical analysis are relatively neglected topics, given their importance in gene mapping. A few inopportunely placed errors, if ignored, can tremendously affect evidence for linkage. The present study takes a fresh look at the calculation of pedigree likelihoods in the presence of genotyping error. To accommodate genotyping error, we present extensions to the Lander-Green-Kruglyak deterministic algorithm for small pedigrees and to the Markov-chain Monte Carlo stochastic algorithm for large pedigrees. These extensions can accommodate a variety of error models and refrain from simplifying assumptions, such as allowing, at most, one error per pedigree. In principle, almost any statistical genetic analysis can be performed taking errors into account, without actually correcting or deleting suspect genotypes. Three examples illustrate the possibilities. These examples make use of the full pedigree data, multiple linked markers, and a prior error model. The first example is the estimation of genotyping error rates from pedigree data. The second-and currently most useful-example is the computation of posterior mistyping probabilities. These probabilities cover both Mendelian-consistent and Mendelian-inconsistent errors. The third example is the selection of the true pedigree structure connecting a group of people from among several competing pedigree structures. Paternity testing and twin zygosity testing are typical applications.

A transmission/disequilibrium test that allows for genotyping errors in the analysis of single-nucleotide polymorphism data

The present study assesses the effects of genotyping errors on the type I error rate of a particular transmission/disequilibrium test (TDT(std)), which assumes that data are errorless, and introduces a new transmission/disequilibrium test (TDT(ae)) that allows for random genotyping errors. We evaluate the type I error rate and power of the TDT(ae) under a variety of simulations and perform a power comparison between the TDT(std) and the TDT(ae), for errorless data. Both the TDT(std) and the TDT(ae) statistics are computed as two times a log-likelihood difference, and both are asymptotically distributed as chi(2) with 1 df. Genotype data for trios are simulated under a null hypothesis and under an alternative (power) hypothesis. For each simulation, errors are introduced randomly via a computer algorithm with different probabilities (called "allelic error rates"). The TDT(std) statistic is computed on all trios that show Mendelian consistency, whereas the TDT(ae) statistic is computed on all trios. The results indicate that TDT(std) shows a significant increase in type I error when applied to data in which inconsistent trios are removed. This type I error increases both with an increase in sample size and with an increase in the allelic error rates. TDT(ae) always maintains correct type I error rates for the simulations considered. Factors affecting the power of the TDT(ae) are discussed. Finally, the power of TDT(std) is at least that of TDT(ae) for simulations with errorless data. Because data are rarely error free, we recommend that researchers use methods, such as the TDT(ae), that allow for errors in genotype data.

A radiation hybrid map of the cat genome: implications for comparative mapping

Ordered gene maps of mammalian species are becoming increasingly valued in assigning gene variants to function in human and animal models, as well as recapitulating the natural history of genome organization. To extend this power to the domestic cat, a radiation hybrid (RH) map of the cat was constructed integrating 424 Type I-coding genes with 176 microsatellite markers, providing coverage over all 20 feline chromosomes. Alignment of parallel RH maps of human and cat reveal 100 conserved segments ordered (CSOs) between the species, nearly three times the number observed with reciprocal chromosome painting analyses. The observed number is equivalent to theoretical predictions of the number of conserved segments to be found between cat and human, implying that 300-400 Type I gene markers is sufficient to reveal nearly all conserved segments for species that exhibit the most frequently observed "slow" rate of genome reorganization. The cat-human RH map comparisons provide a new genomic tool for comparative gene mapping in the cat and related Felidae, and provide confirmation that the cat genome organization is remarkably conserved compared with human. These data demonstrate that ordered RH-based gene maps provide the most precise assessment of comparing genomes, short of contig construction or full-sequence determination.

RHO--radiation hybrid ordering

Radiation hybrid (RH) mapping is a somatic cell technique that is used for ordering markers along a chromosome and estimating the physical distances between them. With the advent of this mapping technique, analyzing the experimental data is becoming a challenging and demanding computational task. In this paper we present the software package RHO (radiation hybrid ordering). The package implements a number of heuristics that attempt to order genomic markers along a chromosome, given as input the results of an RH experiment. The heuristics are based on reducing an appropriate optimization problem to the traveling salesman problem (TSP). The reduced optimization problem is either the nonparametric obligate chromosome breaks (OCBs) or the parametric maximum likelihood estimation (MLE). We tested our package on both simulated and publicly available RH data. For synthetic RH data, the reconstructed markers' permutation is very close to the original permutation, even with fairly high error rates. For real data we used the framework markers' data from the Whitehead Institute maps. For most of the chromosomes (18 out of 23), there is a perfect agreement or nearly perfect agreement (reversal of chromosome arm or arms) between our maps and the Whitehead framework maps. For the remaining five chromosomes, our maps improve on the Whitehead framework maps with respect to both optimization criteria, having higher likelihood and fewer breakpoints. For three chromosomes, the results differ significantly (lod score >1.75), with chromosome 2 having the largest improvement (lod score 3.776).

A radiation hybrid map of bovine chromosome one

Radiation hybrid (RH) mapping has proven to be an extremely powerful approach to constructing high density maps of human chromosomes and is experiencing increased use in other animals, including cattle. A 5000 rad bovine whole-genome radiation hybrid panel was recently constructed in order to integrate existing cattle linkage maps with evolutionarily conserved genes and provide high resolution comparative maps relative to humans and mice. We utilized this panel to construct a 19 marker framework map of bovine chromosome 1 (BTA1), which included 8 Type I loci and 11 Type II loci ordered with at least 1000:1 odds. A 35 marker comprehensive map including 15 Type I loci and 20 Type II loci was also produced. Of the 15 Type I loci ordered on the comprehensive map, three are ordered on HSA3 and five are ordered in three blocks on HSA21 on the human cytogenetic maps.

Fine-mapping of quantitative trait loci by identity by descent in outbred populations: application to milk production in dairy cattle

We previously mapped a quantitative trait locus (QTL) affecting milk production to bovine chromosome 14. To refine the map position of this QTL, we have increased the density of the genetic map of BTA14q11-16 by addition of nine microsatellites and three single nucleotide polymorphisms. Fine-mapping of the QTL was accomplished by a two-tiered approach. In the first phase, we identified seven sires heterozygous "Qq" for the QTL by marker-assisted segregation analysis in a Holstein-Friesian pedigree comprising 1,158 individuals. In a second phase, we genotyped the seven selected sires for the newly developed high-density marker map and searched for a shared haplotype flanking an hypothetical, identical-by-descent QTL allele with large substitution effect. The seven chromosomes increasing milk fat percentage were indeed shown to carry a common chromosome segment with an estimated size of 5 cM predicted to contain the studied QTL. The same haplotype was shown to be associated with increased fat percentage in the general population as well, providing additional support in favor of the location of the QTL within the corresponding interval.

A High-Density Integrated Genetic Linkage and Radiation Hybrid Map of the Laboratory Rat

The laboratory rat (Rattus norvegicus) is a key animal model for biomedical research. However, the genetic infrastructure required for connecting phenotype and genotype in the rat is currently incomplete. Here, we report the construction and integration of two genomic maps: a dense genetic linkage map of the rat and the first radiation hybrid (RH) map of the rat. The genetic map was constructed in two F2 intercrosses (SHRSP × BN and FHH × ACI), containing a total of 4736 simple sequence length polymorphism (SSLP) markers. Allele sizes for 4328 of the genetic markers were characterized in 48 of the most commonly used inbred strains. The RH map is a lod ≥ 3 framework map, including 983 SSLPs, thereby allowing integration with markers on various genetic maps and with markers mapped on the RH panel. Together, the maps provide an integrated reference to >3000 genes and ESTs and >8500 genetic markers (5211 of our SSLPs and >3500 SSLPs developed by other groups). [Bihoreau et al. (1997); James and Tanigami, RHdb (http://www.ebi.ac.uk/RHdb/index.html); Wilder (http://www.nih.gov/niams/scientific/ratgbase); Serikawa et al. (1992); RATMAP server (http://ratmap.gen.gu.se)] RH maps (v. 2.0) have been posted on our web sites at http://goliath.ifrc.mcw.edu/LGR/index.htmlor http://curatools.curagen.com/ratmap. Both web sites provide an RH mapping server where investigators can localize their own RH vectors relative to this map. The raw data have been deposited in the RHdb database. Taken together, these maps provide the basic tools for rat genomics. The RH map provides the means to rapidly localize genetic markers, genes, and ESTs within the rat genome. These maps provide the basic tools for rat genomics. They will facilitate studies of multifactorial disease and functional genomics, allow construction of physical maps, and provide a scaffold for both directed and large-scale sequencing efforts and comparative genomics in this important experimental organism.

Tissue restricted expression of two human Frzbs in preadipocytes and pancreas

Frzb is a newly discovered family of secreted glycoproteins that function to modulate signaling activity of Wnt. Frzb proteins share sequence homology with the extracellular domain of the Wnt receptor (frizzled) and are capable of binding to Wnt. Thus, Frzb functions to antagonize Wnt activity by sequestering Wnt and preventing its binding to the frizzled receptor. Since the initial identification of bovine and human Frzb, several related members of this family have been isolated from rodent and human. In this paper, we describe the cloning and expression of two human frzb homologues termed hFRP-1b and hFRP-2. These human FRPs share significant homology to mouse sFRP-1 and sFRP-2 (55 and 98% identity at amino acid level, respectively). Northern blot experiments revealed that these Frzb homologues have highly restricted tissue distribution. hFRP-1b is exclusively expressed in pancreatic tissue while high levels of hFRP-2 were found in adipose tissue. In addition, low levels of hFRP-2 were also observed in other tissues including heart, pancreas and muscle. Remarkably, FRP-2 is predominantly expressed in un-differentiated preadipocytes in both rodent and man. The expression of FRP-2 is also significantly reduced in fat pads from obese mice. Taken together, these data indicate that distinctive members of the Frzb family exhibit different expression patterns in vivo, suggesting their ability to modulate diverse aspects of Wnt signaling. The expression and dysregulation of sFRP-2 in fat and obesity also suggest a potential roles on the Wnt signaling pathway in the pathology of obesity and related metabolic diseases. Molecular cloning and expression of these Frzbs will allow detailed molecular and biochemical analysis of Wnt-Frzb interaction and their impact on Wnt-Frizzled receptor signal transduction.

Estimation of distances and map construction using radiation hybrids

A method of estimating distances between pairs of genetic markers is described that directly uses their observed joint frequency distribution in a panel of radiation hybrids (RHs). The distance measure is based on the strength of association between marker pairs, which is high for close markers and decays with distance. These distances are then submitted to a previous method that generates linear coordinates for the markers directly from the intermarker distance matrix. This method of map building from RH data is simpler than others, because it uses only the observed joint frequency distributions of markers in the panel, and does not attempt to model unobserved quantities such as the retention of different sized fragments that contain the markers. It also incorporates directly the observed variation in retention of different markers, without needing a model for differential fragment retention dependent on chromosomal location, which is generally not known. Only small, precise distances are used in map construction, thereby reducing any effects of different fragment retention frequencies and local variations in X-ray sensitivity. The method is tested by simulation, and known marker distances and locations are successfully recovered from RH raw data. The method is also applied to publicly available data sets related to the recent transcript map of the human genome.

Building human genome maps with radiation hybrids

Genome maps are crucial tools in human genetic research, providing known landmarks for locating disease genes and frameworks for large-scale sequencing. Radiation hybrid mapping is one technique for building genome maps. In this paper, we describe the methods used to build radiation hybrid maps of the entire human genome. We present the hidden Markov model that we employ to estimate the likelihood of a map despite uncertainty about the data, and we discuss the problem of searching for maximum-likelihood maps. We describe the graph algorithms used to find sparse but reliable initial maps and our methods of extending them. Finally, we show results validating our software on simulated data, and we describe our genome-wide human radiation hybrid maps and the evidence supporting them.

Markov chain Monte Carlo methods for radiation hybrid mapping

The ordering of genetic loci is central to genetic mapping at all levels. Markov chain Monte Carlo (MCMC) techniques can provide estimates of the posterior density of orders while accounting naturally for missing data, data errors, and unknown parameters. MCMC sampling schemes have been proposed for mapping problems such as linkage mapping and radiation hybrid mapping. The sampling schemes tend, however, to suffer from poor mixing caused by strong correlations between the model parameters. The method described here investigates the effect of using a modified sampling scheme, simulated tempering, on the mixing characteristics of the Markov chain. The method is illustrated by the analysis of haploid radiation hybrid mapping data; the principles are, however, applicable to a range of mapping problems. The results demonstrate that simulated tempering greatly improves the performance of the MCMC sampling scheme. For the radiation hybrid problem, the approach is probably not suitable for simultaneously ordering very large number of loci (> 100); it could, however, be useful for fine scale mapping of subsections of chromosomes.

On constructing radiation hybrid maps

Radiation hybrid (RH) mapping is a somatic cell method for obtaining ordering information of markers on a chromosome, using relatively few experiments. Given the results of a typical RH experiment, finding the true order of the markers is a challenging algorithmic problem. In this work we present several simple algorithms for ordering and mapping the markers, where the input is the genomic data obtained from RH experiments. We provide a rigorous analysis of these algorithms. In particular, we show that under the standard statistical model for RH, our algorithms are "statistically consistent." That is, given enough hybrids, the algorithms do reconstruct the true markers' order (with high probability). We also prove a simple lower bound for the number of hybrids required (by any algorithm) to correctly reconstruct the order. We have implemented these algorithms, and tested them on synthetic and real data. These simulations show that for practical input sizes (number of markers and hybrids) our algorithms produce outputs that are very close to the true ordering. The simulations also indicate that the true ordering of the markers is usually not the one which minimizes the number of obligate chromosome breaks.

A radiation hybrid framework map of bovine chromosome 13

In this paper we present a 5000-rad radiation hybrid framework map of bovine chromosome 13 (BTA13) containing 13 loci, including five conserved genes and eight polymorphic microsatellites. All framework markers are ordered with odds greater than 1000:1. Furthermore, we present a comprehensive map of BTA13 integrating 11 genes and 16 microsatellites. The proposed order is in general agreement with the recently published medium-density linkage maps. A model of five blocks of genes with conserved order between human, mouse and cattle is presented.

Zooming in on the human-mouse comparative map: genome conservation re-examined on a high-resolution scale

Over the past decade, conservation of genetic linkage groups has been shown in mammals and used to great advantage, fueling significant exchanges of gene mapping and functional information especially between the genomes of humans and mice. As human physical maps increase in resolution from chromosome bands to nucleotide sequence, comparative alignments of mouse and human regions have revealed striking similarities and surprising differences between the genomes of these two best-mapped mammalian species. Whereas, at present, very few mouse and human regions have been compared on the physical level, existing studies provide intriguing insights to genome evolution, including the observation of recent duplications and deletions of genes that may play significant roles in defining some of the biological differences between the two species. Although high-resolution conserved marker-based maps are currently available only for human and mouse, a variety of new methods and resources are speeding the development of comparative maps of additional organisms. These advances mark the first step toward establishment of the human genome as a reference map for vertebrate species, providing evolutionary and functional annotation to human sequence and vast new resources for genetic analysis of a variety of commercially, medically, and ecologically important animal models.

An STS-based radiation hybrid map of the human genome

We have constructed a physical map of the human genome by using a panel of 83 whole genome radiation hybrids (the Stanford G3 panel) in conjunction with 10,478 sequence-tagged sites (STSs) derived from random genomic DNA sequences, previously mapped genetic markers, and expressed sequences. Of these STSs, 5049 are framework markers that fall into 1766 high-confidence bins. An additional 945 STSs are indistinguishable in their map location from one or more of the framework markers. These 5994 mapped STSs have an average spacing of 500 kb. An additional 4484 STSs are positioned with respect to the framework markers. Comparison of the orders of markers on this map with orders derived from independent meiotic and YAC STS-content maps indicates that the error rate in defining high-confidence bins is < 5%. Analysis of 322 random cDNAs indicates that the map covers the vast majority of the human genome. This STS-based radiation hybrid map of the human genome brings us one step closer to the goal of a physical map containing 30,000 unique ordered landmarks with an average marker spacing of 100 kb.

Studies on locus content mapping

Locus content maps are derived from monosomic or disomic chromosomes broken by radiation, shearing, or other clastogen, the fragments being distributed among clones by dilution or incorporation into the cells of another species and scored for segregation of markers. Locus content maps provide evidence about radiosensitivity of chromosome regions, support for order, and approximate location. Omission of the most aberrant and least informative clones increases efficiency of localization. Correct analysis must allow for preferential retention of certain sequences, monosomy or polysomy of donor chromosomes, and error filtration. Combination of these refinements extracts substantially more information from fewer clones. Because of unmodeled peculiarities in the data, the best analysis does not recover the physical map but roughly localizes markers that may be monomorphic and therefore unsuitable for linkage mapping. As with linkage for polymorphic loci, distance in the composite map should be confirmed by physical methods.

Hybrid selection as a method of increasing mapping power for radiation hybrids

Radiation hybrids have become a widely used tool for physical mapping. A drawback of the technique is that large numbers of hybrids are required to construct robust, high-resolution maps. The information contained within a panel of radiation hybrids is limited by the frequency of retention of chromosomal fragments from the donor cell line. In almost all experiments to date, the retention frequency has been below the optimal level; therefore, many hybrids are needed to produce high-quality maps. Because of the labor-intensive nature of large-scale mapping projects, it is important to make panels as small as possible. One method that has been adopted is to produce initially a large number of hybrids that are all typed with a few loci. Those hybrids showing satisfactorily high retention are admitted to the final panel and the rest are discarded. In this way, a panel of radiation hybrids with higher than expected retention can be created. Methods for conducting such a selection regime are discussed. To investigate the potential advantages of selecting hybrids based on their retention frequency, simulations were run under a variety of conditions. As expected panels with high retention (40%) provided better mapping resources than panels with lower (20%) retention. Beginning with an initial panel of 200 hybrids, comparisons of a random sample of 100 hybrids and the set of those 100 hybrids showing the highest marker retention demonstrated that selection may not be always the best strategy despite the increase in mean retention it yields. The selection of hybrids containing large numbers of fragments leads to an overestimation of the frequency of radiation-induced breaks. When breaks occur with high frequency (for example, when high radiation doses are used), the selection of hybrids leads to a loss of linkage and hence an inability to order the markers. As such, the merits of screening hybrids depends on both the radiation dose and the desired map resolution.

Statistical methods for polyploid radiation hybrid mapping

Radiation hybrid mapping is a somatic cell technique for ordering genetic loci along a chromosome and estimating physical distances between adjacent loci. This paper presents a model of fragment generation and retention for data involving two or more copies of the chromosome of interest per clone. Such polyploid data can be generated by initially irradiating normal diploid cells or by pooling haploid or diploid clones. The current model assumes that fragments are generated in the ancestral cell of a clone according to an independent Poisson breakage process along each chromosome. Once generated, fragments are independently retained in the clone with a common retention probability. On the basis of this and less restrictive retention models, statistical criteria such as minimum obligate breaks, maximum likelihood ratios, and Bayesian posterior probabilities can be used to decide locus order. Distances can be estimated by maximum likelihood. Likelihood computation is particularly challenging, and computing techniques from the theory of hidden Markov chains prove crucial. Within this context it is possible to incorporate typing errors. The statistical tools discussed here are applied to 14 loci on the short arm of human chromosome 4.