Transposon Tc1-derived, sequence-tagged sites in Caenorhabditis elegans as markers for gene mapping

We present an approach to map large numbers of Tc1 transposon insertions in the genome of Caenorhabditis elegans. Strains have been described that contain up to 500 polymorphic Tc1 insertions. From these we have cloned and shotgun sequenced over 2000 Tc1 flanks, resulting in an estimated set of 400 or more distinct Tc1 insertion alleles. Alignment of these sequences revealed a weak Tc1 insertion site consensus sequence that was symmetric around the invariant TA target site and reads CAYATATRTG. The Tc1 flanking sequences were compared with 40 Mbp of a C. elegans genome sequence. We found 151 insertions within the sequenced area, a density of ≈1 Tc1 insertion in every 265 kb. As the rest of the C. elegans genome sequence is obtained, remaining Tc1 alleles will fall into place. These mapped Tc1 insertions can serve two functions: (i) insertions in or near genes can be used to isolate deletion derivatives that have that gene mutated; and (ii) they represent a dense collection of polymorphic sequence-tagged sites. We demonstrate a strategy to use these Tc1 sequence-tagged sites in fine-mapping mutations.

Sequence-tagged microsatellite profiling (STMP): a rapid technique for developing SSR markers

We describe a technique, sequence-tagged microsatellite profiling (STMP), to rapidly generate large numbers of simple sequence repeat (SSR) markers from genomic or cDNA. This technique eliminates the need for library screening to identify SSR-containing clones and provides an ∼25-fold increase in sequencing throughput compared to traditional methods. STMP generates short but characteristic nucleotide sequence tags for fragments that are present within a pool of SSR amplicons. These tags are then ligated together to form concatemers for cloning and sequencing. The analysis of thousands of tags gives rise to a representational profile of the abundance and frequency of SSRs within the DNA pool, from which low copy sequences can be identified. As each tag contains sufficient nucleotide sequence for primer design, their conversion into PCR primers allows the amplification of corresponding full-length fragments from the pool of SSR amplicons. These fragments permit the full characterisation of a SSR locus and provide flanking sequence for the development of a microsatellite marker. Alternatively, sequence tag primers can be used to directly amplify corresponding SSR loci from genomic DNA, thereby reducing the cost of developing a microsatellite marker to the synthesis of just one sequence-specific primer. We demonstrate the utility of STMP by the development of SSR markers in bread wheat.

A high-resolution annotated physical map of the human chromosome 13q12-13 region containing the breast cancer susceptibility locus BRCA2.

Various types of physical mapping data were assembled by developing a set of computer programs (Integrated Mapping Package) to derive a detailed, annotated map of a 4-Mb region of human chromosome 13 that includes the BRCA2 locus. The final assembly consists of a yeast artificial chromosome (YAC) contig with 42 members spanning the 13q12-13 region and aligned contigs of 399 cosmids established by cross-hybridization between the cosmids, which were selected from a chromosome 13-specific cosmid library using inter-Alu PCR probes from the YACs. The end sequences of 60 cosmids spaced nearly evenly across the map were used to generate sequence-tagged sites (STSs), which were mapped to the YACs by PCR. A contig framework was generated by STS content mapping, and the map was assembled on this scaffold. Additional annotation was provided by 72 expressed sequences and 10 genetic markers that were positioned on the map by hybridization to cosmids.

A yeast artificial chromosome-based map of the region of chromosome 20 containing the diabetes-susceptibility gene, MODY1, and a myeloid leukemia related gene.

We have generated a physical map of human chromosome bands 20q11.2-20q13.1, a region containing a gene involved in the development of one form of early-onset, non-insulin-dependent diabetes mellitus, MODY1, as well as a putative myeloid tumor suppressor gene. The yeast artificial chromosome contig consists of 71 clones onto which 71 markers, including 20 genes, 5 expressed sequence tags, 32 simple tandem repeat DNA polymorphisms, and 14 sequence-tagged sites have been ordered. This region spans about 18 Mb, which represents about 40% of the physical length of 20q. Using this physical map, we have refined the location of MODY1 to a 13-centimorgan interval (approximately equal to 7 Mb) between D20S169 and D20S176. The myeloid tumor suppressor gene was localized to an 18-centimorgan interval (approximately equal to 13 Mb) between RPN2 and D20S17. This physical map will facilitate the isolation of MODY1 and the myeloid tumor suppressor gene.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Functional copies of a human gene can be directly isolated by transformation-associated recombination cloning with a small 3′ end target sequence

Unique, small sequences (sequence tag sites) have been identified at the 3′ ends of most human genes that serve as landmarks in genome mapping. We investigated whether a single copy gene could be isolated directly from total human DNA by transformation-associated recombination (TAR) cloning in yeast using a short, 3′ unique target. A TAR cloning vector was constructed that, when linearized, contained a small amount (381 bp) of 3′ hypoxanthine phosphoribosyltransferase (HPRT) sequence at one end and an 189-bp Alu repeat at the other end. Transformation with this vector along with human DNA led to selective isolations of the entire HPRT gene as yeast artificial chromosomes (YACs) that extended from the 3′ end sequence to various Alu positions as much as 600 kb upstream. These YACs were retrofitted with a NeoR and a bacterial artificial chromosome (BAC) sequence to transfer the YACs to bacteria and subsequently the BACs to mouse cells by using a Neo selection. Most of the HPRT isolates were functional, demonstrating that TAR cloning retains the functional integrity of the isolated material. Thus, this modified version of TAR cloning, which we refer to as radial TAR cloning, can be used to isolate large segments of the human genome accurately and directly with only a small amount of sequence information.

A high-resolution physical map of human chromosome 11.

The development of a highly reliable physical map with landmark sites spaced an average of 100 kbp apart has been a central goal of the Human Genome Project. We have approached the physical mapping of human chromosome 11 with this goal as a primary target. We have focused on strategies that would utilize yeast artificial chromosome (YAC) technology, thus permitting long-range coverage of hundreds of kilobases of genomic DNA, yet we sought to minimize the ambiguities inherent in the use of this technology, particularly the occurrence of chimeric genomic DNA clones. This was achieved through the development of a chromosome 11-specific YAC library from a human somatic cell hybrid line that has retained chromosome 11 as its sole human component.To maximize the efficiency of YAC contig assembly and extension, we have employed an Alu-PCR-based hybridization screening system. This system eliminates many of the more costly and time-consuming steps associated with sequence tagged site content mapping such as sequencing, primer production, and hierarchical screening, resulting in greater efficiency with increased throughput and reduced cost. Using these approaches, we have achieved YAC coverage for >90% of human chromosome 11, with an average intermarker distance of <100 kbp. Cytogenetic localization has been determined for each contig by fluorescent in situ hybridization and/or sequence tagged site content. The YAC contigs that we have generated should provide a robust framework to move forward to sequence-ready templates for the sequencing efforts of the Human Genome Project as well as more focused positional cloning on chromosome 11.

GenMapDB: a database of mapped human BAC clones

GenMapDB (http://genomics.med.upenn.edu/genmapdb) is a repository of human bacterial artificial chromosome (BAC) clones mapped by our laboratory to sequence-tagged site markers. Currently, GenMapDB contains over 3000 mapped clones that span 19 chromosomes, chromosomes 2, 4, 5, 9–22, X and Y. This database provides positional information about human BAC clones from the RPCI-11 human male BAC library. It also contains restriction fragment analysis data and end sequences of the clones. GenMapDB is freely available to the public. The main purpose of GenMapDB is to organize the mapping data and to allow the research community to search for mapped BAC clones that can be used in gene mapping studies and chromosomal mutation analysis projects.

Identification of urocortin III, an additional member of the corticotropin-releasing factor (CRF) family with high affinity for the CRF2 receptor

The corticotropin-releasing factor (CRF) family of neuropeptides includes the mammalian peptides CRF, urocortin, and urocortin II, as well as piscine urotensin I and frog sauvagine. The mammalian peptides signal through two G protein-coupled receptor types to modulate endocrine, autonomic, and behavioral responses to stress, as well as a range of peripheral (cardiovascular, gastrointestinal, and immune) activities. The three previously known ligands are differentially distributed anatomically and have distinct specificities for the two major receptor types. Here we describe the characterization of an additional CRF-related peptide, urocortin III, in the human and mouse. In searching the public human genome databases we found a partial expressed sequence tagged (EST) clone with significant sequence identity to mammalian and fish urocortin-related peptides. By using primers based on the human EST sequence, a full-length human clone was isolated from genomic DNA that encodes a protein that includes a predicted putative 38-aa peptide structurally related to other known family members. With a human probe, we then cloned the mouse ortholog from a genomic library. Human and mouse urocortin III share 90% identity in the 38-aa putative mature peptide. In the peptide coding region, both human and mouse urocortin III are 76% identical to pufferfish urocortin-related peptide and more distantly related to urocortin II, CRF, and urocortin from other mammalian species. Mouse urocortin III mRNA expression is found in areas of the brain including the hypothalamus, amygdala, and brainstem, but is not evident in the cerebellum, pituitary, or cerebral cortex; it is also expressed peripherally in small intestine and skin. Urocortin III is selective for type 2 CRF receptors and thus represents another potential endogenous ligand for these receptors.

Structure, tissue distribution, and chromosomal localization of the prepronociceptin gene.

Nociceptin (orphanin FQ), the newly discovered natural agonist of opioid receptor-like (ORL1) receptor, is a neuropeptide that is endowed with pronociceptive activity in vivo. Nociceptin is derived from a larger precursor, prepronociceptin (PPNOC), whose human, mouse, and rat genes we have now isolated. The PPNOC gene is highly conserved in the three species and displays organizational features that are strikingly similar to those of the genes of preproenkephalin, preprodynorphin, and preproopiomelanocortin, the precursors to endogenous opioid peptides, suggesting the four genes belong to the same family-i.e., have a common evolutionary origin. The PPNOC gene encodes a single copy of nociceptin as well as of other peptides whose sequence is strictly conserved across murine and human species; hence it is likely to be neurophysiologically significant. Northern blot analysis shows that the PPNOC gene is predominantly transcribed in the central nervous system (brain and spinal cord) and, albeit weakly, in the ovary, the sole peripheral organ expressing the gene. By using a radiation hybrid cell line panel, the PPNOC gene was mapped to the short arm of human chromosome 8 (8p21), between sequence-tagged site markers WI-5833 and WI-1172, in close proximity of the locus encoding the neurofilament light chain NEFL. Analysis of yeast artificial chromosome clones belonging to the WC8.4 contig covering the 8p21 region did not allow to detect the presence of the gene on these yeast artificial chromosomes, suggesting a gap in the coverage within this contig.

Mini-chromosomes derived from the human Y chromosome by telomere directed chromosome breakage.

We have used telomeric DNA to break two acrocentric derivatives of the human Y chromosome into mini-chromosomes that are small enough to be size- fractionated by pulsed-field gel electrophoresis. One of the mini-chromosomes is about 7 Mb in size and sequence-tagged site analysis of this molecule suggests that it corresponds to a simple truncation of the short arm of the Y chromosome. Five of the mini-chromosomes are derived from the long arm, are all rearranged by more than a simple truncation, and range in size from 4.0 Mb to 9 Mb. We have studied the mitotic stabilities of these mini-chromosomes and shown that they are stably maintained by cells proliferating in culture for about 100 cell divisions.

Genomic cloning of methylthioadenosine phosphorylase: a purine metabolic enzyme deficient in multiple different cancers.

5'-Deoxy-5'-methylthioadenosine phosphorylase (methylthioadeno-sine: ortho-phosphate methylthioribosyltransferase, EC 24.2.28; MTAP) plays a role in purine and polyamine metabolism and in the regulation of transmethylation reactions. MTAP is abundant in normal cells but is deficient in many cancers. Recently, the genes for the cyclin-dependent kinase inhibitors p16 and p15 have been localized to the short arm of human chromosome 9 at band p21, where MTAP and interferon alpha genes (IFNA) also map. Homozygous deletions of p16 and p15 are frequent malignant cell lines. However, the order of the MTAP, p16, p15, and IFNA genes on chromosome 9p is uncertain, and the molecular basis for MTAP deficiency in cancer is unknown. We have cloned the MTAP gene, and have constructed a topologic map of the 9p21 region using yeast artificial chromosome clones, pulse-field gel electrophoresis, and sequence-tagged-site PCR. The MTAP gene consists of eight exons and seven introns. Of 23 malignant cell lines deficient in MTAP protein, all but one had complete or partial deletions. Partial or total deletions of the MTAP gene were found in primary T-cell acute lymphoblastic leukemias (T-ALL). A deletion breakpoint of partial deletions found in cell lines and primary T-ALL was in intron 4. Starting from the centromeric end, the gene order on chromosome 9p2l is p15, p16, MTAP, IFNA, and interferon beta gene (IFNB). These results indicate that MTAP deficiency in cancer is primarily due to codeletion of the MTAP and p16 genes.

Short DNA Fragments without Sequence Similarity Are Initiation Sites for Replication in the Chromosome of the Yeast Yarrowia lipolytica

We have previously shown that both a centromere (CEN) and a replication origin are necessary for plasmid maintenance in the yeast Yarrowia lipolytica (Vernis et al., 1997). Because of this requirement, only a small number of centromere-proximal replication origins have been isolated from Yarrowia. We used a CEN-based plasmid to obtain noncentromeric origins, and several new fragments, some unique and some repetitive sequences, were isolated. Some of them were analyzed by two-dimensional gel electrophoresis and correspond to actual sites of initiation (ORI) on the chromosome. We observed that a 125-bp fragment is sufficient for a functional ORI on plasmid, and that chromosomal origins moved to ectopic sites on the chromosome continue to act as initiation sites. These Yarrowia origins share an 8-bp motif, which is not essential for origin function on plasmids. The Yarrowia origins do not display any obvious common structural features, like bent DNA or DNA unwinding elements, generally present at or near eukaryotic replication origins. Y. lipolytica origins thus share features of those in the unicellular Saccharomyces cerevisiae and in multicellular eukaryotes: they are discrete and short genetic elements without sequence similarity.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis.

Competing hypotheses seek to explain the evolution of oxygenic and anoxygenic processes of photosynthesis. Since chlorophyll is less reduced and precedes bacteriochlorophyll on the modern biosynthetic pathway, it has been proposed that chlorophyll preceded bacteriochlorophyll in its evolution. However, recent analyses of nucleotide sequences that encode chlorophyll and bacteriochlorophyll biosynthetic enzymes appear to provide support for an alternative hypothesis. This is that the evolution of bacteriochlorophyll occurred earlier than the evolution of chlorophyll. Here we demonstrate that the presence of invariant sites in sequence datasets leads to inconsistency in tree building (including maximum-likelihood methods). Homologous sequences with different biological functions often share invariant sites at the same nucleotide positions. However, different constraints can also result in additional invariant sites unique to the genes, which have specific and different biological functions. Consequently, the distribution of these sites can be uneven between the different types of homologous genes. The presence of invariant sites, shared by related biosynthetic genes as well as those unique to only some of these genes, has misled the recent evolutionary analysis of oxygenic and anoxygenic photosynthetic pigments. We evaluate an alternative scheme for the evolution of chlorophyll and bacteriochlorophyll.