Uncovering QTL for resistance and survival time to Philasterides dicentrarchi in turbot (Scophthalmus maximus)

Disease resistance-related traits have received increasing importance in aquaculture breeding programs worldwide. Currently, genomic information offers new possibilities in breeding to address the improvement of this kind of traits. The turbot is one of the most promising European aquaculture species, and Philasterides dicentrarchi is a scuticociliate parasite causing fatal disease in farmed turbot. An appealing approach to fight against disease is to achieve a more robust broodstock, which could prevent or diminish the devastating effects of scuticociliatosis on farmed individuals. In the present study, a genome scan for quantitative trait loci (QTL) affecting resistance and survival time to P. dicentrarchi in four turbot families was carried out. The objectives were to identify QTL using different statistical approaches [linear regression (LR) and maximum likelihood (ML)] and to locate significantly associated markers for their application in genetic breeding strategies. Several genomic regions controlling resistance and survival time to P. dicentrarchi were detected. When analyzing each family separately, significant QTL for resistance were identified by the LR method in two linkage groups (LG1 and LG9) and for survival time in LG1, while the ML methodology identified QTL for resistance in LG9 and LG23 and for survival time in LG6 and LG23. The analysis of the total data set identified an additional significant QTL for resistance and survival time in LG3 with the LR method. Significant association between disease resistance-related traits and genotypes was detected for several markers, a single one explaining up to 22% of the phenotypic variance. Obtained results will be essential to identify candidate genes for resistance and to apply them in marker-assisted selection programs to improve turbot production.

QTL detection for Aeromonas salmonicida resistance related traits in turbot (Scophthalmus maximus)

Background

Interactions between fish and pathogens, that may be harmless under natural conditions, often result in serious diseases in aquaculture systems. This is especially important due to the fact that the strains used in aquaculture are derived from wild strains that may not have had enough time to adapt to new disease pressures. The turbot is one of the most promising European aquaculture species. Furunculosis, caused by the bacterium Aeromonas salmonicida, produces important losses to turbot industry. An appealing solution is to achieve more robust broodstock, which can prevent or diminish the devastating effects of epizooties. Genomics strategies have been developed in turbot to look for candidate genes for resistance to furunculosis and a genetic map with appropriate density to screen for genomic associations has been also constructed. In the present study, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out. The objectives were to identify consistent QTL using different statistical approaches (linear regression and maximum likelihood) and to locate the tightest associated markers for their application in genetic breeding strategies.

Results

Significant QTL for resistance were identified by the linear regression method in three linkage groups (LGs 4, 6 and 9) and for survival in two LGs (6 and 9). The maximum likelihood methodology identified QTL in three LGs (5, 6 and 9) for both traits. Significant association between disease traits and genotypes was detected for several markers, some of them explaining up to 17% of the phenotypic variance. We also identified candidate genes located in the detected QTL using data from previously mapped markers.

Conclusions

Several regions controlling resistance to A. salmonicida in turbot have been detected. The observed concordance between different statistical methods at particular linkage groups gives consistency to our results. The detected associated markers could be useful for genetic breeding strategies. A finer mapping will be necessary at the detected QTL intervals to narrow associations and around the closely associated markers to look for candidate genes through comparative genomics or positional cloning strategies. The identification of associated variants at specific genes will be essential, together with the QTL associations detected in this study, for future marker assisted selection programs.

Detection of growth-related QTL in turbot (Scophthalmus maximus)

Background

The turbot (Scophthalmus maximus) is a highly appreciated European aquaculture species. Growth related traits constitute the main goal of the ongoing genetic breeding programs of this species. The recent construction of a consensus linkage map in this species has allowed the selection of a panel of 100 homogeneously distributed markers covering the 26 linkage groups (LG) suitable for QTL search. In this study we addressed the detection of QTL with effect on body weight, length and Fulton's condition factor.

Results

Eight families from two genetic breeding programs comprising 814 individuals were used to search for growth related QTL using the panel of microsatellites available for QTL screening. Two different approaches, maximum likelihood and regression interval mapping, were used in order to search for QTL. Up to eleven significant QTL were detected with both methods in at least one family: four for weight on LGs 5, 14, 15 and 16; five for length on LGs 5, 6, 12, 14 and 15; and two for Fulton's condition factor on LGs 3 and 16. In these LGs an association analysis was performed to ascertain the microsatellite marker with the highest apparent effect on the trait, in order to test the possibility of using them for marker assisted selection.

Conclusions

The use of regression interval mapping and maximum likelihood methods for QTL detection provided consistent results in many cases, although the high variation observed for traits mean among families made it difficult to evaluate QTL effects. Finer mapping of detected QTL, looking for tightly linked markers to the causative mutation, and comparative genomics are suggested to deepen in the analysis of QTL in turbot so they can be applied in marker assisted selection programs.

Restricting inbreeding while maintaining selection response for weight gain in Mus musculus

An experiment with mice was designed to test the relative efficiency of three selection methods that help to minimize the rate of inbreeding during selection. A common house mice (Mus musculus) population was selected for 17 generations to increase the weight gain between 21 and 42 days. The population was split at random into three lines A, B and C where three selection methods were applied: individual selection and random mating, weighted selection with random mating and individual selection with minimum coancestry mating, respectively. There were three replicates for each line. Cumulated selection response was similar in the three lines, but there were differences in the level of inbreeding attained (in percentage): 31.24 (method A), 24.72 (method B) and 27.88 (method C). As consequence, lines B and C (weighted selection and minimum coancestry) showed a lower value of deterioration of fitness traits (the intrauterine mortality and the mortality at birth) than line A (random mating).

Accuracy of genome-wide evaluation for disease resistance in aquaculture breeding programs

Current aquaculture breeding programs aimed at improving resistance to diseases are based on challenge tests, where performance is recorded on sibs of candidates to selection, and on selection between families. Genome-wide evaluation (GWE) of breeding values offers new opportunities for using variation within families when dealing with such traits. However, up-to-date studies on GWE in aquaculture programs have only considered continuous traits. The objectives of this study were to extend GWE methodology, in particular the Bayes B method, to analyze dichotomous traits such as resistance to disease, and to quantify, through computer simulation, the accuracy of GWE for disease resistance in aquaculture sib-based programs, using the methodology developed. Two heritabilities (0.1 and 0.3) and 2 disease prevalences (0.1 and 0.5) were assumed in the simulations. We followed the threshold liability model, which assumes that there is an underlying variable (liability) with a continuous distribution and assumed a BayesB model for the liabilities. It was shown that the threshold liability model used fits very well with the BayesB model of GWE. The advantage of using the threshold model was clear when dealing with disease resistance dichotomous phenotypes, particularly under the conditions where linear models are less appropriate (low heritability and disease prevalence). In the testing set (where individuals are genotyped but not measured), the increase in accuracy for the simulated schemes when using the threshold model ranged from 4 (for heritability equal to 0.3 and prevalence equal to 0.5) to 16% (for heritability and prevalence equal to 0.1) when compared with the linear model.

A note on mate allocation for dominance handling in genomic selection

Estimation of non-additive genetic effects in animal breeding is important because it increases the accuracy of breeding value prediction and the value of mate allocation procedures. With the advent of genomic selection these ideas should be revisited. The objective of this study was to quantify the efficiency of including dominance effects and practising mating allocation under a whole-genome evaluation scenario. Four strategies of selection, carried out during five generations, were compared by simulation techniques. In the first scenario (MS), individuals were selected based on their own phenotypic information. In the second (GSA), they were selected based on the prediction generated by the Bayes A method of whole-genome evaluation under an additive model. In the third (GSD), the model was expanded to include dominance effects. These three scenarios used random mating to construct future generations, whereas in the fourth one (GSD + MA), matings were optimized by simulated annealing. The advantage of GSD over GSA ranges from 9 to 14% of the expected response and, in addition, using mate allocation (GSD + MA) provides an additional response ranging from 6% to 22%. However, mate selection can improve the expected genetic response over random mating only in the first generation of selection. Furthermore, the efficiency of genomic selection is eroded after a few generations of selection, thus, a continued collection of phenotypic data and re-evaluation will be required.

Microsatellite-based genetic diversity and population structure of domestic sheep in northern Eurasia

Background

Identification of global livestock diversity hotspots and their importance in diversity maintenance is essential for making global conservation efforts. We screened 52 sheep breeds from the Eurasian subcontinent with 20 microsatellite markers. By estimating and weighting differently within- and between-breed genetic variation our aims were to identify genetic diversity hotspots and prioritize the importance of each breed for conservation, respectively. In addition we estimated how important within-species diversity hotspots are in livestock conservation.

Results

Bayesian clustering analysis revealed three genetic clusters, termed Nordic, Composite and Fat-tailed. Southern breeds from close to the region of sheep domestication were more variable, but less genetically differentiated compared with more northern populations. Decreasing weight for within-breed diversity component led to very high representation of genetic clusters or regions containing more diverged breeds, but did not increase phenotypic diversity among the high ranked breeds. Sampling populations throughout 14 regional groups was suggested for maximized total genetic diversity.

Conclusions

During initial steps of establishing a livestock conservation program populations from the diversity hot-spot area are the most important ones, but for the full design our results suggested that approximately equal population presentation across environments should be considered. Even in this case, higher per population emphasis in areas of high diversity is appropriate. The analysis was based on neutral data, but we have no reason to think the general trend is limited to this type of data. However, a comprehensive valuation of populations should balance production systems, phenotypic traits and available genetic information, and include consideration of probability of success.

Objectives, criteria and methods for using molecular genetic data in priority setting for conservation of animal genetic resources

The genetic diversity of the world's livestock populations is decreasing, both within and across breeds. A wide variety of factors has contributed to the loss, replacement or genetic dilution of many local breeds. Genetic variability within the more common commercial breeds has been greatly decreased by selectively intense breeding programmes. Conservation of livestock genetic variability is thus important, especially when considering possible future changes in production environments. The world has more than 7500 livestock breeds and conservation of all of them is not feasible. Therefore, prioritization is needed. The objective of this article is to review the state of the art in approaches for prioritization of breeds for conservation, particularly those approaches that consider molecular genetic information, and to identify any shortcomings that may restrict their application. The Weitzman method was among the first and most well-known approaches for utilization of molecular genetic information in conservation prioritization. This approach balances diversity and extinction probability to yield an objective measure of conservation potential. However, this approach was designed for decision making across species and measures diversity as distinctiveness. For livestock, prioritization will most commonly be performed among breeds within species, so alternatives that measure diversity as co-ancestry (i.e. also within-breed variability) have been proposed. Although these methods are technically sound, their application has generally been limited to research studies; most existing conservation programmes have effectively primarily based decisions on extinction risk. The development of user-friendly software incorporating these approaches may increase their rate of utilization.

Genetic diversity in farm animals--a review

Domestication of livestock species and a long history of migrations, selection and adaptation have created an enormous variety of breeds. Conservation of these genetic resources relies on demographic characterization, recording of production environments and effective data management. In addition, molecular genetic studies allow a comparison of genetic diversity within and across breeds and a reconstruction of the history of breeds and ancestral populations. This has been summarized for cattle, yak, water buffalo, sheep, goats, camelids, pigs, horses, and chickens. Further progress is expected to benefit from advances in molecular technology.

Assessing population genetic structure via the maximisation of genetic distance

Background

The inference of the hidden structure of a population is an essential issue in population genetics. Recently, several methods have been proposed to infer population structure in population genetics.

Methods

In this study, a new method to infer the number of clusters and to assign individuals to the inferred populations is proposed. This approach does not make any assumption on Hardy-Weinberg and linkage equilibrium. The implemented criterion is the maximisation (via a simulated annealing algorithm) of the averaged genetic distance between a predefined number of clusters. The performance of this method is compared with two Bayesian approaches: STRUCTURE and BAPS, using simulated data and also a real human data set.

Results

The simulations show that with a reduced number of markers, BAPS overestimates the number of clusters and presents a reduced proportion of correct groupings. The accuracy of the new method is approximately the same as for STRUCTURE. Also, in Hardy-Weinberg and linkage disequilibrium cases, BAPS performs incorrectly. In these situations, STRUCTURE and the new method show an equivalent behaviour with respect to the number of inferred clusters, although the proportion of correct groupings is slightly better with the new method. Re-establishing equilibrium with the randomisation procedures improves the precision of the Bayesian approaches. All methods have a good precision for F_ST ≥ 0.03, but only STRUCTURE estimates the correct number of clusters for F_ST as low as 0.01. In situations with a high number of clusters or a more complex population structure, MGD performs better than STRUCTURE and BAPS. The results for a human data set analysed with the new method are congruent with the geographical regions previously found.

Conclusion

This new method used to infer the hidden structure in a population, based on the maximisation of the genetic distance and not taking into consideration any assumption about Hardy-Weinberg and linkage equilibrium, performs well under different simulated scenarios and with real data. Therefore, it could be a useful tool to determine genetically homogeneous groups, especially in those situations where the number of clusters is high, with complex population structure and where Hardy-Weinberg and/or linkage equilibrium are present.

Preserving population allele frequencies in ex situ conservation programs

Optimization of contributions of parents to progeny by minimizing the average coancestry of the progeny is an effective strategy for maintaining genetic diversity in ex situ conservation programs, but its application on the basis of molecular markers has the negative collateral effect of homogenizing the allelic frequencies at each locus. Because one of the objectives of a conservation program is to preserve the genetic composition of the original endangered population, we devised a method in which markers are used to maintain the allele frequency distribution at each locus as closely as possible to that of the native population. Contributions of parents were obtained so as to minimize changes in allele frequency for a set of molecular markers in a population of reduced size. We used computer simulations, under a range of scenarios, to assess the effectiveness of the method in comparison with methods in which contributions of minimum coancestry are sought, either making use of molecular markers or genealogical information. Our simulations indicated that the proposed method effectively maintained the original distribution of allele frequencies, particularly under strong linkage, and maintained acceptable levels of genetic diversity in the population. Nevertheless, contributions of minimum coancestry determined from pedigree information but ignoring the genealogy previous to the conservation program, was the most effective method for maintaining allelic frequencies in realistic situations.

A quantitative trait locus genome scan for porcine muscle fiber traits reveals overdominance and epistasis

Muscle histochemical characteristics are decisive determinants of meat quality. The relative percentage and diameters of the different muscular fiber types influence crucial aspects of meat such as color, tenderness, and ultimate pH. Despite its relevance, however, the information on muscle fiber genetic architecture is scant, because histochemical muscle characterization is a laborious task. Here we report a complete QTL scan of muscle fiber traits in 160 animals from a F(2) cross between Iberian and Landrace pigs using 139 markers. We identified 20 genome regions distributed along 15 porcine chromosomes (SSC1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, and X) with direct and(or) epistatic effects. Epistasis was frequent and some interactions were highly significant. Chromosomes 10 and 11 seemed to behave as hubs; they harbored 2 individual QTL, but also 6 epistatic regions. Numerous individual QTL effects had cryptic alleles, with opposite effects to phenotypic pure breed differences. Many of the QTL identified here coincided with previous reports for these traits in the literature, and there was overlapping with potential candidate genes and previously reported meat quality QTL.

Iterated prisoner's dilemma in an asocial world dominated by loners, not by defectors

Cooperation among genetically unrelated individuals can arise when pairs of individuals interact repeatedly in the Prisoner's Dilemma. However, the conditions allowing the evolution of reciprocal cooperation become extremely restrictive as the size of the cooperative group increases, because defectors can exploit cooperators more efficiently in larger groups. Here we consider three strategies: Tit for Tat, defector, and loner. Loner beats defector in a non-cooperative world. However, a cooperative strategy Tit for Tat (TFT(0)) that stops cooperation after the first iteration when there is at least one defector in the group, can invade a world of loners, even in sizable groups, if both the TFT(0) and the defector strategies arise at the same frequency by mutation.

Assortative mating through a mechanism of sexual selection

We propose that assortative mating can arise through a mechanism of sexual selection by active female choice of partners based on a 'self-seeking like' decision rule. Using a mathematical model, we show that a gene can be selected that make females to choose mates that are similar to themselves with respect to an arbitrary tag, even if two independent and unlinked genes determine the preference and the tag. The necessary requisite for this process to apply is an asymmetry between partners, such that the female can choose the male but this one must always accept to mate. The fitness advantage is due to linkage disequilibrium built up between both genes. Simulations have been run to check the algebraic results and to analyse the influence of several factors on the evolution of the system. Any factor that favors linkage disequilibrium also favors the evolution of the preference allele. Moreover, in a large population subdivided in small subpopulations connected by migration, the assortative mating homogenizes the population genotypic structure for the tags in contrast to random mating that maintains most of the variation.

Mutual benefit cooperation and ethnic cultural diversity

We show that mutual benefit cooperation can favor the evolution of a preference to interact with individuals that are similar to themselves with respect to one or several arbitrary cultural tags. One necessary requisite to obtain this result is an asymmetry between partners in such a way that one of them (actor) proposes the cooperation and elects the partner, whereas the other (receiver) never rejects the offer because cooperation always reports benefits. The advantage of individuals possessing allele for preferential assortment is due to receiving more offers of mutually beneficial cooperation when there is linkage disequilibrium between the assortment locus and the cultural tags. An especially favorable scenario for the evolution of such preference is a subdivided metapopulation. In this case, the homogeneity within populations and the divergence between populations is favored, facilitating the existence of ethnic groups.

A new method to estimate relatedness from molecular markers

Four major problems can affect the efficiency of methods developed to estimate relatedness between individuals from information of molecular markers: (i) some of them are dependent on the knowledge of the true allelic frequencies in the base population; (ii) they assume that all loci are unlinked and in Hardy-Weinberg and linkage equilibrium; (iii) pairwise methods can lead to incongruous assignations because they take into account only two individuals at a time; (iv) most are usually constructed for particular structured populations (only consider a few relationship classes, e.g. full-sibs vs. unrelated). We have developed a new approach to estimate relatedness that is free from the above limitations. The method uses a 'blind search algorithm' (actually simulated annealing) to find the genealogy that yield a co-ancestry matrix with the highest correlation with the molecular co-ancestry matrix calculated using the markers. Thus (i and ii) it makes no direct assumptions about allelic frequencies or Hardy-Weinberg and linkage equilibrium; (iii) it always provide congruent relationships, as it considers all individuals at a time; (iv) degrees of relatedness can be as complex as desired just increasing the 'depth' (i.e. number of generations) of the proposed genealogies. Computer simulations have shown that the accuracy and robustness against genotyping errors of this new approach is comparable to that of other proposed methods in those particular situations they were developed for, but it is more flexible and can cope with more complex situations.

Sex ratio variation in Iberian pigs

Within the area of sex allocation, one of the topics that has attracted a lot of attention is the sex ratio problem. Fisher (1930) proposed that equal numbers of males and females have been promoted by natural selection and it has an adaptive significance. But the empirical success of Fisher's theory remains doubtful because a sex ratio of 0.50 is also expected from the chromosomal mechanism of sex determination. Another way of approaching the subject is to consider that Fisher's argument relies on the underlying assumption that offspring inherit their parent's tendency in biased sex ratio and therefore that genetic variance for this trait exists. Here, we analyzed sex ratio data of 56,807 piglets coming from 550 boars and 1893 dams. In addition to classical analysis of heterogeneity we performed analyses fitting linear and threshold animal models in a Bayesian framework using Gibbs sampling techniques. The marginal posterior mean of heritability was 2.63 x 10(-4) under the sire linear model and 9.17 x 10(-4) under the sire threshold model. The probability of the hypothesis p(h(2) = 0) fitting the last model was 0.996. Also, we did not detect any trend in sex ratio related to maternal age. From an evolutionary point of view, the chromosomal sex determination acts as a constraint that precludes control of offspring sex ratio in vertebrates and it should be included in the general theory of sex allocation. From a practical view that means that the sex ratio in domestic species is hardly susceptible to modification by artificial selection.

The effect of genetic drift on the variance/covariance components generated by multilocus additive x additive epistatic systems

The effect of population bottlenecks on the components of the genetic variance/covariance generated by n neutral independent additive x additive loci has been studied theoretically. In its simplest version, this situation can be modelled by specifying the allele frequencies and homozygous effects at each locus, and an additional factor measuring the strength of the n-th order epistatic interaction. The variance/covariance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium over replicates randomly derived from the base population, after t bottlenecks of size N (derived components). Formulae were obtained giving the derived components (and the between-line variance) as functions of the ancestral ones (alternatively, in terms of allele frequencies and effects) and the corresponding inbreeding coefficient F(t). The n-th order derived component of the genetic variance/covariance is continuously eroded by inbreeding, but the remaining components may increase initially until a critical F(t) value is attained, which is inversely related to the order of the pertinent component, and subsequently decline to zero. These changes can be assigned to the between-line variances/covariances of gene substitution and epistatic effects induced by drift. Numerical examples indicate that: (1) the derived additive variance/covariance component will generally exceed its ancestral value unless epistasis is weak; (2) the derived epistatic variance/covariance components will generally exceed their ancestral values unless allele frequencies are extreme; (3) for systems showing equal ancestral additive and total non-additive variance/covariance components, those including a smaller number of epistatic loci may generate a larger excess in additive variance/covariance after bottlenecks than others involving a larger number of loci, provided that F(t) is low. Our results indicate that it is unlikely that the rate of evolution may be significantly accelerated after population bottlenecks, in spite of occasional increments of the derived additive variance over its ancestral value.

Benefits from marker-assisted selection under an additive polygenic genetic model1

This study investigated, through stochastic computer simulation, the extra gains expected from marker-assisted selection (MAS) in an infinitesimal model with linkage. The trait under selection was assumed to be controlled by 2,000 loci of additive small effect and evenly distributed in c chromosomes of one Morgan each (and c = 5, 10, 20, or 30). This approach differs from previous studies on the benefits of MAS that have considered mixed inheritance models. Marker information was used together with pedigree information to compute the relationship matrix used in BLUP genetic evaluations. The MAS schemes were compared with schemes where genetic evaluations were performed using standard BLUP (i.e., the relationship matrix is obtained using pedigree information only). When the number of markers was large enough (approximately one marker every 10 cM), there were increases in the accuracy of selection with MAS, and this led to extra gains compared with standard BLUP for all genome sizes considered. The benefit from MAS increased over generations. At the last generation of selection (Generation 10), the response from MAS was 11, 9, 7, and 5% greater than with standard BLUP for genomes with 5, 10, 20, and 30 chromosomes, respectively. Thus, although small, gains from MAS were nonetheless detectable for genome sizes typical of livestock populations.

Inferences on the role of insertion in a mutation accumulation experiment with Drosophila melanogaster using RAPDs

The genetic variability for RAPDs band pattern was studied in a set of 157 mutation accumulation (MA) lines of Drosophila melanogaster. These MA lines were derived from the same isogenic base population and subsequently maintained by full-sib mating during 132 generations. The ancestral pattern of the original isogenic base can be unambiguously established as the consensus pattern of the MA lines and, because these lines are expected to be homozygous, dominance for band pattern is not a concern. Only repeatable changes in band pattern were considered. The number of ancestral bands detected implies that nine-nucleotide targets are enough for repeatable PCR amplification. Compared with the ancestral pattern, one MA line lost one band and two MA lines gained a new one. These results can be accounted for by the insertion of transposable elements occurring at a rate 0.07 < i < 0.21 per whole haploid genome and generation. This range is typical for Drosophila and consistent with the previously observed mobility for the roo family, supporting the generality of previous estimates of spontaneous mutation rates for morphological and fitness traits based on these MA lines. The sequence of one of the new bands suggests that the Idefix family is also active in the lines.

Mutual benefit can promote the evolution of preferential interactions and in this way can lead to the evolution of true altruism

We analyse the evolution of the assortment of encounters through active choice of companions among individuals that interact cooperatively in a situation of mutual benefit. Using a simple mathematical model, we show that mutual benefit can favour the evolution of a preference to interact with individuals that are similar to themselves with respect to an arbitrary tag even when both the preference and the tag depend on two independent and unlinked genes. Two necessary requisites to obtain this result are: (i) a small population or a large subdivided metapulation and (ii) an asymmetry between partners in such a way that one of them (donor) proposes the cooperation and elects the partner, whereas the other (receiver) never rejects the offer. We also show that mutual benefit can be the starting point for the evolution of altruistic behaviours as long as there are preferential interactions. This requires that the tag used in the election of partners is the altruistic or selfish behaviour itself.

EPISTASIS AND THE TEMPORAL CHANGE IN THE ADDITIVE VARIANCE-COVARIANCE MATRIX INDUCED BY DRIFT

The effect of population bottlenecks on the components of the genetic covariance generated by two neutral independent epistatic loci has been studied theoretically (additive, covA; dominance, covD; additive-by-additive, covAA; additive-by-dominance, covAD; and dominance-by-dominance, covDD). The additive-by-additive model and a more general model covering all possible types of marginal gene action at the single-locus level (additive/dominance epistatic model) were considered. The covariance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium over replicates randomly derived from the base population, after t consecutive bottlenecks of equal size N (derived components). Formulae were obtained in terms of the allele frequencies and effects at each locus, the corresponding epistatic effects and the inbreeding coefficient Ft. These expressions show that the contribution of nonadditive loci to the derived additive covariance (covAt) does not linearly decrease with inbreeding, as in the pure additive case, and may initially increase or even change sign in specific situations. Numerical examples were also analyzed, restricted for simplicity to the case of all covariance components being positive. For additive-by-additive epistasis, the condition covAt > covA only holds for high frequencies of the allele decreasing the metric traits at each locus (negative allele) if epistasis is weak, or for intermediate allele frequencies if it is strong. For the additive/dominance epistatic model, however, covAt > covA applies for low frequencies of the negative alleles at one or both loci and mild epistasis, but this result can be progressively extended to intermediate frequencies as epistasis becomes stronger. Without epistasis the same qualitative results were found, indicating that marginal dominance induced by epistasis can be considered as the primary cause of an increase of the additive covariance after bottlenecks. For all models, the magnitude of the ratio covAt/covA was inversely related to N and t.