The Mitochondrial Contribution to Animal Performance, Adaptation, and Life-History Variation

Animals display tremendous variation in their rates of growth, reproductive output, and longevity. While the physiological and molecular mechanisms that underlie this variation remain poorly understood, the performance of the mitochondrion has emerged as a key player. Mitochondria not only impact the performance of eukaryotes via their capacity to produce ATP, but they also play a role in producing heat and reactive oxygen species and function as a major signaling hub for the cell. The papers included in this special issue emerged from a symposium titled "Inside the Black Box: The Mitochondrial Basis of Life-history Variation and Animal Performance." Based on studies of diverse animal taxa, three distinct themes emerged from these papers. (1) When linking mitochondrial function to components of fitness, it is crucial that mitochondrial assays are performed in conditions as close as the intracellular conditions experienced by the mitochondria in vivo. (2) Functional plasticity allows mitochondria to retain their performance, as well as that of their host, over a range of exogenous conditions, and selection on mitochondrial and nuclear-derived proteins can optimize the match between the environment and the bioenergetic capacity of the mitochondrion. Finally, (3) studies of wild and wild-derived animals suggest that mitochondria play a central role in animal performance and life history strategy. Taken as a whole, we hope that these papers will foster discussion and inspire new hypotheses and innovations that will further our understanding of the mitochondrial processes that underlie variation in life history traits and animal performance.

Patterns of alternative splicing in response to cold acclimation in fish

Phenotypic plasticity is an important aspect of an organism's response to environmental change that often requires the modulation of gene expression. These changes in gene expression can be quantitative, as a result of increases or decreases in the amounts of specific transcripts, or qualitative, as a result of the expression of alternative transcripts from the same gene (e.g. via alternative splicing of pre-mRNAs). Although the role of quantitative changes in gene expression in phenotypic plasticity is well known, relatively few studies have examined the role of qualitative changes. Here, we use skeletal muscle RNA-seq data from Atlantic killifish (Fundulus heteroclitus), threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio) to investigate the extent of qualitative changes in gene expression in response to cold acclimation. Fewer genes demonstrated alternative splicing than differential expression as a result of cold acclimation; however, differences in splicing were detected for 426 to 866 genes depending on species, indicating that large numbers of qualitative changes in gene expression are associated with cold acclimation. Many of these alternatively spliced genes were also differentially expressed, and there was functional enrichment for involvement in muscle contraction among the genes demonstrating qualitative changes in response to cold acclimation. Additionally, there was a common group of 29 genes with cold-acclimation-mediated changes in splicing in all three species, suggesting that there may be a set of genes with expression patterns that respond qualitatively to prolonged exposure to cold temperatures across fishes.

Tolerance traits related to climate change resilience are independent and polygenic

The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish, Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine-tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

The DNA Methylation Landscape of Stickleback Reveals Patterns of Sex Chromosome Evolution and Effects of Environmental Salinity

Epigenetic mechanisms such as DNA methylation are a key component of dosage compensation on sex chromosomes and have been proposed as an important source of phenotypic variation influencing plasticity and adaptive evolutionary processes, yet little is known about the role of DNA methylation in an ecological or evolutionary context in vertebrates. The threespine stickleback (Gasterosteus aculeatus) is an ecological and evolutionary model system that has been used to study mechanisms involved in the evolution of adaptive phenotypes in novel environments as well as the evolution heteromorphic sex chromosomes and dosage compensation in vertebrates. Using whole genome bisulfite sequencing, we compared genome-wide DNA methylation patterns between threespine stickleback males and females and between stickleback reared at different environmental salinities. Apparent hypermethylation of the younger evolutionary stratum of the stickleback X chromosome in females relative to males suggests a potential role of DNA methylation in the evolution of heteromorphic sex chromosomes. We also demonstrate that rearing salinity has genome-wide effects on DNA methylation levels, which has the potential to lead to the accumulation of epigenetic variation between natural populations in different environments.

Persistent and plastic effects of temperature on DNA methylation across the genome of threespine stickleback (Gasterosteus aculeatus)

Epigenetic mechanisms such as changes in DNA methylation have the potential to affect the resilience of species to climate change, but little is known about the response of the methylome to changes in environmental temperature in animals. Using reduced representation bisulfite sequencing, we assessed the effects of development temperature and adult acclimation temperature on DNA methylation levels in threespine stickleback (Gasterosteus aculeatus). Across all treatments, we identified 2130 differentially methylated cytosines distributed across the genome. Both increases and decreases in temperature during development and with thermal acclimation in adults increased global DNA methylation levels. Approximately 25% of the differentially methylated regions (DMRs) responded to both developmental temperature and adult thermal acclimation, and 50 DMRs were common to all treatments, demonstrating a core response of the epigenome to thermal change at multiple time scales. We also identified differentially methylated loci that were specific to a particular developmental or adult thermal response, which could facilitate the accumulation of epigenetic variation between natural populations that experience different thermal regimes. These data demonstrate that thermal history can have long-lasting effects on the epigenome, highlighting the role of epigenetic modifications in the response to temperature change across multiple time scales.

Similarities in temperature-dependent gene expression plasticity across timescales in threespine stickleback (Gasterosteus aculeatus)

Phenotypic plasticity occurs at a variety of timescales, but little is known about the degree to which plastic responses at different timescales are associated with similar underlying molecular processes, which is critical for assessing the effects of plasticity on evolutionary trajectories. To address this issue, we identified differential gene expression in response to developmental temperature in the muscle transcriptome of adult threespine stickleback (Gasterosteus aculeatus) exposed to 12, 18 and 24°C until hatch and then held at 18°C for 9 months and compared these results to differential gene expression in response to adult thermal acclimation in stickleback developed at 18°C and then acclimated to 5 and 25°C as adults. Adult thermal acclimation affected the expression of 7,940 and 7,015 genes in response to cold and warm acclimation, respectively, and 4,851 of these genes responded in both treatments. In contrast, the expression of only 33 and 29 genes was affected by cold and warm development, respectively. The majority of the genes affected by developmental temperature were also affected by adult acclimation temperature. Many genes that were differentially expressed as a result of adult acclimation were associated with previously identified temperature-dependent effects on DNA methylation patterns, suggesting a role of epigenetic mechanisms in regulating gene expression plasticity during acclimation. Taken together, these results demonstrate similarities between the persistent effects of developmental plasticity on gene expression and the effects of adult thermal acclimation, emphasizing the potential for mechanistic links between plasticity acting at these different life stages.

Patterns of mitochondrial membrane remodeling parallel functional adaptations to thermal stress

The effect of temperature on mitochondrial performance is thought to be partly due to its effect on mitochondrial membranes. Numerous studies have shown that thermal acclimation and adaptation can alter the amount of inner-mitochondrial membrane (IMM), but little is known about the capacity of organisms to modulate mitochondrial membrane composition. Using northern and southern subspecies of Atlantic killifish (Fundulus heteroclitus) that are locally adapted to different environmental temperatures, we assessed whether thermal acclimation altered liver mitochondrial respiratory capacity or the composition and amount of IMM. We measured changes in phospholipid headgroups and headgroup-specific fatty acid (FA) remodeling, and used respirometry to assess mitochondrial respiratory capacity. Acclimation to 5°C and 33°C altered mitochondrial respiratory capacity in both subspecies. Northern F. heteroclitus exhibited greater mitochondrial respiratory capacity across acclimation temperatures, consistent with previously observed subspecies differences in whole-organism aerobic metabolism. Mitochondrial phospholipids were altered following thermal acclimation, and the direction of these changes was largely consistent between subspecies. These effects were primarily driven by remodeling of specific phospholipid classes and were associated with shifts in metabolic phenotypes. There were also differences in membrane composition between subspecies that were driven largely by differences in phospholipid classes. Changes in respiratory capacity between subspecies and with acclimation were largely but not completely accounted for by alterations in the amount of IMM. Taken together, these results support a role for changes in liver mitochondrial function in the ectothermic response to thermal stress during both acclimation and adaptation, and implicate lipid remodeling as a mechanism contributing to these changes.

Anatomy of STEM teaching in North American universities

National and local initiatives focused on the transformation of STEM teaching in higher education have multiplied over the last decade. These initiatives often focus on measuring change in instructional practices, but it is difficult to monitor such change without a national picture of STEM educational practices, especially as characterized by common observational instruments. We characterized a snapshot of this landscape by conducting the first large scale observation-based study. We found that lecturing was prominent throughout the undergraduate STEM curriculum, even in classrooms with infrastructure designed to support active learning, indicating that further work is required to reform STEM education. Additionally, we established that STEM faculty’s instructional practices can vary substantially within a course, invalidating the commonly-used teaching evaluations based on a one-time observation.

Although lecture is prominent throughout the undergraduate STEM curriculum, STEM faculty employ varied teaching practices within the same course.

A rainbow trout Oncorhynchus mykiss strain with higher aerobic scope in normoxia also has superior tolerance of hypoxia

This study compared parr from three strains of rainbow trout Oncorhynchus mykiss to examine intraspecific variation in metabolic traits, hypoxia tolerance and upper thermal tolerance in this species. At the strain level, variation in absolute aerobic scope (AAS), critical oxygen level (O_2crit ), incipient lethal oxygen saturation (ILOS) and critical thermal maximum (CT_max ) generally exhibited consistent differences among the strains, suggesting the possibility of functional associations among these traits. This possibility was further supported at the individual level by a positive correlation between ILOS and O_2crit and a negative correlation between O_2crit and AAS. These results indicate that intraspecific differences in hypoxia tolerance among strains of O. mykiss may be primarily determined by differences in the ability to maintain oxygen uptake in hypoxia and that variation in aerobic scope in normoxia probably plays a role in determining the ability of these fish to sustain metabolism aerobically as water oxygen saturation is reduced.

claudin-10 isoform expression and cation selectivity change with salinity in salt-secreting epithelia of Fundulusheteroclitus

To provide insight into claudin (Cldn) tight junction (TJ) protein contributions to branchial salt secretion in marine teleost fishes, this study examined cldn-10 TJ protein isoforms of a euryhaline teleost (mummichog; Fundulus heteroclitus) in association with salinity change and measurements of transepithelial cation selectivity. Mummichogs were transferred from freshwater (FW) to seawater (SW, 35‰) and from SW to hypersaline SW (2SW, 60‰) in a time course with transfer control groups (FW to FW, and SW to SW). FW to SW transfer increased mRNA abundance of cldn-10d and cldn-10e twofold, whilst cldn-10c and cldn-10f transcripts were unchanged. Transfer from SW to 2SW did not alter cldn-10d, and transiently altered cldn-10e abundance, but increased cldn-10c and cldn-10f fourfold. This was coincident with an increased number of single-stranded junctions (observed by transmission electron microscopy). For both salinity transfers, (1) cldn-10e mRNA was acutely responsive (i.e. after 24 h), (2) other responsive cldn-10 isoforms increased later (3-7 days), and (3) cystic fibrosis transmembrane conductance regulator (cftr) mRNA was elevated in accordance with established changes in transcellular Cl^- movement. Changes in mRNA encoding cldn-10c and -10f appeared linked, consistent with the tandem repeat locus in the Fundulus genome, whereas mRNA for tandem cldn-10d and cldn-10e seemed independent of each other. Cation selectivity sequence measured by voltage and conductance responses to artificial SW revealed Eisenman sequence VII: Na⁺>K⁺>Rb⁺∼Cs⁺>Li⁺ Collectively, these data support the idea that Cldn-10 TJ proteins create and maintain cation-selective pore junctions in salt-secreting tissues of teleost fishes.

Subspecies differences in thermal acclimation of mitochondrial function and the role of uncoupling proteins in killifish

Thermal effects on mitochondrial efficiency and ATP production can influence whole-animal thermal tolerance and performance. Thus, organisms may have the capacity to alter mitochondrial processes through acclimation or adaptation to mitigate these effects. One possible mechanism is through the action of uncoupling proteins (UCPs) which can decrease the proton motive force independent of the production of ATP. To test this hypothesis, we examined the mRNA expression patterns of UCP isoforms and characterized the effects of thermal acclimation and putative local thermal adaptation on mitochondrial capacity, proton leak, and P/O ratios in two subspecies of Atlantic killifish (Fundulus heteroclitus). Ucp1 was the dominant isoform in liver and was more highly expressed in northern killifish. We found that cold acclimation increased mitochondrial capacity (state III and maximum substrate oxidation capacity), state II membrane potential, proton leak, and P/O ratios in northern, but not southern killifish liver mitochondria. Palmitate-induced mitochondrial uncoupling was detected in northern, but not southern, killifish liver mitochondria, consistent with the differences in mRNA expression between the subspecies. Taken together, our data suggest that mitochondrial function is more plastic in response to thermal acclimation in northern killifish than southern killifish and that UCP1 may play a role in regulating the proton motive force in northern, but not southern killifish in response to thermal acclimation. These data demonstrate the potential for adaptive variation in mitochondrial plasticity in response to cold.

Functional Annotation of All Salmonid Genomes (FAASG): an international initiative supporting future salmonid research, conservation and aquaculture

We describe an emerging initiative - the 'Functional Annotation of All Salmonid Genomes' (FAASG), which will leverage the extensive trait diversity that has evolved since a whole genome duplication event in the salmonid ancestor, to develop an integrative understanding of the functional genomic basis of phenotypic variation. The outcomes of FAASG will have diverse applications, ranging from improved understanding of genome evolution, to improving the efficiency and sustainability of aquaculture production, supporting the future of fundamental and applied research in an iconic fish lineage of major societal importance.

Intrinsic reproductive isolating mechanisms in the maintenance of a hybrid zone between ecologically divergent subspecies

Understanding factors involved in maintaining stable hybrid zones is important for predicting the ultimate fate of the interacting taxa, but the relative importance of mechanisms such as ecological selection and intrinsic reproductive isolation remains unclear. Most studies of reproductive isolation in hybrid zones have focused either on zones with strongly bimodal patterns in genotype or phenotype frequencies, with relatively strong isolation, or unimodal zones with relatively weak isolation, whereas less is known about more intermediate classes of hybrid zone. Here, we utilize a hybrid zone of this intermediate type occurring between northern and southern subspecies of Atlantic killifish, Fundulus heteroclitus, to identify isolating mechanisms playing a role in maintaining this type of zone. The two subspecies differ in environmental tolerance, and we found some evidence of microhabitat preference between subspecies within a small tidal creek at the centre of the hybrid zone. There was also an association between sex, mitochondrial genotype and habitat within this creek. Fertilization success did not differ between consubspecific and heterosubspecific crosses, but hatching success was significantly lower for crosses involving southern males and northern females, and crosses between southern females and northern males had altered developmental rates. Southern females and northern males showed patterns consistent with positive assortative mating. Together, these results indicate a role for a combination of factors including assortative mating and/or early hybrid inviability in the maintenance of this hybrid zone and suggest that hybrid zones with intermediate levels of reproductive isolation are likely to be maintained by multiple interacting isolating mechanisms.

Warm acclimation improves hypoxia tolerance in Fundulus heteroclitus

Human activities are increasing both the frequency of hypoxic episodes and the mean temperature of aquatic ecosystems, but few studies have considered the possibility that acclimation to one of these stressors could improve the ability to cope with the other stressor. Here, we used Atlantic killifish, Fundulus heteroclitus, to test this hypothesis. Hypoxia tolerance was measured as time to loss of equilibrium in hypoxia (LOEhyp) at 0.4 kPa oxygen. Time to LOEhyp declined from 73.3 ± 6.9 min at 15 °C to 2.6 ± 3.8 min at 23 °C, and at 30 °C no fish could withstand this level of hypoxia. Prior acclimation to warm temperatures significantly increased time to LOEhyp. Hypoxia tolerance of the southern subspecies of killifish, F. heteroclitus heteroclitus, was greater than that of the northern subspecies, F. heteroclitus macrolepidotus, measured both as critical oxygen tension (Pcrit) and as time to LOEhyp. Warm acclimation offset the negative effects of temperature on time to LOEhyp to a similar extent in the two subspecies. Warm acclimation increased total lamellar surface area of the gill in both subspecies as a result of regression of an interlamellar cell mass (ILCM). However, differences in total lamellar surface area could not explain differences in time to LOEhyp between the subspecies, suggesting that the lower time to LOEhyp of northern fish is related to their higher routine metabolic rate. These data suggest that thermal plasticity in gill morphology can improve the capacity of this species to tolerate hypoxia, and shows how existing plasticity may help organisms to cope with the complex interacting stressors that they will encounter with increasing frequency as our climate changes.

Thermal acclimation and subspecies-specific effects on heart and brain mitochondrial performance in a eurythermal teleost (Fundulus heteroclitus)

Mitochondrial performance may play a role in setting whole-animal thermal tolerance limits and their plasticity, but the relative roles of adjustments in mitochondrial performance across different highly aerobic tissues remain poorly understood. We compared heart and brain mitochondrial responses to acute thermal challenges and to thermal acclimation using high-resolution respirometry in two locally adapted subspecies of Atlantic killifish (Fundulus heteroclitus). We predicted that 5°C acclimation to would result in compensatory increases in mitochondrial performance, while 33°C acclimation would cause suppression of mitochondrial function to minimize the effects of high temperature on mitochondrial metabolism. In contrast, acclimation to both 33 and 5°C decreased mitochondrial performance compared to fish acclimated to 15°C. These adjustments could represent an energetic cost saving mechanism at temperature extremes. Acclimation responses were similar in both heart and brain; however, this effect was smaller in the heart which might indicate its importance in maintaining whole-animal thermal performance. Alternatively, larger acclimation effects in the brain might indicate greater thermal sensitivity compared to the heart. We detected only modest differences between subspecies that were dependent on the tissue assayed. These data demonstrate extensive plasticity in mitochondrial performance following thermal acclimation in killifish, and indicate that the extent of these responses differs between tissues, highlighting the importance and complexity of mitochondrial regulation in thermal acclimation in eurytherms.

Maternal stress has divergent effects on gene expression patterns in the brains of male and female threespine stickleback

Maternal stress can have long-term effects on neurodevelopment that can influence offspring performance and population evolutionary trajectories. To examine the mechanistic basis for these neurodevelopmental effects of maternal stress, we used RNA-seq to assess differential gene expression across the brain transcriptome of adult male and female threespine stickleback (Gasterosteus aculeatus) from stressed and unstressed mothers. We identified sexually divergent effects of maternal stress on the brain transcriptome. In males, genes that were upregulated by maternal stress were enriched for processes involved in synaptic function and organization and steroid hormone-mediated signalling pathways, whereas in females genes that were upregulated by maternal stress were enriched for processes involved in protein translation and metabolic functions. The expression of several genes involved in the hypothalamic-pituitary-interrenal response to stress and epigenetic processes such as the regulation of DNA methylation patterns and miRNAs increased in males and not in females. These data suggest that maternal stress has markedly different effects on cellular pathways in the brains of male and female offspring of mothers that are exposed to stress, which could have important implications when assessing the long-term ecological and evolutionary impacts of stress across generations.