A numerical study of turbulence seeded with light particles is presented. We analyze the statistical properties of coherent, small-scale structures by looking at the trapping events of light particles inside vortex filaments. We study the properties of particles attracting set, measuring its fractal dimension and the probability that the separation between two particles remains within the dissipative scale, even for time lapses as long as the large-scale correlation time, T(L). We show how to estimate the vortex lifetime by studying the moment of inertia of bunches of particles, showing the presence of an exponential lifetime distribution, with events up to T(L). (C) 2010 American Institute of Physics. [doi:10.1063/1.3431660]
We present the results of our numerical simulations of the Rayleigh-Taylor turbulence, performed using a recently proposed (Sbragaglia et al 2009 J. Fluid Mech. 628 299, Scagliarini et al 2010 Phys. Fluids 22 055101) lattice Boltzmann method that can describe consistently a thermal compressible flow subjected to an external forcing. The method allowed us to study the system in both the nearly Boussinesq regime and the strongly compressible regime. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.
Community structures are found to exist ubiquitously in real-world complex networks. We address here the problem of community detection in directed networks. Most of the previous literature ignores edge directions and applies methods designed for community detection in undirected networks, which discards valuable information and often fails when different communities are defined on the basis of incoming and outgoing edges. We suggest extracting information about edge directions using a PageRank random walk and translating such information into edge weights. After extraction we obtain a new weighted directed network in which edge directions can then be safely ignored. We thus transform community detection in directed networks into community detection in reweighted undirected networks. Such an approach can benefit directly from the large volume of algorithms for the detection of communities in undirected networks already developed, since it is not obvious how to extend these algorithms to account for directed networks and the procedure is often difficult. Validations on synthetic and real-world networks demonstrate that the proposed framework can effectively detect communities in directed networks.
MANIA: A GENE NETWORK REVERSE ALGORITHM FOR COMPOUNDS MODE-OF-ACTION AND GENES INTERACTIONS INFERENCE
Lai Darong
;
Lu Hongtao
;
Lauria Mario
;
Di Bernardo Diego
;
Nardini Christine
Understanding the complexity of the cellular machinery represents a grand challenge in molecular biology. To contribute to the deconvolution of this complexity, a novel inference algorithm based on linear ordinary differential equations is proposed, based solely on high-throughput gene expression data. The algorithm can infer (i) gene-gene interactions from steady state expression profiles and (ii) mode-of-action of the components that can trigger changes in the system. Results demonstrate that the proposed algorithm can identify both information with high performances, thus overcoming the limitation of current algorithms that can infer reliably only one.
With the expanding availability of sequencing technologies, research previously centered on the human genome can now afford to include the study of humans' internal ecosystem (human microbiome). Given the scale of the data involved in this metagenomic research (two orders of magnitude larger than the human genome) and their importance in relation to human health, it is crucial to guarantee (along with the appropriate data collection and taxonomy) proper tools for data analysis. We propose to adapt the approaches defined for the analysis of gene-expression microarray in order to infer information in metagenomics. In particular, we applied SAM, a broadly used tool for the identification of differentially expressed genes among different samples classes, to a reported dataset on a research model with mice of two genotypes (a high density lipoprotein knockout mouse and its wild-type counterpart). The data contain two different diets (high-fat or normal-chow) to ensure the onset of obesity, prodrome of metabolic syndromes (MS). By using 16S rRNA gene as a genomic diversity marker, we illustrate how this approach can identify bacterial populations differentially enriched among different genetic and dietary conditions of the host. This approach faithfully reproduces highly-relevant results from phylogenetic and standard statistical analyses, used to explain the role of the gut microbiome in relation to obesity. This represents a promising proof-of-principle for using functional genomic approaches in the fast growing area of metagenomics, and warrants the availability of a large body of thoroughly tested and theoretically sound methodologies to this exciting new field.
Community structure has been found to exist ubiquitously in many different kinds of real world complex networks. Most of the previous literature ignores edge directions and applies methods designed for community finding in undirected networks to find communities. Here, we address the problem of finding communities in directed networks. Our proposed method uses PageRank random walk induced network embedding to transform a directed network into an undirected one, where the information on edge directions is effectively incorporated into the edge weights. Starting from this new undirected weighted network, previously developed methods for undirected network community finding can be used without any modification. Moreover, our method improves on recent work in terms of community definition and meaning. We provide two simulated examples, a real social network and different sets of power law benchmark networks, to illustrate how our method can correctly detect communities in directed networks. (C) 2010 Elsevier B.V. All rights reserved.
Directed network
Community
Random walk
Network embedding
Modularity
The representation of real systems with network models is becoming increasingly common and critical to both capture and simplify systems' complexity, notably, via the partitioning of networks into communities. In this respect, the definition of modularity, a common and broadly used quality measure for networks partitioning, has induced a surge of efficient modularity-based community detection algorithms. However, recently, the optimization of modularity has been found to show a resolution limit, which reduces its effectiveness and range of applications. Therefore, one recent trend in this area of research has been related to the definition of novel quality functions, alternative to modularity. In this paper, however, instead of laying aside the important body of knowledge developed so far for modularity-based algorithms, we propose to use a strategy to preprocess networks before feeding them into modularity-based algorithms. This approach is based on the observation that dynamic processes triggered on vertices in the same community possess similar behavior patterns but dissimilar on vertices in different communities. Validations on real-world and synthetic networks demonstrate that network preprocessing can enhance the modularity-based community detection algorithms to find more natural clusters and effectively alleviates the problem of resolution limit.
For the numerical solution of the hypersingular integral equation of a notched half-plane problem we propose collocation methods which look for an approximation of the derivative of the solution of the original equation. This derivative is the solution of a Cauchy singular integral equation with additional fixed singularities. We also give a solvability analysis of the original equation which motivates the suggested numerical methods.
Singular integral equations
fixed singularities
collocation
We have developed a rat brain organotypic culture model, in which tissue slices contain cortex-subventricular zone-striatum regions, to model neuroblast activity in response to in vitro ischemia. Neuroblast activation has been described in terms of two main parameters, proliferation and migration from the subventricular zone into the injured cortex. We observed distinct phases of neuroblast activation as is known to occur after in vivo ischemia. Thus, immediately after oxygen/glucose deprivation (6–24 hours), neuroblasts reduce their proliferative and migratory activity, whereas, at longer time points after the insult (2 to 5 days), they start to proliferate and migrate into the damaged cortex. Antagonism of ionotropic receptors for extracellular ATP during and after the insult unmasks an early activation of neuroblasts in the subventricular zone, which responded with a rapid and intense migration of neuroblasts into the damaged cortex (within 24 hours). The process is further enhanced by elevating the production of the chemoattractant SDf-1α and may also be boosted by blocking the activation of microglia. This organotypic model which we have developed is an excellent in vitro system to study neurogenesis after ischemia and other neurodegenerative diseases. Its application has revealed a SOS response to oxygen/glucose deprivation, which is inhibited by unfavorable conditions due to the ischemic environment. Finally, experimental quantifications have allowed us to elaborate a mathematical model to describe neuroblast activation and to develop a computer simulation which should have promising applications for the screening of drug candidates for novel therapies of ischemia-related pathologies.
We provide a polynomial algorithm that determines for any given undirected graph G = (V, E), positive integer k, and convex functions fv : N -> R (v ? V ) a subgraph H = (V, F ) of k edges that minimizes ?v?V fv (dH (v)), where dH (v) is the degree of v in H. The motivation and at the same time the main application of the results is the problem of finding a subset of k vertices in a line graph that covers as many edges as possible. The latter problem generalizes the vertex cover problem for line graphs, which is in turn equivalent to the maximum matching problem in graphs. Improving paths or walks for factorization problems have to be completed by pairs of such walks for this problem. We provide several solutions leading to different variants of the problem and also show the limits of the methods by proving the NP-completeness of some direct extensions, in particular
to all convex functions.
Edge-Path-Tree (EPT) graphs are intersection graphs of EPT matrices that is matrices whose columns
are incidence vectors of edge-sets of paths in a given tree. EPT graphs have polynomially many cliques
[M.C. Golumbic, R.E. Jamison, The edge intersection graphs of paths in a tree, Journal of Combinational
Theory Series B 38 (1985) 8-22; C.L. Monma, V.K. Wey, Intersection graphs of paths in a tree, Journal
of Combinational Theory Series B 41 (1986) 141-181]. Therefore, the problem of finding a clique of
maximum weight in these graphs is solvable in strongly polynomial time. We extend this result to a proper
superclass of EPT graphs.
An Edge Path Tree (EPT) family is a family whose members are edge sets of paths in a
tree. Relying on the notion of Pie introduced in [M.C. Golumbic, R.E. Jamison, The edge
intersection graphs of paths in a tree, Journal of Combinatorial Theory, Series B 38 (1985)
8-22], we characterize Ideal and Mengerian EPT families. In particular, we show that an EPT
family is Ideal if and only if it is Mengerian. If, in addition, the EPT family is uniform, then
it is Ideal if and only if it is Unimodular. The latter equivalence generalizes the well-known
fact that the edge set of a graph is an Ideal clutter if and only if the graph is bipartite