Genomic and Epigenomic Influences on Resilience Across Scales: lessons from the responses of fishes to environmental stressors.

Understanding the factors that influence the resilience of biological systems to environmental change is an increasingly critical concern in the face of increasing human impacts on ecosystems and the organisms that inhabit them. However, most considerations of biological resilience have focused at the community and ecosystem levels, whereas here we discuss how including consideration of processes occurring at lower levels of biological organization may provide insights into factors that influence resilience at higher levels of biological organization. Specifically, we explore how processes at the genomic and epigenomic levels may cascade up to influence resilience at higher levels and ask how the concepts of 'resistance', or the capacity of a system to minimize change in response to a disturbance, and 'recovery', or the ability of a system to return to its original state following a disturbance and avoid tipping points and resulting regime shifts map to these lower levels of biological organization. Overall, we argue that substantial changes at these lower levels may be required to support resilience at higher levels, using selected examples of genomic and epigenomic responses of fish to climate-change relevant stressors such as high temperature and hypoxia at the levels of the genome, epigenome and organism.

Transposon wave remodeled the epigenomic landscape in the rapid evolution of X-Chromosome dosage compensation.

Sex chromosome dosage compensation is a model to understand the coordinated evolution of transcription; however, the advanced age of the sex chromosomes in model systems makes it difficult to study how the complex regulatory mechanisms underlying chromosome-wide dosage compensation can evolve. The sex chromosomes of Poecilia picta have undergone recent and rapid divergence, resulting in widespread gene loss on the male Y, coupled with complete X Chromosome dosage compensation, the first case reported in a fish. The recent de novo origin of dosage compensation presents a unique opportunity to understand the genetic and evolutionary basis of coordinated chromosomal gene regulation. By combining a new chromosome-level assembly of P. picta with whole-genome bisulfite sequencing and RNA-seq data, we determine that the YY1 transcription factor (YY1) DNA binding motif is associated with male-specific hypomethylated regions on the X, but not the autosomes. These YY1 motifs are the result of a recent and rapid repetitive element expansion on the P. picta X Chromosome, which is absent in closely related species that lack dosage compensation. Taken together, our results present compelling support that a disruptive wave of repetitive element insertions carrying YY1 motifs resulted in the remodeling of the X Chromosome epigenomic landscape and the rapid de novo origin of a dosage compensation system.

Sex chromosome heteromorphism and the Fast-X effect in poeciliids.

Fast-X evolution has been observed in a range of heteromorphic sex chromosomes. However, it remains unclear how early in the process of sex chromosome differentiation the Fast-X effect becomes detectible. Recently, we uncovered an extreme variation in sex chromosome heteromorphism across poeciliid fish species. The common guppy, Poecilia reticulata, Endler's guppy, P. wingei, swamp guppy, P. picta and para guppy, P. parae, appear to share the same XY system and exhibit a remarkable range of heteromorphism. Species outside this group lack this sex chromosome system. We combined analyses of sequence divergence and polymorphism data across poeciliids to investigate X chromosome evolution as a function of hemizygosity and reveal the causes for Fast-X effects. Consistent with the extent of Y degeneration in each species, we detect higher rates of divergence on the X relative to autosomes, a signal of Fast-X evolution, in P. picta and P. parae, species with high levels of X hemizygosity in males. In P. reticulata, which exhibits largely homomorphic sex chromosomes and little evidence of hemizygosity, we observe no change in the rate of evolution of X-linked relative to autosomal genes. In P. wingei, the species with intermediate sex chromosome differentiation, we see an increase in the rate of nonsynonymous substitutions on the older stratum of divergence only. We also use our comparative approach to test for the time of origin of the sex chromosomes in this clade. Taken together, our study reveals an important role of hemizygosity in Fast-X evolution.

Evolutionary History of the Poecilia picta Sex Chromosomes.

The degree of divergence between the sex chromosomes is not always proportional to their age. In poeciliids, four closely related species all exhibit a male heterogametic sex chromosome system on the same linkage group, yet show a remarkable diversity in X and Y divergence. In Poecilia reticulata and P. wingei, the sex chromosomes remain homomorphic, yet P. picta and P. parae have a highly degraded Y chromosome. To test alternative theories about the origin of their sex chromosomes, we used a combination of pedigrees and RNA-seq data from P. picta families in conjunction with DNA-seq data collected from P. reticulata, P. wingei, P. parae, and P. picta. Phylogenetic clustering analysis of X and Y orthologs, identified through segregation patterns, and their orthologous sequences in closely related species demonstrates a similar time of origin for both the P. picta and P. reticulata sex chromosomes. We next used k-mer analysis to identify shared ancestral Y sequence across all four species, suggesting a single origin to the sex chromosome system in this group. Together, our results provide key insights into the origin and evolution of the poeciliid Y chromosome and illustrate that the rate of sex chromosome divergence is often highly heterogenous, even over relatively short evolutionary time frames.

Rapid Evolution of Complete Dosage Compensation in Poecilia.

Dosage compensation balances gene expression between the sexes in systems with diverged heterogametic sex chromosomes. Theory predicts that dosage compensation should rapidly evolve in tandem with the divergence of sex chromosomes to prevent the deleterious effects of dosage imbalances that occur as a result of sex chromosome divergence. Examples of complete dosage compensation, where gene expression of the entire sex chromosome is compensated, are rare, and have only been found in relatively ancient sex chromosome systems. Consequently, very little is known about the evolutionary dynamics of complete dosage compensation systems. Within the family Poeciliidae the subgenus Lebistes share the same sex chromosome system which originated 18.48-26.08 Ma. In Poecilia reticulata and P. wingei, the Y chromosome has been largely maintained, whereas the Y in the closely related species P. picta and P. parae has rapidly degraded. We recently found P. picta to be the first example of complete dosage compensation in a fish. Here, we show that P. parae also has complete dosage compensation, thus complete dosage compensation likely evolved in the short (∼3.7 Myr) interval after the split of the ancestor of these two species from P. reticulata, but before they diverged from each other. These data suggest that novel dosage compensation mechanisms can evolve rapidly, thus supporting the longstanding theoretical prediction that such mechanisms arise in tandem with rapidly diverging sex chromosomes.

Sex Chromosome Evolution: So Many Exceptions to the Rules.

Genomic analysis of many nonmodel species has uncovered an incredible diversity of sex chromosome systems, making it possible to empirically test the rich body of evolutionary theory that describes each stage of sex chromosome evolution. Classic theory predicts that sex chromosomes originate from a pair of homologous autosomes and recombination between them is suppressed via inversions to resolve sexual conflict. The resulting degradation of the Y chromosome gene content creates the need for dosage compensation in the heterogametic sex. Sex chromosome theory also implies a linear process, starting from sex chromosome origin and progressing to heteromorphism. Despite many convergent genomic patterns exhibited by independently evolved sex chromosome systems, and many case studies supporting these theoretical predictions, emerging data provide numerous interesting exceptions to these long-standing theories, and suggest that the remarkable diversity of sex chromosomes is matched by a similar diversity in their evolution. For example, it is clear that sex chromosome pairs are not always derived from homologous autosomes. In addition, both the cause and the mechanism of recombination suppression between sex chromosome pairs remain unclear, and it may be that the spread of recombination suppression is a more gradual process than previously thought. It is also clear that dosage compensation can be achieved in many ways, and displays a range of efficacy in different systems. Finally, the remarkable turnover of sex chromosomes in many systems, as well as variation in the rate of sex chromosome divergence, suggest that assumptions about the inevitable linearity of sex chromosome evolution are not always empirically supported, and the drivers of the birth-death cycle of sex chromosome evolution remain to be elucidated. Here, we concentrate on how the diversity in sex chromosomes across taxa highlights an equal diversity in each stage of sex chromosome evolution.

A solution to Nature's haemoglobin knockout: a plasma-accessible carbonic anhydrase catalyses CO2 excretion in Antarctic icefish gills.

In all vertebrates studied to date, CO2 excretion depends on the enzyme carbonic anhydrase (CA) that catalyses the rapid conversion of HCO3- to CO2 at the gas-exchange organs. The largest pool of CA is present within red blood cells (RBCs) and, in some vertebrates, plasma-accessible CA (paCA) isoforms participate in CO2 excretion. However, teleost fishes typically do not have paCA at the gills and CO2 excretion is reliant entirely on RBC CA - a strategy that is not possible in icefishes. As the result of a natural knockout, Antarctic icefishes (Channichthyidae) are the only known vertebrates that do not express haemoglobin (Hb) as adults, and largely lack RBCs in the circulation (haematocrit <1%). Previous work has indicated the presence of high levels of membrane-bound CA activity in the gills of icefishes, but without determining its cellular orientation. Thus, we hypothesised that icefishes express a membrane-bound CA isoform at the gill that is accessible to the blood plasma. The CA distribution was compared in the gills of two closely related notothenioid species, one with Hb and RBCs (Notothenia rossii) and one without (Champsocephalus gunnari). Molecular, biochemical and immunohistochemical markers indicate high levels of a Ca4 isoform in the gills of the icefish (but not the red-blooded N. rossii), in a plasma-accessible location that is consistent with a role in CO2 excretion. Thus, in the absence of RBC CA, the icefish gill could exclusively provide the catalytic activity necessary for CO2 excretion - a pathway that is unlike that of any other vertebrate.

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.

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.

Gene expression plasticity in response to salinity acclimation in threespine stickleback ecotypes from different salinity habitats.

Phenotypic plasticity is thought to facilitate the colonization of novel environments and shape the direction of evolution in colonizing populations. However, the relative prevalence of various predicted patterns of changes in phenotypic plasticity following colonization remains unclear. Here, we use a whole-transcriptome approach to characterize patterns of gene expression plasticity in the gills of a freshwater-adapted and a saltwater-adapted ecotype of threespine stickleback (Gasterosteus aculeatus) exposed to a range of salinities. The response of the gill transcriptome to environmental salinity had a large shared component common to both ecotypes (2159 genes) with significant enrichment of genes involved in transmembrane ion transport and the restructuring of the gill epithelium. This transcriptional response to freshwater acclimation is induced at salinities below two parts per thousand. There was also differentiation in gene expression patterns between ecotypes (2515 genes), particularly in processes important for changes in the gill structure and permeability. Only 508 genes that differed between ecotypes also responded to salinity and no specific processes were enriched among this gene set, and an even smaller number (87 genes) showed evidence of changes in the extent of the response to salinity acclimation between ecotypes. No pattern of relative expression dominated among these genes, suggesting that neither gains nor losses of plasticity dominated the changes in expression patterns between the ecotypes. These data demonstrate that multiple patterns of changes in gene expression plasticity can occur following colonization of novel habitats.

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.

Conserved effects of salinity acclimation on thermal tolerance and hsp70 expression in divergent populations of threespine stickleback (Gasterosteus aculeatus).

In natural environments, organisms must cope with complex combinations of abiotic stressors. Here, we use threespine stickleback (Gasterosteus aculeatus) to examine how changes in salinity affect tolerance of high temperatures. Threespine stickleback inhabit a range of environments that vary in both salinity and thermal stability making this species an excellent system for investigating interacting stressors. We examined the effects of environmental salinity on maximum thermal tolerance (CTMax) and 70 kDa heat shock protein (hsp70) gene expression using divergent stickleback ecotypes from marine and freshwater habitats. In both ecotypes, the CTMax of fish acclimated to 20 ppt was significantly higher compared to fish acclimated to 2 ppt. The effect of salinity acclimation on the expression of hsp70-1 and hsp70-2 was similar in both the marine and freshwater stickleback ecotype. There were differences in the expression profiles of hsp70-1 and hsp70-2 during heat shock, with hsp70-2 being induced earlier and to a higher level compared to hsp70-1. These data suggest that the two hsp70 isoforms may have functionally different roles in the heat shock response. Lastly, acute salinity challenge coupled with heat shock revealed that the osmoregulatory demands experienced during the heat shock response have a larger effect on the hsp70 expression profile than does the acclimation salinity.

Epigenomics in marine fishes.

Epigenetic mechanisms are an underappreciated and often ignored component of an organism's response to environmental change and may underlie many types of phenotypic plasticity. Recent technological advances in methods for detecting epigenetic marks at a whole-genome scale have launched new opportunities for studying epigenomics in ecologically relevant non-model systems. The study of ecological epigenomics holds great promise to better understand the linkages between genotype, phenotype, and the environment and to explore mechanisms of phenotypic plasticity. The many attributes of marine fish species, including their high diversity, variable life histories, high fecundity, impressive plasticity, and economic value provide unique opportunities for studying epigenetic mechanisms in an environmental context. To provide a primer on epigenomic research for fish biologists, we start by describing fundamental aspects of epigenetics, focusing on the most widely studied and most well understood of the epigenetic marks: DNA methylation. We then describe the techniques that have been used to investigate DNA methylation in marine fishes to date and highlight some new techniques that hold great promise for future studies. Epigenomic research in marine fishes is in its early stages, so we first briefly discuss what has been learned about the establishment, maintenance, and function of DNA methylation in fishes from studies in zebrafish and then summarize the studies demonstrating the pervasive effects of the environment on the epigenomes of marine fishes. We conclude by highlighting the potential for ongoing research on the epigenomics of marine fishes to reveal critical aspects of the interaction between organisms and their environments.

Conserved structure and expression of hsp70 paralogs in teleost fishes.

The cytosolic 70KDa heat shock proteins (Hsp70s) are widely used as biomarkers of environmental stress in ecological and toxicological studies in fish. Here we analyze teleost genome sequences to show that two genes encoding inducible hsp70s (hsp70-1 and hsp70-2) are likely present in all teleost fish. Phylogenetic and synteny analyses indicate that hsp70-1 and hsp70-2 are distinct paralogs that originated prior to the diversification of the teleosts. The promoters of both genes contain a TATA box and conserved heat shock elements (HSEs), but unlike mammalian HSP70s, both genes contain an intron in the 5' UTR. The hsp70-2 gene has undergone tandem duplication in several species. In addition, many other teleost genome assemblies have multiple copies of hsp70-2 present on separate, small, genomic scaffolds. To verify that these represent poorly assembled tandem duplicates, we cloned the genomic region surrounding hsp70-2 in Fundulus heteroclitus and showed that the hsp70-2 gene copies that are on separate scaffolds in the genome assembly are arranged as tandem duplicates. Real-time quantitative PCR of F. heteroclitus genomic DNA indicates that four copies of the hsp70-2 gene are likely present in the F. heteroclitus genome. Comparison of expression patterns in F. heteroclitus and Gasterosteus aculeatus demonstrates that hsp70-2 has a higher fold increase than hsp70-1 following heat shock in gill but not in muscle tissue, revealing a conserved difference in expression patterns between isoforms and tissues. These data indicate that ecological and toxicological studies using hsp70 as a biomarker in teleosts should take this complexity into account.