Rapid analysis of metagenomic data using signature-based clustering.

BACKGROUND: Sequencing highly-variable 16S regions is a common and often effective approach to the study of microbial communities, and next-generation sequencing (NGS) technologies provide abundant quantities of data for analysis. However, the speed of existing analysis pipelines may limit our ability to work with these quantities of data. Furthermore, the limited coverage of existing 16S databases may hamper our ability to characterise these communities, particularly in the context of complex or poorly studied environments. RESULTS: In this article we present the SigClust algorithm, a novel clustering method involving the transformation of sequence reads into binary signatures. When compared to other published methods, SigClust yields superior cluster coherence and separation of metagenomic read data, while operating within substantially reduced timeframes. We demonstrate its utility on published Illumina datasets and on a large collection of labelled wound reads sourced from patients in a wound clinic. The temporal analysis is based on tracking the dominant clusters of wound samples over time. The analysis can identify markers of both healing and non-healing wounds in response to treatment. Prominent clusters are found, corresponding to bacterial species known to be associated with unfavourable healing outcomes, including a number of strains of Staphylococcus aureus. CONCLUSIONS: SigClust identifies clusters rapidly and supports an improved understanding of the wound microbiome without reliance on a reference database. The results indicate a promising use for a SigClust-based pipeline in wound analysis and prediction, and a possible novel method for wound management and treatment.

Gene networks associated with non-syndromic intellectual disability.

Non-syndromic intellectual disability (NS-ID) is a genetically heterogeneous disorder, with more than 200 candidate genes to date. Despite the increasing number of novel mutations detected, a relatively low number of recurrently mutated genes have been identified, highlighting the complex genetic architecture of the disorder. A systematic search of PubMed and Medline identified 245 genes harbouring non-synonymous variants, insertions or deletions, which were identified as candidate NS-ID genes from case reports or from linkage or pedigree analyses. From this list, 33 genes are common to syndromic intellectual disability (S-ID) and 58 genes are common to certain neurological and neuropsychiatric disorders that often include intellectual disability as a clinical feature. We examined the evolutionary constraint and brain expression of these gene sets, and we performed gene network and protein-protein interaction analyses using GeneGO MetaCoreTM and DAPPLE, respectively. The 245 NS-ID candidate genes were over-represented in axon guidance, synaptogenesis, cell adhesion and neurotransmission pathways, all of which are key neurodevelopmental processes for the establishment of mature neuronal circuitry in the brain. These 245 genes exhibit significantly elevated expression in human brain and are evolutionarily constrained, consistent with expectations for a brain disorder such as NS-ID that is associated with reduced fecundity. In addition, we report enrichment of dopaminergic and glutamatergic pathways for those candidate NS-ID genes that are common to S-ID and/or neurological and neuropsychiatric disorders that exhibit intellectual disability. Collectively, this study provides an overview and analysis of gene networks associated with NS-ID and suggests modulation of neurotransmission, particularly dopaminergic and glutamatergic systems as key contributors to synaptic dysfunction in NS-ID.

Programmed ribosomal frameshift alters expression of west nile virus genes and facilitates virus replication in birds and mosquitoes.

West Nile virus (WNV) is a human pathogen of significant medical importance with close to 40,000 cases of encephalitis and more than 1,600 deaths reported in the US alone since its first emergence in New York in 1999. Previous studies identified a motif in the beginning of non-structural gene NS2A of encephalitic flaviviruses including WNV which induces programmed -1 ribosomal frameshift (PRF) resulting in production of an additional NS protein NS1'. We have previously demonstrated that mutant WNV with abolished PRF was attenuated in mice. Here we have extended our previous observations by showing that PRF does not appear to have a significant role in virus replication, virion formation, and viral spread in several cell lines in vitro. However, we have also shown that PRF induces an over production of structural proteins over non-structural proteins in virus-infected cells and that mutation abolishing PRF is present in ∼11% of the wild type virus population. In vivo experiments in house sparrows using wild type and PRF mutant of New York 99 strain of WNV viruses showed some attenuation for the PRF mutant virus. Moreover, PRF mutant of Kunjin strain of WNV showed significant decrease compared to wild type virus infection in dissemination of the virus from the midgut through the haemocoel, and ultimately the capacity of infected mosquitoes to transmit virus. Thus our results demonstrate an important role for PRF in regulating expression of viral genes and consequently virus replication in avian and mosquito hosts.

Whole transcriptome characterization of the effects of dehydration and rehydration on Cladonia rangiferina, the grey reindeer lichen.

BACKGROUND: Lichens are symbiotic organisms with a fungal and an algal or a cyanobacterial partner. Lichens inhabit some of the harshest climates on earth and most lichen species are desiccation-tolerant. Lichen desiccation-tolerance has been studied at the biochemical level and through proteomics, but the underlying molecular genetic mechanisms remain largely unexplored. The objective of our study was to examine the effects of dehydration and rehydration on the gene expression of Cladonia rangiferina. RESULTS: Samples of C. rangiferina were collected at several time points during both the dehydration and rehydration process and the gene expression intensities were measured using a custom DNA microarray. Several genes, which were differentially expressed in one or more time points, were identified. The microarray results were validated using qRT-PCR analysis. Enrichment analysis of differentially expressed transcripts was also performed to identify the Gene Ontology terms most associated with the rehydration and dehydration process. CONCLUSIONS: Our data identify differential expression patterns for hundreds of genes that are modulated during dehydration and rehydration in Cladonia rangiferina. These dehydration and rehydration events clearly differ from each other at the molecular level and the largest changes to gene expression are observed within minutes following rehydration. Distinct changes are observed during the earliest stage of rehydration and the mechanisms not appear to be shared with the later stages of wetting or with drying. Several of the most differentially expressed genes are similar to genes identified in previous studies that have investigated the molecular mechanisms of other desiccation-tolerant organisms. We present here the first microarray experiment for any lichen species and have for the first time studied the genetic mechanisms behind lichen desiccation-tolerance at the whole transcriptome level.

Fontimonas thermophila gen. nov., sp. nov., a moderately thermophilic bacterium isolated from a freshwater hot spring, and proposal of Solimonadaceae fam. nov. to replace Sinobacteraceae Zhou et al. 2008.

A novel bacterial strain designated HA-01(T) was isolated from a freshwater terrestrial hot spring located at Hot Springs National Park, Arkansas, USA. Cells were Gram-negative-staining, rod-shaped, aerobic, chemo-organotrophic, oxidase- and catalase-positive, non-spore-forming and motile by means of a single polar flagellum. Growth occurred at 37-60 °C, with an optimum between 45 and 50 °C, and at pH 6.5-8.5, with an optimum between pH 6.5 and 7.0. Phylogenetic analysis based on 16S rRNA gene sequences indicated that the closest relatives of strain HA-01(T) were Solimonas aquatica NAA-16(T) (93.8 %), Solimonas flava CW-KD 4(T) (94.1 %), Solimonas soli DCY12(T) (93.1 %), Solimonas variicoloris MN28(T) (94.0 %), Nevskia ramosa Soe1(T) (91.2 %) and Hydrocarboniphaga effusa AP103(T) (91.1 %). Major fatty acids consisted of C(16 : 0), iso-C(16 : 0), C(16 : 1)ω5c and summed feature 8 (C(18 : 1)ω9c, C(18 : 1)ω7c and C(18 : 1)ω6c). Polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine, and the major isoprenoid quinone was Q-8. The DNA G+C content was 64.4 mol%. Based on phylogenetic, phenotypic and chemotaxonomic evidence, it is proposed that strain HA-01(T) represents a novel species in a new genus for which the name Fontimonas thermophila gen. nov., sp. nov. is proposed. The type strain of the type species is HA-01(T) (= DSM 23609(T) = CCUG 59713(T)). A new family, Solimonadaceae fam. nov., is also proposed to replace Sinobacteriaceae Zhou et al. 2008.

Thermoanaerobaculum aquaticum gen. nov., sp. nov., the first cultivated member of Acidobacteria subdivision 23, isolated from a hot spring.

A novel bacterium was isolated from a freshwater hot spring, the Hale House Spring, located at Hot Springs National Park, Hot Springs, AR, USA. Cells of strain MP-01(T) stained Gram-negative, were rod-shaped, non-motile, strictly anaerobic and chemo-organotrophic and did not form spores. Growth occurred at 50-65 °C, with an optimum at 60 °C, at pH 6.0-8.0, with an optimum at pH 6.5-7.0, and at NaCl concentrations up to 0.5 % (w/v), with optimum growth in the absence of NaCl. Strain MP-01(T) was capable of fermentative growth on pyruvate or proteinaceous substrates as well as reducing Fe(III) and Mn(IV). Major fatty acids were iso-C15 : 0, iso-C16 : 0, anteiso-C17 : 0 and iso-C17 : 0. The polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine and the major isoprenoid quinone was MK-10. In the polyamine pattern, sym-homospermidine was the predominant compound. The DNA G+C content was 62.7 mol%. Analysis of the 16S rRNA gene sequence of the isolate indicated that strain MP-01(T) represents the first reported cultivated member of subdivision 23 of the Acidobacteria. It is proposed that strain MP-01(T) represents a novel genus and species, for which the name Thermoanaerobaculum aquaticum gen. nov., sp. nov. is proposed. The type strain of Thermoanaerobaculum aquaticum is MP-01(T) ( = DSM 24856(T) = JCM 18256(T)).

Draft genome sequences of Helicobacter pylori isolates from Malaysia, cultured from patients with functional dyspepsia and gastric cancer.

Helicobacter pylori is the main bacterial causative agent of gastroduodenal disorders and a risk factor for gastric adenocarcinoma and mucosa-associated lymphoid tissue (MALT) lymphoma. The draft genomes of 10 closely related H. pylori isolates from the multiracial Malaysian population will provide an insight into the genetic diversity of isolates in Southeast Asia. These isolates were cultured from gastric biopsy samples from patients with functional dyspepsia and gastric cancer. The availability of this genomic information will provide an opportunity for examining the evolution and population structure of H. pylori isolates from Southeast Asia, where the East meets the West.

Characterization of a transcriptome from a non-model organism, Cladonia rangiferina, the grey reindeer lichen, using high-throughput next generation sequencing and EST sequence data.

BACKGROUND: Lichens are symbiotic organisms that have a remarkable ability to survive in some of the most extreme terrestrial climates on earth. Lichens can endure frequent desiccation and wetting cycles and are able to survive in a dehydrated molecular dormant state for decades at a time. Genetic resources have been established in lichen species for the study of molecular systematics and their taxonomic classification. No lichen species have been characterised yet using genomics and the molecular mechanisms underlying the lichen symbiosis and the fundamentals of desiccation tolerance remain undescribed. We report the characterisation of a transcriptome of the grey reindeer lichen, Cladonia rangiferina, using high-throughput next-generation transcriptome sequencing and traditional Sanger EST sequencing data. RESULTS: Altogether 243,729 high quality sequence reads were de novo assembled into 16,204 contigs and 49,587 singletons. The genome of origin for the sequences produced was predicted using Eclat with sequences derived from the axenically grown symbiotic partners used as training sequences for the classification model. 62.8% of the sequences were classified as being of fungal origin while the remaining 37.2% were predicted as being of algal origin. The assembled sequences were annotated by BLASTX comparison against a non-redundant protein sequence database with 34.4% of the sequences having a BLAST match. 29.3% of the sequences had a Gene Ontology term match and 27.9% of the sequences had a domain or structural match following an InterPro search. 60 KEGG pathways with more than 10 associated sequences were identified. CONCLUSIONS: Our results present a first transcriptome sequencing and de novo assembly for a lichen species and describe the ongoing molecular processes and the most active pathways in C. rangiferina. This brings a meaningful contribution to publicly available lichen sequence information. These data provide a first glimpse into the molecular nature of the lichen symbiosis and characterise the transcriptional space of this remarkable organism. These data will also enable further studies aimed at deciphering the genetic mechanisms behind lichen desiccation tolerance.

Microarray analysis identifies candidate genes for key roles in coral development.

BACKGROUND: Anthozoan cnidarians are amongst the simplest animals at the tissue level of organization, but are surprisingly complex and vertebrate-like in terms of gene repertoire. As major components of tropical reef ecosystems, the stony corals are anthozoans of particular ecological significance. To better understand the molecular bases of both cnidarian development in general and coral-specific processes such as skeletogenesis and symbiont acquisition, microarray analysis was carried out through the period of early development - when skeletogenesis is initiated, and symbionts are first acquired. RESULTS: Of 5081 unique peptide coding genes, 1084 were differentially expressed (P

Maintenance of ancestral complexity and non-metazoan genes in two basal cnidarians.

Cnidarians are among the simplest extant animals; however EST analyses reveal that they have a remarkably high level of genetic complexity. In this article, we show that the full diversity of metazoan signaling pathways is represented in this phylum, as are antagonists previously known only in chordates. Many of the cnidarian ESTs match genes previously known only in non-animal kingdoms. At least some of these represent ancient genes lost by all bilaterians examined so far, rather than genes gained by recent lateral gene transfer.

Proteomic screen in the simple metazoan Hydra identifies 14-3-3 binding proteins implicated in cellular metabolism, cytoskeletal organisation and Ca2+ signalling.

BACKGROUND: 14-3-3 proteins have been implicated in many signalling mechanisms due to their interaction with Ser/Thr phosphorylated target proteins. They are evolutionarily well conserved in eukaryotic organisms from single celled protozoans and unicellular algae to plants and humans. A diverse array of target proteins has been found in higher plants and in human cell lines including proteins involved in cellular metabolism, apoptosis, cytoskeletal organisation, secretion and Ca2+ signalling. RESULTS: We found that the simple metazoan Hydra has four 14-3-3 isoforms. In order to investigate whether the diversity of 14-3-3 target proteins is also conserved over the whole animal kingdom we isolated 14-3-3 binding proteins from Hydra vulgaris using a 14-3-3-affinity column. We identified 23 proteins that covered most of the above-mentioned groups. We also isolated several novel 14-3-3 binding proteins and the Hydra specific secreted fascin-domain-containing protein PPOD. In addition, we demonstrated that one of the 14-3-3 isoforms, 14-3-3 HyA, interacts with one Hydra-Bcl-2 like protein in vitro. CONCLUSION: Our results indicate that 14-3-3 proteins have been ubiquitous signalling components since the start of metazoan evolution. We also discuss the possibility that they are involved in the regulation of cell numbers in response to food supply in Hydra.

Spatiotemporal expression control correlates with intragenic scaffold matrix attachment regions (S/MARs) in Arabidopsis thaliana.

Scaffold/matrix attachment regions (S/MARs) are essential for structural organization of the chromatin within the nucleus and serve as anchors of chromatin loop domains. A significant fraction of genes in Arabidopsis thaliana contains intragenic S/MAR elements and a significant correlation of S/MAR presence and overall expression strength has been demonstrated. In this study, we undertook a genome scale analysis of expression level and spatiotemporal expression differences in correlation with the presence or absence of genic S/MAR elements. We demonstrate that genes containing intragenic S/MARs are prone to pronounced spatiotemporal expression regulation. This characteristic is found to be even more pronounced for transcription factor genes. Our observations illustrate the importance of S/MARs in transcriptional regulation and the role of chromatin structural characteristics for gene regulation. Our findings open new perspectives for the understanding of tissue- and organ-specific regulation of gene expression.

openSputnik--a database to ESTablish comparative plant genomics using unsaturated sequence collections.

The public expressed sequence tag collections are continually being enriched with high-quality sequences that represent an ever-expanding range of taxonomically diverse plant species. While these sequence collections provide biased insight into the populations of expressed genes available within individual species and their associated tissues, the information is conceivably of wider relevance in a comparative context. When we consider the available expressed sequence tag (EST) collections of summer 2004, most of the major plant taxonomic clades are at least superficially represented. Investigation of the five million available plant ESTs provides a wealth of information that has applications in modelling the routes of plant genome evolution and the identification of lineage-specific genes and gene families. Over four million ESTs from over 50 distinct plant species have been collated within an EST analysis pipeline called openSputnik. The ESTs were resolved down into approximately one million unigene sequences. These have been annotated using orthology-based annotation transfer from reference plant genomes and using a variety of contemporary bioinformatics methods to assign peptide, structural and functional attributes. The openSputnik database is available at http://sputnik.btk.fi.

Eclair--a web service for unravelling species origin of sequences sampled from mixed host interfaces.

The identification of the genes that participate at the biological interface of two species remains critical to our understanding of the mechanisms of disease resistance, disease susceptibility and symbiosis. The sequencing of complementary DNA (cDNA) libraries prepared from the biological interface between two organisms provides an inexpensive way to identify the novel genes that may be expressed as a cause or consequence of compatible or incompatible interactions. Sequence classification and annotation of species origin typically use an orthology-based approach and require access to large portions of either genome, or a close relative. Novel species- or clade-specific sequences may have no counterpart within existing databases and remain ambiguous features. Here we present a web-service, Eclair, which utilizes support vector machines for the classification of the origin of expressed sequence tags stemming from mixed host cDNA libraries. In addition to providing an interface for the classification of sequences, users are presented with the opportunity to train a model to suit their preferred species pair. Eclair is freely available at http://eclair.btk.fi.

PlantMarkers--a database of predicted molecular markers from plants.

Molecular markers are required in a broad spectrum of gene screening approaches, ranging from gene-mapping within traditional 'forward'-genetics approaches through QTL identification studies to genotyping and haplotyping studies. As we enter the post-genomics era, the need for genetic markers does not diminish, even in the species with fully sequenced genomes. PlantMarkers is a genetic marker database that contains a comprehensive pool of predicted molecular markers. We have adopted contemporary techniques to identify putative single nucleotide polymorphism (SNP), simple sequence repeat (SSR) and conserved orthologue set markers. A systematic approach to identify as broad a range of putative markers has been undertaken by screening the available openSputnik unigene consensus sequences from over 50 plant species. A web presence at http://markers.btk.fi provides functionality so that a user may search for species-specific markers on the basis of many specific criteria not limited to non-synonymous SNPs segregating between different varieties or measured polymorphic SSRs. Feedback forms are provided with all sequence entries to enable inclusion of, for example, map location for markers validated by the research community.

Bioinformatics in the plant genomic and phenomic domain: The German contribution to resources, services and perspectives.

Plant genetic resources are a substantial opportunity for plant breeding, preservation and maintenance of biological diversity. As part of the German Network for Bioinformatics Infrastructure (de.NBI) the German Crop BioGreenformatics Network (GCBN) focuses mainly on crop plants and provides both data and software infrastructure which are tailored to the needs of the plant research community. Our mission and key objectives include: (1) provision of transparent access to germplasm seeds, (2) the delivery of improved workflows for plant gene annotation, and (3) implementation of bioinformatics services that link genotypes and phenotypes. This review introduces the GCBN's spectrum of web-services and integrated data resources that address common research problems in the plant genomics community.

Genetic screen for signal peptides in Hydra reveals novel secreted proteins and evidence for non-classical protein secretion.

We have screened a Hydra cDNA library for sequences encoding N-terminal signal peptides using the yeast invertase secretion vector pSUC [Jacobs et al., 1997. A genetic selection for isolating cDNAs encoding secreted proteins. Gene 198, 289-296]. We isolated and sequenced 907 positive clones; 88% encoded signal peptides; 12% lacked signal peptides. By searching the Hydra EST database we identified full-length sequences for the selected clones. These encoded 37 known proteins with signal peptides and 40 novel Hydra-specific proteins with signal peptides. Localization of two signal peptide-containing sequences, VEGF and ferritin, to the secretory pathway was confirmed with GFP fusion proteins. In addition, we isolated 105 clones which lacked signal peptides but which supported invertase secretion from yeast. Isolation of plasmids from these clones and retransformation in invertase-negative yeast cells confirmed the phenotype. A GFP fusion protein of one such clone encoding the foot morphogen pedibin was localized to the cytoplasm in transfected Hydra cells and did not enter the ER/Golgi secretory pathway. Secretion of pedibin and other proteins lacking signal peptides appears to occur by a non-classical protein secretion route.

A wing expressed sequence tag resource for Bicyclus anynana butterflies, an evo-devo model.

BACKGROUND: Butterfly wing color patterns are a key model for integrating evolutionary developmental biology and the study of adaptive morphological evolution. Yet, despite the biological, economical and educational value of butterflies they are still relatively under-represented in terms of available genomic resources. Here, we describe an Expression Sequence Tag (EST) project for Bicyclus anynana that has identified the largest available collection to date of expressed genes for any butterfly. RESULTS: By targeting cDNAs from developing wings at the stages when pattern is specified, we biased gene discovery towards genes potentially involved in pattern formation. Assembly of 9,903 ESTs from a subtracted library allowed us to identify 4,251 genes of which 2,461 were annotated based on BLAST analyses against relevant gene collections. Gene prediction software identified 2,202 peptides, of which 215 longer than 100 amino acids had no homology to any known proteins and, thus, potentially represent novel or highly diverged butterfly genes. We combined gene and Single Nucleotide Polymorphism (SNP) identification by constructing cDNA libraries from pools of outbred individuals, and by sequencing clones from the 3' end to maximize alignment depth. Alignments of multi-member contigs allowed us to identify over 14,000 putative SNPs, with 316 genes having at least one high confidence double-hit SNP. We furthermore identified 320 microsatellites in transcribed genes that can potentially be used as genetic markers. CONCLUSION: Our project was designed to combine gene and sequence polymorphism discovery and has generated the largest gene collection available for any butterfly and many potential markers in expressed genes. These resources will be invaluable for exploring the potential of B. anynana in particular, and butterflies in general, as models in ecological, evolutionary, and developmental genetics.

Sputnik: a database platform for comparative plant genomics.

Two million plant ESTs, from 20 different plant species, and totalling more than one 1000 Mbp of DNA sequence, represents a formidable transcriptomic resource. Sputnik uses the potential of this sequence resource to fill some of the information gap in the un-sequenced plant genomes and to serve as the foundation for in silicio comparative plant genomics. The complexity of the individual EST collections has been reduced using optimised EST clustering techniques. Annotation of cluster sequences is performed by exploiting and transferring information from the comprehensive knowledgebase already produced for the completed model plant genome (Arabidopsis thaliana) and by performing additional state of-the-art sequence analyses relevant to today's plant biologist. Functional predictions, comparative analyses and associative annotations for 500 000 plant EST derived peptides make Sputnik (http://mips.gsf.de/proj/sputnik/) a valid platform for contemporary plant genomics.