Characterization of whole genome amplified (WGA) DNA for use in genotyping assay development.

BACKGROUND: Genotyping assays often require substantial amounts of DNA. To overcome the problem of limiting amounts of available DNA, Whole Genome Amplification (WGA) methods have been developed. The multiple displacement amplification (MDA) method using Φ29 polymerase has become the preferred choice due to its high processivity and low error rate. However, the uniformity and fidelity of the amplification process across the genome has not been extensively characterized. RESULTS: To assess amplification uniformity, we used array-based comparative genomic hybridization (aCGH) to evaluate DNA copy number variations (CNVs) in DNAs amplified by two MDA kits: GenomiPhi and REPLI-g. The Agilent Human CGH array containing nearly one million probes was used in this study together with DNAs from a normal subject and 2 cystic fibrosis (CF) patients. Each DNA sample was amplified 4 independent times and compared to its native unamplified DNA. Komogorov distances and Phi correlations showed a high consistency within each sample group. Less than 2% of the probes showed more than 2-fold CNV introduced by the amplification process. The two amplification kits, REPLI-g and GenomiPhi, generate very similar amplified DNA samples despite the differences between the unamplified and amplified DNA samples. The results from aCGH analysis indicated that there were no obvious CNVs in the CFTR gene region due to WGA when compared to unamplified DNA. This was confirmed by quantitative real-time PCR copy number assays at 10 locations within the CFTR gene. DNA sequencing analysis of a 2-kb region within the CFTR gene showed no mutations introduced by WGA. CONCLUSION: The relatively high uniformity and consistency of the WGA process, coupled with the low replication error rate, suggests that WGA DNA may be suitable for accurate genotyping. Regions of the genome that were consistently under-amplified were found to contain higher than average GC content. Because of the consistent differences between the WGA DNA and the native unamplified DNA, characterization of the genomic region of interest, as described here, will be necessary to ensure the reliability of genotyping results from WGA DNA.

Genome-wide expression profiling of urinary bladder implicates desmosomal and cytoskeletal dysregulation in the bladder exstrophy-epispadias complex.

The bladder exstrophy-epispadias complex (BEEC) represents a spectrum of urological abnormalities where part, or all, of the distal urinary tract fails to close during development, becoming exposed on the outer abdominal wall. While the etiology of BEEC remains unknown, strong evidence exists that genetic factors are implicated. To understand the pathways regulating embryonic bladder development and to identify high-probability BEEC candidate genes, we conducted a genome-wide expression profiling (GWEP) study using normal and exstrophic human urinary bladders and human and mouse embryologic bladder-precursor tissues. We identified 162 genes differentially expressed in both embryonic and postnatal human samples. Pathway analysis of these genes revealed 11 biological networks with top functions related to skeletal and muscular system development, cellular assembly and development, organ morphology, or connective tissue disorders. The two most down-regulated genes desmin (DES, fold-change, -74.7) and desmuslin (DMN, fold-change, -53.0) are involved in desmosome mediated cell-cell adhesion and cytoskeletal architecture. Intriguingly, the sixth most overexpressed gene was desmoplakin (DSP, fold-change, +48.8), the most abundant desmosomal protein. We found 30% of the candidate genes to be directly associated with desmosome structure/function or cytoskeletal assembly, pointing to desmosomal and/or cytoskeletal deregulation as an etiologic factor for BEEC. Further findings indicate that p63, PERP, SYNPO2 and the Wnt pathway may also contribute to BEEC etiology. This study provides the first expression profile of urogenital genes during bladder development and points to the high-probability candidate genes for BEEC.

Hepatitis B virus replication induces methylation of both host and viral DNA.

Control of viral replication is a major therapeutic goal to reduce morbidity and mortality from chronic hepatitis B virus (HBV) infection. Recently, methylation has been identified as a novel host defense mechanism, and methylation of viral DNA leads to downregulation of HBV gene expression. To better understand the mechanisms of HBV methylation, cell lines were exposed to HBV using a model system that mimics natural infection and the expression of host DNA methyltransferase genes (DNMTs) was measured. DNMT1, DNMT2, and DNMT3 were all significantly upregulated in response to HBV. DNMT3 was further studied because of its known role in the de novo methylation of DNA. Cotransfection experiments with full-length HBV and DNMT3 led to the downregulation of viral protein and pregenomic RNA production. To investigate whether the upregulation of DNMTs could also have an effect on the methylation of host DNA, cell lines were exposed to HBV in two independent model systems, one that mimics natural infection and a second model with temporary transfection. Host DNA methylation was measured by DNA microarray analysis. Increased methylation of host CpG islands was detected in both experimental systems. Two CpG islands, corresponding to genes SUFU and TIRAP, were selected, and the downregulation of these genes in hepatocellular carcinomas was confirmed. In conclusion, hepatocytes respond to HBV infection by upregulating DNMTs. The DNMTs methylate viral DNA, leading to decreased viral gene expression and decreased viral replication. However, virus-induced overexpression of DNMTs also leads to methylation of host CpG islands.

The incredible shrinking world of DNA microarrays.

The efficacy of microarrays in examining gene expression, gene and genome structure, protein-DNA interactions, whole-genome similarities and differences, microRNA expression, methylation (and more) is no longer in question. It is a fast-developing, cutting edge technology that has grown up along with massive sequence databases and is likely to become part of everyday patient care. Many advances have recently expanded the power and utility of microarrays; among them is our development of a new array tiling technique that dramatically increases the scope of coverage of an oligonucleotide tiling array without substantially increasing its cost.

Fibrolamellar carcinomas show overexpression of genes in the RAS, MAPK, PIK3, and xenobiotic degradation pathways.

Fibrolamellar carcinomas (FLC) are a rare type of primary hepatocellular carcinoma found in younger individuals. FLC are known to have relatively few consistent chromosomal alterations, although a gain of chromosome 1q has been reported. The gene expression of 4 FLC (2 primary FLC and 2 metastatic deposits) were studied using Affymetrix DNA microarray technology (Santa Clara, CA). Selected genes were confirmed by real-time polymerase chain reaction. Relatively few genes were significantly overexpressed-447 genes, case 1; 1298 genes, case 2-corresponding to approximately 0.8% and 2.3%, respectively, of the 56000 transcripts present in the arrays. Of these, 155 genes were overexpressed simultaneously by both tumors. The number of significantly overexpressed genes more than doubled in the 2 metastatic deposits (2777 and 2855 genes compared with 1298 in the primary tumor). Proteins involved in the RAS, MAPK, PIK3, and xenobiotic degradation pathways were commonly overexpressed. Because chromosome 1q is thought to contain an important oncogene, additional attention was focused on this region. Of 114 total genes found overexpressed in common among all primary and metastatic tumors, 11 of 114 genes were located on chromosome 1q: ARF1, CD46, CNIH4, ENSA, FH, NICE-3, PSMB4, RGS2, RGS5, TIMM17A, and UFC1. Primary FLC show overexpression of genes involved in the RAS, MAPK, PIK3, and xenobiotic degradation pathways. Eleven common genes were consistently overexpressed on chromosome 1q among all tumors and metastases and warrant further study as potential oncogenes.

Transposon insertion site profiling chip (TIP-chip).

Mobile elements are important components of our genomes, with diverse and significant effects on phenotype. Not only can transposons inactivate genes by direct disruption and shuffle the genome through recombination, they can also alter gene expression subtly or powerfully. Currently active transposons are highly polymorphic in host populations, including, among hundreds of others, L1 and Alu elements in humans and Ty1 elements in yeast. For this reason, we wished to develop a simple genome-wide method for identifying all transposons in any given sample. We have designed a transposon insertion site profiling chip (TIP-chip), a microarray intended for use as a high-throughput technique for mapping transposon insertions. By selectively amplifying transposon flanking regions and hybridizing them to the array, we can locate all transposons present in a sample. We have tested the TIP-chip extensively to map Ty1 retrotransposon insertions in yeast and have achieved excellent results in two laboratory strains as well as in evolved Ty1 high-copy strains. We are able to identify all of the theoretically detectable transposons in the FY2 lab strain, with essentially no false positives. In addition, we mapped many new transposon copies in the high-copy Ty1 strain and determined its Ty1 insertion pattern.

Stacking the deck: double-tiled DNA microarrays.

Microarrays-high-throughput platforms for analyzing the gene expression and features of total genomic DNA, among other applications-are gaining in popularity as researchers discover ever more uses for their unbiased and broad feature sets. At present, microarray analyses are limited by the number of individual features that can be placed on each array. Here we describe a double-tiling method that significantly increases the number of sequences present on an array, and we show that successful transcriptional profiling is possible and straightforward with such arrays. With this method, we and others can save money and precious samples by using fewer arrays to cover a region, or can carry out investigations at significantly higher resolution without incurring prohibitive costs or increasing the amount of sample required for the experiment.

Role of cocaine- and amphetamine-regulated transcript in estradiol-mediated neuroprotection.

Estrogen reduces brain injury after experimental cerebral ischemia in part through a genomic mechanism of action. Using DNA microarrays, we analyzed the genomic response of the brain to estradiol, and we identified a transcript, cocaine- and amphetamine-regulated transcript (CART), that is highly induced in the cerebral cortex by estradiol under ischemic conditions. Using in vitro and in vivo models of neural injury, we confirmed and characterized CART mRNA and protein up-regulation by estradiol in surviving neurons, and we demonstrated that i.v. administration of a rat CART peptide is protective against ischemic brain injury in vivo. We further demonstrated binding of cAMP response element (CRE)-binding protein to a CART promoter CRE site in ischemic brain and rapid activation by CART of ERK in primary cultured cortical neurons. The findings suggest that CART is an important player in estrogen-mediated neuroprotection and a potential therapeutic agent for stroke and other neurodegenerative diseases.

Multiple-laboratory comparison of microarray platforms.

Microarray technology is a powerful tool for measuring RNA expression for thousands of genes at once. Various studies have been published comparing competing platforms with mixed results: some find agreement, others do not. As the number of researchers starting to use microarrays and the number of cross-platform meta-analysis studies rapidly increases, appropriate platform assessments become more important. Here we present results from a comparison study that offers important improvements over those previously described in the literature. In particular, we noticed that none of the previously published papers consider differences between labs. For this study, a consortium of ten laboratories from the Washington, DC-Baltimore, USA, area was formed to compare data obtained from three widely used platforms using identical RNA samples. We used appropriate statistical analysis to demonstrate that there are relatively large differences in data obtained in labs using the same platform, but that the results from the best-performing labs agree rather well.

TGF-beta-dependent pathogenesis of mitral valve prolapse in a mouse model of Marfan syndrome.

Mitral valve prolapse (MVP) is a common human phenotype, yet little is known about the pathogenesis of this condition. MVP can occur in the context of genetic syndromes, including Marfan syndrome (MFS), an autosomal-dominant connective tissue disorder caused by mutations in fibrillin-1. Fibrillin-1 contributes to the regulated activation of the cytokine TGF-beta, and enhanced signaling is a consequence of fibrillin-1 deficiency. We thus hypothesized that increased TGF-beta signaling may contribute to the multisystem pathogenesis of MFS, including the development of myxomatous changes of the atrioventricular valves. Mitral valves from fibrillin-1-deficient mice exhibited postnatally acquired alterations in architecture that correlated both temporally and spatially with increased cell proliferation, decreased apoptosis, and excess TGF-beta activation and signaling. In addition, TGF-beta antagonism in vivo rescued the valve phenotype, suggesting a cause and effect relationship. Expression analyses identified increased expression of numerous TGF-beta-related genes that regulate cell proliferation and survival and plausibly contribute to myxomatous valve disease. These studies validate a novel, genetically engineered murine model of myxomatous changes of the mitral valve and provide critical insight into the pathogenetic mechanism of such changes in MFS and perhaps more common nonsyndromic variants of mitral valve disease.

Regulation of c-Met-dependent gene expression by PTEN.

Receptor tyrosine kinases (RTK) and the tumor suppressor PTEN co-regulate oncogenic cell signaling pathways. How these interactions influence gene transcription is inadequately understood. We used expression microarrays to investigate the effects of PTEN on gene expression changes caused by activating c-Met in human glioblastoma cells. c-Met activation by scatter factor/hepatocyte growth factor (SF/HGF) altered the expression of 27-fold more genes in PTEN-null U-373MG cells than in PTEN homozygous primary normal human astrocytes (523 vs 19 genes). Restoring wt-PTEN in U-373MG cells dramatically altered patterns of c-Met regulated gene expression. This effect was varied depending on the specific gene in question. PTEN reduced the number of c-Met regulated transcripts from 931 to 502, decreased the relative number of genes upregulated by c-Met from 46 to 25%, and increased the relative number of downregulated genes from 54 to 75%. PTEN and c-Met co-regulated many genes involved in cell growth regulation such as oncogenes, growth factors, transcription factors, and constituents of the ubiquitin pathway. c-Met activation in PTEN-null (but not PTEN reconstituted) cells led to upregulation of the EGFR agonist TGFalpha and subsequently to EGFR activation. Using PTEN mutants, we found that PTEN's transcriptional effects were either lipid-phosphatase dependent, protein-phosphatase dependent, or phosphatase-independent. These results show that PTEN has critical and mechanistically complex effects on RTK-regulated gene transcription. These findings expand our understanding of tumor promoter/suppressor inter-relationships and downstream transcriptional effects of PTEN loss and c-Met overexpression in malignant gliomas.

The role of CTLs in persistent viral infection: cytolytic gene expression in CD8+ lymphocytes distinguishes between individuals with a high or low proviral load of human T cell lymphotropic virus type 1.

The proviral load in human T cell lymphotropic virus type 1 (HTLV-1) infection is typically constant in each infected host, but varies by >1000-fold between hosts and is strongly correlated with the risk of HTLV-1-associated inflammatory disease. However, the factors that determine an individual's HTLV-1 proviral load remain uncertain. Experimental evidence from studies of host genetics, viral genetics, and lymphocyte function and theoretical considerations suggest that a major determinant of the equilibrium proviral load is the CD8+ T cell response to HTLV-1. In this study, we tested the hypothesis that the gene expression profile in circulating CD8+ and CD4+ lymphocytes distinguishes between individuals with a low proviral load of HTLV-1 and those with a high proviral load. We show that circulating CD8+ lymphocytes from individuals with a low HTLV-1 proviral load overexpressed a core group of nine genes with strong functional coherence: eight of the nine genes encode granzymes or other proteins involved in cell-mediated lysis or Ag recognition. We conclude that successful suppression of the HTLV-1 proviral load is associated with strong cytotoxic CD8+ lymphocyte activity in the peripheral blood.

Nonsense surveillance regulates expression of diverse classes of mammalian transcripts and mutes genomic noise.

Premature termination codons induce rapid transcript degradation in eukaryotic cells through nonsense-mediated mRNA decay (NMD). This pathway can modulate phenotypes arising from nonsense or frameshift mutations, but little is known about the physiologic role of NMD in higher eukaryotes. To address this issue, we examined expression profiles in mammalian cells depleted of Rent1 (also called hUpf1), a factor essential for NMD. Upregulated transcripts included those with upstream open reading frames in the 5' untranslated region, alternative splicing that introduces nonsense codons or frameshifts, introns in the 3' untranslated region or selenocysteine codons. Transcripts derived from ancient transposons and endogenous retroviruses were also upregulated. These RNAs are unified by the presence of a spliced intron at least 50 nucleotides downstream of a termination codon, a context sufficient to initiate NMD. Consistent with direct regulation by NMD, representative upregulated transcripts decayed more slowly in cells deficient in NMD. In addition, inhibition of NMD induced by amino acid starvation upregulated transcripts that promote amino acid homeostasis. These results document that nonsense surveillance is a crucial post-transcriptional regulatory event that influences the expression of broad classes of physiologic transcripts, has been functionally incorporated into essential homeostatic mechanisms and suppresses expression of evolutionary remnants.

Microarray and serial analysis of gene expression analyses identify known and novel transcripts overexpressed in hematopoietic stem cells.

The human CD34(+)/CD38(-)/Lin(-) cell subset, comprising approximately 1-10% of the CD34(+) cell population, contains few of the less primitive hematopoietic (lineage-committed) progenitor cells (HPCs) but most of the primitive in vivo engrafting (lympho-)hematopoietic stem cells (HSCs). We analyzed gene expression in CD34(+)/CD38(-)/Lin(-) cell populations isolated from normal human adult donor bone marrow, neonatal placental/umbilical cord blood, and mobilized adult donor peripheral blood stem-progenitor cells. As measured by Affymetrix microarrays, 4746 genes were expressed in CD34(+)/CD38(-)/Lin(-) cells from all three tissues. We also determined the transcriptomes of the stem cell-depleted, HPC-enriched CD34(+)/[CD38/Lin](++) cell population from each tissue. Comparison of CD34(+)/CD38(-)/Lin(-) (HSC-enriched) versus CD34(+)/[CD38/Lin](++) (HPC-enriched, HSC-depleted) cells from each tissue yielded 81 genes overrepresented and 90 genes underrepresented, common to all three of the CD34(+)/CD38(-)/Lin(-) cell populations. These transcripts, which are selectively expressed in HSCs from all three tissues, include a number of known genes (e.g., transcription factors, receptors, and signaling molecules) that might play roles in key functions (e.g., survival, self-renewal, differentiation, and/or migration/adhesion) of human HSCs. Many genes/transcripts of unknown function were also detected by microarray analysis. Serial analysis of gene expression of the bone marrow HSC and HPC populations confirmed expression of most of the overrepresented transcripts for which reliable serial analysis of gene expression tags were detected and additionally suggested that current microarrays do not detect as many as 30% of the transcripts expressed in HSCs, including a number of previously unknown transcripts. This work is a step toward full definition of the transcriptome of normal human HSCs and may identify new genes involved in leukemogenesis and cancer stem cells.