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.

Construct and Compare Gene Coexpression Networks with DAPfinder and DAPview

Background: DAPfinder and DAPview are novel BRB-ArrayTools plug-ins to construct gene coexpression networks and identify significant differences in pairwise gene-gene coexpression between two phenotypes. Results: Each significant difference in gene-gene association represents a Differentially Associated Pair (DAP). Our tools include several choices of filtering methods, gene-gene association metrics, statistical testing methods and multiple comparison adjustments. Network results are easily displayed in Cytoscape. Analyses of glioma experiments and microarray simulations demonstrate the utility of these tools. Conclusions: DAPfinder is a new friendly-user tool for reconstruction and comparison of biological networks.

The impact of measurement errors in the identification of regulatory networks

Background: There are several studies in the literature depicting measurement error in gene expression data and also, several others about regulatory network models. However, only a little fraction describes a combination of measurement error in mathematical regulatory networks and shows how to identify these networks under different rates of noise. Results: This article investigates the effects of measurement error on the estimation of the parameters in regulatory networks. Simulation studies indicate that, in both time series (dependent) and non-time series (independent) data, the measurement error strongly affects the estimated parameters of the regulatory network models, biasing them as predicted by the theory. Moreover, when testing the parameters of the regulatory network models, p-values computed by ignoring the measurement error are not reliable, since the rate of false positives are not controlled under the null hypothesis. In order to overcome these problems, we present an improved version of the Ordinary Least Square estimator in independent (regression models) and dependent (autoregressive models) data when the variables are subject to noises. Moreover, measurement error estimation procedures for microarrays are also described. Simulation results also show that both corrected methods perform better than the standard ones (i.e., ignoring measurement error). The proposed methodologies are illustrated using microarray data from lung cancer patients and mouse liver time series data. Conclusions: Measurement error dangerously affects the identification of regulatory network models, thus, they must be reduced or taken into account in order to avoid erroneous conclusions. This could be one of the reasons for high biological false positive rates identified in actual regulatory network models.

Dynamic light scattering and optical absorption spectroscopy study of pH and temperature stabilities of the extracellular hemoglobin of Glossoscolex paulistus

The extracellular hemoglobin of Glossoscolex paulistus (HbGp) is constituted of subunits containing heme groups, monomers and trimers, and nonheme structures, called linkers, and the whole protein has a minimum molecular mass near 3.1 x 10(6) Da. This and other proteins of the same family are useful model systems for developing blood substitutes due to their extracellular nature, large size, and resistance to oxidation. HbGp samples were studied by dynamic light scattering (DLS). In the pH range 6.0-8.0, HbGp is stable and has a monodisperse size distribution with a z-average hydrodynamic diameter (D-h) of 27 +/- 1 nm. A more alkaline pH induced an irreversible dissociation process, resulting in a smaller D-h of 10 +/- 1 nm. The decrease in D-h suggests a complete hemoglobin dissociation. Gel filtration chromatography was used to show unequivocally the oligomeric dissociation observed at alkaline pH. At pH 9.0, the dissociation kinetics is slow, taking a minimum of 24 h to be completed. Dissociation rate constants progressively increase at higher pH, becoming, at pH 10.5, not detectable by DILS. Protein temperature stability was also pH-dependent. Melting curves for HbGp showed oligomeric dissociation and protein denaturation as a function of pH. Dissociation temperatures were lower at higher pH. Kinetic studies were also performed using ultraviolet-visible absorption at the Soret band. Optical absorption monitors the hemoglobin autoxidation while DLS gives information regarding particle size changes in the process of protein dissociation. Absorption was analyzed at different pH values in the range 9.0-9.8 and at two temperatures, 25 degrees C and 38 degrees C. At 25 degrees C, for pH 9.0 and 9.3, the kinetics monitored by ultraviolet-visible absorption presents a monoexponential behavior, whereas for pH 9.6 and 9.8, a biexponential behavior was observed, consistent with heme heterogeneity at more alkaline pH. The kinetics at 38 degrees C is faster than that at 25 degrees C and is biexponential in the whole pH range. DLS dissociation rates are faster than the autoxidation dissociation rates at 25 degrees C. Autoxiclation and dissociation processes are intimately related, so that oligomeric protein dissociation promotes the increase of autoxidation rate and vice versa. The effect of dissociation is to change the kinetic character of the autoxidation of hemes from monoexponential to biexponential, whereas the reverse change is not as effective. This work shows that DLS can be used to follow, quantitatively and in real time, the kinetics of changes in the oligomerization of biologic complex supramolecular systems. Such information is relevant for the development of mimetic systems to be used as blood substitutes.

Purification and characterization of an antimicrobial peptide produced by Pseudomonas sp strain 4B

An antimicrobial peptide produced by a bacterium isolated from the effluent pond of a bovine abattoir was purified and characterized. The strain was characterized by biochemical profiling and 16S rDNA sequencing as Pseudomonas sp. The antimicrobial peptide was purified by ammonium sulfate precipitation, gel filtration, and ion exchange chromatography. Direct activity on sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was observed. A major band on SDS-PAGE suggested that the antimicrobial peptide has a molecular mass of about 30 kDa. The substance was inhibitory to a broad range of indicator strains, including pathogenic and food spoilage bacteria such as Listeria monocytogenes, Bacillus cereus, Staphylococcus aureus, among other. The partially purified antimicrobial substance remained active over a wide temperature range and was resistant to all proteases tested. This substance showed different properties than other antimicrobials from Pseudomonas species, suggesting a novel antimicrobial peptide was characterized.

Highly Toxic Microcystis aeruginosa Strain, Isolated from So Paulo-Brazil, Produce Hepatotoxins and Paralytic Shellfish Poison Neurotoxins

While evaluating several laboratory-cultured cyanobacteria strains for the presence of paralytic shellfish poison neurotoxins, the hydrophilic extract of Microcystis aeruginosa strain SPC777-isolated from Billings`s reservoir, So Paulo, Brazil-was found to exhibit lethal neurotoxic effect in mouse bioassay. The in vivo test showed symptoms that unambiguously were those produced by PSP. In order to identify the presence of neurotoxins, cells were lyophilized, and the extracts were analyzed by HPLC-FLD and HPLC-MS. HPLC-FLD analysis revealed four main Gonyautoxins: GTX4(47.6%), GTX2(29.5%), GTX1(21.9%), and GTX3(1.0%). HPLC-MS analysis, on other hand, confirmed both epimers, with positive Zwitterions M(+) 395.9 m/z for GTX3/GTX2 and M(+) 411 m/z for GTX4/GTX1 epimers. The hepatotoxins (Microcystins) were also evaluated by ELISA and HPLC-MS analyses. Positive immunoreaction was observed by ELISA assay. Alongside, the HPLC-MS analyses revealed the presence of [l-ser(7)] MCYST-RR. The N-methyltransferase (NMT) domain of the microcystin synthetase gene mcyA was chosen as the target sequence to detect the presence of the mcy gene cluster. PCR amplification of the NMT domain, using the genomic DNA of the SPC777 strain and the MSF/MSR primer set, resulted in the expected 1,369 bp product. The phylogenetic analyses grouped the NMT sequence with the NMT sequences of other known Microcystis with high bootstrap support. The taxonomical position of M. aeruginosa SPC777 was confirmed by a detailed morphological description and a phylogenetic analysis of 16S rRNA gene sequence. Therefore, co-production of PSP neurotoxins and microcystins by an isolated M. aeruginosa strain is hereby reported for the first time.

A new monitor for routine thermal and epithermal neutron fluence rate monitoring in k(0) INAA

The Zr-Au set for monitoring the thermal and epithermal neutron fluence rate and the epithermal spectrum parameter a is not always practicable for routine application of INAA in well-thermalized facilities. An alternative set consisting of Cr, Au and Mo provides values for the thermal neutron fluence rate, f and alpha that are not significantly different from those found via the Zr-Au method and the Cd-covered Zr-method. The IRMM standard SMELS-II was analyzed using the (Au-Cr-Mo) monitor and a good agreement was obtained. (C) 2008 Elsevier Ltd. All rights reserved.

Collimated scanning LS-INAA for testing trace elements homogeneity in Brazilian coffee beans

The degree of homogeneity is normally assessed by the variability of the results of independent analyses of several (e.g., 15) normal-scale replicates. Large sample instrumental neutron activation analysis (LS-INAA) with a collimated Ge detector allows inspecting the degree of homogeneity of the initial batch material, using a kilogram-size sample. The test is based on the spatial distributions of induced radioactivity. Such test was applied to samples of Brazilian whole (green) coffee beans (Coffea arabica and Coffea canephora) of approximately I kg in the frame of development of a coffee reference material. Results indicated that the material do not contain significant element composition inhomogeneities between batches of approximately 30-50 g, masses typically forming the starting base of a reference material.

Fermentation of cacao (Theobroma cacao L.) seeds with a hybrid Kluyveromyces marxianus strain improved product quality attributes

Fermentation of Theobroma cacao (cacao) seeds is an absolute requirement for the full development of chocolate flavor precursors. An adequate aeration of the fermenting cacao seed mass is a fundamental prerequisite for a satisfactory fermentation. Here, we evaluated whether a controlled inoculation of cacao seed fermentation using a Kluyveromyces marxianus hybrid yeast strain, with an increased pectinolytic activity, would improve an earlier liquid drainage (`sweatings`) from the fermentation mass, developing a superior final product quality. Inoculation with K. marxianus increased by one third the volume of drained liquid and affected the microorganism population structure during fermentation, which was detectable up to the end of the process. Introduction of the hybrid yeast affected the profile of total seed protein degradation evaluated by polyacrylamide gel electrophoresis, with improved seed protein degradation, and reduction of titrable acidity. Sensorial evaluation of the chocolate obtained from beans fermented with the K. marxianus inoculation was more accepted by analysts in comparison with the one from cocoa obtained through natural fermentation. The increase in mass aeration during the first 24 h seemed to be fundamental for the improvement of fermentation quality, demonstrating the potential application of this improved hybrid yeast strain with superior exogenous pectinolytic activity.

Adjustment of Mineral Elements in the Culture Medium for the Micropropagation of Three Vriesea Bromeliads from the Brazilian Atlantic Forest: The Importance of Calcium

Many different species of Bromeliaceae are endangered and their conservation requires specific knowledge of their growth habits and propagation. In vitro culture of bromeliads is an important method for efficient clonal propagation and ill vitro seed g,germination can be used to maintain genetic variability. The present work aims to evaluate the in vitro growth and nutrient concentration in leaves of the epiphyte bromeliads Vriesea friburguensis Mez, Vriesea hieroglyphica (Carriere) E. Morren, and Vriesea unilateralis Mez, which exhibit slow rates of growth in vivo and in vitro. Initially, we compared the endogenous mineral composition of bromeliad plantlets grown in half-strength Murashige and Skoog (MS) medium and the mineral composition considered adequate in the literature. This approach suggested that calcium (Ca) is a critical nutrient and this was considered for new media formulation. Three new culture media were defined in which the main changes to half-strength MS medium were an increase in Ca, magnesium, sulfur, copper, and chloride and a decrease in iron, maintaining the nitrate: ammonium rate at approximate to 2:1. The main difference among the three new media formulated was Ca concentration, which varied from 1.5 mm in half-strength MS to 3.0, 6.0, and 12 mm in M2, M3, and M4 media, respectively. Consistently, all three species exhibited significantly higher fresh and dry weight on M4, the newly defined medium with the highest level of Ca (12 mm). Leaf nitrogen, potassium, zinc, magnesium and boron concentrations increased as Ca concentration in the medium increased from 1.5 to 12 mm.

Zone trapping/merging zones in flow analysis: A novel approach for rapid assays involving relatively slow chemical reactions

A novel strategy for accomplishing zone trapping in flow analysis is proposed. The sample and the reagent solutions are simultaneously inserted into convergent carrier streams and the established zones merge together before reaching the detector, where the most concentrated portion of the entire sample zone is trapped. The main characteristics, potentialities and limitations of the strategy were critically evaluated in relation to an analogous flow system with zone stopping. When applied to the spectrophotometric determination of nitrite in river waters, the main figures of merit were maintained, exception made for the sampling frequency which was calculated as 189h(-1), about 32% higher relatively to the analogous system with zone stopping. The sample inserted volume can be increased up to 1.0 mL without affecting sampling frequency and no problems with pump heating or malfunctions were noted after 8-h operation of the system. In contrast to zone stopping, only a small portion of the sample zone is halted with zone trapping, leading to these beneficial effects. (C) 2011 Elsevier B.V. All rights reserved.

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.

Influence of high- and low-carbohydrate diet following glycogen-depleting exercise on heart rate variability and plasma catecholamines

The purpose of this study was to investigate the effects of a short-term low-or high-carbohydrate (CHO) diet consumed after exercise on sympathetic nervous system activity. Twelve healthy males underwent a progressive incremental test; a control measurement of plasma catecholamines and heart rate variability (HRV); an exercise protocol to reduce endogenous CHO stores; a low-or high-CHO diet (counterbalanced order) consumed for 2 days, beginning immediately after the exercise protocol; and a second resting plasma catecholamine and HRV measurement. The exercise and diet protocols and the second round of measurements were performed again after a 1-week washout period. The mean (+/- SD) values of the standard deviation of R-R intervals were similar between conditions (control, 899.0 +/- 146.1 ms; low-CHO diet, 876.8 +/- 115.8 ms; and high-CHO diet, 878.7 +/- 127.7 ms). The absolute high-and low-frequency (HF and LF, respectively) densities of the HRV power spectrum were also not different between conditions. However, normalized HF and LF (i.e., relative to the total power spectrum) were lower and higher, respectively, in the low-CHO diet than in the control diet (mean +/- SD, 17 +/- 9 normalized units (NU) and 83 +/- 9 NU vs. 27 +/- 11 NU and 73 +/- 17 NU, respectively; p < 0.05). The LF/HF ratio was higher with the low-CHO diet than with the control diet (mean +/- SD, 7.2 +/- 6.2 and 4.2 +/- 3.2, respectively; p < 0.05). The mean values of plasma catecholamines were not different between diets. These results suggest that the autonomic control of the heart rate was modified after a short-term low-CHO diet, but plasma catecholamine levels were not altered.

Peak torque and rate of torque development in elderly with and without fall history

Background Falls are one of the greatest concerns among the elderly A number of studies have described peak torque as one of the best fall-related predictor. No studies have comprehensively focused on the rate of torque development of the lower limb muscles among elderly fallers. Then, the aim of this study was to determine the relationship between muscle peak torque and rate of torque development of the lower limb joints in elderly with and without fall history It was also aimed to determine whether these parameters of muscle performance (i e, peak torque and rate of torque development) are related to the number of falls. Methods: Thirty-one women volunteered to participate in the study and were assigned in one of the groups according to the number of falls over the 12 months that preceded the present Then, participants with no fall history (Cl; n = 13; 67.6[7.5] years-old), one fall (GII; n = 8, 66 0[4 91 years-old) and two or more falls (GIII, n = 10; 67.8[8.8] years-old) performed a number of lower limb maximal isometric voluntary contractions from which peak torque and rate of torque development were quantified Findings. Primary outcomes indicated no peak torque differences between experimental groups in any lower limb joint. The rate of torque development of the knee flexor muscles observed in the non-fallers (Cl) was greater than that observed in the fallers (P < 0.05) and had a significant relationship with the number of falls (P < 0 05) Interpretation. The greater knee flexor muscles` rate of torque development found in the non-fallers in comparison to the fallers indicated that the ability of the elderly to rapidly reorganise the arrangement of the lower limb may play a significant role in allowing the elderly to recover balance after a trip. Thus, training stimulus aimed to improve the rate of torque development may be more beneficial to prevent falls among the elderly than other training stimulus, which are not specifically designed to improve the ability to rapidly produce large amounts of torque (C) 2010 Published by Elsevier Ltd