Genetic analysis of environmental variation

Statistical tests for detecting variance effects in quantitative trait studies

Motivation

Identifying variants, both discrete and continuous, that are associated with quantitative traits, or QTs, is the primary focus of quantitative genetics. Most current methods are limited to identifying mean effects, or associations between genotype or covariates and the mean value of a quantitative trait. It is possible, however, that a variant may affect the variance of the quantitative trait in lieu of, or in addition to, affecting the trait mean. Here, we develop a general methodology to identify covariates with variance effects on a quantitative trait using a Bayesian heteroskedastic linear regression model (BTH). We compare BTH with existing methods to detect variance effects across a large range of simulations drawn from scenarios common to the analysis of quantitative traits.

Results

We find that BTH and a double generalized linear model (dglm) outperform classical tests used for detecting variance effects in recent genomic studies. We show BTH and dglm are less likely to generate spurious discoveries through simulations and application to identifying methylation variance QTs and expression variance QTs. We identify four variance effects of sex in the Cardiovascular and Pharmacogenetics study. Our work is the first to offer a comprehensive view of variance identifying methodology. We identify shortcomings in previously used methodology and provide a more conservative and robust alternative. We extend variance effect analysis to a wide array of covariates that enables a new statistical dimension in the study of sex and age specific quantitative trait effects.

Availability and implementation

https://github.com/b2du/bth.

Supplementary information

Supplementary data are available at Bioinformatics online.

Bayesian, Likelihood-Free Modelling of Phenotypic Plasticity and Variability in Individuals and Populations

There is a paradigm shift from the traditional focus on the "average" individual towards the definition and analysis of trait variation within individual life-history and among individuals in populations. This is a result of increasing availability of individual phenotypic data. The shift allows the use of genetic and environment-driven variations to assess robustness to challenge, gain greater understanding of organismal biological processes, or deliver individual-targeted treatments or genetic selection. These consequences apply, in particular, to variation in ontogenetic growth. We propose an approach to parameterise mathematical models of individual traits (e.g., reaction norms, growth curves) that address two challenges: 1) Estimation of individual traits while making minimal assumptions about data distribution and correlation, addressed via Approximate Bayesian Computation (a form of nonparametric inference). We are motivated by the fact that available information on distribution of biological data is often less precise than assumed by conventional likelihood functions. 2) Scaling-up to population phenotype distributions while facilitating unbiased use of individual data; this is addressed via a probabilistic framework where population distributions build on separately-inferred individual distributions and individual-trait interpretability is preserved. The approach is tested against Bayesian likelihood-based inference, by fitting weight and energy intake growth models to animal data and normal- and skewed-distributed simulated data. i) Individual inferences were accurate and robust to changes in data distribution and sample size; in particular, median-based predictions were more robust than maximum- likelihood-based curves. These results suggest that the approach gives reliable inferences using few observations and monitoring resources. ii) At the population level, each individual contributed via a specific data distribution, and population phenotype estimates were not disproportionally influenced by outlier individuals. Indices measuring population phenotype variation can be derived for study comparisons. The approach offers an alternative for estimating trait variability in biological systems that may be reliable for various applications, for example, in genetics, health, and individualised nutrition, while using fewer assumptions and fewer empirical observations. In livestock breeding, the potentially greater accuracy of trait estimation (without specification of multitrait variance-covariance parameters) could lead to improved selection and to more decisive estimates of trait heritability.

The Statistical Scale Effect as a Source of Positive Genetic Correlation Between Mean and Variability: A Simulation Study

The selection objective for animal production is the highest income with the lowest production cost, while ensuring the highest animal welfare. A selection experiment for environmental variability of birth weight in mice showed a correlated response in the mean after 20 generations starting from a crossed panmictic population. The relationship between the birth weight and its environmental variability explained the correlated response. The scale effect represents a potential cause of this correlation. The relationship between the mean and the variability implies: the higher the mean, the higher the variability. The study was to quantify by simulation the genetic correlation between a trait and its environmental variability. This can be attributable to the scale effect in a range of coefficients of variation and heritabilities between 0.05 and 0.50. The resulting genetic correlation ranged from 0.1335 to 0.7021 being the highest for the highest heritability and the lowest CV. The scale effect for a trait with heritability between 0.25 and 0.35 and CV between 0.15 and 0.25 generated a genetic correlation between 0.43 and 0.57. The genetic coefficient of variation (GCV) affecting residual variability was modulated by the strength reducing the impact of the scale effect. GCV ranged from 0.0050 to 1.4984. The strength of the scale effect might be in the range between 0 and 1. The scale effect would explain many reported genetic correlation and the additive genetic variance for the variability. This is relevant when increasing the mean of a trait jointly with the reduction of its variability.

Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank

Genotype-by-environment interaction (GEI) is a fundamental component in understanding complex trait variation. However, it remains challenging to identify genetic variants with GEI effects in humans largely because of the small effect sizes and the difficulty of monitoring environmental fluctuations. Here, we demonstrate that GEI can be inferred from genetic variants associated with phenotypic variability in a large sample without the need of measuring environmental factors. We performed a genome-wide variance quantitative trait locus (vQTL) analysis of ~5.6 million variants on 348,501 unrelated individuals of European ancestry for 13 quantitative traits in the UK Biobank and identified 75 significant vQTLs with P < 2.0 × 10-9 for 9 traits, especially for those related to obesity. Direct GEI analysis with five environmental factors showed that the vQTLs were strongly enriched with GEI effects. Our results indicate pervasive GEI effects for obesity-related traits and demonstrate the detection of GEI without environmental data.

Genetic Effects on Dispersion in Urinary Albumin and Creatinine in Three House Mouse (Mus musculus) Cohorts

Conventionally, quantitative genetics concerns the heredity of trait means, but there is growing evidence for the existence of architectures in which certain alleles cause random variance in phenotype, termed ‘phenotypic dispersion’ (PD) or ‘variance QTL’ (vQTL), including in physiological traits like disease signs. However, the structure of this phenomenon is still poorly known. PD for urinary albumin (PD_UAlb) and creatinine (PD_UCrea) was mapped using curated data from two nearly genetically identical F₂ mouse (Mus musculus) cohorts (383 male F₂ C57BL/6J×A/J (97 SNP) and 207 male F₂ C57BL/6J×A/J ApoE knockout mice (144 SNP)) and a related mapping cohort (340 male F₂ DBA/2J×C57BL/6J (83 SNP, 8 microsatellites)). PD_UAlb was associated with markers in regions of Chr 1 (5-64 megabases (MB); 141-158 MB), 3 (∼113 MB), 8 (37-68 MB), 14 (92-117 MB) and 17 (14-24 MB) with several positions and quantitative architectures in common between the two C57BL/6J×A/J cohorts, most of which had a negative dominant construction. One locus for PD_UCrea was detected on Chr 19 (57 MB) in the C57BL/6J×A/J ApoE^−/− cohort. The large number of negative dominant loci for albuminuria dispersion relative to conventional quantitative trait loci suggests that the development of albuminuria may be largely genetically dynamic and that randomization in this development is detrimental.

Opportunities to Improve Resilience in Animal Breeding Programs

Resilience is the capacity of an animal to be minimally affected by disturbances or to rapidly return to the state pertained before exposure to a disturbance. However, indicators for general resilience to environmental disturbances have not yet been defined, and perhaps therefore resilience is not yet included in breeding goals. The current developments on big data collection give opportunities to determine new resilience indicators based on longitudinal data, which can aid to incorporate resilience in animal breeding goals. The objectives of this paper were: (1) to define resilience indicator traits based on big data, (2) to define economic values for resilience, and (3) to show the potential to improve resilience of livestock through inclusion of resilience in breeding goals. Resilience might be measured based on deviations from expected production levels over a period of time. Suitable resilience indicators could be the variance of deviations, the autocorrelation of deviations, the skewness of deviations, and the slope of a reaction norm. These (new) resilience indicators provide opportunity to include resilience in breeding programs. Economic values of resilience indicators in the selection index can be calculated based on reduced costs due to labor and treatments. For example, when labor time is restricted, the economic value of resilience increases with an increasing number of animals per farm, and can become as large as the economic value of production. This shows the importance of including resilience in breeding goals. Two scenarios were described to show the additional benefit of including resilience in breeding programs. These examples showed that it is hard to improve resilience with only production traits in the selection index, but that it is possible to greatly improve resilience by including resilience indicators in the selection index. However, when health-related traits are already present in the selection index, the effect is smaller. Nevertheless, inclusion of resilience indicators in the selection index increases the response in the breeding goal and resilience, which results in less labor-demanding, and thus easier-to-manage livestock. Current developments on massive collection of data, and new phenotypes based on these data, offer exciting opportunities to breed for improved resilience of livestock.

Analysis of the impact of cytoplasmic and mitochondrial inheritance on litter size and carcass in rabbits

<p>The effects of mitogenome variation on economically important traits have been reported in a number of domestic animal species. In this study, the first of its kind on rabbits, we have performed the estimation of the contribution of cytoplasmic and D-loop mitochondrial DNA (mtDNA) sequence effects on the litter size and carcass traits in three Pannon rabbit breeds (Pannon Ka, Pannon Large and Pannon White). The observed effects of both estimates, coming from cytoplasmic or D-loop mtDNA variation, were negligible. The most likely explanation for the results obtained is the lack of mitogenome polymorphism in all three populations, as suggested from the analysis performed on the D-loop mtDNA sequence, here assigned to the two most frequent rabbit haplotypes. The extent of potential benefits of the introduction, or alteration, of mitogenome variation in rabbit breeding remains an open question for future research.</p>

vqtl: An R Package for Mean-Variance QTL Mapping

We present vqtl, an R package for mean-variance QTL mapping. This QTL mapping approach tests for genetic loci that influence the mean of the phenotype, termed mean QTL, the variance of the phenotype, termed variance QTL, or some combination of the two, termed mean-variance QTL. It is unique in its ability to correct for variance heterogeneity arising not only from the QTL itself but also from nuisance factors, such as sex, batch, or housing. This package provides functions to conduct genome scans, run permutations to assess the statistical significance, and make informative plots to communicate results. Because it is inter-operable with the popular qtl package and uses many of the same data structures and input patterns, it will be straightforward for geneticists to analyze future experiments with vqtl as well as re-analyze past experiments, possibly discovering new QTL.

QTL Mapping on a Background of Variance Heterogeneity

Standard QTL mapping procedures seek to identify genetic loci affecting the phenotypic mean while assuming that all individuals have the same residual variance. But when the residual variance differs systematically between groups, perhaps due to a genetic or environmental factor, such standard procedures can falter: in testing for QTL associations, they attribute too much weight to observations that are noisy and too little to those that are precise, resulting in reduced power and and increased susceptibility to false positives. The negative effects of such "background variance heterogeneity" (BVH) on standard QTL mapping have received little attention until now, although the subject is closely related to work on the detection of variance-controlling genes. Here we use simulation to examine how BVH affects power and false positive rate for detecting QTL affecting the mean (mQTL), the variance (vQTL), or both (mvQTL). We compare linear regression for mQTL and Levene's test for vQTL, with tests more recently developed, including tests based on the double generalized linear model (DGLM), which can model BVH explicitly. We show that, when used in conjunction with a suitable permutation procedure, the DGLM-based tests accurately control false positive rate and are more powerful than the other tests. We also find that some adverse effects of BVH can be mitigated by applying a rank inverse normal transform. We apply our novel approach, which we term "mean-variance QTL mapping", to publicly available data on a mouse backcross and, after accommodating BVH driven by sire, detect a new mQTL for bodyweight.

QTL Mapping on a Background of Variance Heterogeneity

Standard QTL mapping procedures seek to identify genetic loci affecting the phenotypic mean while assuming that all individuals have the same residual variance. But when the residual variance differs systematically between groups, perhaps due to a genetic or environmental factor, such standard procedures can falter: in testing for QTL associations, they attribute too much weight to observations that are noisy and too little to those that are precise, resulting in reduced power and and increased susceptibility to false positives. The negative effects of such "background variance heterogeneity" (BVH) on standard QTL mapping have received little attention until now, although the subject is closely related to work on the detection of variance-controlling genes. Here we use simulation to examine how BVH affects power and false positive rate for detecting QTL affecting the mean (mQTL), the variance (vQTL), or both (mvQTL). We compare linear regression for mQTL and Levene's test for vQTL, with tests more recently developed, including tests based on the double generalized linear model (DGLM), which can model BVH explicitly. We show that, when used in conjunction with a suitable permutation procedure, the DGLM-based tests accurately control false positive rate and are more powerful than the other tests. We also find that some adverse effects of BVH can be mitigated by applying a rank inverse normal transform. We apply our novel approach, which we term "mean-variance QTL mapping", to publicly available data on a mouse backcross and, after accommodating BVH driven by sire, detect a new mQTL for bodyweight.

Identifying loci affecting trait variability and detecting interactions in genome-wide association studies

Identification of genetic variants with effects on trait variability can provide insights into the biological mechanisms that control variation and can identify potential interactions. We propose a two-degree-of-freedom test for jointly testing mean and variance effects to identify such variants. We implement the test in a linear mixed model, for which we provide an efficient algorithm and software. To focus on biologically interesting settings, we develop a test for dispersion effects, that is, variance effects not driven solely by mean effects when the trait distribution is non-normal. We apply our approach to body mass index in the subsample of the UK Biobank population with British ancestry (n ~408,000) and show that our approach can increase the power to detect associated loci. We identify and replicate novel associations with significant variance effects that cannot be explained by the non-normality of body mass index, and we provide suggestive evidence for a connection between leptin levels and body mass index variability.

Genomic regions underlying uniformity of yearling weight in Nellore cattle evaluated under different response variables

BACKGROUND

In livestock, residual variance has been studied because of the interest to improve uniformity of production. Several studies have provided evidence that residual variance is partially under genetic control; however, few investigations have elucidated genes that control it. The aim of this study was to identify genomic regions associated with within-family residual variance of yearling weight (YW; N = 423) in Nellore bulls with high density SNP data, using different response variables. For this, solutions from double hierarchical generalized linear models (DHGLM) were used to provide the response variables, as follows: a DGHLM assuming non-null genetic correlation between mean and residual variance (rmv ≠ 0) to obtain deregressed EBV for mean (dEBVm) and residual variance (dEBVv); and a DHGLM assuming rmv = 0 to obtain two alternative response variables for residual variance, dEBVv_r0 and log-transformed variance of estimated residuals (ln_[Formula: see text]).

RESULTS

The dEBVm and dEBVv were highly correlated, resulting in common regions associated with mean and residual variance of YW. However, higher effects on variance than the mean showed that these regions had effects on the variance beyond scale effects. More independent association results between mean and residual variance were obtained when null rmv was assumed. While 13 and 4 single nucleotide polymorphisms (SNPs) showed a strong association (Bayes Factor > 20) with dEBVv and ln_[Formula: see text], respectively, only suggestive signals were found for dEBVv_r0. All overlapping 1-Mb windows among top 20 between dEBVm and dEBVv were previously associated with growth traits. The potential candidate genes for uniformity are involved in metabolism, stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation.

CONCLUSIONS

It is necessary to use a strategy like assuming null rmv to obtain genomic regions associated with uniformity that are not associated with the mean. Genes involved not only in metabolism, but also stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation were the most promising biological candidates for uniformity of YW. Although no clear evidence of using a specific response variable was found, we recommend consider different response variables to study uniformity to increase evidence on candidate regions and biological mechanisms behind it.

Breeding for robustness: investigating the genotype-by-environment interaction and micro-environmental sensitivity of Genetically Improved Farmed Tilapia (Oreochromis niloticus)

Robustness has become a highly desirable breeding goal in the globalized agricultural market. Both genotype‐by‐environment interaction (G × E) and micro‐environmental sensitivity are important robustness components of aquaculture production, in which breeding stock is often disseminated to different environments. The objectives of this study were (i) to quantify the degree of G × E by assessing the growth performance of Genetically Improved Farmed Tilapia (GIFT) across three countries (Malaysia, India and China) and (ii) to quantify the genetic heterogeneity of environmental variance for body weight at harvest (BW) in GIFT as a measure of micro‐environmental sensitivity. Selection for BW was carried out for 13 generations in Malaysia. Subsets of 60 full‐sib families from Malaysia were sent to China and India after five and nine generations respectively. First, a multi‐trait animal model was used to analyse the BW in different countries as different traits. The results indicate a strong G × E. Second, a genetically structured environmental variance model, implemented using Bayesian inference, was used to analyse micro‐environmental sensitivity of BW in each country. The analysis revealed the presence of genetic heterogeneity of both BW and its environmental variance in all environments. The presence of genetic variation in residual variance of BW implies that the residual variance can be modified by selection. Incorporating both G × E and micro‐environmental sensitivity information may help in selecting robust genotypes with high performance across environments and resilience to environmental fluctuations.

Directional selection reduces developmental canalization against genetic and environmental perturbations in Drosophila wings

Natural selection may enhance or weaken the robustness of phenotypes against genetic or environmental perturbations. However, important aspects of the relationship between adaptive evolution and canalization remain unclear. Recent work showed that the evolution of larger wing size in a high altitude natural population of Drosophila melanogaster was accompanied by decanalized wing development--specifically a loss of robustness to genetic perturbation. But this study did not address environmental robustness, and it compared populations that may have numerous biological differences. Here, we perform artificial selection on this same trait in D. melanogaster (larger wing length) and directly test whether this directional selection resulted in decanalization. We find that in general, size-selected replicates show greater frequencies of wing defects than control replicates both after mutagenesis (genetic perturbation) and when subjected to high temperature stress (environmental perturbation), although the increase in defect frequency varies importantly among replicates. These results support the hypothesis that directional selection may result in the loss of both genetic and environmental robustness-offering a rare window into the relationship between adaptation and canalization.

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.

Epistasis and destabilizing mutations shape gene expression variability in humans via distinct modes of action

Increasing evidence shows that phenotypic variance is genetically determined, but the underlying mechanisms of genetic control over the variance remain obscure. Here, we conducted variance-association mapping analyses to identify expression variability QTLs (evQTLs), i.e. genomic loci associated with gene expression variance, in humans. We discovered that common genetic variants may contribute to increasing gene expression variance via two distinct modes of action—epistasis and destabilization. Specifically, epistasis explains a quarter of the identified evQTLs, of which the formation is attributed to the presence of ‘third-party’ eQTLs that influence the gene expression mean in a fraction, rather than the entire set, of sampled individuals. On the other hand, the destabilization model explains the other three-quarters of evQTLs, caused by mutations that disrupt the stability of the transcription process of genes. To show the destabilizing effect, we measured discordant gene expression between monozygotic twins, and estimated the stability of gene expression in single samples using repetitive qRT-PCR assays. The mutations that cause destabilizing evQTLs were found to be associated with more pronounced expression discordance between twin pairs and less stable gene expression in single samples. Together, our results suggest that common genetic variants work either interactively or independently to shape the variability of gene expression in humans. Our findings contribute to the understanding of the mechanisms of genetic control over phenotypic variance and may have implications for the development of variance-centred analytic methods for quantitative trait mapping.

Genotype Influences Day-to-Day Variability in Sleep in Drosophila melanogaster

Patterns of sleep often vary among individuals. But sleep and activity may also vary within an individual, fluctuating in pattern across time. One possibility is that these daily fluctuations in sleep are caused by the underlying genotype of the individual. However, differences attributable to genetic causes are difficult to distinguish from environmental factors in outbred populations such as humans. We therefore employed Drosophila as a model of intra-individual variability in sleep using previously collected sleep and activity data from the Drosophila Genetic Reference Panel, a collection of wild-derived inbred lines. Individual flies had significant daily fluctuations in their sleep patterns, and these fluctuations were heritable. Using the standard deviation of sleep parameters as a metric, we conducted a genome-wide association study. We found 663 polymorphisms in 104 genes associated with daily fluctuations in sleep. We confirmed the effects of 12 candidate genes on the standard deviation of sleep parameters. Our results suggest that daily fluctuations in sleep patterns are due in part to gene activity.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Charting the genotype-phenotype map: lessons from the Drosophila melanogaster Genetic Reference Panel

Understanding the genetic architecture (causal molecular variants, their effects, and frequencies) of quantitative traits is important for precision agriculture and medicine and predicting adaptive evolution, but is challenging in most species. The Drosophila melanogaster Genetic Reference Panel (DGRP) is a collection of 205 inbred strains with whole genome sequences derived from a single wild population in Raleigh, NC, USA. The large amount of quantitative genetic variation, lack of population structure, and rapid local decay of linkage disequilibrium in the DGRP and outbred populations derived from DGRP lines present a favorable scenario for performing genome-wide association (GWA) mapping analyses to identify candidate causal genes, polymorphisms, and pathways affecting quantitative traits. The many GWA studies utilizing the DGRP have revealed substantial natural genetic variation for all reported traits, little evidence for variants with large effects but enrichment for variants with low P-values, and a tendency for lower frequency variants to have larger effects than more common variants. The variants detected in the GWA analyses rarely overlap those discovered using mutagenesis, and often are the first functional annotations of computationally predicted genes. Variants implicated in GWA analyses typically have sex-specific and genetic background-specific (epistatic) effects, as well as pleiotropic effects on other quantitative traits. Studies in the DGRP reveal substantial genetic control of environmental variation. Taking account of genetic architecture can greatly improve genomic prediction in the DGRP. These features of the genetic architecture of quantitative traits are likely to apply to other species, including humans. WIREs Dev Biol 2018, 7:e289. doi: 10.1002/wdev.289 This article is categorized under: Invertebrate Organogenesis > Flies.

Genetic variability of environmental sensitivity revealed by phenotypic variation in body weight and (its) correlations to physiological and behavioral traits

Adaptive phenotypic plasticity is a key component of the ability of organisms to cope with changing environmental conditions. Fish have been shown to exhibit a substantial level of phenotypic plasticity in response to abiotic and biotic factors. In the present study, we investigate the link between environmental sensitivity assessed globally (revealed by phenotypic variation in body weight) and more targeted physiological and behavioral indicators that are generally used to assess the sensitivity of a fish to environmental stressors. We took advantage of original biological material, the rainbow trout isogenic lines, which allowed the disentangling of the genetic and environmental parts of the phenotypic variance. Ten lines were characterized for the changes of body weight variability (weight measurements taken every month during 18 months), the plasma cortisol response to confinement stress (3 challenges) and a set of selected behavioral indicators. This study unambiguously demonstrated the existence of genetic determinism of environmental sensitivity, with some lines being particularly sensitive to environmental fluctuations and others rather insensitive. Correlations between coefficient of variation (CV) for body weight and behavioral and physiological traits were observed. This confirmed that CV for body weight could be used as an indicator of environmental sensitivity. As the relationship between indicators (CV weight, risk-taking, exploration and cortisol) was shown to be likely depending on the nature and intensity of the stressor, the joint use of several indicators should help to investigate the biological complexity of environmental sensitivity.

Selection for long and short sleep duration in Drosophila melanogaster reveals the complex genetic network underlying natural variation in sleep

Why do some individuals need more sleep than others? Forward mutagenesis screens in flies using engineered mutations have established a clear genetic component to sleep duration, revealing mutants that convey very long or short sleep. Whether such extreme long or short sleep could exist in natural populations was unknown. We applied artificial selection for high and low night sleep duration to an outbred population of Drosophila melanogaster for 13 generations. At the end of the selection procedure, night sleep duration diverged by 9.97 hours in the long and short sleeper populations, and 24-hour sleep was reduced to 3.3 hours in the short sleepers. Neither long nor short sleeper lifespan differed appreciably from controls, suggesting little physiological consequences to being an extreme long or short sleeper. Whole genome sequence data from seven generations of selection revealed several hundred thousand changes in allele frequencies at polymorphic loci across the genome. Combining the data from long and short sleeper populations across generations in a logistic regression implicated 126 polymorphisms in 80 candidate genes, and we confirmed three of these genes and a larger genomic region with mutant and chromosomal deficiency tests, respectively. Many of these genes could be connected in a single network based on previously known physical and genetic interactions. Candidate genes have known roles in several classic, highly conserved developmental and signaling pathways—EGFR, Wnt, Hippo, and MAPK. The involvement of highly pleiotropic pathway genes suggests that sleep duration in natural populations can be influenced by a wide variety of biological processes, which may be why the purpose of sleep has been so elusive.

One of the biggest mysteries in biology is the need to sleep. Sleep duration has an underlying genetic basis, suggesting that very long and short sleep times could be bred for experimentally. How far can sleep duration be driven up or down? Here we achieved extremely long and short night sleep duration by subjecting a wild-derived population of Drosophila melanogaster to an experimental breeding program. At the end of the breeding program, long sleepers averaged 9.97 hours more nightly sleep than short sleepers. We analyzed whole-genome sequences from seven generations of the experimental breeding to identify allele frequencies that diverged between long and short sleepers, and verified genes and genomic regions with mutation and deficiency testing. These alleles map to classic developmental and signaling pathways, implicating many diverse processes that potentially affect sleep duration.

Influence of physiological stress on the presence of hypoplasia and fluctuating asymmetry in a medieval population from the village of Sypniewo

This study aims to estimate the levels of physiological stress in the medieval rural population of Sypniewo by evaluating patterns of fluctuating asymmetry (FA) and enamel hypoplasia (EH), and provide information on the influence of physiological stress during the prenatal and perinatal period on early childhood development. Stress is defined as any external or internal condition that challenges homeostasis of an organism. FA is associated with physiological stress occurring mainly during prenatal development and early childhood. The level of FA is thought to reflect the intensity of the stressor(s). EH is caused by physiological stress such as nutritional instability during the first years of life. The studied material consisted of 126 skulls from the village of Sypniewo (Poland). Cranial radiographs were taken in postero-anterior (P-A) and basal views. The images were scanned and calibrated. Measurements of the cranium were used to estimate FA. The presence of EH was assessed using standard anthropological methods The highest levels of FA were observed in the region of the cranial base. EH was observed in 29% of individuals from the rural skeletal series. There was no statistically significant correlation between FA and EH occurrence or between sex and the studied stress indicators.

On the relationship between epistasis and genetic variance heterogeneity

Epistasis and genetic variance heterogeneity are two non-additive genetic inheritance patterns that are often, but not always, related. Here we use theoretical examples and empirical results from earlier analyses of experimental data to illustrate the connection between the two. This includes an introduction to the relationship between epistatic gene action, statistical epistasis, and genetic variance heterogeneity, and a brief discussion about how genetic processes other than epistasis can also give rise to genetic variance heterogeneity.

Artificial Selection to Increase the Phenotypic Variance in gmax Fails

Stabilizing selection is important in evolutionary theories of the maintenance of genetic variance and has been invoked as the key process determining macroevolutionary patterns of trait evolution. However, manipulative evidence for the extent of stabilizing selection, particularly on multivariate traits, is lacking. We used artificial disruptive selection in Drosophila serrata as a tool to determine the relative strength of stabilizing selection experienced by multivariate trait combinations with contrasting levels of genetic and mutational variance. Contrary to expectation, when disruptive selection was applied to the major axis of standing genetic variance, g_max, we observed a significant and repeatable decrease in its phenotypic variance. In contrast, the multivariate trait combination predicted to be under strong stabilizing selection showed a significant and repeatable increase in its phenotypic variance. Correlated responses were observed in all selection treatments, and viability selection operating on extreme phenotypes of traits genetically correlated with those directly selected on limited our ability to increase their phenotypic range. Our manipulation revealed that multivariate trait combinations were subject to stabilizing selection; however, we did not observe a direct relationship between the strength of stabilizing selection and the levels of standing genetic variance in multivariate trait combinations. Contrasting patterns of allele frequencies underlying traits with high versus low levels of standing genetic variance may be implicated in determining the response to artificial selection in multivariate trait combinations.

Fluctuating selection across years and phenotypic variation in food-deceptive orchids

Nectarless flowers that deceive pollinators offer an opportunity to study asymmetric plant-insect interactions. Orchids are a widely used model for studying these interactions because they encompass several thousand species adopting deceptive pollination systems. High levels of intra-specific phenotypic variation have been reported in deceptive orchids, suggesting a reduced consistency of pollinator-mediated selection on their floral traits. Nevertheless, several studies report on widespread directional selection mediated by pollinators even in these deceptive orchids. In this study we test the hypothesis that the observed selection can fluctuate across years in strength and direction thus likely contributing to the phenotypic variability of this orchid group. We performed a three-year study estimating selection differentials and selection gradients for nine phenotypic traits involved in insect attraction in two Mediterranean orchid species, namely Orchis mascula and O. pauciflora, both relying on a well-described food-deceptive pollination strategy. We found weak directional selection and marginally significant selection gradients in the two investigated species with significant intra-specific differences in selection differentials across years. Our data do not link this variation with a specific environmental cause, but our results suggest that pollinator-mediated selection in food-deceptive orchids can change in strength and in direction over time. In perennial plants, such as orchids, different selection differentials in the same populations in different flowering seasons can contribute to the maintenance of phenotypic variation often reported in deceptive orchids.

The Genetic Architecture of Ovariole Number in Drosophila melanogaster: Genes with Major, Quantitative, and Pleiotropic Effects

Ovariole number has a direct role in the number of eggs produced by an insect, suggesting that it is a key morphological fitness trait. Many studies have documented the variability of ovariole number and its relationship to other fitness and life-history traits in natural populations of Drosophila However, the genes contributing to this variability are largely unknown. Here, we conducted a genome-wide association study of ovariole number in a natural population of flies. Using mutations and RNAi-mediated knockdown, we confirmed the effects of 24 candidate genes on ovariole number, including a novel gene, anneboleyn (formerly CG32000), that impacts both ovariole morphology and numbers of offspring produced. We also identified pleiotropic genes between ovariole number traits and sleep and activity behavior. While few polymorphisms overlapped between sleep parameters and ovariole number, 39 candidate genes were nevertheless in common. We verified the effects of seven genes on both ovariole number and sleep: bin3, blot, CG42389, kirre, slim, VAChT, and zfh1 Linkage disequilibrium among the polymorphisms in these common genes was low, suggesting that these polymorphisms may evolve independently.

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.

Selection for environmental variance of litter size in rabbits

Background

In recent years, there has been an increasing interest in the genetic determination of environmental variance. In the case of litter size, environmental variance can be related to the capacity of animals to adapt to new environmental conditions, which can improve animal welfare.

Results

We developed a ten-generation divergent selection experiment on environmental variance. We selected one line of rabbits for litter size homogeneity and one line for litter size heterogeneity by measuring intra-doe phenotypic variance. We proved that environmental variance of litter size is genetically determined and can be modified by selection. Response to selection was 4.5% of the original environmental variance per generation. Litter size was consistently higher in the Low line than in the High line during the entire experiment.

Conclusions

We conclude that environmental variance of litter size is genetically determined based on the results of our divergent selection experiment. This has implications for animal welfare, since animals that cope better with their environment have better welfare than more sensitive animals. We also conclude that selection for reduced environmental variance of litter size does not depress litter size.

Genetic control of residual variance of yearling weight in Nellore beef cattle

There is evidence for genetic variability in residual variance of livestock traits, which offers the potential for selection for increased uniformity of production. Different statistical approaches have been employed to study this topic; however, little is known about the concordance between them. The aim of our study was to investigate the genetic heterogeneity of residual variance on yearling weight (YW; 291.15 ± 46.67) in a Nellore beef cattle population; to compare the results of the statistical approaches, the two-step approach and the double hierarchical generalized linear model (DHGLM); and to evaluate the effectiveness of power transformation to accommodate scale differences. The comparison was based on genetic parameters, accuracy of EBV for residual variance, and cross-validation to assess predictive performance of both approaches. A total of 194,628 yearling weight records from 625 sires were used in the analysis. The results supported the hypothesis of genetic heterogeneity of residual variance on YW in Nellore beef cattle and the opportunity of selection, measured through the genetic coefficient of variation of residual variance (0.10 to 0.12 for the two-step approach and 0.17 for DHGLM, using an untransformed data set). However, low estimates of genetic variance associated with positive genetic correlations between mean and residual variance (about 0.20 for two-step and 0.76 for DHGLM for an untransformed data set) limit the genetic response to selection for uniformity of production while simultaneously increasing YW itself. Moreover, large sire families are needed to obtain accurate estimates of genetic merit for residual variance, as indicated by the low heritability estimates (<0.007). Box-Cox transformation was able to decrease the dependence of the variance on the mean and decreased the estimates of genetic parameters for residual variance. The transformation reduced but did not eliminate all the genetic heterogeneity of residual variance, highlighting its presence beyond the scale effect. The DHGLM showed higher predictive ability of EBV for residual variance and therefore should be preferred over the two-step approach.

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.

Estimation of breeding values for uniformity of growth in Atlantic salmon (Salmo salar) using pedigree relationships or single-step genomic evaluation

Background

In farmed Atlantic salmon, heritability for uniformity of body weight is low, indicating that the accuracy of estimated breeding values (EBV) may be low. The use of genomic information could be one way to increase accuracy and, hence, obtain greater response to selection. Genomic information can be merged with pedigree information to construct a combined relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H matrix) for a single-step genomic evaluation (ssGBLUP), allowing realized relationships of the genotyped animals to be exploited, in addition to numerator pedigree relationships (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A matrix). We compared the predictive ability of EBV for uniformity of body weight in Atlantic salmon, when implementing either the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H matrix in the genetic evaluation. We used double hierarchical generalized linear models (DHGLM) based either on a sire-dam (sire-dam DHGLM) or an animal model (animal DHGLM) for both body weight and its uniformity.

Results

With the animal DHGLM, the use of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H instead of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A significantly increased the correlation between the predicted EBV and adjusted phenotypes, which is a measure of predictive ability, for both body weight and its uniformity (41.1 to 78.1%). When log-transformed body weights were used to account for a scale effect, the use of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H instead of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A produced a small and non-significant increase (1.3 to 13.9%) in predictive ability. The sire-dam DHGLM had lower predictive ability for uniformity compared to the animal DHGLM.

Conclusions

Use of the combined numerator and genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H) significantly increased the predictive ability of EBV for uniformity when using the animal DHGLM for untransformed body weight. The increase was only minor when using log-transformed body weights, which may be due to the lower heritability of scaled uniformity, the lower genetic correlation of transformed body weight with its uniformity compared to the untransformed traits, and the small number of genotyped animals in the reference population. This study shows that ssGBLUP increases the accuracy of EBV for uniformity of body weight and is expected to increase response to selection in uniformity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0308-3) contains supplementary material, which is available to authorized users.

Leger RJ. The genetic basis for variation in resistance to infection in the Drosophila melanogaster genetic reference panel

Individuals vary extensively in the way they respond to disease but the genetic basis of this variation is not fully understood. We found substantial individual variation in resistance and tolerance to the fungal pathogen Metarhizium anisopliae Ma549 using the Drosophila melanogaster Genetic Reference Panel (DGRP). In addition, we found that host defense to Ma549 was correlated with defense to the bacterium Pseudomonas aeruginosa Pa14, and several previously published DGRP phenotypes including oxidative stress sensitivity, starvation stress resistance, hemolymph glucose levels, and sleep indices. We identified polymorphisms associated with differences between lines in both their mean survival times and microenvironmental plasticity, suggesting that lines differ in their ability to adapt to variable pathogen exposures. The majority of polymorphisms increasing resistance to Ma549 were sex biased, located in non-coding regions, had moderately large effect and were rare, suggesting that there is a general cost to defense. Nevertheless, host defense was not negatively correlated with overall longevity and fecundity. In contrast to Ma549, minor alleles were concentrated in the most Pa14-susceptible as well as the most Pa14-resistant lines. A pathway based analysis revealed a network of Pa14 and Ma549-resistance genes that are functionally connected through processes that encompass phagocytosis and engulfment, cell mobility, intermediary metabolism, protein phosphorylation, axon guidance, response to DNA damage, and drug metabolism. Functional testing with insertional mutagenesis lines indicates that 12/13 candidate genes tested influence susceptibility to Ma549. Many candidate genes have homologs identified in studies of human disease, suggesting that genes affecting variation in susceptibility are conserved across species.

Derivation of economic values for production traits in aquaculture species

BACKGROUND

In breeding programs for aquaculture species, breeding goal traits are often weighted based on the desired gains but economic gain would be higher if economic values were used instead. The objectives of this study were: (1) to develop a bio-economic model to derive economic values for aquaculture species, (2) to apply the model to determine the economic importance and economic values of traits in a case-study on gilthead seabream, and (3) to validate the model by comparison with a profit equation for a simplified production system.

METHODS

A bio-economic model was developed to simulate a grow-out farm for gilthead seabream, and then used to simulate gross margin at the current levels of the traits and after one genetic standard deviation change in each trait with the other traits remaining unchanged. Economic values were derived for the traits included in the breeding goal: thermal growth coefficient (TGC), thermal feed intake coefficient (TFC), mortality rate (M), and standard deviation of harvest weight ([Formula: see text]). For a simplified production system, improvement in TGC was assumed to affect harvest weight instead of growing period. Using the bio-economic model and a profit equation, economic values were derived for harvest weight, cumulative feed intake at harvest, and overall survival.

RESULTS

Changes in gross margin showed that the order of economic importance of the traits was: TGC, TFC, M, and [Formula: see text]. Economic values in € (kg production)^-1 (trait unit)^-1 were: 0.40 for TGC, -0.45 for TFC, -7.7 for M, and -0.0011 to -0.0010 for [Formula: see text]. For the simplified production system, similar economic values were obtained with the bio-economic model and the profit equation. The advantage of the profit equation is its simplicity, while that of the bio-economic model is that it can be applied to any aquaculture species, because it can include any limiting factor and/or environmental condition that affects production.

CONCLUSIONS

We confirmed the validity of the bio-economic model. TGC is the most important trait to improve, followed by TFC and M, and the effect of [Formula: see text] on gross margin is small.

Joint estimation of (co) variance components and breeding values for mean and dispersion of days from calving to first service in Holstein cow

The present study explored the possibility of selection for uniformity of days from calving to first service (DFS) in dairy cattle. A double hierarchical generalised linear model with an iterative reweighted least-squares algorithm was used to estimate covariance components for the mean and dispersion of DFS. Data included the records of 27 113 Iranian Holstein cows (parity, 1–6) in 15 herds from 1981 to 2007. The estimated additive genetic variance for the mean and dispersion were 32.25 and 0.0139; both of these values had low standard errors. The genetic standard deviation for dispersion of DFS was 0.117, indicating that decreasing the estimated breeding value of dispersion by one genetic standard deviation can increase the uniformity by 12%. A strong positive genetic correlation (0.689) was obtained between the mean and dispersion of DFS. This genetic correlation is favourable since one of the aims of breeding is to simultaneously decrease the mean and increase the uniformity of DFS. The Spearman rank correlations between estimated breeding values in the mean and dispersion for sires with a different number of daughter observations were 0.907. In the studied population, the genetic trend in the mean of DFS was significant and favourable (–0.063 days/year), but the genetic trend in the dispersion of DFS was not significantly different from zero. The results obtained in the present study indicated that the mean and uniformity of DFS can simultaneously be improved in dairy cows.

What's wrong with evolutionary biology?

There have been periodic claims that evolutionary biology needs urgent reform, and this article tries to account for the volume and persistence of this discontent. It is argued that a few inescapable properties of the field make it prone to criticisms of predictable kinds, whether or not the criticisms have any merit. For example, the variety of living things and the complexity of evolution make it easy to generate data that seem revolutionary (e.g. exceptions to well-established generalizations, or neglected factors in evolution), and lead to disappointment with existing explanatory frameworks (with their high levels of abstraction, and limited predictive power). It is then argued that special discontent stems from misunderstandings and dislike of one well-known but atypical research programme: the study of adaptive function, in the tradition of behavioural ecology. To achieve its goals, this research needs distinct tools, often including imaginary agency, and a partial description of the evolutionary process. This invites mistaken charges of narrowness and oversimplification (which come, not least, from researchers in other subfields), and these chime with anxieties about human agency and overall purpose. The article ends by discussing several ways in which calls to reform evolutionary biology actively hinder progress in the field.

Microenvironmental Gene Expression Plasticity Among Individual Drosophila melanogaster

Differences in phenotype among genetically identical individuals exposed to the same environmental condition are often noted in genetic studies. Despite this commonplace observation, little is known about the causes of this variability, which has been termed microenvironmental plasticity. One possibility is that stochastic or technical sources of variance produce these differences. A second possibility is that this variation has a genetic component. We have explored gene expression robustness in the transcriptomes of 730 individual Drosophila melanogaster of 16 fixed genotypes, nine of which are infected with Wolbachia. Three replicates of flies were grown, controlling for food, day/night cycles, humidity, temperature, sex, mating status, social exposure, and circadian timing of RNA extraction. Despite the use of inbred genotypes, and carefully controlled experimental conditions, thousands of genes were differentially expressed, revealing a unique and dynamic transcriptional signature for each individual fly. We found that 23% of the transcriptome was differentially expressed among individuals, and that the variability in gene expression among individuals is influenced by genotype. This transcriptional variation originated from specific gene pathways, suggesting a plastic response to the microenvironment; but there was also evidence of gene expression differences due to stochastic fluctuations. These observations reveal previously unappreciated genetic sources of variability in gene expression among individuals, which has implications for complex trait genetics and precision medicine.

Disintegrating the fly: A mutational perspective on phenotypic integration and covariation

The structure of environmentally induced phenotypic covariation can influence the effective strength and magnitude of natural selection. Yet our understanding of the factors that contribute to and influence the evolutionary lability of such covariation is poor. Most studies have either examined environmental variation without accounting for covariation, or examined phenotypic and genetic covariation without distinguishing the environmental component. In this study, we examined the effect of mutational perturbations on different properties of environmental covariation, as well as mean shape. We use strains of Drosophila melanogaster bearing well-characterized mutations known to influence wing shape, as well as naturally derived strains, all reared under carefully controlled conditions and with the same genetic background. We find that mean shape changes more freely than the covariance structure, and that different properties of the covariance matrix change independently from each other. The perturbations affect matrix orientation more than they affect matrix eccentricity or total variance. Yet, mutational effects on matrix orientation do not cluster according to the developmental pathway that they target. These results suggest that it might be useful to consider a more general concept of "decanalization," involving all aspects of variation and covariation.

Genomic Selection Improves Response to Selection in Resilience by Exploiting Genotype by Environment Interactions

Genotype by environment interactions (GxE) are very common in livestock and hamper genetic improvement. On the other hand, GxE is a source of genetic variation: genetic variation in response to environment, e.g., environmental perturbations such as heat stress or disease. In livestock breeding, there is tendency to ignore GxE because of increased complexity of models for genetic evaluations and lack of accuracy in extreme environments. GxE, however, creates opportunities to increase resilience of animals toward environmental perturbations. The main aim of the paper is to investigate to which extent GxE can be exploited with traditional and genomic selection methods. Furthermore, we investigated the benefit of reaction norm (RN) models compared to conventional methods ignoring GxE. The questions were addressed with selection index theory. GxE was modeled according to a linear RN model in which the environmental gradient is the contemporary group mean. Economic values were based on linear and non-linear profit equations. Accuracies of environment-specific (G)EBV were highest in intermediate environments and lowest in extreme environments. RN models had higher accuracies of (G)EBV in extreme environments than conventional models ignoring GxE. Genomic selection always resulted in higher response to selection in all environments than sib or progeny testing schemes. The increase in response was with genomic selection between 9 and 140% compared to sib testing and between 11 and 114% compared to progeny testing when the reference population consisted of 1 million animals across all environments. When the aim was to decrease environmental sensitivity, the response in slope of the RN model with genomic selection was between 1.09 and 319 times larger than with sib or progeny testing and in the right direction in contrast to sib and progeny testing that still increased environmental sensitivity. This shows that genomic selection with large reference populations offers great opportunities to exploit GxE to increase resilience of animals.

Genetic parameters for uniformity of harvest weight and body size traits in the GIFT strain of Nile tilapia

Background

Animal breeding programs have been very successful in improving the mean levels of traits through selection. However, in recent decades, reducing the variability of trait levels between individuals has become a highly desirable objective. Reaching this objective through genetic selection requires that there is genetic variation in the variability of trait levels, a phenomenon known as genetic heterogeneity of environmental (residual) variance. The aim of our study was to investigate the potential for genetic improvement of uniformity of harvest weight and body size traits (length, depth, and width) in the genetically improved farmed tilapia (GIFT) strain. In order to quantify the genetic variation in uniformity of traits and estimate the genetic correlations between level and variance of the traits, double hierarchical generalized linear models were applied to individual trait values.

Results

Our results showed substantial genetic variation in uniformity of all analyzed traits, with genetic coefficients of variation for residual variance ranging from 39 to 58 %. Genetic correlation between trait level and variance was strongly positive for harvest weight (0.60 ± 0.09), moderate and positive for body depth (0.37 ± 0.13), but not significantly different from 0 for body length and width.

Conclusions

Our results on the genetic variation in uniformity of harvest weight and body size traits show good prospects for the genetic improvement of uniformity in the GIFT strain. A high and positive genetic correlation was estimated between level and variance of harvest weight, which suggests that selection for heavier fish will also result in more variation in harvest weight. Simultaneous improvement of harvest weight and its uniformity will thus require index selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0218-9) contains supplementary material, which is available to authorized users.

Constraints on geographic variation in fiddler crabs (Ocypodidae: Uca) from the western Atlantic

A key question in evolutionary biology is how intraspecific variation biases the evolution of a population and its divergence from other populations. Such constraints potentially limit the extent to which populations respond to selection, but may endure long enough to have macroevolutionary consequences. Previous studies have focused on the association between covariation patterns and divergence among isolated populations. Few have focused on geographic variation among semi-connected populations, however, even though this may be indicative of early selective pressures that could lead to long-term divergence and speciation. Here, we test whether covariation in the shape of the carapace of fiddler crabs (genus Uca Leach, 1814) is important for structuring geographic variation. We find that morphological divergence among populations is associated with evolvability in the direction of divergence in only a few species. The shape of the ancestral covariation matrix in these species differs from other species in having notably more variation concentrated along fewer directions (i.e. higher eccentricity). For most species, there is some evidence that covariation has constrained the range of directions into which populations have diverged but not the degree of divergence. These results suggest that even though fiddler crab populations have diverged morphologically in directions predicted by covariation, constraints on the extent to which divergence has occurred may only be manifested in species where variation patterns are eccentric enough to limit populations' ability to respond effectively in many directions.

Estimating the purebred-crossbred genetic correlation for uniformity of eggshell color in laying hens

Background

Uniformity of eggs is an important aspect for retailers because consumers prefer homogeneous products. One of these characteristics is the color of the eggshell, especially for brown eggs. Existence of a genetic component in environmental variance would enable selection for uniformity of eggshell color. Therefore, the objective of this study was to quantify the genetic variance in environmental variance of eggshell color in purebred and crossbred laying hens, to estimate the genetic correlation between environmental variance of eggshell color in purebred and crossbred laying hens and to estimate genetic correlations between environmental variance at different times of the laying period.

Methods

We analyzed 167,651 and 79,345 eggshell color records of purebred and crossbred laying hens, respectively. The purebred and crossbred laying hens originated mostly from the same sires. Since eggshell color records of crossbred laying hens were collected per cage, these records could be related only to cage and sire family. A double hierarchical generalized linear sire model was used to estimate the genetic variance of the mean of eggshell color and its environmental variance. Approximate standard errors for heritability and the genetic coefficient of variation for environmental variance were derived.

Results

The genetic variance in environmental variance at the log scale was equal to 0.077 and 0.067, for purebred and crossbred laying hens, respectively. The genetic coefficient of variation for environmental variance was equal to 0.28 and 0.26, for purebred and crossbred laying hens, respectively. A genetic correlation of 0.70 was found between purebred and crossbred environmental variance of eggshell color, which indicates that there is some reranking of sires for environmental variance of eggshell color in purebred and crossbred laying hens. Genetic correlations between environmental variance of eggshell color in different laying periods were generally higher than 0.85, except between early laying and mid or late laying periods.

Conclusions

Our results indicate that genetic selection can be efficient to improve uniformity of eggshell color in purebreds and crossbreds, ideally by applying combined crossbred and purebred selection. This methodology can be used to estimate genetic correlations between purebred and crossbred lines for uniformity of other traits and species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0212-2) contains supplementary material, which is available to authorized users.

Genetic analysis of within-litter variation in piglets' birth weight using genomic or pedigree relationship matrices

The objective of this study was to estimate the genetic variance for within-litter variation of birth weight (BW0) using genomic (GRM) or pedigree relationship matrices (PRM) and to compare the accuracy of estimated breeding values (EBV) for within-litter variation of BW0 using GRM and PRM. The BW0 and residual variance of BW0 were modeled by the double hierarchical generalized linear model using GRM or PRM. Data came from 2 dam lines: Landrace and Large White. After editing, the data set in Landrace consisted of 748 sows with 1,938 litters and 29,430 piglets and in Large White of 989 sows with 3,320 litters and 51,818 piglets. To construct GRM, 46,466 (Landrace) and 44,826 (Large White) single nucleotide polymorphisms were used, whereas to construct PRM, 5 generations of pedigree were used. The accuracy of EBV with GRM was estimated with 8-fold cross-validation and compared to PRM. Estimated variance components were highly similar for GRM and PRM. The maternal genetic variance in residual variance of BW0 in Landrace was 0.05 with GRM and 0.06 with PRM. In Large White these were 0.04 with GRM and 0.05 with PRM. The genetic coefficient of variation (GCV SDe) was about 0.10 in both dam lines. This indicates a change of 10% in residual SD of BW0 when achieving a genetic response of 1 genetic standard deviation. The genetic correlation between birth weight and its residual variance was about 0.6 in both dam lines. The accuracies of selection for within-litter variation of birth weight were 0.35 with GRM and 0.23 with PRM in Landrace and 0.29 with GRM and 0.34 with PRM in Large White. In this case, using GRM did not significantly increase accuracies of selection. Results, however, show good opportunities to select for reduced within-litter variation of BW0. Genomic selection can increase accuracy of selection when reference populations contain at least 2,000 sows.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population

Background

In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.

Results

In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.

Conclusions

To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

The Multi-allelic Genetic Architecture of a Variance-Heterogeneity Locus for Molybdenum Concentration in Leaves Acts as a Source of Unexplained Additive Genetic Variance

Genome-wide association (GWA) analyses have generally been used to detect individual loci contributing to the phenotypic diversity in a population by the effects of these loci on the trait mean. More rarely, loci have also been detected based on variance differences between genotypes. Several hypotheses have been proposed to explain the possible genetic mechanisms leading to such variance signals. However, little is known about what causes these signals, or whether this genetic variance-heterogeneity reflects mechanisms of importance in natural populations. Previously, we identified a variance-heterogeneity GWA (vGWA) signal for leaf molybdenum concentrations in Arabidopsis thaliana. Here, fine-mapping of this association reveals that the vGWA emerges from the effects of three independent genetic polymorphisms that all are in strong LD with the markers displaying the genetic variance-heterogeneity. By revealing the genetic architecture underlying this vGWA signal, we uncovered the molecular source of a significant amount of hidden additive genetic variation or “missing heritability”. Two of the three polymorphisms underlying the genetic variance-heterogeneity are promoter variants for Molybdate transporter 1 (MOT1), and the third a variant located ~25 kb downstream of this gene. A fourth independent association was also detected ~600 kb upstream of MOT1. Use of a T-DNA knockout allele highlights Copper Transporter 6; COPT6 (AT2G26975) as a strong candidate gene for this association. Our results show that an extended LD across a complex locus including multiple functional alleles can lead to a variance-heterogeneity between genotypes in natural populations. Further, they provide novel insights into the genetic regulation of ion homeostasis in A. thaliana, and empirically confirm that variance-heterogeneity based GWA methods are a valuable tool to detect novel associations of biological importance in natural populations.

Most biological traits vary in natural populations, and understanding the genetic basis of this variation remains an important challenge. Genome-wide association (GWA) studies have emerged as a powerful tool to address this challenge by dissecting the genetic architecture of trait variation into the contribution of individual genes. This contribution has traditionally been measured as the difference in the phenotypic means between groups of individuals with alternative genotypes at one, or multiple loci. However, instead of altering the trait mean, certain loci alter the variability of the trait. Here, we describe the genetic dissection of one such variance-controlling locus that drives variation in leaf molybdenum concentrations amongst natural accessions of Arabidopsis thaliana. The variance-controlling locus was found to result from the contributions of multiple alleles at multiple loci that are closely linked on the chromosome and is a major contributor to the “missing heritability” for this trait identified in previous studies. This illustrates that multi-allelic genetic architectures can hide large amounts of additive genetic variation, and that it is possible to uncover this hidden variation using the appropriate experimental designs and statistical methods described here.