Electromechanical heterogeneity of the myocardium

Herein we discuss modem data showing that ventricle's working myocardium is highly heterogeneous. Significant transmural differences in electrophysiological and biomechanical properties of cardiomyocytes are reviewed. The reviewed evidence of myocardial heterogeneity constitutes the basis for modem assessment of segmental kinetics of different regions in intact heart. We used muscle duplexes as condensed models of a heterogeneous myocardial system. Experimental data, presented here were obtained both in biological duplexes formed by isolated myocardial preparations and in mathematical models of muscle duplexes. We showed that specific functional heterogeneity of cardiomyocytes, related to their excitation sequence, allowed the myocardium to optimise its contractile function and smooth dispersion of repolarisation.     

Coexistence through spatio-temporal heterogeneity and species sorting in grassland plant communities

The effect of spatial heterogeneity on species coexistence relies on the degree of niche heterogeneity in the habitat and the ability of species to exploit the available niche opportunities. We studied species coexistence in a perennial grassland, and tested whether small-scale disturbances create environmental heterogeneity that affects coexistence and whether the functional diversity of species in the species pool affects the ability of community composition to reflect heterogeneity through species sorting. We manipulated the spatio-temporal heterogeneity of disturbance and the functional diversity of species added as seed and measured their impact on the spatial turnover of species composition. Disturbance increased environmental heterogeneity and spatial turnover, and the effect of heterogeneity on turnover was greatest in the presence of a functionally diverse species pool, showing the importance of trait variation among species for exploiting environmental heterogeneity, and suggesting that coexistence occurred due to species sorting among heterogeneous niches.     

Matrix heterogeneity affects population size of the harvest mice: Bayesian estimation of matrix resistance and model validation

It has been increasingly recognized that landscape matrices are an important factor determining patch connectivity and hence the population size of organisms living in highly fragmented landscapes. However, most previous studies estimated the effect of matrix heterogeneity using prior information regarding dispersal or habitat preferences of a focal organism. Here we estimated matrix resistance of harvest mice in agricultural landscapes using a novel pattern‐oriented modeling with Bayesian estimation and no prior information, and then conducted model validation using data sets independent from those used for model construction. First, we investigated the distribution patterns of harvest mice for approximately 400 habitat patches, and estimated matrix resistance for different matrix types using statistical models incorporating patch size, patch environment, and patch connectivity. We used Bayesian estimation with a Markov chain Monte Carlo algorithm, and searched for appropriate matrix resistance that best explained the distribution pattern. Patch connectivity as well as patch quality was an important determinant of local population size for the harvest mice. Moreover, matrix resistance was far from uniform, with rice and crop fields exhibiting low resistance and forests, creeks, roads and residential areas showing much higher resistance. The deviance explained by this model (heterogeneous matrix model) was much larger than that obtained by the model with no consideration of matrix heterogeneity (homogeneous matrix model). Second, we obtained distribution data from five additional landscapes that were more fragmented than that used for model construction, and used them for model validation. The heterogeneous matrix model well predicted the population size for four out of five landscapes. In contrast, the homogeneous model considerably overestimated population sizes in all cases. Our approach is widely applicable to species living in fragmented landscapes, especially those for which prior information regarding movement or dispersal is difficult to obtain.     

Metabotypes of breast cancer cell lines revealed by non-targeted metabolomics

We present an analysis of intracellular metabolism by non-targeted, high-throughput metabolomics profiling of 18 breast cell lines. We profiled >900 putatively annotated metabolite ions for >100 samples collected under both normoxic and hypoxic conditions and revealed extensive heterogeneity across all metabolic pathways and cell lines. Cell line–specific metabolome profiles dominated over patterns associated with malignancy or with the clinical nomenclature of breast cancer cells. Such characteristic metabolome profiles were reproducible across different laboratories and experiments and exhibited mild to robust changes with change in experimental conditions. To extract a functional overview of cell line heterogeneity, we devised an unsupervised metabotyping procedure that for each pathway automatically recognized metabolic types from metabolome data and assigned cell lines. Our procedure provided a condensed yet global representation of cell line metabolism, revealing the fine structure of metabolic heterogeneity across all tested pathways and cell lines. In follow-up experiments on selected pathways, we confirmed that different metabolic types correlated to differences in the underlying fluxes and difference sensitivity to gene knockdown or pharmacological inhibition. Thus, the identified metabotypes recapitulated functional differences at the pathway level. Metabotyping provides a powerful compression of multi-dimensional data that preserves functional information and serves as a resource for reconciling or understanding heterogeneous metabolic phenotypes or response to inhibition of metabolic pathways.     

Predicting co-complexed protein pairs from heterogeneous data

Proteins do not carry out their functions alone. Instead, they often act by participating in macromolecular complexes and play different functional roles depending on the other members of the complex. It is therefore interesting to identify co-complex relationships. Although protein complexes can be identified in a high-throughput manner by experimental technologies such as affinity purification coupled with mass spectrometry (APMS), these large-scale datasets often suffer from high false positive and false negative rates. Here, we present a computational method that predicts co-complexed protein pair (CCPP) relationships using kernel methods from heterogeneous data sources. We show that a diffusion kernel based on random walks on the full network topology yields good performance in predicting CCPPs from protein interaction networks. In the setting of direct ranking, a diffusion kernel performs much better than the mutual clustering coefficient. In the setting of SVM classifiers, a diffusion kernel performs much better than a linear kernel. We also show that combination of complementary information improves the performance of our CCPP recognizer. A summation of three diffusion kernels based on two-hybrid, APMS, and genetic interaction networks and three sequence kernels achieves better performance than the sequence kernels or diffusion kernels alone. Inclusion of additional features achieves a still better ROC(50) of 0.937. Assuming a negative-to-positive ratio of 600ratio1, the final classifier achieves 89.3% coverage at an estimated false discovery rate of 10%. Finally, we applied our prediction method to two recently described APMS datasets. We find that our predicted positives are highly enriched with CCPPs that are identified by both datasets, suggesting that our method successfully identifies true CCPPs. An SVM classifier trained from heterogeneous data sources provides accurate predictions of CCPPs in yeast. This computational method thereby provides an inexpensive method for identifying protein complexes that extends and complements high-throughput experimental data.     

Gene Function Prediction from Functional Association Networks Using Kernel Partial Least Squares Regression

With the growing availability of large-scale biological datasets, automated methods of extracting functionally meaningful information from this data are becoming increasingly important. Data relating to functional association between genes or proteins, such as co-expression or functional association, is often represented in terms of gene or protein networks. Several methods of predicting gene function from these networks have been proposed. However, evaluating the relative performance of these algorithms may not be trivial: concerns have been raised over biases in different benchmarking methods and datasets, particularly relating to non-independence of functional association data and test data. In this paper we propose a new network-based gene function prediction algorithm using a commute-time kernel and partial least squares regression (Compass). We compare Compass to GeneMANIA, a leading network-based prediction algorithm, using a number of different benchmarks, and find that Compass outperforms GeneMANIA on these benchmarks. We also explicitly explore problems associated with the non-independence of functional association data and test data. We find that a benchmark based on the Gene Ontology database, which, directly or indirectly, incorporates information from other databases, may considerably overestimate the performance of algorithms exploiting functional association data for prediction.     

Methylation patterns of the E-cadherin 5' CpG island are unstable and reflect the dynamic, heterogeneous loss of E-cadherin expression during metastatic progression

Metastatic progression of most common epithelial tumors involves a heterogeneous, transient loss of expression of the homotypic cell adhesion protein, E-cadherin, rather than the uniform loss of a functional protein resulting from coding region mutation. Indeed, whereas E-cadherin loss may promote invasion, reexpression may facilitate cell survival within metastatic deposits. The mechanisms underlying such plasticity are unclear. We now show that the heterogeneous loss of E-cadherin expression in primary human breast cancers reflects a heterogeneous pattern of promoter region methylation, which begins early prior to invasion. In cultured human tumor cells, such heterogeneous methylation is dynamic, varying from allele to allele and shifting in relation to the tumor microenvironment. Following invasion in vitro, which favors diminished E-cadherin expression, the density of promoter methylation markedly increased. When these cells were cultured as spheroids, which requires homotypic cell adhesion, promoter methylation decreased dramatically, and E-cadherin was reexpressed. These data show that the methylation associated with E-cadherin loss in human breast cancer is heterogeneous and unstable and suggest that such epigenetic plasticity may contribute to the dynamic, phenotypic heterogeneity that drives metastatic progression.     

Incorporation of inter-individual heterogeneity into the multistage carcinogenesis model: approach to the analysis of cancer incidence data

We investigate a multistage carcinogenesis frailty model to incorporate inter-individual heterogeneity into carcinogenic response. Attention is focused on inference concerning the effects of different sources of population heterogeneity on cancer rates. The authors consider unobserved variability arising from either carcinogen exposure or background characteristics. Gamma and Inverse-Gaussian distributions are selected for frailty models, and the baseline hazard function is the generalized Armitage-Doll model (i.e. non-frailty model) in which exposure effects shift the age scale instead of acting multiplicatively on cancer rates. For illustration, we apply the method to solid cancer data from a cohort of atomic bomb survivors to examine some features of proposed models. The results show that the Gamma frailty model for the heterogeneity of baseline rates provides the best goodness-of-fit of the model and a non-zero frailty variance. Parameter estimates are, for the most part, comparable between the Gamma and Inverse-Gaussian frailty models. In a heterogeneous population the exposure effects on young adulthood cancer rates might be underestimated for the non-frailty model. Meaningful information regarding each source of heterogeneity has been provided by the proposed method. Therefore, the multistage carcinogenesis frailty model approach is useful for analyses of epidemiological cancer data to assess population heterogeneity and heterogeneity-influenced exposure effects.     

Conformational heterogeneity and probability distributions from single-particle cryo-electron microscopy

Single-particle cryo-electron microscopy (cryo-EM) is a technique that takes projection images of biomolecules frozen at cryogenic temperatures. A major advantage of this technique is its ability to image single biomolecules in heterogeneous conformations. While this poses a challenge for data analysis, recent algorithmic advances have enabled the recovery of heterogeneous conformations from the noisy imaging data. Here, we review methods for the reconstruction and heterogeneity analysis of cryo-EM images, ranging from linear-transformation-based methods to nonlinear deep generative models. We overview the dimensionality-reduction techniques used in heterogeneous 3D reconstruction methods and specify what information each method can infer from the data. Then, we review the methods that use cryo-EM images to estimate probability distributions over conformations in reduced subspaces or predefined by atomistic simulations. We conclude with the ongoing challenges for the cryo-EM community.     

Functional Knowledge Transfer for High-accuracy Prediction of Under-studied Biological Processes

A key challenge in genetics is identifying the functional roles of genes in pathways. Numerous functional genomics techniques (e.g. machine learning) that predict protein function have been developed to address this question. These methods generally build from existing annotations of genes to pathways and thus are often unable to identify additional genes participating in processes that are not already well studied. Many of these processes are well studied in some organism, but not necessarily in an investigator''s organism of interest. Sequence-based search methods (e.g. BLAST) have been used to transfer such annotation information between organisms. We demonstrate that functional genomics can complement traditional sequence similarity to improve the transfer of gene annotations between organisms. Our method transfers annotations only when functionally appropriate as determined by genomic data and can be used with any prediction algorithm to combine transferred gene function knowledge with organism-specific high-throughput data to enable accurate function prediction.We show that diverse state-of-art machine learning algorithms leveraging functional knowledge transfer (FKT) dramatically improve their accuracy in predicting gene-pathway membership, particularly for processes with little experimental knowledge in an organism. We also show that our method compares favorably to annotation transfer by sequence similarity. Next, we deploy FKT with state-of-the-art SVM classifier to predict novel genes to 11,000 biological processes across six diverse organisms and expand the coverage of accurate function predictions to processes that are often ignored because of a dearth of annotated genes in an organism. Finally, we perform in vivo experimental investigation in Danio rerio and confirm the regulatory role of our top predicted novel gene, wnt5b, in leftward cell migration during heart development. FKT is immediately applicable to many bioinformatics techniques and will help biologists systematically integrate prior knowledge from diverse systems to direct targeted experiments in their organism of study.     

Extracting a cellular hierarchy from high-dimensional cytometry data with SPADE

The ability to analyze multiple single-cell parameters is critical for understanding cellular heterogeneity. Despite recent advances in measurement technology, methods for analyzing high-dimensional single-cell data are often subjective, labor intensive and require prior knowledge of the biological system. To objectively uncover cellular heterogeneity from single-cell measurements, we present a versatile computational approach, spanning-tree progression analysis of density-normalized events (SPADE). We applied SPADE to flow cytometry data of mouse bone marrow and to mass cytometry data of human bone marrow. In both cases, SPADE organized cells in a hierarchy of related phenotypes that partially recapitulated well-described patterns of hematopoiesis. We demonstrate that SPADE is robust to measurement noise and to the choice of cellular markers. SPADE facilitates the analysis of cellular heterogeneity, the identification of cell types and comparison of functional markers in response to perturbations.     

Intermediate growth forms as a model for the study of plant clonality functioning: an example with root sprouters

In this contribution we want to show that growth forms intermediate between non-clonal and clonal plants can be used to ask questions about the functional ecology of clonality. We discuss this idea on plants sprouting adventitiously from roots and accomplishing clonal growth via root spacers. Based on extensive literature dealing with growth forms of root sprouting plants, we characterise forms functionally intermediate between clonal root-sprouters and non-clonal plants. We delimit them according to their potential ability to form adventitious shoots and horizontal roots. By reviewing experimental work with root sprouters, we identify the most important triggering factors and developmental constraints influencing these intermediate forms plant age, life-history mode and life-history stage. Using this information we ask questions about the importance of root sprouting in (1) conditions of unpredictable disturbance, where root-sprouting ability may be viewed as a tool for vegetative regeneration, and in (2) temporarily and spatially heterogeneous environment, where foraging by roots may serve as a way of exploiting patchy resources.     

Recursive cluster elimination based support vector machine for disease state prediction using resting state functional and effective brain connectivity

Background

Brain state classification has been accomplished using features such as voxel intensities, derived from functional magnetic resonance imaging (fMRI) data, as inputs to efficient classifiers such as support vector machines (SVM) and is based on the spatial localization model of brain function. With the advent of the connectionist model of brain function, features from brain networks may provide increased discriminatory power for brain state classification.

Methodology/Principal Findings

In this study, we introduce a novel framework where in both functional connectivity (FC) based on instantaneous temporal correlation and effective connectivity (EC) based on causal influence in brain networks are used as features in an SVM classifier. In order to derive those features, we adopt a novel approach recently introduced by us called correlation-purged Granger causality (CPGC) in order to obtain both FC and EC from fMRI data simultaneously without the instantaneous correlation contaminating Granger causality. In addition, statistical learning is accelerated and performance accuracy is enhanced by combining recursive cluster elimination (RCE) algorithm with the SVM classifier. We demonstrate the efficacy of the CPGC-based RCE-SVM approach using a specific instance of brain state classification exemplified by disease state prediction. Accordingly, we show that this approach is capable of predicting with 90.3% accuracy whether any given human subject was prenatally exposed to cocaine or not, even when no significant behavioral differences were found between exposed and healthy subjects.

Conclusions/Significance

The framework adopted in this work is quite general in nature with prenatal cocaine exposure being only an illustrative example of the power of this approach. In any brain state classification approach using neuroimaging data, including the directional connectivity information may prove to be a performance enhancer. When brain state classification is used for disease state prediction, our approach may aid the clinicians in performing more accurate diagnosis of diseases in situations where in non-neuroimaging biomarkers may be unable to perform differential diagnosis with certainty.     

Probabilistic protein function prediction from heterogeneous genome-wide data

Dramatic improvements in high throughput sequencing technologies have led to a staggering growth in the number of predicted genes. However, a large fraction of these newly discovered genes do not have a functional assignment. Fortunately, a variety of novel high-throughput genome-wide functional screening technologies provide important clues that shed light on gene function. The integration of heterogeneous data to predict protein function has been shown to improve the accuracy of automated gene annotation systems. In this paper, we propose and evaluate a probabilistic approach for protein function prediction that integrates protein-protein interaction (PPI) data, gene expression data, protein motif information, mutant phenotype data, and protein localization data. First, functional linkage graphs are constructed from PPI data and gene expression data, in which an edge between nodes (proteins) represents evidence for functional similarity. The assumption here is that graph neighbors are more likely to share protein function, compared to proteins that are not neighbors. The functional linkage graph model is then used in concert with protein domain, mutant phenotype and protein localization data to produce a functional prediction. Our method is applied to the functional prediction of Saccharomyces cerevisiae genes, using Gene Ontology (GO) terms as the basis of our annotation. In a cross validation study we show that the integrated model increases recall by 18%, compared to using PPI data alone at the 50% precision. We also show that the integrated predictor is significantly better than each individual predictor. However, the observed improvement vs. PPI depends on both the new source of data and the functional category to be predicted. Surprisingly, in some contexts integration hurts overall prediction accuracy. Lastly, we provide a comprehensive assignment of putative GO terms to 463 proteins that currently have no assigned function.     

"Guilt by association" is the exception rather than the rule in gene networks

Gene networks are commonly interpreted as encoding functional information in their connections. An extensively validated principle called guilt by association states that genes which are associated or interacting are more likely to share function. Guilt by association provides the central top-down principle for analyzing gene networks in functional terms or assessing their quality in encoding functional information. In this work, we show that functional information within gene networks is typically concentrated in only a very few interactions whose properties cannot be reliably related to the rest of the network. In effect, the apparent encoding of function within networks has been largely driven by outliers whose behaviour cannot even be generalized to individual genes, let alone to the network at large. While experimentalist-driven analysis of interactions may use prior expert knowledge to focus on the small fraction of critically important data, large-scale computational analyses have typically assumed that high-performance cross-validation in a network is due to a generalizable encoding of function. Because we find that gene function is not systemically encoded in networks, but dependent on specific and critical interactions, we conclude it is necessary to focus on the details of how networks encode function and what information computational analyses use to extract functional meaning. We explore a number of consequences of this and find that network structure itself provides clues as to which connections are critical and that systemic properties, such as scale-free-like behaviour, do not map onto the functional connectivity within networks.     

A systematic prediction of multiple drug-target interactions from chemical, genomic, and pharmacological data

In silico prediction of drug-target interactions from heterogeneous biological data can advance our system-level search for drug molecules and therapeutic targets, which efforts have not yet reached full fruition. In this work, we report a systematic approach that efficiently integrates the chemical, genomic, and pharmacological information for drug targeting and discovery on a large scale, based on two powerful methods of Random Forest (RF) and Support Vector Machine (SVM). The performance of the derived models was evaluated and verified with internally five-fold cross-validation and four external independent validations. The optimal models show impressive performance of prediction for drug-target interactions, with a concordance of 82.83%, a sensitivity of 81.33%, and a specificity of 93.62%, respectively. The consistence of the performances of the RF and SVM models demonstrates the reliability and robustness of the obtained models. In addition, the validated models were employed to systematically predict known/unknown drugs and targets involving the enzymes, ion channels, GPCRs, and nuclear receptors, which can be further mapped to functional ontologies such as target-disease associations and target-target interaction networks. This approach is expected to help fill the existing gap between chemical genomics and network pharmacology and thus accelerate the drug discovery processes.     

Integrating biological information into the statistical analysis and design of microarray experiments

Microarray technology is a powerful tool for animal functional genomics studies, with applications spanning from gene identification and mapping, to function and control of gene expression. Microarray assays, however, are complex and costly, and hence generally performed with relatively small number of animals. Nevertheless, they generate data sets of unprecedented complexity and dimensionality. Therefore, such trials require careful planning and experimental design, in addition to tailored statistical and computational tools for their appropriate data mining. In this review, we discuss experimental design and data analysis strategies, which incorporate prior genomic and biological knowledge, such as genotypes and gene function and pathway membership. We focus the discussion on the design of genetical genomics studies, and on significance testing for detection of differential expression. It is shown that the use of prior biological information can improve the efficiency of microarray experiments.     

Demography when history matters: construction and analysis of second-order matrix population models

History matters when individual prior conditions contain important information about the fate of individuals. We present a general framework for demographic models which incorporates the effects of history on population dynamics. The framework incorporates prior condition into the i-state variable and includes an algorithm for constructing the population projection matrix from information on current state dynamics as a function of prior condition. Three biologically motivated classes of prior condition are included: prior stages, linear functions of current and prior stages, and equivalence classes of prior stages. Taking advantage of the matrix formulation of the model, we show how to calculate sensitivity and elasticity of any demographic outcome. Prior condition effects are a source of inter-individual variation in vital rates, i.e., individual heterogeneity. As an example, we construct and analyze a second-order model of Lathyrus vernus, a long-lived herb. We present population growth rate, the stable population distribution, the reproductive value vector, and the elasticity of λ to changes in the second-order transition rates. We quantify the contribution of prior conditions to the total heterogeneity in the stable population of Lathyrus using the entropy of the stable distribution.     

From Experiment to Theory: What Can We Learn from Growth Curves?

Finding an appropriate functional form to describe population growth based on key properties of a described system allows making justified predictions about future population development. This information can be of vital importance in all areas of research, ranging from cell growth to global demography. Here, we use this connection between theory and observation to pose the following question: what can we infer about intrinsic properties of a population (i.e., degree of heterogeneity, or dependence on external resources) based on which growth function best fits its growth dynamics? We investigate several nonstandard classes of multi-phase growth curves that capture different stages of population growth; these models include hyperbolic–exponential, exponential–linear, exponential–linear–saturation growth patterns. The constructed models account explicitly for the process of natural selection within inhomogeneous populations. Based on the underlying hypotheses for each of the models, we identify whether the population that it best fits by a particular curve is more likely to be homogeneous or heterogeneous, grow in a density-dependent or frequency-dependent manner, and whether it depends on external resources during any or all stages of its development. We apply these predictions to cancer cell growth and demographic data obtained from the literature. Our theory, if confirmed, can provide an additional biomarker and a predictive tool to complement experimental research.