Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

On the resilience of long cycles in random graphs

In this paper we determine the local and global resilience of random graphs G(n,p) (p >> n(-1)) with respect to the property of containing a cycle of length at least (1 - alpha)n. Roughly speaking, given alpha > 0, we determine the smallest r(g) (G, alpha) with the property that almost surely every subgraph of G = G(n,p) having more than r(g) (G, alpha)vertical bar E(G)vertical bar edges contains a cycle of length at least (1 - alpha)n (global resilience). We also obtain, for alpha < 1/2, the smallest r(l) (G, alpha) such that any H subset of G having deg(H) (v) larger than r(l) (G, alpha) deg(G) (v) for all v is an element of V(G) contains a cycle of length at least (1 - alpha)n (local resilience). The results above are in fact proved in the more general setting of pseudorandom graphs.

Survival time of random walk in random environment among soft obstacles

We consider a Random Walk in Random Environment (RWRE) moving in an i.i.d. random field of obstacles. When the particle hits an obstacle, it disappears with a positive probability. We obtain quenched and annealed bounds on the tails of the survival time in the general d-dimensional case. We then consider a simplified one-dimensional model (where transition probabilities and obstacles are independent and the RWRE only moves to neighbour sites), and obtain finer results for the tail of the survival time. In addition, we study also the ""mixed"" probability measures (quenched with respect to the obstacles and annealed with respect to the transition probabilities and vice-versa) and give results for tails of the survival time with respect to these probability measures. Further, we apply the same methods to obtain bounds for the tails of hitting times of Branching Random Walks in Random Environment (BRWRE).

Random perturbations of stochastic processes with unbounded variable length memory

We consider binary infinite order stochastic chains perturbed by a random noise. This means that at each time step, the value assumed by the chain can be randomly and independently flipped with a small fixed probability. We show that the transition probabilities of the perturbed chain are uniformly close to the corresponding transition probabilities of the original chain. As a consequence, in the case of stochastic chains with unbounded but otherwise finite variable length memory, we show that it is possible to recover the context tree of the original chain, using a suitable version of the algorithm Context, provided that the noise is small enough.

Gibbs random graphs on point processes

Consider a discrete locally finite subset Gamma of R(d) and the cornplete graph (Gamma, E), with vertices Gamma and edges E. We consider Gibbs measures on the set of sub-graphs with vertices Gamma and edges E` subset of E. The Gibbs interaction acts between open edges having a vertex in common. We study percolation properties of the Gibbs distribution of the graph ensemble. The main results concern percolation properties of the open edges in two cases: (a) when Gamma is sampled from a homogeneous Poisson process; and (b) for a fixed Gamma with sufficiently sparse points. (c) 2010 American Institute of Physics. [doi:10.1063/1.3514605]

An oracle approach for interaction neighborhood estimation in random fields

We consider the problem of interaction neighborhood estimation from the partial observation of a finite number of realizations of a random field. We introduce a model selection rule to choose estimators of conditional probabilities among natural candidates. Our main result is an oracle inequality satisfied by the resulting estimator. We use then this selection rule in a two-step procedure to evaluate the interacting neighborhoods. The selection rule selects a small prior set of possible interacting points and a cutting step remove from this prior set the irrelevant points. We also prove that the Ising models satisfy the assumptions of the main theorems, without restrictions on the temperature, on the structure of the interacting graph or on the range of the interactions. It provides therefore a large class of applications for our results. We give a computationally efficient procedure in these models. We finally show the practical efficiency of our approach in a simulation study.

Identification of archaeal proteins that affect the exosome function in vitro

Background: The archaeal exosome is formed by a hexameric RNase PH ring and three RNA binding subunits and has been shown to bind and degrade RNA in vitro. Despite extensive studies on the eukaryotic exosome and on the proteins interacting with this complex, little information is yet available on the identification and function of archaeal exosome regulatory factors. Results: Here, we show that the proteins PaSBDS and PaNip7, which bind preferentially to poly-A and AU-rich RNAs, respectively, affect the Pyrococcus abyssi exosome activity in vitro. PaSBDS inhibits slightly degradation of a poly-rA substrate, while PaNip7 strongly inhibits the degradation of poly-A and poly-AU by the exosome. The exosome inhibition by PaNip7 appears to depend at least partially on its interaction with RNA, since mutants of PaNip7 that no longer bind RNA, inhibit the exosome less strongly. We also show that FITC-labeled PaNip7 associates with the exosome in the absence of substrate RNA. Conclusions: Given the high structural homology between the archaeal and eukaryotic proteins, the effect of archaeal Nip7 and SBDS on the exosome provides a model for an evolutionarily conserved exosome control mechanism.

Identification of Novel Immunoregulatory Molecules in Human Thymic Regulatory CD4(+)CD25(+) T Cells by Phage Display

Thymic CD4(+)CD25(+) cells play an important role in immune regulation and are continuously developed in the thymus as an independent lineage. How these cells are generated, what are their multiple pathways of suppressive activity and which are their specific markers are questions that remain unanswered. To identify molecules involved in the function and development of human CD4(+)CD25(+) T regulatory cells we targeted thymic CD4(+)CD25(+) cells by peptide phage display. A phage library containing random peptides was screened ex vivo for binding to human thymic CD4(+)CD25(+) T cells. After four rounds of selection on CD4(+)CD25(+) enriched populations of thymocytes, we sequenced several phage displayed peptides and selected one with identity to the Vitamin D Receptor (VDR). We confirmed the binding of the VDR phage to active Vitamin D in vitro, as well as the higher expression of VDR in CD4(+)CD25(+) cells. We suggest that differential expression of VDR on natural Tregs may be related to the relevance of Vitamin D in function and ontogeny of these cells.

The Essential Nucleolar Yeast Protein Nop8p Controls the Exosome Function during 60S Ribosomal Subunit Maturation

The yeast nucleolar protein Nop8p has previously been shown to interact with Nip7p and to be required for 60S ribosomal subunit formation. Although depletion of Nop8p in yeast cells leads to premature degradation of rRNAs, the biochemical mechanism responsible for this phenotype is still not known. In this work, we show that the Nop8p amino-terminal region mediates interaction with the 5.8S rRNA, while its carboxyl-terminal portion interacts with Nip7p and can partially complement the growth defect of the conditional mutant strain Dnop8/GAL:NOP8. Interestingly, Nop8p mediates association of Nip7p to pre-ribosomal particles. Nop8p also interacts with the exosome subunit Rrp6p and inhibits the complex activity in vitro, suggesting that the decrease in 60S ribosomal subunit levels detected upon depletion of Nop8p may result from degradation of pre-rRNAs by the exosome. These results strongly indicate that Nop8p may control the exosome function during pre-rRNA processing.

Structure-Function Analysis of the HrpB2-HrcU Interaction in the Xanthomonas citri Type III Secretion System

Bacterial type III secretion systems deliver protein virulence factors to host cells. Here we characterize the interaction between HrpB2, a small protein secreted by the Xanthomonas citri subsp. citri type III secretion system, and the cytosolic domain of the inner membrane protein HrcU, a paralog of the flagellar protein FlhB. We show that a recombinant fragment corresponding to the C-terminal cytosolic domain of HrcU produced in E. coli suffers cleavage within a conserved Asn264-Pro265-Thr266-His267 (NPTH) sequence. A recombinant HrcU cytosolic domain with N264A, P265A, T266A mutations at the cleavage site (HrcU(AAAH)) was not cleaved and interacted with HrpB2. Furthermore, a polypeptide corresponding to the sequence following the NPTH cleavage site also interacted with HrpB2 indicating that the site for interaction is located after the NPTH site. Non-polar deletion mutants of the hrcU and hrpB2 genes resulted in a total loss of pathogenicity in susceptible citrus plants and disease symptoms could be recovered by expression of HrpB2 and HrcU from extrachromossomal plasmids. Complementation of the Delta hrcU mutant with HrcU(AAAH) produced canker lesions similar to those observed when complemented with wild-type HrcU. HrpB2 secretion however, was significantly reduced in the Delta hrcU mutant complemented with HrcU(AAAH), suggesting that an intact and cleavable NPTH site in HrcU is necessary for total functionally of T3SS in X. citri subsp. citri. Complementation of the Delta hrpB2 X. citri subsp. citri strain with a series of hrpB2 gene mutants revealed that the highly conserved HrpB2 C-terminus is essential for T3SS-dependent development of citrus canker symptoms in planta.

Multifractality in the random parameter model for multivariate time series

The Random Parameter model was proposed to explain the structure of the covariance matrix in problems where most, but not all, of the eigenvalues of the covariance matrix can be explained by Random Matrix Theory. In this article, we explore the scaling properties of the model, as observed in the multifractal structure of the simulated time series. We use the Wavelet Transform Modulus Maxima technique to obtain the multifractal spectrum dependence with the parameters of the model. The model shows a scaling structure compatible with the stylized facts for a reasonable choice of the parameter values. (C) 2009 Elsevier B.V. All rights reserved.

Somesthetic, gustatory, olfactory function and salivary flow in patients with neuropathic trigeminal pain

OBJECTIVES: To determine somesthetic, olfactory, gustative and salivary abnormalities in patients with burning mouth syndrome (BMS), idiopathic trigeminal neuralgia (ITN) and trigeminal postherpetic neuralgia (PHN). SUBJECTS AND METHODS: Twenty patients from each group (BMS, ITN, PHN) and 60 healthy controls were evaluated with a systematized quantitative approach of thermal (cold and warm), mechanical, pain, gustation, olfaction and salivary flow; data were analyzed with ANOVA, Tukey, Kruskal Wallis and Dunn tests with a level of significance of 5%. RESULTS: There were no salivary differences among the groups with matched ages; the cold perception was abnormal only at the mandibular branch of PHN (P = 0.001) and warm was abnormal in all trigeminal branches of PHN and BMS; mechanical sensitivity was altered at the mandibular branch of PHN and in all trigeminal branches of BMS. The salty, sweet and olfactory thresholds were higher in all studied groups; the sour threshold was lower and there were no differences of bitter. CONCLUSION: All groups showed abnormal thresholds of gustation and olfaction; somesthetic findings were discrete in ITN and more common in PHN and BMS; central mechanisms of balance of sensorial inputs might be underlying these observations. Oral Diseases (2010) 16, 482-487

Exercise prevents the effects of experimental arthritis on the metabolism and function of immune cells

Active lymphocytes (LY) and macrophages (M Phi) are involved in the pathophysiology of rheumatoid arthritis (RA) Due to its anti-inflammatory effect. physical exercise may be beneficial in RA by acting on the immune system (IS) Thus, female Wistar rats with type II collagen-induced arthritis (CIA) were submitted to swimming training (6 weeks. 5 days/week. 60 min/day) and some biochemical and immune parameters, such as the metabolism of glucose and glutamine and function of LY and M. were evaluated In addition, plasma levels of some hormones and of interleukin-2 (IL-2) were also determined Results demonstrate that CIA increased lymphocyte proliferation (1.9- and 1 7-fold, respectively, in response to concanavalin A (ConA) and lipopolysaccharide (LPS)), as well as macrophage H(2)O(2) production (1 6-fold), in comparison to control Exercise training prevented the activation of immune cells, induced by CIA. and established a pattern of substrate utilization similar to that described as normal for these cells. Exercise also promoted an elevation of plasma levels of corticosterone (22 2%), progesterone (1 7-fold) and IL-2 (2 6-fold) Our data suggest that chronic exercise is able to counterbalance the effects of CIA on cells of the IS. reinforcing the proposal that the benefits of exercise may not be restricted to aerobic capacity and/or strength improvement Copyright (C) 2010 John Wiley & Sons, Ltd

Analysis of pelvic floor musculature function during pregnancy and postpartum: a cohort study

Aim. To compare the measurements of women`s pelvic floor musculature strength (PFMS) during pregnancy and postpartum period. Background. Pregnancy and childbirth can have an influence on the muscles and pelvic floor and can cause morbidities of women`s genito-urinary tract. Design. A prospective cohort study. Methods. There were included 226 primigravidae women, attended by community health services in the city of Itapecerica da Serra, Sao Paulo, Brazil. The participants were followed in four stages: (1) within 12 weeks of pregnancy; (2) between 36-40 weeks of pregnancy; (3) within 48 hours after childbirth; (4) 42-60 days after childbirth. Data were collected from February 2007-August 2008. The pelvic floor musculature strength was evaluated by perineometry and digital vaginal palpation in stages 1, 2 and 4. The final sample included 110 women who completed all four stages of the study. Results. The pelvic floor musculature strength of the women did not change significantly during pregnancy or after delivery (anova: p = 0 center dot 78). In all three examined stages, a low-intensity pelvic floor musculature strength was prevalent (in mmHg: stage 1 = 15 center dot 9; stage 2 = 15 center dot 2, stage 4 = 14 center dot 7), with scores from 0-3 on the Oxford scale. The pelvic floor musculature strength did not differ in relation to maternal age, skin colour, conjugal status, dyspareunia, stool characteristics, type of delivery, or conditions of the perineum. An interaction between maternal nutritional state and newborn`s weight may affect the pelvic floor musculature strength (manova: p = 0 center dot 04). Conclusion. Pregnancy and childbirth did not reduce significantly pelvic floor musculature strength. The perineometry and digital vaginal palpation used to assess the pelvic floor musculature strength were well accepted by the women. Relevance to clinical practice. In clinical practice, digital vaginal palpation is effective for supporting the diagnosis of urinary, intestinal and sexual dysfunctions. Perineometry use is particularly important together with the performance of perineal exercises with biofeedback in the treatment these disorders.

Aerobic exercise training improves skeletal muscle function and Ca(2+) handling-related protein expression in sympathetic hyperactivity-induced heart failure

Bueno CR Jr, Ferreira JC, Pereira MG, Bacurau AV, Brum PC. Aerobic exercise training improves skeletal muscle function and Ca(2+) handling-related protein expression in sympathetic hyperactivity-induced heart failure. J Appl Physiol 109: 702-709, 2010. First published July 1, 2010; doi: 10.1152/japplphysiol.00281.2010.-The cellular mechanisms of positive effects associated with aerobic exercise training on overall intrinsic skeletal muscle changes in heart failure (HF) remain unclear. We investigated potential Ca(2+) abnormalities in skeletal muscles comprising different fiber compositions and investigated whether aerobic exercise training would improve muscle function in a genetic model of sympathetic hyperactivity-induced HF. A cohort of male 5-mo-old wild-type (WT) and congenic alpha(2A)/alpha(2C) adrenoceptor knockout (ARKO) mice in a C57BL/6J genetic background were randomly assigned into untrained and trained groups. Exercise training consisted of a 8-wk running session of 60 min, 5 days/wk (from 5 to 7 mo of age). After completion of the exercise training protocol, exercise tolerance was determined by graded treadmill exercise test, muscle function test by Rotarod, ambulation and resistance to inclination tests, cardiac function by echocardiography, and Ca(2+) handling-related protein expression by Western blot. alpha(2A)/alpha(2C)ARKO mice displayed decreased ventricular function, exercise intolerance, and muscle weakness paralleled by decreased expression of sarcoplasmic Ca(2+) release-related proteins [alpha(1)-, alpha(2)-, and beta(1)-subunits of dihydropyridine receptor (DHPR) and ryanodine receptor (RyR)] and Ca(2+) reuptake-related proteins [sarco(endo) plasmic reticulum Ca(2+)-ATPase (SERCA) 1/2 and Na(+)/Ca(2+) exchanger (NCX)] in soleus and plantaris. Aerobic exercise training significantly improved exercise tolerance and muscle function and reestablished the expression of proteins involved in sarcoplasmic Ca(2+) handling toward WT levels. We provide evidence that Ca(2+) handling-related protein expression is decreased in this HF model and that exercise training improves skeletal muscle function associated with changes in the net balance of skeletal muscle Ca(2+) handling proteins.