Epigenetics,heritability and longitudinal analysis

Background

Longitudinal data and repeated measurements in epigenome-wide association studies (EWAS) provide a rich resource for understanding epigenetics. We summarize 7 analytical approaches to the GAW20 data sets that addressed challenges and potential applications of phenotypic and epigenetic data. All contributions used the GAW20 real data set and employed either linear mixed effect (LME) models or marginal models through generalized estimating equations (GEE). These contributions were subdivided into 3 categories: (a) quality control (QC) methods for DNA methylation data; (b) heritability estimates pretreatment and posttreatment with fenofibrate; and (c) impact of drug response pretreatment and posttreatment with fenofibrate on DNA methylation and blood lipids.

Results

Two contributions addressed QC and identified large statistical differences with pretreatment and posttreatment DNA methylation, possibly a result of batch effects. Two contributions compared epigenome-wide heritability estimates pretreatment and posttreatment, with one employing a Bayesian LME and the other using a variance-component LME. Density curves comparing these studies indicated these heritability estimates were similar. Another contribution used a variance-component LME to depict the proportion of heritability resulting from a genetic and shared environment. By including environmental exposures as random effects, the authors found heritability estimates became more stable but not significantly different. Two contributions investigated treatment response. One estimated drug-associated methylation effects on triglyceride levels as the response, and identified 11 significant cytosine-phosphate-guanine (CpG) sites with or without adjusting for high-density lipoprotein. The second contribution performed weighted gene coexpression network analysis and identified 6 significant modules of at least 30 CpG sites, including 3 modules with topological differences pretreatment and posttreatment.

Conclusions

Four conclusions from this GAW20 working group are: (a) QC measures are an important consideration for EWAS studies that are investigating multiple time points or repeated measurements; (b) application of heritability estimates between time points for individual CpG sites is a useful QC measure for DNA methylation studies; (c) drug intervention demonstrated strong epigenome-wide DNA methylation patterns across the 2 time points; and (d) new statistical methods are required to account for the environmental contributions of DNA methylation across time. These contributions demonstrate numerous opportunities exist for the analysis of longitudinal data in future epigenetic studies.

     

An exact form of the breeder's equation for the evolution of a quantitative trait under natural selection

Starting with the Price equation, I show that the total evolutionary change in mean phenotype that occurs in the presence of fitness variation can be partitioned exactly into five components representing logically distinct processes. One component is the linear response to selection, as represented by the breeder's equation of quantitative genetics, but with heritability defined as the linear regression coefficient of mean offspring phenotype on parent phenotype. The other components are identified as constitutive transmission bias, two types of induced transmission bias, and a spurious response to selection caused by a covariance between parental fitness and offspring phenotype that cannot be predicted from parental phenotypes. The partitioning can be accomplished in two ways, one with heritability measured before (in the absence of) selection, and the other with heritability measured after (in the presence of) selection. Measuring heritability after selection, though unconventional, yields a representation for the linear response to selection that is most consistent with Darwinian evolution by natural selection because the response to selection is determined by the reproductive features of the selected group, not of the parent population as a whole. The analysis of an explicitly Mendelian model shows that the relative contributions of the five terms to the total evolutionary change depends on the level of organization (gene, individual, or mated pair) at which the parent population is divided into phenotypes, with each frame of reference providing unique insight. It is shown that all five components of phenotypic evolution will generally have nonzero values as a result of various combinations of the normal features of Mendelian populations, including biparental sex, allelic dominance, inbreeding, epistasis, linkage disequilibrium, and environmental covariances between traits. Additive genetic variance can be a poor predictor of the adaptive response to selection in these models. The narrow-sense heritability sigma2A/sigma2P should be viewed as an approximation to the offspring-parent linear regression rather than the other way around.     

The estimation of genetic relationships using molecular markers and their efficiency in estimating heritability in natural populations

Molecular marker data collected from natural populations allows information on genetic relationships to be established without referencing an exact pedigree. Numerous methods have been developed to exploit the marker data. These fall into two main categories: method of moment estimators and likelihood estimators. Method of moment estimators are essentially unbiased, but utilise weighting schemes that are only optimal if the analysed pair is unrelated. Thus, they differ in their efficiency at estimating parameters for different relationship categories. Likelihood estimators show smaller mean squared errors but are much more biased. Both types of estimator have been used in variance component analysis to estimate heritability. All marker-based heritability estimators require that adequate levels of the true relationship be present in the population of interest and that adequate amounts of informative marker data are available. I review the different approaches to relationship estimation, with particular attention to optimizing the use of this relationship information in subsequent variance component estimation.     

Epigenetics,heritability and longitudinal analysis

Background

Longitudinal data and repeated measurements in epigenome-wide association studies (EWAS) provide a rich resource for understanding epigenetics. We summarize 7 analytical approaches to the GAW20 data sets that addressed challenges and potential applications of phenotypic and epigenetic data. All contributions used the GAW20 real data set and employed either linear mixed effect (LME) models or marginal models through generalized estimating equations (GEE). These contributions were subdivided into 3 categories: (a) quality control (QC) methods for DNA methylation data; (b) heritability estimates pretreatment and posttreatment with fenofibrate; and (c) impact of drug response pretreatment and posttreatment with fenofibrate on DNA methylation and blood lipids.

Results

Two contributions addressed QC and identified large statistical differences with pretreatment and posttreatment DNA methylation, possibly a result of batch effects. Two contributions compared epigenome-wide heritability estimates pretreatment and posttreatment, with one employing a Bayesian LME and the other using a variance-component LME. Density curves comparing these studies indicated these heritability estimates were similar. Another contribution used a variance-component LME to depict the proportion of heritability resulting from a genetic and shared environment. By including environmental exposures as random effects, the authors found heritability estimates became more stable but not significantly different. Two contributions investigated treatment response. One estimated drug-associated methylation effects on triglyceride levels as the response, and identified 11 significant cytosine-phosphate-guanine (CpG) sites with or without adjusting for high-density lipoprotein. The second contribution performed weighted gene coexpression network analysis and identified 6 significant modules of at least 30 CpG sites, including 3 modules with topological differences pretreatment and posttreatment.

Conclusions

Four conclusions from this GAW20 working group are: (a) QC measures are an important consideration for EWAS studies that are investigating multiple time points or repeated measurements; (b) application of heritability estimates between time points for individual CpG sites is a useful QC measure for DNA methylation studies; (c) drug intervention demonstrated strong epigenome-wide DNA methylation patterns across the 2 time points; and (d) new statistical methods are required to account for the environmental contributions of DNA methylation across time. These contributions demonstrate numerous opportunities exist for the analysis of longitudinal data in future epigenetic studies.

     

离散性状QTL区间定位的最大似然方法

采用最大似然区间定位法对阈模型与一般线性模型的QTL定位效率进行了比较,并对影响离散性状QTL检测效率的主要因素（QTL效应、性状的遗传力和表型发生率）进行了模拟研究,实验设计为多个家系的女儿设计．资源群体大小为500头。研究结果表明：在QTL参数估计及检验功效方面,阈模型方法具有较大的优势,对离散性状QTL定位的效率明显高于LM（Linear Model）方法,定位的准确性也较高。另外,性状遗传力、QTL效应的大小和性状表型发生率对QTL定位的准确度也有直接的影响,随着性状遗传力和表型发生率的提高,随着QTL效应的增大,QTL定位的效率也进一步提高。     

Prediction error variance and expected response to selection,when selection is based on the best predictor – for Gaussian and threshold characters,traits following a Poisson mixed model and survival traits

In this paper, we consider selection based on the best predictor of animal additive genetic values in Gaussian linear mixed models, threshold models, Poisson mixed models, and log normal frailty models for survival data (including models with time-dependent covariates with associated fixed or random effects). In the different models, expressions are given (when these can be found – otherwise unbiased estimates are given) for prediction error variance, accuracy of selection and expected response to selection on the additive genetic scale and on the observed scale. The expressions given for non Gaussian traits are generalisations of the well-known formulas for Gaussian traits – and reflect, for Poisson mixed models and frailty models for survival data, the hierarchal structure of the models. In general the ratio of the additive genetic variance to the total variance in the Gaussian part of the model (heritability on the normally distributed level of the model) or a generalised version of heritability plays a central role in these formulas.     

性状遗传力与QTL方差对标记辅助选择效果的影响

在采用动物模型标记辅助最佳线性无偏预测方法对个体育种值进行估计的基础上,模拟了在一个闭锁群体内连续对单个性状选择10个世代的情形,并系统地比较了性状遗传力和QTL方差对标记辅助选择所获得的遗传进展、QTL增效基因频率和群体近交系数变化的影响。结果表明：在对高遗传力和QTL方差较小的性状实施标记辅助选择时,可望获得更大的遗传进展;遗传力越高,QTL方差越大,则QTL增效基因频率的上升速度越快;遗传力较高时,群体近交系数上升的速度较为缓慢,而QTL方差对群体近交系数上升速度的影响则不甚明显。结合前人关于标记辅助选择相对效率的研究结果,可以认为：当选择性状的遗传力和QTL方差为中等水平时,标记辅助选择可望获得理想的效果。     

Sex differences in sand lizard telomere inheritance: paternal epigenetic effects increases telomere heritability and offspring survival

Background

To date, the only estimate of the heritability of telomere length in wild populations comes from humans. Thus, there is a need for analysis of natural populations with respect to how telomeres evolve.

Methodology/Principal Findings

Here, we show that telomere length is heritable in free-ranging sand lizards, Lacerta agilis. More importantly, heritability estimates analysed within, and contrasted between, the sexes are markedly different; son-sire heritability is much higher relative to daughter-dam heritability. We assess the effect of paternal age on Telomere Length (TL) and show that in this species, paternal age at conception is the best predictor of TL in sons. Neither paternal age per se at blood sampling for telomere screening, nor corresponding age in sons impact TL in sons. Processes maintaining telomere length are also associated with negative fitness effects, most notably by increasing the risk of cancer and show variation across different categories of individuals (e.g. males vs. females). We therefore tested whether TL influences offspring survival in their first year of life. Indeed such effects were present and independent of sex-biased offspring mortality and offspring malformations.

Conclusions/Significance

TL show differences in sex-specific heritability with implications for differences between the sexes with respect to ongoing telomere selection. Paternal age influences the length of telomeres in sons and longer telomeres enhance offspring survival.     

Heritability estimates for carcass traits of cattle: a review

We present estimates of heritability for carcass traits of cattle published in the scientific literature. Seventy-two papers published from 1962 to 2004, which reported estimates of heritability for carcass traits, were reviewed. The unweighted means of estimates of heritability for 14 carcass traits by slaughter end point (age, weight, and fat depth) were calculated. Among the three end points, carcass weight, backfat thickness, longissimus muscle area, and marbling score were the carcass traits with the most estimates of heritability (56 相似文献   

The quantitative genetics of sexually selected traits,preferred traits and preference: a review and analysis of the data

The maintenance of variation in sexually selected traits is a puzzle that has received increasing attention in the past several decades. Traits that are related to fitness, such as life‐history or sexually selected traits, are expected to have low additive genetic variance (and hence, heritability) due to the rapid fixation of advantageous alleles. However, previous analyses have suggested that the heritabilities of sexually selected traits are on average higher than nonsexually selected traits. We show that the heritabilities of sexually selected traits are not significantly different from those of nonsexually selected traits overall or when separated into the three trait categories: behavioural, morphological and physiological. In contrast with previous findings, the heritability of preference is quite low (h² = 0.25 ± 0.06) and is in the same range as life‐history traits. We distinguish preferred traits as a category of sexually selected traits and find that the heritability of the former is not significantly different than sexually selected traits overall (0.48 ± 0.04 vs. 0.46 ± 0.03). We test the hypothesis that the heritability of sexually selected traits is negatively correlated with the strength of sexual selection. As predicted, there is a significant negative correlation between the heritabilities of sexually selected traits and the strength of selection. This suggests that heritabilities do indeed decrease as sexual selection increases but sexual selection is not strong enough to cause heritabilities of sexually selected traits to deviate from the same type of nonsexually selected traits.     

Optimal designs for linkage disequilibrium mapping and candidate gene association tests in livestock populations

Simulation was used to evaluate the performance of different selective genotyping strategies when using linkage disequilibrium across large half-sib families to position a QTL within a previously defined genomic region. Strategies examined included standard selective genotyping and different approaches of discordant and concordant sib selection applied to arbitrary or selected families. Strategies were compared as a function of effect and frequency of QTL alleles, heritability, and phenotypic expression of the trait. Large half-sib families were simulated for 100 generations and 2% of the population was genotyped in the final generation. Simple ANOVA was applied and the marker with the greatest F-value was considered the most likely QTL position. For traits with continuous phenotypes, genotyping the most divergent pairs of half-sibs from all families was the best strategy in general, but standard selective genotyping was somewhat more precise when heritability was low. When the phenotype was distributed in ordered categories, discordant sib selection was the optimal approach for positioning QTL for traits with high heritability and concordant sib selection was the best approach when genetic effects were small. Genotyping of a few selected sibs from many families was generally more efficient than genotyping many individuals from a few highly selected sires.     

Heritability in the genomics era--concepts and misconceptions

Heritability allows a comparison of the relative importance of genes and environment to the variation of traits within and across populations. The concept of heritability and its definition as an estimable, dimensionless population parameter was introduced by Sewall Wright and Ronald Fisher nearly a century ago. Despite continuous misunderstandings and controversies over its use and application, heritability remains key to the response to selection in evolutionary biology and agriculture, and to the prediction of disease risk in medicine. Recent reports of substantial heritability for gene expression and new estimation methods using marker data highlight the relevance of heritability in the genomics era.     

Heritability of cortisol response to confinement stress in European sea bass dicentrarchus labrax

Background

In fish, the most studied production traits in terms of heritability are body weight or growth, stress or disease resistance, while heritability of cortisol levels, widely used as a measure of response to stress, is less studied. In this study, we have estimated heritabilities of two growth traits (body weight and length) and of cortisol response to confinement stress in the European sea bass.

Findings

The F1 progeny analysed (n = 922) belonged to a small effective breeding population with contributions from an unbalanced family structure of just 10 males and 2 females. Heritability values ranged from 0.54 (±0.21) for body weight to 0.65 (±0.22) for standard body length and were low for cortisol response i.e. 0.08 (±0.06). Genetic correlations were positive (0.94) between standard body length and body weight and negative between cortisol and body weight and between cortisol and standard body length (−0.60 and −0.55, respectively).

Conclusion

This study confirms that in European sea bass, heritability of growth-related traits is high and that selection on such traits has potential. However, heritability of cortisol response to stress is low in European sea bass and since it is known to vary greatly among species, further studies are necessary to understand the reasons for these differences.     

Heritability estimates in contrasting environments as influenced by the adaptation level of barley germ plasm

Low heritability estimates in marginal or stress environments have often been used as one of the main justification for conducting selection work in environments with optimum or near‐optimum conditions for plant growth and grain yield. In this study, we have examined the relationships between grain yield and broad‐sense heritability in four groups of recombinant inbred lines (RILs) obtained from four barley crosses derived from parents differing in adaptation to stress. The RILs and the parents were grown in 13 combinations of years and locations (environments) in Syria and Jordan. Grain yields ranged from about 30 kg ha^?1 to nearly 4000 kg ha^?1 and genotype × environment interactions explained about half of the total variance for environmentally standardised data. Broad‐sense heritability in the individual year–location combinations varied from 0 to 0.68 and both the simple correlation and the rank correlation coefficients between grain yield and heritability were not significant. Genotype × years within individual locations, which measures the repeatability of a location in discriminating between genotypes, was also independent from the yield level, confirming that low‐yielding locations can be reliable selection environments. Also, there was no relationship between the type of cross and the magnitude of heritability in the various environments, but, as expected, the magnitude of heritability was significantly associated with the genetic distance between the parents. It is concluded that, holding all other factors affecting response to selection constant, concerns about the magnitude of heritability at low‐yielding locations are not justified and should not prevent them from being used as selection sites.     

Scaling threshold characters

A simple method of scaling ordered categorical responses having a joint distribution with an underlying normal variable is presented. Scores are developed that maximize heritability of the observed variate and that in the class of scores based on polychotomies: (1) maximize the correlation between score and the underlying genetic value to be predicted, and (2) minimize mean-square prediction error. Several examples suggest little is lost, in terms of heritability, by using equally spaced scores. The proposed scaling method discriminates among candidates for selection that would be tied if equally spaced scores are used and sometimes yields different rankings of candidates.     

Parental effort and response to nestling begging in the house sparrow: repeatability,heritability and parent–offspring co‐evolution

Parental effort has a direct impact on individual fitness. Theoretical models exploring how parental effort evolves to cope with offspring demand and sexual conflicts may differ in the assumptions they make in respect to the genetic heritability of parental behaviours. Only a few attempts, however, have been made to estimate the heritability of parental behaviours and their possible co‐evolution with offspring solicitation behaviour. Analysing parent and offspring behaviours in four generations of cross‐fostered broods of house sparrows, we found that parental effort (food delivery rate) was repeatable across consecutive broods and heritable across generations. In contrast, parental response to experimentally induced changes in nestling begging was neither repeatable across broods nor heritable across generations or correlated to nestling begging. Thus, the results give no indication for genetic covariance between begging intensity and parental response, but provide the first cross‐fostering‐based evidence for the heritability of parental investment levels across generations.     

Genetic analysis of blood pressure in the Milan hypertensive strain of rat (Rattus norvegicus)

Estimates of heritability (h2) of blood pressure level and the number of loci controlling the trait were derived from two genetic crosses involving the Milan hypertensive strain of rat and its control with normal blood pressure. In the genetic cross involving backcrosses, the estimates were h2 = 64% and the number of loci was two or three; there was some evidence of dominance of the alleles for normal blood pressures. In the other cross with only F2's, the degree of genetic determination (heritability in the broad sense) was 45%, involving at least three loci.     

Relationship between pyrimidine dimers, 6-4 photoproducts, repair synthesis and cell survival: studies using cells from patients with trichothiodystrophy

Trichothiodystrophy is a genetic disease which in the majority of cases studied is associated with a deficiency in the ability to repair UV damage in cellular DNA. Three categories of UV response have been identified. In type 1 the response is completely normal, whereas type 2 cells are deficient in excision-repair, with properties indistinguishable from those of XP complementation group D. Type 3 cells have normal survival following UV-irradiation and normal rates of removal of cyclobutane pyrimidine dimer sites. Nevertheless repair synthesis is reduced by 50% in these cell strains and this is associated with a marked reduction in the repair of 6-4 photoproducts from cellular DNA. The present results show that 50% or more of repair synthesis at early times after irradiation of normal primary human fibroblasts is attributable to repair of 6-4 products. They also suggest that repair of cyclobutane dimers is crucial for cell survival.     

HERITABILITY OF NESTLING BEGGING INTENSITY IN THE HOUSE SPARROW (PASSER DOMESTICUS)

Evolutionary theory of parent–offspring conflict assumes that offspring food solicitation behavior, known as begging, and parental response to begging are subjected to selection and coevolution. This assumption implies that begging intensity should be heritable, at least to some degree. Although some studies have suggested that begging is heritable, the evidence for this is rare and mostly indirect. To assess the heritability of begging we used artificial selection, sibling analysis, and the monitoring of begging intensity in four generations of cross-fostered captive house sparrow nestlings. We also contrasted the heritability of begging with that of morphological traits, known to be heritable in this species. Our results show that adult wing length and body mass were heritable as expected. The heritability estimates of the visual and vocal components of nestling begging (standardized for food deprivation and body mass) were low to moderate, as expected for behavioral traits in general, and lower than previously reported for passerine birds. Our sibling analysis shows that common environment had much greater effect on begging than genetic origin, suggesting that begging evolution may be strongly influenced by gene–environment interaction, probably through the mechanisms that adjust begging response to environmental and social conditions.     

Heritability of brain size and surface features in rhesus macaques (Macaca mulatta)

The extent of heritability for overall brain size and regional cortical surface features such as sulcus lengths is important for demonstrating a genetic component to the observed phenotypic differences among individuals and for evaluating the potential for evolutionary change in response to selection. Although the genetics of brain size has been extensively considered, the detailed morphology of the cortical surface has not previously been subjected to genetic analysis. We estimated the heritability of brain size and cortical sulcus lengths using 438 endocranial casts taken from skeletons of rhesus macaques (Macaca mulatta) from the Cayo Santiago population. Estimates were obtained both by mother-offspring regression and symmetric-differences-squared (SDS) methods. Brain size, measured as cranial capacity, was highly and significantly heritable in this population, confirming results of previous studies with laboratory mice. Overall, cortical sulcus lengths were also heritable, with 35% of the sulci significantly heritable at the 5% level in the mother-offspring analysis. The average mother-offspring heritability estimate, 0.31, was the same as the average heritability obtained previously from a series of 56 cranial metric characters. The SDS analyses generally corresponded to the findings based on mother-offspring regressions, although the significance test appeared more conservative. Both gross and detailed morphology of the brain are heritable.