Partitioning of genetic variation across the genome using multimarker methods in a wild bird population

The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.     

Partial genotyping at polymorphic markers can improve heritability estimates in sibling groups

Accurate estimates of heritability () are necessary to assess adaptive responses of populations and evolution of fitness‐related traits in changing environments. For plants, estimates generally rely on maternal progeny designs, assuming that offspring are either half‐sibs or unrelated. However, plant mating systems often depart from half‐sib assumptions, this can bias estimates. Here, we investigate how to accurately estimate in nonmodel species through the analysis of sibling designs with a moderate genotyping effort. We performed simulations to investigate how microsatellite marker information available for only a subset of offspring can improve estimates based on maternal progeny designs in the presence of nonrandom mating, inbreeding in the parental population or maternal effects. We compared the basic family method, considering or not adjustments based on average relatedness coefficients, and methods based on the animal model. The animal model was used with average relatedness information, or with hybrid relatedness information: associating one‐generation pedigree and family assumptions, or associating one‐generation pedigree and average relatedness coefficients. Our results highlighted that methods using marker‐based relatedness coefficients performed as well as pedigree‐based methods in the presence of nonrandom mating (i.e. unequal male reproductive contributions, selfing), offering promising prospects to investigate in situ heritabilities in natural populations. In the presence of maternal effects, only the use of pairwise relatednesses through pedigree information improved the accuracy of estimates. In that case, the amount of father‐related offspring in the sibling design is the most critical. Overall, we showed that the method using both one‐generation pedigree and average relatedness coefficients was the most robust to various ecological scenarios.     

Resource quality affects weapon and testis size and the ability of these traits to respond to selection in the leaf‐footed cactus bug,Narnia femorata

The size of weapons and testes can be central to male reproductive success. Yet, the expression of these traits is often extremely variable. Studies are needed that take a more complete organism perspective, investigating the sources of variation in both traits simultaneously and using developmental conditions that mimic those in nature. In this study, we investigated the components of variation in weapon and testis sizes using the leaf‐footed cactus bug, Narnia femorata (Hemiptera: Coreidae) on three natural developmental diets. We show that the developmental diet has profound effects on both weapon and testis expression and scaling. Intriguingly, males in the medium‐quality diet express large weapons but have relatively tiny testes, suggesting complex allocation decisions. We also find that heritability, evolvability, and additive genetic variation are highest in the high‐quality diet for testis and body mass. This result suggests that these traits may have an enhanced ability to respond to selection during a small window of time each year when this diet is available. Taken together, these results illustrate that normal, seasonal fluctuations in the nutritional environment may play a large role in the expression of sexually selected traits and the ability of these traits to respond to selection.     

Similar recent selection criteria associated with different behavioural effects in two dog breeds

Selection during the last decades has split some established dog breeds into morphologically and behaviourally divergent types. These breed splits are interesting models for behaviour genetics since selection has often been for few and well‐defined behavioural traits. The aim of this study was to explore behavioural differences between selection lines in golden and Labrador retriever, in both of which a split between a common type (pet and conformation) and a field type (hunting) has occurred. We hypothesized that the behavioural profiles of the types would be similar in both breeds. Pedigree data and results from a standardized behavioural test from 902 goldens (698 common and 204 field) and 1672 Labradors (1023 and 649) were analysed. Principal component analysis revealed six behavioural components: curiosity, play interest, chase proneness, social curiosity, social greeting and threat display. Breed and type affected all components, but interestingly there was an interaction between breed and type for most components. For example, in Labradors the common type had higher curiosity than the field type (F_1,1668 = 18.359; P < 0.001), while the opposite was found in goldens (F_1,897 = 65.201; P < 0.001). Heritability estimates showed considerable genetic contributions to the behavioural variations in both breeds, but different heritabilities between the types within breeds was also found, suggesting different selection pressures. In conclusion, in spite of similar genetic origin and similar recent selection criteria, types behave differently in the breeds. This suggests that the genetic architecture related to behaviour differs between the breeds.     

Evolution,plasticity and evolving plasticity of phenology in the tree species Alnus glutinosa

Both traits and the plasticity of these traits are subject to evolutionary change and therefore affect the long‐term persistence of populations and their role in local communities. We subjected clones from 12 different populations of Alnus glutinosa, located along a latitudinal gradient, to two different temperature treatments, to disentangle the distribution of genetic variation in timing of bud burst and bud burst plasticity within and among genotypes, populations, and regions. We calculated heritability and evolvability estimates for bud burst and bud burst plasticity and assessed the influence of divergent selection relative to neutral drift. We observed higher levels of heritability and evolvability for bud burst than for its plasticity, whereas the total phenological heritability and evolvability (i.e. combining timing of bud burst and bud burst plasticity) suggest substantial evolutionary potential with respect to phenology. Earlier bud burst was observed for the low‐latitudinal populations than for the populations from higher latitudes, whereas the high‐latitudinal populations did not show the expected delayed bud burst. This countergradient variation can be due to evolution towards increased phenological plasticity at higher latitudes. However, because we found little evidence for adaptive differences in phenological plasticity across the latitudinal gradient, we suggest differential frost tolerance as the most likely explanation for the observed phenological patterns in A. glutinosa.     

Multivariate selection and intersexual genetic constraints in a wild bird population

When selection differs between the sexes for traits that are genetically correlated between the sexes, there is potential for the effect of selection in one sex to be altered by indirect selection in the other sex, a situation commonly referred to as intralocus sexual conflict (ISC). While potentially common, ISC has rarely been studied in wild populations. Here, we studied ISC over a set of morphological traits (wing length, tarsus length, bill depth and bill length) in a wild population of great tits (Parus major) from Wytham Woods, UK. Specifically, we quantified the microevolutionary impacts of ISC by combining intra‐ and intersex additive genetic (co)variances and sex‐specific selection estimates in a multivariate framework. Large genetic correlations between homologous male and female traits combined with evidence for sex‐specific multivariate survival selection suggested that ISC could play an appreciable role in the evolution of this population. Together, multivariate sex‐specific selection and additive genetic (co)variance for the traits considered accounted for additive genetic variance in fitness that was uncorrelated between the sexes (cross‐sex genetic correlation = ?0.003, 95% CI = ?0.83, 0.83). Gender load, defined as the reduction in a population's rate of adaptation due to sex‐specific effects, was estimated at 50% (95% CI = 13%, 86%). This study provides novel insights into the evolution of sexual dimorphism in wild populations and illustrates how quantitative genetics and selection analyses can be combined in a multivariate framework to quantify the microevolutionary impacts of ISC.     

Genetic (co)variances between milk mineral concentration and chemical composition in lactating Holstein-Friesian dairy cows

Milk mineral concentration is important from both the perspective of processing milk into dairy products and its nutritive value for human consumption. Precise estimates of genetic parameters for milk mineral concentration are lacking because of the considerable resources required to collect vast phenotypes quantities. The milk concentration of calcium (Ca), potassium (K), magnesium (Mg), sodium (Na) and phosphorus (P) in the present study was quantified from mid-IR spectroscopy on 12 223 test-day records from 1717 Holstein-Friesian cows. (Co)variance components were estimated using random regressions to model both the additive genetic and within-lactation permanent environmental variances of each trait. The coefficient of genetic variation averaged across days-in-milk (DIM) was 6.93%, 3.46%, 6.55%, 5.20% and 6.68% for Ca, K, Mg, Na and P concentration, respectively; heritability estimates varied across lactation from 0.31±0.05 (5 DIM) to 0.67±0.04 (181 DIM) for Ca, from 0.18±0.03 (60 DIM) to 0.24±0.05 (305 DIM) for K, from 0.08±0.03 (15 DIM) to 0.37±0.03 (223 DIM) for Mg, from 0.16±0.03 (30 DIM) to 0.37±0.04 (305 DIM) for Na and from 0.21±0.04 (12 DIM) to 0.57±0.04 (211 DIM) for P. Genetic correlations within the same trait across different DIM were almost unity between adjacent DIM but weakened as the time interval between pairwise compared DIM lengthened; genetic correlations were weaker than 0.80 only when comparing both peripheries of the lactation. The analysis of the geometry of the additive genetic covariance matrix revealed that almost 90% of the additive genetic variation was accounted by the intercept term of the covariance functions for each trait. Milk protein concentration and mineral concentration were, in general, positively genetically correlated with each other across DIM, whereas milk fat concentration was positively genetically correlated throughout the entire lactation with Ca, K and Mg; the genetic correlation with fat concentration changed from negative to positive with Na and P at 243 DIM and 50 DIM, respectively. Genetic correlations between somatic cell score and Na ranged from 0.38±0.21 (5 DIM) to 0.79±0.18 (305 DIM). Exploitable genetic variation existed for all milk minerals, although many national breeding objectives are probably contributing to an indirect positive response to selection in milk mineral concentration.     

Direct and indirect genetic effects on reproductive investment in a grasshopper

A fundamental part of the quantitative genetic theory deals with the partitioning of the phenotypic variance into additive genetic and environmental components. During interaction with conspecifics, the interaction partner becomes a part of the environment from the perspective of the focal individual. If the interaction effects have a genetic basis, they are called indirect genetic effects (IGEs) and can evolve along with direct genetic effects. Sexual reproduction is a classic context where potential conflict between males and females can arise from trade‐offs between current and future investments. We studied five female fecundity traits, egg length and number, egg pod length and number and latency to first egg pod, and estimated the direct and IGEs using a half‐sib breeding design in the grasshopper Chorthippus biguttulus. We found that the male IGEs were an order of magnitude lower than the direct genetic effects and were not significantly different from zero. However, there was some indication that IGEs were larger shortly after mating, consistent with the idea that IGEs fade with time after interaction. Female direct heritabilities were moderate to low. Simulation shows that the variance component estimates can appear larger with less data, calling for care when interpreting variance components estimated with low power. Our results illustrate that the contribution of male IGEs is overall low on the phenotypic variance of female fecundity traits. Thus, even in the relevant context of sexual conflict, the influence of male IGEs on the evolutionary trajectory of female reproductive traits is likely to be small.     

WHAT DO WE KNOW ABOUT THE HERITABILITY OF DEVELOPMENTAL INSTABILITY? ANSWERS FROM A BAYESIAN MODEL

In spite of over half a century of research, little is known about the genetic basis of developmental instability (DI). The estimation of the heritability of DI is seriously hampered by the fact that fluctuating asymmetry-FA, that is, the observable outcome of DI-only poorly reflects DI. This results in an underestimation of the heritability of DI. Current methods transforming heritabilities in FA into those of DI fail to take all sources of variation into account and yield incorrect confidence bands that are usually based on unrealistic assumptions. Therefore, a Bayesian latent variable model is developed and explored. Simulations show that with sample sizes currently applied in empirical studies, extremely wide posterior distributions are obtained and data do not allow to distinguish between high (0.5) and low (0.1) heritabilities of DI at all. Even sample sizes of 5000 result in very wide posteriors in many cases. Furthermore, for smaller samples (250 and 1250), up to 70% of the estimates of the heritability of DI were below the mean expected value because of the high skewness of its distribution. Knowing that in only one study, sample sizes were above 5000, there is a need for larger studies to evaluate the evolutionary potential of DI. Designs with relatively low numbers of sires (1-2% of total number of offspring) appear most efficient. Because such high sample sizes are hard to obtain for many study organisms, more insights are required about how data from different traits can be combined in a single analysis. In addition, new designs and methods, such as QTL analyses and micro-array techniques, should be applied to gain a better understanding of the genetic basis of DI.     

Evolutionary analysis of a key floral trait in aquilegia canadensis (ranunculaceae): genetic variation in herkogamy and its effect on the mating system

The mating system of flowering plant populations evolves through selection on genetically based phenotypic variation in floral traits. The physical separation of anthers and stigmas within flowers (herkogamy) is expected to be an important target of selection to limit self-fertilization. We investigated the pattern of phenotypic and genetic variation in herkogamy and its effect of self-fertilization in a broad sample of natural populations of Aquilegia canadensis, a species that is highly selfing despite strong inbreeding depression. Within natural populations, plants exhibit substantial phenotypic variation in herkogamy caused primarily by variation in pistil length rather than stamen length. Compared to other floral traits, herkogamy is much more variable and a greater proportion of variation is distributed among rather than within individuals. We tested for a genetic component of this marked phenotypic variation by growing naturally pollinated seed families from five populations in a common greenhouse environment. For three populations, we detected a significant variation in herkogamy among families, and a positive regression between parental herkogamy measured in the field and progeny herkogamy in the greenhouse, suggesting that there is often genetic variation in herkogamy within natural populations. We estimated levels of self-fertilization for groups of flowers that differed in herkogamy and show that, as expected, herkogamy was associated with reduced selfing in 13 of 19 populations. In six of these populations, we performed floral emasculations to show that this decrease in selfing is due to decreased autogamy (within-flower selfing), the mode of selfing that herkogamy should most directly influence. Taken together, these results suggest that increased herkogamy should be selected to reduce the production of low-quality selfed seed. The combination of high selfing and substantial genetic variation for herkogamy in A. canadensis is enigmatic, and reconciling this observation will require a more integrated analysis of how herkogamy influences not only self-fertilization, but also patterns of outcross pollen import and export.