Plant protein-coding gene families: emerging bioinformatics approaches

Protein-coding gene families are sets of similar genes with a shared evolutionary origin and, generally, with similar biological functions. In plants, the size and role of gene families has been only partially addressed. However, suitable bioinformatics tools are being developed to cluster the enormous number of sequences currently available in databases. Specifically, comparative genomic databases promise to become powerful tools for gene family annotation in plant clades. In this review, I evaluate the data retrieved from various gene family databases, the ease with which they can be extracted and how useful the extracted information is.

Azotobacter vinelandii nitrogenase: “Kinetics of nif gene expression and insights into the roles of FdxN and NifQ in FeMo-co biosynthesis”

The Molybdenum-nitrogenase is responsible for most biological nitrogen fixation activity (BNF) in the biosphere. Due to its great agronomical importance, it has been the subject of profound genetic and biochemical studies. The Mo nitrogenase carries at its active site a unique iron-molybdenum cofactor (FeMoco) that consists of an inorganic 7 Fe, 1 Mo, 1 C, 9 S core coordinated to the organic acid homocitrate. Biosynthesis of FeMo-co occurs outside nitrogenase through a complex and highly regulated pathway involving proteins acting as molecular scaffolds, metallocluster carriers or enzymes that provide substrates in appropriate chemical forms. Specific expression regulatory factors tightly control the accumulation levels of all these other components. Insertion of FeMo-co into a P-cluster containing apo-NifDK polypeptide results in nitrogenase reconstitution. Investigation of FeMo-co biosynthesis has uncovered new radical chemistry reactions and new roles for Fe-S clusters in biology.

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.

A genomic approach to gene fusion technology

Gene expression profiling provides powerful analyses of transcriptional responses to cellular perturbation. In contrast to DNA array-based methods, reporter gene technology has been underused for this application. Here we describe a genomewide, genome-registered collection of Escherichia coli bioluminescent reporter gene fusions. DNA sequences from plasmid-borne, random fusions of E. coli chromosomal DNA to a Photorhabdus luminescens luxCDABE reporter allowed precise mapping of each fusion. The utility of this collection covering about 30% of the transcriptional units was tested by analyzing individual fusions representative of heat shock, SOS, OxyR, SoxRS, and cya/crp stress-responsive regulons. Each fusion strain responded as anticipated to environmental conditions known to activate the corresponding regulatory circuit. Thus, the collection mirrors E. coli's transcriptional wiring diagram. This genomewide collection of gene fusions provides an independent test of results from other gene expression analyses. Accordingly, a DNA microarray-based analysis of mitomycin C-treated E. coli indicated elevated expression of expected and unanticipated genes. Selected luxCDABE fusions corresponding to these up-regulated genes were used to confirm or contradict the DNA microarray results. The power of partnering gene fusion and DNA microarray technology to discover promoters and define operons was demonstrated when data from both suggested that a cluster of 20 genes encoding production of type I extracellular polysaccharide in E. coli form a single operon.

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.