Something Old and Something New: Wedding Recombinant Inbred Lines with Traditional Line Cross Analysis Increases Power to Describe Gene Interactions

In this paper we present a novel approach to quantifying genetic architecture that combines recombinant inbred lines (RIL) with line cross analysis (LCA). LCA is a method of quantifying directional genetic effects (i.e. summed effects of all loci) that differentiate two parental lines. Directional genetic effects are thought to be critical components of genetic architecture for the long term response to selection and as a cause of inbreeding depression. LCA typically begins with two inbred parental lines that are crossed to produce several generations such as F1, F2, and backcrosses to each parent. When a RIL population (founded from the same P1 and P2 as was used to found the line cross population) is added to the LCA, the sampling variance of several nonadditive genetic effect estimates is greatly reduced. Specifically, estimates of directional dominance, additive x additive, and dominance x dominance epistatic effects are reduced by 92%, 94%, and 56% respectively. The RIL population can be simultaneously used for QTL identification, thus uncovering the effects of specific loci or genomic regions as elements of genetic architecture. LCA and QTL mapping with RIL provide two qualitatively different measures of genetic architecture with the potential to overcome weaknesses of each approach alone. This approach provides cross-validation of the estimates of additive and additive x additive effects, much smaller confidence intervals on dominance, additive x additive and dominance x dominance estimates, qualitatively different measures of genetic architecture, and the potential when used together to balance the weaknesses of LCA or RIL QTL analyses when used alone.     

Genetic Architecture of Two Fitness-related Traits in Drosophila melanogaster: Ovariole Number and Thorax Length

In Drosophila melanogaster, ovariole number and thorax length are morphological characters thought to be associated with fitness. Maximum daily egg production in females is positively correlated with ovariole number, while thorax length is correlated with male reproductive success and female fecundity. Though both traits are related to fitness, ovariole number is likely to be under stabilizing selection, while thorax length appears to be under directional selection. Current research has focused on examining the sources of variation for ovariole number in relation to fitness, with a view towards elucidating how segregating variation is maintained in natural populations. Here, we utilize a diallel design to explore the genetic architecture of ovariole number and thorax length in nine isogenic lines derived from a natural population. The full diallel design allows the estimation of general combining ability (GCA), specific combining ability (SCA), and also describes variation due to reciprocal effects (RGCA and RSCA). Ovariole number and thorax length differed with respect to their genetic architecture, reflective of the independent selective forces acting on the traits. For ovariole number, GCA accounted for the majority (67.3%) of variation segregating between the lines, with no evidence of reciprocal effects or inbreeding depression; SCA accounted for a small percentage (3.9%) of the variance, suggesting dominance variation; no reciprocal effects were observed. In contrast, for thorax length, the majority of the non-error variance was accounted for by SCA (17.9%), with only one third as much variance (6.2%) due to GCA. Interestingly, RSCA (nuclear–extranuclear interactions) accounted for slightly more variation (7.5%) than GCA in these data. Thus, genetic variation for thorax length is largely in accord with predictions for a fitness trait under directional selection: little additive genetic variation and substantial dominance variation (including a suggestion of inbreeding depression); while the mechanisms underlying the maintenance of variation for ovariole number are more complex.     

Phenotypic plasticity for life-history traits in Drosophila melanogaster. III. Effect of the environment on genetic parameters.

We estimated genetic and environmental variance components for developmental time and dry weight at eclosion in Drosophila melanogaster raised in ten different environments (all combinations of 22, 25 and 28 degrees C and 0.5, 1 and 4% yeast concentration, and 0.25% yeast at 25 degrees C). We used six homozygous lines derived from a natural population for complete diallel crosses in each environment. Additive genetic variances were consistently low for both traits (h2 around 10%). The additive genetic variance of developmental time was larger at lower yeast concentrations, but the heritability did not increase because other components were also larger. The additive genetic effects of the six parental lines changed ranks across environments, suggesting a mechanism for the maintenance of genetic variation in heterogenous environments. The variance due to non-directional dominance was small in most environments. However, there was directional dominance in the form of inbreeding depression for both traits. It was pronounced at high yeast levels and temperatures but disappeared when yeast or temperature were decreased. This meant that the heterozygous flies were more sensitive to environmental differences than homozygous flies. Because dominance effects are not heritable, this suggests that the evolution of plasticity can be constrained when dominance effects are important as a mechanism for plasticity.     

Protected polymorphism and evolutionary stability in pleiotropic models with trait-specific dominance

When alleles have pleiotropic effects on a number of quantitative traits, the degree of dominance between a pair of alleles can be different for each trait. Such trait-specific dominance has been studied previously in models for the maintenance of genetic variation by antagonistic effects of an allele on two fitness components. By generalizing these models to an arbitrary number of fitness components or other phenotypic traits with different degrees of dominance, I show that genetic polymorphism is generally impossible without antagonistic fitness effects of different traits and without trait-specific dominance. I also investigate dominance and pleiotropy from a more long-term evolutionary perspective, allowing for the study of general ecological scenarios, and I discuss the effects of trait-specific dominance on evolutionary stability criteria. When selection is mainly directional and only trait-specific dominance and antagonism cause the emergence of polymorphism, then these polymorphisms can be overtaken by single mutants again, such that they are probably short-lived on an evolutionary time scale. Near evolutionarily singular points where directional selection is absent, trait-specific dominance and overdominance facilitate the emergence of polymorphism and cause evolutionary divergence in some cases. An important outcome of these models is that trait-specific dominance allows for the emergence of genetic polymorphisms without a selective disadvantage for heterozygotes. This removes the scope for the evolution of assortative mate choice and affects dominance modification. Sympatric speciation by disruptive ecological selection requires this heterozygote disadvantage in order to evolve, and therefore it becomes less plausible if the emergence of genetic polymorphism usually occurs via trait-specific dominance and antagonistic effects.     

Linking genetic mechanisms of heterozygosity-fitness correlations to footprints of selection at single loci

Investigations of heterozygosity-fitness correlations (HFCs) are central to the understanding how genetic diversity is maintained in natural populations. Advanced genome-wide approaches will enrich the number of functional loci to be tested. We argue that a combined analysis of the genetic mechanisms of HFCs and selection signals at single loci will allow researchers to better understand the micro-evolutionary basis of HFCs. Different dominance relationships among the alleles at the locus can lead to positive, negative or null HFCs depending on the allele frequency distribution. These scenarios differ in the temporal stability of the HFCs and in the patterns of allele frequency changes over time. Here, we describe a simple theoretical framework that links the analyses of heterozygosity-fitness associations (ecological timescale) with tests for selection signals (evolutionary timescale). Different genomic footprints of selection can be expected for the different underlying genetic mechanisms of HFCs, and this information can be independently used for the classification of HFCs. We suggest that in addition to inbreeding and single-locus overdominant effects also loci under directional selection could play a significant role in the development of heterozygosity-fitness effects in large natural populations under recent or fluctuating ecological changes.     

A quantitative model of the relationship between phenotypic variance and heterozygosity at marker loci under partial selfing.

Negative relationships between allozyme heterozygosity and morphological variance have often been observed and interpreted as evidence for increased developmental stability in heterozygotes. However, inbreeding can also generate such relationships by decreasing heterozygosity at neutral loci and redistributing genetic variance at the same time. I here provide a quantitative genetic model of this process by analogy with heterozygosity-fitness relationships. Inbreeding generates negative heterozygosity-variance relationships irrespective of the genetic architecture of the trait. This holds for fitness traits as well as neutral traits, the effect being stronger for fitness traits under directional dominance or overdominance. The order of magnitude of heterozygosity-variance regressions is compatible with empirical data even with very low inbreeding. Although developmental stability effects cannot be excluded, inbreeding is a parsimonious explanation that should be seriously considered to explain correlations between heterozygosity and both mean and variance of phenotypes in natural populations.     

The quantitative genetics of fluctuating asymmetry

Fluctuating asymmetry (subtle departures from identical expression of a trait across an axis of symmetry) in many taxa is under stabilizing selection for reduced asymmetry. However, lack of reliable estimates of genetic parameters for asymmetry variation hampers our ability to predict the evolutionary outcome of this selection. Here we report on a study, based on analysis of variation within and between isofemale lines and of generation means (line-cross analysis), designed to dissect in detail the quantitative genetics of positional fluctuating asymmetry (PFA) in bristle number in natural populations of Drosophila falleni. PFA is defined as the difference between the two sides of the body in the placement or position of components of a meristic trait. Heritability (measured at 25 degrees C) of two related measures of PFA were 13% and 21%, both of which differed significantly from zero. In contrast, heritability estimates for fluctuating asymmetry in the total number of anterior (0.7%) and transverse (2.4%) sternopleural bristles were smaller, not significant, and in quantitative agreement with previously published estimates. Heritabilities for bristle number (trait size) were considerably greater than that for any asymmetry measure. The experimental design controlled for the potentially confounding effects of common familial environment, and repeated testing revealed that PFA differences between lines were genetically stable for up to 16 generations in the laboratory at 25 degrees C. We performed line cross analysis between strains at the extremes of the PFA distribution (highest and lowest values); parental strains, F1, F1r (reciprocal), F2, backcross, and backcross reciprocal generations were represented. The inheritance of PFA was described best by additive and dominance effects localized to the X-chromosomes, whereas autosomal dominance effects were also detected. Epistatic, maternal, and cytoplasmic effects were not detected. The inheritance of trait size was notably more complex and involved significant autosomal additive, dominance, and epistatic effects; maternal dominance effects; and additive and dominance effects localized to the X-chromosomes. The additive genetic correlation between PFA and its associated measure of trait size was negative (-0.049), but not statistically significant, indicating that the loci contributing additive genetic effects to these traits are probably different. It is suggested that PFA may be a sensitive measure of developmental instability because PFA taps the ability of an organism to integrate interconnected developmental pathways.     

Effects of bottlenecks on quantitative genetic variation in the butterfly Bicyclus anynana

The effects of a single population bottleneck of differing severity on heritability and additive genetic variance was investigated experimentally using a butterfly. An outbred laboratory stock was used to found replicate lines with one pair, three pairs and 10 pairs of adults, as well as control lines with approximately 75 effective pairs. Heritability and additive genetic variance of eight wing pattern characters and wing size were estimated using parent-offspring covariances in the base population and in all daughter lines. Individual morphological characters and principal components of the nine characters showed a consistent pattern of treatment effects in which average heritability and additive genetic variance was lower in one pair and three pair lines than in 10 pair and control lines. Observed losses in heritability and additive genetic variance were significantly greater than predicted by the neutral additive model when calculated with coefficients of inbreeding estimated from demographic parameters alone. However, use of molecular markers revealed substantially more inbreeding, generated by increased variance in family size and background selection. Conservative interpretation of a statistical analysis incorporating this previously undetected inbreeding led to the conclusion that the response to inbreeding of the morphological traits studied showed no significant departure from the neutral additive model. This result is consistent with the evidence for minimal directional dominance for these traits. In contrast, egg hatching rate in the same experimental lines showed strong inbreeding depression, increased phenotypic variance and rapid response to selection, highly indicative of an increase in additive genetic variance due to dominance variance conversion.     

Inbreeding uncovers fundamental differences in the genetic load affecting male and female fertility in a butterfly

Inbreeding depression is most pronounced for traits closely associated with fitness. The traditional explanation is that natural selection eliminates deleterious mutations with additive or dominant effects more effectively than recessive mutations, leading to directional dominance for traits subject to strong directional selection. Here we report the unexpected finding that, in the butterfly Bicyclus anynana, male sterility contributes disproportionately to inbreeding depression for fitness (complete sterility in about half the sons from brother-sister matings), while female fertility is insensitive to inbreeding. The contrast between the sexes for functionally equivalent traits is inconsistent with standard selection arguments, and suggests that trait-specific developmental properties and cryptic selection play crucial roles in shaping genetic architecture. There is evidence that spermatogenesis is less developmentally stable than oogenesis, though the unusually high male fertility load in B. anynana additionally suggests the operation of complex selection maintaining male sterility recessives. Analysis of the precise causes of inbreeding depression will be needed to generate a model that reliably explains variation in directional dominance and reconciles the gap between observed and expected genetic loads carried by populations. This challenging evolutionary puzzle should stimulate work on the occurrence and causes of sex differences in fertility load.     

Dominance Genetic Variance for Traits Under Directional Selection in Drosophila serrata

In contrast to our growing understanding of patterns of additive genetic variance in single- and multi-trait combinations, the relative contribution of nonadditive genetic variance, particularly dominance variance, to multivariate phenotypes is largely unknown. While mechanisms for the evolution of dominance genetic variance have been, and to some degree remain, subject to debate, the pervasiveness of dominance is widely recognized and may play a key role in several evolutionary processes. Theoretical and empirical evidence suggests that the contribution of dominance variance to phenotypic variance may increase with the correlation between a trait and fitness; however, direct tests of this hypothesis are few. Using a multigenerational breeding design in an unmanipulated population of Drosophila serrata, we estimated additive and dominance genetic covariance matrices for multivariate wing-shape phenotypes, together with a comprehensive measure of fitness, to determine whether there is an association between directional selection and dominance variance. Fitness, a trait unequivocally under directional selection, had no detectable additive genetic variance, but significant dominance genetic variance contributing 32% of the phenotypic variance. For single and multivariate morphological traits, however, no relationship was observed between trait–fitness correlations and dominance variance. A similar proportion of additive and dominance variance was found to contribute to phenotypic variance for single traits, and double the amount of additive compared to dominance variance was found for the multivariate trait combination under directional selection. These data suggest that for many fitness components a positive association between directional selection and dominance genetic variance may not be expected.     

Genetics of quantitative traits in Arabidopsis thaliana

The genetic control of 22 quantitative traits, including developmental rates and sizes, was examined in generations of Arabidopsis thaliana derived from the cross between the ecotypes, Columbia (Col) and Landsberg erecta (Ler). The data were obtained from three sets of families raised in the same trial: the 16 basic generations, that is, parents, F(1)'s, F(2)'s, backcrosses, recombinant inbred lines (RILs) and a triple test cross (TTC), the latter produced by crossing the RILs to Col, Ler and their F(1). The data were analysed by two approaches. The first (approach A) involved traditional generation mean and variance component analysis and the second (B), based around the RILs and TTC families, involved marker-based QTL analysis.From (A), genetic differences between Col and Ler were detected for all traits with moderate heritabilities. Height at flowering was the only trait to show heterosis. Dominance was partial to complete for all height traits, and there was no overdominance but there was strong evidence for directional dominance. For most other traits, dominance was ambidirectional and incomplete, with average dominance ratios of around 80%. Epistasis, particularly of the duplicate type that opposes dominance, was a common feature of all traits. The presence of epistasis must imply multiple QTL for all traits.The QTL analysis located 38 significant effects in four regions of chromosomes I, II, IV and V, but not III. QTL affecting rosette size and leaf number were identified in all four regions, with days to maturity on chromosomes IV and V. The only QTL for height was located at the expected position of the erecta gene (chromosome II; 50 cM), but the additive and dominance effects of this single QTL did not adequately explain the generation means. The possible involvement of other interacting height QTL is discussed.     

Epistasis in maize (Zea mays L.)

Summary Three flint and three dent maize (Zea mays L.) inbred lines, their possible F₁ crosses, F₂ and backcross progenies, and all possible three-way crosses were evaluated in a three-year experiment for yield, ear moisture, and plant height. The purpose was to estimate genetic parameters in European breeding materials from (i) generation means analysis, (ii) diallel analysis of generation means, and (iii) analysis of F₁ and three-way cross hybrids. Method (i) was based on the F-metric model and methods (ii) and (iii) on the Eberhart-Gardner (1966) genetic model; both models extended for heterotic maternal effects.Differences among generation means for yield and plant height were mainly attributable to dominance effects. Epistatic effects were significantly different from zero in a few crosses and considerably reduced heterosis in both traits. Additive x additive and domiance x dominance effects for yield were consistently positive and negative, respectively. Significant maternal effects were established to the advantage of generations with a heterozygous seed parent. In the diallel analysis, mean squares for dominance effects were greater than for additive effects for yield and plant height but smaller for ear moisture. Though significant for yield and plant height, epistatic variation was small compared to additive and dominance variation. Estimates of additive x additive epistasis for yield were significantly negative in 11 of 15 crosses, suggesting that advantageous gene combinations in the lines had been disrupted by recombination in the segregating generations. The analysis of hybrids supported the above findings regarding the analysis of variance. However, the estimates of additive x additive epistasis for yield were considerably smaller and only minimally correlated with those from the diallel analysis. Use of noninbred materials as opposed to materials with different levels of inbreeding is considered the main reason for the discrepancies in the results.     

Variance component analysis of bristle characters in local populations of Drosophila melanogaster.

The genetic variabilities of sternopleural and abdominal bristle numbers existing in local natural populations were assessed. Using second chromosome lines of Drosophila melanogaster extracted from three natural populations in Japan (the Ishigakijima, Ogasawara and Aomori populations), experiments were conducted to estimate the components of genetic variances, additive and dominance variances. The following results were obtained: For both sternopleural and abdominal bristle numbers, the additive genetic variances (sigma 2A) were much larger than the dominance variances (sigma 2D) especially in the southern populations. For example, in the Ishigakijima population, for females sternopleural bristle numbers of the inversion-free chromosome group, the additive and dominance variances were estimated to be 1.255 +/- 0.2034 and 0.0552 +/- 0.0180, respectively. The magnitudes of the estimates of additive genetic variances were nearly equal from north to south. By comparing the additive genetic variances of the inversion-free chromosome group with those of the In(2L)t-carrying chromosome group, it was inferred that sufficient number of generations to achieve the equilibrium state has not passed since the introduction of a single or a small number of the In(2L)t-carrying chromosomes to the Ishigakijima population.     

A Study of the Genetic Basis of the Sexual Dimorphism for Wing Length in DROSOPHILA MELANOGASTER

The genetic basis of a sexually dimorphic quantitative character in Drosophila melanogaster was investigated by means of two-way directional selection for increased and decreased differences between male and female wing length. The sex dimorphism (SD), defined as the mean wing length difference between the sexes, within families, provided the criterion for selection.-The two lines (High SD, Low SD) diverged rapidly during the 15 generations of selection, indicating the presence of extensive genetic variability for the genotype-sex interaction underlying the observed sexual dimorphism. There was evidence that genetic variability persisted in both lines when selection was relaxed. Most of the divergence between the two lines remained after 10 generations of relaxed selection.-The change in the level of sex dimorphism in the High line was due primarily to a decrease in male wing length; in the Low line most of the change in SD was the result of a decrease in female wing length. An overall reduction in wing length in both sexes in both lines is interpreted as an effect of inbreeding.-The distribution and nature of the genetic control underlying the SD characteristic of the two selection lines was investigated by chromosome substitution between selection lines using a marked inversion technique. The two lines differed by factors located on each of the three major chromosome pairs. Chromosome III had the greatest effect on the difference in SD level between lines, and showed an overall additive effect when present in homozygous versus heterozygous combination. Chromosome II had the least effect, with a significant dominance effect of the High II being evident when heterozygotes were compared with homozygotes. The effect of the X chromosome was intermediate. There was some evidence of interaction between non-homologous chromosomes.     

Joint modeling of additive and non-additive (genetic line) effects in multi-environment trials

A statistical approach for the analysis of multi-environment trials (METs) is presented, in which selection of best performing lines, best parents, and best combination of parents can be determined. The genetic effect of a line is partitioned into additive, dominance and residual non-additive effects. The dominance effects are estimated through the incorporation of the dominance relationship matrix, which is presented under varying levels of inbreeding. A computationally efficient way of fitting dominance effects is presented which partitions dominance effects into between family dominance and within family dominance line effects. The overall approach is applicable to inbred lines, hybrid lines and other general population structures where pedigree information is available.     

The Genetic Basis of Natural Variation in Seed Size and Seed Number and Their Trade-Off Using Arabidopsis thaliana MAGIC Lines

Offspring number and size are key traits determining an individual’s fitness and a crop’s yield. Yet, extensive natural variation within species is observed for these traits. Such variation is typically explained by trade-offs between fecundity and quality, for which an optimal solution is environmentally dependent. Understanding the genetic basis of seed size and number, as well as any possible genetic constraints preventing the maximization of both, is crucial from both an evolutionary and applied perspective. We investigated the genetic basis of natural variation in seed size and number using a set of Arabidopsis thaliana multiparent advanced generation intercross (MAGIC) lines. We also tested whether life history affects seed size, number, and their trade-off. We found that both seed size and seed number are affected by a large number of mostly nonoverlapping QTL, suggesting that seed size and seed number can evolve independently. The allele that increases seed size at most identified QTL is from the same natural accession, indicating past occurrence of directional selection for seed size. Although a significant trade-off between seed size and number is observed, its expression depends on life-history characteristics, and generally explains little variance. We conclude that the trade-off between seed size and number might have a minor role in explaining the maintenance of variation in seed size and number, and that seed size could be a valid target for selection.     

The Expression of Additive and Nonadditive Genetic Variation under Stress

Experimental lines of Drosophila melanogaster derived from a natural population, which had been isolated in the laboratory for ~70 generations, were crossed to determine if the expression of additive, dominance and epistatic genetic variation in development time and viability was associated with the environment. No association was found between the level of additive genetic effects and environmental value for either trait, but nonadditive genetic effects increased at both extremes of the environmental range for development time. The expression of high levels of dominance and epistatic genetic variation at environmental extremes may be a general expectation for some traits. The disruption of the epistatic gene complexes in the parental lines resulted in hybrid breakdown toward faster development and there was some indication of hybrid breakdown toward higher viability. A combination of genetic drift and natural selection had therefore resulted in different epistatic gene complexes being selected after ~70 generations from a common genetic base. After crossing, the hybrid populations were observed for 10 generations. Epistasis contributed on average 12 hr in development time. Fluctuating asymmetry in sternopleural bristle number also evolved in the hybrid populations, decreasing by >18% in the first seven generations after hybridization.     

Analysis of genetic control of maize response in anther culture within a diallelic set

The androgenic ability in anther culture in vitro was examined within a diallel set of five maize lines. The complex analysis of genetic variation components, Hayman's diagram, and genetic parameters showed that the anther response is under the control of an additive-dominant genetic system. The examined lines possessed the different correlation of dominant and recessive alleles controlling androgenesis. And 44 was the line with the biggest number of recessive genes, which determined the increase in the trait. The level of dominance varied in different loci, though in the whole the degree of dominance approached to the complete one.     

Selection on skewed characters and the paradox of stasis

Observed phenotypic responses to selection in the wild often differ from predictions based on measurements of selection and genetic variance. An overlooked hypothesis to explain this paradox of stasis is that a skewed phenotypic distribution affects natural selection and evolution. We show through mathematical modeling that, when a trait selected for an optimum phenotype has a skewed distribution, directional selection is detected even at evolutionary equilibrium, where it causes no change in the mean phenotype. When environmental effects are skewed, Lande and Arnold's (1983) directional gradient is in the direction opposite to the skew. In contrast, skewed breeding values can displace the mean phenotype from the optimum, causing directional selection in the direction of the skew. These effects can be partitioned out using alternative selection estimates based on average derivatives of individual relative fitness, or additive genetic covariances between relative fitness and trait (Robertson–Price identity). We assess the validity of these predictions using simulations of selection estimation under moderate sample sizes. Ecologically relevant traits may commonly have skewed distributions, as we here exemplify with avian laying date — repeatedly described as more evolutionarily stable than expected — so this skewness should be accounted for when investigating evolutionary dynamics in the wild.