U12-type Spliceosome: Localization and Effects of Splicing Efficiency on Gene Expression

The removal of non-coding sequences, introns, is an essential part of messenger RNA processing. In most metazoan organisms, the U12-type spliceosome processes a subset of introns containing highly conserved recognition sequences. U12-type introns constitute less than 0,5% of all introns and reside preferentially in genes related to information processing functions, as opposed to genes encoding for metabolic enzymes. It has previously been shown that the excision of U12-type introns is inefficient compared to that of U2-type introns, supporting the model that these introns could provide a rate-limiting control for gene expression. The low efficiency of U12-type splicing is believed to have important consequences to gene expression by limiting the production of mature mRNAs from genes containing U12-type introns. The inefficiency of U12-type splicing has been attributed to the low abundance of the components of the U12-type spliceosome in cells, but this hypothesis has not been proven. The aim of the first part of this work was to study the effect of the abundance of the spliceosomal snRNA components on splicing. Cells with a low abundance of the U12-type spliceosome were found to inefficiently process U12-type introns encoded by a transfected construct, but the expression levels of endogenous genes were not found to be affected by the abundance of the U12-type spliceosome. However, significant levels of endogenous unspliced U12-type intron-containing pre-mRNAs were detected in cells. Together these results support the idea that U12-type splicing may limit gene expression in some situations. The inefficiency of U12-type splicing has also promoted the idea that the U12-type spliceosome may control gene expression, limiting the mRNA levels of some U12-type intron-containing genes. While the identities of the primary target genes that contain U12-type introns are relatively well known, little has previously been known about the downstream genes and pathways potentially affected by the efficiency of U12-type intron processing. Here, the effects of U12-type splicing efficiency on a whole organism were studied in a Drosophila line with a mutation in an essential U12-type spliceosome component. Genes containing U12-type introns showed variable gene-specific responses to the splicing defect, which points to variation in the susceptibility of different genes to changes in splicing efficiency. Surprisingly, microarray screening revealed that metabolic genes were enriched among downstream effects, and that the phenotype could largely be attributed to one U12-type intron-containing mitochondrial gene. Gene expression control by the U12-type spliceosome could thus have widespread effects on metabolic functions in the organism. The subcellular localization of the U12-type spliceosome components was studied as a response to a recent dispute on the localization of the U12-type spliceosome. All components studied were found to be nuclear indicating that the processing of U12-type introns occurs within the nucleus, thus clarifying a question central to the field. The results suggest that the U12-type spliceosome can limit the expression of genes that contain U12-type introns in a gene-specific manner. Through its limiting role in pre-mRNA processing, the U12-type splicing activity can affect specific genetic pathways, which in the case of Drosophila are involved in metabolic functions.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Gene therapy: Principles, practice, problems and prospects

The remarkable advances made in recombinant DNA technology over the last two decades have paved way for the use of gene transfer to treat human diseases. Several protocols have been developed for the introduction and expression of genes in humans, but the clinical efficacy has not been conclusively demonstrated in any of them. The eventual success of gene therapy for genetic and acquired disorders depends on the development of better gene transfer vectors for sustained, long term expression of foreign genes as well as a better understanding of the pathophysiology of human diseases, it is heartening to note that some of the gene therapy protocols have found other applications such as the genetic immunization or DNA vaccines, which is being heralded as the third vaccine revolution, Gene therapy is yet to become a dream come true, but the light is seen at the end of the tunnel.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

Identification and functional characterization of a cis-acting positive DNA element regulating CYP 2B1/B2 gene transcription in rat liver

A positive cis-acting DNA element in the near 5'-upstream region of the CYP2B1/B2 genes in rat liver was found to play an important role in the transcription of these genes. An oligonucleotide covering -69 to -98 nt mimicked the gel mobility shift pattern given by the fragment -179 to +29 nt, which was earlier found adequate to confer the regulatory features of this gene. Two major complexes were seen, of which the slower and faster moving complexes became intense under uninduced and Phenobarbitone-induced conditions respectively. Minigene cloned DNA plasmid covering -179 to +181 nt in pUC 19 and Bal 31 mutants derived from this parent were transcribed in whole nuclei and cell free transcription extracts and mutants containing only upto -75 nt of the upstream were poorly transcribed. Transcription extracts from phenobarbitone-injected rat liver nuclei were significantly more active than extracts from uninduced rats in transcribing the minigene constructs. Addition of the oligonucleotide (-69 to -98nt) specifically inhibited the transcription of the minigene construct (-179 to +181 nt) in the cell free transcription system. It is therefore, concluded that the region -69 to -98 nt acts as a positive cis-acting element in the transcription of the CYP2B1/B2 genes and in mediating the inductive effects of phenobarbitone.

Use of protein A gene fusions for the analysis of structure-function relationship of the transactivator protein C of bacteriophage Mu

A sensitive dimerization assay for DNA binding proteins has been developed using gene fusion technology. For this purpose, we have engineered a gene fusion using protein A gene of Staphylococcus aureus and C gene, the late gene transactivator of bacteriophage Mu. The C gene was fused to the 3' end of the gene for protein A to generate an A- C fusion. The overexpressed fusion protein was purified in a single step using immunoglobulin affinity chromatography. Purified fusion protein exhibits DNA binding activity as demonstrated by electrophoretic mobility shift assays. When the fusion protein A-C was mixed with C and analyzed for DNA binding, in addition to C and A-C specific complexes, a single intermediate complex comprising of a heterodimer of C and A-C fusion proteins was observed. Further, the protein A moiety in the fusion protein A-C does not contribute to DNA binding as demonstrated by proteolytic cleavage and circular dichroism (CD) analysis. The assay has also been applied to analyze the DNA binding domain of C protein by generating fusions between protein A and N- and C-terminal deletion mutants of C. The results indicate a role for the region towards the carboxy terminal of the protein in DNA binding. The general applicability of this method is discussed.

Involvement of Synthesis and Phosphorylation of Nuclear Protein Factors That Bind to the Positivecis-Acting Element in the Transcriptional Activation of the CYP2B1/B2 Gene by Phenobarbitonein Vivo

The synthesis and phosphorylation of protein factor(s) that bind to the positivecis-acting element (−69 to −98 nt) of the CYP2B1/B2 gene have been examinedin vivoin the rat. Treatment of rats with cycloheximide, a protein synthetic inhibitor, suppresses basal as well as phenobarbitone-induced levels of CYP2B1/B2 mRNA and its run-on transcription. Under these conditions, complex formation of the nuclear extract with the positive element is also inhibited, as judged by gel shift assays. Treatment of rats with 2-aminopurine, a general protein kinase inhibitor, blocks the phenobarbitone-mediated increase in CYP2B1/B2 mRNA, cell-free transcription of a minigene construct containing the positive element, pP450e179DNA, and binding of nuclear proteins to the positive element. Treatment of rats with okadaic acid, a protein phosphatase inhibitor, mimics the effects of phenobarbitone, but only partially. Thus, both phenobarbitone and okadaic acid individually enhance binding of the nuclear protein(s) to the positive element, cell-free transcription of the minigene construct, and phosphorylation of the not, vert, similar26- and 94-kDa proteins binding to the positive element. But unlike phenobarbitone, okadaic acid is not an inducer of CYP2B1/B2 mRNA or its run-on transcription. Thus, phenobarbitone-responsive positive element interactions constitute only a minimal requirement, and okadaic acid is perhaps not able to bring about the total requirement for activation of CYP2B1/B2 gene transcription that should include interaction between the minimal promoter and further upstream elements. An intriguing feature is the antagonistic effect of okadaic acid on phenobarbitone-mediated effects on CYP2B1/B2 mRNA levels, cell-free and run-on transcription, and nuclear protein binding to the positive element. The reason for this antagonism is not clear. It is concluded that phenobarbitone treatment enhancesin vivothe synthesis and phosphorylation of protein factors binding to the positive element and these constitute a minimal requirement for the transcriptional activation of the CYP2B1/B2 gene.

Genetic studies of the PRP17 gene of Saccharomyces cerevisiae: a domain essential for function maps to a nonconserved region of the protein

The PRP17 gene product is required for the second step of pre-mRNA splicing reactions. The C-terminal half of this protein bears four repeat units with homology to the beta transducin repeat. Missense mutations in three temperature-sensitive prp17 mutants map to a region in the N-terminal half of the protein. We have generated, in vitro, 11 missense alleles at the beta transducin repeat units and find that only one affects function in vivo. A phenotypically silent missense allele at the fourth repeat unit enhances the slow-growing phenotype conferred by an allele at the third repeat, suggesting an interaction between these domains. Although many missense mutations in highly conserved amino acids lack phenotypic effects, deletion analysis suggests an essential role for these units. Only mutations in the N-terminal nonconserved domain of PRP17 are synthetically lethal in combination with mutations in PRP16 and PRP18, two other gene products required for the second splicing reaction. A mutually allele-specific interaction between Prp17 and snr7, with mutations in U5 snRNA, was observed. We therefore suggest that the functional region of Prp17p that interacts with Prp18p, Prp16p, and U5 snRNA is the N terminal region of the protein.

Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

Type III Secretion System of Phytopathogenic Bacterium Pseudomonas syringae:From Gene to Function

The type III secretion system (T3SS) is an essential requirement for the virulence of many Gram-negative bacteria which infect plants, animals and men. Pathogens use the T3SS to deliver effector proteins from the bacterial cytoplasm to the eukaryotic host cells, where the effectors subvert host defenses. The best candidates for directing effector protein traffic are the bacterial type III-associated appendages, called needles or pili. In plant pathogenic bacteria, the best characterized example of a T3SS-associated appendage is the HrpA pilus of the plant pathogen Pseudomonas syringae pv. tomato DC3000. The components of the T3SS in plant pathogens are encoded by a cluster of hrp (hypersensitive reaction and pathogenicity) genes. Two major classes of T3SS-secreted proteins are: harpin proteins such as HrpZ which are exported into extracellular space, and avirulence (Avr) proteins such as AvrPto which are translocated directly to the plant cytoplasm. This study deals with the structural and functional characterization of the T3SS-associated HrpA pilus and the T3SS-secreted harpins. By insertional mutagenesis analysis of HrpA, we located the optimal epitope insertion site in the amino-terminus of HrpA, and revealed the potential application of the HrpA pilus as a carrier of antigenic determinants for vaccination. By pulse-expression of proteins combined with immuno-electron microscopy, we discovered the Hrp pilus assembly strategy as addition of HrpA subunits to the distal end of the growing pilus, and we showed for the first time that secretion of HrpZ occurs at the tip of the pilus. The pilus thus functions as a conduit delivering proteins to the extracellular milieu. By using phage-display and scanning-insertion mutagenesis methods we identified a conserved HrpZ-binding peptide and localized the peptide-binding site to the central domain of HrpZ. We also found that the HrpZ specifically interacts with a host bean protein. Taken together, the current results provide deeper insight into the molecular mechanism of T3SS-associated pilus assembly and effector protein translocation, which will be helpful for further studies on the pathogenic mechanisms of Gram-negative bacteria and for developing new strategies to prevent bacterial infection.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

Human parvovirus B19: tissue persistence and prevalence of prototypic and new variants

Human parvovirus B19 is a minute ssDNA virus causing a wide variety of diseases, including erythema infectiosum, arthropathy, anemias, and fetal death. After primary infection, genomic DNA of B19 has been shown to persist in solid tissues of not only symptomatic but also of constitutionally healthy, immunocompetent individuals. In this thesis, the viral DNA was shown to persist as an apparently intact molecule of full length, and without persistence-specific mutations. Thus, although the mere presence of B19 DNA in tissue can not be used as a diagnostic criterion, a possible role in the pathogenesis of diseases e.g. through mRNA or protein production can not be excluded. The molecular mechanism, the host-cell type and the possible clinical significance of B19 DNA tissue persistence are yet to be elucidated. In the beginning of this work, the B19 genomic sequence was considered highly conserved. However, new variants were found: V9 was detected in 1998 in France, in serum of a child with aplastic crisis. This variant differed from the prototypic B19 sequences by ~10 %. In 2002 we found, persisting in skin of constitutionally healthy humans, DNA of another novel B19 variant, LaLi. Genetically this variant differed from both the prototypic sequences and the variant V9 also by ~10%. Simultaneously, B19 isolates with DNA sequences similar to LaLi were introduced by two other groups, in the USA and France. Based on phylogeny, a classification scheme based on three genotypes (B19 types 1-3) was proposed. Although the B19 virus is mainly transmitted via the respiratory route, blood and plasma-derived products contaminated with high levels of B19 DNA have also been shown to be infectious. The European Pharmacopoeia stipulates that, in Europe, from the beginning of 2004, plasma pools for manufacture must contain less than 104 IU/ml of B19 DNA. Quantitative PCR screening is therefore a prerequisite for restriction of the B19 DNA load and obtaining of safe plasma products. Due to the DNA sequence variation among the three B19 genotypes, however, B19 PCR methods might fail to detect the new variants. We therefore examined the suitability of the two commercially available quantitative B19 PCR tests, LightCycler-Parvovirus B19 quantification kit (Roche Diagnostics) and RealArt Parvo B19 LC PCR (Artus), for detection, quantification and differentiation of the three B19 types known, including B19 types 2 and 3. The former method was highly sensitive for detection of the B19 prototype but was not suitable for detection of types 2 and 3. The latter method detected and differentiated all three B19 virus types. However, one of the two type-3 strains was detected at a lower sensitivity. Then, we assessed the prevalence of the three B19 virus types among Finnish blood donors, by screening pooled plasma samples derived from >140 000 blood-donor units: none of the pools contained detectable levels of B19 virus types 2 or 3. According to the results of other groups, B19 type 2 was absent also among Danish blood-donors, and extremely rare among symptomatic European patients. B19 type 3 has been encountered endemically in Ghana and (apparently) in Brazil, and sporadical cases have been detected in France and the UK. We next examined the biological characteristics of these virus types. The p6 promoter regions of virus types 1-3 were cloned in front of a reporter gene, the constructs were transfected into different cell lines, and the promoter activities were measured. As a result, we found that the activities of the three p6 promoters, although differing in sequence by >20%, were of equal strength, and most active in B19-permissive cells. Furthermore, the infectivity of the three B19 types was examined in two B19-permissive cell lines. RT-PCR revealed synthesis of spliced B19 mRNAs, and immunofluorescence verified the production of NS1 and VP proteins in the infected cells. These experiments suggested similar host-cell tropism and showed that the three virus types are strains of the same species, i.e. human parvovirus B19. Last but not least, the sera from subjects infected in the past either with B19 type 1 or type 2 (as evidenced by tissue persistence of the respective DNAs), revealed in VP1/2- and VP2-EIAs a 100 % cross-reactivity between virus types 1 and 2. These results, together with similar studies by others, indicate that the three B19 genotypes constitute a single serotype.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.