The Role of (Co)variation in Shaping the Response to Selection in New World Leaf-Nosed Bats.

AbstractUnderstanding and predicting the evolutionary responses of complex morphological traits to selection remains a major challenge in evolutionary biology. Because traits are genetically correlated, selection on a particular trait produces both direct effects on the distribution of that trait and indirect effects on other traits in the population. The correlations between traits can strongly impact evolutionary responses to selection and may thus impose constraints on adaptation. Here, we used museum specimens and comparative quantitative genetic approaches to investigate whether the covariation among cranial traits facilitated or constrained the response to selection during the major dietary transitions in one of the world's most ecologically diverse mammalian families-the phyllostomid bats. We reconstructed the set of net selection gradients that would have acted on each cranial trait during the major transitions to feeding specializations and decomposed the selection responses into their direct and indirect components. We found that for all transitions, most traits capturing craniofacial length evolved toward adaptive directions owing to direct selection. Additionally, we showed instances of dietary transitions in which the complex interaction between the patterns of covariation among traits and the strength and direction of selection either constrained or facilitated evolution. Our work highlights the importance of considering the within-species covariation estimates to quantify evolvability and to disentangle the relative contribution of variational constraints versus selective causes for observed patterns.

The developmental-genetics of canalization.

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Variation in mouse pelvic morphology maps to locations enriched in Sox9 Class II and Pitx1 regulatory features.

Variation in pelvic morphology has a complex genetic basis and its patterning and specification is governed by conserved developmental pathways. Whether the mechanisms underlying the differentiation and specification of the pelvis also produce the morphological covariation on which natural selection may act, is still an open question in evolutionary developmental biology. We use high-resolution quantitative trait locus (QTL) mapping in the F34 generation of an advanced intercross experiment (LG,SM-G34 ) to characterize the genetic architecture of the mouse pelvis. We test the prediction that genomic features linked to developmental patterning and differentiation of the hind limb and pelvis and the regulation of chondrogenesis are overrepresented in QTL. We find 31 single QTL trait associations at the genome- or chromosome-wise significance level coalescing to 27 pleiotropic loci. We recover further QTL at a more relaxed significance threshold replicating locations found in a previous experiment in an earlier generation of the same population. QTL were more likely than chance to harbor Pitx1 and Sox9 Class II chromatin immunoprecipitation-seq features active during development of skeletal features. There was weak or no support for the enrichment of seven more categories of developmental features drawn from the literature. Our results suggest that genotypic variation is channeled through a subset of developmental processes involved in the generation of phenotypic variation in the pelvis. This finding indicates that the evolvability of complex traits may be subject to biases not evident from patterns of covariance among morphological features or developmental patterning when either is considered in isolation.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse.

BACKGROUND: Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression. RESULTS: Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape. CONCLUSIONS: Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

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.

A high-fat diet alters genome-wide DNA methylation and gene expression in SM/J mice.

BACKGROUND: While the genetics of obesity has been well defined, the epigenetics of obesity is poorly understood. Here, we used a genome-wide approach to identify genes with differences in both DNA methylation and expression associated with a high-fat diet in mice. RESULTS: We weaned genetically identical Small (SM/J) mice onto a high-fat or low-fat diet and measured their weights weekly, tested their glucose and insulin tolerance, assessed serum biomarkers, and weighed their organs at necropsy. We measured liver gene expression with RNA-seq (using 21 total libraries, each pooled with 2 mice of the same sex and diet) and DNA methylation with MRE-seq and MeDIP-seq (using 8 total libraries, each pooled with 4 mice of the same sex and diet). There were 4356 genes with expression differences associated with diet, with 184 genes exhibiting a sex-by-diet interaction. Dietary fat dysregulated several pathways, including those involved in cytokine-cytokine receptor interaction, chemokine signaling, and oxidative phosphorylation. Over 7000 genes had differentially methylated regions associated with diet, which occurred in regulatory regions more often than expected by chance. Only 5-10% of differentially methylated regions occurred in differentially expressed genes, however this was more often than expected by chance (p = 2.2 × 10- 8). CONCLUSIONS: Discovering the gene expression and methylation changes associated with a high-fat diet can help to identify new targets for epigenetic therapies and inform about the physiological changes in obesity. Here, we identified numerous genes with altered expression and methylation that are promising candidates for further study.

Additive genetic variation in the craniofacial skeleton of baboons (genus Papio) and its relationship to body and cranial size.

OBJECTIVES: Determining the genetic architecture of quantitative traits and genetic correlations among them is important for understanding morphological evolution patterns. We address two questions regarding papionin evolution: (1) what effect do body and cranial size, age, and sex have on phenotypic (VP ) and additive genetic (VA ) variation in baboon crania, and (2) how might additive genetic correlations between craniofacial traits and body mass affect morphological evolution? MATERIALS AND METHODS: We use a large captive pedigreed baboon sample to estimate quantitative genetic parameters for craniofacial dimensions (EIDs). Our models include nested combinations of the covariates listed above. We also simulate the correlated response of a given EID due to selection on body mass alone. RESULTS: Covariates account for 1.2-91% of craniofacial VP . EID VA decreases across models as more covariates are included. The median genetic correlation estimate between each EID and body mass is 0.33. Analysis of the multivariate response to selection reveals that observed patterns of craniofacial variation in extant baboons cannot be attributed solely to correlated response to selection on body mass, particularly in males. DISCUSSION: Because a relatively large proportion of EID VA is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual VP patterns. This may conflate direct selection effects on craniofacial variation with those resulting from a correlated response to body mass selection. This shared genetic variation may partially explain how selection for increased body mass in two different papionin lineages produced remarkably similar craniofacial phenotypes.

Modularity: genes, development and evolution.

Modularity has emerged as a central concept for evolutionary biology, providing the field with a theory of organismal structure and variation. This theory has reframed long standing questions and serves as a unified conceptual framework for genetics, developmental biology and multivariate evolution. Research programs in systems biology and quantitative genetics are bridging the gap between these fields. While this synthesis is ongoing, some major themes have emerged and empirical evidence for modularity has become abundant. In this review, we look at modularity from an historical perspective, highlighting its meaning at different levels of biological organization and the different methods that can be used to detect it. We then explore the relationship between quantitative genetic approaches to modularity and developmental genetic studies. We conclude by investigating the dynamic relationship between modularity and the adaptive landscape and how this potentially shapes evolution and can help bridge the gap between micro- and macroevolution.

Gene expression in human meniscal tears has limited association with early degenerative changes in knee articular cartilage.

Purpose/Aim: Meniscus tears are a common injury to the knee associated with the development of osteoarthritis. Gene expression in the injured meniscus may be associated with early degeneration in the articular cartilage. The purpose of this study was to test the hypothesis that gene expression in meniscus tears is associated with early degenerative changes in the articular cartilage at the time of partial meniscectomy. MATERIALS AND METHODS: Torn meniscus was removed at the time of partial meniscectomy in 68 patients without radiographic osteoarthritis. Meniscal mRNA expression was measured by quantitative PCR for multiple molecular markers of osteoarthritis and cartilage homeostasis. The presence of early degenerative changes in the knee was recorded by X-ray (N = 63), magnetic resonance imaging (MRI, N = 48), and arthroscopy (N = 63). Gene expression was tested for correlation with the presence/absence of degenerative changes after adjusting for age, sex, and body mass index. RESULTS: Overall gene expression varied significantly with degenerative changes based on X-ray (P = 0.047) and MRI (P = 0.018). The linear combination of gene variation was also significant. However, only adiponectin (ADIPOQ) (P = 0.015) was expressed at a significantly lower level in patients with chondrosis on MRI, while the expression of ADIPOQ (P = 0.035) and resistin (RETN) (P = 0.017) was higher in patients with early degenerative changes on X-ray. None of the genes varied significantly with presence/absence of chondrosis as measured by arthroscopy. CONCLUSIONS: There is an overall association of gene expression in meniscal tears to early degenerative changes in the knee, but only a limited number of specific genes demonstrate this relationship. The roles of adiponectin and resistin in knee injury and osteoarthritis deserve further study.

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 F16 generation) and developmental (41 QTL identified in an F34 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.

The effect of dietary fat intake on hepatic gene expression in LG/J AND SM/J mice.

BACKGROUND: The liver plays a major role in regulating metabolic homeostasis and is vital for nutrient metabolism. Identifying the genetic factors regulating these processes could lead to a greater understanding of how liver function responds to a high-fat diet and how that response may influence susceptibilities to obesity and metabolic syndrome. In this study we examine differences in hepatic gene expression between the LG/J and SM/J inbred mouse strains and how gene expression in these strains is affected by high-fat diet. LG/J and SM/J are known to differ in their responses to a high-fat diet for a variety of obesity- and diabetes-related traits, with the SM/J strain exhibiting a stronger phenotypic response to diet. RESULTS: Dietary intake had a significant effect on gene expression in both inbred lines. Genes up-regulated by a high-fat diet were involved in biological processes such as lipid and carbohydrate metabolism; protein and amino acid metabolic processes were down regulated on a high-fat diet. A total of 259 unique transcripts exhibited a significant diet-by-strain interaction. These genes tended to be associated with immune function. In addition, genes involved in biochemical processes related to non-alcoholic fatty liver disease (NAFLD) manifested different responses to diet between the two strains. For most of these genes, SM/J had a stronger response to the high-fat diet than LG/J. CONCLUSIONS: These data show that dietary fat impacts gene expression levels in SM/J relative to LG/J, with SM/J exhibiting a stronger response. This supports previous data showing that SM/J has a stronger phenotypic response to high-fat diet. Based upon these findings, we suggest that SM/J and its cross with the LG/J strain provide a good model for examining non-alcoholic fatty liver disease and its role in metabolic syndrome.

Further evidence for phenotypic signatures of hybridization in descendant baboon populations.

Hybridization may have played a substantial role in shaping the diversity of our evolving lineage. Although recent genomic evidence has shown that hybridization occurred between anatomically modern humans (AMHS) and Neanderthals, it remains difficult to pin down precisely where and when this gene flow took place. Investigations of the hybrid phenotype in primates and other mammals are providing models for identifying signatures of hybridization in the fossil record. However, our understanding of intra- and inter-taxon variation in hybrids is still limited. Moreover, there is little evidence from these studies that is pertinent to the question of how long hybrid skeletal traits persist in descendants, and therefore it is not clear whether observed hybrid phenotypes are evidence of recent (e.g., F1) or much earlier hybridization events. Here, we present an analysis updating a previous study of cranial variation in pedigreed olive and yellow baboons and their hybrids. Results suggest that traits previously associated with hybrids in baboons and other mammalian species are also present in this expanded data set; many of these traits are highly heritable, confirming a genetic basis for their variation in this mixed population. While F1 animals - and especially F1 males - still have the highest number of dental anomalies, these and other atypical traits persist into later hybrid generations (such as F2 and B1). Moreover, non-F1 recombinants also show extremely rare trait variations, including reduced canines and rotated teeth. However, these results must be considered in light of the possibility that some founding individuals may have themselves been unrecognized hybrids. Despite this, the data are compelling, and indicate once again that further controlled research remains to be done on primates and other mammals in order to better understand variation in the hybrid phenotype.

Pleiotropic scaling of gene effects and the 'cost of complexity'.

As perceived by Darwin, evolutionary adaptation by the processes of mutation and selection is difficult to understand for complex features that are the product of numerous traits acting in concert, for example the eye or the apparatus of flight. Typically, mutations simultaneously affect multiple phenotypic characters. This phenomenon is known as pleiotropy. The impact of pleiotropy on evolution has for decades been the subject of formal analysis. Some authors have suggested that pleiotropy can impede evolutionary progress (a so-called 'cost of complexity'). The plausibility of various phenomena attributed to pleiotropy depends on how many traits are affected by each mutation and on our understanding of the correlation between the number of traits affected by each gene substitution and the size of mutational effects on individual traits. Here we show, by studying pleiotropy in mice with the use of quantitative trait loci (QTLs) affecting skeletal characters, that most QTLs affect a relatively small subset of traits and that a substitution at a QTL has an effect on each trait that increases with the total number of traits affected. This suggests that evolution of higher organisms does not suffer a 'cost of complexity' because most mutations affect few traits and the size of the effects does not decrease with pleiotropy.

Genetic effects at pleiotropic loci are context-dependent with consequences for the maintenance of genetic variation in populations.

Context-dependent genetic effects, including genotype-by-environment and genotype-by-sex interactions, are a potential mechanism by which genetic variation of complex traits is maintained in populations. Pleiotropic genetic effects are also thought to play an important role in evolution, reflecting functional and developmental relationships among traits. We examine context-dependent genetic effects at pleiotropic loci associated with normal variation in multiple metabolic syndrome (MetS) components (obesity, dyslipidemia, and diabetes-related traits). MetS prevalence is increasing in Western societies and, while environmental in origin, presents substantial variation in individual response. We identify 23 pleiotropic MetS quantitative trait loci (QTL) in an F(16) advanced intercross between the LG/J and SM/J inbred mouse strains (Wustl:LG,SM-G16; n = 1002). Half of each family was fed a high-fat diet and half fed a low-fat diet; and additive, dominance, and parent-of-origin imprinting genotypic effects were examined in animals partitioned into sex, diet, and sex-by-diet cohorts. We examine the context-dependency of the underlying additive, dominance, and imprinting genetic effects of the traits associated with these pleiotropic QTL. Further, we examine sequence polymorphisms (SNPs) between LG/J and SM/J as well as differential expression of positional candidate genes in these regions. We show that genetic associations are different in different sex, diet, and sex-by-diet settings. We also show that over- or underdominance and ecological cross-over interactions for single phenotypes may not be common, however multidimensional synthetic phenotypes at loci with pleiotropic effects can produce situations that favor the maintenance of genetic variation in populations. Our findings have important implications for evolution and the notion of personalized medicine.

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.

On the relationship between ontogenetic and static allometry.

Ontogenetic and static allometries describe how a character changes in size when the size of the organism changes during ontogeny and among individuals measured at the same developmental stage, respectively. Understanding the relationship between these two types of allometry is crucial to understanding the evolution of allometry and, more generally, the evolution of shape. However, the effects of ontogenetic allometry on static allometry remain largely unexplored. Here, we first show analytically how individual variation in ontogenetic allometry and body size affect static allometry. Using two longitudinal data sets on ontogenetic and static allometry, we then estimate variances and covariances for the different parameters of the ontogenetic allometry defined in our model and assess their relative contribution to the static allometric slope. The mean ontogenetic allometry is the main parameter that determines the static allometric slope, while the covariance between the ontogenetic allometric slope and body size generates most of the discrepancies between ontogenetic and static allometry. These results suggest that the apparent evolutionary stasis of the static allometric slope is not generated by internal (developmental) constraints but more likely results from external constraints imposed by selection.

The road to modularity.

A network of interactions is called modular if it is subdivided into relatively autonomous, internally highly connected components. Modularity has emerged as a rallying point for research in developmental and evolutionary biology (and specifically evo-devo), as well as in molecular systems biology. Here we review the evidence for modularity and models about its origin. Although there is an emerging agreement that organisms have a modular organization, the main open problem is the question of whether modules arise through the action of natural selection or because of biased mutational mechanisms.

Heritability of articular cartilage regeneration and its association with ear wound healing in mice.

OBJECTIVE: Emerging evidence suggests that genetic components contribute significantly to cartilage degeneration in osteoarthritis pathophysiology, but little information is available on the genetics of cartilage regeneration. Therefore, this study was undertaken to investigate cartilage regeneration in genetic murine models using common inbred strains and a set of recombinant inbred (RI) lines generated from LG/J (healer of ear wounds) and SM/J (nonhealer) inbred mouse strains. METHODS: An acute full-thickness cartilage injury was introduced in the trochlear groove of 8-week-old mice (n=265) through microsurgery. Mouse knee joints were sagittally sectioned and stained with toluidine blue to evaluate regeneration. For the ear wound phenotype, a bilateral 2-mm through-and-through puncture was created in 6-week-old mice (n=229), and healing outcomes were measured after 30 days. Broad-sense heritability and genetic correlations were calculated for both phenotypes. RESULTS: Time-course analysis of the RI mouse lines showed no significant regeneration until 16 weeks after surgery; at that time, the strains could be segregated into 3 categories: good, intermediate, and poor healers. Analysis of heritability (H2) showed that both cartilage regeneration (H2=26%; P=0.006) and ear wound closure (H2=53%; P<0.00001) were significantly heritable. The genetic correlations between the two healing phenotypes for common inbred mouse strains (r=0.92) and RI mouse lines (r=0.86) were found to be extremely high. CONCLUSION: Our findings indicate that articular cartilage regeneration in mice is heritable, the differences between the mouse lines are due to genetic differences, and a strong genetic correlation between the two phenotypes exists, indicating that they plausibly share a common genetic basis. We therefore surmise that LG/J by SM/J intercross mice can be used to dissect the genetic basis of variation in cartilage regeneration.

Lack of p21 expression links cell cycle control and appendage regeneration in mice.

Animals capable of regenerating multiple tissue types, organs, and appendages after injury are common yet sporadic and include some sponge, hydra, planarian, and salamander (i.e., newt and axolotl) species, but notably such regenerative capacity is rare in mammals. The adult MRL mouse strain is a rare exception to the rule that mammals do not regenerate appendage tissue. Certain commonalities, such as blastema formation and basement membrane breakdown at the wound site, suggest that MRL mice may share other features with classical regenerators. As reported here, MRL fibroblast-like cells have a distinct cell-cycle (G2/M accumulation) phenotype and a heightened basal and wound site DNA damage/repair response that is also common to classical regenerators and mammalian embryonic stem cells. Additionally, a neutral and alkaline comet assay displayed a persistent level of intrinsic DNA damage in cells derived from the MRL mouse. Similar to mouse ES cells, the p53-target p21 was not expressed in MRL ear fibroblasts. Because the p53/p21 axis plays a central role in the DNA damage response and cell cycle control, we directly tested the hypothesis that p21 down-regulation could functionally induce a regenerative response in an appendage of an otherwise nonregenerating mouse strain. Using the ear hole closure phenotype, a genetically mapped and reliable quantitative indicator of regeneration in the MRL mouse, we show that the unrelated Cdkn1a(tmi/Tyj)/J p21(-/-) mouse (unlike the B6129SF2/J WT control) closes ear holes similar to MRL mice, providing a firm link between cell cycle checkpoint control and tissue regeneration.

The Collaborative Cross, a community resource for the genetic analysis of complex traits.

The goal of the Complex Trait Consortium is to promote the development of resources that can be used to understand, treat and ultimately prevent pervasive human diseases. Existing and proposed mouse resources that are optimized to study the actions of isolated genetic loci on a fixed background are less effective for studying intact polygenic networks and interactions among genes, environments, pathogens and other factors. The Collaborative Cross will provide a common reference panel specifically designed for the integrative analysis of complex systems and will change the way we approach human health and disease.