The status, quality, and expansion of the NIH full-length cDNA project: the Mammalian Gene Collection (MGC)

The National Institutes of Health's Mammalian Gene Collection (MGC) project was designed to generate and sequence a publicly accessible cDNA resource containing a complete open reading frame (ORF) for every human and mouse gene. The project initially used a random strategy to select clones from a large number of cDNA libraries from diverse tissues. Candidate clones were chosen based on 5'-EST sequences, and then fully sequenced to high accuracy and analyzed by algorithms developed for this project. Currently, more than 11,000 human and 10,000 mouse genes are represented in MGC by at least one clone with a full ORF. The random selection approach is now reaching a saturation point, and a transition to protocols targeted at the missing transcripts is now required to complete the mouse and human collections. Comparison of the sequence of the MGC clones to reference genome sequences reveals that most cDNA clones are of very high sequence quality, although it is likely that some cDNAs may carry missense variants as a consequence of experimental artifact, such as PCR, cloning, or reverse transcriptase errors. Recently, a rat cDNA component was added to the project, and ongoing frog (Xenopus) and zebrafish (Danio) cDNA projects were expanded to take advantage of the high-throughput MGC pipeline.

Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences

The National Institutes of Health Mammalian Gene Collection (MGC) Program is a multiinstitutional effort to identify and sequence a cDNA clone containing a complete ORF for each human and mouse gene. ESTs were generated from libraries enriched for full-length cDNAs and analyzed to identify candidate full-ORF clones, which then were sequenced to high accuracy. The MGC has currently sequenced and verified the full ORF for a nonredundant set of >9,000 human and >6,000 mouse genes. Candidate full-ORF clones for an additional 7,800 human and 3,500 mouse genes also have been identified. All MGC sequences and clones are available without restriction through public databases and clone distribution networks (see http:mgc.nci.nih.gov).

Pooled library tissue tags for EST-based gene discovery

MOTIVATION

In gene discovery projects based on EST sequencing, effective post-sequencing identification methods are important in determining tissue sources of ESTs within pooled cDNA libraries. In the past, such identification efforts have been characterized by higher than necessary failure rates due to the presence of errors within the subsequence containing the oligo tag intended to define the tissue source for each EST.

RESULTS

A large-scale EST-based gene discovery program at The University of Iowa has led to the creation of a unique software method named UITagCreator usable in the creation of large sets of synthetic tissue identification tags. The identification tags provide error detection and correction capability and, in conjunction with automated annotation software, result in a substantial improvement in the accurate identification of the tissue source in the presence of sequencing and base-calling errors. These identification rates are favorable, relative to past paradigms.

AVAILABILITY

The UITagCreator source code and installation instructions, along with detection software usable in concert with created tag sets, is freely available at http://genome.uiowa.edu/pubsoft/software.html

CONTACT

tomc@eng.uiowa.edu

Automated construction of high-density comparative maps between rat, human, and mouse

Animal models have been used primarily as surrogates for humans, having similar disease-based phenotypes. Genomic organization also tends to be conserved between species, leading to the generation of comparative genome maps. The emergence of radiation hybrid (RH) maps, coupled with the large numbers of available Expressed Sequence Tags (ESTs), has revolutionized the way comparative maps can be built. We used publicly available rat, mouse, and human data to identify genes and ESTs with interspecies sequence identity (homology), identified their UniGene relationships, and incorporated their RH map positions to build integrated comparative maps with >2100 homologous UniGenes mapped in more than one species (approximately 6% of all mammalian genes). The generation of these maps is iterative and labor intensive; therefore, we developed a series of computer tools (not described here) based on our algorithm that identifies anchors between species and produces printable and on-line clickable comparative maps that link to a wide variety of useful tools and databases. The maps were constructed using sequence-based comparisons, thus creating "hooks" for further sequence-based annotation of human, mouse, and rat sequences. Currently, this map enables investigators to link the physiology of the rat with the genetics of the mouse and the clinical significance of the human.

An analysis of allelic variation in the ABCA4 gene

PURPOSE

To assess the allelic variation of the ATP-binding transporter protein (ABCA4).

METHODS

A combination of single-strand conformation polymorphism (SSCP) and automated DNA sequencing was used to systematically screen this gene for sequence variations in 374 unrelated probands with a clinical diagnosis of Stargardt disease, 182 patients with age-related macular degeneration (AMD), and 96 normal subjects.

RESULTS

There was no significant difference in the proportion of any single variant or class of variant between the control and AMD groups. In contrast, truncating variants, amino acid substitutions, synonymous codon changes, and intronic variants were significantly enriched in patients with Stargardt disease when compared with their presence in subjects without Stargardt disease (Kruskal-Wallis P < 0.0001 for each variant group). Overall, there were 2480 instances of 213 different variants in the ABCA4 gene, including 589 instances of 97 amino acid substitutions, and 45 instances of 33 truncating variants.

CONCLUSIONS

Of the 97 amino acid substitutions, 11 occurred at a frequency that made them unlikely to be high-penetrance recessive disease-causing variants (HPRDCV). After accounting for variants in cis, one or more changes that were compatible with HPRDCV were found on 35% of all Stargardt-associated alleles overall. The nucleotide diversity of the ABCA4 coding region, a collective measure of the number and prevalence of polymorphic sites in a region of DNA, was found to be 1.28, a value that is 9 to 400 times greater than that of two other macular disease genes that were examined in a similar fashion (VMD2 and EFEMP1).

Generation of a high-density rat EST map

We have developed a high-density EST map of the rat, consisting of >11,000 ESTs. These ESTs were placed on a radiation hybrid framework map of genetic markers spanning all 20 rat autosomes, plus the X chromosome. The framework maps have a total size of approximately 12,400 cR, giving an average correspondence of 240 kb/cR. The frameworks are all LOD 3 chromosomal maps consisting of 775 radiation-hybrid-mapped genetic markers and ESTs. To date, we have generated radiation-hybrid-mapping data for >14,000 novel ESTs identified by our Rat Gene Discovery and Mapping Project (http://ratEST.uiowa.edu), from which we have placed >11,000 on our framework maps. To minimize mapping errors, ESTs were mapped in duplicate and consensus RH vectors produced for use in the placement procedure. This EST map was then used to construct high-density comparative maps between rat and human and rat and mouse. These maps will be a useful resource for positional cloning of genes for rat models of human diseases and in the creation and verification of a tiling set of map order for the upcoming rat-genome sequencing.

Generation of a high-density rat EST map

We have developed a high-density EST map of the rat, consisting of >11,000 ESTs. These ESTs were placed on a radiation hybrid framework map of genetic markers spanning all 20 rat autosomes, plus the X chromosome. The framework maps have a total size of approximately 12,400 cR, giving an average correspondence of 240 kb/cR. The frameworks are all LOD 3 chromosomal maps consisting of 775 radiation-hybrid-mapped genetic markers and ESTs. To date, we have generated radiation-hybrid-mapping data for >14,000 novel ESTs identified by our Rat Gene Discovery and Mapping Project (http://ratEST.uiowa.edu), from which we have placed >11,000 on our framework maps. To minimize mapping errors, ESTs were mapped in duplicate and consensus RH vectors produced for use in the placement procedure. This EST map was then used to construct high-density comparative maps between rat and human and rat and mouse. These maps will be a useful resource for positional cloning of genes for rat models of human diseases and in the creation and verification of a tiling set of map order for the upcoming rat-genome sequencing.

Anopheles gambiae pilot gene discovery project: identification of mosquito innate immunity genes from expressed sequence tags generated from immune-competent cell lines

Together with AIDS and tuberculosis, malaria is at the top of the list of devastating infectious diseases. However, molecular genetic studies of its major vector, Anopheles gambiae, are still quite limited. We have conducted a pilot gene discovery project to accelerate progress in the molecular analysis of vector biology, with emphasis on the mosquito's antimalarial immune defense. A total of 5,925 expressed sequence tags were determined from normalized cDNA libraries derived from immune-responsive hemocyte-like cell lines. The 3,242 expressed sequence tag-containing cDNA clones were grouped into 2,380 clone clusters, potentially representing unique genes. Of these, 1,118 showed similarities to known genes from other organisms, but only 27 were identical to previously known mosquito genes. We identified 38 candidate genes, based on sequence similarity, that may be implicated in immune reactions including antimalarial defense; 19 of these were shown experimentally to be inducible by bacterial challenge, lending support to their proposed involvement in mosquito immunity.

An autosomal genomic screen for autism. Collaborative linkage study of autism

Autism is a severe neurodevelopmental disorder defined by social and communication deficits and ritualistic-repetitive behaviors that are detectable in early childhood. The etiology of idiopathic autism is strongly genetic, and oligogenic transmission is likely. The first stage of a two-stage genomic screen for autism was carried out by the Collaborative Linkage Study of Autism on individuals affected with autism from 75 families ascertained through an affected sib-pair. The strongest multipoint results were for regions on chromosomes 13 and 7. The highest maximum multipoint heterogeneity LOD (MMLS/het) score is 3.0 at D13S800 (approximately 55 cM from the telomere) under the recessive model, with an estimated 35% of families linked to this locus. The next highest peak is an MMLS/het score of 2.3 at 19 cM, between D13S217 and D13S1229. Our third highest MMLS/het score of 2.2 is on chromosome 7 and is consistent with the International Molecular Genetic Study of Autism Consortium report of a possible susceptibility locus somewhere within 7q31-33. These regions and others will be followed up in the second stage of our study by typing additional markers in both the original and a second set of identically ascertained autism families, which are currently being collected. By comparing results across a number of studies, we expect to be able to narrow our search for autism susceptibility genes to a small number of genomic regions. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 88:609-615, 1999.