The evolutionary role of modularity and integration in the hominoid cranium

Patterns of morphological integration and modularity among shape features emerge from genetic and developmental factors with varying pleiotropic effects. Factors or processes affecting morphology only locally may respond to selection more easily than common factors that may lead to deleterious side effects and hence are expected to be more conserved. We briefly review evidence for such global factors in primate cranial development as well as for local factors constrained to either the face or the neurocranium. In a sample comprising 157 crania of Homo sapiens, Pan troglodytes, and Gorilla gorilla, we statistically estimated common and local factors of shape variation from Procrustes coordinates of 347 landmarks and semilandmarks. Common factors with pleiotropic effects on both the face and the neurocranium account for a large amount of shape variation, but mainly by extension or truncation of otherwise conserved developmental pathways. Local factors (modular shape characteristics) have more degrees of freedom for evolutionary change than mere ontogenetic scaling. Cranial shape is similarly integrated during development in all three species, but human evolution involves dissociation among several characteristics. The dissociation has probably been achieved by evolutionary alterations and by the novel emergence of local factors affecting characteristics that are controlled at the same time by the common factors.     

Distinct genetic architectures for male and female inflorescence traits of maize

We compared the genetic architecture of thirteen maize morphological traits in a large population of recombinant inbred lines. Four traits from the male inflorescence (tassel) and three traits from the female inflorescence (ear) were measured and studied using linkage and genome-wide association analyses and compared to three flowering and three leaf traits previously studied in the same population. Inflorescence loci have larger effects than flowering and leaf loci, and ear effects are larger than tassel effects. Ear trait models also have lower predictive ability than tassel, flowering, or leaf trait models. Pleiotropic loci were identified that control elongation of ear and tassel, consistent with their common developmental origin. For these pleiotropic loci, the ear effects are larger than tassel effects even though the same causal polymorphisms are likely involved. This implies that the observed differences in genetic architecture are not due to distinct features of the underlying polymorphisms. Our results support the hypothesis that genetic architecture is a function of trait stability over evolutionary time, since the traits that changed most during the relatively recent domestication of maize have the largest effects.     

Functional genetics of intraspecific ecological interactions in Arabidopsis thaliana

Studying the genetic basis of traits involved in ecological interactions is a fundamental part of elucidating the connections between evolutionary and ecological processes. Such knowledge allows one to link genetic models of trait evolution with ecological models describing interactions within and between species. Previous work has shown that connections between genetic and ecological processes in Arabidopsis thaliana may be mediated by the fact that quantitative trait loci (QTL) with 'direct' effects on traits of individuals also have pleiotropic 'indirect' effects on traits expressed in neighbouring plants. Here, we further explore these connections by examining functional relationships between traits affected directly and indirectly by the same QTL. We develop a novel approach using structural equation models (SEMs) to determine whether observed pleiotropic effects result from traits directly affected by the QTL in focal individuals causing the changes in the neighbours' phenotypes. This hypothesis was assessed using SEMs to test whether focal plant phenotypes appear to mediate the connection between the focal plants' genotypes and the phenotypes of their neighbours, or alternatively, whether the connection between the focal plants' genotypes and the neighbours' phenotypes is mediated by unmeasured traits. We implement this analysis using a QTL of major effect that maps to the well-characterized flowering locus, FRIGIDA. The SEMs support the hypothesis that the pleiotropic indirect effects of this locus arise from size and developmental timing-related traits in focal plants affecting the expression of developmental traits in their neighbours. Our findings provide empirical insights into the genetics and nature of intraspecific ecological interactions. Our technique holds promise in directing future work into the genetic basis and functional relationship of traits mediating and responding to ecological interactions.     

Evolutionary integration of the frog cranium

Evolutionary integration (covariation) of traits has long fascinated biologists because of its potential to elucidate factors that have shaped morphological evolution. Studies of tetrapod crania have identified patterns of evolutionary integration that reflect functional or developmental interactions among traits, but no studies to date have sampled widely across the species-rich lissamphibian order Anura (frogs). Frogs exhibit a vast range of cranial morphologies, life history strategies, and ecologies. Here, using high-density morphometrics we capture cranial morphology for 172 anuran species, sampling every extant family. We quantify the pattern of evolutionary modularity in the frog skull and compare patterns in taxa with different life history modes. Evolutionary changes across the anuran cranium are highly modular, with a well-integrated “suspensorium” involved in feeding. This pattern is strikingly similar to that identified for caecilian and salamander crania, suggesting replication of patterns of evolutionary integration across Lissamphibia. Surprisingly, possession of a feeding larval stage has no notable influence on cranial integration across frogs. However, late-ossifying bones exhibit higher integration than early-ossifying bones. Finally, anuran cranial modules show diverse morphological disparities, supporting the hypothesis that modular variation allows mosaic evolution of the cranium, but we find no consistent relationship between degree of within-module integration and disparity.     

DEVELOPMENTAL QUANTITATIVE GENETICS AND THE EVOLUTION OF ONTOGENIES

A model is presented which permits integration of developmental information into genetic discussions about evolutionary change in morphology. Development of a trait is described in terms of an ontogenetic trajectory whose properties are defined by a small number of parameters. Some evolutionary aspects of development are examined from the perspective of this quantitative genetic model. Particular attention is given to the developmental origin of pleiotropic effects, developmental constraints, heterochrony, and the growth and morphogenesis of complex morphologies. The role of genetic maternal effects in mammalian development is briefly examined, particularly as it relates to selection on developmental traits.     

Juvenile hormone mediates developmental integration between exaggerated traits and supportive traits in the horned flour beetle Gnatocerus cornutus

Sexually selected exaggerated traits are often coupled with modifications in other nontarget traits. In insects with weapons, enlargements of nontarget characters that functionally support the weapon often occur (i.e. supportive traits). The support of sexual traits requires developmental coordination among functionally related multiple traits—an explicit example of morphological integration. The genetic theory predicts that developmental integration among different body modules, for which development is regulated via different sets of genes, is likely to be coordinated by pleiotropic factors. However, the developmental backgrounds of morphological integrations are largely unknown. We tested the hypothesis that the juvenile hormone (JH), as a pleiotropic factor, mediates the integration between exaggerated and supportive traits in an armed beetle Gnatocerus cornutus. During combat, males of this beetle use exaggerated mandibles to lift up their opponents with the supportive traits, that is, the head and prothoracic body parts. Application of methoprene, a JH analog (JHA), during the larval to prepupal period, induced the formation of large mandibles relative to the body sizes in males. Morphometric examination of nontarget traits elucidated an increase in the relative sizes of supportive traits, including the head and prothoracic body parts. In addition, reductions in the hind wing area and elytra length, which correspond to flight and reproductive abilities, were detected. Our findings are consistent with the genetic theory and support the idea that JH is a key pleiotropic factor that coordinates the developmental integration of exaggerated traits and supportive characters, as well as resource allocation trade‐offs.     

Phenotypic Integration Without Modularity: Testing Hypotheses About the Distribution of Pleiotropic Quantitative Trait Loci in a Continuous Space

Theories of phenotypic integration have relied heavily on the concept of modularity in order to model the ways in which traits in an organism correlate and covary. Recent investigations suggest that, while some functional and developmental processes may be morphologically and ontogenetically localized, and thus modular in a developmental sense, there is a great deal of overlap among these influences on patterns of integration in the adult form. This can result in blurry boundaries between hypothesized modules constructed to test hypotheses about phenotypic integration. This investigation tests hypotheses about the contribution of pleiotropic quantitative trait loci (QTL) to phenotypic integration in the mouse mandible without using a priori categorical hypotheses about which traits constitute a module. We ask two main questions: (1) Are the effects of pleiotropic QTL localized to highly correlated traits or more spread out among traits than one might expect by chance? (2) Does the pattern of trait influence when all pleiotropic QTL are considered together deviate from what we might expect if QTL affect traits without regard for the correlations among traits? We find that a large proportion of pleiotropic QTL affect traits that are more highly correlated than we expect by chance with the remainder having effects that are distributed as if by chance. Furthermore, the overall distribution of the effects of pleiotropic QTL differs significantly from the null distribution of no association between pleiotropic effects on traits and correlations among traits. The main modular hypothesis used by earlier studies often does not predict the distribution of sets of traits sharing a common QTL. These results suggest that there is a clear tendency for pleiotropic effects of QTL to be localized but that the localization may be best thought of as occurring in a continuous space rather being clustered in discrete modules.     

The "eyespot module" and eyespots as modules: development, evolution, and integration of a complex phenotype

Organisms are inherently modular, yet modules also evolve in response to selection for functional integration or functional specialization of traits. For serially repeated homologous traits, there is a clear expectation that selection on the function of individual traits will reduce the integration between traits and subdivide a single ancestral module. The eyespots on butterfly wings are one example of serially repeated morphological traits that share a common developmental mechanism but are subject to natural and sexual selection for divergent functions. Here, I test two hypotheses about the organization of the eyespot pattern into independent dorsal-ventral and anterior-posterior modules, using a graphical modeling technique to examine patterns of eyespot covariation among and within wing surfaces in the butterfly Bicyclus anynana. Although there is a hierarchical and complex pattern of integration among eyespots, the results show a surprising mismatch between patterns of eyespot integration and the developmental and evolutionary eyespot units identified in previous empirical studies. These results are discussed in light of the relationships between developmental, functional, and evolutionary modules, and they suggest that developmental sources of independent trait variation are often masked by developmental sources of trait integration.     

Pleiotropic Effects of methoprene-tolerant(Met), a Gene Involved in Juvenile Hormone Metabolism, on Life History Traits in Drosophila melanogaster

Life history theory assumes that there are alleles with pleiotropic effects on fitness components. Although quantitative genetic data are often consistent with pleiotropy, there are few explicit examples of pleiotropic loci. The Drosophila melanogaster gene Methoprene-tolerant (Met) may be such a locus. The Met gene product, a putative juvenile hormone receptor, facilitates the action of juvenile hormone (JH) and JH analogs; JH affects many life history traits in arthropods. Here we use quantitative complementation to investigate effects of Met mutant and wildtype alleles on female developmental time, onset of reproduction, and fecundity. Whereas the alleles did not differ in their effects on developmental time, we detected allelic variation for the onset of reproduction and for age-specific fecundity. Alleles influenced phenotypic co-variances among traits (developmental time and onset of reproduction; onset of reproduction and both early and late fecundity; early and late fecundity), suggesting that alleles of Met vary in their pleiotropic effects upon life history. Furthermore, the genetic covariance between developmental time and early fecundity attributed to alleles of Met was negative, indicating consistent pleiotropic effects among alleles on these traits. The allelic effects of Met support genetic models where pleiotropy at genes associated with hormone regulation can contribute to the evolution of life history traits.     

Indirect genetic effects from ecological interactions in Arabidopsis thaliana

Indirect genetic effects arise when genes expressed in one individual affect the expression of traits in other individuals. The importance of indirect genetic effects has been recognized for a diversity of evolutionary processes including kin selection, sexual selection, community structure and multilevel selection, but data regarding their genetic architecture and prevalence throughout the genome remain scarce, especially for interactions between unrelated individuals. Using a set of 411 Bay-0 x Shahdara Arabidopsis recombinant inbred lines grown with Landsberg neighbours, we examined quantitative trait loci (QTL) having direct and indirect effects on size, developmental, and fitness related traits. Using an interval mapping approach, we identified 15 QTL with direct effects and found that 13 of these QTL had significant indirect effects on trait expression in neighbouring plants. These results suggest widespread pleiotropy, as nearly all direct effect QTL have associated pleiotropic indirect effects. Paradoxically, most indirect effects were of the same sign as direct effects, creating a pattern of nearly universal positive pleiotropy that makes most covariances between direct and indirect effects positive. These results are consistent with a complex genetic basis for intraspecific interactions, but suggest that interactions between neighbouring plants are largely positive, rather than negative as would be expected for competition. In addition to their evolutionary and ecological importance, these pleiotropic relationships between DGE and IGE loci have implications for quantitative genetic studies of natural populations as well as experimental design considerations. Additionally, studies that ignore IGEs may over- or underestimate quantitative genetic parameters, as well as the effect of and variance contributed by QTL.     

Pleiotropic effects on mandibular morphology I. Developmental morphological integration and differential dominance

Pleiotropy refers to a single genetic locus that affects more than one phenotypic trait. Pleiotropic effects of genetic loci are thought to play an important role in evolution, reflecting functional and developmental relationships among phenotypes. In a previous study, we examined pleiotropic effects displayed by quantitative trait loci (QTLs) on murine mandibular morphology in relation to mandibular structure and function. In replicating most of our previous QTLs and increasing our sample size, this study strengthens and extends our earlier results. As in our previous study, we find that QTL effects tend to be restricted to developmentally or functionally related traits. In addition, we examine patterns of differential dominance for pleiotropic QTL effects. Differential dominance occurs when dominance patterns for a single locus vary among traits. We find that multivariate overdominance is a common and substantial phenomenon, and may potentially provide an explanation for the persistence of heterozygosity in natural populations.     

Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies

Genetic correlations between traits may cause correlated responses to selection. Previous models described the conditions under which genetic correlations are expected to be maintained. Selection, mutation, and migration are all proposed to affect genetic correlations, regardless of whether the underlying genetic architecture consists of pleiotropic or tightly linked loci affecting the traits. Here, we investigate the conditions under which pleiotropy and linkage have different effects on the genetic correlations between traits by explicitly modeling multiple genetic architectures to look at the effects of selection strength, degree of correlational selection, mutation rate, mutational variance, recombination rate, and migration rate. We show that at mutation-selection(-migration) balance, mutation rates differentially affect the equilibrium levels of genetic correlation when architectures are composed of pairs of physically linked loci compared to architectures of pleiotropic loci. Even when there is perfect linkage (no recombination within pairs of linked loci), a lower genetic correlation is maintained than with pleiotropy, with a lower mutation rate leading to a larger decrease. These results imply that the detection of causal loci in multitrait association studies will be affected by the type of underlying architectures, whereby pleiotropic variants are more likely to be underlying multiple detected associations. We also confirm that tighter linkage between nonpleiotropic causal loci maintains higher genetic correlations at the traits and leads to a greater proportion of false positives in association analyses.     

Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution

Biologists have long sought to understand which genes and what kinds of changes in their sequences are responsible for the evolution of morphological diversity. Here, I outline eight principles derived from molecular and evolutionary developmental biology and review recent studies of species divergence that have led to a genetic theory of morphological evolution, which states that (1) form evolves largely by altering the expression of functionally conserved proteins, and (2) such changes largely occur through mutations in the cis-regulatory sequences of pleiotropic developmental regulatory loci and of the target genes within the vast networks they control.     

DIFFERENTIAL DOMINANCE OF PLEIOTROPIC LOCI FOR MOUSE SKELETAL TRAITS

The term "differential dominance" describes the situation in which the dominance effects at a pleiotropic locus vary between traits. Directional selection on the phenotype can lead to balancing selection on differentially dominant pleiotropic loci. Even without any individual overdominant traits, some linear combination of traits will display overdominance at a locus displaying differential dominance. Multivariate overdominance may be responsible, in part, for high levels of heterozygosity found in natural populations. We examine differential dominance of 70 mouse skeletal traits at 92 quantitative trait loci (QTL). Our results indicate moderate to strong additive and dominance effects at pleiotropic loci, low levels of individual-trait overdominance, and universal multivariate overdominance. Multivariate overdominance affects a range of 6% to 81% of morphospace, with a mean of 32%. Multivariate overdominance tends to affect a larger percentage of morphospace at pleiotropic loci with antagonistic effects on multiple traits (42%). We conclude that multivariate overdominance is common and should be considered in models and in empirical studies of the role of genetic variation in evolvability.     

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₁₆ 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.     

Evolution of branched regulatory genetic pathways: directional selection on pleiotropic loci accelerates developmental system drift

Developmental systems are regulated by a web of interacting loci. One common and useful approach in studying the evolution of development is to focus on classes of interacting elements within these systems. Here, we use individual-based simulations to study the evolution of traits controlled by branched developmental pathways involving three loci, where one locus regulates two different traits. We examined the system under a variety of selective regimes. In the case where one branch was under stabilizing selection and the other under directional selection, we observed "developmental system drift": the trait under stabilizing selection showed little phenotypic change even though the loci underlying that trait showed considerable evolutionary divergence. This occurs because the pleiotropic locus responds to directional selection and compensatory mutants are then favored in the pathway under stabilizing selection. Though developmental system drift may be caused by other mechanisms, it seems likely that it is accelerated by the same underlying genetic mechanism as that producing the Dobzhansky-Muller incompatibilities that lead to speciation in both linear and branched pathways. We also discuss predictions of our model for developmental system drift and how different selective regimes affect probabilities of speciation in the branched pathway system.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Pleiotropic effects on mandibular morphology II: differential epistasis and genetic variation in morphological integration

The evolution of morphological modularity through the sequestration of pleiotropy to sets of functionally and developmentally related traits requires genetic variation in the relationships between traits. Genetic variation in relationships between traits can result from differential epistasis, where epistatic relationships for pairs of loci are different for different traits. This study maps relationship quantitative trait loci (QTLs), specifically QTLs that affect the relationship between individual mandibular traits and mandible length, across the genome in an F2 intercross of the LG/J and SM/J inbred mouse strains (N = 1045). We discovered 23 relationship QTLs scattered throughout the genome. All mandibular traits were involved in one or more relationship QTL. When multiple traits were affected at a relationship QTL, the traits tended to come from a developmentally restricted region of the mandible, either the muscular processes or the alveolus. About one-third of the relationship QTLs correspond to previously located trait QTLs affecting the same traits. These results comprise examples of genetic variation necessary for an evolutionary response to selection on the range of pleiotropic effects.     

The genetic architecture of domestication in the chicken: effects of pleiotropy and linkage

The extent of pleiotropy and epistasis in quantitative traits remains equivocal. In the case of pleiotropy, multiple quantitative trait loci are often taken to be pleiotropic if their confidence intervals overlap, without formal statistical tests being used to ascertain if these overlapping loci are statistically significantly pleiotropic. Additionally, the degree to which the genetic correlations between phenotypic traits are reflected in these pleiotropic quantitative trait loci is often variable, especially in the case of antagonistic pleiotropy. Similarly, the extent of epistasis in various morphological, behavioural and life-history traits is also debated, with a general problem being the sample sizes required to detect such effects. Domestication involves a large number of trade-offs, which are reflected in numerous behavioural, morphological and life-history traits which have evolved as a consequence of adaptation to selective pressures exerted by humans and captivity. The comparison between wild and domestic animals allows the genetic analysis of the traits that differ between these population types, as well as being a general model of evolution. Using a large F(2) intercross between wild and domesticated chickens, in combination with a dense SNP and microsatellite marker map, both pleiotropy and epistasis were analysed. The majority of traits were found to segregate in 11 tight 'blocks' and reflected the trade-offs associated with domestication. These blocks were shown to have a pleiotropic 'core' surrounded by more loosely linked loci. In contrast, epistatic interactions were almost entirely absent, with only six pairs identified over all traits analysed. These results give insights both into the extent of such blocks in evolution and the development of domestication itself.     

Shared quantitative trait loci underlying the genetic correlation between continuous traits

We review genetic correlations among quantitative traits in light of their underlying quantitative trait loci (QTL). We derive an expectation of genetic correlation from the effects of underlying loci and test whether published genetic correlations can be explained by the QTL underlying the traits. While genetically correlated traits shared more QTL (33%) on average than uncorrelated traits (11%), the actual number of shared QTL shared was small. QTL usually predicted the sign of the correlation with good accuracy, but the quantitative prediction was poor. Approximately 25% of trait pairs in the data set had at least one QTL with antagonistic effects. Yet a significant minority (20%) of such trait pairs have net positive genetic correlations due to such antagonistic QTL 'hidden' within positive genetic correlations. We review the evidence on whether shared QTL represent single pleiotropic loci or closely linked monotropic genes, and argue that strict pleiotropy can be viewed as one end of a continuum of recombination rates where r=0. QTL studies of genetic correlation will likely be insufficient to predict evolutionary trajectories over long time spans in large panmictic populations, but will provide important insights into the trade-offs involved in population and species divergence.