Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

AF5q31, a newly identified AF4-related gene, is fused to MLL in infant acute lymphoblastic leukemia with ins(5;11)(q31;q13q23)

Infant acute lymphoblastic leukemia (ALL) with MLL gene rearrangements is characterized by early pre-B phenotype (CD10−/CD19+) and poor treatment outcome. The t(4;11), creating MLL-AF4 chimeric transcripts, is the predominant 11q23 chromosome translocation in infant ALL and is associated with extremely poor prognosis as compared with other 11q23 translocations. We analyzed an infant early preB ALL with ins(5;11)(q31;q13q23) and identified the AF5q31 gene on chromosome 5q31 as a fusion partner of the MLL gene. The AF5q31 gene, which encoded a protein of 1,163 aa, was located in the vicinity of the cytokine cluster region of chromosome 5q31 and contained at least 16 exons. The AF5q31 gene was expressed in fetal heart, lung, and brain at relatively high levels and fetal liver at a low level, but the expression in these tissues decreased in adults. The AF5q31 protein was homologous to AF4-related proteins, including AF4, LAF4, and FMR2. The AF5q31 and AF4 proteins had three homologous regions, including the transactivation domain of AF4, and the breakpoint of AF5q31 was located within the region homologous to the transactivation domain of AF4. Furthermore, the clinical features of this patient with the MLL-AF5q31 fusion transcript, characterized by the early pre-B phenotype (CD10−/CD19+) and poor outcome, were similar to those of patients having MLL-AF4 chimeric transcripts. These findings suggest that AF5q31 and AF4 might define a new family particularly involved in the pathogenesis of 11q23-associated-ALL.

Cluster analysis and display of genome-wide expression patterns

A system of cluster analysis for genome-wide expression data from DNA microarray hybridization is described that uses standard statistical algorithms to arrange genes according to similarity in pattern of gene expression. The output is displayed graphically, conveying the clustering and the underlying expression data simultaneously in a form intuitive for biologists. We have found in the budding yeast Saccharomyces cerevisiae that clustering gene expression data groups together efficiently genes of known similar function, and we find a similar tendency in human data. Thus patterns seen in genome-wide expression experiments can be interpreted as indications of the status of cellular processes. Also, coexpression of genes of known function with poorly characterized or novel genes may provide a simple means of gaining leads to the functions of many genes for which information is not available currently.

Hox cluster genomics in the horn shark, Heterodontus francisci

Reconstructing the evolutionary history of Hox cluster origins will lead to insights into the developmental and evolutionary significance of Hox gene clusters in vertebrate phylogeny and to their role in the origins of various vertebrate body plans. We have isolated two Hox clusters from the horn shark, Heterodontus francisci. These have been sequenced and compared with one another and with other chordate Hox clusters. The results show that one of the horn shark clusters (HoxM) is orthologous to the mammalian HoxA cluster and shows a structural similarity to the amphioxus cluster, whereas the other shark cluster (HoxN) is orthologous to the mammalian HoxD cluster based on cluster organization and a comparison with noncoding and Hox gene-coding sequences. The persistence of an identifiable HoxA cluster over an 800-million-year divergence time demonstrates that the Hox gene clusters are highly integrated and structured genetic entities. The data presented herein identify many noncoding sequence motifs conserved over 800 million years that may function as genetic control motifs essential to the developmental process.

Global analysis of gene expression in pulmonary fibrosis reveals distinct programs regulating lung inflammation and fibrosis

The molecular mechanisms of pulmonary fibrosis are poorly understood. We have used oligonucleotide arrays to analyze the gene expression programs that underlie pulmonary fibrosis in response to bleomycin, a drug that causes lung inflammation and fibrosis, in two strains of susceptible mice (129 and C57BL/6). We then compared the gene expression patterns in these mice with 129 mice carrying a null mutation in the epithelial-restricted integrin β6 subunit (β6−/−), which develop inflammation but are protected from pulmonary fibrosis. Cluster analysis identified two distinct groups of genes involved in the inflammatory and fibrotic responses. Analysis of gene expression at multiple time points after bleomycin administration revealed sequential induction of subsets of genes that characterize each response. The availability of this comprehensive data set should accelerate the development of more effective strategies for intervention at the various stages in the development of fibrotic diseases of the lungs and other organs.

Novel role of phosphorylation in Fe–S cluster stability revealed by phosphomimetic mutations at Ser-138 of iron regulatory protein 1

Animals regulate iron metabolism largely through the action of the iron regulatory proteins (IRPs). IRPs modulate mRNA utilization by binding to iron-responsive elements (IRE) in the 5′ or 3′ untranslated region of mRNAs encoding proteins involved in iron homeostasis or energy production. IRP1 is also the cytosolic isoform of aconitase. The activities of IRP1 are mutually exclusive and are modulated through the assembly/disassembly of its [4Fe–4S] cluster, reversibly converting it between an IRE-binding protein and cytosolic aconitase. IRP1 is also phosphoregulated by protein kinase C, but the mechanism by which phosphorylation posttranslationally increases IRE binding activity has not been fully defined. To investigate this, Ser-138 (S138), a PKC phosphorylation site, was mutated to phosphomimetic glutamate (S138E), aspartate (S138D), or nonphosphorylatable alanine (S138A). The S138E IRP1 mutant and, to a lesser extent, the S138D IRP1 mutant were impaired in aconitase function in yeast when grown aerobically but not when grown anaerobically. Purified wild-type and mutant IRP1s could be reconstituted to active aconitases anaerobically. However, when exposed to oxygen, the [4Fe–4S] cluster of the S138D and S138E mutants decayed 5-fold and 20-fold faster, respectively, than was observed for wild-type IRP1. Our findings suggest that stability of the Fe–S cluster of IRP1 can be regulated by phosphorylation and reveal a mechanism whereby the balance between the IRE binding and [4Fe–4S] forms of IRP1 can be modulated independently of cellular iron status. Furthermore, our results show that IRP1 can function as an oxygen-modulated posttranscriptional regulator of gene expression.

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ≤ 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 × 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

The Arabidopsis CBF Gene Family Is Composed of Three Genes Encoding AP2 Domain-Containing Proteins Whose Expression Is Regulated by Low Temperature but Not by Abscisic Acid or Dehydration1

We have identified two genes from Arabidopsis that show high similarity with CBF1, a gene encoding an AP2 domain-containing transcriptional activator that binds to the low-temperature-responsive element CCGAC and induces the expression of some cold-regulated genes, increasing plant freezing tolerance. These two genes, which we have named CBF2 and CBF3, also encode proteins containing AP2 DNA-binding motifs. Furthermore, like CBF1, CBF2 and CBF3 proteins also include putative nuclear-localization signals and potential acidic activation domains. The CBF2 and CBF3 genes are linked to CBF1, constituting a cluster on the bottom arm of chromosome IV. The high level of similarity among the three CBF genes, their tandem organization, and the fact that they have the same transcriptional orientation all suggest a common origin. CBF1, CBF2, and CBF3 show identical expression patterns, being induced very rapidly by low-temperature treatment. However, in contrast to most of the cold-induced plant genes characterized, they are not responsive to abscisic acid or dehydration. Taken together, all of these data suggest that CBF2 and CBF3 may function as transcriptional activators, controlling the level of low-temperature gene expression and promoting freezing tolerance through an abscisic acid-independent pathway.

Analysis of gene expression during myc oncogene-induced lymphomagenesis in the bursa of Fabricius

The transcriptional effects of deregulated myc gene overexpression are implicated in tumorigenesis in a spectrum of experimental and naturally occurring neoplasms. In follicles of the chicken bursa of Fabricius, myc induction of B-cell neoplasia requires a target cell population present during early bursal development and progresses through preneoplastic transformed follicles to metastatic lymphomas. We developed a chicken immune system cDNA microarray to analyze broad changes in gene expression that occur during normal embryonic B-cell development and during myc-induced neoplastic transformation in the bursa. The number of mRNAs showing at least 3-fold change was greater during myc-induced lymphomagenesis than during normal development, and hierarchical cluster analysis of expression patterns revealed that levels of several hundred mRNAs varied in concert with levels of myc overexpression. A set of 41 mRNAs were most consistently elevated in myc-overexpressing preneoplastic and neoplastic cells, most involved in processes thought to be subject to regulation by Myc. The mRNAs for another cluster of genes were overexpressed in neoplasia independent of myc expression level, including a small subset with the expression signature of embryonic bursal lymphocytes. Overexpression of myc, and some of the genes overexpressed with myc, may be important for generation of preneoplastic transformed follicles. However, expression profiles of late metastatic tumors showed a large variation in concert with myc expression levels, and some showed minimal myc overexpression. Therefore, high-level myc overexpression may be more important in the early induction of these lymphomas than in maintenance of late-stage metastases.

Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian leukocyte receptor cluster

An extensive, highly diversified multigene family of novel immune-type receptor (nitr) genes has been defined in Danio rerio (zebrafish). The genes are predicted to encode type I transmembrane glycoproteins consisting of extracellular variable (V) and V-like C2 (V/C2) domains, a transmembrane region and a cytoplasmic tail. All of the genes examined encode immunoreceptor tyrosine-based inhibition motifs in the cytoplasmic tail. Radiation hybrid panel mapping and analysis of a deletion mutant line (b240) indicate that a minimum of ≈40 nitr genes are contiguous in the genome and span ≈0.6 Mb near the top of zebrafish linkage group 7. One flanking region of the nitr gene complex shares conserved synteny with a region of mouse chromosome 7, which shares conserved synteny with human 19q13.3-q13.4 that encodes the leukocyte receptor cluster. Antibody-induced crosslinking of Nitrs that have been introduced into a human natural killer cell line inhibits the phosphorylation of mitogen-activated protein kinase that is triggered by natural killer-sensitive tumor target cells. Nitrs likely represent intermediates in the evolution of the leukocyte receptor cluster.

The chaperone/usher pathways of Pseudomonas aeruginosa: Identification of fimbrial gene clusters (cup) and their involvement in biofilm formation

Pseudomonas aeruginosa, an important opportunistic human pathogen, persists in certain tissues in the form of specialized bacterial communities, referred to as biofilm. The biofilm is formed through series of interactions between cells and adherence to surfaces, resulting in an organized structure. By screening a library of Tn5 insertions in a nonpiliated P. aeruginosa strain, we identified genes involved in early stages of biofilm formation. One class of mutations identified in this study mapped in a cluster of genes specifying the components of a chaperone/usher pathway that is involved in assembly of fimbrial subunits in other microorganisms. These genes, not previously described in P. aeruginosa, were named cupA1–A5. Additional chaperone/usher systems (CupB and CupC) have been also identified in the genome of P. aeruginosa PAO1; however, they do not appear to play a role in adhesion under the conditions where the CupA system is expressed and functions in surface adherence. The identification of these putative adhesins on the cell surface of P. aeruginosa suggests that this organism possess a wide range of factors that function in biofilm formation. These structures appear to be differentially regulated and may function at distinct stages of biofilm formation, or in specific environments colonized by this organism.

The highly recombinogenic bz locus lies in an unusually gene-rich region of the maize genome

The bronze (bz) locus exhibits the highest rate of recombination of any gene in higher plants. To investigate the possible basis of this high rate of recombination, we have analyzed the physical organization of the region around the bz locus. Two adjacent bacterial artificial chromosome clones, comprising a 240-kb contig centered around the Bz-McC allele, were isolated, and 60 kb of contiguous DNA spanning the two bacterial artificial chromosome clones was sequenced. We find that the bz locus lies in an unusually gene-rich region of the maize genome. Ten genes, at least eight of which are shown to be transcribed, are contained in a 32-kb stretch of DNA that is uninterrupted by retrotransposons. We have isolated nearly full length cDNAs corresponding to the five proximal genes in the cluster. The average intertranscript distance between them is just 1 kb, revealing a surprisingly compact packaging of adjacent genes in this part of the genome. At least 11 small insertions, including several previously described miniature inverted repeat transposable elements, were detected in the introns and 3′ untranslated regions of genes and between genes. The gene-rich region is flanked at the proximal and distal ends by retrotransposon blocks. Thus, the maize genome appears to have scattered regions of high gene density similar to those found in other plants. The unusually high rate of intragenic recombination seen in bz may be related to the very high gene density of the region.

The evolution of the vertebrate Dlx gene family.

The vertebrate Dlx gene family consists of homeobox-containing transcription factors distributed in pairs on the same chromosomes as the Hox genes. To investigate the evolutionary history of Dlx genes, we have cloned five new zebrafish family members and have provided additional sequence information for two mouse genes. Phylogenetic analyses of Dlx gene sequences considered in the context of their chromosomal arrangements suggest that an initial tandem duplication produced a linked pair of Dlx genes after the divergence of chordates and arthropods but prior to the divergence of tunicates and vertebrates. This pair of Dlx genes was then duplicated in the chromosomal events that led to the four clusters of Hox genes characteristic of bony fish and tetrapods. It is possible that a pair of Dlx genes linked to the Hoxc cluster has been lost from mammals. We were unable to distinguish between independent duplication and retention of the ancestral state of bony vertebrates to explain the presence of a greater number of Dlx genes in zebrafish than mammals. Determination of the linkage relationship of these additional zebrafish Dlx genes to Hox clusters should help resolve this issue.

Identification of planarian homeobox sequences indicates the antiquity of most Hox/homeotic gene subclasses.

The homeotic gene complex (HOM-C) is a cluster of genes involved in the anteroposterior axial patterning of animal embryos. It is composed of homeobox genes belonging to the Hox/HOM superclass. Originally discovered in Drosophila, Hox/HOM genes have been identified in organisms as distantly related as arthropods, vertebrates, nematodes, and cnidarians. Data obtained in parallel from the organization of the complex, the domains of gene expression during embryogenesis, and phylogenetic relationships allow the subdivision of the Hox/HOM superclass into five classes (lab, pb/Hox3, Dfd, Antp, and Abd-B) that appeared early during metazoan evolution. We describe a search for homologues of these genes in platyhelminths, triploblast metazoans emerging as an outgroup to the great coelomate ensemble. A degenerate PCR screening for Hox/HOM homeoboxes in three species of triclad planarians has revealed 10 types of Antennapedia-like genes. The homeobox-containing sequences of these PCR fragments allowed the amplification of the homeobox-coding exons for five of these genes in the species Polycelis nigra. A phylogenetic analysis shows that two genes are clear orthologues of Drosophila labial, four others are members of a Dfd/Antp superclass, and a seventh gene, although more difficult to classify with certainty, may be related to the pb/Hox3 class. Together with previously identified Hox/HOM genes in other flatworms, our analyses demonstrate the existence of an elaborate family of Hox/HOM genes in the ancestor of all triploblast animals.

Basis for selection of improved carbohydrate-binding single-chain antibodies from synthetic gene libraries.

A technique is described for the simultaneous and controlled random mutation of all three heavy or light chain complementarity-determining regions (CDRs) in a single-chain Fv specific for the O polysaccharide of Salmonella serogroup B. Sense oligonucleotides were synthesized such that the central bases encoding a CDR were randomized by equimolar spiking with A, G, C, and T at a level of 10% while the antisense strands contained inosine in the spiked regions. Phage display of libraries assembled from the spiked oligonucleotides by a synthetic ligase chain reaction demonstrated a bias for selection of mutants that formed dimers and higher oligomers. Kinetic analyses showed that oligomerization increased association rates in addition to slowing dissociation rates. In combination with some contribution from reduced steric clashes with residues in heavy-chain CDR2, oligomerization resulted in functional affinities that were much higher than that of the monomeric form of the wild-type single-chain Fv.