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).     

Molecular breeding of viruses

Genetic recombination is a major force driving the evolution of many viruses. Recombination between two copackaged retroviral genomes may occur at rates as high as 40% per replication cycle. This enables genetic information to be shuffled rapidly, leading to recombinants with new patterns of mutations and phenotypes. The in vitro process of DNA shuffling (molecular breeding) mimics this mechanism on a vastly parallel and accelerated scale. Multiple homologous parental sequences are recombined in parallel, leading to a diverse library of complex recombinants from which desired improvements can be selected. Different proteins and enzymes have been improved using DNA shuffling. We report here the first application of molecular breeding to viruses. A single round of shuffling envelope sequences from six murine leukaemia viruses (MLV) followed by selection yielded a chimaeric clone with a completely new tropism for Chinese Hamster Ovary (CHOK1) cells. The composition and properties of the selected clone indicated that this particular permutation of parental sequences cannot be readily attained by natural retroviral recombination. This example demonstrates that molecular breeding can enhance the inherently high evolutionary potential of retroviruses to obtain desired phenotypes. It can be an effective tool, when information is limited, to optimize viruses for gene therapy and vaccine applications when multiple complex functions must be simultaneously balanced.     

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.     

Control of neurulation by the nucleosome assembly protein-1-like 2

Neurulation is a complex process of histogenesis involving the precise temporal and spatial organization of gene expression. Genes influencing neurulation include proneural genes determining primary cell fate, neurogenic genes involved in lateral inhibition pathways and genes controlling the frequency of mitotic events. This is reflected in the aetiology and genetics of human and mouse neural tube defects, which are of both multifactorial and multigenic origin. The X-linked gene Nap1l2, specifically expressed in neurons, encodes a protein that is highly similar to the nucleosome assembly (NAP) and SET proteins. We inactivated Nap1l2 in mice by gene targeting, leading to embryonic lethality from mid-gestation onwards. Surviving mutant chimaeric embryos showed extensive surface ectoderm defects as well as the presence of open neural tubes and exposed brains similar to those observed in human spina bifida and anencephaly. These defects correlated with an overproduction of neuronal precursor cells. Protein expression studies showed that the Nap1l2 protein binds to condensing chromatin during S phase and in apoptotic cells, but remained cytoplasmic during G1 phase. Nap1l2 therefore likely represents a class of tissue-specific factors interacting with chromatin to regulate neuronal cell proliferation.     

A new type of mutation causes a splicing defect in ATM

Disease-causing splicing mutations described in the literature primarily produce changes in splice sites and, to a lesser extent, variations in exon-regulatory sequences such as the enhancer elements. The gene ATM is mutated in individuals with ataxia-telangiectasia; we have identified the aberrant inclusion of a cryptic exon of 65 bp in one affected individual with a deletion of four nucleotides (GTAA) in intron 20. The deletion is located 12 bp downstream and 53 bp upstream from the 5' and 3' ends of the cryptic exon, respectively. Through analysis of the splicing defect using a hybrid minigene system, we identified a new intron-splicing processing element (ISPE) complementary to U1 snRNA, the RNA component of the U1 small nuclear ribonucleoprotein (snRNP). This element mediates accurate intron processing and interacts specifically with U1 snRNP particles. The 4-nt deletion completely abolished this interaction, causing activation of the cryptic exon. On the basis of this analysis, we describe a new type of U1 snRNP binding site in an intron that is essential for accurate intron removal. Deletion of this sequence is directly involved in the splicing processing defect.     

New genes involved in cancer identified by retroviral tagging

Retroviral insertional mutagenesis in BXH2 and AKXD mice induces a high incidence of myeloid leukemia and B- and T-cell lymphoma, respectively. The retroviral integration sites (RISs) in these tumors thus provide powerful genetic tags for the discovery of genes involved in cancer. Here we report the first large-scale use of retroviral tagging for cancer gene discovery in the post-genome era. Using high throughput inverse PCR, we cloned and analyzed the sequences of 884 RISs from a tumor panel composed primarily of B-cell lymphomas. We then compared these sequences, and another 415 RIS sequences previously cloned from BXH2 myeloid leukemias and from a few AKXD lymphomas, against the recently assembled mouse genome sequence. These studies identified 152 loci that are targets of retroviral integration in more than one tumor (common retroviral integration sites, CISs) and therefore likely to encode a cancer gene. Thirty-six CISs encode genes that are known or predicted to be genes involved in human cancer or their homologs, whereas others encode candidate genes that have not yet been examined for a role in human cancer. Our studies demonstrate the power of retroviral tagging for cancer gene discovery in the post-genome era and indicate a largely unrecognized complexity in mouse and presumably human cancer.     

Construction of a mouse yeast artificial chromosome library in a recombination-deficient strain of yeast.

We have constructed a new generation yeast artificial chromosome (YAC) library from female C57BL/10 mice in a recombination-deficient strain of Saccharomyces cerevisiae carrying a mutation in the RAD52 gene. The YAC library contains 41,568 clones with an average insert size of 240 kilobases, representing a greater than threefold coverage of the mouse genome. Currently, the library can be screened by polymerase chain reaction and we have isolated positive clones at a number of loci in the mouse genome. This rad52 library should enable a long-term assessment of the effect of one of the yeast recombination pathway genes on both, genome-wide YAC clone stability and the frequency of chimaeric clones.     

Insertional mutagenesis in zebrafish rapidly identifies genes essential for early vertebrate development

To rapidly identify genes required for early vertebrate development, we are carrying out a large-scale, insertional mutagenesis screen in zebrafish, using mouse retroviral vectors as the mutagen. We will obtain mutations in 450 to 500 different genes--roughly 20% of the genes that can be mutated to produce a visible embryonic phenotype in this species--and will clone the majority of the mutated alleles. So far, we have isolated more than 500 insertional mutants. Here we describe the first 75 insertional mutants for which the disrupted genes have been identified. In agreement with chemical mutagenesis screens, approximately one-third of the mutants have developmental defects that affect primarily one or a small number of organs, body shape or swimming behavior; the rest of the mutants show more widespread or pleiotropic abnormalities. Many of the genes we identified have not been previously assigned a biological role in vivo. Roughly 20% of the mutants result from lesions in genes for which the biochemical and cellular function of the proteins they encode cannot be deduced with confidence, if at all, from their predicted amino-acid sequences. All of the genes have either orthologs or clearly related genes in human. These results provide an unbiased view of the genetic construction kit for a vertebrate embryo, reveal the diversity of genes required for vertebrate development and suggest that hundreds of genes of unknown biochemical function essential for vertebrate development have yet to be identified.     

Regional loss of imprinting and growth deficiency in mice with a targeted deletion of KvDMR1

Genomic imprinting is an epigenetic modification that results in expression from only one of the two parental copies of a gene. Differences in methylation between the two parental chromosomes are often observed at or near imprinted genes. Beckwith-Wiedemann syndrome (BWS), which predisposes to cancer and excessive growth, results from a disruption of imprinted gene expression in chromosome band 11p15.5. One third of individuals with BWS lose maternal-specific methylation at KvDMR1, a putative imprinting control region within intron 10 of the KCNQ1 gene, and it has been proposed that this epimutation results in aberrant imprinting and, consequently, BWS1, 2. Here we show that paternal inheritance of a deletion of KvDMR1 results in the de-repression in cis of six genes, including Cdkn1c, which encodes cyclin-dependent kinase inhibitor 1C. Furthermore, fetuses and adult mice that inherited the deletion from their fathers were 20-25% smaller than their wildtype littermates. By contrast, maternal inheritance of this deletion had no effect on imprinted gene expression or growth. Thus, the unmethylated paternal KvDMR1 allele regulates imprinted expression by silencing genes on the paternal chromosome. These findings support the hypothesis that loss of methylation in BWS patients activates the repressive function of KvDMR1 on the maternal chromosome, resulting in abnormal silencing of CDKN1C and the development of BWS.     

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.     

Mouse mutants from chemically mutagenized embryonic stem cells

The drive to characterize functions of human genes on a global scale has stimulated interest in large-scale generation of mouse mutants. Conventional germ-cell mutagenesis with N-ethyl-N-nitrosourea (ENU) is compromised by an inability to monitor mutation efficiency, strain and interlocus variation in mutation induction, and extensive husbandry requirements. To overcome these obstacles and develop new methods for generating mouse mutants, we devised protocols to generate germline chimaeric mice from embryonic stem (ES) cells heavily mutagenized with ethylmethanesulphonate (EMS). Germline chimaeras were derived from cultures that underwent a mutation rate of up to 1 in 1,200 at the Hprt locus (encoding hypoxanthine guanine phosphoribosyl transferase). The spectrum of mutations induced by EMS and the frameshift mutagen ICR191 was consistent with that observed in other mammalian cells. Chimaeras derived from ES cells treated with EMS transmitted mutations affecting several processes, including limb development, hair growth, hearing and gametogenesis. This technology affords several advantages over traditional mutagenesis, including the ability to conduct shortened breeding schemes and to screen for mutant phenotypes directly in ES cells or their differentiated derivatives.     

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.     

Closing gaps in the human genome with fosmid resources generated from multiple individuals

The human genome sequence has been finished to very high standards; however, more than 340 gaps remained when the finished genome was published by the International Human Genome Sequencing Consortium in 2004. Using fosmid resources generated from multiple individuals, we targeted gaps in the euchromatic part of the human genome. Here we report 2,488,842 bp of previously unknown euchromatic sequence, 363,114 bp of which close 26 of 250 euchromatic gaps, or 10%, including two remaining euchromatic gaps on chromosome 19. Eight (30.7%) of the closed gaps were found to be polymorphic. These sequences allow complete annotation of several human genes as well as the assignment of mRNAs. The gap sequences are 2.3-fold enriched in segmentally duplicated sequences compared to the whole genome. Our analysis confirms that not all gaps within 'finished' genomes are recalcitrant to subcloning and suggests that the paired-end-sequenced fosmid libraries could prove to be a rich resource for completion of the human euchromatic genome.     

Chimaeric mice deficient in dystroglycans develop muscular dystrophy and have disrupted myoneural synapses

Mutations in the dystrophin gene (DMD) and in genes encoding several dystrophin-associated proteins result in Duchenne and other forms of muscular dystrophy. alpha-Dystroglycan (Dg) binds to laminins in the basement membrane surrounding each myofibre and docks with beta-Dg, a transmembrane protein, which in turn interacts with dystrophin or utrophin in the subplasmalemmal cytoskeleton. alpha- and beta-Dgs are thought to form the functional core of a larger complex of proteins extending from the basement membrane to the intracellular cytoskeleton, which serves as a superstructure necessary for sarcolemmal integrity. Dgs have also been implicated in the formation of synaptic densities of acetylcholine receptors (AChRs) on skeletal muscle. Here we report that chimaeric mice generated with ES cells targeted for both Dg alleles have skeletal muscles essentially devoid of Dgs and develop a progressive muscle pathology with changes emblematic of muscular dystrophies in humans. In addition, many neuromuscular junctions are disrupted in these mice. The ultrastructure of basement membranes and the deposition of laminin within them, however, appears unaffected in Dg-deficient muscles. We conclude that Dgs are necessary for myofibre survival and synapse differentiation or stability, but not for the formation of the muscle basement membrane, and that Dgs may have more than a purely structural function in maintaining muscle integrity.     

Duplication-degeneration as a mechanism of gene fission and the origin of new genes in Drosophila species

Gene fission and fusion, the processes by which a single gene is split into two separate genes and two adjacent genes are fused into a single gene, respectively, are among the primary processes that generate new genes. Despite their seeming reversibility, nothing is known about the mechanism of gene fission. Because the nucleotide sequences of fission genes record little about their origination process, conventional analysis of duplicate genes may not be powerful enough to unravel the underlying mechanism. In a survey for young genes in species of the Drosophila melanogaster subgroup using fluorescence in situ hybridization, we identified a young gene family, monkey king, whose genesis sheds light on the evolutionary process of gene fission. Its members originated 1-2 million years ago as retroposed duplicates and evolved into fission genes that separately encode protein domains from a multidomain ancestor. The mechanism underlying this process is gene duplication with subsequent partial degeneration.     

Homology-based annotation yields 1,042 new candidate genes in the Drosophila melanogaster genome

The approach to annotating a genome critically affects the number and accuracy of genes identified in the genome sequence. Genome annotation based on stringent gene identification is prone to underestimate the complement of genes encoded in a genome. In contrast, over-prediction of putative genes followed by exhaustive computational sequence, motif and structural homology search will find rarely expressed, possibly unique, new genes at the risk of including non-functional genes. We developed a two-stage approach that combines the merits of stringent genome annotation with the benefits of over-prediction. First we identify plausible genes regardless of matches with EST, cDNA or protein sequences from the organism (stage 1). In the second stage, proteins predicted from the plausible genes are compared at the protein level with EST, cDNA and protein sequences, and protein structures from other organisms (stage 2). Remote but biologically meaningful protein sequence or structure homologies provide supporting evidence for genuine genes. The method, applied to the Drosophila melanogaster genome, validated 1,042 novel candidate genes after filtering 19,410 plausible genes, of which 12,124 matched the original 13,601 annotated genes. This annotation strategy is applicable to genomes of all organisms, including human.