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.

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.

Multiclass Hammersley-Aldous-Diaconis process and multiclass-customer queues

In the Hammersley-Aldous-Diaconis process, infinitely many particles sit in R and at most one particle is allowed at each position. A particle at x, whose nearest neighbor to the right is at y, jumps at rate y - x to a position uniformly distributed in the interval (x, y). The basic coupling between trajectories with different initial configuration induces a process with different classes of particles. We show that the invariant measures for the two-class process can be obtained as follows. First, a stationary M/M/1 queue is constructed as a function of two homogeneous Poisson processes, the arrivals with rate, and the (attempted) services with rate rho > lambda Then put first class particles at the instants of departures (effective services) and second class particles at the instants of unused services. The procedure is generalized for the n-class case by using n - 1 queues in tandem with n - 1 priority types of customers. A multi-line process is introduced; it consists of a coupling (different from Liggett's basic coupling), having as invariant measure the product of Poisson processes. The definition of the multi-line process involves the dual points of the space-time Poisson process used in the graphical construction of the reversed process. The coupled process is a transformation of the multi-line process and its invariant measure is the transformation described above of the product measure.

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp aurantifolii

Background: Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C. Results: We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein. Conclusion: We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Contrasting photoelectrochemical behaviour of two isomeric supramolecular dyes based on meso-tetra(pyridyl)porphyrin incorporating four (mu(3)-oxo)- triruthenium(III) clusters

A saddle shaped tetracluster porphyrin species containing four [Ru(3)O(OAc)(6)(py)(2)](+) clusters coordinated to the N-pyridyl atoms of 5,10,15,20-tetra(3-pyridyl)porphyrin, H(2)(3-TCPyP), has been investigated in comparison with the planar tetra(4-pyridyl) porphyrin analogue H(2)(4-TCPyP). The steric effects from the bulky peripheral complexes play a critical role in the H(2)(3-TCPyP) species, determining a non-planar configuration around the porphyrin centre and precluding any significant pi-electronic coupling, in contrast with the less hindered H(2)(4-TCPyP) species. Both systems exhibit a photoelectrochemical response in the presence of nanocrystalline TiO(2) films, involving the porphyrin excitation around 450 nm. However, only in the H(2)(4-TCPyP) case do the cluster moieties also contribute to the photoinduced electron injection process at 670 nm, reflecting the relevance of the electronic coupling between the porphyrin centre and the peripheral complexes.

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.

Soil Carbon Turnover Measurement by Physical Fractionation at a Forest-to-Pasture Chronosequence in the Brazilian Amazon

The effect of conversion from forest-to-pasture upon soil carbon stocks has been intensively discussed, but few studies focus on how this land-use change affects carbon (C) distribution across soil fractions in the Amazon basin. We investigated this in the 20 cm depth along a chronosequence of sites from native forest to three successively older pastures. We performed a physicochemical fractionation of bulk soil samples to better understand the mechanisms by which soil C is stabilized and evaluate the contribution of each C fraction to total soil C. Additionally, we used a two-pool model to estimate the mean residence time (MRT) for the slow and active pool C in each fraction. Soil C increased with conversion from forest-to-pasture in the particulate organic matter (> 250 mu m), microaggregate (53-250 mu m), and d-clay (< 2 mu m) fractions. The microaggregate comprised the highest soil C content after the conversion from forest-to-pasture. The C content of the d-silt fraction decreased with time since conversion to pasture. Forest-derived C remained in all fractions with the highest concentration in the finest fractions, with the largest proportion of forest-derived soil C associated with clay minerals. Results from this work indicate that microaggregate formation is sensitive to changes in management and might serve as an indicator for management-induced soil carbon changes, and the soil C changes in the fractions are dependent on soil texture.

Dynamics of carbon, biomass, and structure in two Amazonian forests

Amazon forests are potentially globally significant sources or sinks for atmospheric carbon dioxide. In this study, we characterize the spatial trends in carbon storage and fluxes in both live and dead biomass (necromass) in two Amazonian forests, the Biological Dynamic of Forest Fragments Project (BDFFP), near Manaus, Amazonas, and the Tapajos National Forest (TNF) near Santarem, Para. We assessed coarse woody debris (CWD) stocks, tree growth, mortality, and recruitment in ground-based plots distributed across the terra firme forest at both sites. Carbon dynamics were similar within each site, but differed significantly between the sites. The BDFFP and the TNF held comparable live biomass (167 +/- 7.6 MgC.ha(-1) versus 149 +/- 6.0 MgC.ha(-1), respectively), but stocks of CWD were 2.5 times larger at TNF (16.2 +/- 1.5 MgC.ha(-1) at BDFFP, versus 40.1 +/- 3.9 MgC.ha(-1) at TNF). A model of current forest dynamics suggests that the BDFFP was close to carbon balance, and its size class structure approximated a steady state. The TNF, by contrast, showed rapid carbon accrual to live biomass (3.24 +/- 0.22 MgC.ha(-1).a(-1) in TNF, 2.59 +/- 0.16 MgC.ha(-1).a(-1) in BDFFP), which was more than offset by losses from large stocks of CWD, as well as ongoing shifts of biomass among size classes. This pattern in the TNF suggests recovery from a significant disturbance. The net loss of carbon from the TNF will likely last 10 - 15 years after the initial disturbance (controlled by the rate of decay of coarse woody debris), followed by uptake of carbon as the forest size class structure and composition continue to shift. The frequency and longevity of forests showing such disequilibruim dynamics within the larger matrix of the Amazon remains an essential question to understanding Amazonian carbon balance.

Functional specialization of Eucalyptus fine roots: contrasting potential uptake rates for nitrogen, potassium and calcium tracers at varying soil depths

1. Little is known about the role of deep roots in the nutrition of forest trees and their ability to provide a safety-net service taking up nutrients leached from the topsoil. 2. To address this issue, we studied the potential uptake of N, K and Ca by Eucalyptus grandis trees (6 years of age - 25 m mean height), in Brazil, as a function of soil depth, texture and water content. We injected NO(3)(-)- (15)N, Rb(+) (analogue of K(+)) and Sr(2+) (analogue of Ca(2+)) tracers simultaneously in a solution through plastic tubes at 10, 50, 150 and 300 cm in depth in a sandy and a clayey Ferralsol soil. A complete randomized design was set up with three replicates of paired trees per injection depth and soil type. Recently expanded leaves were sampled at various times after tracer injection in the summer, and the experiment was repeated in the winter. Soil water contents were continuously monitored at the different depths in the two soils. 3. Determination of foliar Rb and Sr concentrations and (15)N atom % made it possible to estimate the relative uptake potential (RUP) of tracer injections from the four soil depths and the specific RUP (SRUP), defined as RUP, per unit of fine root length density in the corresponding soil layer. 4. The highest tracer uptake rates were found in the topsoil, but contrasting RUP distributions were observed for the three tracers. Whilst the RUP was higher for NO(3)(-)- (15)N than for Rb(+) and Sr(2+) in the upper 50 cm of soil, the highest SRUP values for Sr(2+) and Rb(+) were found at a depth of 300 cm in the sandy soil, as well as in the clayey soil when gravitational solutions reached that depth. 5. Our results suggest that the fine roots of E. grandis trees exhibit contrasting potential uptake rates with depth depending on the nutrient. This functional specialization of roots might contribute to the high growth rates of E. grandis trees, efficiently providing the large amounts of nutrients required throughout the development of these fast-growing plantations.

Radiossensitivity of two propagules of citrus

Studies have shown that the radiossensitivity in plants varies depending on the varieties and the propagules used in mutagenic treatment. The purpose of this study was to evaluate the radiossensitivity of two types of propagules (buds and in vitro segments of epicotyl) in `Murcott` tangor, `Thomas` and `Fremont` mandarins and `Rangpur` lime (only in vitro segments of epicotyl) and to compare the gamma-rays sensitivity among propagules and among varieties. The following doses were used: 0, 10, 20, 30, 40 and 50 Gy of gamma-rays. The parameters shoot mortality and shoot height, for experiment of bud irradiation, and number of regenerated shoots per explant and percentage of responsive explants, for experiment of peicotyl segments irradiation were evaluated after 60 days. The mutagenic doses tested in buds caused significant reduction in shoot height of all varieties. In the experiment of irradiation of in vitro segments of epicotyl, only `Murcott` tangor and `Rangpur` lime showed significant reduction in the number of regenerated shoots per explant, due to the increase of mutagen doses. Results indicate that radiossensitivity of in vitro (segments of epicotyl) and in vivo (buds) propagules is variable depending on the variety. Also, in some cases the in vitro propagules are more sensitive, to irradiation and in other cases, there is no differential sensibility.

Evaluation of grinding methods for pellets preparation aiming at the analysis of plant materials by laser induced breakdown spectrometry

It has been demonstrated that laser induced breakdown spectrometry (LIBS) can be used as an alternative method for the determination of macro (P, K. Ca, Mg) and micronutrients (B, Fe, Cu, Mn, Zn) in pellets of plant materials. However, information is required regarding the sample preparation for plant analysis by LIBS. In this work, methods involving cryogenic grinding and planetary ball milling were evaluated for leaves comminution before pellets preparation. The particle sizes were associated to chemical sample properties such as fiber and cellulose contents, as well as to pellets porosity and density. The pellets were ablated at 30 different sites by applying 25 laser pulses per site (Nd:YAG@1064 nm, 5 ns, 10 Hz, 25J cm(-2)). The plasma emission collected by lenses was directed through an optical fiber towards a high resolution echelle spectrometer equipped with an ICCD. Delay time and integration time gate were fixed at 2.0 and 4.5 mu s, respectively. Experiments carried out with pellets of sugarcane, orange tree and soy leaves showed a significant effect of the plant species for choosing the most appropriate grinding conditions. By using ball milling with agate materials, 20 min grinding for orange tree and soy, and 60 min for sugarcane leaves led to particle size distributions generally lower than 75 mu m. Cryogenic grinding yielded similar particle size distributions after 10 min for orange tree, 20 min for soy and 30 min for sugarcane leaves. There was up to 50% emission signal enhancement on LIBS measurements for most elements by improving particle size distribution and consequently the pellet porosity. (C) 2011 Elsevier B.V. All rights reserved.

ON TWO-DIMENSIONAL QUASITOPOLOGICAL FIELD THEORIES

We study a class of lattice field theories in two dimensions that includes gauge theories. We show that in these theories it is possible to implement a broader notion of local symmetry, based on semisimple Hopf algebras. A character expansion is developed for the quasitopological field theories, and partition functions are calculated with this tool. Expected values of generalized Wilson loops are defined and studied with the character expansion.

High-yield kraft pulping of Eucalyptus grandis Hill ex Maiden biotreated by Ceriporiopsis subvermispora under two different culture conditions

In the present study, it was evaluated how two different culture conditions for the biotreatment of Eucalyptus grandis by Ceriporiopsis subvermispora affect a subsequent high-yield kraft pulping process. Under the varied culture conditions investigated, different extracellular enzyme activities were observed. Manganese-peroxidase (MnP) secretion was 3.7 times higher in cultures supplemented with glucose plus corn-steep liquor (glucose/CSL) as compared to non-supplemented (NS) cultures. The biotreated samples underwent diverse levels of wood component degradation as losses of weight and lignin were increased in glucose/CSL cultures. Mass balances for lignin removal during kraft pulping showed that delignification was facilitated when both biotreated wood samples were cooked. Delignification efficiency did not correlate positively with MnP levels in the cultures. On the other hand, biopulps from NS and glucose/CSL cultures saved 27% and 38% beating time to achieve 288 Schopper-Riegler freeness during refining, respectively. Biopulps disposed of decreased tensile and tear resistances, thus easier refining of the biokraft pulps seems to be a consequence of less resistant fiber walls. Improved beatability of biopulps was tentatively related to short fibers and fines formation during refining. We suggest that to some extent polysaccharide depolymerization occurred during the biotreatment, which also resulted in diminished pulp yields in the case of glucose/CSL cultures.

Aqueous two-phase extraction using thermoseparating copolymer: a new system for phenolic compounds removal from hemicelullosic hydrolysate

BACKGROUND: The hydrolysis of hemicellulosic material can provide liquor with high xylose concentration (which can be used as a fermentation medium) and phenolic compounds (Phs), potentially immunostimulating compounds. However, these hydrolysates must be detoxified in order to remove the Phs that can act as inhibitors in bioconversions. RESULTS: Aqueous two-phase systems composed of thermoseparating copolymers were used for rice straw hydrolysate detoxification. The hydrolysis process was able to promote chemical breakdown of 85% of the total hemicellulose content, 14% of the cellulose, and 2% of the lignin. The hydrolysate obtained contained 19.7 g L-1 of xylose and several phenolic compounds, such as vanillin, vanillic acid, ferullic acid, etc. The phenolics extraction was studied as a function of copolymer molar mass (1100 g mol(-1), 2000 g mol(-1) and 2800 g mol(-1)), their percentages (from 5% to 50%) and Phs initial concentration. Phenolic compounds extraction of around 80% was obtained under the following conditions: 20% (w/w) and 35% (w/w) copolymer 1100 g mol-1, 35% (w/w) copolymer 2000 g mol(-1) and 35% (w/w) copolymer 2800 g mol(-1) at 25 degrees C. CONCLUSIONS: The results demonstrated the viability of this method for the removal of Phs from rice straw hydrolysate, which has potential uses in bioconversion processes. (c) 2007 Society of Chemical Industry.