A survey of expressed genes in Caenorhabditis elegans.

As an adjunct to the genomic sequencing of Caenorhabditis elegans, we have investigated a representative cDNA library of 1,517 clones. A single sequence read has been obtained from the 5' end of each clone, allowing its characterization with respect to the public databases, and the clones are being localized on the genome map. The result is the identification of about 1,200 of the estimated 15,000 genes of C. elegans. More than 30% of the inferred protein sequences have significant similarity to existing sequences in the databases, providing a route towards in vivo analysis of known genes in the nematode. These clones also provide material for assessing the accuracy of predicted exons and splicing patterns and will lead to a more accurate estimate of the total number of genes in the organism than has hitherto been available.     

Mutational hot spot within a new RPGR exon in X-linked retinitis pigmentosa

The gene RPGR was previously identified in the RP3 region of Xp21.1 and shown to be mutated in 10-20% of patients with the progressive retinal degeneration X-linked retinitis pigmentosa (XLRP). The mutations predominantly affected a domain homologous to RCC1, a guanine nucleotide exchange factor for the small GTPase Ran, although they were present in fewer than the 70-75% of XLRP patients predicted from linkage studies. Mutations in the RP2 locus at Xp11.3 were found in a further 10-20% of XLRP patients, as predicted from linkage studies. Because the mutations in the remainder of the XLRP patients may reside in undiscovered exons of RPGR, we sequenced a 172-kb region containing the entire gene. Analysis of the sequence disclosed a new 3' terminal exon that was mutated in 60% of XLRP patients examined. This exon encodes 567 amino acids, with a repetitive domain rich in glutamic acid residues. The sequence is conserved in the mouse, bovine and Fugu rubripes genes. It is preferentially expressed in mouse and bovine retina, further supporting its importance for retinal function. Our results suggest that mutations in RPGR are the only cause of RP3 type XLRP and account for the disease in over 70% of XLRP patients and an estimated 11% of all retinitis pigmentosa patients.     

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.     

Nova regulates brain-specific splicing to shape the synapse

Alternative RNA splicing greatly increases proteome diversity and may thereby contribute to tissue-specific functions. We carried out genome-wide quantitative analysis of alternative splicing using a custom Affymetrix microarray to assess the role of the neuronal splicing factor Nova in the brain. We used a stringent algorithm to identify 591 exons that were differentially spliced in the brain relative to immune tissues, and 6.6% of these showed major splicing defects in the neocortex of Nova2-/- mice. We tested 49 exons with the largest predicted Nova-dependent splicing changes and validated all 49 by RT-PCR. We analyzed the encoded proteins and found that all those with defined brain functions acted in the synapse (34 of 40, including neurotransmitter receptors, cation channels, adhesion and scaffold proteins) or in axon guidance (8 of 40). Moreover, of the 35 proteins with known interaction partners, 74% (26) interact with each other. Validating a large set of Nova RNA targets has led us to identify a multi-tiered network in which Nova regulates the exon content of RNAs encoding proteins that interact in the synapse.     

Open-reading-frame sequence tags (OSTs) support the existence of at least 17,300 genes in C. elegans

The genome sequences of Caenorhabditis elegans, Drosophila melanogaster and Arabidopsis thaliana have been predicted to contain 19,000, 13,600 and 25,500 genes, respectively. Before this information can be fully used for evolutionary and functional studies, several issues need to be addressed. First, the gene number estimates obtained in silico and not yet supported by any experimental data need to be verified. For example, it seems biologically paradoxical that C. elegans would have 50% more genes than Drosophilia. Second, intron/exon predictions need to be tested experimentally. Third, complete sets of open reading frames (ORFs), or "ORFeomes," need to be cloned into various expression vectors. To address these issues simultaneously, we have designed and applied to C. elegans the following strategy. Predicted ORFs are amplified by PCR from a highly representative cDNA library using ORF-specific primers, cloned by Gateway recombination cloning and then sequenced to generate ORF sequence tags (OSTs) as a way to verify identity and splicing. In a sample (n=1,222) of the nearly 10,000 genes predicted ab initio (that is, for which no expressed sequence tag (EST) is available so far), at least 70% were verified by OSTs. We also observed that 27% of these experimentally confirmed genes have a structure different from that predicted by GeneFinder. We now have experimental evidence that supports the existence of at least 17,300 genes in C. elegans. Hence we suggest that gene counts based primarily on ESTs may underestimate the number of genes in human and in other organisms.     

Restriction enzyme-generated siRNA (REGS) vectors and libraries

Small interfering RNA (siRNA) technology facilitates the study of loss of gene function in mammalian cells and animal models, but generating multiple siRNA vectors using oligonucleotides is slow, inefficient and costly. Here we describe a new, enzyme-mediated method for generating numerous functional siRNA constructs from any gene of interest or pool of genes. To test our restriction enzyme-generated siRNA (REGS) system, we silenced a transgene and two endogenous genes and obtained the predicted phenotypes. REGS generated on average 34 unique siRNAs per kilobase of sequence. REGS enabled us to create enzymatically a complex siRNA library (>4 x 10(5) clones) from double-stranded cDNA encompassing known and unknown genes with 96% of the clones containing inserts of the appropriate size.     

An azoospermic man with a de novo point mutation in the Y-chromosomal gene USP9Y

In humans, deletion of any one of three Y-chromosomal regions- AZFa, AZFb or AZFc-disrupts spermatogenesis, causing infertility in otherwise healthy men. Although candidate genes have been identified in all three regions, no case of spermatogenic failure has been traced to a point mutation in a Y-linked gene, or to a deletion of a single Y-linked gene. We sequenced the AZFa region of the Y chromosome and identified two functional genes previously described: USP9Y (also known as DFFRY) and DBY (refs 7,8). Screening of the two genes in 576 infertile and 96 fertile men revealed several sequence variants, most of which appear to be heritable and of little functional consequence. We found one de novo mutation in USP9Y: a 4-bp deletion in a splice-donor site, causing an exon to be skipped and protein truncation. This mutation was present in a man with nonobstructive azoospermia (that is, no sperm was detected in semen), but absent in his fertile brother, suggesting that the USP9Y mutation caused spermatogenic failure. We also identified a single-gene deletion associated with spermatogenic failure, again involving USP9Y, by re-analysing a published study.     

Computational identification of promoters and first exons in the human genome.

The identification of promoters and first exons has been one of the most difficult problems in gene-finding. We present a set of discriminant functions that can recognize structural and compositional features such as CpG islands, promoter regions and first splice-donor sites. We explain the implementation of the discriminant functions into a decision tree that constitutes a new program called FirstEF. By using different models to predict CpG-related and non-CpG-related first exons, we showed by cross-validation that the program could predict 86% of the first exons with 17% false positives. We also demonstrated the prediction accuracy of FirstEF at the genome level by applying it to the finished sequences of human chromosomes 21 and 22 as well as by comparing the predictions with the locations of the experimentally verified first exons. Finally, we present the analysis of the predicted first exons for all of the 24 chromosomes of the human genome.     

Genome-wide analysis of single-nucleotide polymorphisms in human expressed sequences

Single-nucleotide polymorphisms (SNPs) have been explored as a high-resolution marker set for accelerating the mapping of disease genes. Here we report 48,196 candidate SNPs detected by statistical analysis of human expressed sequence tags (ESTs), associated primarily with coding regions of genes. We used Bayesian inference to weigh evidence for true polymorphism versus sequencing error, misalignment or ambiguity, misclustering or chimaeric EST sequences, assessing data such as raw chromatogram height, sharpness, overlap and spacing, sequencing error rates, context-sensitivity and cDNA library origin. Three separate validations-comparison with 54 genes screened for SNPs independently, verification of HLA-A polymorphisms and restriction fragment length polymorphism (RFLP) testing-verified 70%, 89% and 71% of our predicted SNPs, respectively. Our method detects tenfold more true HLA-A SNPs than previous analyses of the EST data. We found SNPs in a large fraction of known disease genes, including some disease-causing mutations (for example, the HbS sickle-cell mutation). Our comprehensive analysis of human coding region polymorphism provides a public resource for mapping of disease genes (available at http://www.bioinformatics.ucla.edu/snp).     

Cloning of the alpha-adducin gene from the Huntington's disease candidate region of chromosome 4 by exon amplification.

We have applied the technique of exon amplification to the isolation of genes from the chromosome 4p16.3 Huntington's disease (HD) candidate region. Exons recovered from cosmid Y24 identified cDNA clones corresponding to the alpha-subunit of adducin, a calmodulin-binding protein that is thought to promote assembly of spectrin-actin complexes in the formation of the membrane cytoskeleton, alpha-adducin is widely expressed and, at least in brain, is encoded by alternatively spliced mRNAs. The alpha-adducin gene maps immediately telomeric to D4S95, in a region likely to contain the HD defect, and must be scrutinized to establish whether it is the site of the HD mutation.     

A negative element in SMN2 exon 7 inhibits splicing in spinal muscular atrophy

Spinal muscular atrophy (SMA) is a relatively common neurodegenerative disease caused by homozygous loss of the survival motor neuron 1 (SMN1) gene. Humans possess a linked, nearly identical gene, SMN2, which produces a functional SMN protein but at levels insufficient to compensate for loss of SMN1 (refs. 1,2). A C/T transition at position +6 in exon 7 is all that differentiates the two genes, but this is sufficient to prevent efficient exon 7 splicing in SMN2 (refs. 2,3). Here we show that the C/T transition functions not to disrupt an exonic splicing enhancer (ESE) in SMN1 (ref. 4), as previously suggested, but rather to create an exonic splicing silencer (ESS) in SMN2. We show that this ESS functions as a binding site for a known repressor protein, hnRNP A1, which binds to SMN2 but not SMN1 exon 7 RNA. We establish the physiological importance of these results by using small interfering RNAs to reduce hnRNP A protein levels in living cells and show that this results in efficient SMN2 exon 7 splicing. Our findings not only define a new mechanism underlying the inefficient splicing of SMN2 exon 7 but also illustrate more generally the remarkable sensitivity and precision that characterizes control of mRNA splicing.     

Large scale cDNA sequencing for analysis of quantitative and qualitative aspects of gene expression.

Large scale sequencing of cDNAs provides a complementary approach to structural analysis of the human genome by generating expressed sequence tags (ESTs). We have initiated the large-scale sequencing of a 3'-directed cDNA library from the human liver cell line HepG2, that is a non-biased representation of the mRNA population. 982 random cDNA clones were sequenced yielding more than 270 kilobases. A significant portion of the identified genes encoded secretable proteins and components for protein-synthesis. The abundance of cDNA species varied from 2.2% to less than 0.004%. Fifty two percent of the mRNA were abundant species consisting of 173 genes and the rest were non-abundant, consisting of about 6,600 genes.