Rapid gene mapping in Caenorhabditis elegans using a high density polymorphism map

Single nucleotide polymorphisms (SNPs) are valuable genetic markers of human disease. They also comprise the highest potential density marker set available for mapping experimentally derived mutations in model organisms such as Caenorhabditis elegans. To facilitate the positional cloning of mutations we have identified polymorphisms in CB4856, an isolate from a Hawaiian island that shows a uniformly high density of polymorphisms compared with the reference Bristol N2 strain. Based on 5.4 Mbp of aligned sequences, we predicted 6,222 polymorphisms. Furthermore, 3,457 of these markers modify restriction enzyme recognition sites ('snip-SNPs') and are therefore easily detected as RFLPs. Of these, 493 were experimentally confirmed by restriction digest to produce a snip-SNP map of the worm genome. A mapping strategy using snip-SNPs and bulked segregant analysis (BSA) is outlined. CB4856 is crossed into a mutant strain, and exclusion of CB4856 alleles of a subset of snip-SNPs in mutant progeny is assessed with BSA. The proximity of a linked marker to the mutation is estimated by the relative proportion of each form of the biallelic marker in populations of wildtype and mutant genomes. The usefulness of this approach is illustrated by the rapid mapping of the dyf-5 gene.     

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.     

A single gene network accurately predicts phenotypic effects of gene perturbation in Caenorhabditis elegans

The fundamental aim of genetics is to understand how an organism's phenotype is determined by its genotype, and implicit in this is predicting how changes in DNA sequence alter phenotypes. A single network covering all the genes of an organism might guide such predictions down to the level of individual cells and tissues. To validate this approach, we computationally generated a network covering most C. elegans genes and tested its predictive capacity. Connectivity within this network predicts essentiality, identifying this relationship as an evolutionarily conserved biological principle. Critically, the network makes tissue-specific predictions-we accurately identify genes for most systematically assayed loss-of-function phenotypes, which span diverse cellular and developmental processes. Using the network, we identify 16 genes whose inactivation suppresses defects in the retinoblastoma tumor suppressor pathway, and we successfully predict that the dystrophin complex modulates EGF signaling. We conclude that an analogous network for human genes might be similarly predictive and thus facilitate identification of disease genes and rational therapeutic targets.     

Evidence for co-evolution of gene order and recombination rate

There is increasing evidence in eukaryotic genomes that gene order is not random, even allowing for tandem duplication. Notably, in numerous genomes, genes of similar expression tend to be clustered. Are there other reasons for clustering of functionally similar genes? If genes are linked to enable genetic, rather than physical clustering, then we also expect that clusters of certain genes might be associated with blocks of reduced recombination rates. Here we show that, in yeast, essential genes are highly clustered and this clustering is independent of clustering of co-expressed genes and of tandem duplications. Adjacent pairs of essential genes are preferentially conserved through evolution. Notably, we also find that clusters of essential genes are in regions of low recombination and that larger clusters have lower recombination rates. These results suggest that selection acts to modify both the fine-scale intragenomic variation in the recombination rate and the distribution of genes and provide evidence for co-evolution of gene order and recombination rate.     

Effective targeted gene 'knockdown' in zebrafish

The sequencing of the zebrafish genome should be completed by the end of 2002. Direct assignment of function on the basis of this information would be facilitated by the development of a rapid, targeted 'knockdown' technology in this model vertebrate. We show here that antisense, morpholino-modified oligonucleotides (morpholinos) are effective and specific translational inhibitors in zebrafish. We generated phenocopies of mutations of the genes no tail (ref. 2), chordin (ref. 3), one-eyed-pinhead (ref. 4), nacre (ref. 5) and sparse (ref. 6), removing gene function from maternal through post-segmentation and organogenesis developmental stages. We blocked expression from a ubiquitous green fluorescent protein (GFP) transgene, showing that, unlike tissue-restricted limitations found with RNA-based interference in the nematode, all zebrafish cells readily respond to this technique. We also developed also morpholino-based zebrafish models of human disease. Morpholinos targeted to the uroporphyrinogen decarboxylase gene result in embryos with hepatoerythropoietic porphyria. We also used morpholinos for the determination of new gene functions. We showed that embryos with reduced sonic hedgehog (ref. 9) signalling and reduced tiggy-winkle hedgehog (ref. 10) function exhibit partial cyclopia and other specific midline abnormalities, providing a zebrafish genetic model for the common human disorder holoprosencephaly. Conserved vertebrate processes and diseases are now amenable to a systematic, in vivo, reverse-genetic paradigm using zebrafish embryos.     

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.     

Disruption of dog-1 in Caenorhabditis elegans triggers deletions upstream of guanine-rich DNA

Genetic integrity is crucial to normal cell function, and mutations in genes required for DNA replication and repair underlie various forms of genetic instability and disease, including cancer. One structural feature of intact genomes is runs of homopolymeric dC/dG. Here we describe an unusual mutator phenotype in Caenorhabditis elegans characterized by deletions that start around the 3' end of polyguanine tracts and terminate at variable positions 5' from such tracts. We observed deletions throughout genomic DNA in about half of polyguanine tracts examined, especially those containing 22 or more consecutive guanine nucleotides. The mutator phenotype results from disruption of the predicted gene F33H2.1, which encodes a protein with characteristics of a DEAH helicase and which we have named dog-1 (for deletions of guanine-rich DNA). Nematodes mutated in dog-1 showed germline as well as somatic deletions in genes containing polyguanine tracts, such as vab-1. We propose that DOG-1 is required to resolve the secondary structures of guanine-rich DNA that occasionally form during lagging-strand DNA synthesis.     

Amplification of siRNA in Caenorhabditis elegans generates a transgenerational sequence-targeted histone H3 lysine 9 methylation footprint

Exogenous double-stranded RNA (dsRNA) has been shown to exert homology-dependent effects at the level of both target mRNA stability and chromatin structure. Using C. elegans undergoing RNAi as an animal model, we have investigated the generality, scope and longevity of dsRNA-targeted chromatin effects and their dependence on components of the RNAi machinery. Using high-resolution genome-wide chromatin profiling, we found that a diverse set of genes can be induced to acquire locus-specific enrichment of histone H3 lysine 9 trimethylation (H3K9me3), with modification footprints extending several kilobases from the site of dsRNA homology and with locus specificity sufficient to distinguish the targeted locus from the other 20,000 genes in the C. elegans genome. Genetic analysis of the response indicated that factors responsible for secondary siRNA production during RNAi were required for effective targeting of chromatin. Temporal analysis revealed that H3K9me3, once triggered by dsRNA, can be maintained in the absence of dsRNA for at least two generations before being lost. These results implicate dsRNA-triggered chromatin modification in C. elegans as a programmable and locus-specific response defining a metastable state that can persist through generational boundaries.     

AAV serotype 2 vectors preferentially integrate into active genes in mice

Recombinant adeno-associated virus serotype 2 (rAAV2) is a promising vector for gene therapy because it can achieve long-term stable transgene expression in animals and human subjects after direct administration of vectors into various target tissues. In the liver, although stable transgene expression primarily results from extrachromosomal vector genomes, a series of experiments has shown that vector genomes integrate into host chromosomes in hepatocytes at a low frequency. Despite the low integration efficiency, recent reports of retroviral insertional mutagenesis in mice and two human subjects have raised concerns about the potential for rAAV2-mediated insertional mutagenesis. Here we characterize rAAV2-targeted chromosomal integration sites isolated from selected or non-selected hepatocytes in vector-injected mouse livers. We document frequent chromosomal deletions of up to 2 kb at integration sites (14 of 14 integrations, 100%; most of the deletions were <0.3 kb) and preferred integration into genes (21 of 29 integrations, 72%). In addition, all of the targeted genes analyzed (20 of 20 targeted genes, 100%) were expressed in the liver. This is the first report to our knowledge on host chromosomal effects of rAAV2 integration in animals, and it provides insights into the nature of rAAV2 vector integration into chromosomes in quiescent somatic cells in animals and human subjects.     

Conservation and diversification of Wnt signaling function during the evolution of nematode vulva development

Cell-fate specification and cell-cell signaling have been well studied during vulva development in Caenorhabditis elegans and provide a paradigm in evolutionary developmental biology. Pristionchus pacificus has been developed as a 'satellite' organism with an integrated physical and genetic map that allows detailed comparisons to C. elegans. A common aspect of vulva formation in both species is the polarization of the P7.p lineage, which is responsible for vulval symmetry. In C. elegans, Wnt signaling is crucial for P7.p cell-fate patterning; nothing is known about vulval symmetry in P. pacificus. We isolated mutations that disrupt polarization of the P7.p lineage in P. pacificus and found that the corresponding gene encodes a Frizzled-like molecule. In addition, mutations in Ppa-lin-17 (encoding Frizzled) and morpholino knock-down of Ppa-lin-44 (encoding Wnt), Ppa-egl-20 (encoding Wnt), Ppa-mig-5 (encoding Dsh), Ppa-apr-1 (encoding APC) and Ppa-bar-1 (encoding beta-catenin) results in gonad-independent vulva differentiation, indicating that these genes have a role in a negative signaling process. In contrast, in C. elegans, Wnt signaling has a positive role in vulva induction, and mutations in bar-1 result in a hypoinduced phenotype. Therefore, whereas the molecular mechanisms that generate vulval symmetry are conserved, the genetic control of vulva induction diversified during evolution.     

A natural allele of Nxf1 suppresses retrovirus insertional mutations

Endogenous retroviruses have shaped the evolution of mammalian genomes. Host genes that control the effects of retrovirus insertions are therefore of great interest. The modifier-of-vibrator-1 locus (Mvb1) controls levels of correctly processed mRNA from genes mutated by endogenous retrovirus insertions into introns, including the Pitpn(vb) tremor mutation and the Eya1(BOR) model of human branchiootorenal syndrome. Positional complementation cloning identifies Mvb1 as the nuclear export factor Nxf1, providing an unexpected link between the mRNA export receptor and pre-mRNA processing. Population structure of the suppressive allele in wild Mus musculus castaneus suggests selective advantage. A congenic Mvb1(CAST) allele is a useful tool for modifying gene expression from existing mutations and could be used to manipulate engineered mutations containing retroviral elements.     

The complete family of genes encoding G proteins of Caenorhabditis elegans

Caenorhabditis elegans is the first animal whose genomic sequence has been determined. One of the new possibilities in post-sequence genetics is the analysis of complete gene families at once. We studied the family of heterotrimeric G proteins. C. elegans has 20 Galpha, 2 Gbeta and 2 Ggamma genes. There is 1 homologue of each of the 4 mammalian classes of Galpha genes, G(i)/G(o)alpha, G(s)alpha , G(q)alpha and G12alpha, and there are 16 new alpha genes. Although the conserved Galpha subunits are expressed in many neurons and muscle cells, GFP fusions indicate that 14 new Galpha genes are expressed almost exclusively in a small subset of the chemosensory neurons of C. elegans. We generated loss-of-function alleles using target-selected gene inactivation. None of the amphid-expressed genes are essential for viability, and only four show any detectable phenotype (chemotaxis defects), suggesting extensive functional redundancy. On the basis of functional analysis, the 20 genes encoding Galpha proteins can be divided into two groups: those that encode subunits affecting muscle activity (homologues of G(i)/G(o)alpha, G(s)alpha and G(q)), and those (14 new genes) that encode proteins most likely involved in perception.     

Identification of acquired somatic mutations in the gene encoding chromatin-remodeling factor ATRX in the alpha-thalassemia myelodysplasia syndrome (ATMDS)

Inherited mutations of specific genes have elucidated the normal roles of the proteins they encode by relating specific mutations to particular phenotypes. But many potentially informative mutations in such genes are lethal early in development. Consequently, inherited mutations may not reflect all the functional roles of such proteins. Acquired, somatic defects should reflect a wider spectrum of mutations because they are not prone to negative selection in development. It has been difficult to identify such mutations so far, but microarray analysis provides a new opportunity to do so. Using this approach, we have shown that in individuals with myelodysplasia associated with alpha-thalassemia (ATMDS), somatic mutations of the gene encoding the chromatin remodeling factor ATRX cause an unexpectedly severe hematological phenotype compared with the wide spectrum of inherited mutations affecting this gene. These findings cast new light on this pleiotropic cofactor, which appears to be an essential component rather than a mere facilitator of globin gene expression.     

C. elegans ORFeome version 1.1: experimental verification of the genome annotation and resource for proteome-scale protein expression

To verify the genome annotation and to create a resource to functionally characterize the proteome, we attempted to Gateway-clone all predicted protein-encoding open reading frames (ORFs), or the 'ORFeome,' of Caenorhabditis elegans. We successfully cloned approximately 12,000 ORFs (ORFeome 1.1), of which roughly 4,000 correspond to genes that are untouched by any cDNA or expressed-sequence tag (EST). More than 50% of predicted genes needed corrections in their intron-exon structures. Notably, approximately 11,000 C. elegans proteins can now be expressed under many conditions and characterized using various high-throughput strategies, including large-scale interactome mapping. We suggest that similar ORFeome projects will be valuable for other organisms, including humans.     

Modular epistasis in yeast metabolism

Epistatic interactions, manifested in the effects of mutations on the phenotypes caused by other mutations, may help uncover the functional organization of complex biological networks. Here, we studied system-level epistatic interactions by computing growth phenotypes of all single and double knockouts of 890 metabolic genes in Saccharomyces cerevisiae, using the framework of flux balance analysis. A new scale for epistasis identified a distinctive trimodal distribution of these epistatic effects, allowing gene pairs to be classified as buffering, aggravating or noninteracting. We found that the ensuing epistatic interaction network could be organized hierarchically into function-enriched modules that interact with each other 'monochromatically' (i.e., with purely aggravating or purely buffering epistatic links). This property extends the concept of epistasis from single genes to functional units and provides a new definition of biological modularity, which emphasizes interactions between, rather than within, functional modules. Our approach can be used to infer functional gene modules from purely phenotypic epistasis measurements.     

A systematic, genome-wide, phenotype-driven mutagenesis programme for gene function studies in the mouse

As the human genome project approaches completion, the challenge for mammalian geneticists is to develop approaches for the systematic determination of mammalian gene function. Mouse mutagenesis will be a key element of studies of gene function. Phenotype-driven approaches using the chemical mutagen ethylnitrosourea (ENU) represent a potentially efficient route for the generation of large numbers of mutant mice that can be screened for novel phenotypes. The advantage of this approach is that, in assessing gene function, no a priori assumptions are made about the genes involved in any pathway. Phenotype-driven mutagenesis is thus an effective method for the identification of novel genes and pathways. We have undertaken a genome-wide, phenotype-driven screen for dominant mutations in the mouse. We generated and screened over 26,000 mice, and recovered some 500 new mouse mutants. Our work, along with the programme reported in the accompanying paper, has led to a substantial increase in the mouse mutant resource and represents a first step towards systematic studies of gene function in mammalian genetics.