Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5364555

Formats

Article sections

Authors

Related links

Sci Rep. 2017; 7: 45191.

Published online 2017 March 24. doi: 10.1038/srep45191

PMCID: PMC5364555

Jingli Zhao,^{1,}^{2,}^{*} Shuling Li,^{3,}^{*} Lijuan Wang,^{4} Li Jiang,^{2} Runqing Yang,^{a,}^{1} and Yuehua Cui^{b,}^{5,}^{6}

Received 2016 June 17; Accepted 2017 February 22.

Copyright © 2017, The Author(s)

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

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 plants^{1} and animals^{2}^{,3}^{,4}^{,5}^{,6}^{,7}^{,8}^{,9} and has been recognized for its role in shaping developmental processes^{10}^{,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 levels^{13}^{,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 scaling^{16}. 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 model^{17} 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 allometry^{18}^{,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 evolve^{21}.

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 squares^{9}^{,22}^{,23}^{,24} or maximum likelihood methods^{25}. Multi-step tests for contrast models have been also proposed to identify the imprinting pattern^{9}^{,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 model^{28}^{,29} and into genotypic effects of genetic model for mapping QTL^{30}^{,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 mammals^{33}^{,34}^{,35}. Imprinting also affects adult traits either through the persistent effects of early growth and development^{36}^{,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 strains^{39}. 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.

Let denote *m* partial body sizes and *y* be the entire body size. The joint static allometry model^{17} is defined as

where *β*_{0} is an intercept, and are partial scaling exponents of the *j*th 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

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

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

where *y*_{i} is the entire body size for the *i*th individual; *b*_{i} represents the *l*th fixed effect in *l* systematic environments, such as sexes; *a*_{j}, *d*_{j} and *i*_{j} are additive, dominance and imprinting genetic effects for the *j*th marker; *h*_{il}, *z*_{ij}, *w*_{ij} and *s*_{ij} are indicator variables corresponding to *b*_{l}, *a*_{j}, *d*_{j} and *i*_{j}, where *z*_{ij}, *w*_{ij} and *s*_{ij} are defined in Mantey *et al*.^{40}; *g*_{i} is the polygenic genetic effect derived from the pedigree information, assuming that with A being a relationship matrix and being the polygenic genetic variance; and *e*_{i} is the residual error with .

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 model^{41}^{,42}:

where with *θ**b, a, d, i* or *g*. Let , with **h**, **z**, **w** or **s** and *h, z, w* or *s*, , with **θ****b**, **a**, **d** or **i** and , then in matrix form, model (3) can be rewritten as,

In this model, genetic effects are fixed for all markers and polygenic effects are random for all individuals. Thus, and , where **A** and **G** are the numerator-relationship matrix and polygenic genetic covariance matrix for multiple allometric scalings, respectively.

Restricted maximum likelihood method^{43} 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

Similarly, another *t* statistic is calculated as

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 (H_{0}: *a*=0 and *d*=0), polar over-dominance (H_{0}: *a*=0 and *d*=*i*) and polar under-dominance (H_{0}: *a*=0 and *d*=−*i*). Following the definition of different imprinting types and the corresponding null hypothesis^{26}^{,27}^{,44}, the imprinting pattern for the detected QTL can be statistically inferred by formula (7).

A set of reciprocal cross families for distinguishing the four genotypes considering parental origin was derived from an F_{2} intercross of large (LG/J) and small (SM/J) inbred mouse strains^{45}^{,46}^{,47}. Ten LG/J females were crossed with 10 SM/J males to produce 54 F_{1} hybrids. These F_{1} hybrids were intercrossed to produce 510 F_{2} animals. The F_{2} males and females were then reciprocally mated to produce 158 full-sibling F_{2} families with a total of 1,632 F_{3} 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 (*x*_{1}) was measured with callipers. The mice were then immediately dissected by necropsy, and reproductive fat depot (*x*_{2}), heart (*x*_{3}), kidneys (*x*_{4}), spleen (*x*_{5}), and liver (*x*_{6}) were weighed to the nearest 0.01g 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 F_{2} animals and their F_{3} offspring. The F_{2} and F_{3} genotypes were used to reconstruct haplotypes using the “block-extension algorithm” in the PedPhase program^{48}. With the inferred haplotype information, it was possible to distinguish all four genotypes in the F_{3} 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 birth^{45} 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

where *x*_{4}, *x*_{5} and *x*_{6} 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.

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*).

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.

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 , residual variance obtained from the real data analysis, and polygenic genetic covariance matrix 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 estimated from the real data analysis. In scenario 2, we investigated the performance of QTN detection with . 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.

Taking the joint static allometry scaling model^{49} 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 traits^{42} and mapping procedure for allometry scalings^{30}^{,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 equation^{16}. 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

where **x**_{i}**p**_{i} 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 genes^{50}. Their findings of imprinting loci, effects and patterns on adult body compositions were strongly supported by genetic evidences of imprinting on chromosomes 7^{51}, 12^{52} and 18^{53}. 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 mice^{54}.

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 operator^{55}^{,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 regression^{57}^{,58} and genomic best linear unbiased prediction^{59}^{,60}^{,61} can be introduced to assess genomic variation of allometries based on our proposed model with high density markers.

**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.

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).

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.

- 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 chromosome18. 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**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |