Cloud BioLinux: pre-configured and on-demand bioinformatics computing for the genomics community.

BACKGROUND: A steep drop in the cost of next-generation sequencing during recent years has made the technology affordable to the majority of researchers, but downstream bioinformatic analysis still poses a resource bottleneck for smaller laboratories and institutes that do not have access to substantial computational resources. Sequencing instruments are typically bundled with only the minimal processing and storage capacity required for data capture during sequencing runs. Given the scale of sequence datasets, scientific value cannot be obtained from acquiring a sequencer unless it is accompanied by an equal investment in informatics infrastructure. RESULTS: Cloud BioLinux is a publicly accessible Virtual Machine (VM) that enables scientists to quickly provision on-demand infrastructures for high-performance bioinformatics computing using cloud platforms. Users have instant access to a range of pre-configured command line and graphical software applications, including a full-featured desktop interface, documentation and over 135 bioinformatics packages for applications including sequence alignment, clustering, assembly, display, editing, and phylogeny. Each tool's functionality is fully described in the documentation directly accessible from the graphical interface of the VM. Besides the Amazon EC2 cloud, we have started instances of Cloud BioLinux on a private Eucalyptus cloud installed at the J. Craig Venter Institute, and demonstrated access to the bioinformatic tools interface through a remote connection to EC2 instances from a local desktop computer. Documentation for using Cloud BioLinux on EC2 is available from our project website, while a Eucalyptus cloud image and VirtualBox Appliance is also publicly available for download and use by researchers with access to private clouds. CONCLUSIONS: Cloud BioLinux provides a platform for developing bioinformatics infrastructures on the cloud. An automated and configurable process builds Virtual Machines, allowing the development of highly customized versions from a shared code base. This shared community toolkit enables application specific analysis platforms on the cloud by minimizing the effort required to prepare and maintain them.

The response of marine picoplankton to ocean acidification.

Since industrialization global CO(2) emissions have increased, and as a consequence oceanic pH is predicted to drop by 0.3-0.4 units before the end of the century - a process coined 'ocean acidification'. Consequently, there is significant interest in how pH changes will affect the ocean's biota and integral processes. We investigated marine picoplankton (0.2-2 µm diameter) community response to predicted end of century CO(2) concentrations, via a 'high-CO(2) ' (∼ 750 ppm) large-volume (11 000 l) contained seawater mesocosm approach. We found little evidence of changes occurring in bacterial abundance or community composition due to elevated CO(2) under both phytoplankton pre-bloom/bloom and post-bloom conditions. In contrast, significant differences were observed between treatments for a number of key picoeukaryote community members. These data suggested a key outcome of ocean acidification is a more rapid exploitation of elevated CO(2) levels by photosynthetic picoeukaryotes. Thus, our study indicates the need for a more thorough understanding of picoeukaryote-mediated carbon flow within ocean acidification experiments, both in relation to picoplankton carbon sources, sinks and transfer to higher trophic levels.

Spatial scaling of arbuscular mycorrhizal fungal diversity is affected by farming practice.

Evidence suggests that microbial communities show patterns of spatial scaling which can be driven by geographical distance and environmental heterogeneity. Here we demonstrate that human management can have a major impact on microbial distribution patterns at both the local and landscape scale. Mycorrhizal fungi are vital components of terrestrial ecosystems, forming a mutualistic symbiosis with plant roots which has a major impact on above ground ecology and productivity. We used contrasting agricultural systems to investigate the spatial scaling of the most widespread mycorrhizal fungus group, the arbuscular mycorrhizal fungi (AMF). Using multiple sampling sites with a maximum separation of 250 km we describe for the first time the roles which land management, environmental heterogeneity and geographical distance play in determining spatial patterns of microbial distribution. Analysis of AMF taxa-area relationships at each sampling site revealed that AMF diversity and spatial turnover was greater under organic relative to conventional farm management. At the regional scale (250 km) distance-decay analyses showed that there was significant change in AMF community composition with distance, and that this was greater under organic relative to conventional management. Environmental heterogeneity was found to be the major factor determining turnover of AMF taxa at the landscape scale. Overall we demonstrate that human management can play a key role in determining the turnover of microbial communities at both the local and regional scales.

Oligonucleotide microarrays for bacteriophage expression studies.

Gene expression microarrays offer the ability to monitor the expression of all phage genes over an infection cycle. However, there are relatively few reports to date of microarrays being used to investigate phage biology. This chapter aims to provide an overview of how to design and implement a microarray experiment to investigate phage biology. Given the nature of microarrays being specific to an organism, each will provide a number of unique issues. In this chapter, we outline the basic theory behind microarrays and provide details on how to implement a microarray experiment from the design of oligonucleotide probes through to the hybridisation of microarrays. The matter of designing oligonucleotide probes will be discussed with regards to how probe length, secondary structure, free energy, probe orientation and amplification all have to be taken into account. As means of an example, the conditions used for the hybridisation of an array designed to be specific to the cyanophage S-PM2 is detailed.

Open software for biologists: from famine to feast.

Developing and deploying specialized computing systems for specific research communities is achievable, cost effective and has wide-ranging benefits.

Annotation of environmental OMICS data: application to the transcriptomics domain.

Researchers working on environmentally relevant organisms, populations, and communities are increasingly turning to the application of OMICS technologies to answer fundamental questions about the natural world, how it changes over time, and how it is influenced by anthropogenic factors. In doing so, the need to capture meta-data that accurately describes the biological "source" material used in such experiments is growing in importance. Here, we provide an overview of the formation of the "Env" community of environmental OMICS researchers and its efforts at considering the meta-data capture needs of those working in environmental OMICS. Specifically, we discuss the development to date of the Env specification, an informal specification including descriptors related to geographic location, environment, organism relationship, and phenotype. We then describe its application to the description of environmental transcriptomic experiments and how we have used it to extend the Minimum Information About a Microarray Experiment (MIAME) data standard to create a domain-specific extension that we have termed MIAME/Env. Finally, we make an open call to the community for participation in the Env Community and its future activities.

The architecture of variant surface glycoprotein gene expression sites in Trypanosoma brucei.

Trypanosoma brucei evades the immune system by switching between Variant Surface Glycoprotein (VSG) genes. The active VSG gene is transcribed in one of approximately 20 telomeric expression sites (ESs). It has been postulated that ES polymorphism plays a role in host adaptation. To gain more insight into ES architecture, we have determined the complete sequence of Bacterial Artificial Chromosomes (BACs) containing DNA from three ESs and their flanking regions. There was variation in the order and number of ES-associated genes (ESAGs). ESAGs 6 and 7, encoding transferrin receptor subunits, are the only ESAGs with functional copies in every ES that has been sequenced until now. A BAC clone containing the VO2 ES sequences comprised approximately half of a 330 kb 'intermediate' chromosome. The extensive similarity between this intermediate chromosome and the left telomere of T. brucei 927 chromosome I, suggests that this previously uncharacterised intermediate size class of chromosomes could have arisen from breakage of megabase chromosomes. Unexpected conservation of sequences, including pseudogenes, indicates that the multiple ESs could have arisen through a relatively recent amplification of a single ES.

Unveiling the transcriptional features associated with coccolithovirus infection of natural Emiliania huxleyi blooms.

Lytic viruses have been implicated in the massive cellular lysis observed during algal blooms, through which they assume a prominent role in oceanic carbon and nutrient flows. Despite their impact on biogeochemical cycling, the transcriptional dynamics of these important oceanic events is still poorly understood. Here, we employ an oligonucleotide microarray to monitor host (Emiliania huxleyi) and virus (coccolithovirus) transcriptomic features during the course of E. huxleyi blooms induced in seawater-based mesocosm enclosures. Host bloom development and subsequent coccolithovirus infection was associated with a major shift in transcriptional profile. In addition to the expected metabolic requirements typically associated with viral infection (amino acid and nucleotide metabolism, as well as transcription- and replication-associated functions), the results strongly suggest that the manipulation of lipid metabolism plays a fundamental role during host-virus interaction. The results herein reveal the scale, so far massively underestimated, of the transcriptional domination that occurs during coccolithovirus infection in the natural environment.

Pyrosequencing of Mytilus galloprovincialis cDNAs: tissue-specific expression patterns.

BACKGROUND: Mytilus species are important in marine ecology and in environmental quality assessment, yet their molecular biology is poorly understood. Molecular aspects of their reproduction, hybridisation between species, mitochondrial inheritance, skewed sex ratios of offspring and adaptation to climatic and pollution factors are priority areas. METHODOLOGY/PRINCIPAL FINDINGS: To start to address this situation, expressed genetic transcripts from M. galloprovincialis were pyrosequenced. Transcripts were isolated from the digestive gland, foot, gill and mantle of both male and female mussels. In total, 175,547 sequences were obtained and for foot and mantle, 90% of the sequences could be assembled into contiguous fragments but this reduced to 75% for the digestive gland and gill. Transcripts relating to protein metabolism and respiration dominated including ribosomal proteins, cytochrome oxidases and NADH dehydrogenase subunits. Tissue specific variation was identified in transcripts associated with mitochondrial energy metabolism, with the digestive gland and gill having the greatest transcript abundance. Using fragment recruitment it was also possible to identify sites of potential small RNAs involved in mitochondrial transcriptional regulation. Sex ratios based on Vitelline Envelop Receptor for Lysin and Vitelline Coat Lysin transcript abundances, indicated that an equal sex distribution was maintained. Taxonomic profiling of the M. galloprovincialis tissues highlighted an abundant microbial flora associated with the digestive gland. Profiling of the tissues for genes involved in intermediary metabolism demonstrated that the gill and digestive gland were more similar to each other than to the other two tissues, and specifically the foot transcriptome was most dissimilar. CONCLUSIONS: Pyrosequencing has provided extensive genomic information for M. galloprovincialis and generated novel observations on expression of different tissues, mitochondria and associated microorganisms. It will also facilitate the much needed production of an oligonucleotide microarray for the organism.

Bias in assessments of marine microbial biodiversity in fosmid libraries as evaluated by pyrosequencing.

On the basis of 16S rRNA gene sequencing, the SAR11 clade of marine bacteria has an almost universal distribution, being detected as abundant sequences in all marine provinces. Yet, SAR11 sequences are rarely detected in fosmid libraries, suggesting that the widespread abundance may be an artefact of PCR cloning and that SAR11 has a relatively low abundance. Here the relative abundance of SAR11 is explored in both a fosmid library and a metagenomic sequence data set from the same biological community taken from fjord surface water from Bergen, Norway. Pyrosequenced data and 16S clone data confirmed an 11-15% relative abundance of SAR11 within the community. In contrast, not a single SAR11 fosmid was identified in a pooled shotgun sequence data set of 100 fosmid clones. This underrepresentation was evidenced by comparative abundances of SAR11 sequences assessed by taxonomic annotation and fragment recruitment. Analysis revealed a similar underrepresentation of low-GC Flavobacteriaceae. We speculate that a contributing factor towards the fosmid bias may be DNA fragmentation during preparation because of the low GC content of SAR11 sequences and other underrepresented taxa. This study suggests that, although fosmid libraries can be extremely useful, caution must be taken when directly inferring community composition from metagenomic fosmid libraries.