A genome-wide scalable SNP genotyping assay using microarray technology

Oligonucleotide probe arrays have enabled massively parallel analysis of gene expression levels from a single cDNA sample. Application of microarray technology to analyzing genomic DNA has been stymied by the sequence complexity of the entire human genome. A robust, single base-resolution direct genomic assay would extend the reach of microarray technology. We developed an array-based whole-genome genotyping assay that does not require PCR and enables effectively unlimited multiplexing. The assay achieves a high signal-to-noise ratio by combining specific hybridization of picomolar concentrations of whole genome-amplified DNA to arrayed probes with allele-specific primer extension and signal amplification. As proof of principle, we genotyped several hundred previously characterized SNPs. The conversion rate, call rate and accuracy were comparable to those of high-performance PCR-based genotyping assays.     

Evidence for an instructive mechanism of de novo methylation in cancer cells

DNA methylation has a role in the regulation of gene expression during normal mammalian development but can also mediate epigenetic silencing of CpG island genes in cancer and other diseases. Many individual genes (including tumor suppressors) have been shown to undergo de novo methylation in specific tumor types, but the biological logic inherent in this process is not understood. To decipher this mechanism, we have adopted a new approach for detecting CpG island DNA methylation that can be used together with microarray technology. Genome-wide analysis by this technique demonstrated that tumor-specific methylated genes belong to distinct functional categories, have common sequence motifs in their promoters and are found in clusters on chromosomes. In addition, many are already repressed in normal cells. These results are consistent with the hypothesis that cancer-related de novo methylation may come about through an instructive mechanism.     

Duplicate genes increase gene expression diversity within and between species

Using microarray gene expression data from several Drosophila species and strains, we show that duplicated genes, compared with single-copy genes, significantly increase gene expression diversity during development. We show further that duplicate genes tend to cause expression divergences between Drosophila species (or strains) to evolve faster than do single-copy genes. This conclusion is also supported by data from different yeast strains.     

Widespread aneuploidy revealed by DNA microarray expression profiling

Expression profiling using DNA microarrays holds great promise for a variety of research applications, including the systematic characterization of genes discovered by sequencing projects. To demonstrate the general usefulness of this approach, we recently obtained expression profiles for nearly 300 Saccharomyces cerevisiae deletion mutants. Approximately 8% of the mutants profiled exhibited chromosome-wide expression biases, leading to spurious correlations among profiles. Competitive hybridization of genomic DNA from the mutant strains and their isogenic parental wild-type strains showed they were aneuploid for whole chromosomes or chromosomal segments. Expression profile data published by several other laboratories also suggest the use of aneuploid strains. In five separate cases, the extra chromosome harboured a close homologue of the deleted gene; in two cases, a clear growth advantage for cells acquiring the extra chromosome was demonstrated. Our results have implications for interpreting whole-genome expression data, particularly from cells known to suffer genomic instability, such as malignant or immortalized cells.     

RNA sequencing shows no dosage compensation of the active X-chromosome

Mammalian cells from both sexes typically contain one active X chromosome but two sets of autosomes. It has previously been hypothesized that X-linked genes are expressed at twice the level of autosomal genes per active allele to balance the gene dose between the X chromosome and autosomes (termed 'Ohno's hypothesis'). This hypothesis was supported by the observation that microarray-based gene expression levels were indistinguishable between one X chromosome and two autosomes (the X to two autosomes ratio (X:AA) ~1). Here we show that RNA sequencing (RNA-Seq) is more sensitive than microarray and that RNA-Seq data reveal an X:AA ratio of ~0.5 in human and mouse. In Caenorhabditis elegans hermaphrodites, the X:AA ratio reduces progressively from ~1 in larvae to ~0.5 in adults. Proteomic data are consistent with the RNA-Seq results and further suggest the lack of X upregulation at the protein level. Together, our findings reject Ohno’s hypothesis, necessitating a major revision of the current model of dosage compensation in the evolution of sex chromosomes.     

Bayesian method to predict individual SNP genotypes from gene expression data

RNA profiling can be used to capture the expression patterns of many genes that are associated with expression quantitative trait loci (eQTLs). Employing published putative cis eQTLs, we developed a Bayesian approach to predict SNP genotypes that is based only on RNA expression data. We show that predicted genotypes can accurately and uniquely identify individuals in large populations. When inferring genotypes from an expression data set using eQTLs of the same tissue type (but from an independent cohort), we were able to resolve 99% of the identities of individuals in the cohort at P(adjusted) ≤ 1 × 10(-5). When eQTLs derived from one tissue were used to predict genotypes using expression data from a different tissue, the identities of 90% of the study subjects could be resolved at P(adjusted) ≤ 1 × 10(-5). We discuss the implications of deriving genotypic information from RNA data deposited in the public domain.     

An epi-allelic series of p53 hypomorphs created by stable RNAi produces distinct tumor phenotypes in vivo

The application of RNA interference (RNAi) to mammalian systems has the potential to revolutionize genetics and produce novel therapies. Here we investigate whether RNAi applied to a well-characterized gene can stably suppress gene expression in hematopoietic stem cells and produce detectable phenotypes in mice. Deletion of the Trp53 tumor suppressor gene greatly accelerates Myc-induced lymphomagenesis, resulting in highly disseminated disease. To determine whether RNAi suppression of Trp53 could produce a similar phenotype, we introduced several Trp53 short hairpin RNAs (shRNAs) into hematopoietic stem cells derived from E(mu)-Myc transgenic mice, and monitored tumor onset and overall pathology in lethally irradiated recipients. Different Trp53 shRNAs produced distinct phenotypes in vivo, ranging from benign lymphoid hyperplasias to highly disseminated lymphomas that paralleled Trp53-/- lymphomagenesis in the E(mu)-Myc mouse. In all cases, the severity and type of disease correlated with the extent to which specific shRNAs inhibited p53 activity. Therefore, RNAi can stably suppress gene expression in stem cells and reconstituted organs derived from those cells. In addition, intrinsic differences between individual shRNA expression vectors targeting the same gene can be used to create an 'epi-allelic series' for dissecting gene function in vivo.     

Loss of genomic methylation causes p53-dependent apoptosis and epigenetic deregulation

Cytosine methylation of mammalian DNA is essential for the proper epigenetic regulation of gene expression and maintenance of genomic integrity. To define the mechanism through which demethylated cells die, and to establish a paradigm for identifying genes regulated by DNA methylation, we have generated mice with a conditional allele for the maintenance DNA methyltransferase gene Dnmt1. Cre-mediated deletion of Dnmt1 causes demethylation of cultured fibroblasts and a uniform p53-dependent cell death. Mutational inactivation of Trp53 partially rescues the demethylated fibroblasts for up to five population doublings in culture. Oligonucleotide microarray analysis showed that up to 10% of genes are aberrantly expressed in demethylated fibroblasts. Our results demonstrate that loss of Dnmt1 causes cell-type-specific changes in gene expression that impinge on several pathways, including expression of imprinted genes, cell-cycle control, growth factor/receptor signal transduction and mobilization of retroelements.     

Mining the human genome using microarrays of open reading frames

To test the hypothesis that the human genome project will uncover many genes not previously discovered by sequencing of expressed sequence tags (ESTs), we designed and produced a set of microarrays using probes based on open reading frames (ORFs) in 350 Mb of finished and draft human sequence. Our approach aims to identify all genes directly from genomic sequence by querying gene expression. We analysed genomic sequence with a suite of ORF prediction programs, selected approximately one ORF per gene, amplified the ORFs from genomic DNA and arrayed the amplicons onto treated glass slides. Of the first 10,000 arrayed ORFs, 31% are completely novel and 29% are similar, but not identical, to sequences in public databases. Approximately one-half of these are expressed in the tissues we queried by microarray. Subsequent verification by other techniques confirmed expression of several of the novel genes. Expressed sequence tags (ESTs) have yielded vast amounts of data, but our results indicate that many genes in the human genome will only be found by genomic sequencing.     

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.     

Epistasis analysis with global transcriptional phenotypes

Classical epistasis analysis can determine the order of function of genes in pathways using morphological, biochemical and other phenotypes. It requires knowledge of the pathway's phenotypic output and a variety of experimental expertise and so is unsuitable for genome-scale analysis. Here we used microarray profiles of mutants as phenotypes for epistasis analysis. Considering genes that regulate activity of protein kinase A in Dictyostelium, we identified known and unknown epistatic relationships and reconstructed a genetic network with microarray phenotypes alone. This work shows that microarray data can provide a uniform, quantitative tool for large-scale genetic network analysis.     

Deficiency in short-chain fatty acid beta-oxidation affects theta oscillations during sleep

In rodents, the electroencephalogram (EEG) during paradoxical sleep and exploratory behavior is characterized by theta oscillations. Here we show that a deficiency in short-chain acyl-coenzyme A dehydrogenase (encoded by Acads) in mice causes a marked slowing in theta frequency during paradoxical sleep only. We found Acads expression in brain regions involved in theta generation, notably the hippocampus. Microarray analysis of gene expression in mice with mutations in Acads indicates overexpression of Glo1 (encoding glyoxylase 1), a gene involved in the detoxification of metabolic by-products. Administration of acetyl-L-carnitine (ALCAR) to mutant mice significantly recovers slow theta and Glo1 overexpression. Thus, an underappreciated metabolic pathway involving fatty acid beta-oxidation also regulates theta oscillations during sleep.