Genetic variation in variability: Phenotypic variability of fledging weight and its evolution in a songbird population

Multivariate analysis of herd structure and genetic resource indicators in seedstock beef cattle herds

Genetic, environmental, technological and financial resources are used differently in cattle herds that participate in the same breeding programme. The percentages of calves sired by sires within their own herd or from external herds vary across herds, as do the intensities of use of reproductive biotechnologies. These divergences may be related to differences in the indicators of genetic performance for economic traits. The aim of this study was to determine the factors related to herd structure and genetic resource utilization that exert the greatest influence on the genetic merit of seedstock herds within a Nellore breeding programme. The database comprised 21 factors, along with genomic-enhanced expected progeny differences (GE-EPDs) for growth, reproductive and carcass traits, as well as a selection index of animals from 128 herds. By combining principal component analysis and cluster analysis, we were able to group the herds. We identified statistically significant differences (p < 0.05) in the mean values of the factors, GE-EPDs and genetic trends among the groups of herds. Differences in the percentage of sires from external herds and in sire age between the groups of herds were the factors most associated with differences in mean GE-EPDs and genetic trends. Using young sires from other herds or lineages is an effective strategy in animal breeding. By enhancing genetic variability, this approach does not only improve the genetic quality of herds but also accelerates genetic progress in desired traits over time. Therefore, to ensure the success of this strategy, it is crucial that seedstock herds undergo a thorough selection process aimed at maximizing the genetic potential of future generations of beef cattle.

Effects of species diversity on trait expression of the clonal herb Taraxacum officinale and its relation to genotype diversity and phenotypic plasticity

Plant species respond to varying plant species diversity and associated changes in their abiotic and biotic environment with changes in their phenotype. However, it is not clear to what degree this phenotypic differentiation is due to genotype diversity within populations or phenotypic plasticity of plant individuals. We studied individuals of 16 populations of the clonal herb Taraxacum officinale grown in plant communities of different species richness in a 17-year-old grassland biodiversity experiment (Jena Experiment). We collected 12 individuals in each population to measure phenotypic traits and identify distinct genotypes using microsatellite DNA markers. Plant species richness did not influence population-level genotype and trait diversity. However, it affected the expression of several phenotypic traits, e.g. leaf and inflorescence number, maximum leaf length and seed mass, which increased with increasing plant species richness. Moreover, population-level trait diversity correlated positively with genotype richness for leaf dry matter content (LDMC) and negatively with inflorescence number. For several traits (i.e. seed mass, germination rate, LDMC, specific leaf area (SLA)), a larger portion of variance was explained by genotype identity, while variance in other traits (i.e. number of inflorescences, leaf nitrogen concentration, leaf number, leaf length) resided within genotypes and thus was mostly due to phenotypic plasticity. Overall, our findings show that plant species richness positively affected the population means of some traits related to whole-plant performance, whose variation was achieved through both phenotypic plasticity and genotype composition of a population.

Genomic prediction for grain yield and micro-environmental sensitivity in winter wheat

Individuals within a common environment experience variations due to unique and non-identifiable micro-environmental factors. Genetic sensitivity to micro-environmental variation (i.e. micro-environmental sensitivity) can be identified in residuals, and genotypes with lower micro-environmental sensitivity can show greater resilience towards environmental perturbations. Micro-environmental sensitivity has been studied in animals; however, research on this topic is limited in plants and lacking in wheat. In this article, we aimed to (i) quantify the influence of genetic variation on residual dispersion and the genetic correlation between genetic effects on (expressed) phenotypes and residual dispersion for wheat grain yield using a double hierarchical generalized linear model (DHGLM); and (ii) evaluate the predictive performance of the proposed DHGLM for prediction of additive genetic effects on (expressed) phenotypes and its residual dispersion. Analyses were based on 2,456 advanced breeding lines tested in replicated trials within and across different environments in Denmark and genotyped with a 15K SNP-Illumina-BeadChip. We found that micro-environmental sensitivity for grain yield is heritable, and there is potential for its reduction. The genetic correlation between additive effects on (expressed) phenotypes and dispersion was investigated, and we observed an intermediate correlation. From these results, we concluded that breeding for reduced micro-environmental sensitivity is possible and can be included within breeding objectives without compromising selection for increased yield. The predictive ability and variance inflation for predictions of the DHGLM and a linear mixed model allowing heteroscedasticity of residual variance in different environments (LMM-HET) were evaluated using leave-one-line-out cross-validation. The LMM-HET and DHGLM showed good and similar performance for predicting additive effects on (expressed) phenotypes. In addition, the accuracy of predicting genetic effects on residual dispersion was sufficient to allow genetic selection for resilience. Such findings suggests that DHGLM may be a good choice to increase grain yield and reduce its micro-environmental sensitivity.

Benchmarking in a beef cattle breeding program: Lessons from the best breeders

Beef cattle breeding programs offer genetic evaluations and consulting services on animal breeding practices to help breeders improve the genetic merit of their herds. Some breeders are more willing to apply best practices and technologies than others. Consequently, the average genetic merit and genetic trends differ across herds. We benchmarked some parameters of an average herd (AVE) and the corresponding parameters of herds with higher genetic merit (TOP), both participating in a commercial Nellore breeding program. Expected progeny differences (EPD) for growth, reproductive and carcass traits and a selection index (SI) of animals born from 2005 to 2019 on 128 farms located in Brazil, Bolivia and Paraguay were used to compute the AVE parameters. The 20 herds with higher mean SI of animals born in the last five birth seasons were classified as TOP herds. The mean SI and EPD of animals born in the last five seasons in the TOP herds were, respectively, 89% and 79% to 206% higher (p ≤ 0.001) than those of animals from the AVE herd. Genetic trends over the entire period were also higher (50% for SI and 31% to 88% separately for each trait, p ≤ 0.006) in the TOP herds compared to the AVE herd. Although the difference in the numbers of cows, bulls and calves between the TOP and AVE herds did not reach statistical significance (p = 0.175, p = 0.273 and p = 0.061, respectively), the numbers of progeny per cow and per bull were 21% (p = 0.012) and 26% (p = 0.047) higher in the TOP herds, respectively. Multiple ovulation and embryo transfer and in vitro fertilization and embryo transfer (MOET/IVF) accounted for a higher percentage of births in the TOP herds compared to AVE (24.6% vs. 12.5%, p = 0.002). The generation interval was 17% shorter (p < 0.001) in the TOP herds compared to AVE. The average inbreeding coefficient of animals from the TOP herds (1.08 ± 0.52%) did not differ (p = 0.78) from that of AVE animals (1.26 ± 0.96%). In general, AVE herds are evolving in the desirable direction but differences in genetic merit between AVE and TOP herds are increasing over time. The more frequent use of MOET/IVF, a lower cow-to-bull ratio, and a larger family size (progeny per cow or per bull) can help achieve larger selection differentials and increase genetic trends and average genetic merits of TOP herds compared to AVE herds.

Development of resilience indicator traits based on daily step count data for dairy cattle breeding

Background

Resilient animals are minimally affected by disturbances, such as diseases and heat stress, and quickly recover. Daily activity data can potentially indicate resilience, because resilient animals likely keep variations due to disturbances that threat animal homeostasis at a low magnitude. We used daily step count of cows to define resilience indicators based on theory, exploratory analysis and literature, and then investigated if they can be used to genetically improve resilience by estimating heritability and repeatability, and genetic associations with other resilience-related traits, i.e. health traits, longevity, fertility, and body condition score (BCS).

Results

Two groups of resilience indicators were defined: indicators describing (1) mean step count at different lactation stages for individual cows, and (2) fluctuations in step count from individual step count curves. Heritability estimates were highest for resilience indicators describing mean step count, from 0.22 for the 2-week period pre-partum to 0.45 for the whole lactation. High mean step count was consistently, but weakly, genetically correlated with good health, fertility, and longevity, and high BCS. Heritability estimates of resilience indicators describing fluctuations ranged from 0.01 for number of step count drops to 0.15 for the mean of negative residuals from individual curves. Genetic correlations with health traits, longevity, fertility, and BCS were mostly weak, but were moderate and favorable for autocorrelation of residuals (− 0.33 to − 0.44) and number of step count drops (− 0.44 to − 0.56) with hoof health, fertility, and BCS. Resilience indicators describing variability of residuals and mean of negative residuals showed strong genetic correlations with mean step count (0.86 to 0.95, absolute), which suggests that adjustment for step count level is needed. After adjustment, ‘mean of negative residuals’ was highly genetically correlated with hoof health, fertility, and BCS.

Conclusions

Mean step count, autocorrelation and mean of negative residuals showed most potential as resilience indicators based on resilience theory, heritability, and genetic associations with health, fertility, and body condition score. Other resilience indicators were heritable, but had unfavorable genetic correlations with several health traits. This study is an important first step in the exploration of the use of activity data to breed more resilient livestock.

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Estimation of macro- and micro-genetic environmental sensitivity in unbalanced datasets

Genotype-by-environment interaction is caused by variation in genetic environmental sensitivity (GES), which can be subdivided into macro- and micro-GES. Macro-GES is genetic sensitivity to macro-environments (definable environments often shared by groups of animals), while micro-GES is genetic sensitivity to micro-environments (individual environments). A combined reaction norm and double hierarchical generalised linear model (RN-DHGLM) allows for simultaneous estimation of base genetic, macro- and micro-GES effects. The accuracy of variance components estimated using a RN-DHGLM has been explicitly studied for balanced data and recommendation of a data size with a minimum of 100 sires with at least 100 offspring each have been made. In the current study, the data size (numbers of sires and progeny) and structure requirements of the RN-DHGLM were investigated for two types of unbalanced datasets. Both datasets had a variable number of offspring per sire, but one dataset also had a variable number of offspring within macro-environments. The accuracy and bias of the estimated macro- and micro-GES effects and the estimated breeding values (EBVs) obtained using the RN-DHGLM depended on the data size. Reasonably accurate and unbiased estimates were obtained with data containing 500 sires with 20 offspring or 100 sires with 50 offspring, regardless of the data structure. Variable progeny group sizes, alone or in combination with an unequal number of offspring within macro-environments, had little impact on the dispersion of the EBVs or the bias and accuracy of variance component estimation, but resulted in lower accuracies of the EBVs. Compared to genetic correlations of zero, a genetic correlation of 0.5 between base genetic, macro- and micro-GES components resulted in a slight decrease in the percentage of replicates that converged out of 100 replicates, but had no effect on the dispersion and accuracy of variance component estimation or the dispersion of the EBVs. The results show that it is possible to apply the RN-DHGLM to unbalanced datasets to obtain estimates of variance due to macro- and micro-GES. Furthermore, the levels of accuracy and bias of variance estimates when analysing macro- and micro-GES simultaneously are determined by average family size, with limited impact from variability in family size and/or cohort size. This creates opportunities for the use of field data from populations with unbalanced data structures when estimating macro- and micro-GES.

Exploring the use of residual variance for uniformity of body weight in meat quail lines using Bayesian inference

1. Uniformity in animal products is an important aspect of the production system. Several studies have reported estimates of genetics on residual variance in different species, indicating that it could be exploited to improve uniformity by selection. Nevertheless, there are no reports about the possibilities of such a selection strategy in meat quail.2. Records of hatching weight (HW) and body weight at 42 days (W42) of female and male birds from two meat quail lines (UFV1 and UFV2) were analysed. A three-step genetic evaluation was used to investigate the effect of genetic variation on residual variance of HW and W42 in both lines. In Step 1, a single-trait model was fitted to the data. In Step 2, log-transformed squared estimated residuals (ln(ê²)) were evaluated for these traits. In Step 3, a multi-trait analysis was performed to estimate the genetic correlation between the additive genetic effects for HW, W42, and their respective ln(ê²).3. The heritability estimates ranged from 0.12 to 0.23 for HW and from 0.22 to 0.35 for W42. The estimated heritabilities for the residual part were low and ranged from 0.0003 to 0.02 for both traits, and the genetic coefficient of variation residual variance estimates ranged from 0.31 to 0.42 for HW and from 0.09 to 0.25 for W42. Genetic correlations between the means (HW and W42) and ln(ê²) values were both positive and did not differ from zero, indicating no association between mean and ln(ê²).4. In conclusion, the uniformity of HW and W42 could be improved by selecting for lower residual variance in both meat quail lines, but the accuracy of selection may be low due to low heritability for uniformity, mainly for W42.

Phenotypic noise and the cost of complexity

Experimental studies demonstrate the existence of phenotypic diversity despite constant genotype and environment. Theoretical models based on a single phenotypic character predict that during an adaptation event, phenotypic noise should be positively selected far from the fitness optimum because it increases the fitness of the genotype, and then be selected against when the population reaches the optimum. It is suggested that because of this fitness gain, phenotypic noise should promote adaptive evolution. However, it is unclear how the selective advantage of phenotypic noise is linked to the rate of evolution, and whether any advantage would hold for more realistic, multidimensional phenotypes. Indeed, complex organisms suffer a cost of complexity, where beneficial mutations become rarer as the number of phenotypic characters increases. Using a quantitative genetics approach, we first show that for a one-dimensional phenotype, phenotypic noise promotes adaptive evolution on plateaus of positive fitness, independently from the direct selective advantage on fitness. Second, we show that for multidimensional phenotypes, phenotypic noise evolves to a low-dimensional configuration, with elevated noise in the direction of the fitness optimum. Such a dimensionality reduction of the phenotypic noise promotes adaptive evolution and numerical simulations show that it reduces the cost of complexity.

Patterns of environmental variance across environments and traits in domestic cattle

The variance in phenotypic trait values is a product of environmental and genetic variation. The sensitivity of traits to environmental variation has a genetic component and is likely to be under selection. However, there are few studies investigating the evolution of this sensitivity, in part due to the challenges of estimating the environmental variance. The livestock literature provides a wealth of studies that accurately partition components of phenotypic variance, including the environmental variance, in well-defined environments. These studies involve breeds that have been under strong selection on mean phenotype in optimal environments for many generations, and therefore represent an opportunity to study the potential evolution of trait sensitivity to environmental conditions. Here, we use literature on domestic cattle to examine the evolution of micro-environmental variance (CVR-the coefficient of residual variance) by testing for differences in expression of CVR in animals from the same breed reared in different environments. Traits that have been under strong selection did not follow a null expectation of an increase in CVR in heterogenous environments (e.g., grazing), a pattern that may reflect evolution of increased uniformity in heterogeneous environments. When comparing CVR across environments of different levels of optimality, here measured by trait mean, we found a reduction in CVR in the more optimal environments for both life history and growth traits. Selection aimed at increasing trait means in livestock breeds typically occurs in the more optimal environments, and we therefore suspect that the decreased CVR is a consequence of evolution of the expression of micro-environmental variance in this environment. Our results highlight the heterogeneity in micro-environmental variance across environments and point to possible connections to the intensity of selection on trait means.

Genetic variance for behavioural 'predictability' of stress response

Genetic factors underpinning phenotypic variation are required if natural selection is to result in adaptive evolution. However, evolutionary and behavioural ecologists typically focus on variation among individuals in their average trait values and seek to characterize genetic contributions to this. As a result, less attention has been paid to if and how genes could contribute towards within-individual variance or trait 'predictability'. In fact, phenotypic 'predictability' can vary among individuals, and emerging evidence from livestock genetics suggests this can be due to genetic factors. Here, we test this empirically using repeated measures of a behavioural stress response trait in a pedigreed population of wild-type guppies. We ask (a) whether individuals differ in behavioural predictability and (b) whether this variation is heritable and so evolvable under selection. Using statistical methodology from the field of quantitative genetics, we find support for both hypotheses and also show evidence of a genetic correlation structure between the behavioural trait mean and individual predictability. We show that investigating sources of variability in trait predictability is statistically tractable and can yield useful biological interpretation. We conclude that, if widespread, genetic variance for 'predictability' will have major implications for the evolutionary causes and consequences of phenotypic variation.

The Evolution of Variance Control

Genetically identical individuals can be phenotypically variable, even in constant environmental conditions. The ubiquity of this phenomenon, known as 'intra-genotypic variability', is increasingly evident and the relevant mechanistic underpinnings are beginning to be understood. In parallel, theory has delineated a number of formal expectations for contexts in which such a feature would be adaptive. Here, we review empirical evidence across biological systems and theoretical expectations, including nonlinear averaging and bet hedging. We synthesize existing results to illustrate the dependence of selection outcomes both on trait characteristics, features of environmental variability, and species' demographic context. We conclude by discussing ways to bridge the gap between empirical evidence of intra-genotypic variability, studies demonstrating its genetic component, and evidence that it is adaptive.

Genetics and genomics of uniformity and resilience in livestock and aquaculture species: A review

Genetic control of residual variance offers opportunities to increase uniformity and resilience of livestock and aquaculture species. Improving uniformity and resilience of animals will improve health and welfare of animals and lead to more homogenous products. Our aims in this review were to summarize the current models and methods to study genetic control of residual variance, genetic parameters and genomic results for residual variance and discuss future research directions. Typically, the genetic coefficient of variation is high (median = 0.27; range 0-0.86) and the heritability of residual variance is low (median = 0.01; range 0-0.10). Higher heritabilities can be achieved when increasing the number of records per animal. Divergent selection experiments have supported the feasibility of selecting for high or low residual variance. Genomic studies have revealed associations in regions related to stress, including those from the heat shock protein family. Although the number of studies is growing, genetic control of residual variance is still poorly understood, but big data and genomics offer great opportunities.

Environmental coupling of heritability and selection is rare and of minor evolutionary significance in wild populations

Predicting the rate of adaptation to environmental change in wild populations is important for understanding evolutionary change. However, predictions may be unreliable if the two key variables affecting the rate of evolutionary change-heritability and selection-are both affected by the same environmental variable. To determine how general such an environmentally induced coupling of heritability and selection is, and how this may influence the rate of adaptation, we made use of freely accessible, open data on pedigreed wild populations to answer this question at the broadest possible scale. Using 16 populations from 10 vertebrate species, which provided data on 50 traits (relating to body mass, morphology, physiology, behaviour and life history), we found evidence for an environmentally induced relationship between heritability and selection in only 6 cases, with weak evidence that this resulted in an increase or decrease in the expected selection response. We conclude that such a coupling of heritability and selection is unlikely to strongly affect evolutionary change, even though both heritability and selection are commonly postulated to be dependent on the environment.

Modelling the co-evolution of indirect genetic effects and inherited variability

When individuals interact, their phenotypes may be affected not only by their own genes but also by genes in their social partners. This phenomenon is known as Indirect Genetic Effects (IGEs). In aquaculture species and some plants, however, competition not only affects trait levels of individuals, but also inflates variability of trait values among individuals. In the field of quantitative genetics, the variability of trait values has been studied as a quantitative trait in itself, and is often referred to as inherited variability. Such studies, however, consider only the genetic effect of the focal individual on trait variability and do not make a connection to competition. Although the observed phenotypic relationship between competition and variability suggests an underlying genetic relationship, the current quantitative genetic models of IGE and inherited variability do not allow for such a relationship. The lack of quantitative genetic models that connect IGEs to inherited variability limits our understanding of the potential of variability to respond to selection, both in nature and agriculture. Models of trait levels, for example, show that IGEs may considerably change heritable variation in trait values. Currently, we lack the tools to investigate whether this result extends to variability of trait values. Here we present a model that integrates IGEs and inherited variability. In this model, the target phenotype, say growth rate, is a function of the genetic and environmental effects of the focal individual and of the difference in trait value between the social partner and the focal individual, multiplied by a regression coefficient. The regression coefficient is a genetic trait, which is a measure of cooperation; a negative value indicates competition, a positive value cooperation, and an increasing value due to selection indicates the evolution of cooperation. In contrast to the existing quantitative genetic models, our model allows for co-evolution of IGEs and variability, as the regression coefficient can respond to selection. Our simulations show that the model results in increased variability of body weight with increasing competition. When competition decreases, i.e., cooperation evolves, variability becomes significantly smaller. Hence, our model facilitates quantitative genetic studies on the relationship between IGEs and inherited variability. Moreover, our findings suggest that we may have been overlooking an entire level of genetic variation in variability, the one due to IGEs.

Maternal Effects in a Wild Songbird Are Environmentally Plastic but Only Marginally Alter the Rate of Adaptation

Despite ample evidence for the presence of maternal effects (MEs) in a variety of traits and strong theoretical indications for their evolutionary consequences, empirical evidence to what extent MEs can influence evolutionary responses to selection remains ambiguous. We tested the degree to which MEs can alter the rate of adaptation of a key life-history trait, clutch size, using an individual-based model approach parameterized with experimental data from a long-term study of great tits (Parus major). We modeled two types of MEs: (i) an environmentally plastic ME, in which the relationship between maternal and offspring clutch size depended on the maternal environment via offspring condition, and (ii) a fixed ME, in which this relationship was constant. Although both types of ME affected the rate of adaptation following an abrupt environmental shift, the overall effects were small. We conclude that evolutionary consequences of MEs are modest at best in our study system, at least for the trait and the particular type of ME we considered here. A closer link between theoretical and empirical work on MEs would hence be useful to obtain accurate predictions about the evolutionary consequences of MEs more generally.

The effects of scaling on age, sex and size relationships in Red-legged Partridges

Wild birds differ in size according to their age and sex, adult birds being larger than juveniles. In the galliforms, males are larger than females, in contrast to some groups, such as the raptors, in which the females are larger. Size generally influences the rank hierarchy within a group of birds, although the age, sex, temperament and behaviour of an individual may override its size related rank order. The scaled size of birds according to age and sex affects their physiology and behaviour. Precise details of body-size differences by age and sex are poorly known in most partridge species. We measured 13,814 wild partridges in a homogenous population over 14 years of study to evaluate size differences within a uniform habitat and population management regime. We show that wild Red-legged Partridges have scaled mass, and body- and wing-lengths consistent with age/sex classes. Power functions between mass and body-length (as a proxy for walking efficiency), and between mass and wing-length (for flight efficiency) differ between juvenile females and males, and adult females and males. We discuss these findings and their physiological, behavioural and ecological implications.

Genetics of Genome-Wide Recombination Rate Evolution in Mice from an Isolated Island

Recombination rate is a heritable quantitative trait that evolves despite the fundamentally conserved role that recombination plays in meiosis. Differences in recombination rate can alter the landscape of the genome and the genetic diversity of populations. Yet our understanding of the genetic basis of recombination rate evolution in nature remains limited. We used wild house mice (Mus musculus domesticus) from Gough Island (GI), which diverged recently from their mainland counterparts, to characterize the genetics of recombination rate evolution. We quantified genome-wide autosomal recombination rates by immunofluorescence cytology in spermatocytes from 240 F2 males generated from intercrosses between GI-derived mice and the wild-derived inbred strain WSB/EiJ. We identified four quantitative trait loci (QTL) responsible for inter-F2 variation in this trait, the strongest of which had effects that opposed the direction of the parental trait differences. Candidate genes and mutations for these QTL were identified by overlapping the detected intervals with whole-genome sequencing data and publicly available transcriptomic profiles from spermatocytes. Combined with existing studies, our findings suggest that genome-wide recombination rate divergence is not directional and its evolution within and between subspecies proceeds from distinct genetic loci.

Heritable Micro-environmental Variance Covaries with Fitness in an Outbred Population of Drosophila serrata

The genetic basis of stochastic variation within a defined environment, and the consequences of such micro-environmental variance for fitness are poorly understood . Using a multigenerational breeding design in Drosophila serrata, we demonstrated that the micro-environmental variance in a set of morphological wing traits in a randomly mating population had significant additive genetic variance in most single wing traits. Although heritability was generally low (<1%), coefficients of additive genetic variance were of a magnitude typical of other morphological traits, indicating that the micro-environmental variance is an evolvable trait. Multivariate analyses demonstrated that the micro-environmental variance in wings was genetically correlated among single traits, indicating that common mechanisms of environmental buffering exist for this functionally related set of traits. In addition, through the dominance genetic covariance between the major axes of micro-environmental variance and fitness, we demonstrated that micro-environmental variance shares a genetic basis with fitness, and that the pattern of selection is suggestive of variance-reducing selection acting on micro-environmental variance.

Genetic basis of between-individual and within-individual variance of docility

Between-individual variation in phenotypes within a population is the basis of evolution. However, evolutionary and behavioural ecologists have mainly focused on estimating between-individual variance in mean trait and neglected variation in within-individual variance, or predictability of a trait. In fact, an important assumption of mixed-effects models used to estimate between-individual variance in mean traits is that within-individual residual variance (predictability) is identical across individuals. Individual heterogeneity in the predictability of behaviours is a potentially important effect but rarely estimated and accounted for. We used 11 389 measures of docility behaviour from 1576 yellow-bellied marmots (Marmota flaviventris) to estimate between-individual variation in both mean docility and its predictability. We then implemented a double hierarchical animal model to decompose the variances of both mean trait and predictability into their environmental and genetic components. We found that individuals differed both in their docility and in their predictability of docility with a negative phenotypic covariance. We also found significant genetic variance for both mean docility and its predictability but no genetic covariance between the two. This analysis is one of the first to estimate the genetic basis of both mean trait and within-individual variance in a wild population. Our results indicate that equal within-individual variance should not be assumed. We demonstrate the evolutionary importance of the variation in the predictability of docility and illustrate potential bias in models ignoring variation in predictability. We conclude that the variability in the predictability of a trait should not be ignored, and present a coherent approach for its quantification.