Trait contributions to fish community assembly emerge from trophic interactions in an individual-based model

Community ecology seeks to understand and predict the characteristics of communities that can develop under different environmental conditions, but most theory has been built on analytical models that are limited in the diversity of species traits that can be considered simultaneously. We address that limitation with an individual-based model to simulate assembly of fish communities characterized by life history and trophic interactions with multiple physiological tradeoffs as constraints on species performance. Simulation experiments were carried out to evaluate the distribution of 6 life history and 4 feeding traits along gradients of resource productivity and prey accessibility. These experiments revealed that traits differ greatly in importance for species sorting along the gradients. Body growth rate emerged as a key factor distinguishing community types and defining patterns of community stability and coexistence, followed by egg size and maximum body size. Dominance by fast-growing, relatively large, and fecund species occurred more frequently in cases where functional responses were saturated (i.e. high productivity and/or prey accessibility). Such dominance was associated with large biomass fluctuations and priority effects, which prevented richness from increasing with productivity and may have limited selection on secondary traits, such as spawning strategies and relative size at maturation. Our results illustrate that the distribution of species traits and the consequences for community dynamics are intimately linked and strictly dependent on how the benefits and costs of these traits are balanced across different conditions. © 2012 Elsevier B.V.

Principal components and factor analytic models for birth to mature weights in Nellore cattle

A total of 61,528 weight records from 22,246 Nellore animals born between 1984 and 2002 were used to compare different multiple-trait analysis methods for birth to mature weights. The following models were used: standard multivarite model (MV), five reduced-rank models fitting the first 1, 2, 3, 4 and 5 genetic principal components, and five models using factor analysis with 1, 2, 3, 4 and 5 factors. Direct additive genetic random effects and residual effects were included in all models. In addition, maternal genetic and maternal permanent environmental effects were included as random effects for birth and weaning weight. The models included contemporary group as fixed effect and age of animal at recording (except for birth weight) and age of dam at calving as linear and quadratic effects (for birth weight and weaning weight). The maternal genetic, maternal permanent environmental and residual (co)variance matrices were assumed to be full rank. According to model selection criteria, the model fitting the three first principal components (PC3) provided the best fit, without the need for factor analysis models. Similar estimates of phenotypic, direct additive and maternal genetic, maternal permanent environmental and residual (co)variances were obtained with models MV and PC3. Direct heritability ranged from 0.21 (birth weight) to 0.45 (weight at 6 years of age). The genetic and phenotypic correlations obtained with model PC3 were slightly higher than those estimated with model MV. In general, the reduced-rank model substantially decreased the number of parameters in the analyses without reducing the goodness-of-fit. © 2013 Elsevier B.V.

Linear and Poisson models for genetic evaluation of tick resistance in cross-bred Hereford x Nellore cattle

Cattle resistance to ticks is measured by the number of ticks infesting the animal. The model used for the genetic analysis of cattle resistance to ticks frequently requires logarithmic transformation of the observations. The objective of this study was to evaluate the predictive ability and goodness of fit of different models for the analysis of this trait in cross-bred Hereford x Nellore cattle. Three models were tested: a linear model using logarithmic transformation of the observations (MLOG); a linear model without transformation of the observations (MLIN); and a generalized linear Poisson model with residual term (MPOI). All models included the classificatory effects of contemporary group and genetic group and the covariates age of animal at the time of recording and individual heterozygosis, as well as additive genetic effects as random effects. Heritability estimates were 0.08 ± 0.02, 0.10 ± 0.02 and 0.14 ± 0.04 for MLIN, MLOG and MPOI models, respectively. The model fit quality, verified by deviance information criterion (DIC) and residual mean square, indicated fit superiority of MPOI model. The predictive ability of the models was compared by validation test in independent sample. The MPOI model was slightly superior in terms of goodness of fit and predictive ability, whereas the correlations between observed and predicted tick counts were practically the same for all models. A higher rank correlation between breeding values was observed between models MLOG and MPOI. Poisson model can be used for the selection of tick-resistant animals. © 2013 Blackwell Verlag GmbH.

Estimates of genetic parameters for total milk yield over multiple ages in Brazilian Murrah buffaloes using different models

The objective of this study was to estimate variance components and genetic parameters for accumulated 305-day milk yield (MY305) over multiple ages, from 24 to 120 months of age, applying random regression (RRM), repeatability (REP) and multi-trait (MT) models. A total of 4472 lactation records from 1882 buffaloes of the Murrah breed were utilized. The contemporary group (herd-year-calving season) and number of milkings (two levels) were considered as fixed effects in all models. For REP and RRM, additive genetic, permanent environmental and residual effects were included as random effects. MT considered the same random effects as did REP and RRM with the exception of permanent environmental effect. Residual variances were modeled by a step function with 1, 4, and 6 classes. The heritabilities estimated with RRM increased with age, ranging from 0.19 to 0.34, and were slightly higher than that obtained with the REP model. For the MT model, heritability estimates ranged from 0.20 (37 months of age) to 0.32 (94 months of age). The genetic correlation estimates for MY305 obtained by RRM (L23.res4) and MT models were very similar, and varied from 0.77 to 0.99 and from 0.77 to 0.99, respectively. The rank correlation between breeding values for MY305 at different ages predicted by REP, MT, and RRM were high. It seems that a linear and quadratic Legendre polynomial to model the additive genetic and animal permanent environmental effects, respectively, may be sufficient to explain more parsimoniously the changes in MY305 genetic variation with age.

Predicting breeding values for milk yield of Guzera (Bos indicus) cows using random regression models

Given the importance of Guzera breeding programs for milk production in the tropics, the objective of this study was to compare alternative random regression models for estimation of genetic parameters and prediction of breeding values. Test-day milk yields records (TDR) were collected monthly, in a maximum of 10 measurements. The database included 20,524 records of first lactation from 2816 Guzera cows. TDR data were analyzed by random regression models (RRM) considering additive genetic, permanent environmental and residual effects as random and the effects of contemporary group (CG), calving age as a covariate (linear and quadratic effects) and mean lactation curve as fixed. The genetic additive and permanent environmental effects were modeled by RRM using Wilmink, All and Schaeffer and cubic B-spline functions as well as Legendre polynomials. Residual variances were considered as heterogeneous classes, grouped differently according to the model used. Multi-trait analysis using finite-dimensional models (FDM) for testday milk records (TDR) and a single-trait model for 305-days milk yields (default) using the restricted maximum likelihood method were also carried out as further comparisons. Through the statistical criteria adopted, the best RRM was the one that used the cubic B-spline function with five random regression coefficients for the genetic additive and permanent environmental effects. However, the models using the Ali and Schaeffer function or Legendre polynomials with second and fifth order for, respectively, the additive genetic and permanent environmental effects can be adopted, as little variation was observed in the genetic parameter estimates compared to those estimated by models using the B-spline function. Therefore, due to the lower complexity in the (co)variance estimations, the model using Legendre polynomials represented the best option for the genetic evaluation of the Guzera lactation records. An increase of 3.6% in the accuracy of the estimated breeding values was verified when using RRM. The ranks of animals were very close whatever the RRM for the data set used to predict breeding values. Considering P305, results indicated only small to medium difference in the animals' ranking based on breeding values predicted by the conventional model or by RRM. Therefore, the sum of all the RRM-predicted breeding values along the lactation period (RRM305) can be used as a selection criterion for 305-day milk production. (c) 2014 Elsevier B.V. All rights reserved.

Random regression models using different functions to model milk flow in dairy cows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic analysis of the temperament of Nellore cattle using linear and threshold models

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genotype x environment interaction for age at first calving, scrotal circumference, and yearling weight in Nellore cattle using reaction norms in multitrait random regression models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Quantitative trait loci with sex-specific effects for internal organs weights and hematocrit value in a broiler-layer cross

Rapid growth in broilers is associated with susceptibility to metabolic disorders such as pulmonary hypertension syndrome (ascites) and sudden death. This study describes a genome search for QTL associated with relative weight of cardio respiratory and metabolically important organs (heart, lungs, liver and gizzard), and hematocrit value in a Brazilian broiler-layer cross. QTL with similar or different effects across sexes were investigated. At 42 days of age after fasted for 6 h, the F2 chickens were weighed and slaughtered. Weights and percentages of the weight relative to BW42 of gizzard, heart, lungs, liver and hematocrit were used in the QTL search. Parental, F1 and F2 individuals were genotyped with 128 genetic markers (127 microsatellites and 1 SNP) covering 22 linkage groups. QTL mapping analyses were carried out using mixed models. A total of 11 genome-wide significant QTL and five suggestive linkages were mapped. Thus, genome-wide significant QTL with similar effects across sexes were mapped to GGA2, 4 and 14 for heart weight, and to GGA2, 8 and 12 for gizzard %. Additionally, five genome-wide significant QTL with different effects across sexes were mapped to GGA 8, 19 and 26 for heart weight; GGA26 for heart % and GGA3 for hematocrit value. Five QTL were detected in chromosomal regions where QTL for similar traits were previously mapped in other F2 chicken populations. Seven novel genome-wide significant QTL are reported here, and 21 positional candidate genes in QTL regions were identified.

Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Simulating the evolution of a clonal trait in plants with sexual and vegetative reproduction

Aims Phenotypic optimality models neglect genetics. However, especially when heterozygous genotypes ire fittest, evolving allele, genotype and phenotype frequencies may not correspond to predicted optima. This was not previously addressed for organisms with complex life histories. Methods Therefore, we modelled the evolution of a fitness-relevant trait of clonal plants, stolon internode length. We explored the likely case of air asymmetric unimodal fitness profile with three model types. In constant selection models (CSMs), which are gametic, but not spatially explicit, evolving allele frequencies in the one-locus and five-loci cases did not correspond to optimum stolon internode length predicted by the spatially explicit, but not gametic, phenotypic model. This deviation was due to the asymmetry of the fitness profile. Gametic, spatially explicit individual-based (SEIB) modeling allowed us relaxing the CSM assumptions of constant selection with exclusively sexual reproduction. Important findings For entirely vegetative or sexual reproduction, predictions. of the gametic SEIB model were close to the ones of spatially explicit CSMs gametic phenotypic models, hut for mixed modes of reproduction they appoximated those of gametic, not spatially explicit CSMs. Thus, in contrast to gametic SEIB models, phenotypic models and, especially for few loci, also CSMs can be very misleading. We conclude that the evolution of trails governed by few quantitative trait loci appears hardly predictable by simple models, that genetic algorithms aiming at technical optimization may actually, miss the optimum and that selection may lead to loci with smaller effects, in derived compared with ancestral lines.

Chapter Six – Individual Trait Variation and Diversity in Food Webs

In recent years, there has been a renewed interest in the ecological consequences of individual trait variation within populations. Given that individual variability arises from evolutionary dynamics, to fully understand eco-evolutionary feedback loops, we need to pay special attention to how standing trait variability affects ecological dynamics. There is mounting empirical evidence that intra-specific phenotypic variation can exceed species-level means, but theoretical models of multi-trophic species coexistence typically neglect individual-level trait variability. What is needed are multispecies datasets that are resolved at the individual level that can be used to discriminate among alternative models of resource selection and species coexistence in food webs. Here, using one the largest individual-based datasets of a food web compiled to date, along with an individual trait-based stochastic model that incorporates Approximate Bayesian computation methods, we document intra-population variation in the strength of prey selection by different classes or predator phenotypes which could potentially alter the diversity and coexistence patterns of food webs. In particular, we found that strongly connected individual predators preferentially consumed common prey, whereas weakly connected predators preferentially selected rare prey. Such patterns suggest that food web diversity may be governed by the distribution of predator connectivity and individual trait variation in prey selection. We discuss the consequences of intra-specific variation in prey selection to assess fitness differences among predator classes (or phenotypes) and track longer term food web patterns of coexistence accounting for several phenotypes within each prey and predator species.

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G x E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G x E interaction is not necessary for the maintenance of more than two alleles in two demes.

Seasonal Variation in the Capacity for Plant Trait Measures to Predict Grassland Carbon and Water Fluxes

There is a need for accurate predictions of ecosystem carbon (C) and water fluxes in field conditions. Previous research has shown that ecosystem properties can be predicted from community abundance-weighted means (CWM) of plant functional traits and measures of trait variability within a community (FDvar). The capacity for traits to predict carbon (C) and water fluxes, and the seasonal dependency of these trait-function relationships has not been fully explored. Here we measured daytime C and water fluxes over four seasons in grasslands of a range of successional ages in southern England. In a model selection procedure, we related these fluxes to environmental covariates and plant biomass measures before adding CWM and FDvar plant trait measures that were scaled up from measures of individual plants grown in greenhouse conditions. Models describing fluxes in periods of low biological activity contained few predictors, which were usually abiotic factors. In more biologically active periods, models contained more predictors, including plant trait measures. Field-based plant biomass measures were generally better predictors of fluxes than CWM and FDvar traits. However, when these measures were used in combination traits accounted for additional variation. Where traits were significant predictors their identity often reflected seasonal vegetation dynamics. These results suggest that database derived trait measures can improve the prediction of ecosystem C and water fluxes. Controlled studies and those involving more detailed flux measurements are required to validate and explore these findings, a worthwhile effort given the potential for using simple vegetation measures to help predict landscape-scale fluxes.

Association of Trait Resilience With Peritraumatic and Posttraumatic Stress in Patients With Myocardial Infarction.

OBJECTIVE Acute myocardial infarction (MI) is a life-threatening condition, leading to immediate fear and distress in many patients. Approximately 18% of patients develop posttraumatic stress disorder in the aftermath of MI. Trait resilience has shown to be a protective factor for the development of posttraumatic stress disorder. However, whether this buffering effect has already an impact on peritraumatic distress and applies to patients with MI is elusive. METHODS We investigated 98 consecutive patients with acute MI within 48 hours after having reached stable circulatory conditions and 3 months thereafter. Peritraumatic distress was assessed retrospectively with three single-item questions about pain, fear, and helplessness during MI. All patients completed the Posttraumatic Diagnostic Scale (PDS) and the Resilience Scale to self-rate posttraumatic stress and trait resilience. RESULTS Multivariate models adjusting for sociodemographic and medical factors showed that trait resilience was not associated with peritraumatic distress, but significantly so with posttraumatic stress. Patients with greater trait resilience showed lower PDS scores (b = -0.06, p < .001). There was no significant relationship between peritraumatic distress scores and PDS scores; resilience did not emerge as a moderator of this relationship. CONCLUSIONS The findings suggest that trait resilience does not buffer the perception of acute MI as stressful per se but may enhance better coping with the traumatic experience in the longer term, thus preventing the development of MI-associated posttraumatic stress. Trait resilience may play an important role in posttraumatic stress symptoms triggered by medical diseases such as acute MI.