Uniform integration of genome mapping data using intersection graphs

MOTIVATION: The methods for analyzing overlap data are distinct from those for analyzing probe data, making integration of the two forms awkward. Conversion of overlap data to probe-like data elements would facilitate comparison and uniform integration of overlap data and probe data using software developed for analysis of STS data. RESULTS: We show that overlap data can be effectively converted to probe-like data elements by extracting maximal sets of mutually overlapping clones. We call these sets virtual probes, since each set determines a site in the genome corresponding to the region which is common among the clones of the set. Finding the virtual probes is equivalent to finding the maximal cliques of a graph. We modify a known maximal-clique algorithm such that it finds all virtual probes in a large dataset within minutes. We illustrate the algorithm by converting fingerprint and Alu-PCR overlap data to virtual probes. The virtual probes are then analyzed using double-linkage intersection graphs and structure graphs to show that methods designed for STS data are also applicable to overlap data represented as virtual probes. Next we show that virtual probes can produce a uniform integration of different kinds of mapping data, in particular STS probe data and fingerprint and Alu-PCR overlap data. The integrated virtual probes produce longer double-linkage contigs than STS probes alone, and in conjunction with structure graphs they facilitate the identification and elimination of anomalies. Thus, the virtual-probe technique provides: (i) a new way to examine overlap data; (ii) a basis on which to compare overlap data and probe data using the same systems and standards; and (iii) a unique and useful way to uniformly integrate overlap data with probe data.     

Chromosomal localization of ESTs obtained from human fetal liver via BAC-mediated FISH mapping.

A total of 55 expressed sequence tags (ESTs) randomly chosen from our collection of fetal liver ESTs were mapped to chromosomes by fluorescence in situ hybridization (FISH) mapping techniques. To generate FISH mapping probes, the genomic DNAs for each EST were selected by screening an arrayed human bacterial artificial chromosome (BAC) library. In total, 73 BACs were used for mapping of the 55 ESTs. Among them, 70 BACs representing 52 ESTs unequivocally mapped to single chromosomal regions. The remaining 3 BACs representing 3 ESTs were localized to multiple regions, suggesting that BACs may have very low chimerism. Our mapping results were compared with EST mapping databases deposited in NCBI. Thirty-six of 55 ESTs corresponded to previously mapped positions of ESTs, 2 ESTs mapped to different positions from previously determined ones, and it was found that 17 ESTs have been mapped on new locations from this study. These mapping data may be used for completing the framework of the human physical map, and also for providing a good starting point for searching disease-related genes.     

Integration of animal linkage and BAC contig maps using overgo hybridization

The alignment of genome linkage maps, defined primarily by segregation of sequence-tagged site (STS) markers, with BAC contig physical maps and full genome sequences requires high throughput mechanisms to identify BAC clones that contain specific STS. A powerful technique for this purpose is multi-dimensional hybridization of "overgo" probes. The probes are chosen from available STS sequence data by selecting unique probe sequences that have a common melting temperature. We have hybridized sets of 216 overgo probes in subset pools of 36 overgos at a time to filter-spotted chicken BAC clone arrays. A four-dimensional pooling strategy, including one degree of redundancy, has been employed. This requires 24 hybridizations to completely assign BACs for all 216 probes. Results to date are consistent with about a 10% failure rate in overgo probe design and a 15-20% false negative detection rate within a group of 216 markers. Three complete rounds of overgo hybridization, each to sets of about 39,000 BACs (either BAMHI or ECORI partial digest inserts) generated a total of 1853 BAC alignments for 517 mapped chicken genome STS markers. These data are publicly available, and they have been used in the assembly of a first generation BAC contig map of the chicken genome.     

Mapping human chromosomes by walking with sequence-tagged sites from end fragments of yeast artificial chromosome inserts.

Sequence-tagged sites (STSs) derived from end fragments of chromosome-specific yeast artificial chromosomes (YACs) can facilitate the assembly of an overlapping YAC/STS map. Contigs form rapidly by iteratively screening YAC collections with end-fragment STSs from YACs that have not yet been detected by any previous STS. The map is rendered rapidly useful during its assembly by incorporating supplementary STSs from genes and genetic linkage probes with known locations. Methods for the systematic development and testing of the end-fragments STSs are given here, and a group of 100 STSs is presented for the X chromosome. The mapping strategy is shown to be successful in simulations with portions of the X chromosome already largely mapped into overlapping YACs by other means.     

Isolation of sequences from Xp22.3 and deletion mapping using sex chromosome rearrangements from human X-Y interchange sex reversals

A repeated DNA element (STIR) interspersed in Xp22.3 and on the Y chromosome has been used as a tag to isolate seven single-copy probes from the human sex chromosomes. The seven probes detect X-specific loci located in Xp22.3. Using a panel of X-chromosomal deletions from X-Y interchange sex reversals (XX males and XY females), these X-specific loci and some additional ones were mapped to four contiguous intervals of Xp22.3, proximal to the pseudoautosomal region and distal to STS. The construction of this deletion map of the terminal part of the human X chromosome can serve as a starting point for a long-range physical map of Xp22.3 and for a more accurate mapping of genetic diseases located in Xp22.3.     

Reconstructing an ultrametric galled phylogenetic network from a distance matrix

Given a distance matrix M that specifies the pairwise evolutionary distances between n species, the phylogenetic tree reconstruction problem asks for an edge-weighted phylogenetic tree that satisfies M, if one exists. We study some extensions of this problem to rooted phylogenetic networks. Our main result is an O(n(2) log n)-time algorithm for determining whether there is an ultrametric galled network that satisfies M, and if so, constructing one. In fact, if such an ultrametric galled network exists, our algorithm is guaranteed to construct one containing the minimum possible number of nodes with more than one parent (hybrid nodes). We also prove that finding a largest possible submatrix M' of M such that there exists an ultrametric galled network that satisfies M' is NP-hard. Furthermore, we show that given an incomplete distance matrix (i.e. where some matrix entries are missing), it is also NP-hard to determine whether there exists an ultrametric galled network which satisfies it.     

A 3347-locus genetic recombination map of sequence-tagged sites reveals features of genome organization, transmission and evolution of cotton (Gossypium)

We report genetic maps for diploid (D) and tetraploid (AtDt) Gossypium genomes composed of sequence-tagged sites (STS) that foster structural, functional, and evolutionary genomic studies. The maps include, respectively, 2584 loci at 1.72-cM ( approximately 600 kb) intervals based on 2007 probes (AtDt) and 763 loci at 1.96-cM ( approximately 500 kb) intervals detected by 662 probes (D). Both diploid and tetraploid cottons exhibit negative crossover interference; i.e., double recombinants are unexpectedly abundant. We found no major structural changes between Dt and D chromosomes, but confirmed two reciprocal translocations between At chromosomes and several inversions. Concentrations of probes in corresponding regions of the various genomes may represent centromeres, while genome-specific concentrations may represent heterochromatin. Locus duplication patterns reveal all 13 expected homeologous chromosome sets and lend new support to the possibility that a more ancient polyploidization event may have predated the A-D divergence of 6-11 million years ago. Identification of SSRs within 312 RFLP sequences plus direct mapping of 124 SSRs and exploration for CAPS and SNPs illustrate the "portability" of these STS loci across populations and detection systems useful for marker-assisted improvement of the world's leading fiber crop. These data provide new insights into polyploid evolution and represent a foundation for assembly of a finished sequence of the cotton genome.     

High-resolution melt analysis to identify and map sequence-tagged site anchor points onto linkage maps: a white lupin (Lupinus albus) map as an exemplar

* The provision of sequence-tagged site (STS) anchor points allows meaningful comparisons between mapping studies but can be a time-consuming process for nonmodel species or orphan crops. * Here, the first use of high-resolution melt analysis (HRM) to generate STS markers for use in linkage mapping is described. This strategy is rapid and low-cost, and circumvents the need for labelled primers or amplicon fractionation. * Using white lupin (Lupinus albus, x = 25) as a case study, HRM analysis was applied to identify 91 polymorphic markers from expressed sequence tag (EST)-derived and genomic libraries. Of these, 77 generated STS anchor points in the first fully resolved linkage map of the species. The map also included 230 amplified fragment length polymorphisms (AFLP) loci, spanned 1916 cM (84.2% coverage) and divided into the expected 25 linkage groups. * Quantitative trait loci (QTL) analyses performed on the population revealed genomic regions associated with several traits, including the agronomically important time to flowering (tf), alkaloid synthesis and stem height (Ph). Use of HRM-STS markers also allowed us to make direct comparisons between our map and that of the related crop, Lupinus angustifolius, based on the conversion of RFLP, microsatellite and single nucleotide polymorphism (SNP) markers into HRM markers.     

Steroid sulfatase and estrogen sulfotransferase in normal human tissue and breast carcinoma

Steroid sulfatase (STS) hydrolyzes inactive estrone sulfate (E1-S) to estrone (E1), while estrogen sulfotransferase (EST; SULT 1E1 or STE gene) sulfonates estrogens to estrogen sulfates. They are considered to play important roles in the regulation of local estrogenic actions in various human tissues, however, their biological significance remains largely unknown. Therefore, we examined the expression of STS and EST in non-pathologic human tissues and breast carcinomas. STS expression was very weak except for the placenta, while EST expression was markedly detected in various tissues examined. In breast carcinoma tissues, STS and EST immunoreactivity was detected in carcinoma cells in 74 and 44% of cases, respectively, and was significantly associated with their mRNA levels and enzymatic activities. STS immunoreactivity was significantly correlated with the tumor size, and an increased risk of recurrence. EST immunoreactivity was inversely correlated with the tumor size or lymph node status. Moreover, EST immunoreactivity was significantly associated with a decreased risk of recurrence or improved prognosis. Our results suggest that EST is involved in protecting various peripheral tissues from excessive estrogenic effects. In the breast carcinoma, STS and EST are suggested to play important roles in the regulation of in situ estrogen production in the breast carcinomas.     

ProbeMatchDB--a web database for finding equivalent probes across microarray platforms and species

SUMMARY: ProbeMatchDB is a web-based database designed to facilitate the search of EST/cDNA sequences or STS markers that can be used to represent the same gene across different microarray platforms and species. It can be used for finding equivalent EST clones in the Research Genetics sequence verified clone set based on results from Affymetirx GeneChips. It will also help to identify probes representing orthologous genes across human, mouse and rat on different microarray platforms. AVAILABILITY: The database is accessible at http://brainarray.mhri.med.umich.edu/MARRAY/BC_ASP/brainarray.htm by clicking the 'Query ProbeMatchDB' link.     

High-resolution, human–bovine comparative mapping based on a closed YAC contig spanning the bovine mh locus

A closed YAC contig spanning the mh locus was assembled by STS content mapping with seven microsatellite markers, eight genes or EST, and nine STS corresponding to YAC ends. The contig comprises 27 YACs, has an average depth of 4.3 YACs, and spans an estimated 1.2 Mb. A linkage map was constructed based on five of the microsatellite markers anchored to the contig and shown to span 7 cM, yielding a ratio of 160 kb/1 cM for the corresponding chromosome region. Comparative mapping data indicate that the constructed contig spans an evolutionary breakpoint connecting two chromosome segments that are syntenic but not adjacent in the human. Consolidation of human gene order by means of whole genome radiation hybrids and its comparison with the bovine order as inferred from the contig confirm conservation of gene order within segments. Received: 6 August 1998 / Accepted: 28 October 1998     

The chicken RH map: current state of progress and microchromosome mapping

The ChickRH6 radiation hybrid panel has been used to construct consensus chromosome radiation hybrid (RH) maps of the chicken genome. Markers genotyped were either from throughout the genome or targeted to specific chromosomes and a large proportion (one third) of data was the result of collaborative efforts. Altogether, 2,531 markers were genotyped, allowing the construction of RH reference maps for 20 chromosomes and linkage groups for four other chromosomes. Amongst the markers, 581 belong to the framework maps, while 1,721 are on the comprehensive maps. Around 800 markers still have to be assigned to linkage groups. Our attempt to assign the supercontigs from the chrun (virtual chromosome containing all the genome sequence that could not be attributed to a chromosome) as well as EST (Expressed Sequence Tag) contigs that do not have a BLAST hit in the genome assembly led to the construction of new maps for microchromosomes either absent or for which very little data is present in the genome assembly. RH data is presented through our ChickRH webserver (http://chickrh.toulouse.inra.fr/), which is a mapping tool as well as the official repository RH database for genotypes. It also displays the RH reference maps and comparison charts with the sequence thus highlighting the possible discrepancies. Future improvements of the RH maps include complete coverage of the sequence assigned to chromosomes, further mapping of the chrun and mapping of EST contigs absent from the assembly. This will help finish the mapping of the smallest gene-rich microchromosomes.     

Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.     

Generation of a 5.5-Mb BAC/PAC contig of pig chromosome 6q1.2 and its integration with existing RH, genetic and comparative maps

We generated a sequence-ready BAC/PAC contig spanning approximately 5.5 Mb on porcine chromosome 6q1.2, which represents a very gene-rich genome region. STS content mapping was used as the main strategy for the assembly of the contig and a total of 6 microsatellite markers, 53 gene-related STS and 116 STS corresponding to BAC and PAC end sequences were analyzed. The contig comprises 316 BAC and PAC clones covering the region between the genes GPI and LIPE. The correct contig assembly was verified by RH-mapping of STS markers and comparative mapping of BAC/PAC end sequences using BLAST searches. The use of microsatellite primer pairs allowed the integration of the physical maps with the genetic map of this region. Comparative mapping of the porcine BAC/PAC contig with respect to the gene-rich region on the human chromosome 19q13.1 map revealed a completely conserved gene order of this segment, however, physical distances differ somewhat between HSA19q13.1 and SSC6q1.2. Three major differences in DNA content between human and pig are found in two large intergenic regions and in one region of a clustered gene family, respectively. While there is a complete conservation of gene order between pig and human, the comparative analysis with respect to the rodent species mouse and rat shows one breakpoint where a genome segment is inverted.     

Few-sweep estimation of evoked potential based on a generalized subspace approach

A generalized subspace approach is proposed for single channel brain evoked potential (EP) extraction from background electroencephalogram (EEG) signal. The method realizes the optimum estimate of EP signal from the observable noisy signal. The underlying principle is to project the signal and noise into signal and noise coefficient subspace respectively by applying projection matrix at first. Secondly, coefficient weighting matrix is achieved based on the autocorrelation matrices of the noise and the noisy signal. With the coefficient weighting matrix, we can remove the noise projection coefficients and estimate the signal ones. EP signal is then obtained by averaging the signals estimated with the reconstruction matrix. Given different signal-to-noise ratio (SNR) conditions, the algorithm can estimate the EP signal with only two sweeps observable noisy signals. Our approach is shown to have excellent capability of estimating EP signal even in poor SNR conditions. The interference of spontaneous EEG has been eliminated with significantly improved SNR. The simulation results have demonstrated the effectiveness and superior performance of the proposed method.     

基于EST的新基因克隆策略

表达序列标签(expressed sequence tags, EST) 是从随机选择的cDNA 克隆进行单向测序获得的短的cDNA序列, 代表一个完整基因的一部分。随着生物信息学和基因定位的迅猛发展, EST已成为基因定位、基因克隆、基因表达分析的有力工具。近年来, 由于EST数据库的迅速扩张, 运用EST来克隆和定位基因, 使得新基因克隆的策略发生了革命性变革。尽管存在一些不足, 实践证明EST可大大加速新基因的发现与研究。本文将就EST技术尤其是它在新基因克隆中的应用策略作详细介绍。     

Theoretical basis for the identification of allelic variants that encode drug efficacy and toxicity

Almost all drugs that produce a favorable response (efficacy) may also produce adverse effects (toxicity). The relative strengths of drug efficacy and toxicity that vary in human populations are controlled by the combined influences of multiple genes and environmental influences. Genetic mapping has proven to be a powerful tool for detecting and identifying specific DNA sequence variants on the basis of the haplotype map (HapMap) constructed from single-nucleotide polymorphisms (SNPs). In this article, we present a novel statistical model for sequence mapping of two different but related drug responses. This model is incorporated by mathematical functions of drug response to varying doses or concentrations and the statistical device used to model the correlated structure of the residual (co)variance matrix. We implement a closed-form solution for the EM algorithm to estimate the population genetic parameters of SNPs and the simplex algorithm to estimate the curve parameters describing the pharmacodynamic changes of different genetic variants and matrix-structuring parameters. Extensive simulations are performed to investigate the statistical properties of our model. The implications of our model in pharmacogenetic and pharmacogenomic research are discussed.     

Biochemical and Structural Studies of the Large Ycf4-Photosystem I Assembly Complex of the Green Alga Chlamydomonas reinhardtii

Ycf4 is a thylakoid protein essential for the accumulation of photosystem I (PSI) in Chlamydomonas reinhardtii. Here, a tandem affinity purification tagged Ycf4 was used to purify a stable Ycf4-containing complex of >1500 kD. This complex also contained the opsin-related COP2 and the PSI subunits PsaA, PsaB, PsaC, PsaD, PsaE, and PsaF, as identified by mass spectrometry (liquid chromatography–tandem mass spectrometry) and immunoblotting. Almost all Ycf4 and COP2 in wild-type cells copurified by sucrose gradient ultracentrifugation and subsequent ion exchange column chromatography, indicating the intimate and exclusive association of Ycf4 and COP2. Electron microscopy revealed that the largest structures in the purified preparation measure 285 × 185 Å; these particles may represent several large oligomeric states. Pulse-chase protein labeling revealed that the PSI polypeptides associated with the Ycf4-containing complex are newly synthesized and partially assembled as a pigment-containing subcomplex. These results indicate that the Ycf4 complex may act as a scaffold for PSI assembly. A decrease in COP2 to 10% of wild-type levels by RNA interference increased the salt sensitivity of the Ycf4 complex stability but did not affect the accumulation of PSI, suggesting that COP2 is not essential for PSI assembly.     

High-resolution mapping of the barley leaf rust resistance gene Rph5 using barley expressed sequence tags (ESTs) and synteny with rice

The rapidly growing expressed sequence tag (EST) resources of species representing the Poacea family and availability of comprehensive sequence information for the rice (Oryza sativa) genome create an excellent opportunity for comparative genome analysis. Extensive synteny between rice chromosome 1 and barley (Hordeum vulgare L.) chromosome 3 has proven extremely useful for saturation mapping of chromosomal regions containing target genes of large-genome barley with conserved orthologous genes from the syntenic regions of the rice genome. Rph5 is a gene conferring resistance to the barley leaf rust pathogen Puccinia hordei. It was mapped to chromosome 3HS, which is syntenic with rice chromosome 1S. The objective of this study was to increase marker density within the sub-centimorgan region around Rph5, using sequence-tagged site (STS) markers that were developed based on barley ESTs syntenic to the phage (P1)-derived artificial chromosome (PAC) clones comprising the distal region of rice chromosome 1S. Five rice PAC clones were used as queries in a blastn search to screen 375,187 barley ESTs. Ninety-four non-redundant EST sequences were identified from the EST database and used as templates to design 174 pairs of primer combinations. As a result, 9 barley EST-based STS markers were incorporated into the ‘Bowman’ × ‘Magnif 102’ high-resolution map of the Rph5 region. More importantly, six markers, including five EST-derived STS sequences, were found to co-segregate with Rph5. The results of this study demonstrate the usefulness of rice genomic resources for efficient deployment of barley ESTs for marker saturation of targeted barley genomic regions.