It is elemental: soil nutrient stoichiometry drives bacterial diversity.

It is well established that resource quantity and elemental stoichiometry play major roles in shaping below and aboveground plant biodiversity, but their importance for shaping microbial diversity in soil remains unclear. Here, we used statistical modeling on a regional database covering 179 locations and six ecosystem types across Scotland to evaluate the roles of total carbon (C), nitrogen (N) and phosphorus (P) availabilities and ratios, together with land use, climate and biotic and abiotic factors, in determining regional scale patterns of soil bacterial diversity. We found that bacterial diversity and composition were primarily driven by variation in soil resource stoichiometry (total C:N:P ratios), itself linked to different land uses, and secondarily driven by other important biodiversity drivers such as climate, soil spatial heterogeneity, soil pH, root influence (plant-soil microbe interactions) and microbial biomass (soil microbe-microbe interactions). In aggregate, these findings provide evidence that nutrient stoichiometry is a strong predictor of bacterial diversity and composition at a regional scale.

The Role of Microbial Community Composition in Controlling Soil Respiration Responses to Temperature.

Rising global temperatures may increase the rates of soil organic matter decomposition by heterotrophic microorganisms, potentially accelerating climate change further by releasing additional carbon dioxide (CO2) to the atmosphere. However, the possibility that microbial community responses to prolonged warming may modify the temperature sensitivity of soil respiration creates large uncertainty in the strength of this positive feedback. Both compensatory responses (decreasing temperature sensitivity of soil respiration in the long-term) and enhancing responses (increasing temperature sensitivity) have been reported, but the mechanisms underlying these responses are poorly understood. In this study, microbial biomass, community structure and the activities of dehydrogenase and ß-glucosidase enzymes were determined for 18 soils that had previously demonstrated either no response or varying magnitude of enhancing or compensatory responses of temperature sensitivity of heterotrophic microbial respiration to prolonged cooling. The soil cooling approach, in contrast to warming experiments, discriminates between microbial community responses and the consequences of substrate depletion, by minimising changes in substrate availability. The initial microbial community composition, determined by molecular analysis of soils showing contrasting respiration responses to cooling, provided evidence that the magnitude of enhancing responses was partly related to microbial community composition. There was also evidence that higher relative abundance of saprophytic Basidiomycota may explain the compensatory response observed in one soil, but neither microbial biomass nor enzymatic capacity were significantly affected by cooling. Our findings emphasise the key importance of soil microbial community responses for feedbacks to global change, but also highlight important areas where our understanding remains limited.

The effect of elevated carbon dioxide on the interaction between Eucalyptus grandis and diverse isolates of Pisolithus sp. is associated with a complex shift in the root transcriptome.

Using the newly available genome for Eucalyptus grandis, we sought to determine the genome-wide traits that enable this host to form mutualistic interactions with ectomycorrhizal (ECM) Pisolithus sp. and to determine how future predicted concentrations of atmospheric carbon dioxide (CO2 ) will affect this relationship. We analyzed the physiological and transcriptomic responses of E. grandis during colonization by different Pisolithus sp. isolates under conditions of ambient (400 ppm) and elevated (650 ppm) CO2 to tease out the gene expression profiles associated with colonization status. We demonstrate that E. grandis varies in its susceptibility to colonization by different Pisolithus isolates in a manner that is not predictable by geographic origin or the internal transcribed spacer (ITS)-based phylogeny of the fungal partner. Elevated concentrations of CO2 alter the receptivity of E. grandis to Pisolithus, a change that is correlated to a dramatic shift in the transcriptomic profile of the root. These data provide a starting point for understanding how future environmental change may alter the signaling between plants and their ECM partners and is a step towards determining the mechanism behind previously observed shifts in Eucalypt-associated fungal communities exposed to elevated concentrations of atmospheric CO2 .

Differential gene expression and alternative splicing between diploid and tetraploid watermelon.

The exploitation of synthetic polyploids for producing seedless fruits is well known in watermelon. Tetraploid progenitors of triploid watermelon plants, compared with their diploid counterparts, exhibit wide phenotypic differences. Although many factors modulate alternative splicing (AS) in plants, the effects of autopolyploidization on AS are still unknown. In this study, we used tissues of leaf, stem, and fruit of diploid and tetraploid sweet watermelon to understand changes in gene expression and the occurrence of AS. RNA-sequencing analysis was performed along with reverse transcription quantitative PCR and rapid amplification of cDNA ends (RACE)-PCR to demonstrate changes in expression and splicing. All vegetative tissues except fruit showed an increased level of AS in the tetraploid watermelon throughout the growth period. The ploidy levels of diploids and the tetraploid were confirmed using a ploidy analyser. We identified 5362 and 1288 genes that were up- and downregulated, respectively, in tetraploid as compared with diploid plants. We further confirmed that 22 genes underwent AS events across tissues, indicating possibilities of generating different protein isoforms with altered functions of important transcription factors and transporters. Arginine biosynthesis, chlorophyllide synthesis, GDP mannose biosynthesis, trehalose biosynthesis, and starch and sucrose degradation pathways were upregulated in autotetraploids. Phloem protein 2, chloroplastic PGR5-like protein, zinc-finger protein, fructokinase-like 2, MYB transcription factor, and nodulin MtN21 showed AS in fruit tissues. These results should help in developing high-quality seedless watermelon and provide additional transcriptomic information related to other cucurbits.

Loss of microbial diversity in soils is coincident with reductions in some specialized functions.

Loss of microbial diversity is considered a major threat because of its importance for ecosystem functions, but there is a lack of conclusive evidence that diversity itself is reduced under anthropogenic stress, and about the consequences of diversity loss. Heavy metals are one of the largest, widespread pollutant types globally, and these represent a significant environmental stressor for terrestrial microbial communities. Using combined metagenomics and functional assays, we show that the compositional and functional response of microbial communities to long-term heavy metal stress results in a significant loss of diversity. Our results indicate that even at a moderate loss of diversity, some key specialized functions (carried out by specific groups) may be compromised. Together with previous work, our data suggest disproportionate impact of contamination on microbes that carry out specialized, but essential, ecosystem functions. Based on these findings, we propose a conceptual framework to explicitly consider diversity of functions and microbial functional groups to test the relationship between biodiversity and soil functions.

Methylome-wide association study of schizophrenia: identifying blood biomarker signatures of environmental insults.

IMPORTANCE: Epigenetic studies present unique opportunities to advance schizophrenia research because they can potentially account for many of its clinical features and suggest novel strategies to improve disease management. OBJECTIVE: To identify schizophrenia DNA methylation biomarkers in blood. DESIGN, SETTING, AND PARTICIPANTS: The sample consisted of 759 schizophrenia cases and 738 controls (N = 1497) collected in Sweden. We used methyl-CpG-binding domain protein-enriched genome sequencing of the methylated genomic fraction, followed by next-generation DNA sequencing. We obtained a mean (SD) number of 68 (26.8) million reads per sample. This massive data set was processed using a specifically designed data analysis pipeline. Critical top findings from our methylome-wide association study (MWAS) were replicated in independent case-control participants using targeted pyrosequencing of bisulfite-converted DNA. MAIN OUTCOMES AND MEASURES: Status of schizophrenia cases and controls. RESULTS: Our MWAS suggested a considerable number of effects, with 25 sites passing the highly conservative Bonferroni correction and 139 sites significant at a false discovery rate of 0.01. Our top MWAS finding, which was located in FAM63B, replicated with P = 2.3 × 10-10. It was part of the networks regulated by microRNA that can be linked to neuronal differentiation and dopaminergic gene expression. Many other top MWAS results could be linked to hypoxia and, to a lesser extent, infection, suggesting that a record of pathogenic events may be preserved in the methylome. Our findings also implicated a site in RELN, one of the most frequently studied candidates in methylation studies of schizophrenia. CONCLUSIONS AND RELEVANCE: To our knowledge, the present study is one of the first MWASs of disease with a large sample size using a technology that provides good coverage of methylation sites across the genome. Our results demonstrated one of the unique features of methylation studies that can capture signatures of environmental insults in peripheral tissues. Our MWAS suggested testable hypotheses about disease mechanisms and yielded biomarkers that can potentially be used to improve disease management.

Ethylene and jasmonic acid act as negative modulators during mutualistic symbiosis between Laccaria bicolor and Populus roots.

The plant hormones ethylene, jasmonic acid and salicylic acid have interconnecting roles during the response of plant tissues to mutualistic and pathogenic symbionts. We used morphological studies of transgenic- or hormone-treated Populus roots as well as whole-genome oligoarrays to examine how these hormones affect root colonization by the mutualistic ectomycorrhizal fungus Laccaria bicolor S238N. We found that genes regulated by ethylene, jasmonic acid and salicylic acid were regulated in the late stages of the interaction between L. bicolor and poplar. Both ethylene and jasmonic acid treatments were found to impede fungal colonization of roots, and this effect was correlated to an increase in the expression of certain transcription factors (e.g. ETHYLENE RESPONSE FACTOR1) and a decrease in the expression of genes associated with microbial perception and cell wall modification. Further, we found that ethylene and jasmonic acid showed extensive transcriptional cross-talk, cross-talk that was opposed by salicylic acid signaling. We conclude that ethylene and jasmonic acid pathways are induced late in the colonization of root tissues in order to limit fungal growth within roots. This induction is probably an adaptive response by the plant such that its growth and vigor are not compromised by the fungus.

A comprehensive family-based replication study of schizophrenia genes.

IMPORTANCE: Schizophrenia (SCZ) is a devastating psychiatric condition. Identifying the specific genetic variants and pathways that increase susceptibility to SCZ is critical to improve disease understanding and address the urgent need for new drug targets. OBJECTIVE: To identify SCZ susceptibility genes. DESIGN: We integrated results from a meta-analysis of 18 genome-wide association studies (GWAS) involving 1,085,772 single-nucleotide polymorphisms (SNPs) and 6 databases that showed significant informativeness for SCZ. The 9380 most promising SNPs were then specifically genotyped in an independent family-based replication study that, after quality control, consisted of 8107 SNPs. SETTING: Linkage meta-analysis, brain transcriptome meta-analysis, candidate gene database, OMIM, relevant mouse studies, and expression quantitative trait locus databases. PATIENTS: We included 11,185 cases and 10,768 control subjects from 6 databases and, after quality control 6298 individuals (including 3286 cases) from 1811 nuclear families. MAIN OUTCOMES AND MEASURES: Case-control status for SCZ. RESULTS: Replication results showed a highly significant enrichment of SNPs with small P values. Of the SNPs with replication values of P.01, the proportion of SNPs that had the same direction of effects as in the GWAS meta-analysis was 89% in the combined ancestry group (sign test, P < 2.20 x 10(-16) and 93% in subjects of European ancestry only (P < 2.20 < 10(-16)). Our results supported the major histocompatibility complex region showing a3.7-fold overall enrichment of replication values of P < .01 in subjects from European ancestry. We replicated SNPs in TCF4 (P = 2.53 x 10(-10)) and NOTCH4 (P = 3.16 x 10(-7)) that are among the most robust SCZ findings. More novel findings included POM121L2 (P = 3.51 x 10(-7)), AS3MT (P = 9.01 x 10(-7)), CNNM2 (P = 6.07 = 10(-7)), and NT5C2(P = 4.09 x 10(-7)). To explore the many small effects, we performed pathway analyses. The most significant pathways involved neuronal function (axonal guidance, neuronal systems, and L1 cell adhesion molecule interaction)and the immune system (antigen processing, cell adhesion molecules relevant to T cells, and translocation to immunological synapse). CONCLUSIONS AND RELEVANCE: We replicated novel SCZ disease genes and pathogenic pathways. Better understanding the molecular and biological mechanisms involved with schizophrenia may improve disease management and may identify new drug targets.

Methylome-wide comparison of human genomic DNA extracted from whole blood and from EBV-transformed lymphocyte cell lines.

DNA from Epstein-Barr virus-transformed lymphocyte cell lines (LCLs) has proven useful for studies of genetic sequence polymorphisms. Whether LCL DNA is suitable for methylation studies is less clear. We conduct a genome-wide methylation investigation using an array set with 45 million probes to investigate the methylome of LCL DNA and technical duplicates of WB DNA from the same 10 individuals. We focus specifically on methylation sites that show variation between individuals and, therefore, are potentially useful as biomarkers. The sample correlations for the methylation variable probes ranged from 0.69 to 0.78 for the WB duplicates and from 0.27 to 0.72 for WB vs LCL. To compare the pattern of the methylation signals, we grouped adjacent probes based on their inter-correlations. These analyses showed ∼29 000 and ∼14 000 blocks in WB and LCL, respectively. Merely 31% of the methylated regions detected in WB were detectable in LCLs. Furthermore, we observed significant differences in mean difference between WB and LCL as compared with duplicates of WB (P-value =2.2 × 10(-16)). Our study shows that there are substantial differences in the DNA methylation patterns between LCL and WB. Thus, LCL DNA should not be used as a proxy for WB DNA in methylome-wide studies.

Genome-wide pharmacogenomic study of neurocognition as an indicator of antipsychotic treatment response in schizophrenia.

Neurocognitive deficits are a core feature of schizophrenia and, therefore, represent potentially critical outcome variables for assessing antipsychotic treatment response. We performed genome-wide association studies (GWAS) with 492K single nucleotide polymorphisms (SNPs) in a sample of 738 patients with schizophrenia from the Clinical Antipsychotic Trials of Intervention Effectiveness study. Outcome variables consisted of a neurocognitive battery administered at multiple time points over an 18-month period, measuring processing speed, verbal memory, vigilance, reasoning, and working memory domains. Genetic mediation of improvements in each of these five domains plus a composite neurocognitive measure was assessed for each of five antipsychotics (olanzapine, perphenazine, quetiapine, risperidone, and ziprasidone). Six SNPs achieved genome-wide significance using a pre-specified threshold that ensures, on average, only 1 in 10 findings is a false discovery. These six SNPs were located within, or in close proximity to, genes EHF, SLC26A9, DRD2, GPR137B, CHST8, and IL1A. The more robust findings, that is those significant across multiple neurocognitive domains and having adjacent SNPs showing evidence for association, were rs286913 at the EHF gene (p-value 6.99 × 10(-8), q-value 0.034, mediating the effects of ziprasidone on vigilance), rs11240594 at SLC26A9 (p-value 1.4 × 10(-7), q-value 0.068, mediating the effects of olanzapine on processing speed), and rs11677416 at IL1A (p-value 6.67 × 10(-7), q-value 0.081, mediating the effects of olanzapine on working memory). This study has generated several novel candidate genes for antipsychotic response. However, our findings will require replication and functional validation. To facilitate replication efforts, we provide all GWAS p-values for download.

Mining mammalian transcript data for functional long non-coding RNAs.

BACKGROUND: The role of long non-coding RNAs (lncRNAs) in controlling gene expression has garnered increased interest in recent years. Sequencing projects, such as Fantom3 for mouse and H-InvDB for human, have generated abundant data on transcribed components of mammalian cells, the majority of which appear not to be protein-coding. However, much of the non-protein-coding transcriptome could merely be a consequence of 'transcription noise'. It is therefore essential to use bioinformatic approaches to identify the likely functional candidates in a high throughput manner. PRINCIPAL FINDINGS: We derived a scheme for classifying and annotating likely functional lncRNAs in mammals. Using the available experimental full-length cDNA data sets for human and mouse, we identified 78 lncRNAs that are either syntenically conserved between human and mouse, or that originate from the same protein-coding genes. Of these, 11 have significant sequence homology. We found that these lncRNAs exhibit: (i) patterns of codon substitution typical of non-coding transcripts; (ii) preservation of sequences in distant mammals such as dog and cow, (iii) significant sequence conservation relative to their corresponding flanking regions (in 50% cases, flanking regions do not have homology at all; and in the remaining, the degree of conservation is significantly less); (iv) existence mostly as single-exon forms (8/11); and, (v) presence of conserved and stable secondary structure motifs within them. We further identified orthologous protein-coding genes that are contributing to the pool of lncRNAs; of which, genes implicated in carcinogenesis are significantly over-represented. CONCLUSION: Our comparative mammalian genomics approach coupled with evolutionary analysis identified a small population of conserved long non-protein-coding RNAs (lncRNAs) that are potentially functional across Mammalia. Additionally, our analysis indicates that amongst the orthologous protein-coding genes that produce lncRNAs, those implicated in cancer pathogenesis are significantly over-represented, suggesting that these lncRNAs could play an important role in cancer pathomechanisms.

Genomic assessment of the evolution of the prion protein gene family in vertebrates.

Prion diseases are devastating neurological disorders caused by the propagation of particles containing an alternative beta-sheet-rich form of the prion protein (PrP). Genes paralogous to PrP, called Doppel and Shadoo, have been identified, that also have neuropathological relevance. To aid in the further functional characterization of PrP and its relatives, we annotated completely the PrP gene family (PrP-GF), in the genomes of 42 vertebrates, through combined strategic application of gene prediction programs and advanced remote homology detection techniques (such as HMMs, PSI-TBLASTN and pGenThreader). We have uncovered several previously undescribed paralogous genes and pseudogenes. We find that current high-quality genomic evidence indicates that the PrP relative Doppel, was likely present in the last common ancestor of present-day Tetrapoda, but was lost in the bird lineage, since its divergence from reptiles. Using the new gene annotations, we have defined the consensus of structural features that are characteristic of the PrP and Doppel structures, across diverse Tetrapoda clades. Furthermore, we describe in detail a transcribed pseudogene derived from Shadoo that is conserved across primates, and that overlaps the meiosis gene, SYCE1, thus possibly regulating its expression. In addition, we analysed the locus of PRNP/PRND for significant conservation across the genomic DNA of eleven mammals, and determined the phylogenetic penetration of non-coding exons. The genomic evidence indicates that the second PRNP non-coding exon found in even-toed ungulates and rodents, is conserved in all high-coverage genome assemblies of primates (human, chimp, orang utan and macaque), and is, at least, likely to have fallen out of use during primate speciation. Furthermore, we have demonstrated that the PRNT gene (at the PRNP human locus) is conserved across at least sixteen mammals, and evolves like a long non-coding RNA, fashioned from fragments of ancient, long, interspersed elements. These annotations and evolutionary analyses will be of further use for functional characterisation of the PrP-GF, and will be updatable in a semi-automated fashion as more genomes accumulate.

Strong association between pseudogenization mechanisms and gene sequence length.

UNLABELLED: Pseudogenes arise from the decay of gene copies following either RNA-mediated duplication (processed pseudogenes) or DNA-mediated duplication (nonprocessed pseudogenes). Here, we show that long protein-coding genes tend to produce more nonprocessed pseudogenes than short genes, whereas the opposite is true for processed pseudogenes. Protein-coding genes longer than 3000 bp are 6 times more likely to produce nonprocessed pseudogenes than processed ones. REVIEWERS: This article was reviewed by Dr. Dan Graur and Dr. Craig Nelson (nominated by Dr. J Peter Gogarten).

Assessing the genomic evidence for conserved transcribed pseudogenes under selection.

BACKGROUND: Transcribed pseudogenes are copies of protein-coding genes that have accumulated indicators of coding-sequence decay (such as frameshifts and premature stop codons), but nonetheless remain transcribed. Recent experimental evidence indicates that transcribed pseudogenes may regulate the expression of homologous genes, through antisense interference, or generation of small interfering RNAs (siRNAs). Here, we assessed the genomic evidence for such transcribed pseudogenes of potential functional importance, in the human genome. The most obvious indicators of such functional importance are significant evidence of conservation and selection pressure. RESULTS: A variety of pseudogene annotations from multiple sources were pooled and filtered to obtain a subset of sequences that have significant mid-sequence disablements (frameshifts and premature stop codons), and that have clear evidence of full-length mRNA transcription. We found 1750 such transcribed pseudogene annotations (TPAs) in the human genome (corresponding to approximately 11.5% of human pseudogene annotations). We checked for syntenic conservation of TPAs in other mammals (rhesus monkey, mouse, rat, dog and cow). About half of the human TPAs are conserved in rhesus monkey, but strikingly, very few in mouse (approximately 3%). The TPAs conserved in rhesus monkey show evidence of selection pressure (relative to surrounding intergenic DNA) on: (i) their GC content, and (ii) their rate of nucleotide substitution. This is in spite of distributions of Ka/Ks (ratios of non-synonymous to synonymous substitution rates), congruent with a lack of protein-coding ability. Furthermore, we have identified 68 human TPAs that are syntenically conserved in at least two other mammals. Interestingly, we observe three TPA sequences conserved in dog that have intermediate character (i.e., evidence of both protein-coding ability and pseudogenicity), and discuss the implications of this. CONCLUSION: Through evolutionary analysis, we have identified candidate sequences for functional human transcribed pseudogenes, and have pinpointed 68 strong candidates for further investigation as potentially functional transcribed pseudogenes across multiple mammal species.

The missing link: Bordetella petrii is endowed with both the metabolic versatility of environmental bacteria and virulence traits of pathogenic Bordetellae.

BACKGROUND: Bordetella petrii is the only environmental species hitherto found among the otherwise host-restricted and pathogenic members of the genus Bordetella. Phylogenetically, it connects the pathogenic Bordetellae and environmental bacteria of the genera Achromobacter and Alcaligenes, which are opportunistic pathogens. B. petrii strains have been isolated from very different environmental niches, including river sediment, polluted soil, marine sponges and a grass root. Recently, clinical isolates associated with bone degenerative disease or cystic fibrosis have also been described. RESULTS: In this manuscript we present the results of the analysis of the completely annotated genome sequence of the B. petrii strain DSMZ12804. B. petrii has a mosaic genome of 5,287,950 bp harboring numerous mobile genetic elements, including seven large genomic islands. Four of them are highly related to the clc element of Pseudomonas knackmussii B13, which encodes genes involved in the degradation of aromatics. Though being an environmental isolate, the sequenced B. petrii strain also encodes proteins related to virulence factors of the pathogenic Bordetellae, including the filamentous hemagglutinin, which is a major colonization factor of B. pertussis, and the master virulence regulator BvgAS. However, it lacks all known toxins of the pathogenic Bordetellae. CONCLUSION: The genomic analysis suggests that B. petrii represents an evolutionary link between free-living environmental bacteria and the host-restricted obligate pathogenic Bordetellae. Its remarkable metabolic versatility may enable B. petrii to thrive in very different ecological niches.

Dynamics of reductive genome evolution in mitochondria and obligate intracellular microbes.

Reductive evolution in mitochondria and obligate intracellular microbes has led to a significant reduction in their genome size and guanine plus cytosine content (GC). We show that genome shrinkage during reductive evolution in prokaryotes follows an exponential decay pattern and provide a method to predict the extent of this decay on an evolutionary timescale. We validated predictions by comparison with estimated extents of genome reduction known to have occurred in mitochondria and Buchnera aphidicola, through comparative genomics and by drawing on available fossil evidences. The model shows how the mitochondrial ancestor would have quickly shed most of its genome, shortly after its incorporation into the protoeukaryotic cell and prior to codivergence subsequent to the split of eukaryotic lineages. It also predicts that the primary rickettsial parasitic event would have occurred between 180 and 425 million years ago (MYA), an event of relatively recent evolutionary origin considering the fact that Rickettsia and mitochondria evolved from a common alphaproteobacterial ancestor. This suggests that the symbiotic events of Rickettsia and mitochondria originated at different time points. Moreover, our model results predict that the ancestor of Wigglesworthia glossinidia brevipalpis, dated around the time of origin of its symbiotic association with the tsetse fly (50-100 MYA), was likely to have been an endosymbiont itself, thus supporting an earlier proposition that Wigglesworthia, which is currently a maternally inherited primary endosymbiont, evolved from a secondary endosymbiont.

Genome sequence of the ubiquitous hydrocarbon-degrading marine bacterium Alcanivorax borkumensis.

Alcanivorax borkumensis is a cosmopolitan marine bacterium that uses oil hydrocarbons as its exclusive source of carbon and energy. Although barely detectable in unpolluted environments, A. borkumensis becomes the dominant microbe in oil-polluted waters. A. borkumensis SK2 has a streamlined genome with a paucity of mobile genetic elements and energy generation-related genes, but with a plethora of genes accounting for its wide hydrocarbon substrate range and efficient oil-degradation capabilities. The genome further specifies systems for scavenging of nutrients, particularly organic and inorganic nitrogen and oligo-elements, biofilm formation at the oil-water interface, biosurfactant production and niche-specific stress responses. The unique combination of these features provides A. borkumensis SK2 with a competitive edge in oil-polluted environments. This genome sequence provides the basis for the future design of strategies to mitigate the ecological damage caused by oil spills.

"Plasmo2D": an ancillary proteomic tool to aid identification of proteins from Plasmodium falciparum.

Bioinformatics tools to aid gene and protein sequence analysis have become an integral part of biology in the post-genomic era. Release of the Plasmodium falciparum genome sequence has allowed biologists to define the gene and the predicted protein content as well as their sequences in the parasite. Using pI and molecular weight as characteristics unique to each protein, we have developed a bioinformatics tool to aid identification of proteins from Plasmodium falciparum. The tool makes use of a Virtual 2-DE generated by plotting all of the proteins from the Plasmodium database on a pI versus molecular weight scale. Proteins are identified by comparing the position of migration of desired protein spots from an experimental 2-DE and that on a virtual 2-DE. The procedure has been automated in the form of user-friendly software called "Plasmo2D". The tool can be downloaded from http://144.16.89.25/Plasmo2D.zip.

Novel hydrolase diversity retrieved from a metagenome library of bovine rumen microflora.

A metagenome expression library of bulk DNA extracted from the rumen content of a dairy cow was established in a phage lambda vector and activity-based screening employed to explore the functional diversity of the microbial flora. Twenty-two clones specifying distinct hydrolytic activities (12 esterases, nine endo-beta-1,4-glucanases and one cyclodextrinase) were identified in the library and characterized. Sequence analysis of the retrieved enzymes revealed that eight (36%) were entirely new and formed deep-branched phylogenetic lineages with no close relatives among known ester- and glycosyl-hydrolases. Bioinformatic analyses of the hydrolase gene sequences, and the sequences and contexts of neighbouring genes, suggested tentative phylogenetic assignments of the rumen organisms producing the retrieved enzymes. The phylogenetic novelty of the hydrolases suggests that some of them may have potential for new applications in biocatalysis.