Novel transcribed regions in the human genome.

We have used genomic tiling arrays to identify transcribed regions throughout the human genome. Analysis of the mapping results of RNA isolated from five cell/tissue types, NB4 cells, NB4 cells treated with retinoic acid (RA), NB4 cells treated with 12-O-tetradecanoylphorbol-13 acetate (TPA), neutrophils, and placenta, throughout the ENCODE region reveals a large number of novel transcribed regions. Interestingly, neutrophils exhibit a great deal of novel expression in several intronic regions. Comparison of the hybridization results of NB4 cells treated with different stimuli relative to untreated cells reveals that many new regions are expressed upon cell differentiation. One such region is the Hox locus, which contains a large number of novel regions expressed in a number of cell types. Analysis of the trinucleotide composition of the novel transcribed regions reveals that it is similar to that of known exons. These results suggest that many of the novel transcribed regions may have a functional role.

Transcriptional landscape of the human and fly genomes: nonlinear and multifunctional modular model of transcriptomes.

Regions of the genome not coding for proteins or not involved in cis-acting regulatory activities are frequently viewed as lacking in functional value. However, a number of recent large-scale studies have revealed significant regulated transcription of unannotated portions of a variety of plant and animal genomes, allowing a new appreciation of the widespread transcription of large portions of the genome. High-resolution mapping of the sites of transcription of the human and fly genomes has provided an alternative picture of the extent and organization of transcription and has offered insights for biological functions of some of the newly identified unannotated transcripts. Considerable portions of the unannotated transcription observed are developmental or cell-type-specific parts of protein-coding transcripts, often serving as novel, alternative 5' transcriptional start sites. These distal 5' portions are often situated at significant distances from the annotated gene and alternatively join with or ignore portions of other intervening genes to comprise novel unannotated protein-coding transcripts. These data support an interlaced model of the genome in which many regions serve multifunctional purposes and are highly modular in their utilization. This model illustrates the underappreciated organizational complexity of the genome and one of the functional roles of transcription from unannotated portions of the genome.

Prokaryotic RNA preparation methods useful for high density array analysis: comparison of two approaches.

High density oligonucleotide arrays have been used extensively for expression studies of eukaryotic organisms. We have designed a prokaryotic high density oligonucleotide array using the complete Escherichia coli genome sequence to monitor expression levels of all genes and intergenic regions in the genome. Because previously described methods for preparing labeled target nucleic acids are not useful for prokaryotic cell analysis using such arrays, a mRNA enrichment and direct labeling protocol was developed together with a cDNA synthesis protocol. The reproducibility of each labeling method was determined using high density oligonucleotide probe arrays as a read-out methodology and the expression results from direct labeling were compared to the expression results from the cDNA synthesis. About 50% of all annotated E.coli open reading frames are observed to be transcribed, as measured by both protocols, when the cells were grown in rich LB medium. Each labeling method individually showed a high degree of concordance in replica experiments (95 and 99%, respectively), but when each sample preparation method was compared to the other, approximately 32% of the genes observed to be expressed were discordant. However, both labeling methods can detect the same relative gene expression changes when RNA from IPTG-induced cells was labeled and compared to RNA from uninduced E.coli cells.

Reprogramming of the macrophage transcriptome in response to interferon-gamma and Mycobacterium tuberculosis: signaling roles of nitric oxide synthase-2 and phagocyte oxidase.

Macrophage activation determines the outcome of infection by Mycobacterium tuberculosis (Mtb). Interferon-gamma (IFN-gamma) activates macrophages by driving Janus tyrosine kinase (JAK)/signal transducer and activator of transcription-dependent induction of transcription and PKR-dependent suppression of translation. Microarray-based experiments reported here enlarge this picture. Exposure to IFN-gamma and/or Mtb led to altered expression of 25% of the monitored genome in macrophages. The number of genes suppressed by IFN-gamma exceeded the number of genes induced, and much of the suppression was transcriptional. Five times as many genes related to immunity and inflammation were induced than suppressed. Mtb mimicked or synergized with IFN-gamma more than antagonized its actions. Phagocytosis of nonviable Mtb or polystyrene beads affected many genes, but the transcriptional signature of macrophages infected with viable Mtb was distinct. Studies involving macrophages deficient in inducible nitric oxide synthase and/or phagocyte oxidase revealed that these two antimicrobial enzymes help orchestrate the profound transcriptional remodeling that underlies macrophage activation.

Flexible use of high-density oligonucleotide arrays for single-nucleotide polymorphism discovery and validation.

A method for identifying and validating single nucleotide polymorphisms (SNPs) with high-density oligonucleotide arrays without the need for locus-specific polymerase chain reactions (PCR) is described in this report. Genomic DNAs were divided into subsets with complexity of ~10 Mb by restriction enzyme digestion and gel-based fragment size resolution, ligated to a common adaptor, and amplified with one primer in a single PCR reaction. As a demonstration of this approach, a total of 124 SNPs were located in 190 kb of genomic sequences distributed across the entire human genome by hybridizing to high-density variant detection arrays (VDA). A set of independent validation experiments was conducted for these SNPs employing bead-based affinity selection followed by hybridization of the affinity-selected SNP-containing fragments to the same VDA that was used to identify the SNPs. A total of 98.7% (74/75) of these SNPs were confirmed using both DNA dideoxynucleotide sequencing and the VDA methodologies. With flexible sample preparation, high-density oligonucleotide arrays can be tailored for even larger scale genome-wide SNP discovery as well as validation.

Temporal gene regulation during HIV-1 infection of human CD4+ T cells.

CD4(+) T-cell depletion is a characteristic of human immunodeficiency virus type 1 (HIV-1) infection. In this study, modulation of mRNA expression of 6800 genes was monitored simultaneously at eight time points in a CD4(+) T-cell line (CEM-GFP) during HIV infection. The responses to infection included: (1) >30% decrease at 72 h after infection in overall host-cell production of monitored mRNA synthesis, with the replacement of host-cell mRNA by viral mRNA, (2) suppression of the expression of selected mitochondrial and DNA repair gene transcripts, (3) increased expression of the proapoptotic gene and its gene p53-induced product Bax, and (4) activation of caspases 2, 3, and 9. The intense HIV-1 transcription resulted in the repression of much cellular RNA expression and was associated with the induction of apoptosis of infected cells but not bystander cells. This choreographed host gene response indicated that the subversion of the cell transcriptional machinery for the purpose of HIV-1 replication is akin to genotoxic stress and represents a major factor leading to HIV-induced apoptosis.

Comparing genomes within the species Mycobacterium tuberculosis.

The study of genetic variability within natural populations of pathogens may provide insight into their evolution and pathogenesis. We used a Mycobacterium tuberculosis high-density oligonucleotide microarray to detect small-scale genomic deletions among 19 clinically and epidemiologically well-characterized isolates of M. tuberculosis. The pattern of deletions detected was identical within mycobacterial clones but differed between different clones, suggesting that this is a suitable genotyping system for epidemiologic studies. An analysis of genomic deletions among an extant population of pathogenic bacteria provided a novel perspective on genomic organization and evolution. Deletions are likely to contain ancestral genes whose functions are no longer essential for the organism's survival, whereas genes that are never deleted constitute the minimal mycobacterial genome. As the amount of genomic deletion increased, the likelihood that the bacteria will cause pulmonary cavitation decreased, suggesting that the accumulation of mutations tends to diminish their pathogenicity. Array-based comparative genomics is a promising approach to exploring molecular epidemiology, microbial evolution, and pathogenesis.

Detection of deleted genomic DNA using a semiautomated computational analysis of GeneChip data.

Genomic diversity within and between populations is caused by single nucleotide mutations, changes in repetitive DNA systems, recombination mechanisms, and insertion and deletion events. The contribution of these sources to diversity, whether purely genetic or of phenotypic consequence, can only be investigated if we have the means to quantitate and characterize diversity in many samples. With the advent of complete sequence characterization of representative genomes of different species, the possibility of developing protocols to screen for genetic polymorphism across entire genomes is actively being pursued. The large numbers of measurements such approaches yield demand that we pay careful attention to the numerical analysis of data. In this paper we present a novel application of an Affymetrix GeneChip to perform genome-wide screens for deletion polymorphism. A high-density oligonucleotide array formatted for mRNA expression and targeted at a fully sequenced 4.4-million-base pair Mycobacterium tuberculosis standard strain genome was adapted to compare genomic DNA. Hybridization intensities to 111,000 probe pairs (perfect complement and mismatch complement) were measured for genomic DNA from a clinical strain and from a vaccine organism. Because individual probe-pair hybridization intensities exhibit limited sensitivity/specificity characteristics to detect deletions, data-analytical methodology to exploit measurements from multiple probes in tandem locations across the genome was developed. The TSTEP (Tandem Set Terminal Extreme Probability) algorithm designed specifically to analyze the tandem hybridization measurements data was applied and shown to discover genomic deletions with high sensitivity. The TSTEP algorithm provides a foundation for similar efforts to characterize deletions in many hybridization measures in similar-sized and larger genomes. Issues relating to the design of genome content screening experiments and the implications of these methods for studying population genomics and the evolution of genomes are discussed.

Analysis of drug pharmacology towards predicting drug behavior by expression profiling using high-density oligonucleotide arrays.

An important aspect of the drug development process is prediction of efficacious and toxic side effects. Profiling of mRNA expression is a powerful approach to analyze the molecular phenotype of cells under various conditions, for example, in response to stimulation by compounds. We attempt to explore the approach of using expression profiling to identify patterns or fingerprints that are correlated with specific drug properties or behaviors. Identification of such expression patterns may also lead to revelation of the potential action mechanism of drugs and fingerprints indicative of certain drug efficacy or side effects. We describe here a strategy that was used to identify a set of genes whose differential expression pattern correlates with activation mode and target specificity of a specific group of drug compounds.

The transcriptional responses of respiratory epithelial cells to Bordetella pertussis reveal host defensive and pathogen counter-defensive strategies.

Bordetella pertussis, the causative agent of whooping cough, has many well-studied virulence factors and a characteristic clinical presentation. Despite this information, it is not clear how B. pertussis interaction with host cells leads to disease. In this study, we examined the interaction of B. pertussis with a human bronchial epithelial cell line (BEAS-2B) and measured host transcriptional profiles by using high-density DNA microarrays. The early transcriptional response to this pathogen is dominated by altered expression of cytokines, DNA-binding proteins, and NFkappaB-regulated genes. This previously unrecognized response to B. pertussis was modified in similar but nonidentical fashions by the antiinflammatory agents dexamethasone and sodium salicylate. Cytokine protein expression was confirmed, as was neutrophil chemoattraction. We show that B. pertussis induces mucin gene transcription by BEAS-2B cells then counters this defense by using mucin as a binding substrate. A set of genes is described for which the catalytic activity of pertussis toxin is both necessary and sufficient to regulate transcription. Host genomic transcriptional profiling, in combination with functional assays to evaluate subsequent biological events, provides insight into the complex interaction of host and pathogen.

Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse.

Single-nucleotide polymorphisms (SNPs) have been the focus of much attention in human genetics because they are extremely abundant and well-suited for automated large-scale genotyping. Human SNPs, however, are less informative than other types of genetic markers (such as simple-sequence length polymorphisms or microsatellites) and thus more loci are required for mapping traits. SNPs offer similar advantages for experimental genetic organisms such as the mouse, but they entail no loss of informativeness because bi-allelic markers are fully informative in analysing crosses between inbred strains. Here we report a large-scale analysis of SNPs in the mouse genome. We characterized the rate of nucleotide polymorphism in eight mouse strains and identified a collection of 2,848 SNPs located in 1,755 sequence-tagged sites (STSs) using high-density oligonucleotide arrays. Three-quarters of these SNPs have been mapped on the mouse genome, providing a first-generation SNP map of the mouse. We have also developed a multiplex genotyping procedure by which a genome scan can be performed with only six genotyping reactions per animal.

Emergence of dual resistance to zidovudine and lamivudine in HIV-1-infected patients treated with zidovudine plus lamivudine as initial therapy.

Presence of mutations associated with resistance to zidovudine or lamivudine was determined in isolates of HIV-1 obtained after long-term follow-up of 64 infected individuals who received zidovudine, lamivudine, or both drugs as initial antiretroviral therapy. Zidovudine resistance mutations were less frequent in isolates from patients treated with combination lamivudine plus zidovudine compared with zidovudine alone, but these mutations accumulated over time. Phenotypic resistance to both drugs was found in isolates from 3 of 23 patients. In 3 other patients, lamivudine-resistant virus detected at week 12 was replaced by wild-type virus after longer follow-up, which correlated with a return to baseline levels of plasma HIV-1 RNA. These results show that dual resistance to zidovudine and lamivudine develops over time despite the delayed emergence of zidovudine-resistant mutations. These results also suggest a selective advantage in vivo for HIV-1 species that are wild-type at RT codon 184.

HIV-1 infectability of CD4+ lymphocytes with relation to beta-chemokines and the CCR5 coreceptor.

CD4+ lymphocytes exhibit variable permissiveness to the replication of HIV-1. A cohort of sexually-exposed-yet-uninfected individuals were previously shown to have CD4+ lymphocytes refractory for M-tropic viral replication. In particular, two individuals from this population whose CD4+ lymphocytes exhibited complete resistance to M-tropic viral replication were later shown to be homozygous for a 32bp (delta32) deletion in the gene encoding for CCR5. In screening diverse populations, HIV-1 infected individuals heterozygous for the delta32 allele were statistically favored in their disease course to harbor lower viral loads and exhibit slower rates of CD4+ cell loss when compared to control CCR5 wild-type individuals. Further comparative analysis between individuals in the exposed but uninfected cohort who demonstrated intermediate levels of in vitro viral replication and CD4+ lymphocytes isolated from uninfected delta32 heterozygous individuals indicate that reduced levels of in vitro M-tropic replication are a CCR5-related phenomenon: CD4+ lymphocytes from these individuals were more sensitive to the HIV-1 blocking effects of recombinant chemokines, displayed lower CCR5 cell surface expression levels and a proportionate increase in the production of RANTES when compared to CD4+ lymphocytes from control individuals. These results suggest that the CCR5 phenotype is important in determining the replicative capacity of M-tropic HIV-1 in vitro. The implications of these results with relation to HIV-1 transmission and disease progression are discussed.

High density synthetic oligonucleotide arrays.

Experimental genomics involves taking advantage of sequence information to investigate and understand the workings of genes, cells and organisms. We have developed an approach in which sequence information is used directly to design high-density, two-dimensional rays of synthetic oligonucleotides. The GeneChipe probe arrays are made using spatially patterned, light-directed combinatorial chemical synthesis and contain up to hundreds of thousands of different oligonucleotides on a small glass surface. The arrays have been designed and used for quantitative and highly parallel measurements of gene expression, to discover polymorphic loci and to detect the presence of thousands of alternative alleles. Here, we describe the fabrication of the arrays, their design and some specific applications to high-throughput genetic and cellular analysis.

Mycobacterium species identification and rifampin resistance testing with high-density DNA probe arrays.

Species identification within the genus Mycobacterium and subsequent antibiotic susceptibility testing still rely on time-consuming, culture-based methods. Despite the recent development of DNA probes, which greatly reduce assay time, there is a need for a single platform assay capable of answering the multitude of diagnostic questions associated with this genus. We describe the use of a DNA probe array based on two sequence databases: one for the species identification of mycobacteria (82 unique 16S rRNA sequences corresponding to 54 phenotypical species) and the other for detecting Mycobacterium tuberculosis rifampin resistance (rpoB alleles). Species identification or rifampin resistance was determined by hybridizing fluorescently labeled, amplified genetic material generated from bacterial colonies to the array. Seventy mycobacterial isolates from 27 different species and 15 rifampin-resistant M. tuberculosis strains were tested. A total of 26 of 27 species were correctly identified as well as all of the rpoB mutants. This parallel testing format opens new perspectives in terms of patient management for bacterial diseases by allowing a number of genetic tests to be simultaneously run.

Cellular gene expression altered by human cytomegalovirus: global monitoring with oligonucleotide arrays.

Mechanistic insights to viral replication and pathogenesis generally have come from the analysis of viral gene products, either by studying their biochemical activities and interactions individually or by creating mutant viruses and analyzing their phenotype. Now it is possible to identify and catalog the host cell genes whose mRNA levels change in response to a pathogen. We have used DNA array technology to monitor the level of approximately 6,600 human mRNAs in uninfected as compared with human cytomegalovirus-infected cells. The level of 258 mRNAs changed by a factor of 4 or more before the onset of viral DNA replication. Several of these mRNAs encode gene products that might play key roles in virus-induced pathogenesis, identifying them as intriguing targets for further study.

Simultaneous genotyping and species identification using hybridization pattern recognition analysis of generic Mycobacterium DNA arrays.

High-density oligonucleotide arrays can be used to rapidly examine large amounts of DNA sequence in a high throughput manner. An array designed to determine the specific nucleotide sequence of 705 bp of the rpoB gene of Mycobacterium tuberculosis accurately detected rifampin resistance associated with mutations of 44 clinical isolates of M. tuberculosis. The nucleotide sequence diversity in 121 Mycobacterial isolates (comprised of 10 species) was examined by both conventional dideoxynucleotide sequencing of the rpoB and 16S genes and by analysis of the rpoB oligonucleotide array hybridization patterns. Species identification for each of the isolates was similar irrespective of whether 16S sequence, rpoB sequence, or the pattern of rpoB hybridization was used. However, for several species, the number of alleles in the 16S and rpoB gene sequences provided discordant estimates of the genetic diversity within a species. In addition to confirming the array's intended utility for sequencing the region of M. tuberculosis that confers rifampin resistance, this work demonstrates that this array can identify the species of nontuberculous Mycobacteria. This demonstrates the general point that DNA microarrays that sequence important genomic regions (such as drug resistance or pathogenicity islands) can simultaneously identify species and provide some insight into the organism's population structure.

Reduced HIV-1 infectability of CD4+ lymphocytes from exposed-uninfected individuals: association with low expression of CCR5 and high production of beta-chemokines.

We examined the human immunodeficiency virus type 1 infectability of CD4+ lymphocytes isolated from CCR5 wild-type individuals, individuals heterozygous for the delta32 allele of CCR5, and HIV-1-exposed but uninfected (EU) individuals who had CD4+ lymphocytes refractory to M-tropic viral replication. None of the EU individuals were found to be heterozygous for the delta32 allele. The CD4+ lymphocytes isolated from CCR5/delta32 and EU individuals were less infectable with an M-tropic viral isolate of HIV-1 than CCR5/CCR5 control individuals but were equally as infectable with a T-tropic viral isolate. The restriction to M-tropic viral isolate replication did not associate with any profound genotypic change in the CCR5 gene. CD4+ lymphocytes from CCR5/delta32 and CCR5/CCR5 EU individuals were more sensitive to the HIV-inhibitory effects of the recombinant beta-chemokines RANTES, MIP-1alpha, and MIP-1beta than were CD4+ lymphocytes from CCR5/CCR5 control individuals. CD4+ lymphocytes from EU individuals also showed increased sensitivity to recombinant beta-chemokines and low surface expression of CCR5. A phenotype of low CCR5 expression and high secretion of beta-chemokines is associated with reduced infectability of cells by M-tropic HIV-1. This phenotype may also be associated with protection against sexual transmission of HIV-1.

Human immunodeficiency virus type 1 reverse transcriptase genotype and drug susceptibility changes in infected individuals receiving dideoxyinosine monotherapy for 1 to 2 years.

The genetic mechanisms of human immunodeficiency virus type 1 (HIV-1) resistance to dideoxyinosine (ddI) in vivo have been described based on data from primary HIV-1 isolates. To better define the spectrum of HIV-1 reverse transcriptase (RT) changes occurring during ddI therapy, we determined the genotype and ddI susceptibility of the RT gene of HIV RNA isolated from the plasma of 23 patients who had received 1 to 2 years (mean, 87 +/- 16 weeks) of ddI monotherapy. Population-based sequencing of plasma virus showed that 12 of 23 (52%) patients developed known ddI resistance mutations: L74V (7 patients), K65R (2 patients), L74V with M184V (3 patients), and L74V with K65R (1 patient). Five patients developed one or more known zidovudine resistance mutations (at codons 41, 67, 70, 215, and/or 219) during the study. Other amino acid substitutions were found, but only S68G and L210W occurred in more than one patient. Studies of sensitivity to ddI were performed on population-based recombinant-virus stocks generated by homologous recombination between a plasmid containing an HXB2 clone with the RT gene deleted and RT-PCR products of the RT genes from patients' plasma RNA. The sequences of the virus stocks produced by this procedure were typically identical to the sequence of the input PCR product from plasma RNA. Both an MT-2 cell-based culture assay and a cell-free virion-associated RT inhibition assay showed that viruses possessing an L74V and/or M184V mutation or a K65R mutation had reduced sensitivity to ddI. Viruses without these specific mutations had no change in sensitivity to ddI. The results presented here show that the spectrum of RT mutations in a population of patients on ddI monotherapy is more complex than previously described. The development of multiple mutational patterns, including those that confer resistance to other nucleoside analogs, highlights the complexity of using the currently available nucleoside analogs for antiretroviral therapy.

Extensive polymorphisms observed in HIV-1 clade B protease gene using high-density oligonucleotide arrays.

Naturally occurring mutations in HIV-1-infected patients have important implications for therapy and the outcome of clinical studies. However, little is known about the prevalence of mutations that confer resistance to HIV-1 protease inhibitors in isolates derived from patients naive for such inhibitors. In the first clinical application of high-density oligonucleotide array sequencing, the sequences of 167 viral isolates from 102 patients have been determined. The DNA sequence of USA HIV-1 clade B proteases was found to be extremely variable and 47.5% of the 99 amino acid positions varied. This level of amino acid diversity is greater than that previously known for all worldwide HIV-1 clades combined (40%). Many of the amino acid changes that are known to contribute to drug resistance occurred as natural polymorphisms in isolates from patients who had never received protease inhibitors.