Multi-trait and random regression mature weight heritability and breeding value estimates in Nelore cattle

Mature weight breeding values were estimated using a multi-trait animal model (MM) and a random regression animal model (RRM). Data consisted of 82 064 weight records from 8 145 animals, recorded from birth to eight years of age. Weights at standard ages were considered in the MM. All models included contemporary groups as fixed effects, and age of dam (linear and quadratic effects) and animal age as covariates. In the RRM, mean trends were modelled through a cubic regression on orthogonal polynomials of animal age and genetic maternal and direct and maternal permanent environmental effects were also included as random. Legendre polynomials of orders 4, 3, 6 and 3 were used for animal and maternal genetic and permanent environmental effects, respectively, considering five classes of residual variances. Mature weight (five years) direct heritability estimates were 0.35 (MM) and 0.38 (RRM). Rank correlation between sires' breeding values estimated by MM and RRM was 0.82. However, selecting the top 2% (12) or 10% (62) of the young sires based on the MM predicted breeding values, respectively 71% and 80% of the same sires would be selected if RRM estimates were used instead. The RRM modelled the changes in the (co) variances with age adequately and larger breeding value accuracies can be expected using this model.

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.

The Lobel-Komlos-Sos conjecture for trees of diameter 5 and for certain caterpillars

Loebl, Komlos, and Sos conjectured that if at least half the vertices of a graph G have degree at least some k is an element of N, then every tree with at most k edges is a subgraph of G. We prove the conjecture for all trees of diameter at most 5 and for a class of caterpillars. Our result implies a bound on the Ramsey number r( T, T') of trees T, T' from the above classes.

Do plant species influence soil CO(2) and N(2)O fluxes in a diverse tropical forest?

To test whether plant species influence greenhouse gas production in diverse ecosystems, we measured wet season soil CO(2) and N(2)O fluxes close to similar to 300 large (>35 cm in diameter at breast height (DBH)) trees of 15 species at three clay-rich forest sites in central Amazonia. We found that soil CO(2) fluxes were 38% higher near large trees than at control sites >10 m away from any tree (P < 0.0001). After adjusting for large tree presence, a multiple linear regression of soil temperature, bulk density, and liana DBH explained 19% of remaining CO(2) flux variability. Soil N(2)O fluxes adjacent to Caryocar villosum, Lecythis lurida, Schefflera morototoni, and Manilkara huberi were 84%-196% greater than Erisma uncinatum and Vochysia maxima, both Vochysiaceae. Tree species identity was the most important explanatory factor for N(2)O fluxes, accounting for more than twice the N(2)O flux variability as all other factors combined. Two observations suggest a mechanism for this finding: (1) sugar addition increased N(2)O fluxes near C. villosum twice as much (P < 0.05) as near Vochysiaceae and (2) species mean N(2)O fluxes were strongly negatively correlated with tree growth rate (P = 0.002). These observations imply that through enhanced belowground carbon allocation liana and tree species can stimulate soil CO(2) and N(2)O fluxes (by enhancing denitrification when carbon limits microbial metabolism). Alternatively, low N(2)O fluxes potentially result from strong competition of tree species with microbes for nutrients. Species-specific patterns in CO(2) and N(2)O fluxes demonstrate that plant species can influence soil biogeochemical processes in a diverse tropical forest.

BORON UPTAKE AND DISTRIBUTION IN FIELD GROWN CITRUS TREES

In low fertility tropical soils, boron (B) deficiency impairs fruit production. However, little information is available on the efficiency of nutrient application and use by trees. Therefore, this work verified the effects of soil and foliar applications of boron in a commercial citrus orchard. An experiment was conducted with fertigated 4-year-old `Valencia` sweet orange trees on `Swingle` citrumelo rootstock. Boron (isotopically-enriched 10B) was supplied to trees once or twice in the growing season, either dripped in the soil or sprayed on the leaves. Trees were sampled at different periods and separated into different parts for total B contents and 10B/11B isotope ratios analyses. Soil B applied via fertigation was more efficient than foliar application for the organs grown after the B fertilization. Recovery of labeled B by fruits was 21% for fertigation and 7% for foliar application. Residual effects of nutrient application in the grove were observed in the year after labeled fertilizer application, which greater proportions derived from the soil supply.

Nickel sorption capacity of ground xylem of Quercus ilex trees and effects of selected ligands present in the xylem sap

In this work the influence of four different ligands present in the xylem sap of Quercus ilex (histidine, citric, oxalic and aspartic acids) on Ni(II) adsorption by xylem was investigated. Grinded xylem was trapped in acrylic columns and solutions of Ni(II), in the absence and presence of the four ligands prepared in KNO(3) 0-1 mol L(-1) at pH 5.5, were percolated through the column. Aliquots of solutions were recovered in the column end for Ni determination by graphite furnace atomic absorption spectrometry (GFAAS). The experimental. data to describe Ni sorption by xylem in both the presence and absence of ligands was better explained by the Freundlich isotherm model. The decreasing affinity order of ligands for Ni was: oxalic acid > citric acid > histidine > aspartic acid. On the other hand, the Ni(II) adsorption by xylem increased following the inverse sequence of ligands. Potentiometric titrations of acidic groups were carried out to elucidate the sorption site groups available in Q. ilex xylem. The potentiometric titration has shown three sorption sites: pK(a) 2.6 (57.7% of the sorption sites), related to monobasic aliphatic carboxylic acids or nitrogen aromatic bases, pK(a) 8.1 (9.6%) and pK(a) 9.9 (32.7%), related to phenolic groups. (C) 2008 Elsevier GmbH. All rights reserved.

Seasonal and diurnal changes in photosynthetic limitation of young sweet orange trees

This study tests the hypothesis that potted sweet orange plants show a significant variation in photosynthesis over seasonal and diurnal cycles. even in well-hydrated conditions. This hypothesis was tested by measuring diurnal variations in leaf gas exchange, chlorophyll fluorescence, leaf water potential, and the responses of CO(2) assimilation to increasing air CO(2) concentrations in 1-year-old `Valencia` sweet orange scions grafted onto `Cleopatra` mandarin rootstocks during the winter and summer seasons in a subtropical climate. In addition, diurnal leaf gas exchange was evaluated under controlled conditions, with constant environmental conditions during both winter and summer. In relation to our hypothesis, a greater rate of photosynthesis is found during the summer compared to the winter. Reduced photosynthesis during winter was induced by cool night conditions, as the diurnal fluctuation of environmental conditions was not limiting. Low air and soil temperatures caused decreases in the stomatal conductance and in the rates of the biochemical reactions underlying photosynthesis (ribulose-1,5-bisphosphate (RuBP) carboxylation and RuBP regeneration) during the winter compared to the values obtained for those markers in the Summer. Citrus photosynthesis during the summer was riot impaired by biochemical or photochemical reactions. as CO(2) assimilation was only limited by stomatal conductance due to high leaf-to-air vapor pressure difference (VPD) during the afternoon. During the winter, the reduction in photosynthesis during the afternoon Was Caused by decreases in RuBP regeneration and stomatal conductance, which are both precipitated by low night temperature. (c) 2009 Elsevier B.V. All rights reserved.

Log-modified Weibull regression models with censored data: Sensitivity and residual analysis

This paper proposes a regression model considering the modified Weibull distribution. This distribution can be used to model bathtub-shaped failure rate functions. Assuming censored data, we consider maximum likelihood and Jackknife estimators for the parameters of the model. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and we also present some ways to perform global influence. Besides, for different parameter settings, sample sizes and censoring percentages, various simulations are performed and the empirical distribution of the modified deviance residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for a martingale-type residual in log-modified Weibull regression models with censored data. Finally, we analyze a real data set under log-modified Weibull regression models. A diagnostic analysis and a model checking based on the modified deviance residual are performed to select appropriate models. (c) 2008 Elsevier B.V. All rights reserved.

On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

Regression models for grouped survival data: Estimation and sensitivity analysis

In this study, regression models are evaluated for grouped survival data when the effect of censoring time is considered in the model and the regression structure is modeled through four link functions. The methodology for grouped survival data is based on life tables, and the times are grouped in k intervals so that ties are eliminated. Thus, the data modeling is performed by considering the discrete models of lifetime regression. The model parameters are estimated by using the maximum likelihood and jackknife methods. To detect influential observations in the proposed models, diagnostic measures based on case deletion, which are denominated global influence, and influence measures based on small perturbations in the data or in the model, referred to as local influence, are used. In addition to those measures, the local influence and the total influential estimate are also employed. Various simulation studies are performed and compared to the performance of the four link functions of the regression models for grouped survival data for different parameter settings, sample sizes and numbers of intervals. Finally, a data set is analyzed by using the proposed regression models. (C) 2010 Elsevier B.V. All rights reserved.

Large Woody Debris Input and Its Influence on Channel Structure in Agricultural Lands of Southeast Brazil

Riparian forests are important for the structure and functioning of stream ecosystems, providing structural components such as large woody debris (LWD). Changes in these forests will cause modifications in the LWD input to streams, affecting their structure. In order to assess the influence of riparian forests changes in LWD supply, 15 catchments (third and fourth order) with riparian forests at different conservation levels were selected for sampling. In each catchment we quantified the abundance, volume and diameter of LWD in stream channels; the number, area and volume of pools formed by LWD and basal area and tree diameter of riparian forest. We found that riparian forests were at a secondary successional stage with predominantly young trees (diameter at breast height < 10 cm) in all studied streams. Results showed that basal area and diameter of riparian forest differed between the stream groups (forested and non-forested), but tree density did not differ between groups. Differences were also observed in LWD abundance, volume, frequency of LWD pools with subunits and area and volume of LWD pools. LWD diameter, LWD that form pools diameter and frequency of LWD pools without subunits did not differ between stream groups. Regression analyses showed that LWD abundance and volume, and frequency of LWD pools (with and without subunits) were positively related with the proportion of riparian forest. LWD diameter was not correlated to riparian tree diameter. The frequency of LWD pools was correlated to the abundance and volume of LWD, but characteristics of these pools (area and volume) were not correlated to the diameter of LWD that formed the pools. These results show that alterations in riparian forest cause modifications in the LWD abundance and volume in the stream channel, affecting mainly the structural complexity of these ecosystems (reduction in the number and structural characteristics of LWD pools). Our results also demonstrate that riparian forest conservation actions must consider not only its extension, but also successional stage to guarantee the quantity and quality of LWD necessary to enable the structuring of stream channels.

Association between decay fungi and subterranean termites in the wood biodeterioration process of Tipuana tipu (Benth.) O. Kuntze trees of Sao Paulo city, SP

In the present paper the process of wood biodeterioration of tipuana trees planted in 7 regions of the city of Sao Paulo, SP was evaluated. On the sidewalks, 1109 trees were analyzed taking into consideration the occurrence and association of the xylophagous organisms (decay fungi and subterranean termites), the wood deterioration and the BHD (breast height diameter). The percentage of wood internal deterioration (%) was obtained by non destructive analysis, using a penetrometer. The results had shown that 75% of the tipuana trees presented BHD superior to 50 cm, characterizing them as adult. Decay fungi in the roots and/or trunk had been observed in 338 trees (30.5%). Subterranean termites of Heterotermes sp. and Coptotermes gestroi species had occurred in 307 trees (27.7%), the latter in high infestation level. The association between the fungi and termites was observed, as well as its relation with the BHD, where a greater value of BHD meant higher wood biodeterioration intensity. For tipuana trees, the BHD was considered an indicative attribute of the internal deterioration intensity, caused by these xylophagous organisms.

Radial variation of the wood anatomical structure of Gmelina arborea trees from different climatic and management conditions in Costa Rica

The evaluations of the effect of the climatic conditions and of the intensity of forest management in the trunk of the Gmelina arborea Linn. Roxb. trees are restricted to its physical-mechanical properties and use. The present work has as objective to study the radial variations of the wood anatomy of the gmelina trees sampled in plantations of 30 sites in Costa Rica, characterized by two climatic conditions (tropical dry and humid) and three intensities of forest management (intensive, moderate and without management). The results of the analyses demonstrated the existence of radial variation of the different anatomical parameters, except for the fiber lumen diameter and multiple vessels in the wood of the gmelina trees. For the wood anatomical elements, fibers (width, lumen diameter, and length), vessels (multiple vessels, diameter and frequency) and radial parenchyma (height) relationships were observed with the climate (tropical humid and dry). The radial variations of the wood anatomical elements were, also, influenced by the management regimes of the gmelina trees.

Why don`t our stands grow even faster? Control of production and carbon cycling in eucalypt plantations

The growth of Eucalyptus stands varies several fold across sites, under the influence of resource availability, stand age and stand structure. We describe a series of related studies that aim to understand the mechanisms that drive this great range in stand growth rates. In a seven-year study in Hawaii of Eucalyptus saligna at a site that was not water limited, we showed that nutrient availability differences led to a two-fold difference in stand wood production. Increasing nutrient supply in mid-rotation raised productivity to the level attained in continuously fertilised plots. Fertility affected the age-related decline in wood and foliage production; production in the intensive fertility treatments declined more slowly than in the minimal fertility treatments. The decline in stem production was driven largely by a decline in canopy photosynthesis. Over time, the fraction of canopy photosynthesis partitioned to below-ground allocation increased, as did foliar respiration, further reducing wood production. The reason for the decline in photosynthesis was uncertain, but it was not caused by nutrient limitation, a decline in leaf area or in photosynthetic capacity, or by hydraulic limitation. Most of the increase in carbon stored from conversion of the sugarcane plantation to Eucalyptus plantation was in the above-ground woody biomass. Soil carbon showed no net change. This study and other studies on carbon allocation showed that resource availability changes the fraction of annual photosynthesis used below-ground and for wood production. High resources (nutrition or water) decrease the partitioning below-ground and increase partitioning to wood production. Annual foliage and wood respiration and foliage production as a fraction of annual photosynthesis was remarkably constant across a wide range of fertility treatments and forest age. In the Brazil Eucalyptus Productivity Project, stand structure was manipulated by planting clonal Eucalyptus all at once or in three groups at three-monthly intervals, producing a stand where trees did not segregate into dominants and one that had strong dominance. The uneven stand structure reduced production 10-15% throughout the rotation.

Variation of wood density and relationship with the tree-ring width of mahogany trees, Swietenia macrophylla, from Amazonian tropical forest of Peru

Mahogany trees, Swietenia macrophylla, occur in open rainforest, semi deciduous and deciduous and dense rainforest of Peruvian Amazonian tropical forest. They occur, preferentially, in areas with a defined dry season, with typical phenology and seasonal variation activity, forming distinct tree-rings. The present work had as aim to determine the wood density radial variation of 14 mahogany trees, of two populations of the Peruvian Amazonian tropical forest, through the X-ray densitometry and to evaluate their application as methodology, compared to the classic method of measurement table, for the determination of the treering width. The radial wood apparent density of the trees profiles rendered it possible to delimit the areas of juvenile-adult wood and of the heartwood-sapwood, relative to the anatomical structure and chemical composition differences, due to the extractives and the vessels obstruction by tyloses. The mean, minimum and maximum wood apparent density of the mahogany trees for the Populations A and B were of 0.70; 0.29; 1.01 g.cm(-3) and 0.81; 0.29; 1.19 g.cm(-3), respectively. The analysis of the variance and mean test indicate differences of mean wood density among the mahogany trees of each population, probably due to the age of the trees. There was no correlation between mean wood density of mahogany trees among the two populations, as well as, between the tree-ring width and the respective mean density. The X-ray densitometry technique is an important tool in the evaluation of the radial variation of wood apparent density and the delimitation of tree-ring boundaries, with correlations of 0.94 and 0.93 in relation to measurement table, for each sampled population.