Evolutionary Divergence of Aggregatibacter actinomycetemcomitans.

Gram-negative facultative Aggregatibacter actinomycetemcomitans is an oral pathogen associated with periodontitis. The genetic heterogeneity among A. actinomycetemcomitans strains has been long recognized. This study provides a comprehensive genomic analysis of A. actinomycetemcomitans and the closely related nonpathogenic Aggregatibacter aphrophilus. Whole genome sequencing by Illumina MiSeq platform was performed for 31 A. actinomycetemcomitans and 2 A. aphrophilus strains. Sequence similarity analysis shows a total of 3,220 unique genes across the 2 species, where 1,550 are core genes present in all genomes and 1,670 are variable genes (accessory genes) missing in at least 1 genome. Phylogenetic analysis based on 397 concatenated core genes distinguished A. aphrophilus and A. actinomycetemcomitans. The latter was in turn divided into 5 clades: clade b (serotype b), clade c (serotype c), clade e/f (serotypes e and f), clade a/d (serotypes a and d), and clade e' (serotype e strains). Accessory genes accounted for 14.1% to 23.2% of the A. actinomycetemcomitans genomes, with a majority belonging to the category of poorly characterized by Cluster of Orthologous Groups classification. These accessory genes were often organized into genomic islands (n = 387) with base composition biases, suggesting their acquisitions via horizontal gene transfer. There was a greater degree of similarity in gene content and genomic islands among strains within clades than between clades. Strains of clade e' isolated from human were found to be missing the genomic island that carries genes encoding cytolethal distending toxins. Taken together, the results suggest a pattern of sequential divergence, starting from the separation of A. aphrophilus and A. actinomycetemcomitans through gain and loss of genes and ending with the divergence of the latter species into distinct clades and serotypes. With differing constellations of genes, the A. actinomycetemcomitans clades may have evolved distinct adaptation strategies to the human oral cavity.

Comparative genomic hybridization and transcriptome analysis with a pan-genome microarray reveal distinctions between JP2 and non-JP2 genotypes of Aggregatibacter actinomycetemcomitans.

It was postulated that the highly virulent JP2 genotype of Aggregatibacter actinomycetemcomitans may possess a constellation of distinct virulence determinants not found in non-JP2 genotypes. This study compared the genome content and the transcriptome of the serotype b JP2 genotype and the closely related serotype b non-JP2 genotype of A. actinomycetemcomitans. A custom-designed pan-genomic microarray of A. actinomycetemcomitans was constructed and validated against a panel of 11 sequenced reference strains. The microarray was subsequently used for comparative genomic hybridization of serotype b strains of JP2 (six strains) and non-JP2 (six strains) genotypes, and for transcriptome analysis of strains of JP2 (three strains) and non-JP2 (two strains). Two JP2-specific and two non-JP2-specific genomic islands were identified. In one instance, distinct genomic islands were found to be inserted into the same locus among strains of different genotypes. Transcriptome analysis identified five operons, including the leukotoxin operon, to have at least two genes with an expression ratio of 2 or greater between genotypes. Two of the differentially expressed operons were members of the membrane-bound nitrate reductase system (nap operon) and the Tol-Pal system of gram-negative bacterial species. This study is the first to demonstrate the differences in the full genome content and gene expression between A. actinomycetemcomitans strains of JP2 and non-JP2 genotypes. The information is essential for designing hypothesis-driven experiments to examine the pathogenic mechanisms of A. actinomycetemcomitans.

Comparison of Lumin LED fluorescent attachment, fluorescent microscopy and Ziehl-Neelsen for AFB diagnosis.

DESIGN: Cross-sectional studies in Russia (n = 502) and Macedonia (n = 205), with fluorochrome-stained sputum examined by 1) the new Lumin light emitting diode (LED) fluorescent attachment on a light microscope, and 2) conventional fluorescent microscope (CFM) available in each laboratory, and compared to 3) Ziehl-Neelsen (ZN) restaining/reading of the same smears. Poor readings of ZN-restained smears in Russia stimulated a retrospective laboratory registry analysis for sensitivity and specificity of directly ZN-stained smears (n = 791) from a previous period. RESULTS: In Macedonia, the sensitivity of the Lumin and CFM were 87.8%, and that of restained ZN smears with conventional light microscope was 78.0%. In Russia, sensitivity was as follows: Lumin 72.8%, CFM 52.5%; re-stained ZN smears 28.5% and directly ZN stained smears 55.6%. CONCLUSION: Fluorescence microscopy is more sensitive than conventional microscopy. The Lumin attachment to conventional light microscopes provided results equal to or better than the CFMs. Smear restaining for ZN showed a 12% advantage for Lumin and CFM in Macedonia, in line with other meta-analyses. Restaining for ZN gave poor results in Russia for unknown reasons. Retrospective analysis of directly ZN-stained smears showed 55.6% sensitivity compared to the Lumin (72.8%), which is also in line with the superiority of fluorescent microscopy reported in literature.

Context-specific infinite mixtures for clustering gene expression profiles across diverse microarray dataset.

MOTIVATION: Identifying groups of co-regulated genes by monitoring their expression over various experimental conditions is complicated by the fact that such co-regulation is condition-specific. Ignoring the context-specific nature of co-regulation significantly reduces the ability of clustering procedures to detect co-expressed genes due to additional 'noise' introduced by non-informative measurements. RESULTS: We have developed a novel Bayesian hierarchical model and corresponding computational algorithms for clustering gene expression profiles across diverse experimental conditions and studies that accounts for context-specificity of gene expression patterns. The model is based on the Bayesian infinite mixtures framework and does not require a priori specification of the number of clusters. We demonstrate that explicit modeling of context-specificity results in increased accuracy of the cluster analysis by examining the specificity and sensitivity of clusters in microarray data. We also demonstrate that probabilities of co-expression derived from the posterior distribution of clusterings are valid estimates of statistical significance of created clusters. AVAILABILITY: The open-source package gimm is available at http://eh3.uc.edu/gimm.

Bayesian mixture model based clustering of replicated microarray data.

MOTIVATION: Identifying patterns of co-expression in microarray data by cluster analysis has been a productive approach to uncovering molecular mechanisms underlying biological processes under investigation. Using experimental replicates can generally improve the precision of the cluster analysis by reducing the experimental variability of measurements. In such situations, Bayesian mixtures allow for an efficient use of information by precisely modeling between-replicates variability. RESULTS: We developed different variants of Bayesian mixture based clustering procedures for clustering gene expression data with experimental replicates. In this approach, the statistical distribution of microarray data is described by a Bayesian mixture model. Clusters of co-expressed genes are created from the posterior distribution of clusterings, which is estimated by a Gibbs sampler. We define infinite and finite Bayesian mixture models with different between-replicates variance structures and investigate their utility by analyzing synthetic and the real-world datasets. Results of our analyses demonstrate that (1) improvements in precision achieved by performing only two experimental replicates can be dramatic when the between-replicates variability is high, (2) precise modeling of intra-gene variability is important for accurate identification of co-expressed genes and (3) the infinite mixture model with the 'elliptical' between-replicates variance structure performed overall better than any other method tested. We also introduce a heuristic modification to the Gibbs sampler based on the 'reverse annealing' principle. This modification effectively overcomes the tendency of the Gibbs sampler to converge to different modes of the posterior distribution when started from different initial positions. Finally, we demonstrate that the Bayesian infinite mixture model with 'elliptical' variance structure is capable of identifying the underlying structure of the data without knowing the 'correct' number of clusters. AVAILABILITY: The MS Windows based program named Gaussian Infinite Mixture Modeling (GIMM) implementing the Gibbs sampler and corresponding C++ code are available at http://homepages.uc.edu/~medvedm/GIMM.htm SUPPLEMENTAL INFORMATION: http://expression.microslu.washington.edu/expression/kayee/medvedovic2003/medvedovic_bioinf2003.html

Gene expression in Pseudomonas aeruginosa biofilms.

Bacteria often adopt a sessile biofilm lifestyle that is resistant to antimicrobial treatment. Opportunistic pathogenic bacteria like Pseudomonas aeruginosa can develop persistent infections. To gain insights into the differences between free-living P. aeruginosa cells and those in biofilms, and into the mechanisms underlying the resistance of biofilms to antibiotics, we used DNA microarrays. Here we show that, despite the striking differences in lifestyles, only about 1% of genes showed differential expression in the two growth modes; about 0.5% of genes were activated and about 0.5% were repressed in biofilms. Some of the regulated genes are known to affect antibiotic sensitivity of free-living P. aeruginosa. Exposure of biofilms to high levels of the antibiotic tobramycin caused differential expression of 20 genes. We propose that this response is critical for the development of biofilm resistance to tobramycin. Our results show that gene expression in biofilm cells is similar to that in free-living cells but there are a small number of significant differences. Our identification of biofilm-regulated genes points to mechanisms of biofilm resistance to antibiotics.

A comprehensive view of regulation of gene expression by double-stranded RNA-mediated cell signaling.

Double-stranded (ds) RNA, a common component of virus-infected cells, is a potent inducer of the type I interferon and other cellular genes. For identifying the full repertoire of human dsRNA-regulated genes, a cDNA microarray hybridization screening was conducted using mRNA from dsRNA-treated GRE cells. Because these cells lack all type I interferon genes, the possibility of gene induction by autocrine actions of interferon was eliminated. Our screen identified 175 dsRNA-stimulated genes (DSG) and 95 dsRNA-repressed genes. A subset of the DSGs was also induced by different inflammatory cytokines and viruses demonstrating interconnections among disparate signaling pathways. Functionally, the DSGs encode proteins involved in signaling, apoptosis, RNA synthesis, protein synthesis and processing, cell metabolism, transport, and structure. Induction of such a diverse family of genes by dsRNA has major implications in host-virus interactions and in the use of RNA(i) technology for functional ablation of specific genes.

Integrated genomic and proteomic analyses of a systematically perturbed metabolic network.

We demonstrate an integrated approach to build, test, and refine a model of a cellular pathway, in which perturbations to critical pathway components are analyzed using DNA microarrays, quantitative proteomics, and databases of known physical interactions. Using this approach, we identify 997 messenger RNAs responding to 20 systematic perturbations of the yeast galactose-utilization pathway, provide evidence that approximately 15 of 289 detected proteins are regulated posttranscriptionally, and identify explicit physical interactions governing the cellular response to each perturbation. We refine the model through further iterations of perturbation and global measurements, suggesting hypotheses about the regulation of galactose utilization and physical interactions between this and a variety of other metabolic pathways.

Global impact of influenza virus on cellular pathways is mediated by both replication-dependent and -independent events.

Influenza virus, the causative agent of the common flu, is a worldwide health problem with significant economic consequences. Studies of influenza virus biology have revealed elaborate mechanisms by which the virus interacts with its host cell as it inhibits the synthesis of cellular proteins, evades the innate antiviral response, and facilitates production of viral RNAs and proteins. With the advent of DNA array technology it is now possible to obtain a large-scale view of how viruses alter the environment within the host cell. In this study, the cellular response to influenza virus infection was examined by monitoring the steady-state mRNA levels for over 4,600 cellular genes. Infections with active and inactivated influenza viruses identified changes in cellular gene expression that were dependent on or independent of viral replication, respectively. Viral replication resulted in the downregulation of many cellular mRNAs, and the effect was enhanced with time postinfection. Interestingly, several genes involved in protein synthesis, transcriptional regulation, and cytokine signaling were induced by influenza virus replication, suggesting that some may play essential or accessory roles in the viral life cycle or the host cell's stress response. The gene expression pattern induced by inactivated viruses revealed induction of the cellular metallothionein genes that may represent a protective response to virus-induced oxidative stress. Genome-scale analyses of virus infections will help us to understand the complexities of virus-host interactions and may lead to the discovery of novel drug targets or antiviral therapies.

Prostate short-chain dehydrogenase reductase 1 (PSDR1): a new member of the short-chain steroid dehydrogenase/reductase family highly expressed in normal and neoplastic prostate epithelium.

Genes regulated by androgenic hormones are of critical importance for the normal physiological function of the human prostate gland, and they contribute to the development and progression of prostate carcinoma. We used cDNA microarrays comprised of prostate-derived cDNAs to profile transcripts regulated by androgens in prostate cancer cells. This study identified a novel gene that we have designated prostate short-chain dehydrogenase/reductase 1 (PSDR1), that exhibits increased expression on exposure to androgens in the LNCaP prostate cancer cell line. Northern analysis demonstrated that PSDR1 is highly expressed in the prostate gland relative to other normal human tissues. The PSDR1 cDNA and putative protein exhibit homology to the family of short-chain dehydrogenase/reductase enzymes and thus identify a new member of this family. Cloning and analysis of the putative PSDR1 promoter region identified a potential androgen-response element. We used a radiation-hybrid panel to map the PSDR1 gene to chromosome 14q23-24.3. In situ hybridization localizes PSDR1 expression to normal and neoplastic prostate epithelium. These results identify a new gene involved in the androgen receptor-regulated gene network of the human prostate that may play a role in the pathogenesis of prostate carcinoma.

Interaction of pseudomonas aeruginosa with epithelial cells: identification of differentially regulated genes by expression microarray analysis of human cDNAs.

Pseudomonas aeruginosa is an opportunistic pathogen that plays a major role in lung function deterioration in cystic fibrosis patients. To identify critical host responses during infection, we have used high-density DNA microarrays, consisting of 1,506 human cDNA clones, to monitor gene expression in the A549 lung pneumocyte cell line during exposure to P. aeruginosa. We have identified host genes that are differentially expressed upon infection, several of which require interaction with P. aeruginosa and the expression of the major subunit of type IV pili, PilA. Differential expression of genes involved in various cellular functions was identified, and we selected the gene encoding the transcription factor interferon regulatory factor 1 (IRF-1) for further analysis. The levels of the IRF-1 transcript increased 3- to 4-fold in A549 cells after adherence by P. aeruginosa. A similar increase of IRF-1 mRNA was observed in A549 cells exposed to wild-type P. aeruginosa when compared to an isogenic, nonpiliated strain. However, this difference was abolished when serum was present during the incubation of bacteria. Exposure of A549 cells to purified P. aeruginosa lipopolysaccharide did not result in a significant increase in IRF-1 mRNA. Although the P. aeruginosa-induced increased IRF-1 expression depends on the presence of bacterial adhesin, our findings do not preclude the possibility that other bacterial products are responsible for IRF-1 activation, which is enhanced by bacterial adherence to cells. These data show that microarray technology can be an important tool for studying the complex interplay between bacterial pathogens and host.

PART-1: a novel human prostate-specific, androgen-regulated gene that maps to chromosome 5q12.

Genes regulated by androgenic hormones are of critical importance for the normal physiological function of the human prostate gland, and they contribute to the development and progression of prostate carcinoma. We used cDNA microarrays containing 1500 prostate-derived cDNAs to profile transcripts regulated by androgens in prostate cancer cells. This study identified a novel gene that we have designated PART-1 (prostate androgen-regulated transcript 1), which exhibited increased expression upon exposure to androgens in the LNCaP prostate cancer cell line. Northern analysis demonstrated that PART-1 is highly expressed in the prostate gland relative to other normal human tissues and is expressed as different transcripts using at least three different polyadenylation signals. The PART-1 cDNA and putative protein are not significantly homologous to any sequences in the nonredundant public sequence databases. Cloning and analysis of the putative PART-1 promoter region identified a potential binding site for the homeobox gene PBX-la, but no consensus androgen response element or sterol-regulatory element binding sites were identified. We used a radiation hybrid panel and fluorescence in situ hybridization to map the PART-1 gene to chromosome 5q12, a region that has been suggested to harbor a prostate tumor suppressor gene. These results identify a new gene involved in the androgen receptor-regulated gene network of the human prostate that may play a role in the etiology of prostate carcinogenesis.

Large-scale monitoring of host cell gene expression during HIV-1 infection using cDNA microarrays.

Human immunodeficiency virus type 1 (HIV-1) infection alters the expression of host cell genes at both the mRNA and protein levels. To obtain a more comprehensive view of the global effects of HIV infection of CD4-positive T-cells at the mRNA level, we performed cDNA microarray analysis on approximately 1500 cellular cDNAs at 2 and 3 days postinfection (p.i.) with HIV-1. Host cell gene expression changed little at 2 days p.i., but at 3 days p.i. 20 cellular genes were identified as differentially expressed. Genes involved in T-cell signaling, subcellular trafficking, and transcriptional regulation, as well as several uncharacterized genes, were among those whose mRNAs were differentially regulated. These results support the hypothesis that HIV-1 infection alters expression of a broad array of cellular genes and provides a framework for future functional studies on the differentially expressed mRNA products.

The prostate expression database (PEDB): status and enhancements in 2000.

The Prostate Expression Database (PEDB) is an online resource designed to access and analyze gene expression information derived from the human prostate. PEDB archives >55 000 expressed sequence tags (ESTs) from 43 cDNA libraries in a curated relational database that provides detailed library information including tissue source, library construction methods, sequence diversity and sequence abundance. The differential expression of each EST species can be viewed across all libraries using a Virtual Expression Analysis Tool (VEAT), a graphical user interface written in Java for intra- and inter-library species comparisons. Recent enhancements to PEDB include: (i) the functional categorization of annotated EST assemblies using a classification scheme developed at The Institute for Genome Research; (ii) catalogs of expressed genes in specific prostate tissue sources designated as transcriptomes; and (iii) the addition of prostate proteome information derived from two-dimensional electrophoreses and mass spectrometry of prostate cancer cell lines. PEDB may be accessed via the WWW at http://www.mbt.washington.edu/PEDB/

Negative selection: a method for obtaining low-abundance cDNAs using high-density cDNA clone arrays.

The identification of the entire complement of genes expressed in a cell, tissue, or organism provides a framework for understanding biological properties and establishes a tool set for subsequent functional studies. The large-scale sequencing of randomly selected clones from cDNA libraries has been successfully employed as a method for identifying a large fraction of these expressed genes. However, this approach is limited by the inherent redundancy of cellular transcripts reflecting widely variant levels of gene transcription. As a result, a high percentage of transcript duplications are encountered as the number of sequenced clones accrues. To address this problem, we have developed a negative hybridization selection method that employs the hybridization of complex cDNA probes to high-density arrays of cDNA clones and the subsequent selection of clones with a null or low hybridization signal. This approach was applied to a cDNA library constructed from normal human prostate tissue and resulted in the reduction of highly expressed prostate cDNAs from 6.8 to 0.57% with an overall decline in clone redundancy from 33 to 11%. The selected clones also reflected a more diverse cDNA population, with 89% of the clones representing distinctly different cDNAs compared with 67% of the randomly selected clones. This method compares favorably with cDNA library re-association normalization approaches and offers several distinct advantages, including the flexibility to use previously prepared libraries, and the ability to employ an iterative screening approach for continued accrual of cDNAs representing rare transcripts.

Comparative hybridization of an array of 21,500 ovarian cDNAs for the discovery of genes overexpressed in ovarian carcinomas.

Comparative hybridization of cDNA arrays is a powerful tool for the measurement of differences in gene expression between two or more tissues. We optimized this technique and employed it to discover genes with potential for the diagnosis of ovarian cancer. This cancer is rarely identified in time for a good prognosis after diagnosis. An array of 21,500 unknown ovarian cDNAs was hybridized with labeled first-strand cDNA from 10 ovarian tumors and six normal tissues. One hundred and thirty-four clones are overexpressed in at least five of the 10 tumors. These cDNAs were sequenced and compared to public sequence databases. One of these, the gene HE4, was found to be expressed primarily in some ovarian cancers, and is thus a potential marker of ovarian carcinoma.

PEDB: the Prostate Expression Database.

The Prostate Expression Database (PEDB) is a curated relational database and suite of analysis tools designed for the study of prostate gene expression in normal and disease states. Expressed Sequence Tags (ESTs) and full-length cDNA sequences derived from more than 40 human prostate cDNA libraries are maintained and represent a wide spectrum of normal and pathological conditions. Detailed library information including tissue source, library construction methods, sequence diversity and abundance are available in a library archive. Prostate ESTs are assembled into distinct species groups using the multiple alignment program CAP2 and are annotated with information from the GenBank, dbEST and Unigene public sequence databases. Annotated sequences in PEDB are searched using the BLAST algorithm. The differential expression of each EST species can be viewed across all libraries using a Virtual Expression Analysis Tool (VEAT), a graphical user interface written in Java for intra- and inter-library species comparisons. PEDB may be accessed via the World Wide Web at http://www.mbt.washington.edu/PEDB/

Snapshot of a large dynamic replicon in a halophilic archaeon: megaplasmid or minichromosome?

Extremely halophilic archaea, which flourish in hypersaline environments, are known to contain a variety of large dynamic replicons. Previously, the analysis of one such replicon, pNRC100, in Halobacterium sp. strain NRC-1, showed that it undergoes high-frequency insertion sequence (IS) element-mediated insertions and deletions, as well as inversions via recombination between 39-kb-long inverted repeats (IRs). Now, the complete sequencing of pNRC100, a 191,346-bp circle, has shown the presence of 27 IS elements representing eight families. A total of 176 ORFs or likely genes of 850-bp average size were found, 39 of which were repeated within the large IRs. More than one-half of the ORFs are likely to represent novel genes that have no known homologs in the databases. Among ORFs with previously characterized homologs, three different copies of putative plasmid replication and four copies of partitioning genes were found, suggesting that pNRC100 evolved from IS element-mediated fusions of several smaller plasmids. Consistent with this idea, putative genes typically found on plasmids, including those encoding a restriction-modification system and arsenic resistance, as well as buoyant gas-filled vesicles and a two-component regulatory system, were found on pNRC100. However, additional putative genes not expected on an extrachromosomal element, such as those encoding an electron transport chain cytochrome d oxidase, DNA nucleotide synthesis enzymes thioredoxin and thioredoxin reductase, and eukaryotic-like TATA-binding protein transcription factors and a chromosomal replication initiator protein were also found. A multi-step IS element-mediated process is proposed to account for the acquisition of these chromosomal genes. The finding of essential genes on pNRC100 and its property of resistance to curing suggest that this replicon may be evolving into a new chromosome.

An expressed-sequence-tag database of the human prostate: sequence analysis of 1168 cDNA clones.

The human prostate is a complex glandular organ with functional development under hormonal regulation. Diseases of the prostate result in significant morbidity and mortality in the form of benign prostatic hypertrophy and prostate adenocarcinoma. The characterization of the molecular framework of the human prostate at the level of expressed genes will facilitate the understanding of normal and pathological prostate biology. The purposes of this study were to acquire an initial assessment of the qualitative and quantitative diversity of gene expression in the normal human prostate and to determine the extent that genes with prostate-restricted expression can be assessed using an expressed sequence tag approach. We have constructed a directional cDNA library from normal adult human prostate tissue and partially sequenced the 5' end of 1168 randomly selected cDNA clones, resulting in more than 400 kb of DNA sequence. Homology searches of the sequenced cDNAs against the GenBank and dbEST databases revealed that 43% of the sequences are identical to human genes whose functions are known, 5% are similar but not identical to known genes in humans or lower organisms, 5% match the mitochondrial genome, 9% are composed of interspersed DNA repeats, 30% are homologous to sequences in the dbEST database without a described function, and 6% are novel sequences. A total of 780 distinct species were identified. In addition to the 74 novel transcripts, 4 genes, prostate-specific antigen (PSA), prostate secretory protein (PSP), prostate acid phosphatase (PAP), and human glandular kallekrein 2 (HK2), have no homologous sequences in the databases that originate from sources other than prostate and thus may represent genes with prostate-restricted expression. Sequences matching PSA, PSP, and PAP each accounted for > 1% of the total ESTs and represent highly abundant transcripts, correlating with the abundance of these proteins in the prostate gland. No novel transcripts were represented by more than one EST and thus are expressed at levels much lower than the known prostate-specific genes.