INTERACTING PHENOTYPES AND THE EVOLUTIONARY PROCESS: I. DIRECT AND INDIRECT GENETIC EFFECTS OF SOCIAL INTERACTIONS

Interacting phenotypes are traits whose expression is affected by interactions with conspecifics. Commonly-studied interacting phenotypes include aggression, courtship, and communication. More extreme examples of interacting phenotypes—traits that exist exclusively as a product of interactions—include social dominance, intraspecific competitive ability, and mating systems. We adopt a quantitative genetic approach to assess genetic influences on interacting phenotypes. We partition genetic and environmental effects so that traits in conspecifics that influence the expression of interacting phenotypes are a component of the environment. When the trait having the effect is heritable, the environmental influence arising from the interaction has a genetic basis and can be incorporated as an indirect genetic effect. However, because it has a genetic basis, this environmental component can evolve. Therefore, to consider the evolution of interacting phenotypes we simultaneously consider changes in the direct genetic contributions to a trait (as a standard quantitative genetic approach would evaluate) as well as changes in the environmental (indirect genetic) contribution to the phenotype. We then explore the ramifications of this model of inheritance on the evolution of interacting phenotypes. The relative rate of evolution in interacting phenotypes can be quite different from that predicted by a standard quantitative genetic analysis. Phenotypic evolution is greatly enhanced or inhibited depending on the nature of the direct and indirect genetic effects. Further, unlike most models of phenotypic evolution, a lack of variation in direct genetic effects does not preclude evolution if there is genetic variance in the indirect genetic contributions. The available empirical evidence regarding the evolution of behavior expressed in interactions, although limited, supports the predictions of our model.     

THE LOCI OF REPEATED EVOLUTION: A CATALOG OF GENETIC HOTSPOTS OF PHENOTYPIC VARIATION

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene (“gene reuse”). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution—defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation—are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain‐of‐function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small‐effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large‐effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.     

VARIATION IN GENETIC ARCHITECTURE OF CALLING SONG AMONG POPULATIONS OF ALLONEMOBIUS SOCIUS,A. FASCIATUS,AND A HYBRID POPULATION: DRIFT OR SELECTION?

Predictions using quantitative genetic models generally assume that the variance-covariance matrices remain constant over time. This assumption is based on the supposition that selection is generally weak and hence variation lost through selection can be replaced by new mutations. Whether this is generally true can only be ascertained from empirical studies. Ideally for such a study we should be able to make a prediction concerning the relative strength of selection versus genetic drift. If the latter force is prevalent then the variance-covariances matrices should be proportional to each other. Previous studies have indicated that females in the two sibling cricket species Allonemobius socius and A. fasciatus do not discriminate between males of the two species by their calling song. Therefore, differences between the calling song of the two males most likely result from drift rather than sexual selection. We test this hypothesis by comparing the genetic architecture of calling song of three populations of A. fasciatus with two populations of A. socius. We found no differences among populations within species, but significant differences in the G (genetic) and P (phenotypic) matrices between species, with the matrices being proportional as predicted under the hypothesis of genetic drift. Because of the proportional change in the (co)variances no differences between species are evident in the heritabilities or genetic correlations. Comparison of the two species with a hybrid population from a zone of overlap showed highly significant nonproportional variation in genetic architecture. This variation is consistent with a general mixture of two separate genomes or selection. Qualitative conclusions reached using the phenotypic matrices are the same as those reached using the genetic matrices supporting the hypothesis that the former may be used as surrogate measures of the latter.     

PHYLOGENETIC ANALYSIS OF THE EVOLUTIONARY CORRELATION USING LIKELIHOOD

Many evolutionary processes can lead to a change in the correlation between continuous characters over time or on different branches of a phylogenetic tree. Shifts in genetic or functional constraint, in the selective regime, or in some combination thereof can influence both the evolution of continuous traits and their relation to each other. These changes can often be mapped on a phylogenetic tree to examine their influence on multivariate phenotypic diversification. We propose a new likelihood method to fit multiple evolutionary rate matrices (also called evolutionary variance–covariance matrices) to species data for two or more continuous characters and a phylogeny. The evolutionary rate matrix is a matrix containing the evolutionary rates for individual characters on its diagonal, and the covariances between characters (of which the evolutionary correlations are a function) elsewhere. To illustrate our approach, we apply the method to an empirical dataset consisting of two features of feeding morphology sampled from 28 centrarchid fish species, as well as to data generated via phylogenetic numerical simulations. We find that the method has appropriate type I error, power, and parameter estimation. The approach presented herein is the first to allow for the explicit testing of how and when the evolutionary covariances between characters have changed in the history of a group.     

THE EVOLUTIONARY STABILITY OF CROSS‐SEX,CROSS‐TRAIT GENETIC COVARIANCES

Although knowledge of the selective agents behind the evolution of sexual dimorphism has advanced considerably in recent years, we still lack a clear understanding of the evolutionary durability of cross‐sex genetic covariances that often constrain its evolution. We tested the relative stability of cross‐sex genetic covariances for a suite of homologous contact pheromones of the fruit fly Drosophila serrata, along a latitudinal gradient where these traits have diverged in mean. Using a Bayesian framework, which allowed us to account for uncertainty in all parameter estimates, we compared divergence in the total amount and orientation of genetic variance across populations, finding divergence in orientation but not total variance. We then statistically compared orientation divergence of within‐sex ( G ) to cross‐sex ( B ) covariance matrices. In line with a previous theoretical prediction, we find that the cross‐sex covariance matrix, B , is more variable than either within‐sex G matrix. Decomposition of B matrices into their symmetrical and nonsymmetrical components revealed that instability is linked to the degree of asymmetry. We also find that the degree of asymmetry correlates with latitude suggesting a role for spatially varying natural selection in shaping genetic constraints on the evolution of sexual dimorphism.     

THE QUANTITATIVE AND MOLECULAR GENETIC ARCHITECTURE OF A SUBDIVIDED SPECIES

In an effort to elucidate the evolutionary mechanisms that determine the genetic architecture of a species, we have analyzed 17 populations of the microcrustacean Daphnia pulex for levels of genetic variation at the level of life-history characters and molecular markers in the nuclear and mitochondrial genomes. This species is highly subdivided, with approximately 30% of the variation for nuclear molecular markers and 50% of the variation for mitochondrial markers being distributed among populations. The average level of genetic subdivision for quantitative traits is essentially the same as that for nuclear markers, which superficially suggests that the life-history characters are diverging at the neutral rate. However, the existence of a strong correlation between the levels of population subdivision and broadsense heritabilities of individual traits argues against this interpretation, suggesting instead that the among-population divergence of some quantitative traits (most notably body size) is being driven by local adaptation to different environments. The fact that the mean phenotypes of the individual populations are also strongly correlated with local levels of homozygosity indicates that variation in local inbreeding plays a role in population differentiation. Rather than being a passive consequence of local founder effects, levels of homozygosity may be selected for directly for their effects on the phenotype (adaptive inbreeding depression). There is no relationship between the levels of variation within populations for molecular markers and quantitative characters, and this is explained by the fact that the average standing genetic variation for life-history characters in this species is equivalent to only 33 generations of variation generated by mutation.     

GENETIC COVARIANCE OF FITNESS CORRELATES: WHAT GENETIC CORRELATIONS ARE MADE OF AND WHY IT MATTERS

The genetic variance-covariance matrix, G, is determined in part by functional architecture, the pathways by which variation in genotype influences phenotype. I develop a simple architectural model for G for two traits under directional selection constrained by their dependence on a common limiting resource. I assume that genetic variance is maintained by mutation-selection balance. The relative numbers of loci that play a role in acquiring versus allocating a limiting resource play a crucial role in determining genetic covariance. If many loci are involved in acquiring a resource, genetic covariance may be either negative or positive at equilibrium, depending on the fitness function and the input of mutational variance. The form of G does not necessarily reveal the constraint on resource acquisition inherent in the system, and therefore studies estimating G do not test for the existence of life-history tradeoffs. Characters may evolve in patterns that are unpredictable from G. Experiments are suggested that would indicate if this model could explain observations of positive genetic covariance.     

GENETIC CONSTRAINTS ON MACROEVOLUTION: THE EVOLUTION OF HOST AFFILIATION IN THE LEAF BEETLE GENUS OPHRAELLA

We hypothesize that the evolution of an ecologically important character, the host associations of specialized phytophagous insects, has been influenced by limitations on genetic variation. Using as a historical framework a phylogenetic reconstruction of the history of host associations in the beetle genus Ophraella (Chrysomelidae), we have employed quantitative-genetic methods to screen four species for genetic variation in larval survival, oviposition (in one species only), and feeding responses to their congeners' host plants, in the Asteraceae. We here report results of studies of one species and evaluate the results from all four. Analysis of half-sib/full-sib families and of progenies of wild females of O. notulata, a specialist on Iva (Ambrosiinae), provided evidence of genetic variation in larval consumption of five of six test plants and in adult consumption of four of six. Larval mortality was complete on five plants; only on Ambrosia, a close relative of the natural host, was there appreciable, and genetically variable, survival. Oviposition on Ambrosia showed marginally significant evidence of genetic variation; a more distantly related plant elicited no oviposition at all. In compiling results from four Ophraella species, reported in this and two other papers, we found no evidence of genetic variation in 18 of 39 tests of feeding responses and 14 of 16 tests of larval survival on congeners' hosts. This result is consistent with the hypothesis that absence or paucity of genetic variation may constrain or at least bias the evolution of host associations. The lower incidence of genetic variation in survival than in feeding behavior may imply, according to recent models, that avoidance is a more common evolutionary response to novel plants than adaptation. The usually great disparity between mean performance on congeners' hosts and the species' natural hosts, and an almost complete lack of evidence for negative genetic correlations, argue against the likelihood that speciation has occurred by sympatric host shift. The presence versus apparent absence of genetic variation in consumption was correlated with the propinquity of relationship between the beetle species tested and the species that normally feeds on the test plant, suggesting that the history of host shifts in Ophraella has been guided in part by restrictions on genetic variation. It was also correlated with the propinquity of relationship between a test plant and the beetle's natural host. The contributions of plant relationships and insect relationships, themselves correlated in part, to the pattern of genetic variation, are not readily distinguishable, but together accord with phylogenetic evidence that these and other phytophagous insects adapt most readily to related plants. In this instance, therefore, the macroevolution of an ecologically important character appears to have been influenced by genetic constraints. We hypothesize that absence of the structural prerequisites for genetic variation in complex characters may affect genetic variation and the trajectory of evolution.     

HOMOGENEITY OF THE GENETIC VARIANCE-COVARIANCE MATRIX FOR ANTIPREDATOR TRAITS IN TWO NATURAL POPULATIONS OF THE GARTER SNAKE THAMNOPHIS ORDINOIDES

Quantitative genetic models of evolution rely on the genetic variance-covariance matrix to predict the phenotypic response to selection. Both prospective and retrospective studies of phenotypic evolution across generations rely on assumptions about the constancy of patterns of genetic covariance through time. In the absence of robust theoretical predictions about the stability of genetic covariances, this assumption must be tested with empirical comparisons of genetic parameters among populations and species. Genetic variance-covariance matrices were estimated for a suite of antipredator traits in two populations of the northwestern garter snake, Thamnophis ordinoides. The characters studied include color pattern and antipredator behaviors that interact to facilitate escape from predators. Significant heritabilities for all traits were detected in both populations. Genetic correlations and covariances were found among behaviors in both populations and between color pattern and behavior in one of the populations. Phenotypic means differed among populations, but pairwise comparisons revealed no heterogeneity of genetic parameters between the populations. The structure of the genetic variance-covariance matrix has apparently not changed significantly during the divergence of these two populations.     

THE DUAL ROLE OF SELECTION AND EVOLUTIONARY HISTORY AS REFLECTED IN GENETIC CORRELATIONS

The patterns of genetic correlations between a series of eye and antenna characters were compared among two sets of spring-dwelling and cave-dwelling populations of Gammarus minus. The two sets of populations originate from different drainages and represent two separate invasions of cave habitats from surface-dwelling populations. Matrix correlations, using permutation tests, indicated significant correlations both between populations in the same basin and from the same habitat. The technique of biplot, which allows for the simultaneous consideration of relationships between different genetic correlations and different populations, was used to further analyze the correlation structure. A rank-3 biplot indicated that spring and cave populations were largely differentiated by eye-antennal correlations, whereas basins were differentiated by both eye-antennal and antennal-antennal correlations. Eye-antennal correlations, which are likely to be subject to selection, were most similar within habitats, which are likely to have similar selective regimes.     

QUANTIFYING EVOLUTIONARY POTENTIAL OF MARINE FISH LARVAE: HERITABILITY,SELECTION, AND EVOLUTIONARY CONSTRAINTS

For many marine fish, intense larval mortality may provide considerable opportunity for selection, yet much less is known about the evolutionary potential of larval traits. We combined field demographic studies and manipulative experiments to estimate quantitative genetic parameters for both larval size and swimming performance for a natural population of a common coral‐reef fish, the bicolor damselfish (Stegastes partitus). We also examined selection on larval size by synthesizing information from published estimates of selective mortality. We introduce a method that uses the Lande–Arnold framework for examining selection on quantitative traits to empirically reconstruct adaptive landscapes. This method allows the relationship between phenotypic value and fitness components to be described across a broad range of trait values. Our results suggested that despite strong viability selection for large larvae and moderate heritability (h²= 0.29), evolutionary responses of larvae would likely be balanced by reproductive selection favoring mothers that produce more, smaller offspring. Although long‐term evolutionary responses of larval traits may be constrained by size‐number trade‐offs, our results suggest that phenotypic variation in larval size may be an ecologically important source of variability in population dynamics through effects on larval survival and recruitment to benthic populations.     

ON GENETIC SEGREGATION AND THE EVOLUTION OF SEX

It has recently been argued that because the genetic load borne by an asexual species resulting from segregation, relative to a comparable sexual population, is greater than two, sex can overcome its twofold disadvantage and succeed. We evaluate some of the assumptions underlying this argument and discuss alternative assumptions. Further, we simulate the dynamics of competition between sexual and asexual types. We find that for populations of size 100 and 500 the advantages of segregation do not outweigh the cost of producing males. We conclude that, at least for small populations, drift and the cost of sex govern the evolution of sexuality, not selection or segregation. We believe, however, that if sexual and asexual populations were isolated for a sufficiently long period, segregation might impart a fitness advantage upon sexuals that could compensate for the cost of sex and allow sexuals to outcompete asexuals upon their reunion.     

THE EVOLUTIONARY GENETICS AND DEVELOPMENTAL BASIS OF WING PATTERN VARIATION IN THE BUTTERFLY BICYCLUS ANYNANA

We have studied interactions between developmental processes and genetic variation for the eyespot color pattern on the adult dorsal forewing of the nymphalid butterfly, Bicyclus anynana. Truncation selection was applied in both an upward and a downward direction to the size of a single eyespot consisting of rings with wing scales of differing color pigments. High heritabilities resulted in rapid responses to selection yielding divergent lines with very large or very small eyespots. Strong correlated responses occurred in most of the other eyespots on each wing surface. The cells at the center of a presumptive eyespot (the “focus”) act in the early pupal stage to establish the adult wing pattern. The developmental fate of the scale cells within an eyespot is specified by the “signaling” properties of the focus and the “response” thresholds of the epidermis. The individual eyespots can be envisaged as developmental homologues. Grafting experiments performed with the eyespot foci of the selected lines showed that additive genetic variance exists for both the response and, in particular, the signaling components of the developmental system. The results are discussed in the context of how constraints on the evolution of this wing pattern may be related to the developmental organization.     

THE MEASUREMENT OF SELECTION ON QUANTITATIVE TRAITS: BIASES DUE TO ENVIRONMENTAL COVARIANCES BETWEEN TRAITS AND FITNESS

The use of regression techniques for estimating the direction and magnitude of selection from measurements on phenotypes has become widespread in field studies. A potential problem with these techniques is that environmental correlations between fitness and the traits examined may produce biased estimates of selection gradients. This report demonstrates that the phenotypic covariance between fitness and a trait, used as an estimate of the selection differential in estimating selection gradients, has two components: a component induced by selection itself and a component due to the effect of environmental factors on fitness. The second component is shown to be responsible for biases in estimates of selection gradients. The use of regressions involving genotypic and breeding values instead of phenotypic values can yield estimates of selection gradients that are not biased by environmental covariances. Statistical methods for estimating the coefficients of such regressions, and for testing for biases in regressions involving phenotypic values, are described.     

TRANSLATING BETWEEN MICROEVOLUTIONARY PROCESS AND MACROEVOLUTIONARY PATTERNS: THE CORRELATION STRUCTURE OF INTERSPECIFIC DATA

As species evolve along a phylogenetic tree, we expect closely related species to retain some phenotypic similarities due to their shared evolutionary histories. The amount of expected similarity depends both on the hierarchical phylogenetic structure, and on the specific magnitude and types of evolutionary changes that accumulate during each generation. In this study, we show how models of microevolutionary change can be translated into the resulting macroevolutionary patterns. We illustrate how the structure of phenotypic covariances expected in interspecific measurements can be derived, and how this structure depends on the microevolutionary forces guiding phenotypic change at each generation. We then explore the covariance structure expected from several simple microevolutionary models of phenotypic evolution, including various combinations of random genetic drift, directional selection, stabilizing selection, and environmental change, as well as models of punctuated or burst-like evolution. We find that stabilizing selection leads to patterns of exponential decrease of between species covariance with phylogenetic distance. This is different from the usual linear patterns of decrease assumed in most comparative and systematic methods. Nevertheless, linear patterns of decrease can result from many processes in addition to random genetic drift, such as directional and fluctuating selection as well as modes of punctuated change. Our framework can be used to develop methods for (1) phylogenetic reconstruction; (2) inference of the evolutionary process from comparative data; and (3) conducting or evaluating statistical analyses of comparative data while taking phylogenetic history into account.     

The G matrix under fluctuating correlational mutation and selection

Theoretical quantitative genetics provides a framework for reconstructing past selection and predicting future patterns of phenotypic differentiation. However, the usefulness of the equations of quantitative genetics for evolutionary inference relies on the evolutionary stability of the additive genetic variance-covariance matrix (G matrix). A fruitful new approach for exploring the evolutionary dynamics of G involves the use of individual-based computer simulations. Previous studies have focused on the evolution of the eigenstructure of G. An alternative approach employed in this paper uses the multivariate response-to-selection equation to evaluate the stability of G. In this approach, I measure similarity by the correlation between response-to-selection vectors due to random selection gradients. I analyze the dynamics of G under several conditions of correlational mutation and selection. As found in a previous study, the eigenstructure of G is stabilized by correlational mutation and selection. However, over broad conditions, instability of G did not result in a decreased consistency of the response to selection. I also analyze the stability of G when the correlation coefficients of correlational mutation and selection and the effective population size change through time. To my knowledge, no prior study has used computer simulations to investigate the stability of G when correlational mutation and selection fluctuate. Under these conditions, the eigenstructure of G is unstable under some simulation conditions. Different results are obtained if G matrix stability is assessed by eigenanalysis or by the response to random selection gradients. In this case, the response to selection is most consistent when certain aspects of the eigenstructure of G are least stable and vice versa.     

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Abstract Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with Q _ST generally exceeding G _ST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of Q _ST and G _ST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria , regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex .     

THE EFFECTS OF SELECTION AND GENETIC DRIFT ON THE GENOMIC DISTRIBUTION OF SEXUALLY ANTAGONISTIC ALLELES

Sexual antagonism (SA) occurs when an allele that is beneficial to one sex, is detrimental to the other. This conflict can result in balancing, directional, or disruptive selection acting on SA alleles. A body of theory predicts the conditions under which sexually antagonistic mutants will invade and be maintained in stable polymorphism under balancing selection. There remains, however, considerable debate over the distribution of SA genetic variation across autosomes and sex chromosomes, with contradictory evidence coming from data and theory. In this article, we investigate how the interplay between selection and genetic drift will affect the genomic distribution of sexually antagonistic alleles. The effective population sizes can differ between the autosomes and the sex chromosomes due to a number of ecological factors and, consequently, the distribution of SA genetic variation in genomes. In general, we predict the interplay of SA selection and genetic drift should lead to the accumulation of SA alleles on the X in male heterogametic (XY) species and, on the autosomes in female heterogametic (ZW) species, especially when sexual competition is strong among males.     

MATERNAL INHERITANCE AND ITS EFFECT ON ADAPTIVE EVOLUTION: A QUANTITATIVE GENETIC ANALYSIS OF MATERNAL EFFECTS IN A NATURAL PLANT POPULATION

A mother can influence a trait in her offspring both by the genes she transmits (Mendelian inheritance) and by maternal attributes that directly affect that trait in her offspring (maternal inheritance). Maternal inheritance can alter the direction, rate, and duration of adaptive evolution from standard Mendelian models and its impact on adaptive evolution is virtually unexplored in natural populations. In a hierarchical quantitative genetic analysis to determine the magnitude and structure of maternal inheritance in the winter annual plant, Collinsia verna, I consider three potential models of inheritance. These range from a standard Mendelian model estimating only direct (i.e., Mendelian) additive and environmental variance components to a maternal inheritance model estimating six additive and environmental variance components: direct additive and environmental variances; maternal additive and environmental variances; and the direct-maternal additive () and environmental covariances. The structure of maternal inheritance differs among the 10 traits considered at four stages in the life cycle. Early in the life cycle, seed weight and embryo weight display substantial , a negative , and a positive . Subsequently, cotyledon diameter displays and of roughly the same magnitude and negative . For fall rosettes, leaf number and length are best described by a Mendelian model. In the spring, leaf length displays maternal inheritance with significant and and a negative . All maternally inherited traits show significant negative . Predicted response to selection under maternal inheritance depends on and as well as . Negative results in predicted responses in the opposite direction to selection for seed weight and embryo weight and predicted responses near zero for all subsequent maternally inherited traits. Maternal inheritance persists through the life cycle of this annual plant for a number of size-related traits and will alter the direction and rate of evolutionary response in this population.