Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a QST - FST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.Peer reviewe

Influence of interspecific interference competition on the genetic structure of Calopteryx splendens populations

Understanding the effects of interspecific competition on genetic diversity will deepen our knowledge on species evolution. In the case of Calopteryx splendens and C. virgo, sympatric damselfly species, interspecific interference competition by C. virgo has remarkable effects on territoriality of C. splendens resulting in reproductive character displacement. Since territoriality is correlated with phenotype and mating success, we investigated the effects of interspecific interference competition on genetic diversity of C. splendens populations. Using amplified fragment length polymorphisms (AFLP), we determined the population genetic structure of 12 C. splendens populations and used the genetic diversity information to relate heterozygosity of C. splendens to abundance of C. virgo in sympatric populations. We found that heterozygosity of C. splendens males decreased with increasing abundance of C. virgo males. This result most likely reflects changes in effective population size due to interspecific interference competition and shows an effect on genetic structure in damselfly populations.</p

Cytosine methylation patterns suggest a role of methylation in plastic and adaptive responses to temperature in European grayling (Thymallus thymallus) populations

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.Peer reviewe

Elasmobranch qPCR reference genes: a case study of hypoxia preconditioned epaulette sharks

<p>Abstract</p> <p>Background</p> <p>Elasmobranch fishes are an ancient group of vertebrates which have high potential as model species for research into evolutionary physiology and genomics. However, no comparative studies have established suitable reference genes for quantitative PCR (qPCR) in elasmobranchs for any physiological conditions. Oxygen availability has been a major force shaping the physiological evolution of vertebrates, especially fishes. Here we examined the suitability of 9 reference candidates from various functional categories after a single hypoxic insult or after hypoxia preconditioning in epaulette shark (<it>Hemiscyllium ocellatum</it>).</p> <p>Results</p> <p>Epaulette sharks were caught and exposed to hypoxia. Tissues were collected from 10 controls, 10 individuals with single hypoxic insult and 10 individuals with hypoxia preconditioning (8 hypoxic insults, 12 hours apart). We produced sequence information for reference gene candidates and monitored mRNA expression levels in four tissues: cerebellum, heart, gill and eye. The stability of the genes was examined with analysis of variance, geNorm and NormFinder. The best ranking genes in our study were <it>eukaryotic translation elongation factor 1 beta </it>(<it>eef1b</it>), <it>ubiquitin </it>(<it>ubq</it>) and <it>polymerase (RNA) II (DNA directed) polypeptide F </it>(<it>polr2f</it>). The performance of the <it>ribosomal protein L6 </it>(<it>rpl6</it>) was tissue-dependent. Notably, in one tissue the analysis of variance indicated statistically significant differences between treatments for genes that were ranked as the most stable candidates by reference gene software.</p> <p>Conclusions</p> <p>Our results indicate that <it>eef1b </it>and <it>ubq </it>are generally the most suitable reference genes for the conditions and tissues in the present epaulette shark studies. These genes could also be potential reference gene candidates for other physiological studies examining stress in elasmobranchs. The results emphasise the importance of inter-group variation in reference gene evaluation.</p

Measuring telomere length and telomere dynamics in evolutionary biology and ecology

Telomeres play a fundamental role in the protection of chromosomal DNA and in the regulation of cellular senescence. Recent work in human epidemiology and evolutionary ecology suggests adult telomere length (TL) may reflect past physiological stress and predict subsequent morbidity and mortality, independent of chronological age. Several different methods have been developed to measure TL, each offering its own technical challenges. The aim of this review is to provide an overview of the advantages and drawbacks of each method for researchers, with a particular focus on issues that are likely to face ecologists and evolutionary biologists collecting samples in the field or in organisms that may never have been studied in this context before. We discuss the key issues to consider and wherever possible try to provide current consensus view regarding best practice with regard to sample collection and storage, DNA extraction and storage, and the five main methods currently available to measure TL. Decisions regarding which tissues to sample, how to store them, how to extract DNA, and which TL measurement method to use cannot be prescribed, and are dependent on the biological question addressed and the constraints imposed by the study system. What is essential for future studies of telomere dynamics in evolution and ecology is that researchers publish full details of their methods and the quality control thresholds they employ

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Transe-QTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.Peer reviewe

Immigrant reproductive dysfunction facilitates ecological speciation

The distributions of species are not only determined by where they can survive – they must also be able to reproduce. Although immigrant inviability is a well-established concept, the fact that immigrants also need to be able to effectively reproduce in foreign environments has not been fully appreciated in the study of adaptive divergence and speciation. Fertilization and reproduction are sensitive life-history stages that could be detrimentally affected for immigrants in non-native habitats. We propose that “immigrant reproductive dysfunction” is a hitherto overlooked aspect of reproductive isolation caused by natural selection on immigrants. This idea is supported by results from experiments on an externally fertilizing fish (sand goby, Pomatoschistus minutus). Growth and condition of adults were not affected by non-native salinity whereas males spawning as immigrants had lower sperm motility and hatching success than residents. We interpret these results as evidence for local adaptation or acclimation of sperm, and possibly also components of paternal care. The resulting loss in fitness, which we call “immigrant reproductive dysfunction,” has the potential to reduce gene flow between populations with locally adapted reproduction, and it may play a role in species distributions and speciation

Sperm adaptation in relation to salinity in three goby species

In externally fertilizing species, the gametes of both males and females are exposed to the influences of the environment into which they are released. Sperm are sensitive to abiotic factors such as salinity, but they are also affected by biotic factors such as sperm competition. In this study, the authors compared the performance of sperm of three goby species, the painted goby, Pomatoschistus pictus, the two-spotted goby, Pomatoschistus flavescens, and the sand goby, Pomatoschistus minutus. These species differ in their distributions, with painted goby having the narrowest salinity range and sand goby the widest. Moreover, data from paternity show that the two-spotted goby experiences the least sperm competition, whereas in the sand goby sperm competition is ubiquitous. The authors took sperm samples from dissected males and exposed them to high salinity water (31 PSU) representing the North Sea and low salinity water (6 PSU) representing the brackish Baltic Sea Proper. They then used computer-assisted sperm analysis to measure the proportion of motile sperm and sperm swimming speed 10 min and 20 h after sperm activation. The authors found that sperm performance depended on salinity, but there seemed to be no relationship to the species' geographical distribution in relation to salinity range. The species differed in the proportion of motile sperm, but there was no significant decrease in sperm motility during 20 h. The sand goby was the only species with motile sperm after 72 h