Exact expressions and numerical evaluation of average evolvability measures for characterizing and comparing [Formula: see text] matrices

Theory predicts that the additive genetic covariance ([Formula: see text]) matrix determines a population's short-term (in)ability to respond to directional selection-evolvability in the Hansen-Houle sense-which is typically quantified and compared via certain scalar indices called evolvability measures. Often, interest is in obtaining the averages of these measures across all possible selection gradients, but explicit formulae for most of these average measures have not been known. Previous authors relied either on approximations by the delta method, whose accuracy is generally unknown, or Monte Carlo evaluations (including the random skewers analysis), which necessarily involve random fluctuations. This study presents new, exact expressions for the average conditional evolvability, average autonomy, average respondability, average flexibility, average response difference, and average response correlation, utilizing their mathematical structures as ratios of quadratic forms. The new expressions are infinite series involving top-order zonal and invariant polynomials of matrix arguments, and can be numerically evaluated as their partial sums with, for some measures, known error bounds. Whenever these partial sums numerically converge within reasonable computational time and memory, they will replace the previous approximate methods. In addition, new expressions are derived for the average measures under a general normal distribution for the selection gradient, extending the applicability of these measures into a substantially broader class of selection regimes.

Morphological integration and evolutionary potential of the primate shoulder: Variation among taxa and implications for genetic covariances with the basicranium, pelvis, and arm

Within the primate order, the morphology of the shoulder girdle is immensely variable and has been shown to reflect the functional demands of the upper limb. The observed morphological variation among extant primate taxa consequently has been hypothesized to be driven by selection for different functional demands. Evolutionary analyses of the shoulder girdle often assess this anatomical region, and its traits, individually, therefore implicitly assuming independent evolution of the shoulder girdle. However, the primate shoulder girdle has developmental and functional covariances with the basicranium and pelvic girdle that have been shown to potentially influence its evolution. It is unknown whether these relationships are similar or even present across primate taxa, and how they may affect morphological variation among primates. This study evaluates the strength of covariance and evolutionary potential across four anatomical regions: shoulder girdle, basicranium, pelvis, and distal humerus. Measures of morphological integration and evolutionary potential (conditioned covariance and evolutionary flexibility) are assessed across eight anthropoid primate taxa. Results demonstrate a consistent pattern of morphological constraint within paired anatomical regions across primates. Differences in evolutionary flexibility are observed among primate genera, with humans having the highest evolutionary potential overall. This pattern does not follow functional differences, but rather a separation between monkeys and apes. Therefore, evolutionary hypotheses of primate shoulder girdle morphological variation that evaluate functional demands alone may not account for the effect of these relationships. Collectively, our findings suggest differences in genetic covariance among anatomical regions may have contributed to the observable morphological variation among taxa.

Evolutionary constraints on physiology confound range shift predictions of two nacellid limpets

Physiological comparisons are fundamental to quantitative assessments of the capacity of species to persist within their current distribution and to predict their rates of redistribution in response to climate change. Yet, the degree to which physiological traits are conserved through evolutionary history may fundamentally constrain the capacity for species to adapt and shift their geographic range. Taxa that straddle major climate transitions provide the opportunity to test the mechanisms underlying evolutionary constraints and how such constraints may influence range shift predictions. Here we focus on two abundant and shallow water nacellid limpets which have representative species on either side of the Polar front. We test the thermal thresholds of the Southern Patagonian limpet, Nacella deaurata and show that its optimal temperatures for growth (4 °C), activity (-1.2 to -0.2 °C) and survival (1 to 8 °C) are mismatched to its currently experienced annual sea surface temperature range (5.9 to 10 °C). Comparisons with the congeneric Antarctic limpet, N. concinna, reveal an evolutionary constraint on N. deaurata physiology, with overlapping thermal capacities, suggesting that a cold climate legacy has been maintained through the evolution of these species. These physiological assessments predict that the South American range of N. deaurata will likely decline with continued warming. It is, however, one of the first species with demonstrated physiological capacity to successfully colonize the cold Southern Ocean. With the expected increase in opportunities for transport within high southern latitudes, N. deaurata has the potential to establish and drive ecological change within the shallow Southern Ocean.

Evolutionary constraints on flicker fusion frequency in Lepidoptera

Flying insects occupy both diurnal and nocturnal niches, and their visual systems encounter distinct challenges in both conditions. Visual adaptations, such as superposition eyes of moths, enhance sensitivity to low light levels but trade off with spatial and temporal resolution. Conversely, apposition eyes of butterflies enable high spatial resolution but are poorly sensitive in dim light. Although diel activity patterns of insects influence visual processing, their role in evolution of visual systems is relatively unexplored. Lepidopteran insects present an excellent system to study how diel activity patterns and phylogenetic position influence the visual transduction system. We addressed this question by comparing electroretinography measurements of temporal response profiles of diverse Lepidoptera to light stimuli that were flickering at different frequencies. Our data show that the eyes of diurnal butterflies are sensitive to visual stimuli of higher temporal frequencies than nocturnal moths. Hesperiid skippers, which are typically diurnal or crepuscular, exhibit intermediate phenotypes with peak sensitivity across broader frequency range. Across all groups, species within families exhibited similar phenotypes irrespective of diel activity. Thus, Lepidopteran photoreceptors may have diversified under phylogenetic constraints, and shifts in their sensitivity to higher temporal frequencies occurred concomitantly with the evolution of diurnal lifestyles.

Evaluating testosterone as a phenotypic integrator: From tissues to individuals to species

Hormones have the potential to bring about rapid phenotypic change; however, they are highly conserved over millions of years of evolution. Here, we examine the evolution of hormone-mediated phenotypes, and the extent to which regulation is achieved via independence or integration of the many components of endocrine systems. We focus on the sex steroid testosterone (T), its cognate receptor (androgen receptor) and related endocrine components. We pose predictions about the mechanisms underlying phenotypic integration, including coordinated sensitivity to T within and among tissues and along the HPG axis. We then assess these predictions with case studies from wild birds, asking whether gene expression related to androgenic signaling naturally co-varies among individuals in ways that would promote phenotypic integration. Finally, we review how mechanisms of integration and independence vary over developmental or evolutionary time, and we find limited support for integration.

Y-chromosomes can constrain adaptive evolution via epistatic interactions with other chromosomes

Background

Variation in the non-coding regions of Y-chromosomes have been shown to influence gene regulation throughout the genome in some systems; a phenomenon termed Y-linked regulatory variation (YRV). This type of sex-specific genetic variance could have important implications for the evolution of male and female traits. If YRV contributes to the additive genetic variation of an autosomally coded trait shared between the sexes (e.g. body size), then selection could facilitate sexually dimorphic evolution via the Y-chromosome. In contrast, if YRV is entirely non-additive (i.e. interacts epistatically with other chromosomes), then Y-chromosomes could constrain trait evolution in both sexes whenever they are selected for the same trait value. The ability for this phenomenon to influence such fundamental evolutionary dynamics remains unexplored.

Results

Here we address the evolutionary contribution of Y-linked variance by selecting for improved male geotaxis in populations possessing multiple Y-chromosomes (i.e. possessed Y-linked additive and/or epistatic variation) or a single Y-chromosome variant (i.e. possessed no Y-linked variation). We found that males from populations possessing Y-linked variation did not significantly respond to selection; however, males from populations with no Y-linked variation did respond. These patterns suggest the presence of a large quantity of Y-linked epistatic variance in the multi-Y population that dramatically slowed its response.

Conclusions

Our results imply that YRV is unlikely to facilitate the evolution of sexually dimorphic traits (at least for the trait examined here), but can interfere with the rate of trait evolution in both males and females. This result could have real biological implications as it suggests that YRV can affect how quickly a population responds to new selective pressures (e.g. invasive species, novel pathogens, or climate change). Considering that YRV influences hundreds of genes and is likely typical of other independently-evolved hemizygous chromosomes, YRV-like phenomena may represent common and significant costs to hemizygous sex determination.

Morphological homeostasis in the fossil record

Morphological homeostasis limits the extent to which genetic and/or environmental variation is translated into phenotypic variation, providing generation-to-generation fitness advantage under a stabilizing selection regime. Depending on its lability, morphological homeostasis might also have a longer-term impact on evolution by restricting the variation-and thus the response to directional selection-of a trait. The fossil record offers an inviting opportunity to investigate whether and how morphological homeostasis constrained trait evolution in lineages or clades on long timescales (thousands to millions of years) that are not accessible to neontological studies. Fossils can also reveal insight into the nature of primitive developmental systems that might not be predictable from the study of modern organisms. The ability to study morphological homeostasis in fossils is strongly limited by taphonomic processes that can destroy, blur, or distort the original biological signal: genetic data are unavailable; phenotypic data can be modified by tectonic or compaction-related deformation; time-averaging limits temporal resolution; and environmental variation is hard to study and impossible to control. As a result of these processes, neither allelic sensitivity (and thus genetic canalization) nor macroenvironmental sensitivity (and thus environmental canalization) can be unambiguously assessed in the fossil record. However, homeorhesis-robustness against microenvironmental variation (developmental noise)-can be assessed in ancient developmental systems by measuring the level of fluctuating asymmetry (FA) in a nominally symmetric trait. This requires the analysis of multiple, minimally time-averaged samples of exquisite preservational quality. Studies of FA in fossils stand to make valuable contributions to our understanding of the deep-time significance of homeorhesis. Few empirical studies have been conducted to date, and future paleontological research focusing on how homeorhesis relates to evolutionary rate (including stasis), species survivorship, and purported macroevolutionary trends in evolvability would reap high reward.

Nuclear genetic codes with a different meaning of the UAG and the UAA codon

BACKGROUND

Departures from the standard genetic code in eukaryotic nuclear genomes are known for only a handful of lineages and only a few genetic code variants seem to exist outside the ciliates, the most creative group in this regard. Most frequent code modifications entail reassignment of the UAG and UAA codons, with evidence for at least 13 independent cases of a coordinated change in the meaning of both codons. However, no change affecting each of the two codons separately has been documented, suggesting the existence of underlying evolutionary or mechanistic constraints.

RESULTS

Here, we present the discovery of two new variants of the nuclear genetic code, in which UAG is translated as an amino acid while UAA is kept as a termination codon (along with UGA). The first variant occurs in an organism noticed in a (meta)transcriptome from the heteropteran Lygus hesperus and demonstrated to be a novel insect-dwelling member of Rhizaria (specifically Sainouroidea). This first documented case of a rhizarian with a non-canonical genetic code employs UAG to encode leucine and represents an unprecedented change among nuclear codon reassignments. The second code variant was found in the recently described anaerobic flagellate Iotanema spirale (Metamonada: Fornicata). Analyses of transcriptomic data revealed that I. spirale uses UAG to encode glutamine, similarly to the most common variant of a non-canonical code known from several unrelated eukaryotic groups, including hexamitin diplomonads (also a lineage of fornicates). However, in these organisms, UAA also encodes glutamine, whereas it is the primary termination codon in I. spirale. Along with phylogenetic evidence for distant relationship of I. spirale and hexamitins, this indicates two independent genetic code changes in fornicates.

CONCLUSIONS

Our study documents, for the first time, that evolutionary changes of the meaning of UAG and UAA codons in nuclear genomes can be decoupled and that the interpretation of the two codons by the cytoplasmic translation apparatus is mechanistically separable. The latter conclusion has interesting implications for possibilities of genetic code engineering in eukaryotes. We also present a newly developed generally applicable phylogeny-informed method for inferring the meaning of reassigned codons.

Measuring Accelerated Rates of Insertions and Deletions Independent of Rates of Nucleotide Substitution

Evolutionary constraint for insertions and deletions (indels) is not necessarily equal to constraint for nucleotide substitutions for any given region of a genome. Knowing the variation in indel-specific evolutionary rates across the sequence will aid our understanding of evolutionary constraints on indels, and help us infer how indels have contributed to the evolution of the sequence. However, unlike for nucleotide substitutions, there has been no phylogenetic method that can statistically infer significantly different rates of indels across the sequence space independent of substitution rates. Here, we have developed a software that will find sites with accelerated evolutionary rates specific to indels, by introducing a scaling parameter that only applies to the indel rates and not to the nucleotide substitution rates. Using the software, we show that we can find regions of accelerated rates of indels in the protein alignments of primate genomes. We also confirm that the sites that have high rates of indels are different from the sites that have high rates of nucleotide substitutions within the protein sequences. By identifying regions with accelerated rates of indels independent of nucleotide substitutions, we will be able to better understand the impact of indel mutations on protein sequence evolution.

Intrasegmental recombination does not contribute to the long-term evolution of group A rotavirus

Rotavirus is a genetically diverse pathogen with an eleven-segmented, double-stranded RNA genome. Intrasegmental recombination has been proposed as a potential mechanism to generate antigenic diversity and a possible route of escape from vaccine-imposed selective pressure. Here intrasegmental recombination was studied by performing a genome-wide scan across the eleven genome segments of 797 publically available rotavirus strains. Sixty-two sequences, or 0.7% of sequences analyzed, have evidence of intrasegmental homologous recombination. None of the specific recombination events is seen in more than one sequence. This uniqueness is consistent with either a spurious finding of recombination or the possibility that recombinant sequences arise naturally but are rapidly purged from the rotavirus population through selection. Arguments for the former explanation are presented. This analysis finds no evidence that intrasegmental recombination leads to ongoing transmission or plays a constructive role in rotavirus evolution. These results have practical implications for phylogenetic analyses and suggest a fundamental constraint that may have shaped rotavirus genome structure and evolution.