Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

A quantitative genetic signature of senescence in a short-lived perennial plant.

The evolution of senescence (the physiological decline of organisms with age) poses an apparent paradox because it represents a failure of natural selection to increase the survival and reproductive performance of organisms. The paradox can be resolved if natural selection becomes less effective with age, because the death of postreproductive individuals should have diminished effects on Darwinian fitness [1, 2]. A substantial body of empirical work is consistent with this prediction for animals, which transmit their genes to progeny via an immortal germline. However, such evidence is still lacking in plants, which lack a germline and whose reproduction is diffuse and modular across the soma. Here, we provide experimental evidence for a genetic basis of senescence in the short-lived perennial plant Silene latifolia. Our pedigree-based analysis revealed a marked increase with age in the additive genetic variance of traits closely associated with fitness. This result thus extends to plants the quantitative genetic support for the evolutionary theory of senescence.

Heritability and quantitative genetic divergence of serotiny, a fire-persistence plant trait.

BACKGROUND AND AIMS: Although it is well known that fire acts as a selective pressure shaping plant phenotypes, there are no quantitative estimates of the heritability of any trait related to plant persistence under recurrent fires, such as serotiny. In this study, the heritability of serotiny in Pinus halepensis is calculated, and an evaluation is made as to whether fire has left a selection signature on the level of serotiny among populations by comparing the genetic divergence of serotiny with the expected divergence of neutral molecular markers (QST-FST comparison). METHODS: A common garden of P. halepensis was used, located in inland Spain and composed of 145 open-pollinated families from 29 provenances covering the entire natural range of P. halepensis in the Iberian Peninsula and Balearic Islands. Narrow-sense heritability (h(2)) and quantitative genetic differentiation among populations for serotiny (QST) were estimated by means of an 'animal model' fitted by Bayesian inference. In order to determine whether genetic differentiation for serotiny is the result of differential natural selection, QST estimates for serotiny were compared with FST estimates obtained from allozyme data. Finally, a test was made of whether levels of serotiny in the different provenances were related to different fire regimes, using summer rainfall as a proxy for fire regime in each provenance. KEY RESULTS: Serotiny showed a significant narrow-sense heritability (h(2)) of 0·20 (credible interval 0·09-0·40). Quantitative genetic differentiation among provenances for serotiny (QST = 0·44) was significantly higher than expected under a neutral process (FST = 0·12), suggesting adaptive differentiation. A significant negative relationship was found between the serotiny level of trees in the common garden and summer rainfall of their provenance sites. CONCLUSIONS: Serotiny is a heritable trait in P. halepensis, and selection acts on it, giving rise to contrasting serotiny levels among populations depending on the fire regime, and supporting the role of fire in generating genetic divergence for adaptive traits.

Quantitative genetic analysis of methylxanthines and phenolic compounds in mate progenies

The objective of this work was to determine the contents of methylxanthines, caffeine and theobromine, and phenolic compounds, chlorogenic and caffeic acids, in 51 mate progenies (half-sib families) and estimate the heritability of genetic parameters. Mate progenies were from five Brazilian municipalities: Pinhão, Ivaí, Barão de Cotegipe, Quedas do Iguaçu, and Cascavel. The progenies were grown in the Ivaí locality. The contents of the compounds were obtained by high performance liquid chromatography (HPLC). The estimation of genetic parameters by the restricted maximum likelihood (REML) and the prediction of genotypic values via best linear unbiased prediction (BLUP) were obtained by the Selegen - REML/BLUP software. Caffeine (0.248-1.663%) and theobromine (0.106-0.807%) contents were significantly different (p<0.05) depending on the region of origin, with high individual heritability (ĥ²>0.5). The two different progeny groups determined for chlorogenic (1.365-2.281%) and caffeic (0.027-0.037%) acid contents were not significantly different (p<0.05) depending on the locality of origin. Individual heritability values were low to medium for chlorogenic (ĥ²<0.4) and caffeic acid (ĥ²<0.3). The content of the compounds and the values of genetic parameters could support breeding programs for mate.

Quantitative genetic divergence and standing genetic (co)variance in thermal reaction norms along latitude

Although the potential to adapt to warmer climate is constrained by genetic trade-offs, our understanding of how selection and mutation shape genetic (co)variances in thermal reaction norms is poor. Using 71 isofemale lines of the fly Sepsis punctum, originating from northern, central, and southern European climates, we tested for divergence in juvenile development rate across latitude at five experimental temperatures. To investigate effects of evolutionary history in different climates on standing genetic variation in reaction norms, we further compared genetic (co)variances between regions. Flies were reared on either high or low food resources to explore the role of energy acquisition in determining genetic trade-offs between different temperatures. Although the latter had only weak effects on the strength and sign of genetic correlations, genetic architecture differed significantly between climatic regions, implying that evolution of reaction norms proceeds via different trajectories at high latitude versus low latitude in this system. Accordingly, regional genetic architecture was correlated to region-specific differentiation. Moreover, hot development temperatures were associated with low genetic variance and stronger genetic correlations compared to cooler temperatures. We discuss the evolutionary potential of thermal reaction norms in light of their underlying genetic architectures, evolutionary histories, and the materialization of trade-offs in natural environments.

Quantitative genetic study of age at subsequent rebreeding in Nellore cattle by using survival analysis

The continuous trait age at subsequent rebreeding (ASR) was evaluated using survival analysis in Nellore breed cows that conceived for the first time at approximately 14 months of age. This methodology was chosen because the restricted breeding season produces censored data. The dataset contained 2885 records of ASR (in days). Records of females that did not produce calves in the following year after being exposed to a sire were considered censored (48.3% of the total). The statistical model used was a Weibull mixed survival model, which included fixed effects of contemporary groups (CG) and period and a random effect of individual animal. The effect of contemporary groups on ASR was significant (P < 0.01). Heritabilities obtained for ASR were 0.03 and 0.04 in logarithmic and original scales, respectively. These results indicate that the genetic selection response for subsequent reproduction of 2-year-old Nellore breed females is not expected to be effective based on survival analysis. Furthermore, these results suggest that environmental improvement is fundamental to this important trait. It should be highlighted that an increase in the average date of birth can produce an adverse effect in the future, since this cannot be compensated by genetic improvement.

The quantitative genetic basis of offspring solicitation and parental response in a passerine bird with biparental care

The coevolution of parental investment and offspring solicitation is driven by partly different evolutionary interests of genes expressed in parents and their offspring. In species with biparental care, the outcome of this conflict ma!: be influenced by the sexual conflict over parental investment, Models for the resolution of such family conflicts have made so far untested assumptions about genetic variation and covariation in the parental resource provisioning response and the level of offspring solicitation. Using a combination of cross-fostering and begging playback experiments, we show that, in the great tit (Parus major), (i) the begging call intensity of nestlings depends on their common origin, suggesting genetic variation for this begging display, (ii) only mothers respond to begging calls by increased food provisioning, and (iii! the size of the parental response is positively related to the begging call intensity of nestlings in the maternal but not paternal line. This study indicates that genetic covariation, its differential expression in the maternal and paternal lines and/or early environmental and parental effects need to be taken into account when predicting the phenotypic outcome of the conflict over investment between genes expressed in each parent and the offspring. [References: 36]

Testing for biases in selection on avian reproductive traits and partitioning direct and indirect selection using quantitative genetic models

Key life history traits such as breeding time and clutch size are frequently both heritable and under directional selection, yet many studies fail to document micro-evolutionary responses. One general explanation is that selection estimates are biased by the omission of correlated traits that have causal effects on fitness, but few valid tests of this exist. Here we show, using a quantitative genetic framework and six decades of life-history data on two free-living populations of great tits Parus major, that selection estimates for egg-laying date and clutch size are relatively unbiased. Predicted responses to selection based on the Robertson-Price Identity were similar to those based on the multivariate breeder’s equation, indicating that unmeasured covarying traits were not missing from the analysis. Changing patterns of phenotypic selection on these traits (for laying date, linked to climate change) therefore reflect changing selection on breeding values, and genetic constraints appear not to limit their independent evolution. Quantitative genetic analysis of correlational data from pedigreed populations can be a valuable complement to experimental approaches to help identify whether apparent associations between traits and fitness are biased by missing traits, and to parse the roles of direct versus indirect selection across a range of environments.

Quantitative genetic analysis of grain yield in an Australian Brassica napus doubled-haploid population

Genetic influence on the variance in P3 amplitude and latency

The P3(00) event-related potential (ERP) component is widely used as a measure of cognitive functioning and provides a sensitive electrophysiological index of the attentional and working memory demands of a task. This study investigated what proportion of the variance in the amplitude and latency of the P3, elicited in a delayed response working memory task, could be attributed to genetic factors. In 335 adolescent twin pairs and 48 siblings, the amplitude and latency of the P3 were examined at frontal, central, and parietal sites. Additive genetic factors accounted for 48% to 61% of the variance in P3 amplitude. Approximately one-third of the genetic variation at frontal sites was mediated by a common genetic factor that also influenced the genetic variation at parietal and central sites. Familial resemblance in P3 latency was due to genetic influence that accounted for 44% to 50% of the variance. Genetic covariance in P3 latency across sites was substantial, with a large part of the variance found at parietal, central, and frontal sites attributed to a common genetic factor. The findings provide further evidence that the P3 is a promising phenotype of neural activity of the brain and has the potential to be used in linkage and association analysis in the search for quantitative trait loci (QTLs) influencing cognition.

Genetic epidemiology of fecal egg excretion during Schistosoma mansoni infection in an endemic area in Minas Gerais, Brazil

There is considerable variation in the level of fecal egg excretion during Schistosoma mansoni infections. Within a single endemic area, the distribution of egg counts is typically overdispersed, with the majority of eggs excreted coming from a minority of residents. The purpose of this study was to quantify the influence of genetic factors on patterns of fecal egg excretion in a rural study sample in Brazil. Individual fecal egg excretions, expressed in eggs per gram of feces, were determined by the Kato-Katz method on stool samples collected on three different days. Detailed genealogic information was gathered at the time of sampling, which allowed assignment of 461 individuals to 14 pedigrees containing between 3 and 422 individuals. Using a maximum likelihood variance decomposition approach, we performed quantitative genetic analyses to determine if genetic factors could partially account for the observed pattern of fecal egg excretion. The quantitative genetic analysis indicated that between 21-37% of the variation in S. mansoni egg counts was attributable to additive genetic factors and that shared environment, as assessed by common household, accounted for a further 12-21% of the observed variation. A maximum likelihood heritability (h²) estimate of 0.44 ± 0.14 (mean ± SE) was found for the 9,604 second- and higher-degree pairwise relationships in the study sample, which is consistent with the upper limit (37%) of the genetic factor determined in the variance decomposition analysis. These analyses point to the significant influence of additive host genes on the pattern of S. mansoni fecal egg excretion in this endemic area.

High genetic variability and low local diversity in a population of arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF) are ecologically important root symbionts of most terrestrial plants. Ecological studies of AMF have concentrated on differences between species; largely assuming little variability within AMF species. Although AMF are clonal, they have evolved to contain a surprisingly high within-species genetic variability, and genetically different nuclei can coexist within individual spores. These traits could potentially lead to within-population genetic variation, causing differences in physiology and symbiotic function in AMF populations, a consequence that has been largely neglected. We found highly significant genetic and phenotypic variation among isolates of a population of Glomus intraradices but relatively low total observed genetic diversity. Because we maintained the isolated population in a constant environment, phenotypic variation can be considered as variation in quantitative genetic traits. In view of the large genetic differences among isolates by randomly sampling two individual spores, <50% of the total observed population genetic diversity is represented. Adding an isolate from a distant population did not increase total observed genetic diversity. Genetic variation exceeded variation in quantitative genetic traits, indicating that selection acted on the population to retain similar traits, which might be because of the multigenomic nature of AMF, where considerable genetic redundancy could buffer the effects of changes in the genetic content of phenotypic traits. These results have direct implications for ecological research and for studying AMF genes, improving commercial AMF inoculum, and understanding evolutionary mechanisms in multigenomic organisms.