A statistical model for dissecting genomic imprinting through genetic mapping

As a result of nonequivalent genetic contribution of maternal and paternal genomes to offsprings, genomic imprinting or called parent-of-origin effect, has been broadly identified in plants, animals and humans. Its role in shaping organism’s development has been unanimously recognized. However, statistical methods for identifying imprinted quantitative trait loci (iQTL) and estimating the imprinted effect have not been well developed. In this article, we propose an efficient statistical procedure for genomewide estimating and testing the effects of significant iQTL underlying the quantitative variation of interested traits. The developed model can be applied to two different genetic cross designs, backcross and F₂ families derived from inbred lines. The proposed procedure is built within the maximum likelihood framework and implemented with the EM algorithm. Extensive simulation studies show that the proposed model is well performed in a variety of situations. To demonstrate the usefulness of the proposed approach, we apply the model to a published data in an F₂ family derived from LG/S and SM/S mouse stains. Two partially maternal imprinting iQTL are identified which regulate the growth of body weight. Our approach provides a testable framework for identifying and estimating iQTL involved in the genetic control of complex traits.     

A statistical framework for genome-wide scanning and testing of imprinted quantitative trait loci

Non-equivalent expression of alleles at a locus results in genomic imprinting. In this article, a statistical framework for genome-wide scanning and testing of imprinted quantitative trait loci (iQTL) underlying complex traits is developed based on experimental crosses of inbred line species in backcross populations. The joint likelihood function is composed of four component likelihood functions with each of them derived from one of four backcross families. The proposed approach models genomic imprinting effect as a probability measure with which one can test the degree of imprinting. Simulation results show that the model is robust for identifying iQTL with various degree of imprinting ranging from no imprinting, partial imprinting to complete imprinting. Under various simulation scenarios, the proposed model shows consistent parameter estimation with reasonable precision and high power in testing iQTL. When a QTL shows Mendelian effect, the proposed model also outperforms traditional Mendelian model. Extension to incorporate maternal effect is also given. The developed model, built within the maximum likelihood framework and implemented with the EM algorithm, provides a quantitative framework for testing and estimating iQTL involved in the genetic control of complex traits.     

Functional mapping imprinted quantitative trait loci underlying developmental characteristics

Background

Genomic imprinting, a phenomenon referring to nonequivalent expression of alleles depending on their parental origins, has been widely observed in nature. It has been shown recently that the epigenetic modification of an imprinted gene can be detected through a genetic mapping approach. Such an approach is developed based on traditional quantitative trait loci (QTL) mapping focusing on single trait analysis. Recent studies have shown that most imprinted genes in mammals play an important role in controlling embryonic growth and post-natal development. For a developmental character such as growth, current approach is less efficient in dissecting the dynamic genetic effect of imprinted genes during individual ontology.

Results

Functional mapping has been emerging as a powerful framework for mapping quantitative trait loci underlying complex traits showing developmental characteristics. To understand the genetic architecture of dynamic imprinted traits, we propose a mapping strategy by integrating the functional mapping approach with genomic imprinting. We demonstrate the approach through mapping imprinted QTL controlling growth trajectories in an inbred F₂ population. The statistical behavior of the approach is shown through simulation studies, in which the parameters can be estimated with reasonable precision under different simulation scenarios. The utility of the approach is illustrated through real data analysis in an F₂ family derived from LG/J and SM/J mouse stains. Three maternally imprinted QTLs are identified as regulating the growth trajectory of mouse body weight.

Conclusion

The functional iQTL mapping approach developed here provides a quantitative and testable framework for assessing the interplay between imprinted genes and a developmental process, and will have important implications for elucidating the genetic architecture of imprinted traits.     

Parametric proportional hazards model for mapping genomic imprinting of survival traits

A number of imprinted genes have been observed in plants, animals and humans. They not only control growth and developmental traits, but may also be responsible for survival traits. Based on the Cox proportional hazards (PH) model, we constructed a general parametric model for dissecting genomic imprinting, in which a baseline hazard function is selectable for fitting the effects of imprinted quantitative trait loci (iQTL) genotypes on the survival curve. The expectation–maximisation (EM) algorithm is derived for solving the maximum likelihood estimates of iQTL parameters. The imprinting patterns of the detected iQTL are statistically tested under a series of null hypotheses. The Bayesian information criterion (BIC) model selection criterion is employed to choose an optimal baseline hazard function with maximum likelihood and parsimonious parameterisation. We applied the proposed approach to analyse the published data in an F₂ population of mice and concluded that, among five commonly used survival distributions, the log-logistic distribution is the optimal baseline hazard function for the survival time of hyperoxic acute lung injury (HALI). Under this optimal model, five QTL were detected, among which four are imprinted in different imprinting patterns.     

Bayesian inference for genomic imprinting underlying developmental characteristics

The identification of imprinted genes is becoming a standard procedure in searching for quantitative trait loci (QTL) underlying complex traits. When a developmental characteristic such as growth or drug response is observed at multiple time points, understanding the dynamics of gene function governing the underlying feature should provide more biological information regarding the genetic control of an organism. Recognizing that differential imprinting can be development-specific, mapping imprinted genes considering the dynamic imprinting effect can provide additional biological insights into the epigenetic control of a complex trait. In this study, we proposed a Bayesian imprinted QTL (iQTL) mapping framework considering the dynamics of imprinting effects and model multiple iQTLs with an efficient Bayesian model selection procedure. The method overcomes the limitation of likelihood-based mapping procedure, and can simultaneously identify multiple iQTLs with different gene action modes across the whole genome with high computational efficiency. An inference procedure using Bayes factors to distinguish different imprinting patterns of iQTL was proposed. Monte Carlo simulations were conducted to evaluate the performance of the method. The utility of the approach was illustrated through an analysis of a body weight growth data set in an F(2) family derived from LG/J and SM/J mouse stains. The proposed Bayesian mapping method provides an efficient and computationally feasible framework for genome-wide multiple iQTL inference with complex developmental traits.     

A random model for mapping imprinted quantitative trait loci in a structured pedigree: an implication for mapping canine hip dysplasia

Genetic imprinting may have played a more notable role in shaping embryonic development of plants, animals, and humans than previously appreciated. Quantitative trait loci that are imprinted (iQTL) exert monoallelic effects, depending on the parent of origin, which is an exception to the laws of Mendelian genetics. In this article, we present a modified random effect-based mapping model to use in a genome-wide scan for the distribution of iQTL that contribute to genetic variance for a complex trait in a structured pedigree. This model, implemented with the maximum likelihood method, capitalizes on a network of relatedness for maternally and paternally derived alleles through identical-by-descent sharing, thus allowing for the discrimination of the genetic variances due to alleles derived from maternal and paternal parents. The model was employed to map iQTL responsible for canine hip dysplasia in a multihierarchical canine pedigree, founded with seven greyhounds and six Labrador retrievers. Of eight significant QTL detected, three, located on CFA1, CFA8, and CF28, were found to trigger significant parent-of-origin effects on the age of femoral capital ossification measured at the left and right hips of a canine. The detected iQTL provide important candidate regions for fine-mapping of imprinted genes and for studying their structure and function in the control of complex traits.     

Bayesian mapping of genome-wide epistatic imprinted loci for quantitative traits

Genomic imprinting, an epigenetic phenomenon of parent-of-origin-specific gene expression, has been widely observed in plants, animals, and humans. To detect imprinting genes influencing quantitative traits, the least squares and maximum likelihood approaches for fitting a single quantitative trait locus (QTL) and Bayesian methods for simultaneously modeling multiple QTL have been adopted, respectively, in various studies. However, most of these studies have only estimated imprinting main effects and thus ignored imprinting epistatic effects. In the presence of extremely complex genomic imprinting architectures, we introduce a Bayesian model selection method to analyze the multiple interacting imprinted QTL (iQTL) model. This approach will greatly enhance the computational efficiency through setting the upper bound of the number of QTLs and performing selective sampling for QTL parameters. The imprinting types of detected main-effect QTLs can be estimated from the Bayes factor statistic formulated by the posterior probabilities for the genetic effects being compared. The performance of the proposed method is demonstrated by several simulation experiments. Moreover, this method is applied to dissect the imprinting genetic architecture for body weight in mouse and fruit weight in tomato. Matlab code for implementing this approach will be available from the authors upon request.     

A statistical design for testing transgenerational genomic imprinting in natural human populations

Genomic imprinting is a phenomenon in which the same allele is expressed differently, depending on its parental origin. Such a phenomenon, also called the parent-of-origin effect, has been recognized to play a pivotal role in embryological development and pathogenesis in many species. Here we propose a statistical design for detecting imprinted loci that control quantitative traits based on a random set of three-generation families from a natural population in humans. This design provides a pathway for characterizing the effects of imprinted genes on a complex trait or disease at different generations and testing transgenerational changes of imprinted effects. The design is integrated with population and cytogenetic principles of gene segregation and transmission from a previous generation to next. The implementation of the EM algorithm within the design framework leads to the estimation of genetic parameters that define imprinted effects. A simulation study is used to investigate the statistical properties of the model and validate its utilization. This new design, coupled with increasingly used genome-wide association studies, should have an immediate implication for studying the genetic architecture of complex traits in humans.     

A Model for Transgenerational Imprinting Variation in Complex Traits

Despite the fact that genetic imprinting, i.e., differential expression of the same allele due to its different parental origins, plays a pivotal role in controlling complex traits or diseases, the origin, action and transmission mode of imprinted genes have still remained largely unexplored. We present a new strategy for studying these properties of genetic imprinting with a two-stage reciprocal F mating design, initiated with two contrasting inbred lines. This strategy maps quantitative trait loci that are imprinted (i.e., iQTLs) based on their segregation and transmission across different generations. By incorporating the allelic configuration of an iQTL genotype into a mixture model framework, this strategy provides a path to trace the parental origin of alleles from previous generations. The imprinting effects of iQTLs and their interactions with other traditionally defined genetic effects, expressed in different generations, are estimated and tested by implementing the EM algorithm. The strategy was used to map iQTLs responsible for survival time with four reciprocal F populations and test whether and how the detected iQTLs inherit their imprinting effects into the next generation. The new strategy will provide a tool for quantifying the role of imprinting effects in the creation and maintenance of phenotypic diversity and elucidating a comprehensive picture of the genetic architecture of complex traits and diseases.     

Characterization of genomic imprinting effects and patterns with parametric accelerated failure time model

Genomic imprinting, a genetic phenomenon of non-equivalent allele expression that depends on parental origins, has been ubiquitously observed in nature. It does not only control the traits of growth and development but also may be responsible for survival traits. Based on the accelerated failure time model, we construct a general parametric model for mapping the imprinted QTL (iQTL). Within the framework of interval mapping, maximum likelihood estimation of iQTL parameters is implemented via EM algorithm. The imprinting patterns of the detected iQTL are statistically tested according to a series of null hypotheses. BIC model selection criterion is employed to choose an optimal baseline hazard function with maximum likelihood and parsimonious parameters. Simulations are used to validate the proposed mapping procedure. A published dataset from a mouse model system was used to illustrate the proposed framework. Results show that among the five commonly used survival distributions, Log-logistic distribution is the optimal baseline hazard function for mapping QTL of hyperoxic acute lung injury (HALI) survival; under the log-logistic distribution, four QTLs were identified, in which only one QTL was inherited in Mendelian fashion, whereas others were imprinted in different imprinting patterns.     

Functional mapping of reaction norms to multiple environmental signals

Whether there are different genes involved in response to different environmental signals and how these genes interact to determine the final expression of the trait are of fundamental importance in agricultural and biological research. We present a statistical framework for mapping environment-induced genes (or quantitative trait loci, QTLs) of major effects on the expression of a trait that respond to changing environments. This framework is constructed with a maximum-likelihood-based mixture model, in which the mean and covariance structure of environment-induced responses is modelled. The means for responses to continuous environmental states, referred to as reaction norms, are approximated for different QTL genotypes by mathematical equations that were derived from fundamental biological principles or based on statistical goodness-of-fit to observational data. The residual covariance between different environmental states was modelled by autoregressive processes. Such an approach to studying the genetic control of reaction norms can be expected to be advantageous over traditional mapping approaches in which no biological principles and statistical structures are considered. We demonstrate the analytical procedure and power of this approach by modelling the photosynthetic rate process as a function of temperature and light irradiance. Our approach allows for testing how a QTL affects the reaction norm of photosynthetic rate to a specific environment and whether there exist different QTLs to mediate photosynthetic responses to temperature and light irradiance, respectively.     

Modelling epistatic effects of embryo and endosperm QTL on seed quality traits

Coordinated expression of embryo and endosperm tissues is required for proper seed development. The coordination among these two tissues is controlled by the interaction between multiple genes expressed in the embryo and endosperm genomes. In this article, we present a statistical model for testing whether quantitative trait loci (QTL) active in different genomes, diploid embryo and triploid endosperm, epistatically affect a trait expressed on the endosperm tissue. The maximum likelihood approach, implemented with the EM algorithm, was derived to provide the maximum likelihood estimates of the locations of embryo- and endosperm-specific QTL and their main effects and epistatic effects. This model was used in a real example for rice in which two QTL, one from the embryo genome and the other from the endosperm genome, exert a significant interaction effect on gel consistency on the endosperm. Our model has successfully detected Waxy, a candidate gene in the embryo genome known to regulate one of the major steps of amylose biosynthesis in the endosperm. This model will have great implications for agricultural and evolutionary genetic research.     

Functional mapping - how to map and study the genetic architecture of dynamic complex traits

The development of any organism is a complex dynamic process that is controlled by a network of genes as well as by environmental factors. Traditional mapping approaches for analysing phenotypic data measured at a single time point are too simple to reveal the genetic control of developmental processes. A general statistical mapping framework, called functional mapping, has been proposed to characterize, in a single step, the quantitative trait loci (QTLs) or nucleotides (QTNs) that underlie a complex dynamic trait. Functional mapping estimates mathematical parameters that describe the developmental mechanisms of trait formation and expression for each QTL or QTN. The approach provides a useful quantitative and testable framework for assessing the interplay between gene actions or interactions and developmental changes.     

A robust multiple-locus method for quantitative trait locus analysis of non-normally distributed multiple traits

Linear regression-based quantitative trait loci/association mapping methods such as least squares commonly assume normality of residuals. In genetics studies of plants or animals, some quantitative traits may not follow normal distribution because the data include outlying observations or data that are collected from multiple sources, and in such cases the normal regression methods may lose some statistical power to detect quantitative trait loci. In this work, we propose a robust multiple-locus regression approach for analyzing multiple quantitative traits without normality assumption. In our method, the objective function is least absolute deviation (LAD), which corresponds to the assumption of multivariate Laplace distributed residual errors. This distribution has heavier tails than the normal distribution. In addition, we adopt a group LASSO penalty to produce shrinkage estimation of the marker effects and to describe the genetic correlation among phenotypes. Our LAD-LASSO approach is less sensitive to the outliers and is more appropriate for the analysis of data with skewedly distributed phenotypes. Another application of our robust approach is on missing phenotype problem in multiple-trait analysis, where the missing phenotype items can simply be filled with some extreme values, and be treated as outliers. The efficiency of the LAD-LASSO approach is illustrated on both simulated and real data sets.     

Cloning quantitative trait loci by insertional mutagenesis

Summary This study explores the theoretical potential of insertional mutagenesis (i.e., mutagenesis as a result of integration of novel DNA sequences into the germ line), as a means of cloning quantitative trait loci (QTL). The approach presented is based on a direct search for mutagenic effects of a quantitative nature, and makes no assumptions as to the nature of the loci affecting quantitative trait value. Since there are a very large number of potential insertion sites in the genome but only a limited number of target sites that can affect any particular trait, large numbers of inserts will have to be generated and screened. The effects of allelic variants at any single QTL on phenotype value are expected to be small relative to sampling variation. Thus two of three stages of replicate testing will be required for each insert in order to bring overall Type I error down to negligible proportions and yet maintain good statistical power for inserts with true effects on the quantitative traits under consideration. The overall effort involved will depend on the spectrum of mutagenic effects produced by insertional mutagenesis. This spectrum is presently unknown, but using reasonable estimates, about 10,000 inserts would have to be tested, at reasonable replicate numbers (n 30) and Type I error (=0.01) in the first testing stage, to provide a high likelihood of detecting at least one insert with a true effect on a given quantitative trait of interest. Total offspring numbers required per true quantitative mutagenic effect detected decrease strongly with increased number of traits scored and increased number of inserts per initial transformed parent. In fact, it would appear that successful implementation of experiments of this sort will require the introduction of multiple independent inserts in the original parent individuals, by means of repeated transformation, or use of transposable elements as inserts. When biologically feasible, selfing would appear to be the method of choice for insert replication, and in all cases the experiments must be carried out in inbred lines to reduce error variation due to genetic segregation, and avoid confounding mutational effects of the insert with effects due to linkage with nearby segregating QTL. The special qualifications of Arabidopsis thaliana for studies of this sort are emphasized, and problems raised by somaclonal variation are discussed.Contribution from the Agricultural Research Organization, The Volcani Center, Bet Dagan, Israel. No. 1834-E, 1986 series     

Quantitative trait nucleotide analysis using Bayesian model selection

Although much attention has been given to statistical genetic methods for the initial localization and fine mapping of quantitative trait loci (QTLs), little methodological work has been done to date on the problem of statistically identifying the most likely functional polymorphisms using sequence data. In this paper we provide a general statistical genetic framework, called Bayesian quantitative trait nucleotide (BQTN) analysis, for assessing the likely functional status of genetic variants. The approach requires the initial enumeration of all genetic variants in a set of resequenced individuals. These polymorphisms are then typed in a large number of individuals (potentially in families), and marker variation is related to quantitative phenotypic variation using Bayesian model selection and averaging. For each sequence variant a posterior probability of effect is obtained and can be used to prioritize additional molecular functional experiments. An example of this quantitative nucleotide analysis is provided using the GAW12 simulated data. The results show that the BQTN method may be useful for choosing the most likely functional variants within a gene (or set of genes). We also include instructions on how to use our computer program, SOLAR, for association analysis and BQTN analysis.     

Maternal effects as the cause of parent-of-origin effects that mimic genomic imprinting

Epigenetic effects are increasingly recognized as an important source of variation in complex traits and have emerged as the focus of a rapidly expanding area of research. Principle among these effects is genomic imprinting, which has generally been examined in analyses of complex traits by testing for parent-of-origin-dependent effects of alleles. However, in most of these analyses maternal effects are confounded with genomic imprinting because they can produce the same patterns of phenotypic variation expected for various forms of imprinting. Distinguishing between the two is critical for genetic and evolutionary studies because they have entirely different patterns of gene expression and evolutionary dynamics. Using a simple single-locus model, we show that maternal genetic effects can result in patterns that mimic those expected under genomic imprinting. We further demonstrate how maternal effects and imprinting effects can be distinguished using genomic data from parents and offspring. The model results are applied to a genome scan for quantitative trait loci (QTL) affecting growth- and weight-related traits in mice to illustrate how maternal effects can mimic imprinting. This genome scan revealed five separate maternal-effect loci that caused a diversity of patterns mimicking those expected under various modes of genomic imprinting. These results demonstrate that the appearance of parent-of-origin-dependent effects (POEs) of alleles at a locus cannot be taken as direct evidence that the locus is imprinted. Moreover, they show that, in gene mapping studies, genetic data from both parents and offspring are required to successfully differentiate between imprinting and maternal effects as the cause of apparent parent-of-origin effects of alleles.     

Mitochondrial genetic effects on latent class variables associated with susceptibility to alcoholism

We report the results of statistical genetic analyses of data from the Collaborative Study on the Genetics of Alcoholism prepared for the Genetic Analysis Workshop 14 to detect and characterize maternally inherited mitochondrial genetic effects on variation in latent class psychiatric/behavioral variables employed in the diagnosis of alcoholism. Using published extensions to variance decomposition methods for statistical genetic analysis of continuous and discrete traits we: 1) estimated the proportion of the variance in each trait due to the effects of mitochondrial DNA (mtDNA), 2) tested for pleiotropy, both mitochondrial genetic and residual additive genetic, between trait pairs, and 3) evaluated whether the simultaneous estimation of mitochondrial genetic effects on these traits improves our ability to detect and localize quantitative trait loci (QTL) in the nuclear genome. After correction for multiple testing, we find significant (p < 0.009) mitochondrial genetic contributions to the variance for two latent class variables. Although we do detect significant residual additive genetic correlations between the two traits, there is no evidence of a residual mitochondrial genetic correlation between them. Evidence for autosomal QTL for these traits is improved when linkage screens are conditioned on significant mitochondrial genetic effects. We conclude that mitochondrial genes may contribute to variation in some latent class psychiatric/behavioral variables associated with alcoholism.