NeatMap - non-clustering heat map alternatives in R

<p>Abstract</p> <p>Background</p> <p>The clustered heat map is the most popular means of visualizing genomic data. It compactly displays a large amount of data in an intuitive format that facilitates the detection of hidden structures and relations in the data. However, it is hampered by its use of cluster analysis which does not always respect the intrinsic relations in the data, often requiring non-standardized reordering of rows/columns to be performed post-clustering. This sometimes leads to uninformative and/or misleading conclusions. Often it is more informative to use dimension-reduction algorithms (such as Principal Component Analysis and Multi-Dimensional Scaling) which respect the topology inherent in the data. Yet, despite their proven utility in the analysis of biological data, they are not as widely used. This is at least partially due to the lack of user-friendly visualization methods with the visceral impact of the heat map.</p> <p>Results</p> <p>NeatMap is an R package designed to meet this need. NeatMap offers a variety of novel plots (in 2 and 3 dimensions) to be used in conjunction with these dimension-reduction techniques. Like the heat map, but unlike traditional displays of such results, it allows the entire dataset to be displayed while visualizing relations between elements. It also allows superimposition of cluster analysis results for mutual validation. NeatMap is shown to be more informative than the traditional heat map with the help of two well-known microarray datasets.</p> <p>Conclusions</p> <p>NeatMap thus preserves many of the strengths of the clustered heat map while addressing some of its deficiencies. It is hoped that NeatMap will spur the adoption of non-clustering dimension-reduction algorithms.</p

Reordering Hierarchical Tree Based on Bilateral Symmetric Distance

BACKGROUND: In microarray data analysis, hierarchical clustering (HC) is often used to group samples or genes according to their gene expression profiles to study their associations. In a typical HC, nested clustering structures can be quickly identified in a tree. The relationship between objects is lost, however, because clusters rather than individual objects are compared. This results in a tree that is hard to interpret. METHODOLOGY/PRINCIPAL FINDINGS: This study proposes an ordering method, HC-SYM, which minimizes bilateral symmetric distance of two adjacent clusters in a tree so that similar objects in the clusters are located in the cluster boundaries. The performance of HC-SYM was evaluated by both supervised and unsupervised approaches and compared favourably with other ordering methods. CONCLUSIONS/SIGNIFICANCE: The intuitive relationship between objects and flexibility of the HC-SYM method can be very helpful in the exploratory analysis of not only microarray data but also similar high-dimensional data

Going beyond Clustering in MD Trajectory Analysis: An Application to Villin Headpiece Folding

Recent advances in computing technology have enabled microsecond long all-atom molecular dynamics (MD) simulations of biological systems. Methods that can distill the salient features of such large trajectories are now urgently needed. Conventional clustering methods used to analyze MD trajectories suffer from various setbacks, namely (i) they are not data driven, (ii) they are unstable to noise and changes in cut-off parameters such as cluster radius and cluster number, and (iii) they do not reduce the dimensionality of the trajectories, and hence are unsuitable for finding collective coordinates. We advocate the application of principal component analysis (PCA) and a non-metric multidimensional scaling (nMDS) method to reduce MD trajectories and overcome the drawbacks of clustering. To illustrate the superiority of nMDS over other methods in reducing data and reproducing salient features, we analyze three complete villin headpiece folding trajectories. Our analysis suggests that the folding process of the villin headpiece is structurally heterogeneous

Effects of host genetics and environment on egg-associated microbiotas in brown trout (Salmo trutta).

Recent studies found fish egg-specific bacterial communities that changed over the course of embryogenesis, suggesting an interaction between the developing host and its microbiota. Indeed, single-strain infections demonstrated that the virulence of opportunistic bacteria is influenced by environmental factors and host immune genes. However, the interplay between a fish embryo host and its microbiota has not been studied yet at the community level. In order to test whether host genetics affects the assemblage of egg-associated bacteria, adult brown trout (Salmo trutta) were sampled from a natural population. Their gametes were used for full-factorial in vitro fertilizations to separate sire from dam effects. In total 2,520 embryos were singly raised under experimental conditions that differently support microbial growth. High-throughput 16S rRNA amplicon sequencing was applied to characterize bacterial communities on milt and fertilized eggs across treatments. Dam and sire identity influenced embryo mortality, time until hatching, and composition of egg-associated microbiotas, but no link between bacterial communities on milt and on fertilized eggs could be found. Elevated resources increased embryo mortality and modified bacterial communities with a shift in their putative functional potential. Resource availability did not significantly affect any parental effects on embryo performance. Sire identity affected bacterial diversity that turned out to be a significant predictor of hatching time: embryos associated with high bacterial diversity hatched later. We conclude that both host genetics and the availability of resources define diversity and composition of egg-associated bacterial communities that then affect the life-history of their hosts. This article is protected by copyright. All rights reserved

Global gene expression profiling of Plasmodium falciparum in response to the anti-malarial drug pyronaridine

<p>Abstract</p> <p>Background</p> <p>Pyronaridine (PN) and chloroquine (CQ) are structurally related anti-malarial drugs with primarily the same mode of action. However, PN is effective against several multidrug-resistant lines of <it>Plasmodium falciparum</it>, including CQ resistant lines, suggestive of important operational differences between the two drugs.</p> <p>Methods</p> <p>Synchronized trophozoite stage cultures of <it>P. falciparum </it>strain K1 (CQ resistant) were exposed to 50% inhibitory concentrations (IC₅₀) of PN and CQ, and parasites were harvested from culture after 4 and 24 hours exposure. Global transcriptional changes effected by drug treatment were investigated using DNA microarrays.</p> <p>Results</p> <p>After a 4 h drug exposure, PN induced a greater degree of transcriptional perturbation (61 differentially expressed features) than CQ (10 features). More genes were found to respond to 24 h treatments with both drugs, and 461 features were found to be significantly responsive to one or both drugs across all treatment conditions.</p> <p>Filtering was employed to remove features unrelated to primary drug action, specifically features representing genes developmentally regulated, secondary stress/death related processes and sexual stage development. The only significant gene ontologies represented among the 46 remaining features after filtering relate to host exported proteins from multi-gene families.</p> <p>Conclusions</p> <p>The malaria parasite's molecular responses to PN and CQ treatment are similar in terms of the genes and pathways affected. However, PN appears to exert a more rapid response than CQ. The faster action of PN may explain why PN is more efficacious than CQ, particularly against CQ resistant isolates. In agreement with several other microarray studies of drug action on the parasite, it is not possible, however, to discern mechanism of drug action from the drug-responsive genes.</p

Increased diversity of egg-associated bacteria on brown trout (Salmo trutta) at elevated temperatures.

The taxonomic composition of egg-associated microbial communities can play a crucial role in the development of fish embryos. In response, hosts increasingly influence the composition of their associated microbial communities during embryogenesis, as concluded from recent field studies and laboratory experiments. However, little is known about the taxonomic composition and the diversity of egg-associated microbial communities within ecosystems; e.g., river networks. We sampled late embryonic stages of naturally spawned brown trout at nine locations within two different river networks and applied 16S rRNA pyrosequencing to describe their bacterial communities. We found no evidence for a significant isolation-by-distance effect on the composition of bacterial communities, and no association between neutral genetic divergence of fish host (based on 11 microsatellites) and phylogenetic distances of the composition of their associated bacterial communities. We characterized core bacterial communities on brown trout eggs and compared them to corresponding water samples with regard to bacterial composition and its presumptive function. Bacterial diversity was positively correlated with water temperature at the spawning locations. We discuss this finding in the context of the increased water temperatures that have been recorded during the last 25 years in the study area

Transcriptome analysis of the spalax hypoxia survival response includes suppression of apoptosis and tight control of angiogenesis

BACKGROUND: The development of complex responses to hypoxia has played a key role in the evolution of mammals, as inadequate response to this condition is frequently associated with cardiovascular diseases, developmental disorders, and cancers. Though numerous studies have used mice and rats in order to explore mechanisms that contribute to hypoxia tolerance, these studies are limited due to the high sensitivity of most rodents to severe hypoxia. The blind subterranean mole rat Spalax is a hypoxia tolerant rodent, which exhibits unique longevity and therefore has invaluable potential in hypoxia and cancer research. RESULTS: Using microarrays, transcript abundance was measured in brain and muscle tissues from Spalax and rat individuals exposed to acute and chronic hypoxia for varying durations. We found that Spalax global gene expression response to hypoxia differs from that of rat and is characterized by the activation of functional groups of genes that have not been strongly associated with the response to hypoxia in hypoxia sensitive mammals. Using functional enrichment analysis of Spalax hypoxia induced genes we found highly significant overrepresentation of groups of genes involved in anti apoptosis, cancer, embryonic/sexual development, epidermal growth factor receptor binding, coordinated suppression and activation of distinct groups of transcription factors and membrane receptors, in addition to angiogenic related processes. We also detected hypoxia induced increases of different critical Spalax hub gene transcripts, including antiangiogenic genes associated with cancer tolerance in Down syndrome human individuals. CONCLUSIONS: This is the most comprehensive study of Spalax large scale gene expression response to hypoxia to date, and the first to use custom Spalax microarrays. Our work presents novel patterns that may underlie mechanisms with critical importance to the evolution of hypoxia tolerance, with special relevance to medical research

Cultivation of stable, reproducible microbial communities from different fecal donors using minibioreactor arrays (MBRAs)

Background: Continuous-flow culture models are one tool for studying complex interactions between members of human fecal microbiotas because they allow studies to be completed during an extended period of time under conditions where pH, nutrient availability, and washout of waste products and dead cells can be controlled. Because many of the existing well-validated continuous-flow models are large and complex, we were interested in developing a simpler continuous-flow system that would allow microbial community dynamics to be examined in higher throughput while still maintaining complex microbial communities. To this end, we developed minibioreactor arrays (MBRAs), small volume bioreactors (15 ml) that allow simultaneous cultivation of up to 48 microbial communities in a single anaerobic chamber. Results: We used MBRA to characterize the microbial community dynamics of replicate reactors inoculated from three different human fecal donors and reactors seeded with feces pooled from these three donors. We found that MBRA could be used to efficiently cultivate complex microbial communities that were a subset of the initial fecal inoculum (15–25 % of fecal OTUs initially observed). After an initial acclimation period of approximately 1 week, communities in each reactor stabilized and exhibited day-to-day variation similar to that observed in stable mouse fecal communities. Replicate reactors were predominately populated by shared core microbial communities; variation between replicate reactors was primarily driven by shifts in abundance of shared operational taxonomic units (OTUs). Consistent with differences between fecal donors, MBRA communities present in reactors seeded with different fecal samples had distinct composition and structure. Conclusions: From these analyses, we conclude that MBRAs can be used to cultivate communities that recapitulate key features of human fecal communities and are a useful tool to facilitate higher-throughput studies of the dynamics of these communities

A global transcriptional analysis of Plasmodium falciparum malaria reveals a novel family of telomere-associated lncRNAs

Background: Mounting evidence suggests a major role for epigenetic feedback in Plasmodium falciparum transcriptional regulation. Long non-coding RNAs (lncRNAs) have recently emerged as a new paradigm in epigenetic remodeling. We therefore set out to investigate putative roles for lncRNAs in P. falciparum transcriptional regulation. Results: We used a high-resolution DNA tiling microarray to survey transcriptional activity across 22.6% of the P. falciparum strain 3D7 genome. We identified 872 protein-coding genes and 60 putative P. falciparum lncRNAs under developmental regulation during the parasite's pathogenic human blood stage. Further characterization of lncRNA candidates led to the discovery of an intriguing family of lncRNA telomere-associated repetitive element transcripts, termed lncRNA-TARE. We have quantified lncRNA-TARE expression at 15 distinct chromosome ends and mapped putative transcriptional start and termination sites of lncRNA-TARE loci. Remarkably, we observed coordinated and stage-specific expression of lncRNA-TARE on all chromosome ends tested, and two dominant transcripts of approximately 1.5 kb and 3.1 kb transcribed towards the telomere. Conclusions: We have characterized a family of 22 telomere-associated lncRNAs in P. falciparum. Homologous lncRNA-TARE loci are coordinately expressed after parasite DNA replication, and are poised to play an important role in P. falciparum telomere maintenance, virulence gene regulation, and potentially other processes of parasite chromosome end biology. Further study of lncRNA-TARE and other promising lncRNA candidates may provide mechanistic insight into P. falciparum transcriptional regulation.Organismic and Evolutionary BiologyStem Cell and Regenerative BiologyOther Research Uni