Comparative analyses of multi-species sequences from targeted genomic regions.

The systematic comparison of genomic sequences from different organisms represents a central focus of contemporary genome analysis. Comparative analyses of vertebrate sequences can identify coding and conserved non-coding regions, including regulatory elements, and provide insight into the forces that have rendered modern-day genomes. As a complement to whole-genome sequencing efforts, we are sequencing and comparing targeted genomic regions in multiple, evolutionarily diverse vertebrates. Here we report the generation and analysis of over 12 megabases (Mb) of sequence from 12 species, all derived from the genomic region orthologous to a segment of about 1.8 Mb on human chromosome 7 containing ten genes, including the gene mutated in cystic fibrosis. These sequences show conservation reflecting both functional constraints and the neutral mutational events that shaped this genomic region. In particular, we identify substantial numbers of conserved non-coding segments beyond those previously identified experimentally, most of which are not detectable by pair-wise sequence comparisons alone. Analysis of transposable element insertions highlights the variation in genome dynamics among these species and confirms the placement of rodents as a sister group to the primates.

The NIEHS Xenopus maternal EST project: interim analysis of the first 13,879 ESTs from unfertilized eggs.

The sequencing of expressed sequence tags (ESTs) from Xenopus laevis has lagged behind efforts on many other common experimental organisms and man, partly because of the pseudotetraploid nature of the Xenopus genome. Nonetheless, large collections of Xenopus ESTs would be useful in gene discovery, oligonucleotide-based knockout studies, gene chip analyses of normal and perturbed development, mapping studies in the related diploid frog X. tropicalis, and for other reasons. We have created a normalized library of cDNAs from unfertilized Xenopus eggs. These cells contain all of the information necessary for the first several cell divisions in the early embryo, as well as much of the information needed for embryonic pattern formation and cell fate determination. To date, we have successfully sequenced 13,879 ESTs out of 16,607 attempts (83.6% success rate), with an average sequence read length of 508 bp. Using a fragment assembly program, these ESTs were assembled into 8,985 'contigs' comprised of up to 11 ESTs each. When these contigs were used to search publicly available databases, 46.2% bore no relationship to protein or DNA sequences in the database at the significance level of 1e-6. Examination of a sample of 100 of the assembled contigs revealed that most ( approximately 87%) were comprised of two apparent allelic variants. Expression profiles of 16 of the most prominent contigs showed that 12 exhibited some degree of zygotic expression. These findings have implications for sequence-specific applications for Xenopus ESTs, particularly the use of allele-specific oligonucleotides for knockout studies, differential hybridization techniques such as gene chip analysis, and the establishment of accurate nomenclature and databases for this species.

Gene discovery using computational and microarray analysis of transcription in the Drosophila melanogaster testis.

Identification and annotation of all the genes in the sequenced Drosophila genome is a work in progress. Wild-type testis function requires many genes and is thus of potentially high value for the identification of transcription units. We therefore undertook a survey of the repertoire of genes expressed in the Drosophila testis by computational and microarray analysis. We generated 3141 high-quality testis expressed sequence tags (ESTs). Testis ESTs computationally collapsed into 1560 cDNA set used for further analysis. Of those, 11% correspond to named genes, and 33% provide biological evidence for a predicted gene. A surprising 47% fail to align with existing ESTs and 16% with predicted genes in the current genome release. EST frequency and microarray expression profiles indicate that the testis mRNA population is highly complex and shows an extended range of transcript abundance. Furthermore, >80% of the genes expressed in the testis showed onefold overexpression relative to ovaries, or gonadectomized flies. Additionally, >3% showed more than threefold overexpression at p <0.05. Surprisingly, 22% of the genes most highly overexpressed in testis match Drosophila genomic sequence, but not predicted genes. These data strongly support the idea that sequencing additional cDNA libraries from defined tissues, such as testis, will be important tools for refined annotation of the Drosophila genome. Additionally, these data suggest that the number of genes in Drosophila will significantly exceed the conservative estimate of 13,601.

Mutation of a gene encoding a putative chaperonin causes McKusick-Kaufman syndrome.

McKusick-Kaufman syndrome (MKKS, MIM 236700) is a human developmental anomaly syndrome comprising hydrometrocolpos (HMC), postaxial polydactyly (PAP) and congenital heart disease (CHD). MKKS has been mapped in the Old Order Amish population to 20p12, between D20S162 and D20S894 (ref. 3). Here we describe the identification of a gene mutated in MKKS. We analysed the approximately 450-kb candidate region by sample sequencing, which revealed the presence of several known genes and EST clusters. We evaluated candidate transcripts by northern-blot analysis of adult and fetal tissues. We selected one transcript with widespread expression, MKKS, for analysis in a patient from the Amish pedigree and a sporadic, non-Amish case. The Old Order Amish patient was found to be homozygous for an allele that had two missense substitutions and the non-Amish patient was a compound heterozygote for a frameshift mutation predicting premature protein truncation and a distinct missense mutation. The MKKS predicted protein shows amino acid similarity to the chaperonin family of proteins, suggesting a role for protein processing in limb, cardiac and reproductive system development. We believe that this is the first description of a human disorder caused by mutations affecting a putative chaperonin molecule.

The genomic region encompassing the nephropathic cystinosis gene (CTNS): complete sequencing of a 200-kb segment and discovery of a novel gene within the common cystinosis-causing deletion.

Nephropathic cystinosis is an autosomal recessive disorder caused by the defective transport of cystine out of lysosomes. Recently, the causative gene (CTNS) was identified and presumed to encode an integral membrane protein called cystinosin. Many of the disease-associated mutations in CTNS are deletions, including one >55 kb in size that represents the most common cystinosis allele encountered to date. In an effort to determine the precise genomic organization of CTNS and to gain sequence-based insight about the DNA within and flanking cystinosis-associated deletions, we mapped and sequenced the region of human chromosome 17p13 encompassing CTNS. Specifically, a bacterial artificial chromosome (BAC)-based physical map spanning CTNS was constructed by sequence-tagged site (STS)-content mapping. The resulting BAC contig provided the relative order of 43 STSs. Two overlapping BACs, which together contain all of the CTNS exons as well as extensive amounts of flanking DNA, were selected and subjected to shotgun sequencing. A total of 200,237 bp of contiguous, high-accuracy sequence was generated. Analysis of the resulting data revealed a number of interesting features about this genomic region, including the long-range organization of CTNS, insight about the breakpoints and intervening DNA associated with the common cystinosis-causing deletion, and structural information about five genes neighboring CTNS (human ortholog of rat vanilloid receptor subtype 1 gene, CARKL, TIP-1, P2X5, and HUMINAE). In particular, sequence analysis detected the presence of a novel gene (CARKL) residing within the most common cystinosis-causing deletion. This gene encodes a previously unknown protein that is predicted to function as a carbohydrate kinase. Interestingly, both CTNS and CARKL are absent in nearly half of all cystinosis patients (i.e., those homozygous for the common deletion). [The sequence data described in this paper have been submitted to the GenBank data library under accession nos. AF168787 and AF163573.]

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.

Gene expression in proliferating human erythroid cells.

A complete understanding of human erythropoiesis will require a robust description of transcriptional activity in hematopoietic cells that proliferate and differentiate in response to erythropoietin (EPO). For this purpose, we cultured peripheral blood mononuclear cells in the presence or in the absence of EPO and examined the transcriptional profile of those cells arising only in response to EPO. A distinct population of CD71( +) cells that demonstrated an average of six additional doublings in suspension culture and erythroid colony formation in methylcellulose was isolated. Suppression subtractive hybridization of mRNA isolated from those cells permitted the identification of transcribed genes. A summary of 719 expressed sequence tags (ESTs) describing 505 independent transcripts is provided here with a full analysis of each EST available at http://hembase.niddk.nih.gov. Several transcripts that matched genes previously reported in the context of erythroid differentiation including 4 cell surface proteins were expressed at this developmental stage. Active chromatin remodeling was suggested by the identification of 4 histone proteins, 4 high-mobility group proteins, 13 transcription factors, and 6 genes involved in DNA recombination and repair. Numerous genes associated with leukemic translocations were also recognized including topoisomerases I and II, nucleophosmin, Translin, EGR1, dek, pim-1, TFG, and MLL. In addition to known transcripts, 44 novel EST were discovered. This transcriptional profile provides the first genomic-scale description of gene activity in erythroid progenitor cells.

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.

2006 expressed-sequence tags derived from human chromosome 7-enriched cDNA libraries.

The establishment and mapping of gene-specific DNA sequences greatly complement the ongoing efforts to map and sequence all human chromosomes. To facilitate our studies of human chromosome 7, we have generated and analyzed 2006 expressed-sequence tags (ESTs) derived from a collection of direct selection cDNA libraries that are highly enriched for human chromosome 7 gene sequences. Similarity searches indicate that approximately two-thirds of the ESTs are not represented by sequences in the public databases, including those in dbEST. In addition, a large fraction (68%) of the ESTs do not have redundant or overlapping sequences within our collection. Human DNA-specific sequence-tagged sites (STSs) have been developed from 190 of the ESTs. Remarkably, 180 (96%) of these STSs map to chromosome 7, demonstrating the robustness of chromosome enrichment in constructing the direct selection cDNA libraries. Thus far, 140 of these EST-specific STSs have been assigned unequivocally to YAC contigs that are distributed across the chromosome. Together, these studies provide > 2000 ESTs highly enriched for chromosome 7 gene sequences, 180 new chromosome 7 STSs corresponding to ESTs, and a definitive demonstration of the ability to enrich for chromosome-specific cDNAs by direct selection. Furthermore, the libraries, sequence data, and mapping information will contribute to the construction of a chromosome 7 transcript map.

A collection of 1814 human chromosome 7-specific STSs.

An established goal of the ongoing Human Genome Project is the development and mapping of sequence-tagged sites (STSs) every 100 kb, on average, across all human chromosomes. En route to constructing such a physical map of human chromosome 7, we have generated 1814 chromosome 7-specific STSs. The corresponding PCR assays were designed by the use of DNA sequence determined in our laboratory (79%) or generated elsewhere (21%) and were demonstrated to be suitable for screening yeast artificial chromosome (YAC) libraries. This collection provides the requisite landmarks for constructing a physical map of chromosome 7 at < 100-kb average spacing of STSs.

Mapping the RP10 locus for autosomal dominant retinitis pigmentosa on 7q: refined genetic positioning and localization within a well-defined YAC contig.

Retinitis pigmentosa is a genetically heterogeneous disease that has autosomal dominant, autosomal recessive and X-linked forms. Autosomal dominant retinitis pigmentosa (adRP) has thus far been associated with eight distinct loci, including the rhodopsin and peripherin/RDS genes as well as unidentified genes on chromosomes 7p, 7q, 8q, 17p, 17q, and 19q. The RP10 locus for adRP on chromosome 7q was first mapped in a Spanish family; later, an unrelated American family was identified that also showed linkage to 7q. By combining the linkage results from both families, we are able to assign the disease gene to a 5-cM interval on 7q. Based on extensive physical mapping of this region, the genetic interval is now fully contained within a approximately 5-Mb segment on a well-defined YAC contig. These studies significantly reduce the size of the RP10 critical region, exclude a number of possible candidate genes, and provide the necessary cloned DNA for the positional cloning of the RP10 gene.

Dependence of lactose metabolism upon mutarotase encoded in the gal operon in Escherichia coli.

A new gene (galM) has been identified as the fourth cistron of the gal operon, encoding enzymes for the metabolism of galactose and lactose in Escherichia coli. Induction of the gal operon either from the gal promoters or from a neighboring prophage lambda promoter expresses the galM gene as well. The new structure of the gal operon from the promoter end is galE-galT-galK-galM in counter-clockwise orientation on the chromosome. Genetic and biochemical analyses have revealed that the galM gene product has mutarotase activity, which converts alpha-aldose to the beta-anomer. Unlike mutarotase from other bacteria in which the enzyme is primarily processed for export and secretion, the mutarotase from E. coli does not appear to be processed and yet is still found in periplasm (and culture media when overexpressed) in significant amounts. Although the interconversion of the sugar anomers occurs spontaneously in pure water in vitro, the in vivo formation of alpha-D-galactopyranose (the substrate for phosphorylation) from beta-D-galactopyranose (generated by beta-galactosidase hydrolysis of lactose) is largely dependent upon the presence of the mutarotase. This shows that efficient lactose metabolism requires mutarotase. These results give credence to the idea that the activity of intracellular water is not high enough to permit a simple extrapolation of observed in vitro reactions to in vivo situations in every case.