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/.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

Exceptional times for the dynamical discrete web

The dynamical discrete web (DyDW), introduced in the recent work of Howitt and Warren, is a system of coalescing simple symmetric one-dimensional random walks which evolve in an extra continuous dynamical time parameter tau. The evolution is by independent updating of the underlying Bernoulli variables indexed by discrete space-time that define the discrete web at any fixed tau. In this paper, we study the existence of exceptional (random) values of tau where the paths of the web do not behave like usual random walks and the Hausdorff dimension of the set of such exceptional tau. Our results are motivated by those about exceptional times for dynamical percolation in high dimension by Haggstrom, Peres and Steif, and in dimension two by Schramm and Steif. The exceptional behavior of the walks in the DyDW is rather different from the situation for the dynamical random walks of Benjamini, Haggstrom, Peres and Steif. For example, we prove that the walk from the origin S(0)(tau) violates the law of the iterated logarithm (LIL) on a set of tau of Hausdorff dimension one. We also discuss how these and other results should extend to the dynamical Brownian web, the natural scaling limit of the DyDW. (C) 2009 Elsevier B.V. All rights reserved.

A PHASE TRANSITION FOR COMPETITION INTERFACES

We study the competition interface between two growing clusters in a growth model associated to last-passage percolation. When the initial unoccupied set is approximately a cone, we show that this interface has an asymptotic direction with probability 1. The behavior of this direction depends on the angle theta of the cone: for theta >= 180 degrees, the direction is deterministic, while for theta < 180 degrees, it is random, and its distribution can be given explicitly in certain cases. We also obtain partial results on the fluctuations of the interface around its asymptotic direction. The evolution of the competition interface in the growth model can be mapped onto the path of a second-class particle in the totally asymmetric simple exclusion process; from the existence of the limiting direction for the interface, we obtain a new and rather natural proof of the strong law of large numbers (with perhaps a random limit) for the position of the second-class particle at large times.

POLLING SYSTEMS WITH PARAMETER REGENERATION, THE GENERAL CASE

We consider a polling model with multiple stations, each with Poisson arrivals and a queue of infinite capacity. The service regime is exhaustive and there is Jacksonian feedback of served customers. What is new here is that when the server comes to a station it chooses the service rate and the feedback parameters at random; these remain valid during the whole stay of the server at that station. We give criteria for recurrence, transience and existence of the sth moment of the return time to the empty state for this model. This paper generalizes the model, when only two stations accept arriving jobs, which was considered in [Ann. Appl. Probab. 17 (2007) 1447-1473]. Our results are stated in terms of Lyapunov exponents for random matrices. From the recurrence criteria it can be seen that the polling model with parameter regeneration can exhibit the unusual phenomenon of null recurrence over a thick region of parameter space.

Polarity-Dependent Transcranial Direct Current Stimulation Effects on Central Auditory Processing

Given the polarity dependent effects of transcranial direct current stimulation (tDCS) in facilitating or inhibiting neuronal processing, and tDCS effects on pitch perception, we tested the effects of tDCS on temporal aspects of auditory processing. We aimed to change baseline activity of the auditory cortex using tDCS as to modulate temporal aspects of auditory processing in healthy subjects without hearing impairment. Eleven subjects received 2mA bilateral anodal, cathodal and sham tDCS over auditory cortex in a randomized and counterbalanced order. Subjects were evaluated by the Random Gap Detection Test (RGDT), a test measuring temporal processing abilities in the auditory domain, before and during the stimulation. Statistical analysis revealed a significant interaction effect of time vs. tDCS condition for 4000 Hz and for clicks. Post-hoc tests showed significant differences according to stimulation polarity on RGDT performance: anodal improved 22.5% and cathodal decreased 54.5% subjects' performance, as compared to baseline. For clicks, anodal also increased performance in 29.4% when compared to baseline. tDCS presented polarity-dependent effects on the activity of the auditory cortex, which results in a positive or negative impact in a temporal resolution task performance. These results encourage further studies exploring tDCS in central auditory processing disorders.

Placebo Response of Non-Pharmacological and Pharmacological Trials in Major Depression: A Systematic Review and Meta-Analysis

Background: Although meta-analyses have shown that placebo responses are large in Major Depressive Disorder (MDD) trials; the placebo response of devices such as repetitive transcranial magnetic stimulation (rTMS) has not been systematically assessed. We proposed to assess placebo responses in two categories of MDD trials: pharmacological (antidepressant drugs) and non-pharmacological (device-rTMS) trials. Methodology/Principal Findings: We performed a systematic review and meta-analysis of the literature from April 2002 to April 2008, searching MEDLINE, Cochrane, Scielo and CRISP electronic databases and reference lists from retrieved studies and conference abstracts. We used the keywords placebo and depression and escitalopram for pharmacological studies; and transcranial magnetic stimulation and depression and sham for non-pharmacological studies. All randomized, double-blinded, placebo-controlled, parallel articles on major depressive disorder were included. Forty-one studies met our inclusion criteria-29 in the rTMS arm and 12 in the escitalopram arm. We extracted the mean and standard values of depression scores in the placebo group of each study. Then, we calculated the pooled effect size for escitalopram and rTMS arm separately, using Cohen's d as the measure of effect size. We found that placebo response are large for both escitalopram (Cohen's d-random-effects model-1.48; 95% C.I. 1.26 to 1.6) and rTMS studies (0.82; 95% C.I. 0.63 to 1). Exploratory analyses show that sham response is associated with refractoriness and with the use of rTMS as an add-on therapy, but not with age, gender and sham method utilized. Conclusions/Significance: We confirmed that placebo response in MDD is large regardless of the intervention and is associated with depression refractoriness and treatment combination (add-on rTMS studies). The magnitude of the placebo response seems to be related with study population and study design rather than the intervention itself.

Environmental stresses inhibit splicing in the aquatic fungus Blastocladiella emersonii

Background: Exposure of cells to environmental stress conditions can lead to the interruption of several intracellular processes, in particular those performed by macromolecular complexes such as the spliceosome. Results: During nucleotide sequencing of cDNA libraries constructed using RNA isolated from B. emersonii cells submitted to heat shock and cadmium stress, a large number of ESTs with retained introns was observed. Among the 6,350 ESTs obtained through sequencing of stress cDNA libraries, 181 ESTs presented putative introns (2.9%), while sequencing of cDNA libraries from unstressed B. emersonii cells revealed only 0.2% of ESTs containing introns. These data indicate an enrichment of ESTs with introns in B. emersonii stress cDNA libraries. Among the 85 genes corresponding to the ESTs that retained introns, 19 showed more than one intron and three showed three introns, with intron length ranging from 55 to 333 nucleotides. Canonical splicing junctions were observed in most of these introns, junction sequences being very similar to those found in introns from genes previously characterized in B. emersonii, suggesting that inhibition of splicing during stress is apparently a random process. Confirming our observations, analyses of gpx3 and hsp70 mRNAs by Northern blot and S1 protection assays revealed a strong inhibition of intron splicing in cells submitted to cadmium stress. Conclusion: In conclusion, data indicate that environmental stresses, particularly cadmium treatment, inhibit intron processing in B. emersonii, revealing a new adaptive response to cellular exposure to this heavy metal.

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.

Disparity between Multilocus Enzyme Electrophoresis, Microsatellite Markers and Pulsed-Field Gel Electrophoresis in epidemiological tracking of Candida albicans

Various molecular systems are available for epidemiological, genetic, evolutionary, taxonomic and systematic studies of innumerable fungal infections, especially those caused by the opportunistic pathogen C. albicans. A total of 75 independent oral isolates were selected in order to compare Multilocus Enzyme Electrophoresis (MLEE), Electrophoretic Karyotyping (EK) and Microsatellite Markers (Simple Sequence Repeats - SSRs), in their abilities to differentiate and group C. albicans isolates (discriminatory power), and also, to evaluate the concordance and similarity of the groups of strains determined by cluster analysis for each fingerprinting method. Isoenzyme typing was performed using eleven enzyme systems: Adh, Sdh, M1p, Mdh, Idh, Gdh, G6pdh, Asd, Cat, Po, and Lap (data previously published). The EK method consisted of chromosomal DNA separation by pulsed-field gel electrophoresis using a CHEF system. The microsatellite markers were investigated by PCR using three polymorphic loci: EF3, CDC3, and HIS3. Dendrograms were generated by the SAHN method and UPGMA algorithm based on similarity matrices (S(SM)). The discriminatory power of the three methods was over 95%, however a paired analysis among them showed a parity of 19.7-22.4% in the identification of strains. Weak correlation was also observed among the genetic similarity matrices (S(SM)(MLEE) x S(SM)(EK) x S(SM)(SSRs)). Clustering analyses showed a mean of 9 +/- 12.4 isolates per cluster (3.8 +/- 8 isolates/taxon) for MLEE, 6.2 +/- 4.9 isolates per cluster (4 +/- 4.5 isolates/taxon) for SSRs, and 4.1 +/- 2.3 isolates per cluster (2.6 +/- 2.3 isolates/taxon) for EK. A total of 45 (13%), 39(11.2%), 5 (1.4%) and 3 (0.9%) clusters pairs from 347 showed similarity (Si) of 0.1-10%, 10.1-20%, 20.1-30% and 30.1-40%, respectively. Clinical and molecular epidemiological correlation involving the opportunistic pathogen C. albicans may be attributed dependently of each method of genotyping (i.e., MLEE, EK, and SSRs) supplemented with similarity and grouping analysis. Therefore, the use of genotyping systems that give results which offer minimum disparity, or the combination of the results of these systems, can provide greater security and consistency in the determination of strains and their genetic relationships. (C) 2010 Elsevier B.V. All rights reserved.

Characterization of unisexual flower development in the endangered mahogany tree Swietenia macrophylla King. (Meliaceae)

The selection of candidate plus trees of desirable phenotypes from tropical forest trees and the rapid devastation of the natural environments in which these trees are found have created the need for a more detailed knowledge of the floral and reproductive biology of tropical tree species. In this article, the organogenic processes related to unisexual flower development in tropical mahogany, Swietenia macrophylla, are described. Mahogany inflorescences at different developmental stages were evaluated using scanning electron microscopy or optical microscopy of histological sections. The unisexual flowers of S. macrophylla are usually formed in a thyrse, in which the positions of the female and male flowers are not random. Differences between male and female flowers arise late during development. Both female and male flowers can only be structurally distinguished after stage 9, where ovule primordia development is arrested in male flowers and microspore development is aborted in female flower anthers. After this stage, male and female flowers can be distinguished by the naked eye as a result of differences in the dimensions of the gynoecium. The floral characteristics of S. macrophylla (distribution of male and female flowers within the inflorescence, and the relative number of male to female flowers) have practical implications for conservation strategies of this endangered species. (c) 2008 The Linnean Society of London, Botanical Journal of the Linnean Society, 2008, 156, 529-535.

Stretching and folding mechanism in foams

We have described the stretching and folding of foams in a vertical Hele-Shaw cell containing air and a surfactant solution, from a sequence of upside-down flips. Besides the firactal dimension of the foam, we have observed the logistic growth for the soap film length. The stretching and folding mechanism is present during the foam formation, and this mechanism is observed even after the foam has reached its respective maximum fractal dimension. Observing the motion of bubbles inside the foam, large bubbles present power spectrum associated with random walk motion in both directions, while the small bubbles are scattered like balls in a Galton board. (C) 2008 Published by Elsevier B.V.

Enrollment in Physical Education Is Associated With Health-Related Behavior Among High School Students

METHODS: A total of 4210 students attending public high schools in Pernambuco (northeast of Brazil) were selected using random 2-stage cluster sampling. Data were collected by using the Global School-based Student Health Survey. The independent variable was the frequency of participation in PE classes, whereas physical activity, television viewing, smoking, and alcohol, fruit, vegetables and soda consumption were dependent variables. Logistic regressions were carried out to perform crude and adjusted analysis of the association between enrollment in PE classes and health-related behaviors. RESULTS: Sixty-five percent of students do not take part in PE classes, with a significantly higher proportion among females (67.8%). It was observed that enrollment in PE classes was positively associated with physical activity, TV viewing, and fruit consumption, but was negatively associated with soda drinking. The likelihood of reporting being active and eating fruit on a daily basis was 27% and 45% higher, respectively, among those who participate in at least 2 classes per week in comparison with those who do not. Students who participate in PE classes had 28-30% higher likelihood of reporting lower TV viewing during week days. CONCLUSIONS: Findings suggest that higher levels of enrollment in PE classes could play a role in the promotion of health-related behaviors among high school students.

Atenolol blunts blood pressure increase during dynamic resistance exercise in hypertensives

center dot Dynamic resistance exercise promotes a sizeable increase in blood pressure during its execution in non medicated hypertensives. WHAT THIS STUDY ADDS center dot Atenolol not only decreases blood pressure level but also mitigates the increase of blood pressure during dynamic resistance exercise in hypertensive patients. An increase in blood pressure during resistance exercise might be at least in part attributed to an increase in cardiac output. AIMS This study was conducted to determine whether atenolol was able to decrease BP level and mitigate BP increase during dynamic resistance exercise performed at three different intensities in hypertensives. METHODS Ten essential hypertensives (systolic/diastolic BP between 140/90 and 160/105 mmHg) were blindly studied after 6 weeks of placebo and atenolol. In each phase, volunteers executed, in a random order, three protocols of knee-extension exercises to fatigue: (i) one set at 100% of 1 RM; (ii) three sets at 80% of 1 RM; and (iii) three sets at 40% of 1 RM. Intra-arterial radial blood pressure was measured throughout the protocols. RESULTS Atenolol decreased systolic BP maximum values achieved during the three exercise protocols (100% = 186 +/- 4 vs. 215 +/- 7, 80% = 224 +/- 7 vs. 247 +/- 9 and 40% = 223 +/- 7 vs. 252 +/- 16 mmHg, P < 0.05). Atenolol also mitigated an increase in systolic BP in the first set of exercises (100% = +38 +/- 5 vs. +54 +/- 9; 80% = +68 +/- 11 vs. +84 +/- 13 and 40% = +69 +/- 7 vs. +84 +/- 14, mmHg, P < 0.05). Atenolol decreased diastolic BP values and mitigated its increase during exercise performed at 100% of 1 RM (126 +/- 6 vs. 145 +/- 6 and +41 +/- 6 vs. +52 +/- 6, mmHg, P < 0.05), but not at the other exercise intensities. CONCLUSIONS Atenolol was effective in both reducing systolic BP maximum values and mitigating BP increase during resistance exercise performed at different intensities in hypertensive subjects.

The rating of perceived exertion predicts intermittent vertical jump demand and performance

The aims of this study were (a) to assess the ability of the rating of perceived exertion (RPE) to predict performance (i.e. number of vertical jumps performed to a fixed jump height) of an intermittent vertical jump exercise, and (b) to determine the ability of RPE to describe the physiological demand of such exercise. Eight healthy men performed intermittent vertical jumps with rest periods of 4, 5, and 6s until fatigue. Heart rate and RPE were recorded every five jumps throughout the sessions. The number of vertical jumps performed was also recorded. Random coefficient growth curve analysis identified relationships between the number of vertical jumps and both RPE and heart rate for which there were similar slopes. In addition, there were no differences between individual slopes and the mean slope for either RPE or heart rate. Moreover, RPE and number of jumps were highly correlated throughout all sessions (r=0.97-0.99; P0.001), as were RPE and heart rate (r=0.93-0.97; P0.001). The findings suggest that RPE can both predict the performance of intermittent vertical jump exercise and describe the physiological demands of such exercise.