Around the genomes: the Drosophila genome project

Tracking microRNA Processing Signals by Degradome Sequencing Data Analysis

Degradome sequencing (degradome-seq) was widely used for cleavage site mapping on the microRNA (miRNA) targets. Here, the application value of degradome-seq data in tracking the miRNA processing intermediates was reported. By adopting the parameter “signal/noise” ratio, prominent degradome signals on the miRNA precursors were extracted. For the 15 species analyzed, the processing of many miRNA precursors were supported by the degradome-seq data. We found that the supporting ratio of the “high-confidence” miRNAs annotated in miRBase was much higher than that of the “low-confidence.” For a specific species, the percentage of the miRNAs with degradome-supported processing signals was elevated by the increment of degradome sampling diversity. More interestingly, the tissue- or cell line-specific processing patterns of the miRNA precursors partially contributed to the accumulation patterns of the mature miRNAs. In this study, we also provided examples to show the value of the degradome-seq data in miRNA annotation. Based on the distribution of the processing signals, a renewed model was proposed that the stems of the miRNA precursors were diced through a “single-stranded cropping” mode, and “loop-to-base” processing was much more prevalent than previously thought. Together, our results revealed the remarkable capacity of degradome-seq in tracking miRNA processing signals.

Molecular Population Genetics

Molecular population genetics aims to explain genetic variation and molecular evolution from population genetics principles. The field was born 50 years ago with the first measures of genetic variation in allozyme loci, continued with the nucleotide sequencing era, and is currently in the era of population genomics. During this period, molecular population genetics has been revolutionized by progress in data acquisition and theoretical developments. The conceptual elegance of the neutral theory of molecular evolution or the footprint carved by natural selection on the patterns of genetic variation are two examples of the vast number of inspiring findings of population genetics research. Since the inception of the field, Drosophila has been the prominent model species: molecular variation in populations was first described in Drosophila and most of the population genetics hypotheses were tested in Drosophila species. In this review, we describe the main concepts, methods, and landmarks of molecular population genetics, using the Drosophila model as a reference. We describe the different genetic data sets made available by advances in molecular technologies, and the theoretical developments fostered by these data. Finally, we review the results and new insights provided by the population genomics approach, and conclude by enumerating challenges and new lines of inquiry posed by increasingly large population scale sequence data.

Immune-directed support of rich microbial communities in the gut has ancient roots

The animal gut serves as a primary location for the complex host-microbe interplay that is essential for homeostasis and may also reflect the types of ancient selective pressures that spawned the emergence of immunity in metazoans. In this review, we present a phylogenetic survey of gut host-microbe interactions and suggest that host defense systems arose not only to protect tissue directly from pathogenic attack but also to actively support growth of specific communities of mutualists. This functional dichotomy resulted in the evolution of immune systems much more tuned for harmonious existence with microbes than previously thought, existing as dynamic but primarily cooperative entities in the present day. We further present the protochordate Ciona intestinalis as a promising model for studying gut host-bacterial dialogue. The taxonomic position, gut physiology and experimental tractability of Ciona offer unique advantages in dissecting host-microbe interplay and can complement studies in other model systems.

Drosophila RNAi screening in a postgenomic world

Drosophila melanogaster has a long history as a model organism with several unique features that make it an ideal research tool for the study of the relationship between genotype and phenotype. Importantly fundamental genetic principles as well as key human disease genes have been uncovered through the use of Drosophila. The contribution of the fruit fly to science and medicine continues in the postgenomic era as cell-based Drosophila RNAi screens are a cost-effective and scalable enabling technology that can be used to quantify the contribution of different genes to diverse cellular processes. Drosophila high-throughput screens can also be used as integral part of systems-level approaches to describe the architecture and dynamics of cellular networks.

Drosophila melanogaster as a model for studying the function of animal viral proteins

Studies in which Drosophila melanogaster individuals carrying transgenes of animal viruses were used to analyze the action of animal viral proteins on the cell are reviewed. The data presented suggest that host specificity of viruses is determined by their proteins responsible for the penetration of the virus into the cell, while viral proteins responsible for interactions with the host cell are much less host-specific. Due to this, the model of Drosophila with its developed system of searching for genetic interactions can be used to find intracellular targets for the action of viral proteins of the second group.

On the mechanics of cardiac function of Drosophila embryo

The heart is a vital organ that provides essential circulation throughout the body. Malfunction of cardiac pumping, thus, leads to serious and most of the times, to fatal diseases. Mechanics of cardiac pumping is a complex process, and many experimental and theoretical approaches have been undertaken to understand this process. We have taken advantage of the simplicity of the embryonic heart of an invertebrate, Drosophila melanogaster, to understand the fundamental mechanics of the beating heart. We applied a live imaging technique to the beating embryonic heart combined with analytical imaging tools to study the dynamic mechanics of the pumping. Furthermore, we have identified one mutant line that exhibits aberrant pumping mechanics. The Drosophila embryonic heart consists of only 104 cardiac cells forming a simple straight tube that can be easily accessed for real-time imaging. Therefore, combined with the wealth of available genetic tools, the embryonic Drosophila heart may serve as a powerful model system for studies of human heart diseases, such as arrhythmia and congenital heart diseases. We, furthermore, believe our mechanistic data provides important information that is useful for our further understanding of the design of biological structure and function and for engineering the pumps for medical uses.

graal: a Drosophila gene coding for several mosaic serine proteases

Serine proteases play vital roles in several biological processes such as development and immunity. We have characterized Graal, a large multi-domain serine protease from Drosophila. Graal is spliced in at least three transcripts that are present throughout development. The domains found in Graal proteins are: chitin-binding domains (CBD), scavenger receptor cysteine-rich (SRCR) domains, low density lipoprotein receptor cysteine-rich (LDLR-CR) domains, histidine and proline-rich domains, a NGGYQPP-repeat domain and a serine protease domain. The last 2370 nucleotides of these RNAs are identical and encode a His-rich domain, two SRCR domains, two LDLR-CR domains and a protease domain. The transcription of graal is upregulated after fungal or bacterial infection. Analysis of the Iso1 (y;cn,sp,bw) strain shows that graal transcription is impaired in this fly line due to the insertion of a retrotransposon in the sixth exon. However, no phenotype could be observed consecutive to the absence of graal full length transcripts, particularly in the context of an immune challenge.

ATIDB: Arabidopsis thaliana insertion database

Insertional mutagenesis techniques, including transposon- and T-DNA-mediated mutagenesis, are key resources for systematic identification of gene function in the model plant species Arabidopsis thaliana. We have developed a database (http://atidb.cshl.org/) for archiving, searching and analyzing insertional mutagenesis lines. Flanking sequences from approximately 10 500 insertion lines (including transposon and T-DNA insertions) from several tagging programs in Arabidopsis were mapped to the genome sequence through our annotation system before being entered into the database. The database front end provides World Wide Web searching and analyzing interfaces for genome researchers and other biologists. Users can search the database to identify insertions in a particular gene or perform genome-wide analysis to study the distribution and preference of insertions. Tools integrated with the database include a graphical genome browser, a protein search function, a graphical representation of the insertion distribution and a Blast search function. The database is based on open source components and is available under an open source license.

Integrated graphical analysis of protein sequence features predicted from sequence composition

Several protein sequence analysis algorithms are based on properties of amino acid composition and repetitiveness. These include methods for prediction of secondary structure elements, coiled-coils, transmembrane segments or signal peptides, and for assignment of low-complexity, nonglobular, or intrinsically unstructured regions. The quality of such analyses can be greatly enhanced by graphical software tools that present predicted sequence features together in context and allow judgment to be focused simultaneously on several different types of supporting information. For these purposes, we describe the SFINX package, which allows many different sets of segmental or continuous-curve sequence feature data, generated by individual external programs, to be viewed in combination alongside a sequence dot-plot or a multiple alignment of database matches. The implementation is currently based on extensions to the graphical viewers Dotter and Blixem and scripts that convert data from external programs to a simple generic data definition format called SFS. We describe applications in which dot-plots and flanking database matches provide valuable contextual information for analyses based on compositional and repetitive sequence features. The system is also useful for comparing results from algorithms run with a range of parameters to determine appropriate values for defaults or cutoffs for large-scale genomic analyses.

Map position and expression of the genes in the 38 region of Drosophila

With the completion of the Drosophila genome sequence, an important next step is to extract its biological information by systematic functional analysis of genes. We have produced a high-resolution genetic map of cytological region 38 of Drosophila using 41 deficiency stocks that provide a total of 54 breakpoints within the region. Of a total of 45 independent P-element lines that mapped by in situ hybridization to the region, 14 targeted 7 complementation groups within the 38 region. Additional EMS, X-ray, and spontaneous mutations define a total of 17 complementation groups. Because these two pools partially overlap, the completed analysis revealed 21 distinct complementation groups defined by point mutations. Seven additional functions were defined by trans-heterozygous combinations of deficiencies, resulting in a total of 28 distinct functions. We further produced a developmental expression profile for the 760 kb from 38B to 38E. Of 135 transcription units predicted by GENSCAN, 22 have at least partial homology to mobile genetic elements such as transposons and retroviruses and 17 correspond to previously characterized genes. We analyzed the developmental expression pattern of the remaining genes using poly(A)(+) RNA from ovaries, early and late embryos, larvae, males, and females. We discuss the correlation between GENSCAN predictions and experimentally confirmed transcription units, the high number of male-specific transcripts, and the alignment of the genetic and physical maps in cytological region 38.

From first base: the sequence of the tip of the X chromosome of Drosophila melanogaster, a comparison of two sequencing strategies

We present the sequence of a contiguous 2.63 Mb of DNA extending from the tip of the X chromosome of Drosophila melanogaster. Within this sequence, we predict 277 protein coding genes, of which 94 had been sequenced already in the course of studying the biology of their gene products, and examples of 12 different transposable elements. We show that an interval between bands 3A2 and 3C2, believed in the 1970s to show a correlation between the number of bands on the polytene chromosomes and the 20 genes identified by conventional genetics, is predicted to contain 45 genes from its DNA sequence. We have determined the insertion sites of P-elements from 111 mutant lines, about half of which are in a position likely to affect the expression of novel predicted genes, thus representing a resource for subsequent functional genomic analysis. We compare the European Drosophila Genome Project sequence with the corresponding part of the independently assembled and annotated Joint Sequence determined through "shotgun" sequencing. Discounting differences in the distribution of known transposable elements between the strains sequenced in the two projects, we detected three major sequence differences, two of which are probably explained by errors in assembly; the origin of the third major difference is unclear. In addition there are eight sequence gaps within the Joint Sequence. At least six of these eight gaps are likely to be sites of transposable elements; the other two are complex. Of the 275 genes in common to both projects, 60% are identical within 1% of their predicted amino-acid sequence and 31% show minor differences such as in choice of translation initiation or termination codons; the remaining 9% show major differences in interpretation.

Constitutive expression of a complement-like protein in toll and JAK gain-of-function mutants of Drosophila

We show that Drosophila expresses four genes encoding proteins with significant similarities with the thiolester-containing proteins of the complement C3/alpha(2)-macroglobulin superfamily. The genes are transcribed at a low level during all stages of development, and their expression is markedly up-regulated after an immune challenge. For one of these genes, which is predominantly expressed in the larval fat body, we observe a constitutive expression in gain-of-function mutants of the Janus kinase (JAK) hop and a reduced inducibility in loss-of-function hop mutants. We also observe a constitutive expression in gain-of-function Toll mutants. We discuss the possible roles of these novel complement-like proteins in the Drosophila host defense.

Function and structure of Drosophila glycans

Through the application of classic organismal genetic strategies, such as mutagenesis and interaction screens, Drosophila melanogaster provides opportunities to understand glycan function. For instance, screens for Drosophila genes that establish dorsal-ventral polarity in the embryo or that influence cellular differentiation through signal modulation have identified putative glycan modifying enzymes. Other genetic and molecular approaches have demonstrated the existence of phylogenetically conserved and novel oligosaccharide processing activities and carbohydrate binding proteins. While the structural characterization of Drosophila oligosaccharide diversity has lagged behind the elucidation of glycan function, landmarks are becoming apparent in the carbohydrate terrain. For instance, O-linked GlcNAc and mucins, spatially and temporally regulated N-linked oligosaccharide expression, glycosphingolipids, heparan sulfate, chondroitin sulfate and polysialic acid have all been described. A major challenge for Drosophila glycobiology is to expand the oligosaccharide structural database while endeavoring to link glycan characterization to functional analysis. The completion of the Drosophila genome sequencing project will yield a broad portfolio of glycosyltransferases, glycan modifying enzymes and lectins requiring characterization. To this end, the great range of genetic tools that allow the controlled spatial and temporal expression of transgenes in Drosophila will permit unprecedented manipulation of glycosylation in a whole organism.

Zwitterionic and acidic glycosphingolipids of the Drosophila melanogaster embryo

Defining glycosphingolipid structures in species amenable to genetic manipulation, such as Drosophila melanogaster, provides a foundation for investigating mechanisms that regulate glycolipid expression. Therefore, eight of the 12 major glycosphingolipids, accounting for 64% of lipid-linked carbohydrate in Drosophila embryos, were purified after separation into acidic and zwitterionic pools. The zwitterionic lipids possess phosphoethanolamine (PEtn) linked to one or more GlcNAc residues and comprise a family of serially related structures. The longest characterized glycolipid, an octaosylceramide, designated Nz28, has the structure: GalNAcbeta, 4(PEtn-6)GlcNAcbeta,3Galbeta,3GalNAcalpha,4Ga lNAcbeta, 4(PEtn-6)GlcNAcbeta,3Manbeta,4GlcbetaCer. Heptaosyl (Nz7), hexaosyl (Nz6), pentaosyl (Nz5) and tetraosyl (Nz4) forms of Nz28, sequentially truncated from the nonreducing terminus, possess only one PEtn moiety. The major acidic lipid, designated Az29, possesses two PEtn moieties and a glucuronic acid linked to a Gal-extended Nz28. Two other acidic glycolipids, Az9 and Az6, exhibit one PEtn moiety and the same hexose and N-acetylhexosamine composition as Az29 and Nz6, respectively. The fully extended Drosophila core oligosaccharide differs from that of other dipterans in the linkage at a single glycosidic bond, a distinction with significant structural and biosynthetic consequences. Furthermore, acidic species account for a larger proportion of total glycosphingolipid, and PEtn substitution of GlcNAc is more complete in the Drosophila embryo. Divergent characteristics may reflect interspecies variation or stage-specific glycosphingolipid expression in dipterans.

A genetic analysis of cyromazine resistance in Drosophila melanogaster (Diptera: Drosophilidae)

Flies resistant to the insect growth regulator cyromazine were selected in the F1 generation from a cyromazine-susceptible strain of Drosophila melanogaster (Meigen) treated with ethyl methanesulfonate. Four resistant strains were isolated by screening with cyromazine at a concentration > LC100 of susceptibles. In each strain, resistance is conferred by a single gene mutation. Cyromazine resistance in two of the mutants (rst(1a)cyr1 and rst(1a)cyr2) localizes to map position 17 of the X chromosome. Evidence is presented that these mutations are alleles of the gene rst(1a)cyr. Cyromazine resistance in another of the mutants (Rst(1b)Cyr) is also X-linked, and localizes to map position 49 of the X chromosome. The location of the gene conferring cyromazine resistance in the other mutant (Rst(2b)Cyr) is map position 66 of chromosome II. This is possibly an allele of a previously characterized cyromazine resistance gene, Rst(2)Cyr. Dosage-mortality analyses demonstrate a low level of cyromazine resistance is conferred in all strains.

X-linked late-onset sensorineural deafness caused by a deletion involving OA1 and a novel gene containing WD-40 repeats

We have identified a novel gene, transducin (beta)-like 1 (TBL1), in the Xp22.3 genomic region, that shows high homology with members of the WD-40-repeat protein family. The gene contains 18 exons spanning approximately 150 kb of the genomic region adjacent to the ocular albinism gene (OA1) on the telomeric side. However, unlike OA1, TBL1 is transcribed from telomere to centromere. Northern analysis indicates that TBL1 is ubiquitously expressed, with two transcripts of approximately 2.1 kb and 6.0 kb. The open reading frame encodes a 526-amino acid protein, which shows the presence of six beta-transducin repeats (WD-40 motif) in the C-terminal domain. The homology with known beta-subunits of G proteins and other WD-40-repeat containing proteins is restricted to the WD-40 motif. Genomic analysis revealed that the gene is either partly or entirely deleted in patients carrying Xp22.3 terminal deletions. The complexity of the contiguous gene-syndrome phenotype shared by these patients depends on the number of known disease genes involved in the deletions. Interestingly, one patient carrying a microinterstitial deletion involving the 3' portion of both TBL1 and OA1 shows the OA1 phenotype associated with X-linked late-onset sensorineural deafness. We postulate an involvement of TBL1 in the pathogenesis of the ocular albinism with late-onset sensorineural deafness phenotype.

Genetic Localization of a Drosophila melanogaster Resistance Gene to a Parasitoid Wasp and Physical Mapping of the Region

Drosophila melanogaster larvae usually react against eggs of the parasitoid wasp Leptopilina boulardi by surrounding them with a multicellular melanotic capsule. The genetic determinism of this response has been studied previously using susceptible (non-capsule-forming) and resistant (capsule-forming) strains. The results suggest that differences in their encapsulation response involve a single gene, resistance to Leptopilina boulardi(Rlb), with two alleles, the resistant one being dominant.Rlb confers specific protection against Leptopilina boulardi and is thus probably involved in parasitoid recognition. Recent studies have localized this gene on the right arm of the second chromosome and our aim was to precisely determine its genetic and molecular location. Using strains bearing deletions, we demonstrated that resistance to Leptopilina boulardi is conferred by the55C; 55F3 region and that the 55E2–E6; F3 region is particularly involved. A physical map of the 55C;56A region was then constructed, based on a set of overlapping cosmid and P1 phage clones. Using single and double digests, cross hybridization of restriction fragments, and location of genetically mapped genes and STSs, a complete, five-enzyme restriction map of this 830-kb region was obtained.

Computer-aided target selection-prioritizing targets for antifungal drug discovery

The entire DNA sequence of the Saccharomyces cerevisiae genome was completed in 1996 and represents the first entirely decoded eukaryotic genome. Because major human pathogenic fungi such as Candida albicans are closely related to S. cerevisiae on a molecular level, the question arises as to how this new information can be used to identify and prioritize those genes that are most suitable as targets for antimycotic drug discovery. To tackle this challenge, a software tool called CATS (computer-aided target selection) was developed. The authors describe how it allows an automated and periodically updated assessment of all S. cerevisiae genes to be carried out with regard to their suitability as antifungal targets.

A computer program for aligning a cDNA sequence with a genomic DNA sequence

We address the problem of efficiently aligning a transcribed and spliced DNA sequence with a genomic sequence containing that gene, allowing for introns in the genomic sequence and a relatively small number of sequencing errors. A freely available computer program, described herein, solves the problem for a 100-kb genomic sequence in a few seconds on a workstation.

Large-scale taxonomic profiling of eukaryotic model organisms: a comparison of orthologous proteins encoded by the human, fly, nematode, and yeast genomes

Comparisons of DNA and protein sequences between humans and model organisms, including the yeast Saccharomyces cerevisiae, the nematode Caenorhabditis elegans, and the fruit fly Drosophila melanogaster, are a significant source of information about the function of human genes and proteins in both normal and disease states. Important questions regarding cross-species sequence comparison remain unanswered, including (1) the fraction of the metabolic, signaling, and regulatory pathways that is shared by humans and the various model organisms; and (2) the validity of functional inferences based on sequence homology. We addressed these questions by analyzing the available fractions of human, fly, nematode, and yeast genomes for orthologous protein-coding genes, applying strict criteria to distinguish between candidate orthologous and paralogous proteins. Forty-two quartets of proteins could be identified as candidate orthologs. Twenty-four Drosophila protein sequences were more similar to their human orthologs than the corresponding nematode proteins. Analysis of sequence substitutions and evolutionary distances in this data set revealed that most C. elegans genes are evolving more rapidly than Drosophila genes, suggesting that unequal evolutionary rates may contribute to the differences in similarity to human protein sequences. The available fraction of Drosophila proteins appears to lack representatives of many protein families and domains, reflecting the relative paucity of genomic data from this species.

Mammalian BUB1 protein kinases: map positions and in vivo expression

The spindle assembly checkpoint modulates the timing of anaphase initiation in mitotic cells containing improperly aligned chromosomes and increases the probability of successful delivery of a euploid chromosome set to each daughter cell. We have characterized cDNA sequences from several organisms with highly significant predicted protein sequence homologies to Saccharomyces cerevisiae Bub1p, a protein required for function of the spindle assembly checkpoint in budding yeast. The localization of mouse and human orthologs is in agreement with known conservation of synteny. Mouse backcross mapping data indicate that the murine gene resides on chromosome 2 near IL1A, 73 cM from the mouse centromere. Radiation hybrid mapping data indicate that the human locus exhibits linkage to microsatellite marker D2S176, which is located within 10 cM of human IL1A. Multiple-tissue Northern analysis indicates conservation of expression pattern in mouse and human with markedly high mRNA levels in testis. Northern analysis of two different spindle assembly checkpoint protein gene products from human, BUB1 and MAD2, reveals an expression pattern with common tissue distribution consistent with roles in a common pathway. In addition, we demonstrate that an mRNA found to accumulate in a rat fibroblast cell transformation system encodes rat BUB1, and we find that rat BUB1 mRNA accumulation correlates with the proliferation status of cells in culture.

Analysis of protein domain families in Caenorhabditis elegans

The Caenorhabditis elegans genome sequencing project has completed over half of this nematode's 100-Mb genome. Proteins predicted in the finished sequence have been compiled and released in the data-base Wormpep. Presented here is a comprehensive analysis of protein domain families in Wormpep 11, which comprises 7299 proteins. The relative abundance of common protein domain families was counted by comparing all Wormpep proteins to the Pfam collection of protein families, which is based on recognition by hidden Markov models. This analysis also identified a number of previously unannotated domains. To investigate new apparently nematode-specific protein families, Wormpep was clustered into domain families on the basis of sequence similarity using the Domainer program. The largest clusters that lacked clear homology to proteins outside Nematoda were analyzed in further detail, after which some could be assigned a putative function. We compared all proteins in Wormpep 11 to proteins in the human, Saccharomyces cerevisiae, and Haemophilus influenzae genomes. Among the results are the estimation that over two-thirds of the currently known human proteins are likely to have a homologue in the whole C. elegans genome and that a significant number of proteins are well conserved between C. elegans and H. influenzae, that are not found in S. cerevisiae.

Sequence analysis of a 24-kb contiguous genomic region at the Arabidopsis thaliana PFL locus on chromosome 1

As part of the European Union program of European Scientist Sequencing Arabidopsis (ESSA), the DNA sequence of a 24.053-bp insert of cosmid clone CC17J13 was determined. The cosmid is located on chromosome 1 at the PFL locus (position 30 cM). Analysis of the sequence and comparison to public databases predicts seven genes in this area, thus approximately one gene every 3.3 kb. Three cDNAs corresponding to genes in this region were also sequenced. The homologies and/or possible functions of the (putative) genes are discussed. Proteins encoded by genes in this region include a polyadenylate-binding protein (PAB-3) and a GTP-binding protein (Rab7) as well as a novel protein, possibly involved in double-stranded RNA unwinding and apoptosis. Intriguingly, the gene encoding the PAB-3 protein, which is very specifically expressed, is flanked by putative matrix attachment regions.

Improved splice site detection in Genie

We present an improved splice site predictor for the genefinding program Genie. Genie is based on a generalized Hidden Markov Model (GHMM) that describes the grammar of a legal parse of a multi-exon gene in a DNA sequence. In Genie, probabilities are estimated for gene features by using dynamic programming to combine information from multiple content and signal sensors, including sensors that integrate matches to homologous sequences from a database. One of the hardest problems in genefinding is to determine the complete gene structure correctly. The splice site sensors are the key signal sensors that address this problem. We replaced the existing splice site sensors in Genie with two novel neural networks based on dinucleotide frequencies. Using these novel sensors, Genie shows significant improvements in the sensitivity and specificity of gene structure identification. Experimental results in tests using a standard set of annotated genes showed that Genie identified 86% of coding nucleotides correctly with a specificity of 85%, versus 80% and 84% in the older system. In further splice site experiments, we also looked at correlations between splice site scores and intron and exon lengths, as well as at the effect of distance to the nearest splice site on false positive rates.

A physical map of human chromosome 7: an integrated YAC contig map with average STS spacing of 79 kb

The construction of highly integrated and annotated physical maps of human chromosomes represents a critical goal of the ongoing Human Genome Project. Our laboratory has focused on developing a physical map of human chromosome 7, a approximately 170-Mb segment of DNA that corresponds to an estimated 5% of the human genome. Using a yeast artificial chromosome (YAC)-based sequence-tagged site (STS)-content mapping strategy, 2150 chromosome 7-specific STSs have been established and mapped to a collection of YACs highly enriched for chromosome 7 DNA. The STSs correspond to sequences generated from a variety of DNA sources, with particular emphasis placed on YAC insert ends, genetic markers, and genes. The YACs include a set of relatively nonchimeric clones from a human-hamster hybrid cell line as well as clones isolated from total genomic libraries. For map integration, we have localized 260 STSs corresponding to Genethon genetic markers and 259 STSs corresponding to markers orders by radiation hybrid (RH) mapping on our YAC contigs. Analysis of the data with the program SEGMAP results in the assembly of 22 contigs that are "anchored" on the Genethon genetic map, the RH map, and/or the cytogenetic map. These 22 contigs are ordered relative to one another, are (in all but 3 cases) oriented relative to the centromere and telomeres, and contain > 98% of the mapped STSs. The largest anchored YAC contig, accounting for most of 7p, contains 634 STSs and 1260 YACs. An additional 14 contigs, accounting for approximately 1.5% of the mapped STSs, are assembled but remain unanchored on either the genetic or RH map. Therefore, these 14 "orphan" contigs are not ordered relative to other contigs. In our contig maps, adjacent STSs are connected by two or more YACs in > 95% of cases. With 2150 mapped STSs, our map provides an average STS spacing of approximately 79 kb. The physical map we report here exceeds the goal of 100-kb average STS spacing and should provide an excellent framework for systematic sequencing of the chromosome.