Identifying muscle regulatory elements and genes in the nematode Caenorhabditis elegans.

We report the identification of several putative muscle-specific regulatory elements, and genes which are expressed preferentially in the muscle of the nematode Caenorhabditis elegans. We used computational pattern finding methods to identify cis-regulatory motifs from promoter regions of a set of genes known to express preferentially in muscle; each motif describes the potential binding sites for an unknown regulatory factor. The significance and specificity of the identified motifs were evaluated using several different control sequence sets. Using the motifs, we searched the entire C. elegans genome for genes whose promoter regions have a high probability of being bound by the putative regulatory factors. Genes that met this criterion and were not included in our initial set were predicted to be good candidates for muscle expression. Some of these candidates are additional, known muscle expressed genes and several others are shown here to be preferentially expressed in muscle cells by using GFP (green fluorescent protein) constructs. The methods described here can be used to predict the spatial expression pattern of many uncharacterized genes.

The AZFc region of the Y chromosome features massive palindromes and uniform recurrent deletions in infertile men.

Deletions of the AZFc (azoospermia factor c) region of the Y chromosome are the most common known cause of spermatogenic failure. We determined the complete nucleotide sequence of AZFc by identifying and distinguishing between near-identical amplicons (massive repeat units) using an iterative mapping-sequencing process. A complex of three palindromes, the largest spanning 3 Mb with 99.97% identity between its arms, encompasses the AZFc region. The palindromes are constructed from six distinct families of amplicons, with unit lengths of 115-678 kb, and may have resulted from tandem duplication and inversion during primate evolution. The palindromic complex contains 11 families of transcription units, all expressed in testis. Deletions of AZFc that cause infertility are remarkably uniform, spanning a 3.5-Mb segment and bounded by 229-kb direct repeats that probably served as substrates for homologous recombination.

Surveying Saccharomyces genomes to identify functional elements by comparative DNA sequence analysis.

Comparative sequence analysis has facilitated the discovery of protein coding genes and important functional sequences within proteins, but has been less useful for identifying functional sequence elements in nonprotein-coding DNA because the relatively rapid rate of change of nonprotein-coding sequences and the relative simplicity of non-coding regulatory sequence elements necessitates the comparison of sequences of relatively closely related species. We tested the use of comparative DNA sequence analysis to aid identification of promoter regulatory elements, nonprotein-coding RNA genes, and small protein-coding genes by surveying random DNA sequences of several Saccharomyces yeast species, with the goal of learning which species are best suited for comparisons with S. cerevisiae. We also determined the DNA sequence of a few specific promoters and RNA genes of several Saccharomyces species to determine the degree of conservation of known functional elements within the genome. Our results lead us to conclude that comparative DNA sequence analysis will enable identification of functionally conserved elements within the yeast genome, and suggest a path for obtaining this information.

Rapid gene mapping in Caenorhabditis elegans using a high density polymorphism map.

Single nucleotide polymorphisms (SNPs) are valuable genetic markers of human disease. They also comprise the highest potential density marker set available for mapping experimentally derived mutations in model organisms such as Caenorhabditis elegans. To facilitate the positional cloning of mutations we have identified polymorphisms in CB4856, an isolate from a Hawaiian island that shows a uniformly high density of polymorphisms compared with the reference Bristol N2 strain. Based on 5.4 Mbp of aligned sequences, we predicted 6,222 polymorphisms. Furthermore, 3,457 of these markers modify restriction enzyme recognition sites ('snip-SNPs') and are therefore easily detected as RFLPs. Of these, 493 were experimentally confirmed by restriction digest to produce a snip-SNP map of the worm genome. A mapping strategy using snip-SNPs and bulked segregant analysis (BSA) is outlined. CB4856 is crossed into a mutant strain, and exclusion of CB4856 alleles of a subset of snip-SNPs in mutant progeny is assessed with BSA. The proximity of a linked marker to the mutation is estimated by the relative proportion of each form of the biallelic marker in populations of wildtype and mutant genomes. The usefulness of this approach is illustrated by the rapid mapping of the dyf-5 gene.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms.

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

A physical map of the human genome.

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.

A physical map of the human Y chromosome.

The non-recombining region of the human Y chromosome (NRY), which comprises 95% of the chromosome, does not undergo sexual recombination and is present only in males. An understanding of its biological functions has begun to emerge from DNA studies of individuals with partial Y chromosomes, coupled with molecular characterization of genes implicated in gonadal sex reversal, Turner syndrome, graft rejection and spermatogenic failure. But mapping strategies applied successfully elsewhere in the genome have faltered in the NRY, where there is no meiotic recombination map and intrachromosomal repetitive sequences are abundant. Here we report a high-resolution physical map of the euchromatic, centromeric and heterochromatic regions of the NRY and its construction by unusual methods, including genomic clone subtraction and dissection of sequence family variants. Of the map's 758 DNA markers, 136 have multiple locations in the NRY, reflecting its unusually repetitive sequence composition. The markers anchor 1,038 bacterial artificial chromosome clones, 199 of which form a tiling path for sequencing.

Theories and applications for sequencing randomly selected clones.

Theory is developed for the process of sequencing randomly selected large-insert clones. Genome size, library depth, clone size, and clone distribution are considered relevant properties and perfect overlap detection for contig assembly is assumed. Genome-specific and nonrandom effects are neglected. Order of magnitude analysis indicates library depth is of secondary importance compared to the other variables, especially as clone size diminishes. In such cases, the well-known Poisson coverage law is a good approximation. Parameters derived from these models are used to examine performance for the specific case of sequencing random human BAC clones. We compare coverage and redundancy rates for libraries possessing uniform and nonuniform clone distributions. Results are measured against data from map-based human-chromosome-2 sequencing. We conclude that the map-based approach outperforms random clone sequencing, except early in a project. However, simultaneous use of both strategies can be beneficial if a performance-based estimate for halting random clone sequencing is made. Results further show that the random approach yields maximum effectiveness using nonbiased rather than biased libraries.

A region of the myosin rod important for interaction with paramyosin in Caenorhabditis elegans striated muscle.

The precise arrangement of molecules within the thick filament, as well as the mechanisms by which this arrangement is specified, remains unclear. In this article, we have exploited a unique genetic interaction between one isoform of myosin heavy chain (MHC) and paramyosin in Caenorhabditis elegans to probe the molecular interaction between MHC and paramyosin in vivo. Using chimeric myosin constructs, we have defined a 322-residue region of the MHC A rod critical for suppression of the structural and motility defects associated with the unc-15(e73) allele. Chimeric constructs lacking this region of MHC A either fail to suppress, or act as dominant enhancers of, the e73 phenotype. Although the 322-residue region is required for suppression activity, our data suggest that sequences along the length of the rod also play a role in the isoform-specific interaction between MHC A and paramyosin. Our genetic and cell biological analyses of construct behavior suggest that the 322-residue region of MHC A is important for thick filament stability. We present a model in which this region mediates an avid interaction between MHC A and paramyosin in parallel arrangement in formation of the filament arms.

Comparative genomic sequence analysis of the human and mouse cystic fibrosis transmembrane conductance regulator genes.

The identification of the cystic fibrosis transmembrane conductance regulator gene (CFTR) in 1989 represents a landmark accomplishment in human genetics. Since that time, there have been numerous advances in elucidating the function of the encoded protein and the physiological basis of cystic fibrosis. However, numerous areas of cystic fibrosis biology require additional investigation, some of which would be facilitated by information about the long-range sequence context of the CFTR gene. For example, the latter might provide clues about the sequence elements responsible for the temporal and spatial regulation of CFTR expression. We thus sought to establish the sequence of the chromosomal segments encompassing the human CFTR and mouse Cftr genes, with the hope of identifying conserved regions of biologic interest by sequence comparison. Bacterial clone-based physical maps of the relevant human and mouse genomic regions were constructed, and minimally overlapping sets of clones were selected and sequenced, eventually yielding approximately 1.6 Mb and approximately 358 kb of contiguous human and mouse sequence, respectively. These efforts have produced the complete sequence of the approximately 189-kb and approximately 152-kb segments containing the human CFTR and mouse Cftr genes, respectively, as well as significant amounts of flanking DNA. Analyses of the resulting data provide insights about the organization of the CFTR/Cftr genes and potential sequence elements regulating their expression. Furthermore, the generated sequence reveals the precise architecture of genes residing near CFTR/Cftr, including one known gene (WNT2/Wnt2) and two previously unknown genes that immediately flank CFTR/Cftr.

High-throughput plasmid DNA purification for 3 cents per sample.

To accommodate the increasingly rapid rates of DNA sequencing we have developed and implemented an inexpensive, expeditious method for the purification of double-stranded plasmid DNA clones. The robust nature, high throughput, low degree of technical difficulty and extremely low cost have made it the plasmid DNA preparation method of choice in both our expressed sequence tag (EST) and genome sequencing projects. Here we report the details of the method and describe its application in the generation of more than 700 000 ESTs at a rate exceeding 16 000 per week.

Myotactin, a novel hypodermal protein involved in muscle-cell adhesion in Caenorhabditis elegans.

In C. elegans, assembly of hypodermal hemidesmosome-like structures called fibrous organelles is temporally and spatially coordinated with the assembly of the muscle contractile apparatus, suggesting that signals are exchanged between these cell types to position fibrous organelles correctly. Myotactin, a protein recognized by monoclonal antibody MH46, is a candidate for such a signaling molecule. The antigen, although expressed by hypodermis, first reflects the pattern of muscle elements and only later reflects the pattern of fibrous organelles. Confocal microscopy shows that in adult worms myotactin and fibrous organelles show coincident localization. Further, cell ablation studies show the bodywall muscle cells are necessary for normal myotactin distribution. To investigate myotactin's role in muscle-hypodermal signaling, we characterized the myotactin locus molecularly and genetically. Myotactin is a novel transmembrane protein of approximately 500 kd. The extracellular domain contains at least 32 fibronectin type III repeats and the cytoplasmic domain contains unique sequence. In mutants lacking myotactin, muscle cells detach when embryonic muscle contraction begins. Later in development, fibrous organelles become delocalized and are not restricted to regions of the hypodermis previously contacted by muscle. These results suggest myotactin helps maintain the association between the muscle contractile apparatus and hypodermal fibrous organelles.

Representation of cloned genomic sequences in two sequencing vectors: correlation of DNA sequence and subclone distribution.

Representation of subcloned Caenorhabditis elegans and human DNA sequences in both M13 and pUC sequencing vectors was determined in the context of large scale genomic sequencing. In many cases, regions of subclone under-representation correlated with the occurrence of repeat sequences, and in some cases the under-representation was orientation specific. Factors which affected subclone representation included the nature and complexity of the repeat sequence, as well as the length of the repeat region. In some but not all cases, notable differences between the M13 and pUC subclone distributions existed. However, in all regions lacking one type of subclone (either M13 or pUC), an alternate subclone was identified in at least one orientation. This suggests that complementary use of M13 and pUC subclones would provide the most comprehensive subclone coverage of a given genomic sequence.

High throughput fingerprint analysis of large-insert clones.

As part of the Human Genome Project, the Washington University Genome Sequencing Center has commenced systematic sequencing of human chromsome 7. To organize and supply the effort, we have undertaken the construction of sequence-ready physical maps for defined chromosomal intervals. Map construction is a serial process composed of three main activities. First, candidate STS-positive large-insert PAC and BAC clones are identified. Next, these candidate clones are subjected to fingerprint analysis. Finally, the fingerprint data are used to assemble sequence-ready maps. The fingerprinting method we have devised is key to the success of the overall approach. We present here the details of the method and show that the fingerprints are of sufficient quality to permit the construction of megabase-size contigs in defined regions of the human genome. We anticipate that the high throughput and precision characteristic of our fingerprinting method will make it of general utility.