A search for quantitative trait loci exhibiting imprinting effects on mouse mandible size and shape

Leveraging Descriptor Learning and Functional Map-based Shape Matching for Automatic Landmark Acquisition

Geometric morphometrics is widely employed across the biological sciences for the quantification of morphological traits. However, the scalability of these methods to large datasets is hampered by the requisite placement of landmarks, which can be laborious and time consuming if done manually. Additionally, the selected landmarks embody a particular hypothesis regarding the critical geometry pertinent to the biological inquiry at hand. Modifying this hypothesis lacks flexibility, necessitating the acquisition of an entirely new set of landmarks on the entire dataset to reflect any theoretical adjustments. In our research, we investigate the precision and accuracy of landmarks derived from the comprehensive set of functional correspondences acquired through the functional map framework of geometry processing. We use a deep functional map network to learn shape descriptors that effectively yield functional map-based and point-to-point correspondences between the specimens in our dataset. We then interrogate these maps to identify corresponding landmarks given manually placed landmarks from the entire dataset. We assess our method by automating the landmarking process on a dataset comprising mandibles from various rodent species, comparing its efficacy against MALPACA, a cutting-edge technique for automatic landmark placement. Compared to MALPACA, our model is notably faster and maintains competitive accuracy. The Root Mean Square Error (RMSE) analysis reveals that while MALPACA generally exhibits the lowest RMSE, our models perform comparably, especially with smaller training datasets, suggesting strong generalizability. Visual evaluations confirm the precision of our landmark placements, with deviations remaining within an acceptable range. These findings underscore the potential of unsupervised learning models in anatomical landmark placement, providing a viable and efficient alternative to traditional methods.

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

A complex genetic architecture underlies mandibular evolution in big mice from Gough Island

Some of the most compelling examples of morphological evolution come from island populations. Alterations in the size and shape of the mandible have been repeatedly observed in murid rodents following island colonization. Despite this pattern and the significance of the mandible for dietary adaptation, the genetic basis of island-mainland divergence in mandibular form remains uninvestigated. To fill this gap, we examined mandibular morphology in 609 F2s from a cross between Gough Island mice, the largest wild house mice on record, and mice from a mainland reference strain (WSB). Univariate genetic mapping identifies 3 quantitative trait loci (QTL) for relative length of the temporalis lever arm and 2 distinct QTL for relative condyle length, 2 traits expected to affect mandibular function that differ between Gough Island mice and WSB mice. Multivariate genetic mapping of coordinates from geometric morphometric analyses identifies 27 QTL contributing to overall mandibular shape. Quantitative trait loci show a complex mixture of modest, additive effects dispersed throughout the mandible, with landmarks including the coronoid process and the base of the ascending ramus frequently modulated by QTL. Additive effects of most shape quantitative trait loci do not align with island-mainland divergence, suggesting that directional selection played a limited role in the evolution of mandibular shape. In contrast, Gough Island mouse alleles at QTL for centroid size and QTL for jaw length increase these measures, suggesting selection led to larger mandibles, perhaps as a correlated response to the evolution of larger bodies.

Does Size Matters? Relationships among Social Dominance and Some Morphometric Traits upon Out-of-Season Reproductive Outcomes in Anestrus Dairy Goats Treated with P4 + eCG

Simple Summary

The possible role of the social rank (R) (i.e., low-LSR, middle-MSR, or high-HSR) in anestrus goats exposed to a P4 + eCG (D) (i.e., 100 or 350 IU) estrus induction protocol (EIP) upon some reproductive response variables was evaluated. Results confirmed that the high social ranked goats amalgamated some fundamental factors to be successful: augmented live weight and corporal measurements, aggressiveness, primacy to food access, and enhanced reproductive outcomes. Such morphometric, behavioral, growth-related, and reproductive advantages shown by the HSR-goats gave evidence to emphasize the need to better comprehend the biological foundation of relevant animal traits, and to be able to define future balanced management and breeding programs. While we still have a fragmentary knowledge regarding the role that social rank, live weight, and morphometric traits play in reproductive success, this study contributes to understanding how social dominance, aligned to morphological and growth related traits, modulates and even determines out-of-season reproductive success.

Abstract

The possible role of the social rank (R) (i.e., low-LSR, middle-MSR, or high-HSR) in anestrus goats exposed to a P4 + eCG (D) (i.e., 100 or 350 IU) estrus induction protocol (EIP) was evaluated. Goats (Alpine-Saanen-Nubian × Criollo; n = 70; 25° North) managed under stall-fed conditions were all ultrasound evaluated to confirm anestrous status, while the social rank was determined 30 d prior to the EIP. The response variables included estrus induction (EI, %), duration of estrus (DUR, h), ovulation rate (OR, n), live weight (LW, kg), thoracic perimeter (TP, cm), thoracic diameter (TD, cm), body length (BL, cm), height at withers (HW, cm), beard length (BEA, cm), compactness index (COM, cm), and anamorphosis index (ANA, cm), as affected by R, D, and the R × D interaction were evaluated, while the correlation coefficients among reproductive and morphometric variables were quantified. An R × D interaction (p < 0.05) affected the response variables EI, DUR, and OR. While the largest (p < 0.05) EI% occurred in the HSR goats, irrespective of eCG (i.e., 100 or 350 IU), both the shortest estrus duration (DUR, h) and the lowest ovulation rate (OR, n) occurred in the LSR + D100 combination, with no differences among HSR and MSR either with D100 or D350. Regarding the LW and morphometric response variables, (i.e., LW, TP, TD, BL, HW, BEA, COM, and ANA) all of them favored either the HSR and MSR groups, with the lowest phenotypic values occurring in the LSR-goats. The EI% was observed to be positively correlated (p < 0.05) with DUR (0.71), LW (0.28), TP (0.31), TD (0.34), BL (0.33), HW (0.35), COM (0.23), and ANA (0.23). While DUR was correlated (p < 0.05) with TP (0.26) and ANA (0.24), OR demonstrated no-correlation (p > 0.05) with any response variable, either reproductive or morphometric. As expected, LW had high correlation coefficients (p < 0.01) with TP (0.86), TD (0.88), BL (0.82), HW (0.75), BEA (0.51), COM (0.97), and ANA (0.75). In general, the morphometric variables as a whole demonstrated important correlation coefficients among them (p < 0.01), ranging from 0.38 up to 0.84. To estimate the importance of the morphometric differences between social rank upon estrus induction, a principal component (PC) analysis was performed based on the correlation matrix derived from the corporal measurements. The PC1 and PC2 explained 70.3% and 17.6% of the morphometric variation, respectively. The PC1 was a measure of the goat size (i.e., small, medium, large) and its association with estrus occurrence was evaluated using a logistic regression model; the bigger the goat, the increased probability of being in estrus, by up to five times compared to small goats. Our results confirm that the higher social ranked, larger goats amalgamated some fundamental factors to be successful: aggressiveness, primacy to food access, augmented live weight, and corporal size; all of these were able to modulate out-of-season reproductive success in crossbred dairy goats subjected to an estrus induction protocol and managed under stall-fed conditions in Northern Mexico.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Conserved but flexible modularity in the zebrafish skull: implications for craniofacial evolvability

Morphological variation is the outward manifestation of development and provides fodder for adaptive evolution. Because of this contingency, evolution is often thought to be biased by developmental processes and functional interactions among structures, which are statistically detectable through forms of covariance among traits. This can take the form of substructures of integrated traits, termed modules, which together comprise patterns of variational modularity. While modularity is essential to an understanding of evolutionary potential, biologists currently have little understanding of its genetic basis and its temporal dynamics over generations. To address these open questions, we compared patterns of craniofacial modularity among laboratory strains, defined mutant lines and a wild population of zebrafish (Danio rerio). Our findings suggest that relatively simple genetic changes can have profound effects on covariance, without greatly affecting craniofacial shape. Moreover, we show that instead of completely deconstructing the covariance structure among sets of traits, mutations cause shifts among seemingly latent patterns of modularity suggesting that the skull may be predisposed towards a limited number of phenotypes. This new insight may serve to greatly increase the evolvability of a population by providing a range of 'preset' patterns of modularity that can appear readily and allow for rapid evolution.

Effects of the Y-chromosome and the dominant hemimelia mutation on the morphology of the mouse mandible

The aims of this study were to test whether the Y-chromosome and the autosomal dominant hemimelia (Dh) mutation can affect mandible morphology in mice. I analyzed mandible size and shape using landmark-based geometric morphometrics in 16 DH-Chr Y^@ -+/+ (@ represents one of the inbred strain names) strains and observed significant differences in mandible size. The largest mandible was identified in strain DH-Chr Y^C3H and the smallest in strain DH-Chr Y^KK . Canonical variate and discriminant function analyses suggested that the mandible shapes of strains DH-Chr Y^C3H and DH-Chr Y^KK differed from those of the other strains. Because seven of the DH-Chr Y^@ -+/+ strains were maintained with dominant hemimelia, I also analyzed the potential influence of dominant hemimelia on mandible morphology because dominant hemimelia is known to cause various skeletal malformations. There were no significant differences in mandible size in seven sets of DH-Chr Y^@ -+/+ and DH-Chr Y^@ -Dh/+ strains. However, canonical variate analysis mapped strains DH-Chr Y^CAS -Dh/+ and DH-Chr Y^CBA -Dh/+ mapped distantly from the rest. Additionally, I observed similar patterns of shape change between DH-Chr Y^CAS -+/+ and DH-Chr Y^CAS -Dh/+, and between DH-Chr Y^CBA -+/+ and DH-Chr Y^CBA -Dh/+. These data indicate that the Y-chromosome affects the size and shape of the mouse mandible. Dominant hemimelia affects mandible shape but not size, and its effects emerge depending on the kinds of Y-chromosomes.

A population level atlas of Mus musculus craniofacial skeleton and automated image-based shape analysis

Laboratory mice are staples for evo/devo and genetics studies. Inbred strains provide a uniform genetic background to manipulate and understand gene-environment interactions, while their crosses have been instrumental in studies of genetic architecture, integration and modularity, and mapping of complex biological traits. Recently, there have been multiple large-scale studies of laboratory mice to further our understanding of the developmental basis, evolution, and genetic control of shape variation in the craniofacial skeleton (i.e. skull and mandible). These experiments typically use micro-computed tomography (micro-CT) to capture the craniofacial phenotype in 3D and rely on manually annotated anatomical landmarks to conduct statistical shape analysis. Although the common choice for imaging modality and phenotyping provides the potential for collaborative research for even larger studies with more statistical power, the investigator (or lab-specific) nature of the data collection hampers these efforts. Investigators are rightly concerned that subtle differences in how anatomical landmarks were recorded will create systematic bias between studies that will eventually influence scientific findings. Even if researchers are willing to repeat landmark annotation on a combined dataset, different lab practices and software choices may create obstacles for standardization beyond the underlying imaging data. Here, we propose a freely available analysis system that could assist in the standardization of micro-CT studies in the mouse. Our proposal uses best practices developed in biomedical imaging and takes advantage of existing open-source software and imaging formats. Our first contribution is the creation of a synthetic template for the adult mouse craniofacial skeleton from 25 inbred strains and five F1 crosses that are widely used in biological research. The template contains a fully segmented cranium, left and right hemi-mandibles, endocranial space, and the first few cervical vertebrae. We have been using this template in our lab to segment and isolate cranial structures in an automated fashion from a mixed population of mice, including craniofacial mutants, aged 4-12.5 weeks. As a secondary contribution, we demonstrate an application of nearly automated shape analysis, using symmetric diffeomorphic image registration. This approach, which we call diGPA, closely approximates the popular generalized Procrustes analysis (GPA) but negates the collection of anatomical landmarks. We achieve our goals by using the open-source advanced normalization tools (ANT) image quantification library, as well as its associated R library (ANTsR) for statistical image analysis. Finally, we make a plea to investigators to commit to using open imaging standards and software in their labs to the extent possible to increase the potential for data exchange and improve the reproducibility of findings. Future work will incorporate more anatomical detail (such as individual cranial bones, turbinals, dentition, middle ear ossicles) and more diversity into the template.

Genetics of Skeletal Evolution in Unusually Large Mice from Gough Island

Organisms on islands often undergo rapid morphological evolution, providing a platform for understanding mechanisms of phenotypic change. Many examples of evolution on islands involve the vertebrate skeleton. Although the genetic basis of skeletal variation has been studied in laboratory strains, especially in the house mouse Mus musculus domesticus, the genetic determinants of skeletal evolution in natural populations remain poorly understood. We used house mice living on the remote Gough Island-the largest wild house mice on record-to understand the genetics of rapid skeletal evolution in nature. Compared to a mainland reference strain from the same subspecies (WSB/EiJ), the skeleton of Gough Island mice is considerably larger, with notable expansions of the pelvis and limbs. The Gough Island mouse skeleton also displays changes in shape, including elongations of the skull and the proximal vs. distal elements in the limbs. Quantitative trait locus (QTL) mapping in a large F₂ intercross between Gough Island mice and WSB/EiJ reveals hundreds of QTL that control skeletal dimensions measured at 5, 10, and/or 16 weeks of age. QTL exhibit modest, mostly additive effects, and Gough Island alleles are associated with larger skeletal size at most QTL. The QTL with the largest effects are found on a few chromosomes and affect suites of skeletal traits. Many of these loci also colocalize with QTL for body weight. The high degree of QTL colocalization is consistent with an important contribution of pleiotropy to skeletal evolution. Our results provide a rare portrait of the genetic basis of skeletal evolution in an island population and position the Gough Island mouse as a model system for understanding mechanisms of rapid evolution in nature.

Does 3D Phenotyping Yield Substantial Insights in the Genetics of the Mouse Mandible Shape?

We describe the application of high-resolution 3D microcomputed tomography, together with 3D landmarks and geometric morphometrics, to validate and further improve previous quantitative genetic studies that reported QTL responsible for variation in the mandible shape of laboratory mice using a new backcross between C57BL/6J and A/J inbred strains. Despite the increasing availability of 3D imaging techniques, artificial flattening of the mandible by 2D imaging techniques seems at first an acceptable compromise for large-scale phenotyping protocols, thanks to an abundance of low-cost digital imaging systems such as microscopes or digital cameras. We evaluated the gain of information from considering explicitly this additional third dimension, and also from capturing variation on the bone surface where no precise anatomical landmark can be marked. Multivariate QTL mapping conducted with different landmark configurations (2D vs. 3D; manual vs. semilandmarks) broadly agreed with the findings of previous studies. Significantly more QTL (23) were identified and more precisely mapped when the mandible shape was captured with a large set of semilandmarks coupled with manual landmarks. It appears that finer phenotypic characterization of the mandibular shape with 3D landmarks, along with higher density genotyping, yields better insights into the genetic architecture of mandibular development. Most of the main variation is, nonetheless, preferentially embedded in the natural 2D plane of the hemi-mandible, reinforcing the results of earlier influential investigations.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.

Genetics of murine craniofacial morphology: diallel analysis of the eight founders of the Collaborative Cross

Using eight inbred founder strains of the mouse Collaborative Cross (CC) project and their reciprocal F1 hybrids, we quantified variation in craniofacial morphology across mouse strains, explored genetic contributions to craniofacial variation that distinguish the founder strains, and tested whether specific or summary measures of craniofacial shape display stronger additive genetic contributions. This study thus provides critical information about phenotypic diversity among CC founder strains and about the genetic contributions to this phenotypic diversity, which is relevant to understanding the basis of variation in standard laboratory strains and natural populations. Craniofacial shape was quantified as a series of size-adjusted linear dimensions (RDs) and by principal components (PC) analysis of morphological landmarks captured from computed tomography images from 62 of the 64 reciprocal crosses of the CC founder strains. We first identified aspects of skull morphology that vary between these phenotypically 'normal' founder strains and that are defining characteristics of these strains. We estimated the contributions of additive and various non-additive genetic factors to phenotypic variation using diallel analyses of a subset of these strongly differing RDs and the first eight PCs of skull shape variation. We find little difference in the genetic contributions to RD measures and PC scores, suggesting fundamental similarities in the magnitude of genetic contributions to both specific and summary measures of craniofacial phenotypes. Our results indicate that there are stronger additive genetic effects associated with defining phenotypic characteristics of specific founder strains, suggesting these distinguishing measures are good candidates for use in genotype-phenotype association studies of CC mice. Our results add significantly to understanding of genotype-phenotype associations in the skull, which serve as a foundation for modeling the origins of medically and evolutionarily relevant variation.

Performance of single and multi-atlas based automated landmarking methods compared to expert annotations in volumetric microCT datasets of mouse mandibles

BACKGROUND

Here we present an application of advanced registration and atlas building framework DRAMMS to the automated annotation of mouse mandibles through a series of tests using single and multi-atlas segmentation paradigms and compare the outcomes to the current gold standard, manual annotation.

RESULTS

Our results showed multi-atlas annotation procedure yields landmark precisions within the human observer error range. The mean shape estimates from gold standard and multi-atlas annotation procedure were statistically indistinguishable for both Euclidean Distance Matrix Analysis (mean form matrix) and Generalized Procrustes Analysis (Goodall F-test). Further research needs to be done to validate the consistency of variance-covariance matrix estimates from both methods with larger sample sizes.

CONCLUSION

Multi-atlas annotation procedure shows promise as a framework to facilitate truly high-throughput phenomic analyses by channeling investigators efforts to annotate only a small portion of their datasets.

Mapping of Craniofacial Traits in Outbred Mice Identifies Major Developmental Genes Involved in Shape Determination

The vertebrate cranium is a prime example of the high evolvability of complex traits. While evidence of genes and developmental pathways underlying craniofacial shape determination is accumulating, we are still far from understanding how such variation at the genetic level is translated into craniofacial shape variation. Here we used 3D geometric morphometrics to map genes involved in shape determination in a population of outbred mice (Carworth Farms White, or CFW). We defined shape traits via principal component analysis of 3D skull and mandible measurements. We mapped genetic loci associated with shape traits at ~80,000 candidate single nucleotide polymorphisms in ~700 male mice. We found that craniofacial shape and size are highly heritable, polygenic traits. Despite the polygenic nature of the traits, we identified 17 loci that explain variation in skull shape, and 8 loci associated with variation in mandible shape. Together, the associated variants account for 11.4% of skull and 4.4% of mandible shape variation, however, the total additive genetic variance associated with phenotypic variation was estimated in ~45%. Candidate genes within the associated loci have known roles in craniofacial development; this includes 6 transcription factors and several regulators of bone developmental pathways. One gene, Mn1, has an unusually large effect on shape variation in our study. A knockout of this gene was previously shown to affect negatively the development of membranous bones of the cranial skeleton, and evolutionary analysis shows that the gene has arisen at the base of the bony vertebrates (Eutelostomi), where the ossified head first appeared. Therefore, Mn1 emerges as a key gene for both skull formation and within-population shape variation. Our study shows that it is possible to identify important developmental genes through genome-wide mapping of high-dimensional shape features in an outbred population.

Formation of the face, mandible, and skull is determined in part by genetic factors, but the relationship between genetic variation and craniofacial development is not well understood. We demonstrate how recent advances in mouse genomics and statistical methods can be used to identify genes involved in craniofacial development. We use outbred mice together with a dense panel of genetic markers to identify genetic loci affecting craniofacial shape. Some of the loci we identify are also known from past studies to contribute to craniofacial development and bone formation. For example, the top candidate gene identified in this study, Mn1, is a gene that appeared at a time when animals started to form bony skulls, suggesting that it may be a key gene in this evolutionary innovation. This further suggests that Mn1 and other genes involved in head formation are also responsible for more fine-grained regulation of its shape. Our results confirm that the outbred mouse population used in this study is suitable to identify single genetic factors even under conditions where many genes cooperate to generate a complex phenotype.

Using whole-genome sequences of the LG/J and SM/J inbred mouse strains to prioritize quantitative trait genes and nucleotides

Background

The laboratory mouse is the most commonly used model for studying variation in complex traits relevant to human disease. Here we present the whole-genome sequences of two inbred strains, LG/J and SM/J, which are frequently used to study variation in complex traits as diverse as aging, bone-growth, adiposity, maternal behavior, and methamphetamine sensitivity.

Results

We identified small nucleotide variants (SNVs) and structural variants (SVs) in the LG/J and SM/J strains relative to the reference genome and discovered novel variants in these two strains by comparing their sequences to other mouse genomes. We find that 39% of the LG/J and SM/J genomes are identical-by-descent (IBD). We characterized amino-acid changing mutations using three algorithms: LRT, PolyPhen-2 and SIFT. We also identified polymorphisms between LG/J and SM/J that fall in regulatory regions and highly informative transcription factor binding sites (TFBS). We intersected these functional predictions with quantitative trait loci (QTL) mapped in advanced intercrosses of these two strains. We find that QTL are both over-represented in non-IBD regions and highly enriched for variants predicted to have a functional impact. Variants in QTL associated with metabolic (231 QTL identified in an F₁₆ generation) and developmental (41 QTL identified in an F₃₄ generation) traits were interrogated and we highlight candidate quantitative trait genes (QTG) and nucleotides (QTN) in a QTL on chr13 associated with variation in basal glucose levels and in a QTL on chr6 associated with variation in tibia length.

Conclusions

We show how integrating genomic sequence with QTL reduces the QTL search space and helps researchers prioritize candidate genes and nucleotides for experimental follow-up. Additionally, given the LG/J and SM/J phylogenetic context among inbred strains, these data contribute important information to the genomic landscape of the laboratory mouse.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1592-3) contains supplementary material, which is available to authorized users.

Quantitative trait loci affecting the 3D skull shape and size in mouse and prioritization of candidate genes in-silico

We describe the first application of high-resolution 3D micro-computed tomography, together with 3D landmarks and geometric morphometrics, to map QTL responsible for variation in skull shape and size using a backcross between C57BL/6J and A/J inbred strains. Using 433 animals, 53 3D landmarks, and 882 SNPs from autosomes, we identified seven QTL responsible for the skull size (SCS.qtl) and 30 QTL responsible for the skull shape (SSH.qtl). Size, sex, and direction-of-cross were all significant factors and included in the analysis as covariates. All autosomes harbored at least one SSH.qtl, sometimes up to three. Effect sizes of SSH.qtl appeared to be small, rarely exceeding 1% of the overall shape variation. However, they account for significant amount of variation in some specific directions of the shape space. Many QTL have stronger effect on the neurocranium than expected from a random vector that will parcellate uniformly across the four cranial regions. On the contrary, most of QTL have an effect on the palate weaker than expected. Combined interval length of 30 SSH.qtl was about 315 MB and contained 2476 known protein coding genes. We used a bioinformatics approach to filter these candidate genes and identified 16 high-priority candidates that are likely to play a role in the craniofacial development and disorders. Thus, coupling the QTL mapping approach in model organisms with candidate gene enrichment approaches appears to be a feasible way to identify high-priority candidates genes related to the structure or tissue of interest.

Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter-subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.

Widespread differential maternal and paternal genome effects on fetal bone phenotype at mid-gestation

Parent-of-origin-dependent (epi)genetic factors are important determinants of prenatal development that program adult phenotype. However, data on magnitude and specificity of maternal and paternal genome effects on fetal bone are lacking. We used an outbred bovine model to dissect and quantify effects of parental genomes, fetal sex, and nongenetic maternal effects on the fetal skeleton and analyzed phenotypic and molecular relationships between fetal muscle and bone. Analysis of 51 bone morphometric and weight parameters from 72 fetuses recovered at day 153 gestation (54% term) identified six principal components (PC1-6) that explained 80% of the variation in skeletal parameters. Parental genomes accounted for most of the variation in bone wet weight (PC1, 72.1%), limb ossification (PC2, 99.8%), flat bone size (PC4, 99.7%), and axial skeletal growth (PC5, 96.9%). Limb length showed lesser effects of parental genomes (PC3, 40.8%) and a significant nongenetic maternal effect (gestational weight gain, 29%). Fetal sex affected bone wet weight (PC1, p < 0.0001) and limb length (PC3, p < 0.05). Partitioning of variation explained by parental genomes revealed strong maternal genome effects on bone wet weight (74.1%, p < 0.0001) and axial skeletal growth (93.5%, p < 0.001), whereas paternal genome controlled limb ossification (95.1%, p < 0.0001). Histomorphometric data revealed strong maternal genome effects on growth plate height (98.6%, p < 0.0001) and trabecular thickness (85.5%, p < 0.0001) in distal femur. Parental genome effects on fetal bone were mirrored by maternal genome effects on fetal serum 25-hydroxyvitamin D (96.9%, p < 0.001) and paternal genome effects on alkaline phosphatase (90.0%, p < 0.001) and their correlations with maternally controlled bone wet weight and paternally controlled limb ossification, respectively. Bone wet weight and flat bone size correlated positively with muscle weight (r = 0.84 and 0.77, p < 0.0001) and negatively with muscle H19 expression (r = -0.34 and -0.31, p < 0.01). Because imprinted maternally expressed H19 regulates growth factors by miRNA interference, this suggests muscle-bone interaction via epigenetic factors.

Lines of least resistance and genetic architecture of house mouse (Mus musculus) mandible shape

Evolution along "lines of least resistance (LLR)" is a well-established phenomenon, which is ultimately based on variational properties of the developmental system. However, it has remained unclear which aspects of development are responsible for observed distributions of variation. This article suggests that house mouse mandible shape may present an opportunity for future research in this field. The genetic architecture of mouse mandible shape has been investigated in the laboratory, and its evolution has been surveyed in studies on natural variation. Here, I ask whether evolutionary diversification of mandible shape follows "LLR" by comparing principal directions of genetic and evolutionary variation, and I assess the potential contribution of specific genomic regions to evolutionary divergence along these directions. The role of two aspects of development, QTL number (presumably reflecting mutational target size) and canalization, is comparatively assessed. Overall, our results demonstrate a strong tendency of genetic and evolutionary systems to vary along similar directions, suggesting that mandible shape evolves along LLR at the level of populations, subspecies and species. At the level of genetic factors, effects bearing similarity to directions of evolution are significantly overrepresented, providing support for the mutational target size hypothesis. A role for canalization is not supported by a clear correlation between directions of evolution and size of genetic effects; however, the evidence for canalization remains ambiguous. These results provide some insights into how the developmental system may shape the variational properties of genetic systems and thus influence the direction of evolution.

Body composition and gene expression QTL mapping in mice reveals imprinting and interaction effects

BACKGROUND

Shifts in body composition, such as accumulation of body fat, can be a symptom of many chronic human diseases; hence, efforts have been made to investigate the genetic mechanisms that underlie body composition. For example, a few quantitative trait loci (QTL) have been discovered using genome-wide association studies, which will eventually lead to the discovery of causal mutations that are associated with tissue traits. Although some body composition QTL have been identified in mice, limited research has been focused on the imprinting and interaction effects that are involved in these traits. Previously, we found that Myostatin genotype, reciprocal cross, and sex interacted with numerous chromosomal regions to affect growth traits.

RESULTS

Here, we report on the identification of muscle, adipose, and morphometric phenotypic QTL (pQTL), translation and transcription QTL (tQTL) and expression QTL (eQTL) by applying a QTL model with additive, dominance, imprinting, and interaction effects. Using an F2 population of 1000 mice derived from the Myostatin-null C57BL/6 and M16i mouse lines, six imprinted pQTL were discovered on chromosomes 6, 9, 10, 11, and 18. We also identified two IGF1 and two Atp2a2 eQTL, which could be important trans-regulatory elements. pQTL, tQTL and eQTL that interacted with Myostatin, reciprocal cross, and sex were detected as well. Combining with the additive and dominance effect, these variants accounted for a large amount of phenotypic variation in this study.

CONCLUSIONS

Our study indicates that both imprinting and interaction effects are important components of the genetic model of body composition traits. Furthermore, the integration of eQTL and traditional QTL mapping may help to explain more phenotypic variation than either alone, thereby uncovering more molecular details of how tissue traits are regulated.

Phenotypic variation across chromosomal hybrid zones of the common shrew (Sorex araneus) indicates reduced gene flow

Sorex araneus, the Common shrew, is a species with more than 70 karyotypic races, many of which form parapatric hybrid zones, making it a model for studying chromosomal speciation. Hybrids between races have reduced fitness, but microsatellite markers have demonstrated considerable gene flow between them, calling into question whether the chromosomal barriers actually do contribute to genetic divergence. We studied phenotypic clines across two hybrid zones with especially complex heterozygotes. Hybrids between the Novosibirsk and Tomsk races produce chains of nine and three chromosomes at meiosis, and hybrids between the Moscow and Seliger races produce chains of eleven. Our goal was to determine whether phenotypes show evidence of reduced gene flow at hybrid zones. We used maximum likelihood to fit tanh cline models to geometric shape data and found that phenotypic clines in skulls and mandibles across these zones had similar centers and widths as chromosomal clines. The amount of phenotypic differentiation across the zones is greater than expected if it were dissipating due to unrestricted gene flow given the amount of time since contact, but it is less than expected to have accumulated from drift during allopatric separation in glacial refugia. Only if heritability is very low, Ne very high, and the time spent in allopatry very short, will the differences we observe be large enough to match the expectation of drift. Our results therefore suggest that phenotypic differentiation has been lost through gene flow since post-glacial secondary contact, but not as quickly as would be expected if there was free gene flow across the hybrid zones. The chromosomal tension zones are confirmed to be partial barriers that prevent differentiated races from becoming phenotypically homogenous.

Mapping shape quantitative trait loci using a radius-centroid-contour model

As the consequence of complex interactions between different parts of an organ, shape can be used as a predictor of structural-functional relationships implicated in changing environments. Despite such importance, however, it is no surprise that little is known about the genetic detail involved in shape variation, because no approach is currently available for mapping quantitative trait loci (QTLs) that control shape. Here, we address this problem by developing a statistical model that integrates the principle of shape analysis into a mixture-model-based likelihood formulated for QTL mapping. One state-of-the-art approach for shape analysis is to identify and analyze the polar coordinates of anatomical landmarks on a shape measured in terms of radii from the centroid to the contour at regular intervals. A procrustes analysis is used to align shapes to filter out position, scale and rotation effects on shape variation. To the end, the accurate and quantitative representation of a shape is produced with aligned radius-centroid-contour (RCC) curves, that is, a function of radial angle at the centroid. The high dimensionality of the RCC data, crucial for a comprehensive description of the geometric feature of a shape, is reduced by principal component (PC) analysis, and the resulting PC axes are treated as phenotypic traits, allowing specific QTLs for global and local shape variability to be mapped, respectively. The usefulness and utilization of the new model for shape mapping in practice are validated by analyzing a mapping data collected from a natural population of poplar, Populus szechuanica var tibetica, and identifying several QTLs for leaf shape in this species. The model provides a powerful tool to compute which genes determine biological shape in plants, animals and humans.

Empirical testing of hypotheses about the evolution of genomic imprinting in mammals

The close interaction between mother and offspring in mammals is thought to contribute to the evolution of genomic imprinting or parent-of-origin dependent gene expression. Empirical tests of theories about the evolution of imprinting have been scant for several reasons. Models make different assumptions about the traits affected by imprinted genes and the scenarios in which imprinting is predicted to have been selected for. Thus, competing hypotheses cannot readily be tested against each other. Further, it is far from clear how predictions about expression patterns of genes with specific phenotypic effects can be tested given current methodology of assaying gene expression levels, be it in the brain or in other tissues. We first set out a scenario for testing competing hypotheses and delineate the different assumptions and predictions of models. We then outline how predictions may be tested using mouse models such as intercrosses or recombinant inbred (RI) systems that can be phenotyped for traits relevant to imprinting theories. Further, we briefly discuss different molecular approaches that may be used in conjunction with experiments to ascertain expression patterns of imprinted genes and thus the testing of predictions.

Exploring the effects of gene dosage on mandible shape in mice as a model for studying the genetic basis of natural variation

Mandible shape in the mouse is a complex trait that is influenced by many genetic factors. However, little is known about the action of single genes on adult mandible shape so far, since most developmentally relevant genes are already required during embryogenesis, i.e., knockouts lead to embryonic death or severe deformations, before the mandible is fully formed. We employ here a geometric morphometric approach to identify subtle phenotypic differences caused by dosage effects of candidate genes. We use mouse strains with specific gene modifications (knockouts and knockins) to compare heterozygous animals with controls from the same stock, which is expected to be equivalent to a change of gene expression of the respective locus. Such differences in expression level are also likely to occur as part of the natural variation. We focus on Bmp pathway genes (Bmp4, its antagonist Noggin, and combinations of Bmp5-7 genotypes), but include also two other developmental control genes suspected to affect mandible development in some way (Egfr and Irf6). In addition, we study the effects of Hoxd13, as well as an extracellular matrix constituent (Col2a1). We find that subtle but significant shape differences are caused by differences in gene dosage of several of these genes. The changes seen for Bmp4 and Noggin are partially compatible with the action of these genes known from birds and fish. We find significant shape changes also for Hoxd13, although this gene has so far only been implicated in skeletal patterning processes of the limbs. Comparing the effect sizes of gene dosage changes to the variation found in natural populations of mice as well as quantitative trait loci (QTL) effects on mandible shape, we find that the effect sizes caused by gene dosage changes are at the lower end of the spectrum of natural variation, but larger than the average additive effects found in QTL studies. We conclude that studying gene dosage effects have the potential to provide new insights into aspects of craniofacial development, variation, and evolution.

Exploration of the genetic organization of morphological modularity on the mouse mandible using a set of interspecific recombinant congenic strains between C57BL/6 and mice of the Mus spretus species

Morphological integration and modularity within semi-autonomous modules are essential mechanisms for the evolution of morphological traits. However, the genetic makeup responsible for the control of variational modularity is still relatively unknown. In our study, we tested the hypothesis that the genetic variation for mandible shape clustered into two morphogenetic components: the alveolar group and the ascending ramus. We used the mouse as a model system to investigate genetics determinants of mandible shape. To do this, we used a combination of geometric morphometric tools and a set of 18 interspecific recombinant congenic strains (IRCS) derived from the distantly related species, Mus spretus SEG/Pas and Mus musculus C57BL/6. Quantitative trait loci (QTL) analysis comparing mandible morphometry between the C57BL/6 and the IRCSs identified 42 putative SEG/Pas segments responsible for the genetic variation. The magnitude of the QTL effects was dependent on the proportion of SEG/Pas genome inherited. Using a multivariate correlation coefficient adapted for modularity assessment and a two-block partial least squares analysis to explore the morphological integration, we found that these QTL clustered into two well-integrated morphogenetic groups, corresponding to the ascending ramus and the alveolar region. Together, these results provide evidence that the mouse mandible is subjected to genetic coordination in a modular manner.

The Measurement of Local Variation in Shape

Geometric morphometrics comprises tools for measuring and analyzing shape as captured by an entire set of landmark configurations. Many interesting questions in evolutionary, genetic, and developmental research, however, are only meaningful at a local level, where a focus on "parts" or "traits" takes priority over properties of wholes. To study variational properties of such traits, current approaches partition configurations into subsets of landmarks which are then studied separately. This approach is unable to fully capture both variational and spatial characteristics of these subsets because interpretability of shape differences is context-dependent. Landmarks omitted from a partition usually contain information about that partition's shape. We present an interpolation-based approach that can be used to model shape differences at a local, infinitesimal level as a function of information available globally. This approach belongs in a large family of methods that see shape differences as continuous "fields" spanning an entire structure, for which landmarks serve as reference parameters rather than as data. We show, via analyses of simulated and real data, how interpolation models provide a more accurate representation of regional shapes than partitioned data. A key difference of this interpolation approach from current morphometric practice is that one must assume an explicit interpolation model, which in turn implies a particular kind of behavior of the regions between landmarks. This choice presents novel methodological challenges, but also an opportunity to incorporate and test biomechanical models that have sought to explain tissue-level processes underlying the generation of morphological shape.

Genome-wide association study of three-dimensional facial morphology identifies a variant in PAX3 associated with nasion position

Craniofacial morphology is highly heritable, but little is known about which genetic variants influence normal facial variation in the general population. We aimed to identify genetic variants associated with normal facial variation in a population-based cohort of 15-year-olds from the Avon Longitudinal Study of Parents and Children. 3D high-resolution images were obtained with two laser scanners, these were merged and aligned, and 22 landmarks were identified and their x, y, and z coordinates used to generate 54 3D distances reflecting facial features. 14 principal components (PCs) were also generated from the landmark locations. We carried out genome-wide association analyses of these distances and PCs in 2,185 adolescents and attempted to replicate any significant associations in a further 1,622 participants. In the discovery analysis no associations were observed with the PCs, but we identified four associations with the distances, and one of these, the association between rs7559271 in PAX3 and the nasion to midendocanthion distance (n-men), was replicated (p = 4 × 10(-7)). In a combined analysis, each G allele of rs7559271 was associated with an increase in n-men distance of 0.39 mm (p = 4 × 10(-16)), explaining 1.3% of the variance. Independent associations were observed in both the z (nasion prominence) and y (nasion height) dimensions (p = 9 × 10(-9) and p = 9 × 10(-10), respectively), suggesting that the locus primarily influences growth in the yz plane. Rare variants in PAX3 are known to cause Waardenburg syndrome, which involves deafness, pigmentary abnormalities, and facial characteristics including a broad nasal bridge. Our findings show that common variants within this gene also influence normal craniofacial development.

A comparative assessment of mandible shape in a consomic strain panel of the house mouse (Mus musculus)--implications for epistasis and evolvability of quantitative traits

Background

Expectations of repeatedly finding associations between given genes and phenotypes have been borne out by studies of parallel evolution, especially for traits involving absence or presence of characters. However, it has rarely been asked whether the genetic basis of quantitative trait variation is conserved at the intra- or even at the interspecific level. This question is especially relevant for shape, where the high dimensionality of variation seems to require a highly complex genetic architecture involving many genes.

Results

We analyse here the genetic effects of chromosome substitution strains carrying M. m. musculus chromosomes in a largely M. m. domesticus background on mandible shape and compare them to the results of previously published QTL mapping data between M. m. domesticus strains. We find that the distribution of genetic effects and effect sizes across the genome is consistent between the studies, while the specific shape changes associated with the chromosomes are different. We find also that the sum of the effects from the different M. m. musculus chromosomes is very different from the shape of the strain from which they were derived, as well as all known wild type shapes.

Conclusions

Our results suggest that the relative chromosome-wide effect sizes are comparable between the long separated subspecies M. m. domesticus and M. m. musculus, hinting at a relative stability of genes involved in this complex trait. However, the absolute effect sizes and the effect directions may be allele-dependent, or are context dependent, i.e. epistatic interactions appear to play an important role in controlling shape.

Micro-evolutionary divergence patterns of mandible shapes in wild house mouse (Mus musculus) populations

BACKGROUND

Insights into the micro-evolutionary patterns of morphological traits require an assessment of the natural variation of the trait within and between populations and closely related species. The mouse mandible is a particularly suitable morphological trait for such an analysis, since it has long been used as a model to study the quantitative genetics of shape. In addition, many distinct populations, sub-species and closely related species are known for the house mouse. However, morphological comparisons among wild caught animals require an assessment in how far environmental and technical factors could interfere with the shape change measurements.

RESULTS

Using geometric morphometrics, we have surveyed mandible shapes in 15 natural populations of the genus Mus, with a focus on the subspecies Mus musculus domesticus. In parallel we have carefully assessed possibly confounding technical and biological factors. We find that there are distinct differences on average between populations, subspecies and species, but these differences are smaller than differences between individuals within populations. Populations from summer-dry regions, although more ancestral, are less distinct from each other than are populations from the more recently colonized northern areas. Populations with especially distinct shapes occur in an area of sympatry of M. m. domesticus and M. spretus and on recently colonized sub-antarctic islands. We have also studied a number of inbred strains to assess in how far their mandible shapes resemble those from the wild. We find that they fall indeed into the shape space of natural variation between individuals in populations.

CONCLUSIONS

Although mandible shapes in natural populations can be influenced by environmental variables, these influences are insufficient to explain the average extent of shape differences between populations, such that evolutionary processes must be invoked to explain this level of diversity. We discuss that adaptive evolution may contribute to shape changes between populations, in particular in newly colonized areas. A comparison between inbred strains and wild mice suggests that the laboratory environment has no major systematic effect on the mandible shape and that such strains can be used as representatives of the natural shape differences between individuals.

A flexible estimating equations approach for mapping function-valued traits

In genetic studies, many interesting traits, including growth curves and skeletal shape, have temporal or spatial structure. They are better treated as curves or function-valued traits. Identification of genetic loci contributing to such traits is facilitated by specialized methods that explicitly address the function-valued nature of the data. Current methods for mapping function-valued traits are mostly likelihood-based, requiring specification of the distribution and error structure. However, such specification is difficult or impractical in many scenarios. We propose a general functional regression approach based on estimating equations that is robust to misspecification of the covariance structure. Estimation is based on a two-step least-squares algorithm, which is fast and applicable even when the number of time points exceeds the number of samples. It is also flexible due to a general linear functional model; changing the number of covariates does not necessitate a new set of formulas and programs. In addition, many meaningful extensions are straightforward. For example, we can accommodate incomplete genotype data, and the algorithm can be trivially parallelized. The framework is an attractive alternative to likelihood-based methods when the covariance structure of the data is not known. It provides a good compromise between model simplicity, statistical efficiency, and computational speed. We illustrate our method and its advantages using circadian mouse behavioral data.

A statistical model for mapping morphological shape

BACKGROUND

Living things come in all shapes and sizes, from bacteria, plants, and animals to humans. Knowledge about the genetic mechanisms for biological shape has far-reaching implications for a range spectrum of scientific disciplines including anthropology, agriculture, developmental biology, evolution and biomedicine.

RESULTS

We derived a statistical model for mapping specific genes or quantitative trait loci (QTLs) that control morphological shape. The model was formulated within the mixture framework, in which different types of shape are thought to result from genotypic discrepancies at a QTL. The EM algorithm was implemented to estimate QTL genotype-specific shapes based on a shape correspondence analysis. Computer simulation was used to investigate the statistical property of the model.

CONCLUSION

By identifying specific QTLs for morphological shape, the model developed will help to ask, disseminate and address many major integrative biological and genetic questions and challenges in the genetic control of biological shape and function.

Paternal genetic contribution influences fetal vulnerability to maternal alcohol consumption in a rat model of fetal alcohol spectrum disorder

BACKGROUND

Fetal alcohol exposure causes in the offspring a collection of permanent physiological and neuropsychological deficits collectively termed Fetal Alcohol Spectrum Disorder (FASD). The timing and amount of exposure cannot fully explain the substantial variability among affected individuals, pointing to genetic influences that mediate fetal vulnerability. However, the aspects of vulnerability that depend on the mother, the father, or both, are not known.

METHODOLOGY/PRINCIPAL FINDINGS

Using the outbred Sprague-Dawley (SD) and inbred Brown Norway (BN) rat strains as well as their reciprocal crosses, we administered ethanol (E), pair-fed (PF), or control (C) diets to the pregnant dams. The dams' plasma levels of free thyroxine (fT4), triiodothyronine (T3), free T3 (fT3), and thyroid stimulating hormone (TSH) were measured to elucidate potential differences in maternal thyroid hormonal environment, which affects specific aspects of FASD. We then compared alcohol-exposed, pair fed, and control offspring of each fetal strain on gestational day 21 (G21) to identify maternal and paternal genetic effects on bodyweight and placental weight of male and female fetuses.

CONCLUSIONS

SD and BN dams exhibited different baseline hypothalamic-pituitary-thyroid function. Moreover, the thyroid function of SD dams was more severely affected by alcohol consumption while that of BN dams was relatively resistant. This novel finding suggests that genetic differences in maternal thyroid function are one source of maternal genetic effects on fetal vulnerability to FASD. The fetal vulnerability to decreased bodyweight after alcohol exposure depended on the genetic contribution of both parents, not only maternal contribution as previously thought. In contrast, the effect of maternal alcohol consumption on placental weight was consistent and not strain-dependent. Interestingly, placental weight in fetuses with different paternal genetic contributions exhibited opposite responses to caloric restriction (pair feeding). In summary, these novel findings demonstrate both maternal and paternal genetic contributions to in utero vulnerability to alcohol, refining our understanding of the genetically-based heterogeneity seen in human FASD.

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

Genomic imprinting, where the effects of alleles depend on their parent-of-origin, can be an important component of the genetic architecture of complex traits. Although there has been a rapidly increasing number of studies of genetic architecture that have examined imprinting effects, none have examined whether imprinting effects depend on genetic background. Such effects are critical for the evolution of genomic imprinting because they allow the imprinting state of a locus to evolve as a function of genetic background. Here we develop a two-locus model of epistasis that includes epistatic interactions involving imprinting effects and apply this model to scan the mouse genome for loci that modulate the imprinting effects of quantitative trait loci (QTL). The inclusion of imprinting leads to nine orthogonal forms of epistasis, five of which do not appear in the usual two-locus decomposition of epistasis. Each form represents a change in the imprinting status of one locus across different classes of genotypes at the other locus. Our genome scan identified two different locus pairs that show complex patterns of epistasis, where the imprinting effect at one locus changes across genetic backgrounds at the other locus. Thus, our model provides a framework for the detection of genetic background-dependent imprinting effects that should provide insights into the background dependence and evolution of genomic imprinting. Our application of the model to a genome scan supports this assertion by identifying pairs of loci that show reciprocal changes in their imprinting status as the background provided by the other locus changes.

Relative contribution of additive, dominance, and imprinting effects to phenotypic variation in body size and growth between divergent selection lines of mice

Epigenetic effects attributed to genomic imprinting are increasingly recognized as an important source of variation in quantitative traits. However, little is known about their relative contribution to phenotypic variation compared to those of additive and dominance effects, and almost nothing about their role in phenotypic evolution. Here we address these questions by investigating the relative contribution of additive, dominance, and imprinting effects of quantitative trait loci (QTL) to variation in "early" and "late" body weight in an intercross of mice selected for divergent adult body weight. We identified 18 loci on 13 chromosomes; additive effects accounted for most of the phenotypic variation throughout development, and imprinting effects were always small. Genetic effects on early weight showed more dominance, less additive, and, surprisingly, less imprinting variation than that of late weight. The predominance of additivity of QTL effects on body weight follows the expectation that additive effects account for the evolutionary divergence between selection lines. We hypothesize that the appearance of more imprinting effects on late body weight may be a consequence of divergent selection on adult body weight, which may have indirectly selected for alleles showing partial imprinting effects due to their associated additive effects, highlighting a potential role of genomic imprinting in the response to selection.