Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

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

TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

Genome-Wide Detection of Serpentine Receptor-Like Proteins in Malaria Parasites

Serpentine receptors comprise a large family of membrane receptors distributed over diverse organisms, such as bacteria, fungi, plants and all metazoans. However, the presence of serpentine receptors in protozoan parasites is largely unknown so far. In the present study we performed a genome-wide search for proteins containing seven transmembrane domains (7TM) in the human malaria parasite Plasmodium falciparum and identified four serpentine receptor-like proteins. These proteins, denoted PfSR1, PfSR10, PfSR12 and PfSR25, show membrane topologies that resemble those exhibited by members belonging to different families of serpentine receptors. Expression of the pfsrs genes was detected by Real Time PCR in P. falciparum intraerythrocytic stages, indicating that they potentially code for functional proteins. We also found corresponding homologues for the PfSRs in five other Plasmodium species, two primate and three rodent parasites. PfSR10 and 25 are the most conserved receptors among the different species, while PfSR1 and 12 are more divergent. Interestingly, we found that PfSR10 and PfSR12 possess similarity to orphan serpentine receptors of other organisms. The identification of potential parasite membrane receptors raises a new perspective for essential aspects of malaria parasite host cell infection.

A Morphological and Histological Characterization of Bisexual and Male Flower Types in Pomegranate

Pomegranate [Punica granatum (Punicaceae)] is characterized by having two types of flowers on the same tree: hermaphroditic bisexual flowers and functionally male flowers. This condition, defined as functional andromonoecy, can result in decreased yields resulting from the inability of male flowers to set fruit. Morphological and histological analyses of bisexual and male flowers were conducted using light and scanning electron microscopy (SEM) to characterize the different flower types observed in pomegranate plants and to better understand their developmental differences. Bisexual flowers had a discoid stigma covered with copious exudate, elongated stigmatic papillae, a single elongate style, and numerous stamens inserted on the inner wall of the calyx tube. Using fluorescence staining, high numbers of pollen tubes were observed growing through a central stylar canal. Ovules were numerous, elliptical, and anatropous. In contrast, male flowers had reduced female parts and exhibited shortened pistils of variable heights. Stigmatic papillae of male flowers had little exudate yet supported pollen germination. However, pollen tubes were rarely observed in styles. Ovules in male flowers were rudimentary and exhibited various stages of degeneration. Pollen from both types of flowers was of similar size, approximate to 20 mu m, and exhibited similar percent germination using in vitro germination assays. Pollen germination was strongly influenced by temperature. Maximal germination (greater than 74%) was obtained at 25 and 35 degrees C; pollen germination was significantly lower at 15 degrees C (58%) and 5 degrees C (10%).

Sequential injections as an alternative to gradient exploitation for implementing differential kinetic analysis in a flow injection system

A novel flow-based strategy for implementing simultaneous determinations of different chemical species reacting with the same reagent(s) at different rates is proposed and applied to the spectrophotometric catalytic determination of iron and vanadium in Fe-V alloys. The method relies on the influence of Fe(II) and V(IV) on the rate of the iodide oxidation by Cr(VI) under acidic conditions, the Jones reducing agent is then needed Three different plugs of the sample are sequentially inserted into an acidic KI reagent carrier stream, and a confluent Cr(VI) solution is added downstream Overlap between the inserted plugs leads to a complex sample zone with several regions of maximal and minimal absorbance values. Measurements performed on these regions reveal the different degrees of reaction development and tend to be more precise Data are treated by multivariate calibration involving the PLS algorithm The proposed system is very simple and rugged Two latent variables carried out ca 95% of the analytical information and the results are in agreement with ICP-OES. (C) 2010 Elsevier B V. All rights reserved.

Analysis of Single-Nucleotide Polymorphisms in the crt-o and mdr1 Genes of Plasmodium vivax among Chloroquine-Resistant Isolates from the Brazilian Amazon Region

Plasmodium vivax parasites with chloroquine resistance (CQR) are already circulating in the Brazilian Amazon. Complete single-nucleotide polymorphism (SNP) analyses of coding and noncoding sequences of the pvmdr1 and pvcrt-o genes revealed no associations with CQR, even if some mutations had not been randomly selected. In addition, striking differences in the topologies and numbers of SNPs in these transporter genes between P. vivax and P. falciparum reinforce the idea that mechanisms other than mutations may explain this virulent phenotype in P. vivax.

Skewness, splitting number and vertex deletion of some toroidal meshes

The skewness sk(G) of a graph G = (V, E) is the smallest integer sk(G) >= 0 such that a planar graph can be obtained from G by the removal of sk(C) edges. The splitting number sp(G) of C is the smallest integer sp(G) >= 0 such that a planar graph can be obtained from G by sp(G) vertex splitting operations. The vertex deletion vd(G) of G is the smallest integer vd(G) >= 0 such that a planar graph can be obtained from G by the removal of vd(G) vertices. Regular toroidal meshes are popular topologies for the connection networks of SIMD parallel machines. The best known of these meshes is the rectangular toroidal mesh C(m) x C(n) for which is known the skewness, the splitting number and the vertex deletion. In this work we consider two related families: a triangulation Tc(m) x c(n) of C(m) x C(n) in the torus, and an hexagonal mesh Hc(m) x c(n), the dual of Tc(m) x c(n) in the torus. It is established that sp(Tc(m) x c(n)) = vd(Tc(m) x c(n) = sk(Hc(m) x c(n)) = sp(Hc(m) x c(n)) = vd(Hc(m) x c(n)) = min{m, n} and that sk(Tc(m) x c(n)) = 2 min {m, n}.

Treatment of Obese Adolescents: The Influence of Periodization Models and ACE Genotype

The aims of the present study were to compare the effects of two periodization models on metabolic syndrome risk factors in obese adolescents and verify whether the angiotensin-converting enzyme (ACE) genotype is important in establishing these effects. A total of 32 postpuberty obese adolescents were submitted to aerobic training (AT) and resistance training (RT) for 14 weeks. The subjects were divided into linear periodization (LP, n = 16) or daily undulating periodization (DUP, n = 16). Body composition, visceral and subcutaneous fat, glycemia, insulinemia, homeostasis model assessment of insulin resistance (HOMA-IR), lipid profiles, blood pressure, maximal oxygen consumption (VO(2max)), resting metabolic rate (RMR), muscular endurance were analyzed at baseline and after intervention. Both groups demonstrated a significant reduction in body mass, BMI, body fat, visceral and subcutaneous fat, total and low-density lipoprotein cholesterol, blood pressure and an increase in fat-free mass, VO(2max), and muscular endurance. However, only DUP promoted a reduction in insulin concentrations and HOMA-IR. It is important to emphasize that there was no statics difference between LP and DUP groups; however, it appears that there may be bigger changes in the DUP than LP group in some of the metabolic syndrome risk factors in obese adolescents with regard to the effect size (ES). Both periodization models presented a large effect on muscular endurance. Despite the limitation of sample size, our results suggested that the ACE genotype may influence the functional and metabolic characteristics of obese adolescents and may be considered in the future strategies for massive obesity control.

Cardiac and peripheral adjustments induced by early exercise training intervention were associated with autonomic improvement in infarcted rats: role in functional capacity and mortality

Aims To test the effects of early exercise training (ET) on left ventricular (LV) and autonomic functions, haemodynamics, tissues blood flows (BFs), maximal oxygen consumption (VO(2) max), and mortality after myocardial infarction (MI) in rats. Methods and results Male Wistar rats were divided into: control (C), sedentary-infarcted (SI), and trained-infarcted (TI). One week after MI, TI group underwent an ET protocol (90 days, 50-70% VO2 max). Left ventricular function was evaluated noninvasively and invasively. Baroreflex sensitivity, heart rate variability, and pulse interval were measured. Cardiac output (CO) and regional BFs were determined using coloured microspheres. Infarcted area was reduced in TI (19 +/- 6%) compared with SI (34 +/- 5%) after ET. Exercise training improved the LV and autonomic functions, the CO and regional BF changes induced by MI, as well as increased SERCA2 expression and mRNA vascular endothelial growth factor levels. These changes brought about by ET resulted in mortality rate reduction in the TI (13%) group compared with the SI (54%) group. Conclusion Early aerobic ET reduced cardiac and peripheral dysfunctions and preserved cardiovascular autonomic control after MI in trained rats. Consequently, these ET-induced changes resulted in improved functional capacity and survival after MI.

Exercise intensity optimization for men with high cardiorespiratory fitness

Exercise intensity is a key parameter for exercise prescription but the optimal range for individuals with high cardiorespiratory fitness is unknown. The aims of this study were (1) to determine optimal heart rate ranges for men with high cardiorespiratory fitness based on percentages of maximal oxygen consumption (%VO(2max)) and reserve oxygen consumption (%VO(2reserve)) corresponding to the ventilatory threshold and respiratory compensation point, and ( 2) to verify the effect of advancing age on the exercise intensities. Maximal cardiorespiratory testing was performed on 210 trained men. Linear regression equations were calculated using paired data points between percentage of maximal heart rate (%HR(max)) and %VO(2max) and between percentage of heart rate reserve (%HRR) and %VO(2reserve) attained at each minute during the test. Values of %VO(2max) and %VO(2reserve) at the ventilatory threshold and respiratory compensation point were used to calculate the corresponding values of %HRmax and %HRR, respectively. The ranges of exercise intensity in relation to the ventilatory threshold and respiratory compensation point were achieved at 78-93% of HR(max) and 70-93% of HRR, respectively. Although absolute heart rate decreased with advancing age, there were no age-related differences in %HR(max) and %HRR at the ventilatory thresholds. Thus, in men with high cardiorespiratory fitness, the ranges of exercise intensity based on %HR(max) and %HRR regarding ventilatory threshold were 78-93% and 70-93% respectively, and were not influenced by advancing age.

Exercise training associated with diet improves heart rate recovery and cardiac autonomic nervous system activity in obese children

The purpose of this study was to test the hypotheses that in obese children: 1) hypocaloric diet (D) improves both heart rate recovery at 1 min (Delta HRR1) cfter an exercise test, and cardiac autonomic nervous system activity (CANSA) in obese children; 2) Diet and exercise training (DET) combined leads to greater improvement in both Delta HRR1 after an exercise test and in CANSA, than D alone. Moreover, we examined the relationships among Delta HRR1, CANSA, cardiorespiratory fitness and anthropometric variables (AV) in obese children submitted to D and to DET. 33 obese children (10 +/- 0.2 years; body mass index (BMI) >95(th) percentile) were divided into 2 groups: D (n = 15; BMI = 31 +/- 1 kg/m(2)) and DET (n = 18; 29 +/- 1 kg/m(2)). All children performed a maximal cardiopulmonary exercise test on a treadmill. The Delta HRR1 was defined as the difference between heart rate at peak and at 1-min post-exercise. CANSA was assessed using power spectral analysis of heart rate variability at rest. The sympathovagal balance (low frequency and high frequency ratio, LF/HF) was measured. After interventions, all obese children showed reduced body weight (P < 0.05). The D group did not improve in terms of peak VO(2), Delta HRR1 or LF/HF ratio (P > 0.05). In contrast, the DET group showed increased peak VO(2) (P = 0.01) and improved Delta HRR1 (Delta HRR1 = 37.3 +/- 2.6; P = 0.01) and LF/HF ratio (P = 0.001). The DET group demonstrated significant relationships among Delta HRR1, peak VO(2) and CANSA (P < 0.05). In conclusion, DET, in contrast to D, promoted improved Delta HRR1 and CANSA in obese children, suggesting a positive influence of increased levels of cardiorespiratory fitness by exercise training on cardiac autonomic activity.

Heat and exercise acclimation increases intracellular levels of Hsp72 and inhibits exercise-induced increase in intracellular and plasma Hsp72 in humans

In order to verify the effects of heat and exercise acclimation (HA) on resting and exercise-induced expression of plasma and leukocyte heat shock protein 72 (Hsp72) in humans, nine healthy young male volunteers (25.0 +/- 0.7 years; 80.5 +/- 2.0 kg; 180 +/- 2 cm, mean +/- SE) exercised for 60 min in a hot, dry environment (40 +/- 0A degrees C and 45 A +/- 0% relative humidity) for 11 days. The protocol consisted of running on a treadmill using a controlled hyperthermia technique in which the work rate was adjusted to elevate the rectal temperature by 1A degrees C in 30 min and maintain it elevated for another 30 min. Before and after the HA, the volunteers performed a heat stress test (HST) at 50% of their individual maximal power output for 90 min in the same environment. Blood was drawn before (REST), immediately after (POST) and 1 h after (1 h POST) HST, and plasma and leukocytes were separated and stored. Subjects showed expected adaptations to HA: reduced exercise rectal and mean skin temperatures and heart rate, and augmented sweat rate and exercise tolerance. In HST1, plasma Hsp72 increased from REST to POST and then returned to resting values 1 h POST (REST: 1.11 A +/- 0.07, POST: 1.48 A +/- 0.10, 1 h POST: 1.22 A +/- 0.11 ng mL(-1); p < 0.05). In HST2, there was no change in plasma Hsp72 (REST: 0.94 A +/- 0.08, POST: 1.20 A +/- 0.15, 1 h POST: 1.17 A +/- 0.16 ng mL(-1); p > 0.05). HA increased resting levels of intracellular Hsp72 (HST1: 1 A +/- 0.02 and HST2: 4.2 A +/- 1.2 density units, p < 0.05). Exercise-induced increased intracellular Hsp72 expression was observed on HST1 (HST1: REST, 1 A +/- 0.02 vs. POST, 2.9 A +/- 0.9 density units, mean +/- SE, p < 0.05) but was inhibited on HST2 (HST2: REST, 4.2 +/- 1.2 vs. POST, 4.4 +/- 1.1 density units, p > 0.05). Regression analysis showed that the lower the pre-exercise expression of intracellular Hsp72, the higher the exercise-induced increase (R = -0.85, p < 0.05). In conclusion, HA increased resting leukocyte Hsp72 levels and inhibited exercise-induced expression. This intracellular adaptation probably induces thermotolerance. In addition, the non-increase in plasma Hsp72 after HA may be related to lower stress at the cellular level in the acclimated individuals.

Finger blood pressure during leg resistance exercise

Blood pressure (BP) assessment during resistance exercise can be useful to avoid high BP, reducing cardiovascular risk, especially in hypertensive individuals. However, non-invasive accurate technique for this purpose is not available. The aim of this study was to compare finger photoplethysmographic (FPP) and intra-arterial BP values and responses during resistance exercise. Eight non-medicated hypertensive subjects (5 males, 30-60 years) were evaluated during pre-exercise resting period and during three sets of the knee extension exercise performed at 80% of 1RM until fatigue. BP was measured simultaneously by FPP and intra-arterial methods. Data are mean +/- SD. Systolic BP was significantly higher with FPP than with intra-arterial: at pre-exercise (157 +/- 13 vs. 152 +/- 10 mmHg; p < 0.01) and the mean (202 +/- 29 vs. 198 +/- 26 mmHg; p < 0.01), and the maximal (240 +/- 26 vs. 234 +/- 16 mmHg; p < 0.05) values achieved during exercise. The increase in systolic BP during resistance exercise was similar between FPP and intra-arterial (+ 73 +/- 29 vs. + 71 +/- 18 mmHg; p = 0.59). Diastolic BP values and increases were lower with FPP. In conclusion, FPP provides similar values of BP increment during resistance exercise than intra-arterial method. However, it overestimates by 2.6 +/- 6.1% the maximal systolic BP achieved during this mode of exercise and underestimates by 8.8 +/- 5.8% the maximal diastolic BP.

Role of beta-alanine supplementation on muscle carnosine and exercise performance

ARTIOLI, G. G., B. GUALANO, A. SMITH, J. STOUT, and A. H. LANCHA, JR. Role of beta-Alanine Supplementation on Muscle Carnosine and Exercise Performance. Med. Sci. Sports Exerc., Vol. 42, No. 6, pp. 1162-1173, 2010. In this narrative review, we present and discuss the current knowledge available on carnosine and beta-alanine metabolism as well as the effects of beta-alanine supplementation on exercise performance. Intramuscular acidosis has been attributed to be one of the main causes of fatigue during intense exercise. Carnosine has been shown to play a significant role in muscle pH regulation. Carnosine is synthesized in skeletal muscle from the amino acids L-histidine and beta-alanine. The rate-limiting factor of carnosine synthesis is beta-alanine availability. Supplementation with beta-alanine has been shown to increase muscle carnosine content and therefore total muscle buffer capacity, with the potential to elicit improvements in physical performance during high-intensity exercise. Studies on beta-alanine supplementation and exercise performance have demonstrated improvements in performance during multiple bouts of high-intensity exercise and in single bouts of exercise lasting more than 60 s. Similarly, beta-alanine supplementation has been shown to delay the onset of neuromuscular fatigue. Although beta-alanine does not improve maximal strength or (V) over dotO(2max), some aspects of endurance performance, such as anaerobic threshold and time to exhaustion, can be enhanced. Symptoms of paresthesia may be observed if a single dose higher than 800 mg is ingested. The symptoms, however, are transient and related to the increase in plasma concentration. They can be prevented by using controlled release capsules and smaller dosing strategies. No important side effect was related to the use of this amino acid so far. In conclusion, beta-alanine supplementation seems to be a safe nutritional strategy capable of improving high-intensity anaerobic performance.