Destabilized host-parasite dynamics in newly founded populations.

When species disperse into previously unoccupied habitats, new populations encounter unfamiliar species interactions such as altered parasite loads. Theory predicts that newly founded populations should exhibit destabilized eco-evolutionary fluctuations in infection rates and immune traits. However, to understand founder effects biologists typically rely on retrospective studies of range expansions, missing early-generation infection dynamics. To remedy this, we experimentally founded whole-lake populations of threespine stickleback. Infection rates were temporally stable in native source lakes. In contrast, newly founded populations exhibit destabilized host-parasite dynamics: high starting infection rates led to increases in a heritable immune trait (peritoneal fibrosis), suppressing infection rates. The resulting temporal auto-correlation between infection and immunity suggest that newly founded populations can exhibit rapid host-parasite eco-evolutionary dynamics.

Designing eco-evolutionary experiments for restoration projects: Opportunities and constraints revealed during stickleback introductions.

Eco-evolutionary experiments are typically conducted in semi-unnatural controlled settings, such as mesocosms; yet inferences about how evolution and ecology interact in the real world would surely benefit from experiments in natural uncontrolled settings. Opportunities for such experiments are rare but do arise in the context of restoration ecology-where different "types" of a given species can be introduced into different "replicate" locations. Designing such experiments requires wrestling with consequential questions. (Q1) Which specific "types" of a focal species should be introduced to the restoration location? (Q2) How many sources of each type should be used-and should they be mixed together? (Q3) Which specific source populations should be used? (Q4) Which type(s) or population(s) should be introduced into which restoration sites? We recently grappled with these questions when designing an eco-evolutionary experiment with threespine stickleback (Gasterosteus aculeatus) introduced into nine small lakes and ponds on the Kenai Peninsula in Alaska that required restoration. After considering the options at length, we decided to use benthic versus limnetic ecotypes (Q1) to create a mixed group of colonists from four source populations of each ecotype (Q2), where ecotypes were identified based on trophic morphology (Q3), and were then introduced into nine restoration lakes scaled by lake size (Q4). We hope that outlining the alternatives and resulting choices will make the rationales clear for future studies leveraging our experiment, while also proving useful for investigators considering similar experiments in the future.

A spectral framework to map QTLs affecting joint differential networks of gene co-expression.

Studying the mechanisms underlying the genotype-phenotype association is crucial in genetics. Gene expression studies have deepened our understanding of the genotype → expression → phenotype mechanisms. However, traditional expression quantitative trait loci (eQTL) methods often overlook the critical role of gene co-expression networks in translating genotype into phenotype. This gap highlights the need for more powerful statistical methods to analyze genotype → network → phenotype mechanism. Here, we develop a network-based method, called snQTL, to map quantitative trait loci affecting gene co-expression networks. Our approach tests the association between genotypes and joint differential networks of gene co-expression via a tensor-based spectral statistics, thereby overcoming the ubiquitous multiple testing challenges in existing methods. We demonstrate the effectiveness of snQTL in the analysis of three-spined stickleback (Gasterosteus aculeatus) data. Compared to conventional methods, our method snQTL uncovers chromosomal regions affecting gene co-expression networks, including one strong candidate gene that would have been missed by traditional eQTL analyses. Our framework suggests the limitation of current approaches and offers a powerful network-based tool for functional loci discoveries.

Estimating the magnitude and sensitivity of energy fluxes for stickleback hosts and Schistocephalus solidus parasites using the metabolic theory of ecology.

Parasites are ubiquitous, yet their effects on hosts are difficult to quantify and generalize across ecosystems. One promising metric of parasitic impact uses the metabolic theory of ecology (MTE) to calculate energy flux, an estimate of energy lost to parasites. We investigated the feasibility of using metabolic scaling rules to compare the energetic burden of parasitism among individuals. Specifically, we found substantial sensitivity of energy flux estimates to input parameters used in the MTE equation when using available data from a model host-parasite system (Gasterosteus aculeatus and Schistocephalus solidus). Using literature values, size data from parasitized wild fish, and a respirometry experiment, we estimate that a single S. solidus tapeworm may extract up to 32% of its stickleback host's baseline metabolic energy requirement, and that parasites in multiple infections may collectively extract up to 46%. The amount of energy siphoned from stickleback to tapeworms is large but did not instigate an increase in respiration rate in the current study. This emphasizes the importance of future work focusing on how parasites influence ecosystem energetics. The approach of using the MTE to calculate energy flux provides great promise as a quantitative foundation for such estimates and provides a more concrete metric of parasite impact on hosts than parasite abundance alone.

Evolution of horn length and lifting strength in the Japanese rhinoceros beetle Trypoxylus dichotomus.

What limits the size of nature's most extreme structures? For weapons like beetle horns, one possibility is a tradeoff associated with mechanical levers: as the output arm of the lever system-the beetle horn-gets longer, it also gets weaker. This "paradox of the weakening combatant" could offset reproductive advantages of additional increases in weapon size. However, in contemporary populations of most heavily weaponed species, males with the longest weapons also tend to be the strongest, presumably because selection drove the evolution of compensatory changes to these lever systems that ameliorated the force reductions of increased weapon size. Therefore, we test for biomechanical limits by reconstructing the stages of weapon evolution, exploring whether initial increases in weapon length first led to reductions in weapon force generation that were later ameliorated through the evolution of mechanisms of mechanical compensation. We describe phylogeographic relationships among populations of a rhinoceros beetle and show that the "pitchfork" shaped head horn likely increased in length independently in the northern and southern radiations of beetles. Both increases in horn length were associated with dramatic reductions to horn lifting strength-compelling evidence for the paradox of the weakening combatant-and these initial reductions to horn strength were later ameliorated in some populations through reductions to horn length or through increases in head height (the input arm for the horn lever system). Our results reveal an exciting geographic mosaic of weapon size, weapon force, and mechanical compensation, shedding light on larger questions pertaining to the evolution of extreme structures.

A Call For More Ecologically And Evolutionarily Relevant Studies of Immune Costs.

What are the relative costs and benefits of mounting immune responses? Practitioners of ecoimmunology have grappled with this central question since the field's inception with the main tension being how to make tractable methodological choices that maintain the ecological relevance of induced and measured immune costs. Here, we point out two methodological approaches that we feel are underrepresented in the field, describe risks associated with neglecting these methods, and suggest modern techniques that maximize both the diversity and ecological relevance of collected data. First, it is commonly assumed that frequently used and experimentally convenient immune stimulants will induce ecologically relevant immune responses in study organisms. This can be a dangerous assumption. Even if a stimulant's general immune response properties are well characterized, it is critical to also measure the type and scale of immune responses induced by live pathogens. Second, patterns of immune defenses evolve like other traits, thus a comparative approach is essential to understand what forces shape immune variation. Finally, we describe modern genetic and immunological approaches that will soon become essential tools for ecoimmunologists, and present case studies that exemplify the utility of our recommendations.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation.

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

Evolutionary gain and loss of a pathological immune response to parasitism.

Parasites impose fitness costs on their hosts. Biologists often assume that natural selection favors infection-resistant hosts. Yet, when the immune response itself is costly, theory suggests that selection may sometimes favor loss of resistance, which may result in alternative stable states where some populations are resistant and others are tolerant. Intraspecific variation in immune costs is rarely surveyed in a manner that tests evolutionary patterns, and there are few examples of adaptive loss of resistance. Here, we show that when marine threespine stickleback colonized freshwater lakes, they gained resistance to the freshwater-associated cestode Schistocephalus solidus. Extensive peritoneal fibrosis and inflammation are a commonly observed phenotype that contributes to suppression of cestode growth and viability but also imposes a substantial cost on fecundity. Combining genetic mapping and population genomics, we find that opposing selection generates immune system differences between tolerant and resistant populations, consistent with divergent optimization.

Interspecific variation in cooperative burrowing behavior by Peromyscus mice.

Animals often adjust their behavior according to social context, but the capacity for such behavioral flexibility can vary among species. Here, we test for interspecific variation in behavioral flexibility by comparing burrowing behavior across three species of deer mice (genus Peromyscus) with divergent social systems, ranging from promiscuous (Peromyscus leucopus and Peromyscus maniculatus) to monogamous (Peromyscus polionotus). First, we compared the burrows built by individual mice to those built by pairs of mice in all three species. Although burrow length did not differ in P. leucopus or P. maniculatus, we found that P. polionotus pairs cooperatively constructed burrows that were nearly twice as long as those built by individuals and that opposite-sex pairs dug longer burrows than same-sex pairs. Second, to directly observe cooperative digging behavior in P. polionotus, we designed a burrowing assay in which we could video-record active digging in narrow, transparent enclosures. Using this novel assay, we found, unexpectedly, that neither males nor females spent more time digging with an opposite-sex partner. Rather, we demonstrate that opposite-sex pairs are more socially cohesive and thus more efficient digging partners than same-sex pairs. Together, our study demonstrates how social context can modulate innate behavior and offers insight into how differences in behavioral flexibility may evolve among closely related species.

Between-population differences in constitutive and infection-induced gene expression in threespine stickleback.

Vertebrate immunity is a complex system consisting of a mix of constitutive and inducible defences. Furthermore, host immunity is subject to selective pressure from a range of parasites and pathogens which can produce variation in these defences across populations. As populations evolve immune responses to parasites, they may adapt via a combination of (1) constitutive differences, (2) shared inducible responses, or (3) divergent inducible responses. Here, we leverage a powerful natural host-parasite model system (Gasterosteus aculeatus and Schistochephalus solidus) to tease apart the relative contributions of these three types of adaptations to among-population divergence in response to parasites. Gene expression analyses revealed limited evidence of significant divergence in constitutive expression of immune defence, and strong signatures of conserved inducible responses to the parasite. Furthermore, our results highlight a handful of immune-related genes which show divergent inducible responses which may contribute disproportionately to functional differences in infection success or failure. In addition to investigating variation in evolutionary adaptation to parasite selection, we also leverage this unique data set to improve understanding of cellular mechanisms underlying a putative resistance phenotype (fibrosis). Combined, our results provide a case study in evolutionary immunology showing that a very small number of genes may contribute to genotype differences in infection response.

Immune Gene Expression Covaries with Gut Microbiome Composition in Stickleback.

Commensal microbial communities have immense effects on their vertebrate hosts, contributing to a number of physiological functions, as well as host fitness. In particular, host immunity is strongly linked to microbiota composition through poorly understood bi-directional links. Gene expression may be a potential mediator of these links between microbial communities and host function. However, few studies have investigated connections between microbiota composition and expression of host immune genes in complex systems. Here, we leverage a large study of laboratory-raised fish from the species Gasterosteus aculeatus (three-spined stickleback) to document correlations between gene expression and microbiome composition. First, we examined correlations between microbiome alpha diversity and gene expression. Our results demonstrate robust positive associations between microbial alpha diversity and expression of host immune genes. Next, we examined correlations between host gene expression and abundance of microbial taxa. We identified 15 microbial families that were highly correlated with host gene expression. These families were all tightly correlated with host expression of immune genes and processes, falling into one of three categories-those positively correlated, negatively correlated, and neutrally related to immune processes. Furthermore, we highlight several important immune processes that are commonly associated with the abundance of these taxa, including both macrophage and B cell functions. Further functional characterization of microbial taxa will help disentangle the mechanisms of the correlations described here. In sum, our study supports prevailing hypotheses of intimate links between host immunity and gut microbiome composition.IMPORTANCE Here, we document associations between host gene expression and gut microbiome composition in a nonmammalian vertebrate species. We highlight associations between expression of immune genes and both microbiome diversity and abundance of specific microbial taxa. These findings support other findings from model systems which have suggested that gut microbiome composition and host immunity are intimately linked. Furthermore, we demonstrate that these correlations are truly systemic; the gene expression detailed here was collected from an important fish immune organ (the head kidney) that is anatomically distant from the gut. This emphasizes the systemic impact of connections between gut microbiota and host immune function. Our work is a significant advancement in the understanding of immune-microbiome links in nonmodel, natural systems.

Landscape genomics of an obligate mutualism: Concordant and discordant population structures between the leafcutter ant Atta texana and its two main fungal symbiont types.

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of host-symbiont interactions between the northernmost fungus-farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome-wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant-fungus genome-genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population-genetic structure was concordant between the ants and one cultivar type (M-fungi, concordant clines) but discordant for the other cultivar type (T-fungi). Discordance in population-genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between-nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant-fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna.

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present-day tropical biota. Using genome-wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post-Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as "isolation by instability," highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.

Gene Expression Contributes to the Recent Evolution of Host Resistance in a Model Host Parasite System.

Heritable population differences in immune gene expression following infection can reveal mechanisms of host immune evolution. We compared gene expression in infected and uninfected threespine stickleback (Gasterosteus aculeatus) from two natural populations that differ in resistance to a native cestode parasite, Schistocephalus solidus. Genes in both the innate and adaptive immune system were differentially expressed as a function of host population, infection status, and their interaction. These genes were enriched for loci controlling immune functions known to differ between host populations or in response to infection. Coexpression network analysis identified two distinct processes contributing to resistance: parasite survival and suppression of growth. Comparing networks between populations showed resistant fish have a dynamic expression profile while susceptible fish are static. In summary, recent evolutionary divergence between two vertebrate populations has generated population-specific gene expression responses to parasite infection, affecting parasite establishment and growth.

Contrasting effects of environment and genetics generate a continuum of parallel evolution.

Parallel evolution of similar traits by independent populations in similar environments is considered strong evidence for adaptation by natural selection. Often, however, replicate populations in similar environments do not all evolve in the same way, thus deviating from any single, predominant outcome of evolution. This variation might arise from non-adaptive, population-specific effects of genetic drift, gene flow or limited genetic variation. Alternatively, these deviations from parallel evolution might also reflect predictable adaptation to cryptic environmental heterogeneity within discrete habitat categories. Here, we show that deviations from parallel evolution are the consequence of environmental variation within habitats combined with variation in gene flow. Threespine stickleback (Gasterosteus aculeatus) in adjoining lake and stream habitats (a lake-stream 'pair') diverge phenotypically, yet the direction and magnitude of this divergence is not always fully parallel among 16 replicate pairs. We found that the multivariate direction of lake-stream morphological divergence was less parallel between pairs whose environmental differences were less parallel. Thus, environmental heterogeneity among lake-stream pairs contributes to deviations from parallel evolution. Additionally, likely genomic targets of selection were more parallel between environmentally more similar pairs. In contrast, variation in the magnitude of lake-stream divergence (independent of direction) was better explained by differences in lake-stream gene flow; pairs with greater lake-stream gene flow were less morphologically diverged. Thus, both adaptive and non-adaptive processes work concurrently to generate a continuum of parallel evolution across lake-stream stickleback population pairs.

Recent evolution of extreme cestode growth suppression by a vertebrate host.

Parasites can be a major cause of natural selection on hosts, which consequently evolve a variety of strategies to avoid, eliminate, or tolerate infection. When ecologically similar host populations present disparate infection loads, this natural variation can reveal immunological strategies underlying adaptation to infection and population divergence. For instance, the tapeworm Schistocephalus solidus persistently infects 0-80% of threespine stickleback (Gasterosteus aculeatus) in lakes on Vancouver Island. To test whether these heterogeneous infection rates result from evolved differences in immunity, we experimentally exposed laboratory-reared fish from ecologically similar high-infection and no-infection populations to controlled doses of Schistocephalus We observed heritable between-population differences in several immune traits: Fish from the naturally uninfected population initiated a stronger granulocyte response to Schistocephalus infection, and their granulocytes constitutively generate threefold more reactive oxygen species in cell culture. Despite these immunological differences, Schistocephalus was equally successful at establishing initial infections in both host populations. However, the no-infection fish dramatically suppressed tapeworm growth relative to high-infection fish, and parasite size was intermediate in F1 hybrid hosts. Our results show that stickleback recently evolved heritable variation in their capacity to suppress helminth growth by two orders of magnitude. Data from many natural populations indicate that growth suppression is widespread but not universal and, when present, is associated with reduced infection prevalence. Host suppression of helminth somatic growth may be an important immune strategy that aids in parasite clearance or in mitigating the fitness costs of persistent infection.

Resist Globally, Infect Locally: A Transcontinental Test of Adaptation by Stickleback and Their Tapeworm Parasite.

Parasite infections are a product of both ecological processes affecting host-parasite encounter rates and evolutionary dynamics affecting host susceptibility. However, few studies examine natural infection variation from both ecological and evolutionary perspectives. Here, we describe the ecological and evolutionary factors generating variation in infection rates by a tapeworm (Schistocephalus solidus) in a vertebrate host, the threespine stickleback (Gasterosteus aculeatus). To explore ecological aspects of infection, we measured tapeworm prevalence in Canadian stickleback inhabiting two distinct environments: marine and freshwater. Consistent with ecological control of infection, the tapeworm is very rare in marine environments, even though marine fish are highly susceptible. Conversely, commonly infected freshwater stickleback exhibit substantial resistance in controlled laboratory trials, suggesting that high exposure risk overwhelms their recently evolved resistance. We also tested for parasite adaptation to its host by performing transcontinental reciprocal infections, using stickleback and tapeworm populations from Europe and western Canada. More infections occurred in same-continent host-parasite combinations, indicating parasite "local" adaptation, at least on the scale of continents. However, the recently evolved immunity of freshwater hosts applies to both local and foreign parasites. The pattern of adaptation described here is not wholly compatible with either of the common models of host-parasite coevolution (i.e., matching infection or targeted recognition). Instead, we propose a hybrid, eco-evolutionary model to explain the remarkable pattern of global host resistance and local parasite infectivity.

Partitioning the effects of isolation by distance, environment, and physical barriers on genomic divergence between parapatric threespine stickleback.

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation-by-distance (IBD) and isolation-by-environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake-stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake-stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake-stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.

Evaluation of TagSeq, a reliable low-cost alternative for RNAseq.

RNAseq is a relatively new tool for ecological genetics that offers researchers insight into changes in gene expression in response to a myriad of natural or experimental conditions. However, standard RNAseq methods (e.g., Illumina TruSeq® or NEBNext® ) can be cost prohibitive, especially when study designs require large sample sizes. Consequently, RNAseq is often underused as a method, or is applied to small sample sizes that confer poor statistical power. Low cost RNAseq methods could therefore enable far greater and more powerful applications of transcriptomics in ecological genetics and beyond. Standard mRNAseq is costly partly because one sequences portions of the full length of all transcripts. Such whole-mRNA data are redundant for estimates of relative gene expression. TagSeq is an alternative method that focuses sequencing effort on mRNAs' 3' end, reducing the necessary sequencing depth per sample, and thus cost. We present a revised TagSeq library construction procedure, and compare its performance against NEBNext® , the 'gold-standard' whole mRNAseq method. We built both TagSeq and NEBNext® libraries from the same biological samples, each spiked with control RNAs. We found that TagSeq measured the control RNA distribution more accurately than NEBNext® , for a fraction of the cost per sample (~10%). The higher accuracy of TagSeq was particularly apparent for transcripts of moderate to low abundance. Technical replicates of TagSeq libraries are highly correlated, and were correlated with NEBNext® results. Overall, we show that our modified TagSeq protocol is an efficient alternative to traditional whole mRNAseq, offering researchers comparable data at greatly reduced cost.