Differential gene expression and mitonuclear incompatibilities in fast- and slow-developing interpopulation Tigriopus californicus hybrids.

Mitochondrial functions are intimately reliant on proteins and RNAs encoded in both the nuclear and mitochondrial genomes, leading to inter-genomic coevolution within taxa. Hybridization can break apart coevolved mitonuclear genotypes, resulting in decreased mitochondrial performance and reduced fitness. This hybrid breakdown is an important component of outbreeding depression and early-stage reproductive isolation. However, the mechanisms contributing to mitonuclear interactions remain poorly resolved. Here, we scored variation in developmental rate (a proxy for fitness) among reciprocal F2 interpopulation hybrids of the intertidal copepod Tigriopus californicus and used RNA sequencing to assess differences in gene expression between fast- and slow-developing hybrids. In total, differences in expression associated with developmental rate were detected for 2925 genes, whereas only 135 genes were differentially expressed as a result of differences in mitochondrial genotype. Upregulated expression in fast developers was enriched for genes involved in chitin-based cuticle development, oxidation-reduction processes, hydrogen peroxide catabolic processes and mitochondrial respiratory chain complex I. In contrast, upregulation in slow developers was enriched for DNA replication, cell division, DNA damage and DNA repair. Eighty-four nuclear-encoded mitochondrial genes were differentially expressed between fast- and slow-developing copepods, including 12 subunits of the electron transport system (ETS) which all had higher expression in fast developers than in slow developers. Nine of these genes were subunits of ETS complex I. Our results emphasize the major roles that mitonuclear interactions within the ETS, particularly in complex I, play in hybrid breakdown, and resolve strong candidate genes for involvement in mitonuclear interactions.

The role of mitonuclear incompatibilities in allopatric speciation.

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Strong selective effects of mitochondrial DNA on the nuclear genome.

Oxidative phosphorylation, the primary source of cellular energy in eukaryotes, requires gene products encoded in both the nuclear and mitochondrial genomes. As a result, functional integration between the genomes is essential for efficient adenosine triphosphate (ATP) generation. Although within populations this integration is presumably maintained by coevolution, the importance of mitonuclear coevolution in key biological processes such as speciation and mitochondrial disease has been questioned. In this study, we crossed populations of the intertidal copepod Tigriopus californicus to disrupt putatively coevolved mitonuclear genotypes in reciprocal F2 hybrids. We utilized interindividual variation in developmental rate among these hybrids as a proxy for fitness to assess the strength of selection imposed on the nuclear genome by alternate mitochondrial genotypes. Developmental rate varied among hybrid individuals, and in vitro ATP synthesis rates of mitochondria isolated from high-fitness hybrids were approximately two-fold greater than those of mitochondria isolated from low-fitness individuals. We then used Pool-seq to compare nuclear allele frequencies for high- or low-fitness hybrids. Significant biases for maternal alleles were detected on 5 (of 12) chromosomes in high-fitness individuals of both reciprocal crosses, whereas maternal biases were largely absent in low-fitness individuals. Therefore, the most fit hybrids were those with nuclear alleles that matched their mitochondrial genotype on these chromosomes, suggesting that mitonuclear effects underlie individual-level variation in developmental rate and that intergenomic compatibility is critical for high fitness. We conclude that mitonuclear interactions can have profound impacts on both physiological performance and the evolutionary trajectory of the nuclear genome.

Recovery from hybrid breakdown reveals a complex genetic architecture of mitonuclear incompatibilities.

Reproductive isolation is often achieved when genes that are neutral or beneficial in their genomic background become functionally incompatible in a foreign genomic background, causing inviability, sterility or other forms of low fitness in hybrids. Recent studies suggest that mitonuclear interactions are among the initial incompatibilities to evolve at early stages of population divergence across taxa. Yet, the genomic architecture of mitonuclear incompatibilities has rarely been elucidated. We employ an experimental evolution approach starting with low-fitness F2 interpopulation hybrids of the copepod Tigriopus californicus, in which frequencies of compatible and incompatible nuclear alleles change in response to an alternative mitochondrial background. After about nine generations, we observe a generalized increase in population size and in survivorship, suggesting efficiency of selection against maladaptive phenotypes. Whole genome sequencing of evolved populations showed some consistent allele frequency changes across three replicates of each reciprocal cross, but markedly different patterns between mitochondrial backgrounds. In only a few regions (~6.5% of the genome), the same parental allele was overrepresented irrespective of the mitochondrial background. About 33% of the genome showed allele frequency changes consistent with divergent selection, with the location of these genomic regions strongly differing between mitochondrial backgrounds. In 87% and 89% of these genomic regions, the dominant nuclear allele matched the associated mitochondrial background, consistent with mitonuclear co-adaptation. These results suggest that mitonuclear incompatibilities have a complex polygenic architecture that differs between populations, potentially generating genome-wide barriers to gene flow between closely related taxa.

Consequences of HSF knockdown on gene expression during the heat shock response in Tigriopus californicus.

Although the existence of a cellular heat shock response is nearly universal, its relationship to organismal thermal tolerance is not completely understood. Many of the genes involved are known to be regulated by the highly conserved heat shock transcription factor-1 (HSF-1), yet the regulatory network is not fully characterized. Here, we investigated the role of HSF-1 in gene expression following thermal stress using knockdown of HSF-1 by RNA interference in the intertidal copepod Tigriopus californicus We observed some evidence for decreased transcription of heat shock protein genes following knockdown, supporting the widely acknowledged role of HSF-1 in the heat shock response. However, the majority of differentially expressed genes between the control and HSF-1 knockdown groups were upregulated, suggesting that HSF-1 normally functions to repress their expression. Differential expression observed in genes related to chitin and cuticle formation lends support to previous findings that these processes are highly regulated following heat stress. We performed a genome scan and identified a set of 396 genes associated with canonical heat shock elements. RNA-seq data did not find those genes to be more highly represented in our HSF-1 knockdown treatment, indicating that requirements for binding and interaction of HSF-1 with a given gene are not simply predicted by the presence of HSF-1 binding sites. Further study of the pathways implicated by these results and future comparisons among populations of T. californicus may help us understand the role and importance of HSF-1 in the heat shock response and, more broadly, in organismal thermal tolerance.

Allele-Specific Expression and Evolution of Gene Regulation Underlying Acute Heat Stress Response and Local Adaptation in the Copepod Tigriopus californicus.

Geographic variation in environmental temperature can select for local adaptation among conspecific populations. Divergence in gene expression across the transcriptome is a key mechanism for evolution of local thermal adaptation in many systems, yet the genetic mechanisms underlying this regulatory evolution remain poorly understood. Here we examine gene expression in 2 locally adapted Tigriopus californicus populations (heat tolerant San Diego, SD, and less tolerant Santa Cruz, SC) and their F1 hybrids during acute heat stress response. Allele-specific expression (ASE) in F1 hybrids was used to determine cis-regulatory divergence. We found that the number of genes showing significant allelic imbalance increased under heat stress compared to unstressed controls. This suggests that there is significant population divergence in cis-regulatory elements underlying heat stress response. Specifically, the number of genes showing an excess of transcripts from the more thermal tolerant (SD) population increased with heat stress while that number of genes with an SC excess was similar in both treatments. Inheritance patterns of gene expression also revealed that genes displaying SD-dominant expression phenotypes increase in number in response to heat stress; that is, across loci, gene expression in F1's following heat stress showed more similarity to SD than SC, a pattern that was absent in the control treatment. The observed patterns of ASE and inheritance of gene expression provide insight into the complex processes underlying local adaptation and thermal stress response.

Multiple Modes of Adaptation: Regulatory and Structural Evolution in a Small Heat Shock Protein Gene.

Thermal tolerance is a key determinant of species distribution. Despite much study, the genetic basis of adaptive evolution of thermal tolerance, including the relative contributions of transcriptional regulation versus protein evolution, remains unclear. Populations of the intertidal copepod Tigriopus californicus are adapted to local thermal regimes across their broad geographic range. Upon thermal stress, adults from a heat tolerant southern population, San Diego (SD), upregulate several heat shock proteins (HSPs) to higher levels than those from a less tolerant northern population, Santa Cruz (SC). Suppression of a specific HSP, HSPB1, significantly reduces T. californicus survival following acute heat stress. Sequencing of HSPB1 revealed population specific nucleotide substitutions in both promoter and coding regions of the gene. HSPB1 promoters from heat tolerant populations contain two canonical heat shock elements (HSEs), the binding sites for heat shock transcription factor (HSF), whereas less tolerant populations have mutations in these conserved motifs. Allele specific expression of HSPB1 in F1 hybrids between tolerant and less tolerant populations showed significantly biased expression favoring alleles from tolerant populations and supporting the adaptive divergence in these cis-regulatory variants. The functional impact of population-specific nonsynonymous substitutions in HSPB1 coding sequences was tested by assessing the thermal stabilization properties of SD versus SC HSPB1 protein variants. Recombinant HSPB1 from the southern SD population showed greater capacity for protecting protein structure under elevated temperature. Our results indicate that both regulatory and protein coding sequence evolution within a single gene appear to contribute to thermal tolerance phenotypes and local adaptation among conspecific populations.

Variation in developmental temperature alters adulthood plasticity of thermal tolerance in Tigriopus californicus.

In response to environmental change, organisms rely on both genetic adaptation and phenotypic plasticity to adjust key traits that are necessary for survival and reproduction. Given the accelerating rate of climate change, plasticity may be particularly important. For organisms in warming aquatic habitats, upper thermal tolerance is likely to be a key trait, and many organisms express plasticity in this trait in response to developmental or adulthood temperatures. Although plasticity at one life stage may influence plasticity at another life stage, relatively little is known about this possibility for thermal tolerance. Here, we used locally adapted populations of the copepod Tigriopus californicus to investigate these potential effects in an intertidal ectotherm. We found that low latitude populations had greater critical thermal maxima (CTmax) than high latitude populations, and variation in developmental temperature altered CTmax plasticity in adults. After development at 25°C, CTmax was plastic in adults, whereas no adulthood plasticity in this trait was observed after development at 20°C. This pattern was identical across four populations, suggesting that local thermal adaptation has not shaped this effect among these populations. Differences in the capacities to maintain ATP synthesis rates and to induce heat shock proteins at high temperatures, two likely mechanisms of local adaptation in this species, were consistent with changes in CTmax owing to phenotypic plasticity, which suggests that there is likely mechanistic overlap between the effects of plasticity and adaptation. Together, these results indicate that developmental effects may have substantial impacts on upper thermal tolerance plasticity in adult ectotherms.

Adaptation to a latitudinal thermal gradient within a widespread copepod species: the contributions of genetic divergence and phenotypic plasticity.

Understanding how populations adapt to heterogeneous thermal regimes is essential for comprehending how latitudinal gradients in species diversification are formed, and how taxa will respond to ongoing climate change. Adaptation can occur by innate genetic factors, by phenotypic plasticity, or by a combination of both mechanisms. Yet, the relative contribution of such mechanisms to large-scale latitudinal gradients of thermal tolerance across conspecific populations remains unclear. We examine thermal performance in 11 populations of the intertidal copepod Tigriopus californicus, ranging from Baja California Sur (Mexico) to British Columbia (Canada). Common garden experiments show that survivorship to acute heat-stress differs between populations (by up to 3.8°C in LD50 values), reflecting a strong genetic thermal adaptation. Using a split-brood experiment with two rearing temperatures, we also show that developmental phenotypic plasticity is beneficial to thermal tolerance (by up to 1.3°C), and that this effect differs across populations. Although genetic divergence in heat tolerance strongly correlates with latitude and temperature, differences in the plastic response do not. In the context of climate warming, our results confirm the general prediction that low-latitude populations are most susceptible to local extinction because genetic adaptation has placed physiological limits closer to current environmental maxima, but our results also contradict the prediction that phenotypic plasticity is constrained at lower latitudes.

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.

Genomic evidence for ecological divergence against a background of population homogeneity in the marine snail Chlorostoma funebralis.

The balance between natural selection, gene flow and genetic drift is difficult to resolve in marine invertebrates with extensive dispersal and fluctuating population sizes. The intertidal snail Chlorostoma funebralis has planktonic larvae and previous work using mtDNA polymorphism reported no genetic population structure. Nevertheless, recent studies have documented differences in thermal tolerance and transcriptomic responses to heat stress between northern and southern California, USA, populations. To gain insight into the dynamics influencing adaptive divergence, we used double-digest restriction site-associated DNA (ddRAD) sequencing to identify 1861 genomewide, quality-filtered single-nucleotide polymorphism (SNP) loci for C. funebralis collected from three northern and three southern California sites (15 individuals per population). Considering all SNPs, there was no evidence for genetic differentiation among populations or regions (average FST  = 0.0042). However, outlier tests revealed 34 loci putatively under divergent selection between northern and southern populations, and structure and SNP tree analyses based on these outliers show clear genetic differentiation between geographic regions. Three of these outliers are known or hypothesized to be involved in stress granule formation, a response to environmental stress such as heat. Combined with previous work that found thermally tolerant southern populations show high baseline expression of stress response genes, these results further suggest that thermal stress is a strong selective pressure across C. funebralis populations. Overall, this study increases our understanding of the factors constraining local adaptation in marine organisms, while suggesting that ecologically driven, strong differentiation can occur at relevant loci in a species with planktonic larvae.

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.

RNA-seq reveals regional differences in transcriptome response to heat stress in the marine snail Chlorostoma funebralis.

To investigate the role of gene expression in adaptation of marine ectotherms to different temperatures, we examined the transcriptome-wide thermal stress response in geographically separated populations of the intertidal snail Chlorostoma funebralis. Snails from two southern (heat tolerant) and two northern (heat sensitive) populations were acclimated to a common thermal environment, exposed to an environmentally relevant thermal stress and analysed using RNA-seq. Pooling across all populations revealed 306 genes with differential expression between control and heat-stressed samples, including 163 significantly upregulated and 143 significantly downregulated genes. When considered separately, regional differences in response were widely apparent. Heat shock proteins (Hsps) were upregulated in both regions, but the magnitude of response was significantly greater in northern populations for most Hsp70s, while the southern populations showed greater upregulation for approximately half of the Hsp40s. Of 177 stress-responsive genes in northern populations, 55 responded to heat stress only in northern populations. Several molecular chaperones and antioxidant genes that were not differentially expressed in southern populations showed higher expression under control conditions compared with northern populations. This suggests that evolution of elevated expression of these genes under benign conditions preadapts the southern populations to frequent heat stress and contributes to their higher thermal tolerance. These results indicate that evolution has resulted in different transcriptome responses across populations, including upregulation of genes in response to stress and preadaptation of genes in anticipation of stress (based on evolutionary history of frequent heat exposure). The relative importance of the two mechanisms differs among gene families and among 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.

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.

Evaluating the performance of captive breeding techniques for conservation hatcheries: a case study of the delta smelt captive breeding program.

The delta smelt, an endangered fish species endemic to the San Francisco Bay-Delta, California, United States, was recently brought into captivity for species preservation. This study retrospectively evaluates the implementation of a genetic management plan for the captive delta smelt population. The captive genetic management plan entails tagging fish, molecular data collection, pedigree reconstruction, relatedness estimation, and recommending fish crosses annually in an effort to minimize the average coancestry in the population and limit inbreeding. We employed 12 microsatellite DNA markers to examine temporal genetic diversity in consecutive, discrete generations to determine the effects of intensive genetic management on the population and to quantify the amount of wild genetic diversity present within each captive generation. Wild fish are incorporated into the captive population each generation to minimize genetic drift, and 91% of the original founders are still represented in the F(3) generation. The average mean kinship in the third generation in captivity was 0.0035. There was no evidence of significant genetic divergence of the captive population from the wild population. The results of this study yield management insights into the practical application of genetic management plans for captive populations and conservation hatcheries, in an attempt to preserve the genetic integrity of endangered species.

Loss of mitochondrial performance at high temperatures is correlated with upper thermal tolerance among populations of an intertidal copepod.

Environmental temperatures have pervasive effects on the performance and tolerance of ectothermic organisms, and thermal tolerance limits likely play key roles underlying biogeographic ranges and responses to environmental change. Mitochondria are central to metabolic processes in eukaryotic cells, and these metabolic functions are thermally sensitive; however, potential relationships between mitochondrial function, thermal tolerance limits and local thermal adaptation in general remain unresolved. Loss of ATP synthesis capacity at high temperatures has recently been suggested as a mechanistic link between mitochondrial function and upper thermal tolerance limits. Here we use a common-garden experiment with seven locally adapted populations of intertidal copepods (Tigriopus californicus), spanning approximately 21.5° latitude, to assess genetically based variation in the thermal performance curves of maximal ATP synthesis rates in isolated mitochondria. These thermal performance curves displayed substantial variation among populations with higher ATP synthesis rates at lower temperatures (20-25 °C) in northern populations than in southern populations. In contrast, mitochondria from southern populations maintained ATP synthesis rates at higher temperatures than the temperatures that caused loss of ATP synthesis capacity in mitochondria from northern populations. Additionally, there was a tight correlation between the thermal limits of ATP synthesis and previously determined variation in upper thermal tolerance limits among populations. This suggests that mitochondria may play an important role in latitudinal thermal adaptation in T. californicus, and supports the hypothesis that loss of mitochondrial performance at high temperatures is linked to whole-organism thermal tolerance limits in this ectotherm.

Prezygotic reproductive barriers in precopulatory behavior of tidepool copepod species.

Complexity in prezygotic mating behavior can contribute to the emergence of sexual incompatibility and reproductive isolation. In this study, we performed behavioral tests with two tidepool copepod species of the genus Tigriopus to explore the possibility of precopulatory behavioral isolation. We found that interspecific mating attempts failed prior to genital contact, and that this failure occurred at different behavioral steps between reciprocal pairings. Our results suggest that prezygotic barriers may exist at multiple points of the behavioral process on both male and female sides, possibly due to interspecific differences in mate-recognition cues used at those "checkpoints." While many copepod species are known to show unique precopulatory mate-guarding behavior, the potential contribution of prezygotic behavioral factors to their isolation is not widely recognized. The pattern of sequential mate-guarding behaviors may have allowed the diversification of precopulatory communication and contributed to the evolutionary diversity of the Tigriopus copepods.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence.

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Developmental rate displays effects of inheritance but not of sex in interpopulation hybrids of Tigriopus californicus.

Coevolved genetic interactions within populations can be disrupted by hybridization resulting in loss of fitness in hybrid individuals (i.e., hybrid breakdown). However, the extent to which variation in fitness-related traits among hybrids is inherited across generations remains unclear, and variation in these traits may be sex-specific in hybrids due to differential effects of genetic incompatibilities in females and males. Here we present two experiments investigating variation in developmental rate among reciprocal interpopulation hybrids of the intertidal copepod Tigriopus californicus. Developmental rate is a fitness-related trait in this species that is affected by interactions between mitochondrial-encoded and nuclear-encoded genes in hybrids that result in variation in mitochondrial ATP synthesis capacities. First, we show that F2 -hybrid developmental rate is equivalent in two reciprocal crosses and is unaffected by sex, suggesting that breakdown of developmental rate is likely experienced equally by females and males. Second, we demonstrate that variation in developmental rate among F3 hybrids is heritable; times to copepodid metamorphosis of F4 offspring of fast-developing F3 parents (12.25 ± 0.05 days, µ ± SEM) were significantly faster than those of F4 offspring of slow-developing parents (14.58 ± 0.05 days). Third, we find that ATP synthesis rates in these F4 hybrids are unaffected by the developmental rates of their parents, but that mitochondria from females synthesize ATP at faster rates than mitochondria from males. Taken together, these results suggest that sex-specific effects vary among fitness-related traits in these hybrids, and that effects likely associated with hybrid breakdown display substantial inheritance across hybrid generations.