Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

The genetic regulatory architecture and epigenomic basis for age-related changes in rattlesnake venom

Developmental phenotypic changes can evolve under selection imposed by age- and size-related ecological differences. Many of these changes occur through programmed alterations to gene expression patterns, but the molecular mechanisms and gene-regulatory networks underlying these adaptive changes remain poorly understood. Many venomous snakes, including the eastern diamondback rattlesnake (Crotalus adamanteus), undergo correlated changes in diet and venom expression as snakes grow larger with age, providing models for identifying mechanisms of timed expression changes that underlie adaptive life history traits. By combining a highly contiguous, chromosome-level genome assembly with measures of expression, chromatin accessibility, and histone modifications, we identified cis-regulatory elements and trans-regulatory factors controlling venom ontogeny in the venom glands of C. adamanteus. Ontogenetic expression changes were significantly correlated with epigenomic changes within genes, immediately adjacent to genes (e.g., promoters), and more distant from genes (e.g., enhancers). We identified 37 candidate transcription factors (TFs), with the vast majority being up-regulated in adults. The ontogenetic change is largely driven by an increase in the expression of TFs associated with growth signaling, transcriptional activation, and circadian rhythm/biological timing systems in adults with corresponding epigenomic changes near the differentially expressed venom genes. However, both expression activation and repression contributed to the composition of both adult and juvenile venoms, demonstrating the complexity and potential evolvability of gene regulation for this trait. Overall, given that age-based trait variation is common across the tree of life, we provide a framework for understanding gene-regulatory-network-driven life-history evolution more broadly.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Cell size, body size and Peto's paradox

Carcinogenesis is one of the leading health concerns afflicting presumably every single animal species, including humans. Currently, cancer research expands considerably beyond medicine, becoming a focus in other branches of natural science. Accumulating evidence suggests that a proportional scale of tumor deaths involves domestic and wild animals and poses economical or conservation threats to many species. Therefore, understanding the genetic and physiological mechanisms of cancer initiation and its progression is essential for our future action and contingent prevention. From this perspective, I used an evolutionary-based approach to re-evaluate the baseline for debate around Peto's paradox. First, I review the background of information on which current understanding of Peto's paradox and evolutionary concept of carcinogenesis have been founded. The weak points and limitations of theoretical modeling or indirect reasoning in studies based on intraspecific, comparative studies of carcinogenesis are highlighted. This is then followed by detail discussion of an effect of the body mass in cancer research and the importance of cell size in consideration of body architecture; also, I note to the ambiguity around cell size invariance hypothesis and hard data for variability of cell size across species are provided. Finally, I point to the new research area that is driving concepts to identify exact molecular mechanisms promoting the process of tumorigenesis, which in turn may provide a proximate explanation of Peto's paradox. The novelty of the approach proposed therein lies in intraspecies testing of the effect of differentiation of cell size/number on the probability of carcinogenesis while controlling for the confounding effect of body mass/size.

Evolved high aerobic capacity has context-specific effects on gut microbiota

Gut microbiota is expected to coevolve with the host's physiology and may play a role in adjusting the host's energy metabolism to suit the host's environment. To evaluate the effects of both evolved host metabolism and the environmental context in shaping the gut microbiota, we used a unique combination of (1) experimental evolution to create selection lines for a fast metabolism and (2) a laboratory-to-field translocation study. Mature bank voles Myodes glareolus from lines selected for high aerobic capacity (A lines) and from unselected control (C lines) were released into large (0.2 ha) outdoor enclosures for longitudinal monitoring. To examine whether the natural environment elicited a similar or more pronounced impact on the gut microbiota of the next generation, we also sampled the field-reared offspring. The gut microbiota were characterized using 16S rRNA amplicon sequencing of fecal samples. The artificial selection for fast metabolism had minimal impact on the gut microbiota in laboratory conditions but in field conditions, there were differences between the selection lines (A lines vs. C lines) in the diversity, community, and resilience of the gut microbiota. Notably, the selection lines differed in the less abundant bacteria throughout the experiment. The lab-to-field transition resulted in an increase in alpha diversity and an altered community composition in the gut microbiota, characterized by a significant increase in the relative abundance of Actinobacteria and a decrease of Patescibacteria. Also, the selection lines showed different temporal patterns in changes in microbiota composition, as the average gut microbiota alpha diversity of the C lines, but not A lines, was temporarily reduced during the initial transition to the field. In surviving young voles, the alpha diversity of gut microbiota was significantly higher in A-line than C-line voles. These results indicate that the association of host metabolism and gut microbiota is context-specific, likely mediated by behavioral or physiological modifications in response to the environment.

Does selection for behavioral and physiological performance traits alter glucocorticoid responsiveness in bank voles?

One of the key elements of an animal's Darwinian fitness is its ability to adequately respond to and cope with challenging situations. Glucocorticoid hormones, such as corticosterone, affect an organism's ability to overcome such challenges. We hypothesized that changes in the glucocorticoid response curve contribute to the evolution of increased performance during challenging conditions, and tested it on bank voles (Myodes glareolus) from a multidirectional artificial selection experiment, which involves lines selected for high aerobic exercise metabolism achieved during swimming (A - Aerobic), predatory behavior towards a cricket (P - Predatory) and ability to maintain body mass on a low-quality herbivorous diet (H - Herbivorous), as well as unselected control lines (C - Control). We elicited a glucocorticoid response either by restraining the animal or by maximum pharmacological stimulation, and measured plasma corticosterone levels at baseline, during the response and during the recovery phase. Response-level corticosterone was higher in females, and recovery from maximal level was faster than that of males. Selection did not affect baseline or stress-induced corticosterone levels, but it decreased the maximum corticosterone level in Aerobic and Predatory lines, reducing the difference between stress-induced and maximum levels. Recovery from restraint-induced corticosterone level tended to be slower in the Herbivorous than in the other lines, an effect that was stronger in females than in males. In conclusion, successful selection for increased performance in challenging conditions was not associated with changes in absolute values of the glucocorticoid response to stress, but can affect other characteristics of the glucocorticoid response curve.

Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication

The olive (Olea europaea L. subsp. europaea) is one of the oldest and most socio-economically important cultivated perennial crop in the Mediterranean region. Yet, its origins are still under debate and the genetic bases of the phenotypic changes associated with its domestication are unknown. We generated RNA-sequencing data for 68 wild and cultivated olive trees to study the genetic diversity and structure both at the transcription and sequence levels. To localize putative genes or expression pathways targeted by artificial selection during domestication, we employed a two-step approach in which we identified differentially expressed genes and screened the transcriptome for signatures of selection. Our analyses support a major domestication event in the eastern part of the Mediterranean basin followed by dispersion towards the West and subsequent admixture with western wild olives. While we found large changes in gene expression when comparing cultivated and wild olives, we found no major signature of selection on coding variants and weak signals primarily affected transcription factors. Our results indicated that the domestication of olives resulted in only moderate genomic consequences and that the domestication syndrome is mainly related to changes in gene expression, consistent with its evolutionary history and life history traits.

The effect of monoamines reuptake inhibitors on aerobic exercise performance in bank voles from a selection experiment

Exercise performance depends on both physiological abilities (e.g., muscle strength) and behavioral characteristics (e.g., motivation). We tested the hypothesis that evolution of increased aerobic exercise performance can be facilitated by evolution of neuropsychological mechanisms responsible for motivation to undertake physical activity. We used a unique model system: lines of bank voles Myodes glareolus selected for high swim-induced aerobic metabolism ("aerobic" A lines). In generation 21, voles from the 4 A lines achieved a 57% higher "voluntary maximum" swim-induced aerobic metabolism (VO2swim) than voles from 4 unselected, "control" C lines. In C lines, VO2swim was 9% lower than the maximum forced-exercise aerobic metabolism (VO2run; P = 0.007), while in A lines it was even higher than VO2run, although not significantly (4%, P = 0.15). Thus, we hypothesized that selection changed both the aerobic capacity and the neuronal mechanisms behind motivation to undertake activity. We investigated the influence of reuptake inhibitors of dopamine (DARI), serotonin (SSRI), and norepinephrine (NERI) on VO2swim. The drugs decreased VO2swim both in C and A lines (% decrease compared with saline: DARI 8%, P < 0.001; SSRI 6%, P < 0.001; NERI 8%, P < 0.001), but the proportional response differed between selection directions only for NERI (stronger effect in C lines: P = 0.008) and the difference was marginally non-significant for SSRI (P = 0.07) and DARI (P = 0.06). Thus, the results suggest that all the 3 monoamines are involved in signaling pathways controlling the motivation to be active and that norepinephrine could have played a role in the evolution of increased aerobic exercise performance in our animal model.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis

Background

Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes.

Results

Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes.

Conclusions

Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.

The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii)

Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers.

Preliminary insights into the genetics of bank vole tolerance to Puumala hantavirus in Sweden

Natural reservoirs of zoonotic pathogens generally seem to be capable of tolerating infections. Tolerance and its underlying mechanisms remain difficult to assess using experiments or wildlife surveys. High-throughput sequencing technologies give the opportunity to investigate the genetic bases of tolerance, and the variability of its mechanisms in natural populations. In particular, population genomics may provide preliminary insights into the genes shaping tolerance and potentially influencing epidemiological dynamics. Here, we addressed these questions in the bank vole Myodes glareolus, the specific asymptomatic reservoir host of Puumala hantavirus (PUUV), which causes nephropathia epidemica (NE) in humans. Despite the continuous spatial distribution of M. glareolus in Sweden, NE is endemic to the northern part of the country. Northern bank vole populations in Sweden might exhibit tolerance strategies as a result of coadaptation with PUUV. This may favor the circulation and maintenance of PUUV and lead to high spatial risk of NE in northern Sweden. We performed a genome-scan study to detect signatures of selection potentially correlated with spatial variations in tolerance to PUUV. We analyzed six bank vole populations from Sweden, sampled from northern NE-endemic to southern NE-free areas. We combined candidate gene analyses (Tlr4, Tlr7, and Mx2 genes) and high-throughput sequencing of restriction site-associated DNA (RAD) markers. Outlier loci showed high levels of genetic differentiation and significant associations with environmental data including variations in the regional number of NE human cases. Among the 108 outliers that matched to mouse protein-coding genes, 14 corresponded to immune-related genes. The main biological pathways found to be significantly enriched corresponded to immune processes and responses to hantavirus, including the regulation of cytokine productions, TLR cascades, and IL-7, VEGF, and JAK-STAT signaling. In the future, genome-scan replicates and functional experimentations should enable to assess the role of these biological pathways in M. glareolus tolerance to PUUV.

Experimental evolution of personality traits: open-field exploration in bank voles from a multidirectional selection experiment

Evolution of complex physiological adaptations could be driven by natural selection acting on behavioral traits. Consequently, animal personality traits and their correlation with physiological traits have become an engaging research area. Here, we applied a unique experimental evolution model-lines of bank voles selected for (A) high exercise-induced aerobic metabolism, (H) ability to cope with low-quality herbivorous diet, and (P) intensity of predatory behavior, that is, traits shaping evolutionary path and diversity of mammals-and asked how the selection affected the voles' personality traits, assessed in an open field test. The A- and P-line voles were more active, whereas the H-line voles were less active, compared those from unselected control lines (C). H-line voles moved slower but on more meandering trajectories, which indicated a more thorough exploration, whereas the A- and P-line voles moved faster and on straighter trajectories. A-line voles showed also an increased escape propensity, whereas P-line voles tended to be bolder. The remarkable correlated responses to the selection indicate a common genetic underlying mechanism of behavioral and physiological traits, and support the paradigm of evolutionary physiology built around the concept of correlated evolution of behavior and physiology.

Quantity, Not Quality: Rapid Adaptation in a Polygenic Trait Proceeded Exclusively through Expression Differentiation

A trait's genomic architecture can affect the rate and mechanism of adaptation, and although many ecologically-important traits are polygenic, most studies connecting genotype, phenotype, and fitness in natural populations have focused on traits with relatively simple genetic bases. To understand the genetic basis of polygenic adaptation, we must integrate genomics, phenotypic data, ecology, and fitness effects for a genetically tractable, polygenic trait; snake venoms provide such a system for studying polygenic adaptation because of their genetic tractability and vital ecological role in feeding and defense. We used a venom transcriptome-proteome map, quantitative proteomics, genomics, and fitness assays in sympatric prey to construct a genotype-phenotype-fitness map for the venoms of an island-mainland pair of rattlesnake populations. Reciprocal fitness experiments demonstrated that each population was locally adapted to sympatric prey. We identified significant expression differentiation with little to no coding-sequence variation across populations, demonstrating that expression differentiation was exclusively the genetic basis of polygenic adaptation. Previous research on the genetics of adaptation, however, has largely been biased toward investigating protein-coding regions because of the complexity of gene regulation. Our results showed that biases at the molecular level can be in the opposite direction, highlighting the need for more systematic comparisons of different molecular mechanisms underlying rapid, adaptive evolution in polygenic traits.

A Shift in the Thermoregulatory Curve as a Result of Selection for High Activity-Related Aerobic Metabolism

According to the “aerobic capacity model,” endothermy in birds and mammals evolved as a result of natural selection favoring increased persistent locomotor activity, fuelled by aerobic metabolism. However, this also increased energy expenditure even during rest, with the lowest metabolic rates occurring in the thermoneutral zone (TNZ) and increasing at ambient temperatures (T_a) below and above this range, depicted by the thermoregulatory curve. In our experimental evolution system, four lines of bank voles (Myodes glareolus) have been selected for high swim-induced aerobic metabolism and four unselected lines have been maintained as a control. In addition to a 50% higher rate of oxygen consumption during swimming, the selected lines have also evolved a 7.3% higher mass-adjusted basal metabolic rate. Therefore, we asked whether voles from selected lines would also display a shift in the thermoregulatory curve and an increased body temperature (T_b) during exposure to high T_a. To test these hypotheses we measured the RMR and T_b of selected and control voles at T_a from 10 to 34°C. As expected, RMR within and around the TNZ was higher in selected lines. Further, the T_b of selected lines within the TNZ was greater than the T_b of control lines, particularly at the maximum measured T_a of 34°C, suggesting that selected voles are more prone to hyperthermia. Interestingly, our results revealed that while the slope of the thermoregulatory curve below the lower critical temperature (LCT) is significantly lower in the selected lines, the LCT (26.1°C) does not differ. Importantly, selected voles also evolved a higher maximum thermogenesis, but thermal conductance did not increase. As a consequence, the minimum tolerated temperature, calculated from an extrapolation of the thermoregulatory curve, is 8.4°C lower in selected (−28.6°C) than in control lines (−20.2°C). Thus, selection for high aerobic exercise performance, even though operating under thermally neutral conditions, has resulted in the evolution of increased cold tolerance, which, under natural conditions, could allow voles to inhabit colder environments. Further, the results of the current experiment support the assumptions of the aerobic capacity model of the evolution of endothermy.

Age-related changes of physiological performance and survivorship of bank voles selected for high aerobic capacity

Variation in lifespans is an intriguing phenomenon, but how metabolic rate influence this variation remains unclear. High aerobic capacity can result in health benefits, but also in increased oxidative damage and accelerated ageing. We tested these contradictory predictions using bank voles (Myodes=Clethrionomys glareolus) from lines selected for high swim-induced aerobic metabolism (A), which had about 50% higher maximum metabolic rate and a higher basal and routine metabolic rates, than those from unselected control lines (C). We measured sprint speed (VSmax), forced-running maximum metabolic rate (VO₂run), maximum long-distance running speed (VLmax), running speed at VO₂run (VVO₂), and respiratory quotient at VO₂run (RQ) at three age classes (I: 3-5, II: 12-14, III: 17-19months), and analysed survivorship. We asked if ageing, understood as the age-related decline of the performance traits, differs between the A and C lines. At age class I, voles from A lines had 19% higher VO₂run, and 12% higher VLmax, but tended to have 19% lower VSmax, than those from C lines. RQ was nearly 1.0 for both A and C lines. The pattern of age-related changes differed between the lines mainly between age classes I and II, but not in older animals. VSmax increased by 27% in A lines and by 10% in C lines between age class I and II, but between classes II and III, it increased by 16% in both selection directions. VO₂run decreased by 7% between age class I and II in A lines only, but in C lines it remained constant across all age classes. VLmax decreased by 8% and VVO₂ by 12% between age classes II and III, but similarly in both selection directions. Mortality was higher in A than in C lines only between the age of 1 and 4months. The only trait for which the changes in old animals differed between the lines was RQ. In A lines, RQ increased between age classes II and III, whereas in C lines such an increase occurred between age classes I and II. Thus, we did not find obvious effects of selection on the pattern of ageing. However, the physiological performance and mortality of bank voles remained surprisingly robust to ageing, at least until the age of 17-19months, similar to the maximum lifespan under natural conditions. Therefore, it is possible that the selection could affect the pattern of ageing in even older individuals when symptoms of senility might be more profound.

Genetic and Genomic Response to Selection for Food Consumption in Drosophila melanogaster

Food consumption is an essential component of animal fitness; however, excessive food intake in humans increases risk for many diseases. The roles of neuroendocrine feedback loops, food sensing modalities, and physiological state in regulating food intake are well understood, but not the genetic basis underlying variation in food consumption. Here, we applied ten generations of artificial selection for high and low food consumption in replicate populations of Drosophila melanogaster. The phenotypic response to selection was highly asymmetric, with significant responses only for increased food consumption and minimal correlated responses in body mass and composition. We assessed the molecular correlates of selection responses by DNA and RNA sequencing of the selection lines. The high and low selection lines had variants with significantly divergent allele frequencies within or near 2081 genes and 3526 differentially expressed genes in one or both sexes. A total of 519 genes were both genetically divergent and differentially expressed between the divergent selection lines. We performed functional analyses of the effects of RNAi suppression of gene expression and induced mutations for 27 of these candidate genes that have human orthologs and the strongest statistical support, and confirmed that 25 (93 %) affected the mean and/or variance of food consumption.

Genomic Response to Selection for Predatory Behavior in a Mammalian Model of Adaptive Radiation

If genetic architectures of various quantitative traits are similar, as studies on model organisms suggest, comparable selection pressures should produce similar molecular patterns for various traits. To test this prediction, we used a laboratory model of vertebrate adaptive radiation to investigate the genetic basis of the response to selection for predatory behavior and compare it with evolution of aerobic capacity reported in an earlier work. After 13 generations of selection, the proportion of bank voles (Myodes [=Clethrionomys] glareolus) showing predatory behavior was five times higher in selected lines than in controls. We analyzed the hippocampus and liver transcriptomes and found repeatable changes in allele frequencies and gene expression. Genes with the largest differences between predatory and control lines are associated with hunger, aggression, biological rhythms, and functioning of the nervous system. Evolution of predatory behavior could be meaningfully compared with evolution of high aerobic capacity, because the experiments and analyses were performed in the same methodological framework. The number of genes that changed expression was much smaller in predatory lines, and allele frequencies changed repeatably in predatory but not in aerobic lines. This suggests that more variants of smaller effects underlie variation in aerobic performance, whereas fewer variants of larger effects underlie variation in predatory behavior. Our results thus contradict the view that comparable selection pressures for different quantitative traits produce similar molecular patterns. Therefore, to gain knowledge about molecular-level response to selection for complex traits, we need to investigate not only multiple replicate populations but also multiple quantitative traits.

A dopamine and noradrenaline reuptake inhibitor (bupropion) does not alter exercise performance of bank voles

Physical performance is determined both by biophysical and physiological limitations and behavioral characteristic, specifically motivation. We applied an experimental evolution approach combined with pharmacological manipulation to test the hypothesis that evolution of increased aerobic exercise performance can be triggered by evolution of motivation to undertake physical activity. We used a unique model system: bank voles from A lines, selected for high swim-induced aerobic metabolism (VO₂swim), which achieved a 61% higher mass-adjusted VO₂swim than those from unselected C lines. Because the voles could float on the water surface with only a minimum activity, the maximum rate of metabolism achieved in that test depended not only on their aerobic capacity, but also on motivation to undertake intensive activity. Therefore, we hypothesized that signaling of neurotransmitters putatively involved in regulating physical activity (dopamine and noradrenaline) had changed in response to selection. We measured VO₂swim after intraperitoneal injections of saline or the norepinephrine and dopamine reuptake inhibitor bupropion (20 mg/kg or 30 mg/kg). Additionally, we measured forced-exercise VO₂ (VO₂max). In C lines, VO₂swim (mass-adjusted mean ± standard error (SE): 4.0 ± 0.1 mLO₂/min) was lower than VO₂max (5.0 ± 0.1 mLO₂/min), but in A lines VO2swim (6.0 ± 0.1 mLO₂/min) was as high as VO₂max (6.0 ± 0.1 mLO₂/min). Thus, the selection effectively changed both the physiological-physical performance limit and mechanisms responsible for the willingness to undertake vigorous locomotor activity. Surprisingly, the drug had no effect on the achieved level of VO₂swim. Thus, the results did not allow firm conclusions concerning involvement of these neurotransmitters in evolution of increased aerobic exercise performance in the experimental evolution model system.

Experimental Evolution on a Wild Mammal Species Results in Modifications of Gut Microbial Communities

Comparative studies have shown that diet, life history, and phylogeny interact to determine microbial community structure across mammalian hosts. However, these studies are often confounded by numerous factors. Selection experiments offer unique opportunities to validate conclusions and test hypotheses generated by comparative studies. We used a replicated, 15-generation selection experiment on bank voles (Myodes glareolus) that have been selected for high swim-induced aerobic metabolism, predatory behavior toward crickets, and the ability to maintain body mass on a high-fiber, herbivorous diet. We predicted that selection on host performance, mimicking adaptive radiation, would result in distinct microbial signatures. We collected foregut and cecum samples from animals that were all fed the same nutrient-rich diet and had not been subjected to any performance tests. We conducted microbial inventories of gut contents by sequencing the V4 region of the 16S rRNA gene. We found no differences in cecal microbial community structure or diversity between control lines and the aerobic or predatory lines. However, the cecal chambers of voles selected for herbivorous capability harbored distinct microbial communities that exhibited higher diversity than control lines. The foregut communities of herbivorous-selected voles were also distinct from control lines. Overall, this experiment suggests that differences in microbial communities across herbivorous mammals may be evolved, and not solely driven by current diet or other transient factors.

Evolution of basal metabolic rate in bank voles from a multidirectional selection experiment

A major theme in evolutionary and ecological physiology of terrestrial vertebrates encompasses the factors underlying the evolution of endothermy in birds and mammals and interspecific variation of basal metabolic rate (BMR). Here, we applied the experimental evolution approach and compared BMR in lines of a wild rodent, the bank vole (Myodes glareolus), selected for 11 generations for: high swim-induced aerobic metabolism (A), ability to maintain body mass on a low-quality herbivorous diet (H) and intensity of predatory behaviour towards crickets (P). Four replicate lines were maintained for each of the selection directions and an unselected control (C). In comparison to C lines, A lines achieved a 49% higher maximum rate of oxygen consumption during swimming, H lines lost 1.3 g less mass in the test with low-quality diet and P lines attacked crickets five times more frequently. BMR was significantly higher in A lines than in C or H lines (60.8, 56.6 and 54.4 ml O2 h(-1), respectively), and the values were intermediate in P lines (59.0 ml O2 h(-1)). Results of the selection experiment provide support for the hypothesis of a positive association between BMR and aerobic exercise performance, but not for the association of adaptation to herbivorous diet with either a high or low BMR.

Limits to sustained energy intake. XXIII. Does heat dissipation capacity limit the energy budget of lactating bank voles?

Understanding factors limiting sustained metabolic rate (SusMR) is a central issue in ecological physiology. According to the heat dissipation limit (HDL) theory, the SusMR at peak lactation is constrained by the maternal capacity to dissipate body heat. To test that theory, we shaved lactating bank voles (Myodes glareolus) to experimentally elevate their capacity for heat dissipation. The voles were sampled from lines selected for high aerobic exercise metabolism (A; characterized also by increased basal metabolic rate) and unselected control lines (C). Fur removal significantly increased the peak-lactation food intake (mass-adjusted least square means ± s.e.; shaved: 16.3 ± 0.3 g day(-1), unshaved: 14.4 ± 0.2 g day(-1); P<0.0001), average daily metabolic rate (shaved: 109 ± 2 kJ day(-1), unshaved: 97 ± 2 kJ day(-1); P<0.0001) and metabolisable energy intake (shaved: 215 ± 4 kJ day(-1), unshaved: 185 ± 4 kJ day(-1); P<0.0001), as well as the milk energy output (shaved: 104 ± 4 kJ day(-1); unshaved: 93 ± 4 kJ day(-1); P=0.021) and litter growth rate (shaved: 9.4 ± 0.7 g 4 days(-1), unshaved: 7.7 ± 0.7 g 4 days(-1); P=0.028). Thus, fur removal increased both the total energy budget and reproductive output at the most demanding period of lactation, which supports the HDL theory. However, digestive efficiency was lower in shaved voles (76.0 ± 0.3%) than in unshaved ones (78.5 ± 0.2%; P<0.0001), which may indicate that a limit imposed by the capacity of the alimentary system was also approached. Shaving similarly affected the metabolic and reproductive traits in voles from the A and C lines. Thus, the experimental evolution model did not reveal a difference in the limiting mechanism between animals with inherently different metabolic rates.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

Reproduction is not costly in terms of oxidative stress

One of the core assumptions of life-history theory is the negative trade-off between current and future reproduction. Investment in current reproduction is expected to decrease future reproductive success or survival, but the physiological mechanisms underlying these costs are still obscure. To test for a role of oxidative stress, we measured oxidative damage to lipids and proteins in liver, heart, kidneys, and muscles, as well as the level of antioxidants (total glutathione and catalase), in breeding and non-breeding bank voles. We used females from lines selected for high aerobic metabolism and non-selected control lines and manipulated their reproductive investment by decreasing or increasing litter size. Unlike in most previous studies, the females reared four consecutive litters (the maximum possible during a breeding season). Contrary to predictions, oxidative damage in reproducing females was decreased or not changed, and did not differ between the selected and control lines. Oxidative damage to lipids and proteins in liver was lower in females that weaned enlarged litters than in non-breeding ones, and was intermediate in those with reduced litters. Oxidative damage to proteins in the heart also tended to be lower in breeding females than in non-breeding ones. A negative relationship between the level of oxidative damage and activity of catalase in kidneys indicated a protective action of antioxidants. In conclusion, our study falsified the hypothesis that oxidative stress is a part of the proximate physiological mechanism underlying the fundamental life-history trade-off between current and future reproduction.