Analysis of covariance: an alternative to nutritional indices

Some statistical problems are added to the growing list of cautionary tales regarding the use of the conventional, ratio-based nutritional indices (RCR, RGR, ECI, AD and ECD). Analysis of ratios is based on the, probably unrealistic, assumption of an isometric relationship between denominator and numerator variables. Analysis of covariance (ANCOVA) makes less restrictive assumptions, and additionally provides important information about the data which is lost by using ratio variables. We demonstrate, using computer-generated data sets, some of the pitfalls of statistical analysis of ratios and illustrate how these may be avoided using ANCOVA. Some possible consequences of such statistical iniquities for biological interpretations are discussed.     

Factor analysis models for structuring covariance matrices of additive genetic effects: a Bayesian implementation

Multivariate linear models are increasingly important in quantitative genetics. In high dimensional specifications, factor analysis (FA) may provide an avenue for structuring (co)variance matrices, thus reducing the number of parameters needed for describing (co)dispersion. We describe how FA can be used to model genetic effects in the context of a multivariate linear mixed model. An orthogonal common factor structure is used to model genetic effects under Gaussian assumption, so that the marginal likelihood is multivariate normal with a structured genetic (co)variance matrix. Under standard prior assumptions, all fully conditional distributions have closed form, and samples from the joint posterior distribution can be obtained via Gibbs sampling. The model and the algorithm developed for its Bayesian implementation were used to describe five repeated records of milk yield in dairy cattle, and a one common FA model was compared with a standard multiple trait model. The Bayesian Information Criterion favored the FA model.     

Application of a self-organizing map to select representative species in multivariate analysis: A case study determining diatom distribution patterns across France

Ecological communities consist of a large number of species. Most species are rare or have low abundance, and only a few are abundant and/or frequent. In quantitative community analysis, abundant species are commonly used to interpret patterns of habitat disturbance or ecosystem degradation. Rare species cause many difficulties in quantitative analysis by introducing noises and bulking datasets, which is worsened by the fact that large datasets suffer from difficulties of data handling. In this study we propose a method to reduce the size of large datasets by selecting the most ecologically representative species using a self organizing map (SOM) and structuring index (SI). As an example, we used diatom community data sampled at 836 sites with 941 species throughout the French hydrosystem. Out of the 941 species, 353 were selected. The selected dataset was effectively classified according to the similarities of community assemblages in the SOM map. Compared to the SOM map generated with the original dataset, the community pattern gave a very similar representation of ecological conditions of the sampling sites, displaying clear gradients of environmental factors between different clusters. Our results showed that this computational technique can be applied to preprocessing data in multivariate analysis. It could be useful for ecosystem assessment and management, helping to reduce both the list of species for identification and the size of datasets to be processed for diagnosing the ecological status of water courses.     

Mathematical hazard models of mortality: an alternative to model life tables

A five-parameter competing hazard model of the age pattern of mortality is described, and methods of fitting it to survivorship, death rate, and age structure data are developed and presented. The methods are then applied to published life table and census data to construct life tables for a Late Woodland population, a Christian period Nubian population, and the Yanomama. The advantage of this approach over the use of model life tables is that the hazard model facilitates life-table construction without imposing a particular age pattern of mortality on the data. This development makes it possible to use anthropological data to extend the study of human variation in mortality patterns to small populations.     

Dynamical entropy via entropy of non-random matrices: Application to stability and complexity in modelling ecosystems

In the present paper we have first introduced a measure of dynamical entropy of an ecosystem on the basis of the dynamical model of the system. The dynamical entropy which depends on the eigenvalues of the community matrix of the system leads to a consistent measure of complexity of the ecosystem to characterize the dynamical behaviours such as the stability, instability and periodicity around the stationary states of the system. We have illustrated the theory with some model ecosystems.     

Selecting models to predict the timing of flowering of temperate trees: implications for tree phenology modelling

Classical budburst models (Spring Warming, Sequential, Parallel and Alternating) are unable to fully predict external data, partly because of the methods of optimization used to adjust them. The purpose of this study was to examine different assumptions of budburst models and select those which are best supported by the data, defining new models able to predict external data. Eight models, each differing in one assumption, were fitted and tested using external data. The dataset used to test the models was deduced from aeropalynological data at two stations in France. The results show that some of the models proposed are able to accurately predict external dates of flowering of most of the studied species. The assumptions of those models have been individually tested and shown to improve the models accuracy. Robust estimates of the best predictor models of 12 tree species are presented. The analysis of hypothetical provenance transfer of two species, Buxus sempervirens and Platanus acerifolia, between the two study sites, shows that P. acerifolia estimates are similar in both environments whereas B. sempervirens estimates are variable. This result, which agrees with the genetic characteristics of both species, shows that local adaptation of phenology can also be studied through modelling approaches.     

Prediction of biogas potentials using quick laboratory analyses: Upgrading previous models for application to heterogeneous organic matrices

This study presents an upgrading of the mathematical models to predict anaerobic biogasification potential (ABP) through quick laboratory analyses that have been presented in an earlier study. The aim is to widen the applicability of the models to heterogeneous organic substrates and to improve their reliability through a deeper statistical approach.Three multiple-step linear regressions were obtained using biomass oxygen demand in 20 h (OD₂₀) plus the volatile solids content (VS) of 23 new samples of heterogeneous organic matrices, of 46 samples presented in the earlier work and of the data set comprising all the 69 samples. The two variables chosen were found to be suitable for very heterogeneous materials. To judge the prediction quality, a validation procedure was performed with 12 new samples using model efficiency indexes. The proposed model had good prediction ability for a large variety of organic substrates, and allows the calculation of the ABP value within only 2-day’s laboratory work instead of the 60–90 days required to obtain ABP by anaerobic test.     

Infection of organotypic slice cultures from rat central nervous tissue with Neospora caninum: an alternative approach to study host-parasite interactions

Neospora caninum is an apicomplexan parasite which has emerged as an important cause of bovine abortion worldwide. Abortion is usually triggered by reactivation of dormant bradyzoites during pregnancy and subsequent congenital infection of the foetus, where the central nervous system appears to be most frequently affected. We here report on an organotypic tissue culture model for Neospora infection which can be used to study certain aspects of the cerebral phase of neosporosis within the context of a three-dimensionally organised neuronal network. Organotypic slice cultures of rat cortical tissue were infected with N. caninum tachyzoites, and the kinetics of parasite proliferation, as well as the proliferation-inhibitory effect of interferon-gamma (IFN-gamma), were monitored by either immunofluorescence, transmission electron microscopy, and a quantitative PCR-assay using the LightCycler instrument, respectively. In addition, the neuronal cytoskeletal elements, namely glial acidic protein filaments as well as actin microfilament bundles were shown to be largely colocalising with the pseudocyst periphery. This organotypic culture model for cerebral neosporosis provides a system, which is useful to study the proliferation, ultrastructural characteristics, development, and the interactions of N. caninum within the context of neuronal tissue, which at the same time can be modulated and influenced under controlled conditions, and will be useful in the future to gain more information on the cerebral phase of neosporosis.     

Rapid immunofluorescent determination of cells in the S phase in pea root meristems: An alternative to autoradiography

Photosystem II (PS II) activity and the localization of ribulose-l,5-bisphosphate (RuBP) carboxylase (EC 4.1.1.39) were studied in primary leaves of young maize plants ( Zea mays L. cv. Fronica) by tetra-nitro-blue-tetrazoliumchloride reduction and immunolocalization, respectively. In tissue of 3-day-old plants all chloroplasts were structurally identical. From day 4 they developed into their typical appearance of mesophyll and bundle sheath chloroplasts. First PS II-activity was present in both types of chloroplasts. From day 4 it disappeared in bundle sheath chloroplasts concomitant with the loss of grana. RuBP carboxylase on the other hand was only present in bundle sheath chloroplasts at all stages of development. Thus, the control of the development of the photosystems and the Calvin cycle enzymes seem to differ.     

Aerosol infectivity of a Baculovirus to Trichoplusia ni larvae: An alternative larval inoculation strategy for recombinant protein production

The baculovirus–insect expression system is a popular tool for recombinant protein production. The standard method for infecting insect larvae with recombinant baculovirus for protein production involves either feeding occlusion bodies or injecting budded virus into the cuticle. In this study, we showed that the recombinant Autographa californica multiple nucleopolyhedrovirus (AcMNPV) at titers >10⁸ pfu/mL efficiently infected Trichoplusia ni (T. ni) larvae through aerosol inoculation of budded virus at a pressure of 5.5 × 10⁴ Pa. The dipping T. ni larvae in virus‐containing solution efficiently infected them. These results indicate that surface contamination, either by aerosol or dipping, lead to infection via spiracles. The aerosol infection route for AcMNPV was restricted to T. ni and Plutella xylostella larvae, whereas Spodoptera litura and Helicoverpa armigera larvae were resistant to this inoculation process. The yields of the reporter proteins DsRed and EGFP from T. ni larvae following aerosol infection were nearly identical to those following oral feeding or injection. This alternative baculovirus infection strategy facilitates recombinant protein and virus production by insect larvae. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2009     

Different Ks values for hydrogen of methanogenic bacteria and sulfate reducing bacteria: An explanation for the apparent inhibition of methanogenesis by sulfate

Desulfovibrio vulgaris (Marburg) and Methanobrevibacter arboriphilus (AZ) are anaerobic sewage sludge bacteria which grow on H₂ plus sulfate and H₂ plus CO₂ as sole energy sources, respectively. Their apparent K_s values for H₂ were determined and found to be approximately 1 M for the sulfate reducing bacterium and 6 M for the methanogenic bacterium. In mixed cell suspensions of the two bacteria (adjusted to equal V _max) the rate of H₂ consumption by D. vulgaris was five times that of M. arboriphilus, when the hydrogen supply was rate limiting. The apparent inhibition of methanogenesis was of the same order as expected from the different K_s values for H₂. Difference in substrate affinities can thus account for the inhibition of methanogenesis from H₂ and CO₂ in sulfate rich environments, where the H₂ concentration is well below 5 M.     

An ensemble of models of the acute inflammatory response to bacterial lipopolysaccharide in rats: results from parameter space reduction

In previous work, we developed an 8-state nonlinear dynamic model of the acute inflammatory response, including activated phagocytic cells, pro- and anti-inflammatory cytokines, and tissue damage, and calibrated it to data on cytokines from endotoxemic rats. In the interest of parsimony, the present work employed parametric sensitivity and local identifiability analysis to establish a core set of parameters predominantly responsible for variability in model solutions. Parameter optimization, facilitated by varying only those parameters belonging to this core set, was used to identify an ensemble of parameter vectors, each representing an acceptable local optimum in terms of fit to experimental data. Individual models within this ensemble, characterized by their different parameter values, showed similar cytokine but diverse tissue damage behavior. A cluster analysis of the ensemble of models showed the existence of a continuum of acceptable models, characterized by compensatory mechanisms and parameter changes. We calculated the direct correlations between the core set of model parameters and identified three mechanisms responsible for the conversion of the diverse damage time courses to similar cytokine behavior in these models. Given that tissue damage level could be an indicator of the likelihood of mortality, our findings suggest that similar cytokine dynamics could be associated with very different mortality outcomes, depending on the balance of certain inflammatory elements.     

Cascading Trophic Interactions from Fish to Bacteria and Nutrients after Reduced Sewage Loading: An 18-Year Study of a Shallow Hypertrophic Lake

The effects of major reductions in organic matter, total phosphorus (TP), and total nitrogen (TN) loading on the chemical environment, trophic structure, and dynamics of the hypertrophic, shallow Lake Søbygård were followed for 18 years. After the reduction in organic matter loading in 1976, the lake initially shifted from a summer clear-water state, most likely reflecting high grazing pressure by large Daphniaspecies, to a turbid state with extremely high summer mean chlorophyll a (up to 1400 μg L^{? 1}), high pH (up to 10.2), and low zooplankton grazing. Subsequently, a more variable state with periodically high grazing rates on phytoplankton and bacteria was established. Changes in zooplankton abundance and grazing could be attributed to variations in cyprinid abundance due to a fish kill (probably as a consequence of oxygen depletion) and pH-induced variations in fish recruitment and fry survival. The results suggest strong cascading effects of fish on the abundance and size of zooplankton and phytoplankton and on phytoplankton production. A comparatively weak cascading effect on ciliates and bacterioplankton is suggested. Due to high internal loading, only minor changes were observed in lake-water TP after a reduction in external TP loading of approximately 80% in 1982; net retention of TP was still negative 13 years after the loading reduction, despite a short hydraulic retention time of a few weeks. TN, however, decreased proportionally to the TN-loading reduction in 1987, suggesting a fast N equilibration. Only minor improvement in the environmental state of the lake has been observed. We suggest that another decade will be required before the lake is in equilibrium with present external P loading.