3,400 new expressed sequence tags identify diversity of transcripts in human brain.

We present the results of the partial sequencing of over 3,400 expressed sequence tags (ESTs) from human brain cDNA clones, which increases the number of distinct genes expressed in the brain, that are represented by ESTs, to about 6,000. By choosing clones in an unbiased manner, it is possible to construct a profile of the transcriptional activity of the brain at different stages. Proteins that comprise the cytoskeleton are the most abundant; however, a large variety of regulatory proteins are also seen. About half of the ESTs predicted to contain a protein-coding region have no matches in the public peptide databases and may represent new gene families.

Rapid cDNA sequencing (expressed sequence tags) from a directionally cloned human infant brain cDNA library.

A human infant brain cDNA library, made specifically for production of expressed sequence tags (ESTs) was evaluated by partial sequencing of over 1,600 clones. Advantages of this library, constructed for EST sequencing, include the use of directional cloning, size selection, very low numbers of mitochondrial and ribosomal transcripts, short polyA tails, few non-recombinants and a broad representation of transcripts. 37% of the clones were identified, based on matches to over 320 different genes in the public databases. Of these, two proteins similar to the Alzheimer's disease amyloid precursor protein were identified.

Caenorhabditis elegans expressed sequence tags identify gene families and potential disease gene homologues.

A database containing mapped partial cDNA sequences from Caenorhabditis elegans will provide a ready starting point for identifying nematode homologues of important human genes and determining their functions in C. elegans. A total of 720 expressed sequence tags (ESTs) have been generated from 585 clones randomly selected from a mixed-stage C. elegans cDNA library. Comparison of these ESTs with sequence databases identified 422 new C. elegans genes, of which 317 are not similar to any sequences in the database. Twenty-six new genes have been mapped by YAC clone hybridization. Members of several gene families, including cuticle collagens, GTP-binding proteins, and RNA helicases were discovered. Many of the new genes are similar to known or potential human disease genes, including CFTR and the LDL receptor.

Biochemistry and regulation of pre-mRNA splicing.

During the past year, significant advances have been made in the field of pre-mRNA splicing. It is now clear that members of the serine-arginine-rich protein family are key players in exon definition and function at multiple steps in the spliceosome cycle. Novel findings have been made concerning the role of exon sequences, which function as both constitutive and regulated enhancers of splicing, in trans-splicing and as targets for tissue-specific control of splicing patterns. By combining biochemical approaches in human and yeast extracts with genetic analysis, much has been learned about the RNA-RNA and RNA-protein interactions that are necessary to assemble the various complexes that are found along the pathway to the catalytically active spliceosome.

Mutation of a mutL homolog in hereditary colon cancer.

Some cases of hereditary nonpolyposis colorectal cancer (HNPCC) are due to alterations in a mutS-related mismatch repair gene. A search of a large database of expressed sequence tags derived from random complementary DNA clones revealed three additional human mismatch repair genes, all related to the bacterial mutL gene. One of these genes (hMLH1) resides on chromosome 3p21, within 1 centimorgan of markers previously linked to cancer susceptibility in HNPCC kindreds. Mutations of hMLH1 that would disrupt the gene product were identified in such kindreds, demonstrating that this gene is responsible for the disease. These results suggest that defects in any of several mismatch repair genes can cause HNPCC.

Complementary DNA sequencing: expressed sequence tags and human genome project.

Automated partial DNA sequencing was conducted on more than 600 randomly selected human brain complementary DNA (cDNA) clones to generate expressed sequence tags (ESTs). ESTs have applications in the discovery of new human genes, mapping of the human genome, and identification of coding regions in genomic sequences. Of the sequences generated, 337 represent new genes, including 48 with significant similarity to genes from other organisms, such as a yeast RNA polymerase II subunit; Drosophila kinesin, Notch, and Enhancer of split; and a murine tyrosine kinase receptor. Forty-six ESTs were mapped to chromosomes after amplification by the polymerase chain reaction. This fast approach to cDNA characterization will facilitate the tagging of most human genes in a few years at a fraction of the cost of complete genomic sequencing, provide new genetic markers, and serve as a resource in diverse biological research fields.

Whole-genome random sequencing and assembly of Haemophilus influenzae Rd.

An approach for genome analysis based on sequencing and assembly of unselected pieces of DNA from the whole chromosome has been applied to obtain the complete nucleotide sequence (1,830,137 base pairs) of the genome from the bacterium Haemophilus influenzae Rd. This approach eliminates the need for initial mapping efforts and is therefore applicable to the vast array of microbial species for which genome maps are unavailable. The H. influenzae Rd genome sequence (Genome Sequence DataBase accession number L42023) represents the only complete genome sequence from a free-living organism.

Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii.

The complete 1.66-megabase pair genome sequence of an autotrophic archaeon, Methanococcus jannaschii, and its 58- and 16-kilobase pair extrachromosomal elements have been determined by whole-genome random sequencing. A total of 1738 predicted protein-coding genes were identified; however, only a minority of these (38 percent) could be assigned a putative cellular role with high confidence. Although the majority of genes related to energy production, cell division, and metabolism in M. jannaschii are most similar to those found in Bacteria, most of the genes involved in transcription, translation, and replication in M. jannaschii are more similar to those found in Eukaryotes.

Chromosome 2 sequence of the human malaria parasite Plasmodium falciparum.

Chromosome 2 of Plasmodium falciparum was sequenced; this sequence contains 947,103 base pairs and encodes 210 predicted genes. In comparison with the Saccharomyces cerevisiae genome, chromosome 2 has a lower gene density, introns are more frequent, and proteins are markedly enriched in nonglobular domains. A family of surface proteins, rifins, that may play a role in antigenic variation was identified. The complete sequencing of chromosome 2 has shown that sequencing of the A+T-rich P. falciparum genome is technically feasible.

The minimal gene complement of Mycoplasma genitalium.

The complete nucleotide sequence (580,070 base pairs) of the Mycoplasma genitalium genome, the smallest known genome of any free-living organism, has been determined by whole-genome random sequencing and assembly. A total of only 470 predicted coding regions were identified that include genes required for DNA replication, transcription and translation, DNA repair, cellular transport, and energy metabolism. Comparison of this genome to that of Haemophilus influenzae suggests that differences in genome content are reflected as profound differences in physiology and metabolic capacity between these two organisms.

A whole-genome assembly of Drosophila.

We report on the quality of a whole-genome assembly of Drosophila melanogaster and the nature of the computer algorithms that accomplished it. Three independent external data sources essentially agree with and support the assembly's sequence and ordering of contigs across the euchromatic portion of the genome. In addition, there are isolated contigs that we believe represent nonrepetitive pockets within the heterochromatin of the centromeres. Comparison with a previously sequenced 2.9- megabase region indicates that sequencing accuracy within nonrepetitive segments is greater than 99. 99% without manual curation. As such, this initial reconstruction of the Drosophila sequence should be of substantial value to the scientific community.

The sequence of the human genome.

A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.

The genome sequence of Drosophila melanogaster.

The fly Drosophila melanogaster is one of the most intensively studied organisms in biology and serves as a model system for the investigation of many developmental and cellular processes common to higher eukaryotes, including humans. We have determined the nucleotide sequence of nearly all of the approximately 120-megabase euchromatic portion of the Drosophila genome using a whole-genome shotgun sequencing strategy supported by extensive clone-based sequence and a high-quality bacterial artificial chromosome physical map. Efforts are under way to close the remaining gaps; however, the sequence is of sufficient accuracy and contiguity to be declared substantially complete and to support an initial analysis of genome structure and preliminary gene annotation and interpretation. The genome encodes approximately 13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

RNA binding activity of heterodimeric splicing factor U2AF: at least one RS domain is required for high-affinity binding.

The pre-mRNA splicing factor U2AF (U2 small nuclear ribonucleoprotein particle [snRNP] auxiliary factor) plays a critical role in 3' splice site selection. U2AF binds site specifically to the intron pyrimidine tract between the branchpoint and the 3' splice site and targets U2 snRNP to the branch site at an early step in spliceosome assembly. Human U2AF is a heterodimer composed of large (hU2AF65) and small (hU2AF35) subunits. hU2AF65 contains an arginine-serine-rich (RS) domain and three RNA recognition motifs (RRMs). hU2AF35 has a degenerate RRM and a carboxyl-terminal RS domain. Genetic studies have recently shown that the RS domains on the Drosophila U2AF subunit homologs are each inessential and might have redundant functions in vivo. The site-specific pyrimidine tract binding activity of the U2AF heterodimer has previously been assigned to hU2AF65. While the requirement for the three RRMs on hU2AF65 is firmly established, a role for the large-subunit RS domain in RNA binding remains unresolved. We have analyzed the RNA binding activity of the U2AF heterodimer in vitro. When the Drosophila small-subunit homolog (dU2AF38) was complexed with the large-subunit (dU2AF50) pyrimidine tract, RNA binding activity increased 20-fold over that of free dU2AF50. We detected a similar increase in RNA binding activity when we compared the human U2AF heterodimer and hU2AF65. Surprisingly, the RS domain on dU2AF38 was necessary for the increased binding activity of the dU2AF heterodimer. In addition, removal of the RS domain from the Drosophila large-subunit monomer (dU2AF50DeltaRS) severely impaired its binding activity. However, if the dU2AF38 RS domain was supplied in a complex with dU2AF50DeltaRS, high-affinity binding was restored. These results suggest that the presence of one RS domain of U2AF, on either the large or small subunit, promotes high-affinity pyrimidine tract RNA binding activity, consistent with redundant roles for the U2AF RS domains in vivo.

A 500-kb region on chromosome 16p13.1 contains the pseudoxanthoma elasticum locus: high-resolution mapping and genomic structure.

We have recently mapped the genetic defect underlying pseudoxanthoma elasticum (PXE), an inherited disorder characterized by progressive calcification of elastic fibers in skin, eye, and cardiovascular system, to chromosome 16p 13.1. Here we report further data on the fine-mapping and genomic structure of this locus. Haplotype analysis of informative PXE families narrowed the locus to an interval of less than 500 kb located between markers D16B9621 and D16S764. Three overlapping YAC clones were found to cover this region through YAC-STS content mapping. An overlapping BAC contig was then constructed to cover this interval and the surrounding region. About 80% of this chromosomal region has been fully sequenced using the BAC shotgun technique. Gene content and sequence analysis predicted four genes (MRP1, MRP6, PM5, and a novel transcript) and two pseudogenes (ARA and PKDI) within this interval. By screening a somatic cell hybrid panel we were able to precision-map the breakpoint of Cy185 and the starting point of a chromosomal duplication within 20 kb of BAC A962B4. The present data further refine the localization of PXE, provide additional physical cloning resources, and will aid in the eventual identification of the genetic defect causing PXE.

Identification of new Schistosoma mansoni genes by the EST strategy using a directional cDNA library.

A directional size-selected cDNA library constructed from Schistosoma mansoni (Sm) adult worm RNA was used for the generation of expressed sequence tags (EST). From one or both ends of 429 distinct cDNA clones 607 EST were obtained. Of these, only 16% were previously known Sm genes. More than 22% of the clones had matches with entries for other organisms in the databases. These new Sm genes constituted a broad range of transcripts distributed among cytoplasmic structural and regulatory proteins, enzymes, membrane, nuclear and secretory proteins, and proteins with other functions. Almost 33% of the clones had no significant database matches and thus potentially represent Sm-specific genes. Among the latter, several clones, as judged by their redundancy in the library, appear to represent abundant transcripts. The data, taken as a whole, more than double the number of Sm genes identified by nucleotide sequencing and indicate the potential value of the adoption of genome sequencing strategies for the rapid increase in knowledge of complex disease-causing organisms.

The malaria genome sequencing project.

An international consortium of genome centres, advanced development teams and funding agencies has begun the task of sequencing the genome of the parasite Plasmodium falciparum, the most important cause of human malaria. Sequencing is proceeding chromosome by chromosome, and the annotated sequence of chromosome 2 is nearly finished. With the continual release of sequence data as they are generated, malaria researchers have access to a steady stream of genomic sequences and will soon have the complete annotation of all of the estimated 5000-7000 P. falciparum genes. The task will then be how to best apply these data to the development of new anti-malarial drugs, vaccines and diagnostic tests. This review provides a brief overview of the Malaria Genome Sequencing Project and suggests potential directions for future malaria research.

Optimization of Applied Biosystems 373 sequencer stretch modification.

We have optimized the conditions for using the Stretch modification for the Applied Biosystems 373 Automated DNA Sequencers for sequencing double-stranded DNA using 34-cm well-to-read and 48-cm well-to-read configurations. With the manufacturer's recommended settings, uneven spacing within the first 100 bases was observed, which led to miscalls, insertions and deletions in the analyzed data. A significant decrease in accuracy for reads greater than 400 bases was also observed. Various gel concentrations were tested to improve the base spacing for the first 100 bases while maintaining accuracy and usable length of data. A longer average usable length and better resolution of smaller fragments were achieved by increased acrylamide concentration coupled with increased wattage. Using the Applied Biosystems CATALYST 800 Molecular Biology LabStation, Taq dye primer cycle sequencing reactions were optimized for -21 M13 and M13RP1 primers to produce a more even distribution of dye-labeled fragments that increased the overall signal strengths and decreased background signal. These reaction products, run on the Stretch sequencers using the new gel conditions, provided longer reads with increased reliability and accuracy of the data.

Multipurpose vectors for peptide expression on the M13 viral surface.

We have developed a set of three cloning vectors for the expression of polypeptides on the surface of the M13 viral coat. The M13mp8 genome has been engineered for expression of foreign protein sequences near the NH2-terminus of the mature pIII protein, which is present in five copies on the outside of each M13 viral particle. All three of the vectors carry the same two useful restriction sites for directed cloning of inserts in the pIII coding region; in addition, one vector carries the bacterial gene conferring resistance to the antibiotic tetracycline, and another expresses the lacZ' polypeptide that allows functional complementation of beta-galactosidase activity within the host bacterial cell. All of these vectors propagate well in E. coli DH5 alpha F' cells and do not require helper phage. We demonstrate that a bacteriophage, expressing an eleven amino acid epitope (from human c-myc) at the NH2-terminus of pIII in one of our vectors, can be purified from a vast mixture of other M13 phage through panning techniques. In particular, we find that the c-myc-expressing viral particles can be easily recovered from phage mixtures with the biotinylated form of the monoclonal antibody, 9E10, and streptavidin-coated MagneSphere beads.