Enzymatic production of RNAi libraries from cDNAs

RNA interference (RNAi) induced by small interfering (siRNA) or short hairpin RNA (shRNA) is an important research approach in mammalian genetics. Here we describe a technology called enzymatic production of RNAi library (EPRIL) by which cDNAs are converted by a sequence of enzymatic treatments into an RNAi library consisting of a vast array of different shRNA expression constructs. We applied EPRIL to a single cDNA source and prepared an RNAi library consisting of shRNA constructs with various RNAi efficiencies. High-throughput screening allowed us to rapidly identify the best shRNA constructs from the library. We also describe a new selection scheme using the thymidine kinase gene for obtaining efficient shRNA constructs. Furthermore, we show that EPRIL can be applied to constructing an RNAi library from a cDNA library, providing a basis for future whole-genome phenotypic screening of genes.     

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.     

Genomic screening and replication using the same data set in family-based association testing

The Human Genome Project and its spin-offs are making it increasingly feasible to determine the genetic basis of complex traits using genome-wide association studies. The statistical challenge of analyzing such studies stems from the severe multiple-comparison problem resulting from the analysis of thousands of SNPs. Our methodology for genome-wide family-based association studies, using single SNPs or haplotypes, can identify associations that achieve genome-wide significance. In relation to developing guidelines for our screening tools, we determined lower bounds for the estimated power to detect the gene underlying the disease-susceptibility locus, which hold regardless of the linkage disequilibrium structure present in the data. We also assessed the power of our approach in the presence of multiple disease-susceptibility loci. Our screening tools accommodate genomic control and use the concept of haplotype-tagging SNPs. Our methods use the entire sample and do not require separate screening and validation samples to establish genome-wide significance, as population-based designs do.     

Mutations in a new gene encoding a thiamine transporter cause thiamine-responsive megaloblastic anaemia syndrome.

Thiamine-responsive megaloblastic anaemia syndrome (TRMA; MIM 249270) is an autosomal recessive disorder with features that include megaloblastic anaemia, mild thrombocytopenia and leucopenia, sensorineural deafness and diabetes mellitus. Treatment with pharmacologic doses of thiamine ameliorates the megaloblastic anaemia and diabetes mellitus. A defect in the plasma membrane transport of thiamine has been demonstrated in erythrocytes and cultured skin fibroblasts from TRMA patients. The gene causing TRMA was assigned to 1q23.2-q23.3 by linkage analysis. Here we report the cloning of a new gene, SLC19A2, identified from high-through-put genomic sequences due to homology with SLC19A1, encoding reduced folate carrier 1 (refs 8-10). We cloned the entire coding region by screening a human fetal brain cDNA library. SLC19A2 encodes a protein (of 497 aa) predicted to have 12 transmembrane domains. We identified 2 frameshift mutations in exon 2. a 1-bp insertion and a 2-bp deletion, among four Iranian families with TRMA. The sequence homology and predicted structure of SLC19A2, as well as its role in TRMA, suggest that its gene product is a thiamine carrier, the first to be identified in complex eukaryotes.     

Mitotic recombination is suppressed by chromosomal divergence in hybrids of distantly related mouse strains

Mitotic recombination occurs with high frequency in humans and mice. It leads to loss of heterozygosity (LOH) at important gene loci and can cause disease. However, the genetic modulators of mitotic recombination are not well understood. As recombination depends on a high level of nucleotide sequence homology, we postulate that the frequency of somatic variants derived from mitotic recombination should be diminished in progeny from crosses between strains of mice in which nucleotide sequences have diverged. Here we report that mitotic recombination is suppressed, to various degrees in different tissues, in hybrids of distantly related mouse strains. Reintroduction of greater chromosomal homology by backcrossing restores mitotic recombination in offspring. Thus, chromosomal divergence inhibits mitotic recombination and, consequently, may act as a modifier of cancer susceptibility by limiting the rate of LOH. The suppression of mitotic recombination in some F1 hybrids in which meiotic recombination persists indicates that these processes are differentially affected by chromosomal divergence.     

Isolation of chromosome 21-specific yeast artificial chromosomes from a total human genome library.

A new approach for the isolation of chromosome-specific subsets from a human genomic yeast artificial chromosome (YAC) library is described. It is based on the hybridization with an Alu polymerase chain reaction (PCR) probe. We screened a 1.5 genome equivalent YAC library of megabase insert size with Alu PCR products amplified from hybrid cell lines containing human chromosome 21, and identified a subset of 63 clones representative of this chromosome. The majority of clones were assigned to chromosome 21 by the presence of specific STSs and in situ hybridization. Twenty-nine of 36 STSs that we tested were detected in the subset, and a contig spanning 20 centimorgans in the genetic map and containing 8 STSs in 4 YACs was identified. The proposed approach can greatly speed efforts to construct physical maps of the human genome.     

Fine-scale recombination patterns differ between chimpanzees and humans

Recombination rates seem to vary extensively along the human genome. Pedigree analysis suggests that rates vary by an order of magnitude when measured at the megabase scale, and at a finer scale, sperm typing studies point to the existence of recombination hotspots. These are short regions (1-2 kb) in which recombination rates are 10-1,000 times higher than the background rate. Less is known about how recombination rates change over time. Here we determined to what degree recombination rates are conserved among closely related species by estimating recombination rates from 14 Mb of linkage disequilibrium data in central chimpanzee and human populations. The results suggest that recombination hotspots are not conserved between the two species and that recombination rates in larger (50 kb) genomic regions are only weakly conserved. Therefore, the recombination landscape has changed markedly between the two species.     

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.     

Toward simpler and faster genome-wide mutagenesis in mice

Here we describe a practical Cre-loxP and piggyBac transposon-based mutagenesis strategy to systematically mutate coding sequences and/or the vast noncoding regions of the mouse genome for large-scale functional genomic analysis. To illustrate this approach, we first created loxP-containing loss-of-function alleles in the protocadherin alpha, beta and gamma gene clusters (Pcdha, Pcdhb and Pcdhg). Using these alleles, we show that, under proper guidance, Cre-loxP site-specific recombination can mediate efficient trans-allelic recombination in vivo, facilitating the generation of large germline deletions and duplications including deletions of Pcdha, and Pcdha to Pcdhb, simply by breeding (that is, at frequencies of 5.5%-21.6%). The same breeding method can also generate designed germline translocations between nonhomologous chromosomes at unexpected frequencies of greater than 1%. By incorporating a piggyBac transposon to insert and to distribute loxP sites randomly throughout the mouse genome, we present a simple but comprehensive method for generating genome-wide deletions and duplications, in addition to insertional loss-of-function and conditional rescue alleles, again simply by breeding.     

A complete BAC-based physical map of the Arabidopsis thaliana genome.

Arabidopsis thaliana is a small flowering plant that serves as the major model system in plant molecular genetics. The efforts of many scientists have produced genetic maps that provide extensive coverage of the genome (http://genome-www. stanford.edu/Arabidopsis/maps.html). Recently, detailed YAC, BAC, P1 and cosmid-based physical maps (that is, representations of genomic regions as sets of overlapping clones of corresponding libraries) have been established that extend over wide genomic areas ranging from several hundreds of kilobases to entire chromosomes. These maps provide an entry to gain deeper insight into the A. thaliana genome structure. A. thaliana has been chosen as the subject of the first large-scale project intended to determine the full genome sequence of a plant. This sequencing project, together with the increasing interest in map-based gene cloning, has highlighted the requirement for a complete and accurate physical map of this plant species. To supply the scientific community with a high-quality resource, we present here a complete physical map of A. thaliana using essentially the IGF BAC library. The map consists of 27 contigs that cover the entire genome, except for the presumptive centromeric regions, nucleolar organization regions (NOR) and telomeric areas. This is the first reported map of a complex organism based entirely on BAC clones and it represents the most homogeneous and complete physical map established to date for any plant genome. Furthermore, the analysis performed here serves as a model for an efficient physical mapping procedure using BAC clones that can be applied to other complex genomes.     

Cancer predisposition caused by elevated mitotic recombination in Bloom mice

Bloom syndrome is a disorder associated with genomic instability that causes affected people to be prone to cancer. Bloom cell lines show increased sister chromatid exchange, yet are proficient in the repair of various DNA lesions. The underlying cause of this disease are mutations in a gene encoding a RECQ DNA helicase. Using embryonic stem cell technology, we have generated viable Bloom mice that are prone to a wide variety of cancers. Cell lines from these mice show elevations in the rates of mitotic recombination. We demonstrate that the increased rate of loss of heterozygosity (LOH) resulting from mitotic recombination in vivo constitutes the underlying mechanism causing tumour susceptibility in these mice.     

A Robertsonian translocation suppresses a somatic recombination pathway to loss of heterozygosity

In mammals, loss of APC/Apc gatekeeper function initiates intestinal tumorigenesis. Several different mechanisms have been shown or proposed to mediate functional loss of APC/Apc: mutation in APC/Apc, non-disjunction, homologous somatic recombination and epigenetic silencing. The demonstration that, in the C57BL/6 (B6) Apc(Min/+) mouse model of inherited intestinal cancer, loss of Apc function can occur by loss of heterozygosity (LOH) through somatic recombination between homologs presents an opportunity to search for polymorphisms in the homologous somatic recombination pathway. We report that the Robertsonian translocation Rb(7.18)9Lub (Rb9) suppresses the multiplicity of intestinal adenomas in this mouse model. As the copy number of Rb9 increases, the association with the interphase nucleolus of the rDNA repeats centromeric to the Apc locus on Chromosome 18 is increasingly disrupted. Our analysis shows that homologous somatic recombination is the principal pathway for LOH in adenomas in B6 Apc(Min/+) mice. These studies provide additional evidence that neoplastic growth can initiate in the complete absence of canonical genomic instability.     

Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors

Aberrant methylation of CpG islands and genomic deletion are two predominant mechanisms of gene inactivation in tumorigenesis, but the extent to which they interact is largely unknown. The lack of an integrated approach to study these mechanisms has limited the understanding of tumor genomes and cancer genes. Restriction landmark genomic scanning (RLGS; ref. 1) is useful for global analysis of aberrant methylation of CpG islands, but has not been amenable to alignment with deletion maps because the identity of most RLGS fragments is unknown. Here, we determined the nucleotide sequence and exact chromosomal position of RLGS fragments throughout the genome using the whole chromosome of origin of the fragments and in silico restriction digestion of the human genome sequence. To study the interaction of these gene-inactivation mechanisms in primary brain tumors, we integrated RLGS-based methylation analysis with high-resolution deletion maps from microarray-based comparative genomic hybridization (array CGH; ref. 3). Certain subsets of gene-associated CpG islands were preferentially affected by convergent methylation and deletion, including genes that exhibit tumor-suppressor activity, such as CISH1 (encoding SOCS1; ref. 4), as well as genes such as COE3 that have been missed by traditional non-integrated approaches. Our results show that most aberrant methylation events are focal and independent of deletions, and the rare convergence of these mechanisms can pinpoint biallelic gene inactivation without the use of positional cloning.     

A YAC-based physical map of the mouse genome.

A physical map of the mouse genome is an essential tool for both positional cloning and genomic sequencing in this key model system for biomedical research. Indeed, the construction of a mouse physical map with markers spaced at an average interval of 300 kb is one of the stated goals of the Human Genome Project. Here we report the results of a project at the Whitehead Institute/MIT Center for Genome Research to construct such a physical map of the mouse. We built the map by screening sequenced-tagged sites (STSs) against a large-insert yeast artificial chromosome (YAC) library and then integrating the STS-content information with a dense genetic map. The integrated map shows the location of 9,787 loci, providing landmarks with an average spacing of approximately 300 kb and affording YAC coverage of approximately 92% of the mouse genome. We also report the results of a project at the MRC UK Mouse Genome Centre targeted at chromosome X. The project produced a YAC-based map containing 619 loci (with 121 loci in common with the Whitehead map and 498 additional loci), providing especially dense coverage of this sex chromosome. The YAC-based physical map directly facilitates positional cloning of mouse mutations by providing ready access to most of the genome. More generally, use of this map in addition to a newly constructed radiation hybrid (RH) map provides a comprehensive framework for mouse genomic studies.     

Germline rates of de novo meiotic deletions and duplications causing several genomic disorders

Meiotic recombination between highly similar duplicated sequences (nonallelic homologous recombination, NAHR) generates deletions, duplications, inversions and translocations, and it is responsible for genetic diseases known as 'genomic disorders', most of which are caused by altered copy number of dosage-sensitive genes. NAHR hot spots have been identified within some duplicated sequences. We have developed sperm-based assays to measure the de novo rate of reciprocal deletions and duplications at four NAHR hot spots. We used these assays to dissect the relative rates of NAHR between different pairs of duplicated sequences. We show that (i) these NAHR hot spots are specific to meiosis, (ii) deletions are generated at a higher rate than their reciprocal duplications in the male germline and (iii) some of these genomic disorders are likely to have been underascertained clinically, most notably that resulting from the duplication of 7q11, the reciprocal of the deletion causing Williams-Beuren syndrome.     

The extent of linkage disequilibrium in Arabidopsis thaliana.

Linkage disequilibrium (LD), the nonrandom occurrence of alleles in haplotypes, has long been of interest to population geneticists. Recently, the rapidly increasing availability of genomic polymorphism data has fueled interest in LD as a tool for fine-scale mapping, in particular for human disease loci. The chromosomal extent of LD is crucial in this context, because it determines how dense a map must be for associations to be detected and, conversely, limits how finely loci may be mapped. Arabidopsis thaliana is expected to harbor unusually extensive LD because of its high degree of selfing. Several polymorphism studies have found very strong LD within individual loci, but also evidence of some recombination. Here we investigate the pattern of LD on a genomic scale and show that in global samples, LD decays within approximately 1 cM, or 250 kb. We also show that LD in local populations may be much stronger than that of global populations, presumably as a result of founder events. The combination of a relatively high level of polymorphism and extensive haplotype structure bodes well for developing a genome-wide LD map in A. thaliana.     

Systematic generation of high-resolution deletion coverage of the Drosophila melanogaster genome

In fruit fly research, chromosomal deletions are indispensable tools for mapping mutations, characterizing alleles and identifying interacting loci. Most widely used deletions were generated by irradiation or chemical mutagenesis. These methods are labor-intensive, generate random breakpoints and result in unwanted secondary mutations that can confound phenotypic analyses. Most of the existing deletions are large, have molecularly undefined endpoints and are maintained in genetically complex stocks. Furthermore, the existence of haplolethal or haplosterile loci makes the recovery of deletions of certain regions exceedingly difficult by traditional methods, resulting in gaps in coverage. Here we describe two methods that address these problems by providing for the systematic isolation of targeted deletions in the D. melanogaster genome. The first strategy used a P element-based technique to generate deletions that closely flank haploinsufficient genes and minimize undeleted regions. This deletion set has increased overall genomic coverage by 5-7%. The second strategy used FLP recombinase and the large array of FRT-bearing insertions described in the accompanying paper to generate 519 isogenic deletions with molecularly defined endpoints. This second deletion collection provides 56% genome coverage so far. The latter methodology enables the generation of small custom deletions with predictable endpoints throughout the genome and should make their isolation a simple and routine task.     

The value of data

Data citation and the derivation of semantic constructs directly from datasets have now both found their place in scientific communication. The social challenge facing us is to maintain the value of traditional narrative publications and their relationship to the datasets they report upon while at the same time developing appropriate metrics for citation of data and data constructs.     

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.     

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.