A framework for detecting and characterizing genetic background-dependent imprinting effects

Parent-of-origin effects propagate through networks to shape metabolic traits

Parent-of-origin effects are unexpectedly common in complex traits, including metabolic and neurological traits. Parent-of-origin effects can be modified by the environment, but the architecture of these gene-by-environmental effects on phenotypes remains to be unraveled. Previously, quantitative trait loci (QTL) showing context-specific parent-of-origin effects on metabolic traits were mapped in the F16 generation of an advanced intercross between LG/J and SM/J inbred mice. However, these QTL were not enriched for known imprinted genes, suggesting another mechanism is needed to explain these parent-of-origin effects phenomena. We propose that non-imprinted genes can generate complex parent-of-origin effects on metabolic traits through interactions with imprinted genes. Here, we employ data from mouse populations at different levels of intercrossing (F0, F1, F2, F16) of the LG/J and SM/J inbred mouse lines to test this hypothesis. Using multiple populations and incorporating genetic, genomic, and physiological data, we leverage orthogonal evidence to identify networks of genes through which parent-of-origin effects propagate. We identify a network comprised of three imprinted and six non-imprinted genes that show parent-of-origin effects. This epistatic network forms a nutritional responsive pathway and the genes comprising it jointly serve cellular functions associated with growth. We focus on two genes, Nnat and F2r, whose interaction associates with serum glucose levels across generations in high-fat-fed females. Single-cell RNAseq reveals that Nnat expression increases and F2r expression decreases in pre-adipocytes along an adipogenic trajectory, a result that is consistent with our observations in bulk white adipose tissue.

Fruit quality and DNA methylation are affected by parental order in reciprocal crosses of tomato

Reciprocal effects were found for tomato fruit quality and DNA methylation. The epigenetic identity of reciprocal hybrids indicates that DNA methylation might be one of the mechanisms involved in POEs. Crosses between different genotypes and even between different species are commonly used in plant breeding programs. Reciprocal hybrids are obtained by changing the cross direction (or the sexual role) of parental genotypes in a cross. Phenotypic differences between these hybrids constitute reciprocal effects (REs). The aim of this study was to evaluate phenotypic differences in tomato fruit traits and DNA methylation profiles in three inter- and intraspecific reciprocal crosses. REs were detected for 13 of the 16 fruit traits analyzed. The number of traits with REs was the lowest in the interspecific cross, whereas the highest was found in the cross between recombinant inbred lines (RILs) derived from the same interspecific cross. An extension of gene action analysis was proposed to incorporate parent-of-origin effects (POEs). Maternal and paternal dominance were found in four fruit traits. REs and paternal inheritance were found for epiloci located at coding and non-coding regions. The epigenetic identity displayed by the reciprocal hybrids accounts for the phenotypic differences among them, indicating that DNA methylation might be one of the mechanisms involved in POEs.

Joint modeling of eQTLs and parent-of-origin effects using an orthogonal framework with RNA-seq data

Extensive studies have been conducted on the analysis of genome function, especially on the expression quantitative trait loci (eQTL). These studies offered promising results for characterization of the functional sequencing variation and understanding of the basic processes of gene regulation. Parent of origin effect (POE) is an important epigenetic phenomenon describing that the expression of certain genes depends on their allelic parent-of-origin and it is known to play important roles in human complex diseases. However, traditional eQTL mapping approaches do not allow for the detection of imprinting, or they focus on modeling the additive genetic effect thereby ignoring the estimation of the dominance genetic effect. In this study, we proposed a statistical framework to test the additive and dominance genetic effects of the candidate eQTLs along with detection of the POE with a functional model and an orthogonal model for RNA-seq data. We demonstrated the desirable power and preserved Type I errors of the methods in most scenarios, especially the orthogonal model with un-biased estimation of the genetic effects and over-dispersion of the RNA-seq data. The application to a HapMap project trio dataset validated existing imprinting genes and discovered two novel imprinting genes with potential dominance genetic effect and RB1 and IGF1R genes. This study provides new insights into the next generation statistical modeling of eQTL mapping for better understanding of the genetic architecture underlying the mechanisms of gene expression regulation.

Evolutionary Quantitative Genetics of Genomic Imprinting

Genomic imprinting shapes the genotype-phenotype relationship by creating an asymmetry between the influences of paternally and maternally inherited gene copies. Consequently, imprinting can impact heritable and nonheritable variation, resemblance of relatives, and evolutionary dynamics. Although previous analyses have identified some of the quantitative genetic consequences of imprinting, we lack a framework that cleanly separates the influence of imprinting from other components of variation, particularly dominance. Here we apply a simple orthogonal genetic model to evaluate the roles of genetic (additive and dominance) and epigenetic (imprinting) effects. Imprinting increases the resemblance of relatives who share the expressed allele, and therefore increases variance among families of full or half-siblings. However, only part of this increased variance is heritable and contributes to selection responses. When selection is within, or among, families sharing only a single parent (half-siblings), which is common in selective breeding programs, imprinting can alter overall responses. Selection is more efficient when it acts among families sharing the expressed parent, or within families sharing the parent with lower expression. Imprinting also affects responses to sex-specific selection. When selection is on the sex whose gene copy has lower expression, the response is diminished or delayed the next generation, although the long-term response is unaffected. Our findings have significant implications for understanding patterns of variation, interpretation of short-term selection responses, and the efficacy of selective breeding programs, demonstrating the importance of considering the independent influence of genomic imprinting in quantitative genetics.

Offspring genes indirectly influence sibling and maternal behavioural strategies over resource share

Family members show behavioural strategies predicted to maximize individual fitness. These behaviours depend directly on genes expressed in focal individuals but also indirectly on genes expressed in other family members. However, how sibling and parental behavioural strategies are modified by genes expressed in family members, and to what degree, remains unclear. To answer this question, we have used a split litter design in an experimental population of genetically variable mouse families, and identified loci that indirectly affected sibling and maternal behaviour simultaneously. These loci map to genomic regions that also show a direct effect on offspring behaviour. Directly and indirectly affected traits were significantly correlated at the phenotypic level, illustrating how indirect effects are caused. Genetic variants in offspring that influence solicitation also impacted their siblings' and maternal behaviour. However, in contrast to predictions from sibling competition, unrelated litter mates benefited from increased solicitation. Overall, such indirect genetic effects explained a large proportion of variation seen in behaviours, with candidate genes involved in metabolism to neuronal development. These results reveal that we need to view behavioural strategies as the result of conjoint selection on genetic variation in all interacting family members.

Genetic structure of phenotypic robustness in the collaborative cross mouse diallel panel

Developmental stability and canalization describe the ability of developmental systems to minimize phenotypic variation in the face of stochastic micro-environmental effects, genetic variation and environmental influences. Canalization is the ability to minimize the effects of genetic or environmental effects, whereas developmental stability is the ability to minimize the effects of micro-environmental effects within individuals. Despite much attention, the mechanisms that underlie these two components of phenotypic robustness remain unknown. We investigated the genetic structure of phenotypic robustness in the collaborative cross (CC) mouse reference population. We analysed the magnitude of fluctuating asymmetry (FA) and among-individual variation of cranial shape in reciprocal crosses among the eight parental strains, using geometric morphometrics and a diallel analysis based on a Bayesian approach. Significant differences among genotypes were found for both measures, although they were poorly correlated at the level of individuals. An overall positive effect of inbreeding was found for both components of variation. The strain CAST/EiJ exerted a positive additive effect on FA and, to a lesser extent, among-individual variance. Sex- and other strain-specific effects were not significant. Neither FA nor among-individual variation was associated with phenotypic extremeness. Our results support the existence of genetic variation for both developmental stability and canalization. This finding is important because robustness is a key feature of developmental systems. Our finding that robustness is not related to phenotypic extremeness is consistent with theoretical work that suggests that its relationship to stabilizing selection is not straightforward.

Epistasis and maternal effects in experimental adaptation to chronic nutritional stress in Drosophila

Based on ecological and metabolic arguments, some authors predict that adaptation to novel, harsh environments should involve alleles showing negative (diminishing return) epistasis and/or that it should be mediated in part by evolution of maternal effects. Although the first prediction has been supported in microbes, there has been little experimental support for either prediction in multicellular eukaryotes. Here we use a line-cross design to study the genetic architecture of adaptation to chronic larval malnutrition in a population of Drosophila melanogaster that evolved on an extremely nutrient-poor larval food for 84 generations. We assayed three fitness-related traits (developmental rate, adult female weight and egg-to-adult viability) under the malnutrition conditions in 14 crosses between this selected population and a nonadapted control population originally derived from the same base population. All traits showed a pattern of negative epistasis between alleles improving performance under malnutrition. Furthermore, evolutionary changes in maternal traits accounted for half of the 68% increase in viability and for the whole of 8% reduction in adult female body weight in the selected population (relative to unselected controls). These results thus support both of the above predictions and point to the importance of nonadditive effects in adaptive microevolution.

Genomic imprinting and parent-of-origin effects on complex traits

Parent-of-origin effects occur when the phenotypic effect of an allele depends on whether it is inherited from the mother or the father. Several phenomena can cause parent-of-origin effects, but the best characterized is parent-of-origin-dependent gene expression associated with genomic imprinting. The development of new mapping approaches applied to the growing abundance of genomic data has demonstrated that imprinted genes can be important contributors to complex trait variation. Therefore, to understand the genetic architecture and evolution of complex traits, including complex diseases and traits of agricultural importance, it is crucial to account for these parent-of-origin effects. Here, we discuss patterns of phenotypic variation associated with imprinting, evidence supporting its role in complex trait variation and approaches for identifying its molecular signatures.

Family-based study of HTR2A in suicide attempts: observed gene, gene × environment and parent-of-origin associations

While suicidal behavior is frequently accompanied by serotonergic system alterations, specific associations with genetic variation in the serotonin 2A receptor (HTR2A) gene have been inconsistent. Using a family-based study design of 660 offspring who have made a suicide attempt (SA) and both parents, we conducted an association and linkage analysis using single-nucleotide polymorphisms (SNPs) with extensive gene coverage, and included the study of parent-of-origin (POE) and gene-environment interaction (G × E), also using previously unstudied exposures. The main finding was a G × E between the exon 1 SNP rs6313 and exposure to cumulative types of lifetime stressful life events (SLEs), driven by overtransmission of CT and undertransmission of TT, both in relation to other genotypes. Further exploratory analysis revealed a significant POE in this G × E in female subjects, which followed a polar overdominant inheritance pattern. In addition, rs6310 and rs6305 were found to significantly associate with SA in the total sample. A G × E in female subjects (rs7322347 × physical assault in childhood/adolescence) confirmed features of a previously observed association with SA. Other potentially interesting nominally significant findings were observed, but like the G × E of rs7322347 did not pass a false-discovery rate cutoff. Taken together, this study found multiple associations of HTR2A SNPs on SA, with strongest statistical evidence for a G × E involving rs6313, and further suggested the importance of taking into account different inheritance patterns and G × Es with regard to HTR2A.

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 two-step approach to map quantitative trait loci for meat quality in connected porcine F(2) crosses considering main and epistatic effects

The aim of this study was to map QTL for meat quality traits in three connected porcine F(2) crosses comprising around 1000 individuals. The three crosses were derived from the founder breeds Chinese Meishan, European Wild Boar and Pietrain. The animals were genotyped genomewide for approximately 250 genetic markers, mostly microsatellites. They were phenotyped for seven meat quality traits (pH at 45 min and 24 h after slaughter, conductivity at 45 min and 24 h after slaughter, meat colour, drip loss and rigour). QTL mapping was conducted using a two-step procedure. In the first step, the QTL were mapped using a multi-QTL multi-allele model that was tailored to analyse multiple connected F(2) crosses. It considered additive, dominance and imprinting effects. The major gene RYR1:g.1843C>T affecting the meat quality on SSC6 was included as a cofactor in the model. The mapped QTL were tested for pairwise epistatic effects in the second step. All possible epistatic effects between additive, dominant and imprinting effects were considered, leading to nine orthogonal forms of epistasis. Numerous QTL were found. The most interesting chromosome was SSC6. Not all genetic variance of meat quality was explained by RYR1:g.1843C>T. A small confidence interval was obtained, which facilitated the identification of candidate genes underlying the QTL. Epistasis was significant for the pairwise QTL on SSC12 and SSC14 for pH24 and for the QTL on SSC2 and SSC5 for rigour. Some evidence for additional pairwise epistatic effects was found, although not significant. Imprinting was involved in epistasis.

Diet-dependent genetic and genomic imprinting effects on obesity in mice

Although the current obesity epidemic is of environmental origin, there is substantial genetic variation in individual response to an obesogenic environment. In this study, we perform a genome-wide scan for quantitative trait loci (QTLs) affecting obesity per se, or an obese response to a high-fat diet in mice from the LG/J by SM/J Advanced Intercross (AI) Line (Wustl:LG,SM-G16). A total of 1,002 animals from 78 F₁₆ full sibships were weaned at 3 weeks of age and half of each litter placed on high- and low-fat diets. Animals remained on the diet until 20 weeks of age when they were necropsied and the weights of the reproductive, kidney, mesenteric, and inguinal fat depots were recorded. Effects on these phenotypes, along with total fat depot weight and carcass weight at necropsy, were mapped across the genome using 1,402 autosomal single-nucleotide polymorphism (SNP) markers. Haplotypes were reconstructed and additive, dominance, and imprinting genotype scores were derived every 1 cM along the F₁₆ map. Analysis was performed using a mixed model with additive, dominance, and imprinting genotype scores, their interactions with sex, diet, and with sex-by-diet as fixed effects and with family and its interaction with sex, diet, and sex-by-diet as random effects. We discovered 95 trait-specific QTLs mapping to 40 locations. Most QTLs had additive effects with dominance and imprinting effects occurring at two-thirds of the loci. Nearly every locus interacted with sex and/or diet in important ways demonstrating that gene effects are primarily context dependent, changing depending on sex and/or diet.