A comparison of alternative methods to compute conditional genotype probabilities for genetic evaluation with finite locus models

An increased availability of genotypes at marker loci has prompted the development of models that include the effect of individual genes. Selection based on these models is known as marker-assisted selection (MAS). MAS is known to be efficient especially for traits that have low heritability and non-additive gene action. BLUP methodology under non-additive gene action is not feasible for large inbred or crossbred pedigrees. It is easy to incorporate non-additive gene action in a finite locus model. Under such a model, the unobservable genotypic values can be predicted using the conditional mean of the genotypic values given the data. To compute this conditional mean, conditional genotype probabilities must be computed. In this study these probabilities were computed using iterative peeling, and three Markov chain Monte Carlo (MCMC) methods – scalar Gibbs, blocking Gibbs, and a sampler that combines the Elston Stewart algorithm with iterative peeling (ESIP). The performance of these four methods was assessed using simulated data. For pedigrees with loops, iterative peeling fails to provide accurate genotype probability estimates for some pedigree members. Also, computing time is exponentially related to the number of loci in the model. For MCMC methods, a linear relationship can be maintained by sampling genotypes one locus at a time. Out of the three MCMC methods considered, ESIP, performed the best while scalar Gibbs performed the worst.     

Estimating Additive and Non-Additive Genetic Variances and Predicting Genetic Merits Using Genome-Wide Dense Single Nucleotide Polymorphism Markers

Non-additive genetic variation is usually ignored when genome-wide markers are used to study the genetic architecture and genomic prediction of complex traits in human, wild life, model organisms or farm animals. However, non-additive genetic effects may have an important contribution to total genetic variation of complex traits. This study presented a genomic BLUP model including additive and non-additive genetic effects, in which additive and non-additive genetic relation matrices were constructed from information of genome-wide dense single nucleotide polymorphism (SNP) markers. In addition, this study for the first time proposed a method to construct dominance relationship matrix using SNP markers and demonstrated it in detail. The proposed model was implemented to investigate the amounts of additive genetic, dominance and epistatic variations, and assessed the accuracy and unbiasedness of genomic predictions for daily gain in pigs. In the analysis of daily gain, four linear models were used: 1) a simple additive genetic model (MA), 2) a model including both additive and additive by additive epistatic genetic effects (MAE), 3) a model including both additive and dominance genetic effects (MAD), and 4) a full model including all three genetic components (MAED). Estimates of narrow-sense heritability were 0.397, 0.373, 0.379 and 0.357 for models MA, MAE, MAD and MAED, respectively. Estimated dominance variance and additive by additive epistatic variance accounted for 5.6% and 9.5% of the total phenotypic variance, respectively. Based on model MAED, the estimate of broad-sense heritability was 0.506. Reliabilities of genomic predicted breeding values for the animals without performance records were 28.5%, 28.8%, 29.2% and 29.5% for models MA, MAE, MAD and MAED, respectively. In addition, models including non-additive genetic effects improved unbiasedness of genomic predictions.     

The Genetic Correlation between Height and IQ: Shared Genes or Assortative Mating?

Traits that are attractive to the opposite sex are often positively correlated when scaled such that scores increase with attractiveness, and this correlation typically has a genetic component. Such traits can be genetically correlated due to genes that affect both traits (“pleiotropy”) and/or because assortative mating causes statistical correlations to develop between selected alleles across the traits (“gametic phase disequilibrium”). In this study, we modeled the covariation between monozygotic and dizygotic twins, their siblings, and their parents (total N = 7,905) to elucidate the nature of the correlation between two potentially sexually selected traits in humans: height and IQ. Unlike previous designs used to investigate the nature of the height–IQ correlation, the present design accounts for the effects of assortative mating and provides much less biased estimates of additive genetic, non-additive genetic, and shared environmental influences. Both traits were highly heritable, although there was greater evidence for non-additive genetic effects in males. After accounting for assortative mating, the correlation between height and IQ was found to be almost entirely genetic in nature. Model fits indicate that both pleiotropy and assortative mating contribute significantly and about equally to this genetic correlation.     

Genome-wide detection of non-additive quantitative trait loci for semen production traits in beef and dairy bulls

Semen production traits are important aspects of bull fertility, because semen quantity leads to direct profits for artificial insemination centres, and semen quality is associated with the probability of achieving a pregnancy. Most genome-wide association studies (GWASs) for semen production traits have assumed that each quantitative trait locus (QTL) has an additive effect. However, GWASs that account for non-additive effects are also important in fitness traits, such as bull fertility. Here, we performed a GWAS using models that accounted for additive and non-additive effects to evaluate the importance of non-additive effects on five semen production traits in beef and dairy bulls. A total of 65 463 records for 615 Japanese Black bulls (JB) and 50 734 records for 873 Holstein bulls (HOL), which were previously genotyped using the Illumina BovineSNP50 BeadChip, were used to estimate genetic parameters and perform GWAS. The heritability estimates were low (ranged from 0.11 to 0.23), and the repeatability estimates were low to moderate (ranged from 0.28 to 0.45) in both breeds. The estimated repeatability was approximately twice as high as the estimated heritability for all traits. In this study, only one significant region with an additive effect was detected in each breed, but multiple significant regions with non-additive effects were detected for each breed. In particular, the region at approximately 64 Mbp on Bos taurus autosome 17 had the highest significant non-additive effect on four semen production traits in HOL. The rs41843851 single nucleotide polymorphism (SNP) in the region had a much lower P-value for the non-additive effect (P-value = 1.1 × 10^?31) than for the additive effect (P-value = 1.1 × 10^?8) in sperm motility. The AA and AB genotypes on the SNP had a higher phenotype than the BB genotype in HOL, and there was no bull with the BB genotype in JB. Our results showed that non-additive QTLs affect semen production traits, and a novel QTL accounting for non-additive effects could be detected by GWAS. This study provides new insights into non-additive QTLs that affect fitness traits, such as semen production traits in beef and dairy bulls.     

Genetically based population divergence in overwintering energy mobilization in brook charr (Salvelinus fontinalis)

Investigating the nature of physiological traits potentially related to fitness is important towards a better understanding of how species and/or populations may respond to selective pressures imposed by contrasting environments. In northern species in particular, the ability to mobilize energy reserves to compensate for the low external energy intake during winter is crucial. However, the phenotypic and genetic bases of energy reserve accumulation and mobilization have rarely been investigated, especially pertaining to variation in strategy adopted by different populations. In the present study, we documented variation in several energy reserve variables and estimated their quantitative genetic basis to test the null hypothesis of no difference in variation at those traits among three strains of brook charr (Salvelinus fontinalis) and their reciprocal hybrids. Our results indicate that the strategy of winter energy preparation and mobilization was specific to each strain, whereby (1) domestic fish accumulated a higher amount of energy reserves before winter and kept accumulating liver glycogen during winter despite lower feeding; (2) Laval fish used liver glycogen and lipids during winter and experienced a significant decrease in condition factor; (3) Rupert fish had relatively little energy reserves accumulated at the end of fall and preferentially mobilized visceral fat during winter. Significant heritability for traits related to the accumulation and use of energy reserves was found in the domestic and Laval but not in the Rupert strain. Genetic and phenotypic correlations also varied among strains, which suggested population-specific genetic architecture underlying the expression of these traits. Hybrids showed limited evidence of non-additive effects. Overall, this study provides the first evidence of a genetically based—and likely adaptive—population-specific strategy for energy mobilization related to overwinter survival.     

Heritability and evolvability of fitness and nonfitness traits: Lessons from livestock

Data from natural populations have suggested a disconnection between trait heritability (variance standardized additive genetic variance, V_A) and evolvability (mean standardized V_A) and emphasized the importance of environmental variation as a determinant of trait heritability but not evolvability. However, these inferences are based on heterogeneous and often small datasets across species from different environments. We surveyed the relationship between evolvability and heritability in >100 traits in farmed cattle, taking advantage of large sample sizes and consistent genetic approaches. Heritability and evolvability estimates were positively correlated (r = 0.37/0.54 on untransformed/log scales) reflecting a substantial impact of V_A on both measures. Furthermore, heritabilities and residual variances were uncorrelated. The differences between this and previously described patterns may reflect lower environmental variation experienced in farmed systems, but also low and heterogeneous quality of data from natural populations. Similar to studies on wild populations, heritabilities for life‐history and behavioral traits were lower than for other traits. Traits having extremely low heritabilities and evolvabilities (17% of the studied traits) were almost exclusively life‐history or behavioral traits, suggesting that evolutionary constraints stemming from lack of genetic variability are likely to be most common for classical “fitness” (cf. life‐history) rather than for “nonfitness” (cf. morphological) traits.     

Genomic rearrangements and gene copy-number alterations as a cause of nervous system disorders

Genomic disorders are a group of human genetic diseases caused by genomic rearrangements resulting in copy-number variation (CNV) affecting a dosage-sensitive gene or genes critical for normal development or maintenance. These disorders represent a wide range of clinically distinct entities but include many diseases affecting nervous system function. Herein, we review selected neurodevelopmental, neurodegenerative, and psychiatric disorders either known or suggested to be caused by genomic rearrangement and CNV. Further, we emphasize the cause-and-effect relationship between gene CNV and complex disease traits. We also discuss the prevalence and heritability of CNV, the correlation between CNV and higher-order genome architecture, and the heritability of personality, behavioral, and psychiatric traits. We speculate that CNV could underlie a significant proportion of normal human variation including differences in cognitive, behavioral, and psychological features.     

Inheritance of vertebral number in the three-spined stickleback (Gasterosteus aculeatus)

Intraspecific variation in the number of vertebrae is taxonomically widespread, and both genetic and environmental factors are known to contribute to this variation. However, the relative importance of genetic versus environmental influences on variation in vertebral number has seldom been investigated with study designs that minimize bias due to non-additive genetic and maternal influences. We used a paternal half-sib design and animal model analysis to estimate heritability and causal components of variance in vertebral number in three-spined sticklebacks (Gasterosteus aculeatus). We found that both the number of vertebrae (h² = 0.36) and body size (h² = 0.42) were moderately heritable, whereas the influence of maternal effects was estimated to be negligible. While the number of vertebrae had a positive effect on body size, no evidence for a genetic correlation between body size and vertebral number was detected. However, there was a significant positive environmental correlation between these two traits. Our results support the generalization-in accordance with results from a review of heritability estimates for vertebral number in fish, reptiles and mammals-that the number of vertebrae appears to be moderately to highly heritable in a wide array of species. In the case of the three-spined stickleback, independent evolution of body size and number of vertebrae should be possible given the low genetic correlation between the two traits.     

Sex-specific genetic variances in life-history and morphological traits of the seed beetle Callosobruchus maculatus

Knowledge of heritability and genetic correlations are of central importance in the study of adaptive trait evolution and genetic constraints. We use a paternal half-sib-full-sib breeding design to investigate the genetic architecture of three life-history and morphological traits in the seed beetle, Callosobruchus maculatus. Heritability was significant for all traits under observation and genetic correlations between traits (r(A)) were low. Interestingly, we found substantial sex-specific genetic effects and low genetic correlations between sexes (r(MF)) in traits that are only moderately (weight at emergence) to slightly (longevity) sexually dimorphic. Furthermore, we found an increased sire ([Formula: see text]) compared to dam ([Formula: see text]) variance component within trait and sex. Our results highlight that the genetic architecture even of the same trait should not be assumed to be the same for males and females. Furthermore, it raises the issue of the presence of unnoticed environmental effects that may inflate estimates of heritability. Overall, our study stresses the fact that estimates of quantitative genetic parameters are not only population, time, environment, but also sex specific. Thus, extrapolation between sexes and studies should be treated with caution.     

Assumptions and Properties of Limiting Pathway Models for Analysis of Epistasis in Complex Traits

For most complex traits, results from genome-wide association studies show that the proportion of the phenotypic variance attributable to the additive effects of individual SNPs, that is, the heritability explained by the SNPs, is substantially less than the estimate of heritability obtained by standard methods using correlations between relatives. This difference has been called the “missing heritability”. One explanation is that heritability estimates from family (including twin) studies are biased upwards. Zuk et al. revisited overestimation of narrow sense heritability from twin studies as a result of confounding with non-additive genetic variance. They propose a limiting pathway (LP) model that generates significant epistatic variation and its simple parametrization provides a convenient way to explore implications of epistasis. They conclude that over-estimation of narrow sense heritability from family data (‘phantom heritability’) may explain an important proportion of missing heritability. We show that for highly heritable quantitative traits large phantom heritability estimates from twin studies are possible only if a large contribution of common environment is assumed. The LP model is underpinned by strong assumptions that are unlikely to hold, including that all contributing pathways have the same mean and variance and are uncorrelated. Here, we relax the assumptions that underlie the LP model to be more biologically plausible. Together with theoretical, empirical, and pragmatic arguments we conclude that in outbred populations the contribution of additive genetic variance is likely to be much more important than the contribution of non-additive variance.     

Genetic architecture of ovary size and asymmetry in European honeybee workers

The molecular basis of complex traits is increasingly understood but a remaining challenge is to identify their co-regulation and inter-dependence. Pollen hoarding (pln) in honeybees is a complex trait associated with a well-characterized suite of linked behavioral and physiological traits. In European honeybee stocks bidirectionally selected for pln, worker (sterile helper) ovary size is pleiotropically affected by quantitative trait loci that were initially identified for their effect on foraging behavior. To gain a better understanding of the genetic architecture of worker ovary size in this model system, we analyzed a series of crosses between the selected strains. The crossing results were heterogeneous and suggested non-additive effects. Three significant and three suggestive quantitative trait loci of relatively large effect sizes were found in two reciprocal backcrosses. These loci are not located in genome regions of known effects on foraging behavior but contain several interesting candidate genes that may specifically affect worker-ovary size. Thus, the genetic architecture of this life history syndrome may be comprised of pleiotropic, central regulators that influence several linked traits and other genetic factors that may be downstream and trait specific.     

The heritability of human longevity: A population-based study of 2872 Danish twin pairs born 1870–1900

The aim of this study was to explore, in a large and non-censored twin cohort, the nature (i.e., additive versus non-additive) and magnitude (i.e., heritability) of genetic influences on inter-individual differences in human longevity. The sample comprised all identified and traced non-emigrant like-sex twin pairs born in Denmark during the period 1870–1900 with a zygosity diagnosis and both members of the pairs surviving the age of 15 years. A total of 2872 pairs were included. Age at death was obtained from the Danish Central Person Register, the Danish Cause-of-Death Register and various other registers. The sample was almost non-censored on the date of the last follow-up (May 1, 1994), all but 0.6% had died, leaving a total of 2872 pairs for analysis. Proportions of variance attributable to genetic and environmental factors were assessed from variance-covariance matrices using the structural equation model approach. The most parsimonious explanation of the data was provided by a model that included genetic dominance (non-additive genetic effects caused by interaction within gene loci) and non-shared environmental factors (environmental factors that are individual-specific and not shared in a family). The heritability of longevity was estimated to be 0.26 for males and 0.23 for females. The small sex-difference was caused by a greater impact of non-shared environmental factors in the females. Heritability was found to be constant over the three 10-year birth cohorts included. Thus, longevity seems to be only moderately heritable. The nature of genetic influences on longevity is probably non-additive and environmental influences non-shared. There is no evidence for an impact of shared (family) environment.     

Quantitative genetic variation of leaf size and shape in a mixed diploid and triploid population of Populus

In the interspecific cross of Populus trichocarpa x P. deltoides, unexpected simultaneous occurrence of diploid hybrids and triploid hybrids (with two alleles from the female parent and one from the male parent at each locus) led us to examine the evolutionary genetic significance of this phenomenon. As expected, leaf size and shape of the triploid progeny are closer to the female P. trichocarpa than male P. deltoides parent. Although the pure triploid progeny population did not have higher genetic variance in leaf traits than the pure diploid population, the former appears to hide much non-additive genetic variance and display strong genetic control over the phenotypic plasticity of leaf traits. It is suggested that the cryptic non-additive variance, especially epistasis, can be released when a population is disturbed by changes in the environment. A mixed diploid and triploid progeny population combines phenotypic and genetic characteristics of both pure hybrids and is considered to be of adaptive significance for populars to survive and evolve in a fluctuating environment. The significant effect due to general and specific combining ability differences at the population level suggests that the population divergence of these two species is under additive and non-additive genetic control.     

Repeatability and heritability of reproductive traits in free-ranging snakes

The underlying genetic basis of life-history traits in free-ranging animals is critical to the effects of selection on such traits, but logistical constraints mean that such data are rarely available. Our long-term ecological studies on free-ranging oviparous snakes (keelbacks, Tropidonophis mairii (Gray, 1841), Colubridae) on an Australian floodplain provide the first such data for any tropical reptile. All size-corrected reproductive traits (egg mass, clutch size, clutch mass and post-partum maternal mass) were moderately repeatable between pairs of clutches produced by 69 female snakes after intervals of 49-1152 days, perhaps because maternal body condition was similar between clutches. Parent-offspring regression of reproductive traits of 59 pairs of mothers and daughters revealed high heritability for egg mass (h2= 0.73, SE=0.24), whereas heritability for the other three traits was low (< 0.37). The estimated heritability of egg mass may be inflated by maternal effects such as differential allocation of yolk steroids to different-sized eggs. High heritability of egg size may be maintained (rather than eroded by stabilizing selection) because selection acts on a trait (hatchling size) that is determined by the interaction between egg size and incubation substrate rather than by egg size alone. Variation in clutch size was mainly because of environmental factors (h2=0.04), indicating that one component of the trade-off between egg size and clutch size is under much tighter genetic control than the other. Thus, the phenotypic trade-off between egg size and egg number in keelback snakes occurs because each female snake must allocate a finite amount of energy into eggs of a genetically determined size.     

Sex differences and non-additivity in the effects of genes on personality.

New large-sample data show that non-additive genetic effects, probably epistatic interactions between loci, and sex-limited gene expression are significant features of the genetic architecture of human personality as measured by questionnaire scales of extraversion and neuroticism. Three large data sets--new data on large samples (n = 20,554) of US twins, their spouses, parents, siblings and children, correlations for Australian twins (n = 7,532), and previously published twin data from Finland (n = 14,288)--are subjected to an integrated analysis to test alternative hypotheses about the genetic causes of family resemblance in personality. When allowance is made for differences in reliability of the scales, the combined data are consistent with the same model for variation. There are significant amounts of genetic non-additivity for both dimensions of personality. The evidence favours additive x additive epistatic interactions rather than dominance. In the case of neuroticism, there is especially strong evidence of sex differences in genetic architecture favouring a greater relative contribution of non-additive genetic effects in males. The data confirm previous claims to find no major contribution of the shared environment of twins and siblings to these dimensions of personality. Correlations between spouses are zero, and the correlations for very large samples of siblings and non-identical twins do not differ significantly.     

Lifetime Reproductive Success and Heritability in Nature

The observation that traits closely related to fitness ("fitness traits") have lower heritabilities than traits more distantly associated with fitness has traditionally been framed in terms of Fisher's fundamental theorem of natural selection-fitness traits are expected to have low levels of additive genetic variance due to rapid fixation of alleles conferring highest fitness. Subsequent treatments have challenged this view by pointing out that high environmental and nonadditive genetic contributions to phenotypic variation may also explain the low heritability of fitness traits. Analysis of a large data set from the collared flycatcher Ficedula albicollis confirmed a previous finding that traits closely associated with fitness tend to have lower heritability. However, analysis of coefficients of additive genetic variation (CVA) revealed that traits closely associated with fitness had higher levels of additive genetic variation (VA) than traits more distantly associated with fitness. Hence, the negative relationship between a trait's association with fitness and its heritability was not due to lower levels of VA in fitness traits but was due to their higher residual variance. However, whether the high residual variance was mainly due to higher levels of environmental variance or due to higher levels of nonadditive genetic variance remains a challenge to be addressed by further studies. Our results are consistent with earlier suggestions that fitness-related traits may have more complex genetic architecture than traits more distantly associated with fitness.     

The Relationship of DNA Methylation with Age,Gender and Genotype in Twins and Healthy Controls

Cytosine-5 methylation within CpG dinucleotides is a potentially important mechanism of epigenetic influence on human traits and disease. In addition to influences of age and gender, genetic control of DNA methylation levels has recently been described. We used whole blood genomic DNA in a twin set (23 MZ twin-pairs and 23 DZ twin-pairs, N = 92) as well as healthy controls (N = 96) to investigate heritability and relationship with age and gender of selected DNA methylation profiles using readily commercially available GoldenGate bead array technology. Despite the inability to detect meaningful methylation differences in the majority of CpG loci due to tissue type and locus selection issues, we found replicable significant associations of DNA methylation with age and gender. We identified associations of genetically heritable single nucleotide polymorphisms with large differences in DNA methylation levels near the polymorphism (cis effects) as well as associations with much smaller differences in DNA methylation levels elsewhere in the human genome (trans effects). Our results demonstrate the feasibility of array-based approaches in studies of DNA methylation and highlight the vast differences between individual loci. The identification of CpG loci of which DNA methylation levels are under genetic control or are related to age or gender will facilitate further studies into the role of DNA methylation and disease.     

Genetic parameter estimates in a clonally replicated progeny test of teak (Tectona grandis Linn. f.)

Teak (Tectona grandis Linn. f.) has been planted extensively in the tropics for its highly valued timber. We analysed data from a 3.5-year-old teak progeny test with clonal replication located in northern Australia. Additive and non-additive genetic variances were estimated for commercially important traits. Trees originating from seedlings were on average 2% taller and 4% straighter than those of the same genotype originating from cuttings. Non-additive genetic variance represented 35–50% of total genetic variance for growth traits and 63% of total genetic variance for incidence of flowering. Narrow-sense heritability was 0.22 for diameter, 0.18 for height and volume, 0.07 for stem straightness, 0.05 for insect defoliation, 0.03 for epicormic sprouts and 0.30 for incidence of flowering (estimated on an assumed underlying continuous scale). Broad-sense heritability was 0.37 for diameter, 0.28 for height, 0.35 for volume, 0.12 for stem straightness, 0.06 for insect defoliation, 0.12 for epicormic sprouts and 0.71 for incidence of flowering. Positive correlations were found between tree volume and flowering and between tree volume and stem straightness. The presence of sizeable non-additive variance supports the selection and deployment of clones to capture the full extent of genetic variation in commercially important traits.     

Tail colour signals performance in blue tit nestlings

Indirect sexual selection arises when reproductive individuals choose their mates based on heritable ornaments that are genetically correlated to fitness. Evidence for genetic associations between ornamental colouration and fitness remains scarce. In this study, we investigate the quantitative genetic relationship between different aspects of tail structural colouration (brightness, hue and UV chroma) and performance (cell‐mediated immunity, body mass and wing length) in blue tit (Cyanistes caeruleus) nestlings. In line with previous studies, we find low heritability for structural colouration and moderate heritability for performance measures. Multivariate animal models show positive genetic correlations between the three measures of performance, indicating quantitative genetic variation for overall performance, and tail brightness and UV chroma, two genetically independent colour measures, are genetically correlated with performance (positively and negatively, respectively). Our results suggest that mate choice based on independent aspects of tail colouration can have fitness payoffs in blue tits and provide support for the indirect benefits hypothesis. However, low heritability of tail structural colouration implies that indirect sexual selection on mate choice for this ornament will be a weak evolutionary force.     

Does artificial selection on diapause incidence cause correlational changes in other life‐history traits? Case study in a spider mite population

Theoretical studies suggest that the timing of entering hibernation by arthropods has large effects on long-term fitness, incurring strong selection pressure on diapause attributes every year. On the other hand, diapause attributes are often genetically correlated with other important life-history traits such as fecundity or development time. To understand the evolutionary process of life cycle formation, there is a need to investigate not only diapause attributes themselves but also their genetic association with other life-history traits. The Kanzawa spider mite, Tetranychus kanzawai Kishida (Acari: Tetranychidae), is a small herbivore that lives on the undersurface of host plant leaves. This mite has been investigated for the mode of inheritance of diapause attributes, but scarcely for genetic correlations with other life-history traits. Here, I investigated whether diapause proneness, measured as the proportion of diapausing females under short-day conditions, is genetically correlated with fecundity or development time under long-day conditions using artificial selection experiments. Diapause incidence responded to the selection for both increasing and decreasing directions, suggesting that high genetic variance in diapause proneness is maintained in the study population. However, the change in proportion of diapausing females during the selection period was not associated with responses in fecundity or development time. These results suggest that diapause proneness and other life-history traits have different genetic backgrounds, and thus diapause proneness may freely evolve without being constrained by changes in other life-history traits.