Optimal outbreeding is shaped during larval life-history in the splash pool copepod Tigriopus californicus.

Inbreeding and outbreeding depression are dynamic forms of selection critical to mating system evolution and the efficacy of conservation biology. Most evidence on how the relative severity and timing of these forces are shaped is confined to self-fertilization, distant outcrossing, and intermediate 'optimal outcrossing' in hermaphrodites. We tested the notion that closed population demographics may reduce and delay the costs of inbreeding relative to distant outbreeding in an intertidal copepod with separate sexes and a biphasic larval / post-metamorphic life-history (Tigriopus californicus). At three lifecycle stages (fecundity, metamorphosis, and post-metamorphosis), we quantified the effects of inbreeding and outbreeding in crosses with varying degrees of recent common ancestry. Although inbreeding and outbreeding depression have distinct genetic mechanisms, both manifested the same stage-specific consequences for fitness. Inbreeding and outbreeding depression were not apparent for fecundity, post-metamorphic survival, sex ratio, or the ability to acquire mates, but inbreeding between full siblings and outbreeding between interpopulation hybrids reduced the fraction of offspring that completed metamorphosis by 32% and 47%, respectively. On average, the effects of inbreeding on metamorphic rate were weaker and nearly twice as variable among families than those of outbreeding, suggesting genetic load was less pervasive than the incompatibilities accrued between divergent populations. Overall, our results indicate the transition from larval to juvenile life stages is markedly susceptible to both inbreeding and outbreeding depression in T. californicus. We suggest stage-specific selection acting concurrently with the timing of metamorphosis may be an instrumental factor shaping reproductive optima in species with complex life-histories.

A hybrid beachgrass (Ammophila arenaria × A. breviligulata) is more productive and outcompetes its non-native parent species.

The ability of non-native species to successfully invade new ecosystems sometimes involves evolutionary processes such as hybridization. Hybridization can produce individuals with superior traits that give them a competitive advantage over their parent species, allowing for rapid spread. Here we assess growth, functional morphology, and species interactions between two non-native beachgrass species (Ammophila arenaria and A. breviligulata) and their recently discovered hybrid (A. arenaria × A. breviligulata) on the U.S. Pacific Northwest coast. We asked whether the hybrid beachgrass differs from its parent species in morphology and growth, whether it competes with its parent species, and, if so, what are the potential mechanisms of competition. Plant taxa were grown in low- and high-density monocultures and in two-way interactions in a common garden environment. We show that the hybrid grew taller and more densely, with greater total biomass, than either parent species. The hybrid was also the better competitor, resulting in the model prediction of competitive exclusion against A. breviligulata and, depending on its relative abundance, A. arenaria. The hybrid displays a mixed 'guerilla-phalanx' growth form that allows it to spread laterally and achieve high shoot densities, giving it a competitive advantage. Given the current dominance of A. breviligulata compared to A. arenaria in most of the region where these taxa co-occur, we suggest that the hybrid will grow, compete, and spread quickly with potentially widespread consequences for the two non-native Ammophila congeners and the dunes they build.

Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod.

Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.

Genomic Architecture of Hybrid Male Sterility in a Species Without Sex Chromosomes (Tigriopus californicus, Copepoda: Harpacticoida).

Sterility among hybrids is one of the most prevalent forms of reproductive isolation delineating species boundaries and is expressed disproportionately in heterogametic XY males. While hybrid male sterility (HMS) due to the "large X effect" is a well-recognized mechanism of reproductive isolation, it is less clear how HMS manifests in species that lack heteromorphic sex chromosomes. We evaluated differences in allele frequencies at approximately 460,000 SNPs between fertile and sterile F2 interpopulation male hybrids to characterize the genomic architecture of HMS in a species without sex chromosomes (Tigriopus californicus). We tested associations between HMS and mitochondrial-nuclear and/or nuclear-nuclear signatures of incompatibility. Genomic regions associated with HMS were concentrated on a single chromosome with the same primary 2-Mbp regions identified in one pair of reciprocal crosses. Gene Ontology analysis revealed that annotations associated with spermatogenesis were the most overrepresented within the implicated region, with nine protein-coding genes connected with this process found in the quantitative trait locus of chromosome 2. Our results indicate that a narrow genomic region was associated with the sterility of male hybrids in T. californicus and suggest that incompatibilities among select nuclear loci may replace the large X effect when sex chromosomes are absent.

Myxozoans (Cnidaria) do not Retain Key Oxygen-Sensing and Homeostasis Toolkit Genes.

For aerobic organisms, both the hypoxia-inducible factor pathway and the mitochondrial genomes are key players in regulating oxygen homeostasis. Recent work has suggested that these mechanisms are not as highly conserved as previously thought, prompting more surveys across animal taxonomic levels, which would permit testing of hypotheses about the ecological conditions facilitating evolutionary loss of such genes. The Phylum Cnidaria is known to harbor wide variation in mitochondrial chromosome morphology, including an extreme example, in the Myxozoa, of mitochondrial genome loss. Because myxozoans are obligate endoparasites, frequently encountering hypoxic environments, we hypothesize that variation in environmental oxygen availability could be a key determinant in the evolution of metabolic gene networks associated with oxygen-sensing, hypoxia-response, and energy production. Here, we surveyed genomes and transcriptomes across 46 cnidarian species for the presence of HIF pathway members, as well as for an assortment of hypoxia, mitochondrial, and stress-response toolkit genes. We find that presence of the HIF pathway, as well as number of genes associated with mitochondria, hypoxia, and stress response, do not vary in parallel to mitochondrial genome morphology. More interestingly, we uncover evidence that myxozoans have lost the canonical HIF pathway repression machinery, potentially altering HIF pathway functionality to work under the specific conditions of their parasitic lifestyles. In addition, relative to other cnidarians, myxozoans show loss of large proportions of genes associated with the mitochondrion and involved in response to hypoxia and general stress. Our results provide additional evidence that the HIF regulatory machinery is evolutionarily labile and that variations in the canonical system have evolved in many animal groups.

Little evidence for genetic variation associated with susceptibility to sea star wasting syndrome in the keystone species Pisaster ochraceus.

The keystone species Pisaster ochraceus suffered mass mortalities along the northeast Pacific Ocean from Sea Star Wasting Syndrome (SSWS) outbreaks in 2013-2016. SSWS causation remains of debate, leading to concerns as to whether outbreaks will continue to impact this species. Considering the apparent link between ocean temperature and SSWS, the future of this species and intertidal communities remains uncertain. Surveys of co-occurring apparently normal and wasting P. ochraceus along the central Oregon coast in 2016 allowed us to address whether variation in disease status showed genetic variation that may be associated with differences in susceptibility to SSWS. We performed restriction site-associated DNA sequencing (2bRAD-seq) to genotype ~72,000 single nucleotide polymorphism (SNP) loci across apparently normal and wasting sea stars. Locus-specific analyses of differentiation (FST ) between disease-status groups revealed no signal of genetic differences separating the two groups. Using a multivariate approach, we observed weak separation between the groups, but identified 18 SNP loci showing highest discriminatory power between the groups and scanned the genome annotation for linked genes. A total of 34 protein-coding genes were found to be located within 15 kb (measured by linkage disequilibrium decay) of at least one of the 18 SNPs, and 30 of these genes had homologies to annotated protein databases. Our results suggest that the likelihood of developing SSWS symptoms does not have a strong genetic basis. The few genomic regions highlighted had only modest levels of differentiation, but the genes associated with these regions may form the basis for functional studies aiming to understand disease progression.

Pervasive Mitonuclear Coadaptation Underlies Fast Development in Interpopulation Hybrids of a Marine Crustacean.

Cellular energy production requires coordinated interactions between genetic components from the nuclear and mitochondrial genomes. This coordination results in coadaptation of interacting elements within populations. Interbreeding between divergent gene pools can disrupt coadapted loci and result in hybrid fitness breakdown. While specific incompatible loci have been detected in multiple eukaryotic taxa, the extent of the nuclear genome that is influenced by mitonuclear coadaptation is not clear in any species. Here, we used F2 hybrids between two divergent populations of the copepod Tigriopus californicus to examine mitonuclear coadaptation across the nuclear genome. Using developmental rate as a measure of fitness, we found that fast-developing copepods had higher ATP synthesis capacity than slow developers, suggesting variation in developmental rates is at least partly associated with mitochondrial dysfunction. Using Pool-seq, we detected strong biases for maternal alleles across 7 (of 12) chromosomes in both reciprocal crosses in high-fitness hybrids, whereas low-fitness hybrids showed shifts toward the paternal population. Comparison with previous results on a different hybrid cross revealed largely different patterns of strong mitonuclear coadaptation associated with developmental rate. Our findings suggest that functional coadaptation between interacting nuclear and mitochondrial components is reflected in strong polygenic effects on this life-history phenotype, and reveal that molecular coadaptation follows independent evolutionary trajectories among isolated populations.

Phenotypic Variation in Growth and Gene Expression Under Different Photoperiods in Allopatric Populations of the Copepod Tigriopus californicus.

Daylength is a major environmental condition that varies seasonally and predictably along a latitudinal cline, where higher latitudes exhibit greater ranges in total daylengths. Generally, the circadian clock acts as a network of genes whose expression dynamics are known to control daily rhythms in response to daylength, and it enables the control of many physiological processes such as growth and development. While well studied in many model animals, the influence of daylength variation on phenotypic evolution is poorly examined in marine species. In this study we demonstrate that two allopatric populations of the intertidal crustacean Tigriopus californicus exhibit plastic and divergent phenotypic responses to changes in daylength. Using common-garden experiments, we discovered that shorter daylengths promoted decreased adult body size and faster growth rates in the two divergent populations, suggesting a plastic response to shortened days. In addition, the higher-latitude population exhibited a faster growth rate at any daylength condition, indicating a fixed response, possibly as a result of adaptation to respective natural light regimes. Gene expression profiles of several circadian clock genes, monitored throughout the day by quantitative polymerase chain reaction, revealed that the key core clock genes reach higher daily transcription maxima in the southern population compared to the northern population, pointing to divergent strategies used to respond to changes in daylength. Many modifier genes to the circadian clock showed similar plastic responses to the different daylengths, supporting the existence of at least some conserved gene expression across both populations. Ultimately, our results suggest that photoperiod and daylength exert a potent selective pressure underexplored in marine systems and warranting further future research.

Independent Losses of the Hypoxia-Inducible Factor (HIF) Pathway within Crustacea.

Metazoans respond to hypoxic stress via the hypoxia-inducible factor (HIF) pathway, a mechanism thought to be extremely conserved due to its importance in monitoring cellular oxygen levels and regulating responses to hypoxia. However, recent work revealed that key members of the HIF pathway have been lost in specific lineages (a tardigrade and a copepod), suggesting that this pathway is not as widespread in animals as previously assumed. Using genomic and transcriptomic data from 70 different species across 12 major crustacean groups, we assessed the degree to which the gene HIFα, the master regulator of the HIF pathway, was conserved. Mining of protein domains, followed by phylogenetic analyses of gene families, uncovered group-level losses of HIFα, including one across three orders within Cirripedia, and in three orders within Copepoda. For these groups, additional assessment showed losses of HIF repression machinery (EGLN and VHL). These results suggest the existence of alternative mechanisms for cellular response to low oxygen and highlight these taxa as models useful for probing these evolutionary outcomes.

Loss of the HIF pathway in a widely distributed intertidal crustacean, the copepod Tigriopus californicus.

Hypoxia is a major physiological constraint for which multicellular eukaryotes have evolved robust cellular mechanisms capable of addressing dynamic changes in O2 availability. In animals, oxygen sensing and regulation is primarily performed by the hypoxia-inducible factor (HIF) pathway, and the key components of this pathway are thought to be highly conserved across metazoans. Marine intertidal habitats are dynamic environments, and their inhabitants are known to tolerate wide fluctuations in salinity, temperature, pH, and oxygen. In this study, we show that an abundant intertidal crustacean, the copepod Tigriopus californicus, has lost major genetic components of the HIF pathway, but still shows robust survivorship and transcriptional response to hypoxia. Mining of protein domains across the genome, followed by phylogenetic analyses of gene families, did not identify two key regulatory elements of the metazoan hypoxia response, namely the transcription factor HIF-α and its oxygen-sensing prolyl hydroxylase repressor, EGLN Despite this loss, phenotypic assays revealed that this species is tolerant to extremely low levels of available O2 for at least 24 h at both larval and adult stages. RNA-sequencing (seq) of copepods exposed to nearly anoxic conditions showed differential expression of over 400 genes, with evidence for induction of glycolytic metabolism without a depression of oxidative phosphorylation. Moreover, genes involved in chitin metabolism and cuticle reorganization show categorically a consistent pattern of change during anoxia, highlighting this pathway as a potential solution to low oxygen availability in this small animal with no respiratory structures or pigment.

Rapid metabolic compensation in response to temperature change in the intertidal copepod, Tigriopus californicus.

Animals living in the intertidal zone must adapt to thermal variability, including adjustments in metabolism. We examined metabolic responses to temperature in the copepod, Tigriopus californicus, which inhabits supratidal splash pools along the Pacific coast of North America. We maintained three populations of T. californicus at 20 °C, one from southern California (San Diego, "SD") and two from Oregon (Fogarty Creek, "FCN", Boiler Bay, "BOB") and examined possible geographic patterns in metabolism. We measured oxygen consumption rate (V̇o2) at 20 °C and following 48 h (chronic) acclimation to 25, 27.5 and 30 °C. V̇o2 was temperature-independent, with temperature quotients (Q10) values ≤1 in all populations, indicative of metabolic compensation. We detected no variation in V̇o2 or survival between populations. To explore the time course of metabolic compensation, we performed an acute acclimation experiment in which V̇o2 was measured at 20 °C, following immediate exposure to 25 °C, and following 2 h, 4 h and 6 h exposure to 25 °C. This acute acclimation experiment revealed that V̇o2 increased immediately in SD and FCN, but was no longer different than 20 °C levels by 2 h and 4 h at 25 °C, respectively. BOB showed no significant change in V̇o2, which may indicate complete temperature-independent metabolism or different mechanisms of compensation between populations. This study demonstrates a time course of rapid metabolic compensation in response to temperature that occurs in a small intertidal animal, and suggests intertidal invertebrates can thermally acclimate within a few hours of a significant temperature change.

Novel microRNAs are associated with population divergence in transcriptional response to thermal stress in an intertidal copepod.

The role of gene expression in adaptation to differing thermal environments has been assayed extensively. Yet, in most natural systems, analyses of gene expression reveal only one level of the complexity of regulatory machineries. MicroRNAs (miRNAs) are small noncoding RNAs which are key components of many gene regulatory networks, and they play important roles in a variety of cellular pathways by modulating post-transcriptional quantities of mRNA available for protein synthesis. The characterization of miRNA loci and their regulatory dynamics in nonmodel systems are still largely understudied. In this study, we examine the role of miRNAs in response to high thermal stress in the intertidal copepod Tigriopus californicus. Allopatric populations of this species show varying levels of local adaptation with respect to thermal regimes, and previous studies showed divergence in gene expression between populations from very different thermal environments. We examined the transcriptional response to temperature stress in two populations separated by only 8 km by utilizing RNA-seq to quantify both mRNA and miRNA levels. Using the currently available genome sequence, we first describe the repertoire of miRNAs in T. californicus and assess the degree to which transcriptional response to temperature stress is governed by miRNA activity. The two populations showed large differences in the number of genes involved, the magnitude of change in commonly used genes and in the number of miRNAs involved in transcriptional modulation during stress. Our results suggest that an increased level of regulatory network complexity may underlie improved survivorship under thermal stress in one of the populations.

Genomic signatures of mitonuclear coevolution across populations of Tigriopus californicus.

The copepod Tigriopus californicus shows extensive population divergence and is becoming a model for understanding allopatric differentiation and the early stages of speciation. Here, we report a high-quality reference genome for one population (~190 megabases across 12 scaffolds, and ~15,500 protein-coding genes). Comparison with other arthropods reveals 2,526 genes presumed to be specific to T. californicus, with an apparent proliferation of genes involved in ion transport and receptor activity. Beyond the reference population, we report re-sequenced genomes of seven additional populations, spanning the continuum of reproductive isolation. Populations show extreme mitochondrial DNA divergence, with higher levels of amino acid differentiation than observed in other taxa. Across the nuclear genome, we find elevated protein evolutionary rates and positive selection in genes predicted to interact with mitochondrial DNA and the proteins and RNA it encodes in multiple pathways. Together, these results support the hypothesis that rapid mitochondrial evolution drives compensatory nuclear evolution within isolated populations, thereby providing a potentially important mechanism for causing intrinsic reproductive isolation.

Transcriptome-wide patterns of divergence during allopatric evolution.

Recent studies have revealed repeated patterns of genomic divergence associated with species formation. Such patterns suggest that natural selection tends to target a set of available genes, but is also indicative that closely related taxa share evolutionary constraints that limit genetic variability. Studying patterns of genomic divergence among populations within the same species may shed light on the underlying evolutionary processes. Here, we examine transcriptome-wide divergence and polymorphism in the marine copepod Tigriopus californicus, a species where allopatric evolution has led to replicate sets of populations with varying degrees of divergence and hybrid incompatibility. Our analyses suggest that relatively small effective population sizes have resulted in an exponential decline of shared polymorphisms during population divergence and also facilitated the fixation of slightly deleterious mutations within allopatric populations. Five interpopulation comparisons at three different stages of divergence show that nonsynonymous mutations tend to accumulate in a specific set of proteins. These include proteins with central roles in cellular metabolism, such as those encoded in mtDNA, but also include an additional set of proteins that repeatedly show signatures of positive selection during allopatric divergence. Although our results are consistent with a contribution of nonadaptive processes, such as genetic drift and gene expression levels, generating repeatable patterns of genomic divergence in closely related taxa, they also indicate that adaptive evolution targeting a specific set of genes contributes to this pattern. Our results yield insights into the predictability of evolution at the gene level.

Reverse genetics in the tide pool: knock-down of target gene expression via RNA interference in the copepod Tigriopus californicus.

Reverse genetic tools are essential for characterizing phenotypes of novel genes and testing functional hypotheses generated from next-generation sequencing studies. RNA interference (RNAi) has been a widely used technique for describing or quantifying physiological, developmental or behavioural roles of target genes by suppressing their expression. The marine intertidal copepod Tigriopus californicus has become an emerging model for evolutionary and physiological studies, but this species is not amenable to most genetic manipulation approaches. As crustaceans are susceptible to RNAi-mediated gene knock-down, we developed a simple method for delivery of gene-specific double-stranded RNA that results in significant suppression of target gene transcription levels. The protocol was examined on five genes of interest, and for each, at least 50% knock-down in expression was achieved. While knock-down levels did not reach 100% in any trial, a well-controlled experiment with one heat-shock gene showed unambiguously that such partial gene suppression may cause dramatic changes in phenotype. Copepods with suppressed expression of heat-shock protein beta 1 (hspb1) exhibited dramatically decreased tolerance to high temperatures, validating the importance of this gene during thermal stress, as proposed by a previous study. The application of this RNAi protocol in T. californicus will be invaluable for examining the role of genes putatively involved in reproductive isolation, mitochondrial function and local adaptation.

Hybrid dysfunction and physiological compensation in gene expression.

The formation of new species is often a consequence of genetic incompatibilities accumulated between populations during allopatric divergence. When divergent taxa interbreed, these incompatibilities impact physiology and have a direct cost resulting in reduced hybrid fitness. Recent surveys of gene regulation in interspecific hybrids have revealed anomalous expression across large proportions of the genome, with 30-70% of all genes exhibiting transgressive expression (i.e., higher or lower levels compared with both parental taxa), and these were mostly in the direction of downregulation. However, as most of these studies have focused on pairs of species exhibiting high degrees of reproductive isolation, the association between regulatory disruption and reduced hybrid fitness prior to species formation remains unclear. Within the copepod species Tigriopus californicus, interpopulation hybrids at F2 or later generations show reduced fitness associated with mitochondrial dysfunction. Here we show that in contrast to studies of interspecific hybrids, only 1.2% of the transcriptome is transgressively expressed in F3+ interpopulation hybrids of T. californicus, and nearly 80% of these genes are overexpressed rather than underexpressed; remarkably, none of these genes are among those showing divergent expression between parentals, nor is magnitude of transgressive gene expression in hybrids dependent on levels of protein sequence divergence. Moreover, many genes with transgressive expression are components of functional pathways impacted by mitonuclear incompatibilities in hybrid T. californicus (e.g., oxidative phosphorylation and antioxidant response). Our results suggest that hybrid breakdown at early stages of speciation may result from initial incompatibilities amplified by the cost of compensatory physiological responses.

Ecological novelty by hybridization: experimental evidence for increased thermal tolerance by transgressive segregation in Tigriopus californicus.

Early generations of hybrids can express both genetic incompatibilities and phenotypic novelty. Insights into whether these conflicting interactions between intrinsic and extrinsic selection persist after a few generations of recombination require experimental studies. To address this question, we use interpopulation crosses and recombinant inbred lines (RILs) of the copepod Tigriopus californicus, and focus on two traits that are relevant for the diversification of this species: survivorship during development and tolerance to thermal stress. Experimental crosses between two population pairs show that most RILs between two heat-tolerant populations show enhanced tolerance to temperatures that are lethal to the respective parentals, whereas RILs between a heat-tolerant and a heat-sensitive population are intermediate. Although interpopulation crosses are affected by intrinsic selection at early generational hybrids, most of the sampled F9 RILs have recovered fitness to the level of their parentals. Together, these results suggest that a few generations of recombination allows for an independent segregation of the genes underlying thermal tolerance and cytonuclear incompatibilities, permitting certain recombinant lineages to survive in niches previously unused by parental taxa (i.e., warmer thermal environments) without incurring intrinsic selection.

Elevated oxidative damage is correlated with reduced fitness in interpopulation hybrids of a marine copepod.

Aerobic energy production occurs via the oxidative phosphorylation pathway (OXPHOS), which is critically dependent on interactions between the 13 mitochondrial DNA (mtDNA)-encoded and approximately 70 nuclear-encoded protein subunits. Disruptive mutations in any component of OXPHOS can result in impaired ATP production and exacerbated oxidative stress; in mammalian systems, such mutations are associated with ageing as well as numerous diseases. Recent studies have suggested that oxidative stress plays a role in fitness trade-offs in life-history evolution and functional ecology. Here, we show that outcrossing between populations with divergent mtDNA can exacerbate cellular oxidative stress in hybrid offspring. In the copepod Tigriopus californicus, we found that hybrids that showed evidence of fitness breakdown (low fecundity) also exhibited elevated levels of oxidative damage to DNA, whereas those with no clear breakdown did not show significantly elevated damage. The extent of oxidative stress in hybrids appears to be dependent on the degree of genetic divergence between their respective parental populations, but this pattern requires further testing using multiple crosses at different levels of divergence. Given previous evidence in T. californicus that hybridization disrupts nuclear/mitochondrial interactions and reduces hybrid fitness, our results suggest that such negative intergenomic epistasis may also increase the production of damaging cellular oxidants; consequently, mtDNA evolution may play a significant role in generating postzygotic isolating barriers among diverging populations.

Evidence for compensatory evolution of ribosomal proteins in response to rapid divergence of mitochondrial rRNA.

Rapid evolution of mitochondrial DNA (mtDNA) places intrinsic selective pressures on many nuclear genes involved in mitochondrial functions. Mitochondrial ribosomes, for example, are composed of mtDNA-encoded ribosomal RNAs (rRNAs) and a set of more than 60 nuclear-encoded ribosomal proteins (mRP) distinct from the cytosolic RPs (cRP). We hypothesized that the rapid divergence of mt-rRNA would result in rapid evolution of mRPs relative to cRPs, which respond to slowly evolving nuclear-encoded rRNA. In comparisons of rates of nonsynonymous and synonymous substitutions between a pair of divergent populations of the copepod Tigriopus californicus, we found that mRPs showed elevated levels of amino acid changes relative to cRPs. This pattern was equally strong at the interspecific level, between three pairs of sister species (Nasonia vitripennis vs. N. longicornis, Drosophila melanogaster vs. D. simulans, and Saccharomyces cerevisae vs. S. paradoxus). This high rate of mRP evolution may result in intergenomic incompatibilities between taxonomic lineages, and such incompatibilities could lead to dysfunction of mitochondrial ribosomes and the loss of fitness observed among interpopulation hybrids in T. californicus and interspecific hybrids in other species.

Investigating the molecular basis of local adaptation to thermal stress: population differences in gene expression across the transcriptome of the copepod Tigriopus californicus.

BACKGROUND: Geographic variation in the thermal environment impacts a broad range of biochemical and physiological processes and can be a major selective force leading to local population adaptation. In the intertidal copepod Tigriopus californicus, populations along the coast of California show differences in thermal tolerance that are consistent with adaptation, i.e., southern populations withstand thermal stresses that are lethal to northern populations. To understand the genetic basis of these physiological differences, we use an RNA-seq approach to compare genome-wide patterns of gene expression in two populations known to differ in thermal tolerance. RESULTS: Observed differences in gene expression between the southern (San Diego) and the northern (Santa Cruz) populations included both the number of affected loci as well as the identity of these loci. However, the most pronounced differences concerned the amplitude of up-regulation of genes producing heat shock proteins (Hsps) and genes involved in ubiquitination and proteolysis. Among the hsp genes, orthologous pairs show markedly different thermal responses as the amplitude of hsp response was greatly elevated in the San Diego population, most notably in members of the hsp70 gene family. There was no evidence of accelerated evolution at the sequence level for hsp genes. Among other sets of genes, cuticle genes were up-regulated in SD but down-regulated in SC, and mitochondrial genes were down-regulated in both populations. CONCLUSIONS: Marked changes in gene expression were observed in response to acute sub-lethal thermal stress in the copepod T. californicus. Although some qualitative differences were observed between populations, the most pronounced differences involved the magnitude of induction of numerous hsp and ubiquitin genes. These differences in gene expression suggest that evolutionary divergence in the regulatory pathway(s) involved in acute temperature stress may offer at least a partial explanation of population differences in thermal tolerance observed in Tigriopus.