PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45191.
Published online 2017 March 24. doi:  10.1038/srep45191
PMCID: PMC5364555

Genome-wide random regression analysis for parent-of-origin effects of body composition allometries in mouse

Abstract

Genomic imprinting underlying growth and development traits has been recognized, with a focus on the form of absolute or pure growth. However, little is known about the effect of genomic imprinting on relative growth. In this study, we proposed a random regression model to estimate genome-wide imprinting effects on the relative growth of multiple tissues and organs to body weight in mice. Joint static allometry scaling equation as sub-model is nested within the genetic effects of markers and polygenic effects caused by a pedigree. Both chromosome-wide and genome-wide statistical tests were conducted to identify imprinted quantitative trait nucleotides (QTNs) associated with relative growth of individual tissues and organs to body weight. Real data analysis showed that three of six analysed tissues and organs are significantly associated with body weight in terms of phenotypic relative growth. At the chromosome-wide level, a total 122 QTNs were associated with allometries of kidney, spleen and liver weights to body weight, 36 of which were imprinted with different imprinting fashions. Further, only two imprinted QTNs responsible for relative growth of spleen and liver were verified by genome-wide test. Our approach provides a general framework for statistical inference of genomic imprinting underlying allometry scaling in animals.

Genomic imprinting, an epigenetic phenomenon of parent-of-origin-specific gene expression, has been widely observed in plants1 and animals2,3,4,5,6,7,8,9 and has been recognized for its role in shaping developmental processes10,11,12. Genomic imprinting is a highly complex process that is involved in a number of growth axes operating coordinately at different development stages and showing a time-dependent effect during development. Most imprinted genes play important roles in controlling embryonic and post-natal growth and development in mammals. Based on the parental origin of the expressed allele, imprinting is classified into paternal imprinting and maternal imprinting. Imprinting can be further categorized as complete imprinting, when only one allele is expressed, or partial imprinting, when both alleles are expressed but at different levels13,14,15. To date, imprinting quantitative trait loci (iQTL) for growth and development traits have been identified in which the traits are measured in absolute growth. Little is known about how genomic imprinting affects relative growth, partially due to a lack of efficient statistical modelling and inference procedures.

As a measure of relative growth, allometry scaling describes the relationship between the entire body size and partial body size or between certain two biological traits. A “simple equation of allometry” has been initially developed to quantify allometry scaling16. In a simple allometry equation, the two variables include not only body size measured in scale of length and weight, but also body shape, density and volume. By taking into account the correlations among multiple body parts, the joint static allometry scaling model17 was proposed to simultaneously evaluate allometry scalings of multiple body parts to the entire body. Allometry scaling relationships between different biological traits contain three terms of allometries: static allometry, ontogenetic allometry and evolutionary allometry18,19,20. Static allometry refers to the relative growth between two different traits in adult or at a particular developmental stage. Ontogenetic allometry is the growth trajectory of one trait relative to the other in ontogeny. Evolutionary allometry is the relative growth between traits across species. The differentiation in allometries among traits has been thought to be a driving force by which morphology and structure evolve21.

Methods for detecting imprinting loci have been adapted from methods of interval mapping for Mendelian quantitative trait loci (QTL). Imprinting effects can be estimated using either least squares9,22,23,24 or maximum likelihood methods25. Multi-step tests for contrast models have been also proposed to identify the imprinting pattern9,22,23,24. With Bayesian model selection, Yang, et al.26 estimated genomic imprinting effects and inferred more genomic imprinting patterns than those summarized by Cheverud, et al.27. In the aspect of relative growth, genetic analyses for allometry scalings between biological traits have been carried out by embedding a simple allometry equation into additive genetic effects of the mixed linear model28,29 and into genotypic effects of genetic model for mapping QTL30,31,32. However, all these approaches to map allometry scalings are based on a single QTL model. When allometry scalings are controlled by multiple QTLs, they perform low power to detect QTLs. So far, no statistical method has been proposed to search for imprinted QTL for multiple allometry scalings.

Most current reports on genomic imprinting have been focused on brain and placenta in mammals33,34,35. Imprinting also affects adult traits either through the persistent effects of early growth and development36,37,38 or through direct effects on adult physiology. In absolute growth, a genome-wide mapping has been undertaken to identify quantitative trait nucleotides (QTNs) that have an imprinting effect on adult body composition using a three-generation intercross between inbred mouse strains39. To identify imprinted QTN for relative growth of tissues and organs to body weight in mice, here we developed a random regression model in which a joint static allometry scaling model is nested into genetic effects of markers and polygenic effects. The model was derived under the random regression framework, in which genetic effects of markers and genomic imprinting patterns are statistically inferred using chromosome-wide and genome-wide statistical tests.

Methods

Joint allometric scaling model

Let An external file that holds a picture, illustration, etc.
Object name is srep45191-m1.jpg denote m partial body sizes and y be the entire body size. The joint static allometry model17 is defined as

An external file that holds a picture, illustration, etc.
Object name is srep45191-m2.jpg

where β0 is an intercept, and An external file that holds a picture, illustration, etc.
Object name is srep45191-m3.jpg are partial scaling exponents of the jth partial component to the entire body size. Nonlinear least squares method is generally used to get unbiased estimators of the partial scaling exponents. Under the linear mixed model framework, a linear transformation of model (1) facilitates genetic analysis of allometric scalings. Taking the natural logarithm on both sides of model (1), we have

An external file that holds a picture, illustration, etc.
Object name is srep45191-m4.jpg

With the transformation, the joint static allometric scaling model can be optimized through a stepwise regression analysis.

Random regression model for multiple static allometries

For a family-based population, suppose that n individuals are genotyped for q markers and observed for m partial body sizes. Four possible genotypes, denoted as QQ, Qq, qQ and qq, are distinguished at each marker. The first allele of each genotype is inherited from the paternal parent when considering the parental origin of alleles. In addition to the additive and dominance effects, the two reciprocal heterozygotes Qq and qQ carry the genetic imprinting effect that reflects the difference in allele expression derived from the two parents. Polygenic effects can also be estimated using data from multiple full-sibling and half-sibling families. Following a simple animal model for a single trait, the relationship between the logarithm of entire body size and markers can be modelled as

An external file that holds a picture, illustration, etc.
Object name is srep45191-m5.jpg

where yi is the entire body size for the ith individual; bi represents the lth fixed effect in l systematic environments, such as sexes; aj, dj and ij are additive, dominance and imprinting genetic effects for the jth marker; hil, zij, wij and sij are indicator variables corresponding to bl, aj, dj and ij, where zij, wij and sij are defined in Mantey et al.40; gi is the polygenic genetic effect derived from the pedigree information, assuming that An external file that holds a picture, illustration, etc.
Object name is srep45191-m6.jpg with A being a relationship matrix and An external file that holds a picture, illustration, etc.
Object name is srep45191-m7.jpg being the polygenic genetic variance; and ei is the residual error with An external file that holds a picture, illustration, etc.
Object name is srep45191-m8.jpg.

To genetically analyse allometric scalings of multiple partial to entire body sizes, we embed model (2) into systematic environments, genetic effects of markers and polygenic effects in model (3), yielding the following random regression animal model41,42:

An external file that holds a picture, illustration, etc.
Object name is srep45191-m9.jpg

where An external file that holds a picture, illustration, etc.
Object name is srep45191-m10.jpg with θ [set membership] b, a, d, i or g. Let An external file that holds a picture, illustration, etc.
Object name is srep45191-m11.jpg, with [var phi] [set membership] h, z, w or s and [var phi] [set membership] h, z, w or s, An external file that holds a picture, illustration, etc.
Object name is srep45191-m12.jpg, An external file that holds a picture, illustration, etc.
Object name is srep45191-m13.jpg with θ [set membership] b, a, d or i and An external file that holds a picture, illustration, etc.
Object name is srep45191-m14.jpg, then in matrix form, model (3) can be rewritten as,

An external file that holds a picture, illustration, etc.
Object name is srep45191-m15.jpg

In this model, genetic effects are fixed for all markers and polygenic effects are random for all individuals. Thus,An external file that holds a picture, illustration, etc.
Object name is srep45191-m16.jpg and An external file that holds a picture, illustration, etc.
Object name is srep45191-m17.jpg, where A and G are the numerator-relationship matrix and polygenic genetic covariance matrix for multiple allometric scalings, respectively.

Statistical inference for genomic imprinting

Restricted maximum likelihood method43 is implemented to estimate parameters in the random regression model described above, which gives maximum likelihood estimates and standard errors for fixed genetic effects at each marker, in addition to polygenic and residual variances. A student t statistic is formulated to statistically infer significance of marker genetic effect, which is calculated as

An external file that holds a picture, illustration, etc.
Object name is srep45191-m18.jpg

Similarly, another t statistic is calculated as

An external file that holds a picture, illustration, etc.
Object name is srep45191-m19.jpg

to statistically infer imprinting patterns.

The two test statistics asymptotically follow a standard normal distribution if the sample size is large enough so that the difference between the sample size and the number of estimated parameters exceeds 120. Under a standard normal distribution, the critical value is taken as 1.96 or 1.301 for −log(p) (p is the probability of the test statistic greater than 1.96) at the significance level of 5%. Thus, the iQTLs can be identified at chromosome-wide and genome-wide levels: by chromosome-wide test, the markers with significant genetic effects are screened from each chromosome, and then by genome-wide test, iQTLs are detected from all the significant markers screened by chromosome-wide tests.

In fact, we need to test whether or not a QTL exists, in which fashion the detected QTL inherits and what imprinting pattern the iQTL carries. A marker is identified as a QTL if one of the genetic effects is significantly different from zero by formula (6). In the case of no significant imprinting effect, the detected QTL is defined as a Mendelian QTL; otherwise, it is defined as an iQTL. The imprinting pattern can be further classified as either additive imprinting or dominance imprinting. Classification of imprinting patterns depends on the value and sign of i relative to a and d. The additive imprinting is composed of four subtypes: the complete or partial paternal additive imprinting corresponding to hypothesis d = 0 and a = i or d = 0 and a  i, respectively; and the complete or partial maternal additive imprinting corresponding to hypothesis d = 0 and a = i or d = 0 and a  i, respectively. The dominance imprinting is further classified into bipolar dominance (H0: a = 0 and d = 0), polar over-dominance (H0: a = 0 and d = i) and polar under-dominance (H0: a = 0 and d = i). Following the definition of different imprinting types and the corresponding null hypothesis26,27,44, the imprinting pattern for the detected QTL can be statistically inferred by formula (7).

Case analysis

A set of reciprocal cross families for distinguishing the four genotypes considering parental origin was derived from an F2 intercross of large (LG/J) and small (SM/J) inbred mouse strains45,46,47. Ten LG/J females were crossed with 10 SM/J males to produce 54 F1 hybrids. These F1 hybrids were intercrossed to produce 510 F2 animals. The F2 males and females were then reciprocally mated to produce 158 full-sibling F2 families with a total of 1,632 F3 progenies. Animals were sacrificed after 70 days of age. They were weighed to obtain an overall measure of body size (y). The length of the tail (x1) was measured with callipers. The mice were then immediately dissected by necropsy, and reproductive fat depot (x2), heart (x3), kidneys (x4), spleen (x5), and liver (x6) were weighed to the nearest 0.01 g with a digital scale.

A total of 353 single nucleotide polymorphism markers were chosen from the 4,200 polymorphic markers scored as part of the CTC/Oxford genotyping consortium. These makers were genotyped for all F2 animals and their F3 offspring. The F2 and F3 genotypes were used to reconstruct haplotypes using the “block-extension algorithm” in the PedPhase program48. With the inferred haplotype information, it was possible to distinguish all four genotypes in the F3 population at each marker locus with the paternal allele listed first and the maternal allele second.

Focusing on the form of absolute growth, Cheverud, et al.27 have mapped imprinting effects on the six tissues and organs in mice to know about the contribution of imprinting to quantitative variation in trait expression. By reanalysing the real dataset, we will estimate genomic imprinting effects and infer their patterns for relative growth of the adult body composition to body weight using the random regression model. Before gene mapping, these observations of traits were adjusted for the effects of ages at necropsy and litter sizes at birth45 and the residuals with population mean were used in the following analysis.

Stepwise regression analysis showed that not all six tissues and organs are significantly associated with body weight in relative growth. Thus, we dropped non-significant partial allometry exponents for model (2), and reached the following phenotypic joint static allometry scaling model by

An external file that holds a picture, illustration, etc.
Object name is srep45191-m20.jpg

where x4, x5 and x6 refer to allometry scaling variables for kidney, spleen and liver, respectively. The joint static allometry scaling model of kidney, spleen and liver were chosen for a genome-wide random regression analysis to infer imprinting allometries; fatpad, tail and heart were excluded from the model due to their non-significant phenotypic partial allometry scalings effects.

The gender variable was considered as fixed in the final random regression model of multiple static allometries. Each fixed effect and marker genetic effects were estimated for the partial allometry exponents, along with the covariance matrix for polygenic effects and residual variance in the random regression model were estimated using REML via the DMU package. The initial values were defaulted as zero for each fixed effect, as identity matrix for the additive genetic covariance matrix and as one for the residual variance. Convergence precision for REML was set to 10−6.

Results showed that a total of 122 QTNs were detected using the chromosome-wide tests. These QTNs were distributed on all chromosomes, and 48, 54 and 44 QTNs were associated with the relative growth of kidney, spleen and liver to body weight, respectively (results are shown in Table 1, Table 1S and Table 2S of the Supplementary file). Eleven QTNs were simultaneously associated with two of the three organs, showing pleiotropic effects. The profiles of test statistics for imprinting allometries are depicted in Fig. 1 for kidney (upper), spleen (middle) and liver (bottom). A total of 13, 15 and 11 markers were identified to be imprinted for relative weights of kidney, spleen and liver, respectively, because their test statistics exceeded the critical value of 1.301 at the 5% significance level, as displayed in Fig. 1. Using the chromosome-wide test, Table 1 tabulates the imprinted QTNs for relative growth of kidney, spleen and liver to body weight in mouse. For these imprinted QTNs, three show pleiotropy effects; that is, marker 103 (on chromosome 6) is associated with the relative weight of both spleen and liver, while marker 118 (on chromosome 8) and 206 (on chromosome 11) have effects on the relative weight of both kidney and liver. Genome-wide tests further verified that only 18 QTNs are inferred to be statistically significant at the 5% significance level, among which two are imprinted: marker 151 (on chromosome 15) regulates the relative weight of spleen and marker 184 (on chromosome 12) regulates the relative weight of liver.

Figure 1
The profiles of test statistics of imprinted QTNs for relative growth of kidney (a), spleen (b) and liver (c) to body weight in mouse. The horizontal line in each plot represents the critical value of 1.301 for −log(p).
Table 1
Chromosome-wide imprinted QTNs for relative growth of kidney, spleen and liver to body weight in mouse.

Most of the imprinted QTNs inherit in bipolar dominance fashion with no significant additive and dominant effects. A total of 5 QTNs with over-dominance pattern control relative growths of three organs to body weight, among which SNPs rs3683086, rs13482635 and rs3713033 on chromosome 11, 15 and 19, respectively, were correlated with the kidney trait, while rs6296621 on chromosome 13 and rs3683086 on chromosome 11 were associated with spleen and liver traits, respectively. In particular, allometry of kidney to body weight was found to be regulated by imprinted SNP rs3688854 (on chromosome 2) in an under-dominance fashion, and allometry of spleen was regulated by SNP rs13480638 (on chromosome 10) with a complete maternal additive imprinting pattern. SNP rs13475748 (on chromosome 1) was inferred to be imprinted for the relative growth of liver to body weight, but its imprinting pattern was not defined since both the additive and dominant effects were significant.

Table 2 provides the significant QTNs for relative growth of kidney, spleen and liver to body weight based on the genome-wide test results. Eighteen out of 122 QTNs passed the genome-wide test, among which 4, 9 and 5 were responsible for kidney, spleen and liver, respectively. Some QTNs showed an additive effect on allometries of three organs but no dominance effect (e.g., on liver). Two imprinted QTNs, rs13482486 (on chromosome 15) and rs3662939 (on chromosome 12) on spleen and liver, respectively, were further verified through the genome-wide test. Among the significant QTNs and polygene effects by genome-wide test, the two imprinted QTNs contributed 17.12% and 37.12% of the total genetic variances for the relative growth of spleen and liver, respectively.

Table 2
Genome-wide QTNs for relative growth of kidney, spleen and liver to body weight in mouse.

Simulation study

The purpose of simulation was to investigate the statistical behaviour of detecting QTNs inherited in different patterns with the genome-wide random regression analysis. Based on the results from real data analysis, five additive/dominant QTNs (on chromosome 2) and two imprinting QTNs (on chromosomes 12 and 15) were chosen for simulation analysis. In the simulation, the phenotypic values for kidney, spleen and liver were retained as well as the genotypes of genetic markers. Logarithm of body weights was generated by the estimated fix regression effects for sex An external file that holds a picture, illustration, etc.
Object name is srep45191-m21.jpg, residual variance An external file that holds a picture, illustration, etc.
Object name is srep45191-m22.jpg obtained from the real data analysis, and polygenic genetic covariance matrix An external file that holds a picture, illustration, etc.
Object name is srep45191-m23.jpg for multiple allometry scalings. Two simulation scenarios were considered. In scenario 1, we evaluated the precision of parameter estimation and the power of QTN detection by generating logarithm of body weights with An external file that holds a picture, illustration, etc.
Object name is srep45191-m24.jpg estimated from the real data analysis. In scenario 2, we investigated the performance of QTN detection with An external file that holds a picture, illustration, etc.
Object name is srep45191-m25.jpg. Simulations were repeated 100 times in each scenario to assess the power of QTN detection and precision of parameter estimation. Parameter estimates and statistical powers of QTL detection with genome-wide random regression analysis are shown in Table 3 for the simulated datasets. As can be seen, the higher the relative contribution of the simulated QTN, the greater the power to detect the QTN, which meets the general statistical behaviour in QTL mapping. When a QTN was detected, its genetic patterns can be further accurately inferred. In addition, we found that the power to detect the dominant and imprinting QTNs was low, as compared with that needed to identify additive QTNs. This observation is consistent with a regular QTL mapping study in which a QTL with an additive effect is easier to detect compared with the one with a dominant effect.

Table 3
Parameter estimates (standard deviations) and statistical powers of QTL detection with genome-wide random regression analysis for simulated datasets.

Discussion

Taking the joint static allometry scaling model49 as sub-model, we constructed the random regression model to statistically infer genomic parent-of-origin effects on the relative growth of body composition to body weight in mice. With comparison to the random regression model for growth and developmental traits42 and mapping procedure for allometry scalings30,31,32, there are three major advantages of our analysis method. First, the joint static allometry scaling model can more accurately estimate allometry scalings of multiple tissues and organs to body weights than the simple allometry equation16. More importantly, it facilitates the comparison and genetic analysis of multiple allometry scalings. Second, when analysing growth and developmental traits with the random regression model, it is required for each individual to repeatedly measure the traits in growth and developmental duration. However, such repeat measurement is not a necessity in our study, because for each individual, body compositions were measured only once at necropsy. Third, our method incorporates polygenic effects derived from a pedigree into a random regression model, improving the estimation accuracy of marker effects.

It should be noted that our method can provide the estimates of iQTNs’ effects and polygenic effects but not heritabilities of iQTNs; that is, it cannot answer how much iQTNs contribute to phenotypic variation. If experimental individuals are sacrificed at different ages, the residual covariance matrix for multiple allometries can be estimated by nesting the joint static allometry model into permanent environmental effects caused by multiple ages at necropsy. A complete random regression model can therefore be constructed to identify genomic imprinting for joint allometries, denoted as

An external file that holds a picture, illustration, etc.
Object name is srep45191-m26.jpg

where xipi are random family and permanent environmental effect on multiple allometries. This allows us to successfully evaluate the genetic variation of multiple allometries, such as estimation of heritability for allometry scaling of each tissue and organ to body weight. In addition, a complete pedigree of many families and more recorded individuals are required to stably estimate the parameters in such a complex model. When separating real dataset with gender in case analysis, sub-dataset from male population is even not convergent in REML, this limits discussion about the effect of gender on imprinting status.

With the same dataset, Cheverud, et al.27 examined the contribution of imprinting to quantitative variation in trait expression by estimating imprinting effects on absolute growths of adult body composition traits. Of the eight pleiotropic iQTL the authors identified, only those on chromosomes 7, 12, and centromeric 18 were located in regions previously reported containing imprinted genes50. Their findings of imprinting loci, effects and patterns on adult body compositions were strongly supported by genetic evidences of imprinting on chromosomes 751, 1252 and 1853. In our study, most iQTL for relative growth of body compositions are in new locations that have not previously been associated with imprinting effects on the absolute growth. Only six iQTLs for relative growth overlap with the four iQTLs identified for absolute growth on chromosomes 7, 12, and 18. In particular, the iQTLs for the relative growth on chromosomes 7, 12, and 18 have the same imprinting patterns as those for absolute growth, suggesting the existence of imprinting genes controlling both absolute and relative growth. By the chromosome-wide test, although the detected QTNs with additive, dominance or imprinting effect almost distributed on all chromosomes (see Tables 1S and 2S), only two iQTNs on chromosomes 12 and 15 were found to be associated with the relative growth of spleen and liver by the genome-wide test. Notably, the Pref-1/Dlk1 gene regulating growth retardation and accelerated adiposity is located on chromosome 12 in mice54.

If high density genetic markers are available that can distinguish the four genotypes, the method proposed here can be improved by doing an efficient marker selection, so that zero genetic effects can be first shrunk to zero by the least absolute shrinkage and selection operator55,56 for a sparse oversaturated regression model. Then, chromosome-wide or genome-wide non-zero genetic effects can be statistically inferred within the framework of the proposed model (5). If the detected QTNs only contribute to a small proportion of genetic variation for allometries, genome selection, for instance, whole genome regression57,58 and genomic best linear unbiased prediction59,60,61 can be introduced to assess genomic variation of allometries based on our proposed model with high density markers.

Additional Information

How to cite this article: Zhao, J. et al. Genome-wide random regression analysis for parent-of-origin effects of body composition allometries in mouse. Sci. Rep. 7, 45191; doi: 10.1038/srep45191 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:

Acknowledgments

We sincerely thank Prof. Jason B. Wolf for providing us with a real dataset. This study is partially supported by Special Scientific Research Funds for Central Non-profit Institutes, Chinese Academy of Fishery Sciences (2014B008 and 2014A04XK01) and the Chinese National Natural Science Foundation Grants (30972077, 31110103065 and 31371336).

Footnotes

The authors declare no competing financial interests.

Author Contributions R.Q.Y. proposed the method and supervised J.L.Z. and L.J.W. to write computer codes; S.L.L. and L.J. conducted real data analysis; and R.Q.Y. and Y.H.C. wrote the manuscript. All authors read and approved the manuscript.

References

  • Alleman M. & Doctor J. Genomic imprinting in plants: observations and evolutionary implications. Plant Mol Biol 43, 147–161 (2000). [PubMed]
  • McKeigue P W. S. Association of insulin gene VNTR polymorphism with polycystic ovary syndrome. The Lancet 349, 1771 (1997). [PubMed]
  • Falls J. G., Pulford D. J., Wylie A. A. & Jirtle R. L. Genomic Imprinting: Implications for Human Disease. Am J Pathol 154, 635–647 (1999). [PubMed]
  • Jeon J. T. et al. . A paternally expressed QTL affecting skeletal and cardiac muscle mass in pigs maps to the IGF2 locus. Nat. Genet. 21, 157–158 (1999). [PubMed]
  • Nezer C. et al. . An imprinted QTL with major effect on muscle mass and fat deposition maps to the IGF2 locus in pigs. Nat. Genet. 21, 155–156 (1999). [PubMed]
  • Paterson A. D., Naimark D. M. & Petronis A. The analysis of parental origin of alleles may detect susceptibility loci for complex disorders. Hum. Hered. 49, 197–204 (1999). [PubMed]
  • Mcinnis M. G. et al. . Genome-wide scan of bipolar disorder in 65 pedigrees: supportive evidence for linkage at 8q24, 18q22, 4q32, 2p12, and 13q12. Mol. Psychiatry 8, 288–298 (2003). [PubMed]
  • Van Laere A. S. et al. . A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. Nature 425, 832–836 (2003). [PubMed]
  • Tuiskulahaavisto M. et al. . Quantitative trait loci with parent-of-origin effects in chicken. Genetics Research 84, 57–66 (2004). [PubMed]
  • Tycko B. & Morison I. M. Physiological functions of imprinted genes. Physiol Rev 192, 245–258 (2002). [PubMed]
  • Constância M., Kelsey G. & Reik W. Resourceful imprinting. Nature 432, 53–57 (2004). [PubMed]
  • Isles A. R. & Holland A. J. Imprinted genes and mother–offspring interactions. Early Hum. Dev. 81, 73–77 (2005). [PubMed]
  • Sandovici I. et al. . Familial aggregation of abnormal methylation of parental alleles at the IGF2/H19 and IGF2R differentially methylated regions. Hum. Mol. Genet. 12, 1569–1578 (2003). [PubMed]
  • Naumova A. K. & Croteau S. Mechanisms of Epigenetic Variation: Polymorphic Imprinting. Curr. Genomics 5, 417–429 (2004).
  • Sandovici I. et al. . Interindividual variability and parent of origin DNA methylation differences at specific human Alu elements. Hum. Mol. Genet. 14, 2135–2143 (2005). [PubMed]
  • Huxley J. S. Problems of relative growth. (Methuen & Co., 1932).
  • Gao H., Liu Y., Zhang T., Yang R. & Yang H. Statistical models for jointly analyzing multiple allometries. Journal of theoretical biology 318, 205–209, doi: (2013).10.1016/j.jtbi.2012.11.001 [PubMed] [Cross Ref]
  • Cheverud J. M. Relationships among ontogenetic, static, and evolutionary allometry. Am J Phys Anthropol 59, 139–149 (1982). [PubMed]
  • Klingenberg C. P. & Zimmermann M. Static, Ontogenetic, and Evolutionary Allometry: A Multivariate Comparison in Nine Species of Water Striders. Am Nat 140, 601–620 (1992).
  • Stern D. L. & Emlen D. J. The developmental basis for allometry in insects. Development 126, 1091–1101 (1999). [PubMed]
  • Frankino W. A., Zwaan B. J., Stern D. L. & Brakefield P. M. Natural selection and developmental constraints in the evolution of allometries. Science 307, 718–720 (2005). [PMC free article] [PubMed]
  • Knott S. A., Elsen J. M. & Haley C. S. Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theoretical & Applied Genetics 93, 71–80 (1996). [PubMed]
  • de Koning D. J. et al. . Genome-wide scan for body composition in pigs reveals important role of imprinting. Proceedings of the National Academy of Sciences of the United States of America 97, 7947–7950 (2000). [PubMed]
  • de Koning D. J., Bovenhuis H. & van Arendonk J. A. On the detection of imprinted quantitative trait loci in experimental crosses of outbred species. Biomed. Chromatogr. 161, 931–938 (2002). [PubMed]
  • Cui Y., Lu Q., Cheverud J. M., Littell R. C. & Wu R. Model for mapping imprinted quantitative trait loci in an inbred F 2 design. Genomics 87, 543–551 (2006). [PubMed]
  • Yang R., Wang X., Wu Z., Prows D. R. & Lin M. Bayesian model selection for characterizing genomic imprinting effects and patterns. Bioinformatics 26, 235–241 (2010). [PMC free article] [PubMed]
  • Cheverud J. M. et al. . Genomic imprinting effects on adult body composition in mice. Proceedings of the National Academy of Sciences 105, 4253–4258 (2008). [PubMed]
  • Lande R. On Comparing Coefficients of Variation. Syst Zool 26, 214–217 (1977).
  • Mcguigan K., Nishimura N., Currey M., Dan H. & Cresko W. A. Quantitative genetic variation in static allometry in the threespine stickleback. Integr. Comp. Biol. 50, 1067–1080 (2010). [PubMed]
  • Rongling W. U., Chang-Xing M. A., Lou X. Y. & Casella G. Molecular Dissection of Allometry, Ontogeny, and Plasticity: A Genomic View of Developmental Biology. Bioscience 53, 1041–1047 (2009).
  • Wu R. & Hou W. A hyperspace model to decipher the genetic architecture of developmental processes: allometry meets ontogeny. Genetics 172, 627–637 (2006). [PubMed]
  • Li H. et al. . A conceptual framework for mapping quantitative trait Loci regulating ontogenetic allometry. PLoS ONE 2, e1245–e1245 (2007). [PMC free article] [PubMed]
  • Moore T. & Reik W. Genetic conflict in early development: parental imprinting in normal and abnormal growth. Rev. Reprod. 1, 73–77 (1996). [PubMed]
  • Itier J. M. et al. . Imprinted gene in postnatal growth role. Nature 393, 125–126 (1998). [PubMed]
  • Li L. et al. . Regulation of maternal behavior and offspring growth by paternally expressed Peg3. Science 284, 330–333 (1999). [PubMed]
  • Barker D. J. P. Obesity and early life. Obes. Rev. 8, 45–49 (2007). [PubMed]
  • Hager R. & Johnstone R. A. Early experience and parent-of-origin-specific effects influence female reproductive success in mice. Biol Lett 2, 253–256 (2006). [PMC free article] [PubMed]
  • Meaney M. J. & Szyf M. Maternal care as a model for experience-dependent chromatin plasticity? Trends Neurosci 28, 456–463 (2005). [PubMed]
  • Ehrich T. H. et al. . Pleiotropic effects on mandibular morphology I. developmental morphological integration and differential dominance. Journal of Experimental Zoology Part B Molecular & Developmental Evolution 296, 58–79 (2003). [PubMed]
  • Mantey C., Brockmann G. A., Kalm E. & Reinsch N. Mapping and exclusion mapping of genomic imprinting effects in mouse F2 families. J Hered 96, 329–338(310) (2005). [PubMed]
  • Henderson C. R. Analysis of covariance in the mixed model: higher-level, nonhomogeneous, and random regressions. Biometrics 38, 623–640 (1982). [PubMed]
  • Schaeffer L. R. Application of random regression models in animal breeding. Livestock Production Science 86, 35–45 (2004).
  • Patterson H. D. & Thompson R. Recovery of Inter-Block Information when Block Sizes are Unequal. Biometrika 58, 545–554 (1971).
  • Wolf J. B., Cheverud J. M., Roseman C. & Hager R. Genome-wide analysis reveals a complex pattern of genomic imprinting in mice. PLoS Genet 4 (2008). [PMC free article] [PubMed]
  • Kramer et al. . Genetic variation in body weight gain and composition in the intercross of Large (LG/J) and Small (SM/J) inbred strains of mice. Genet. Mol. Biol. 21, 706–712 (1998).
  • Vaughn T. T. et al. . Mapping quantitative trait loci for murine growth: a closer look at genetic architecture. Genet Res 74, 313–322 (1999). [PubMed]
  • Wolf J. B., Vaughn T. T., Pletscher L. S. & Cheverud J. M. Contribution of maternal effect QTL to genetic architecture of early growth in mice. Heredity 89, 300–310 (2002). [PubMed]
  • Li J. & Jiang T. Efficient inference of haplotypes from genotypes on a pedigree. Journal of Bioinformatics & Computational Biology 1, 41–69 (2009). [PubMed]
  • Gao H., Liu Y., Zhang T., Yang R. & Prows D. R. Parametric proportional hazards model for mapping genomic imprinting of survival traits. Journal of Applied Genetics 54, 79–88 (2013). [PubMed]
  • Luedi P. P., Hartemink A. J. & Jirtle R. L. Genome-wide prediction of imprinted murine genes. Genome Res. 15, 875–884 (2005). [PubMed]
  • Dong C. et al. . Possible genomic imprinting of three human obesity-related genetic loci. Am. J. Hum. Genet. 76, 427–437 (2005). [PubMed]
  • Watanabe T., Tarttelin E., Neubüser A., Kimura M. & Solter D. Fine genetic mapping defines the genetic order of Pax9, Tcf3a, and Acrodysplasia (Adp). Mamm Genome 5, 768–770 (1994). [PubMed]
  • Hagiwara Y. et al. . Screening for imprinted genes by allelic message display: Identification of a paternally expressed gene Impact on mouse chromosome 18. Proc Natl Acad Sci USA 94, 9249–9254 (1997). [PubMed]
  • Yang S. M. et al. . Mice Lacking Paternally Expressed Pref-1/Dlk1 Display Growth Retardation and Accelerated Adiposity. Molecular & Cellular Biology 22, 5585–5592 (2002). [PMC free article] [PubMed]
  • Friedman J., Hastie T. & Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent. Journal of Statistical Software 33, 1–22 (2010). [PMC free article] [PubMed]
  • Tibshirani R. Regression shrinkage and selection via the lasso: a retrospective. Journal of the Royal Statistical Society 73, 273–282 (2011).
  • Meuwissen T. H., Hayes B. J. & Goddard M. E. Prediction of total genetic value using genome-wide dense marker maps. Genetics 157, 1819–1829 (2001). [PubMed]
  • Gianola D. Priors in whole-genome regression: the bayesian alphabet returns. Genetics 194, 573–596 (2013). [PubMed]
  • Visscher P. M. et al. . Assumption-Free Estimation of Heritability from Genome-Wide Identity-by-Descent Sharing between Full Siblings. PLoS Genet 2, 316–325 (2006). [PMC free article] [PubMed]
  • Vanraden P. M. Efficient Methods to Compute Genomic Predictions. J Dairy Sci 91, 4414–4423 (2008). [PubMed]
  • Hayes B. J., Bowman P. J., Daetwyler H. D., Kijas J. W. & Van d. W. J. H. J. Accuracy of genotype imputation in sheep breeds. Anim Genet 43, 72–80 (2012). [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group