Tigriopus Californicus Research Articles

Reproductive and phylogenetic divergence of tidepool copepod populations across a narrow geographical boundary in Baja California

AbstractAimPrevious work on the tidepool copepodTigriopus californicusrevealed a curious case of incipient speciation at the southern end of the species' range in Baja California, Mexico. The present study expands on the geography of this pattern and tests for congruence between reproductive and phylogenetic patterns.LocationThe Pacific coast of North America, from central Baja California to south‐eastern Alaska (27–57° N), including the full range ofT. californicus.MethodsPrimary techniques included mating experiments (> 4000 crosses), phylogeny reconstruction (mitochondrial cytochromecoxidase subunit I) and screening of single nucleotide polymorphisms (SNPs, 42 loci). Analyses used > 8000 copepods for the mating experiments, 86 copepods for the phylogeny and 41 copepods for theSNPassays. Phylogenies were constructed using Bayesian, maximum likelihood and maximum parsimony methods.ResultsPopulations were found to fall into three reproductive groups: northern and southern groups that were reproductively isolated from each other, and an intermediate group that could serve as a conduit for gene flow. The northern and intermediate populations fell into one clade while all southern populations fell into a second clade. These two clades are now separated by less than 12 km at latitude 29.35° N. NuclearSNPdata for a subset of locations confirmed striking divergence between populations on either side of this boundary. The second (southern) clade was further subdivided into two clades separated by the lagoon region of Guerrero Negro (latitude 28° N).Main conclusionsReproductive assays and molecular data (both mitochondrial and nuclear) reveal a sharp break at 29.35° N, a region with no obvious barriers to dispersal, with no evidence for mixing across this narrow transition zone. Results also showed a milder break at the Guerrero Negro Lagoon (28° N), a location where breaks have been reported for other taxa.

Gene conversion yields novel gene combinations in paralogs of GOT1 in the copepod Tigriopus californicus

BackgroundGene conversion of duplicated genes can slow the divergence of paralogous copies over time but can also result in other interesting evolutionary patterns. Islands of genetic divergence that persist in the face of gene conversion can point to gene regions undergoing selection for new functions. Novel combinations of genetic variation that differ greatly from the original sequence can result from the transfer of genetic variation between paralogous genes by rare gene conversion events. Genetically divergent populations of the copepod Tigriopus californicus provide an excellent model to look at the patterns of divergence among paralogs across multiple independent evolutionary lineages.ResultsIn this study the evolution of a set of paralogous genes encoding putative aspartate transaminase proteins (called GOT1 here) are examined in populations of the copepod T. californicus. One pair of duplicated genes, GOT1p1 and GOT1p2, has regions of high divergence between the copies in the face of apparent on-going gene conversion. The GOT1p2 gene also has unique haplotypes in two populations that appear to have resulted from a transfer of genetic variation via inter-paralog gene conversion. A second pair of duplicated genes GOT1Sr and GOT1Sd also shows evidence of gene conversion, but this gene conversion does not appear to have maintained each as a functional copy in all populations.ConclusionsThe patterns of conservation and sequence divergence across this set of paralogous genes among populations of T. californicus suggest that some interesting evolutionary patterns are occurring at these loci. The results for the GOT1p1/GOT1p2 paralogs illustrate how gene conversion can factor in the creation of a mosaic pattern of regions of high divergence and low divergence. When coupled with rare gene conversion events of divergent regions, this pattern can result in the formation of novel proteins differing substantially from either original protein. The evolutionary patterns across these paralogs show how gene conversion can both constrain and facilitate diversification of genetic sequences.

Fitness and morphological outcomes of many generations of hybridization in the copepod Tigriopus californicus

Hybridization between genetically divergent populations is an important evolutionary process, with an outcome that is difficult to predict. We used controlled crosses and freely mating hybrid swarms, followed for up to 30 generations, to examine the morphological and fitness consequences of interpopulation hybridization in the copepod Tigriopus californicus. Patterns of fitness in two generations of controlled crosses were partly predictive of long-term trajectories in hybrid swarms. For one pair of populations, controlled crosses revealed neutral or beneficial effects of hybridization after the F1 generation, and hybrid swarm fitness almost always equalled or exceeded that of the midparent. For a second pair, controlled crosses showed F2 hybrid breakdown, but increased fitness in backcrosses, and hybrid swarm fitness deviated both above and below that of the parentals. Nevertheless, individual swarm replicates exhibited different fitness trajectories over time that were not related in a simple manner to their hybrid genetic composition, and fixation of fitter hybrid phenotypes was not observed. Hybridization did not increase overall morphological variation, and underlying genetic changes may have been masked by phenotypic plasticity. Nevertheless, one type of hybrid swarm exhibited a repeatable pattern of transgressively large eggsacs, indicating a positive effect of hybridization on individual fecundity. Additionally, both parental and hybrid swarms exhibited common phenotypic trends over time, indicating common selective pressures in the laboratory environment. Our results suggest that, in a system where much work has focused on F2 hybrid breakdown, the long-term fitness consequences of interpopulation hybridization are surprisingly benign.

Habitat characteristics and metapopulation dynamics of the copepod Tigriopus californicus

Many organisms live in spatially delimited habitat patches, where local extinctions are compensated by dispersal and recolonization of empty habitat patches. To predict the dynam- ics of such metapopulations in an era of large environmental changes, it is essential to understand the key abiotic factors affecting local occurrence and temporal variation in patch occupancy. Here, we investigated the metapopulation dynamics of the marine copepod Tigriopus californicus, which is restricted to high intertidal and supralittoral rock pools. We monitored populations of T. californicus in northern California for almost 2 yr and observed pronounced seasonal changes in patch occupancy with >80% population turnover annually. Copepod occupancy depended on the area, depth and desiccation of pools and shore height of pools. Most local extinctions were due to pool desiccation, and the pools with an intermediate desiccation risk and low shore heights were the most frequently occupied. We highlight the significance of both the spatial position of pools (shore height) and their ephemerality in explaining occupancy pattern.

THE GENOMIC TRAJECTORY OF HYBRID SWARMS: OUTCOMES OF REPEATED CROSSES BETWEEN POPULATIONS OFTIGRIOPUS CALIFORNICUS

Introgressive hybridization between genetically divergent populations is an important evolutionary process. The degree to which repeated hybridization events between the same parental taxa lead to similar genomic outcomes is unknown. This study addressed this question by following genomic trajectories of replicate hybrid swarms of the copepod Tigriopus californicus over many generations of free mating. Swarm composition was determined both by differential reproductive success of founder individuals and subsequent selection on hybrid genotypes. For one cross, between two populations showing differential fitness in the laboratory and no hybrid breakdown, the genetic trajectory was highly repeatable: replicates rapidly became dominated by alleles from the fitter parent. In a second cross, between two populations showing similar fitness and significant F2 hybrid breakdown, alleles from alternative populations dominated different replicates. Swarms exhibited a general temporal trend of decreasing cytonuclear mismatch. Some patterns of differential introgression across the genome were strikingly congruent amongst swarm replicates, both within and between cross types, and reflected patterns of segregation distortion previously observed within controlled crosses between the same parental populations. Differences in heterozygosity between the sexes, and evidence for a previously suspected sex-distortion locus, suggest that complex interactions between sex and genotype influence hybrid swarm outcome.

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.

The minute brain of the copepod Tigriopus californicus supports a complex ancestral ground pattern of the tetraconate cerebral nervous systems

Copepods are a diverse and ecologically crucial group of minute crustaceans that are relatively neglected in terms of studies on nervous system organization. Recently, morphological neural characters have helped clarify evolutionary relationships within Arthropoda, particularly among Tetraconata (i.e., crustaceans and hexapods), and indicate that copepods occupy an important phylogenetic position relating to both Malacostraca and Hexapoda. This taxon therefore provides the opportunity to evaluate those neural characters common to these two clades likely to be results of shared ancestry (homology) versus convergence (homoplasy). Here we present an anatomical characterization of the brain and central nervous system of the well-studied harpacticoid copepod species Tigriopus californicus. We show that this species is endowed with a complex brain possessing a central complex comprising a protocerebral bridge and central body. Deutocerebral glomeruli are supplied by the antennular nerves, and a lateral protocerebral olfactory neuropil corresponds to the malacostracan hemiellipsoid body. Glomeruli contain synaptic specializations comparable to the presynaptic "T-bars" typical of dipterous insects, including Drosophila melanogaster. Serotonin-like immunoreactivity pervades the brain and ventral nervous system, with distinctive deutocerebral distributions. The present observations suggest that a suite of morphological characters typifying the Tigriopus brain reflect a ground pattern organization of an ancestral Tetraconata, which possessed an elaborate and structurally differentiated nervous system.

Diversifying Selection Underlies the Origin of Allozyme Polymorphism at the Phosphoglucose Isomerase Locus in Tigriopus californicus

The marine copepod Tigriopus californicus lives in intertidal rock pools along the Pacific coast, where it exhibits strong, temporally stable population genetic structure. Previous allozyme surveys have found high frequency private alleles among neighboring subpopulations, indicating that there is limited genetic exchange between populations. Here we evaluate the factors responsible for the diversification and maintenance of alleles at the phosphoglucose isomerase (Pgi) locus by evaluating patterns of nucleotide variation underlying previously identified allozyme polymorphism. Copepods were sampled from eleven sites throughout California and Baja California, revealing deep genetic structure among populations as well as genetic variability within populations. Evidence of recombination is limited to the sample from Pescadero and there is no support for linkage disequilibrium across the Pgi locus. Neutrality tests and codon-based models of substitution suggest the action of natural selection due to elevated non-synonymous substitutions at a small number of sites in Pgi. Two sites are identified as the charge-changing residues underlying allozyme polymorphisms in T. californicus. A reanalysis of allozyme variation at several focal populations, spanning a period of 26 years and over 200 generations, shows that Pgi alleles are maintained without notable frequency changes. Our data suggest that diversifying selection accounted for the origin of Pgi allozymes, while McDonald-Kreitman tests and the temporal stability of private allozyme alleles suggests that balancing selection may be involved in the maintenance of amino acid polymorphisms within populations.

Quantifying the Elevation of Mitochondrial DNA Evolutionary Substitution Rates Over Nuclear Rates in the Intertidal Copepod Tigriopus californicus

Mitochondrial DNA (mtDNA) genomes generally evolve rapidly in animals, but considerable variation in the rates of evolution of mtDNA occurs among taxa. Higher levels of mutation will tend to increase the amount of polymorphism, which should also scale with population size, but there are mixed signals from previous studies on the evolutionary outcomes of the interactions of these processes. The copepod Tigriopus californicus provides an interesting model in which to study the evolution of mtDNA because it has high levels of divergence among populations and there is the suggestion that this divergence could be involved in reproductive isolation. This species also appears to have an elevated mtDNA substitution rate, but previous studies did not provide an accurate measurement. This article examines the rate of mtDNA substitution versus nuclear substitution in T. californicus and finds that the mtDNA rate for synonymous sites averages 55-fold higher, a level that exceeds the rates found in most other invertebrates. Levels of polymorphism are also examined in both mtDNA and nuclear genes, and it is shown that the effective population size of mtDNA genes is much lower than that of nuclear genes. In addition, no correlation between polymorphism in mtDNA and nuclear genes is found across populations, which suggest factors other than demography may shape polymorphism in this species. The results from this study suggest that mtDNA is evolving at a very rapid rate in this copepod species, and this could increase the likelihood that mtDNA evolution is involved in the generation of reproductive isolation.

Long-term experimental hybrid swarms between nearly incompatible Tigriopus californicus populations: persistent fitness problems and assimilation by the superior population

For the intertidal copepod Tigriopus californicus, outbreeding depression for a variety of fitness measures is typically observed in early-generation interpopulation hybrids. We examined both controlled crosses and long-term, freely mating experimental hybrid swarms composed of individuals from Baja California (Mexico) populations Playa Altamira and Punta Morro. In controlled crosses, F1 and F2 hybrids showed large and significant declines in hatching numbers compared to parentals, while reciprocal backcrosses produced no offspring at all. For long-term studies, four treatment groups were initiated: 100%PA, 100%PM, 50%PA: 50%PM, and 80%PA: 20%PM. Replicates were surveyed at 3-month intervals for morphometric, census and fitness measures. The PA and 80PA:20PM treatments had initial fitness below the PM treatment, and went extinct within the first 12 months of the experiment. The 50:50 treatment had fitness below the PM parent at the 3- and 6-month time points, recovered to equivalent or superior fitness from months 9 to 18, and dropped again below PM at month 21. Limited genotyping of diagnostic microsatellites was consistent with PM alleles going nearly to fixation in hybrid replicates and male morphological data were concordant with a shift toward PM values. Results were strikingly different from a recent study of a different pair of populations showing extensive introgression and superior fitness in hybrid populations. This demonstrates how long-term consequences of population mixing depend on the relative fitness and level of compatibility between hybridizing populations.

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

BackgroundGeographic 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.ResultsObserved 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.ConclusionsMarked 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.

A gene-based SNP resource and linkage map for the copepod Tigriopus californicus.

BackgroundAs yet, few genomic resources have been developed in crustaceans. This lack is particularly evident in Copepoda, given the extraordinary numerical abundance, and taxonomic and ecological diversity of this group. Tigriopus californicus is ideally suited to serve as a genetic model copepod and has been the subject of extensive work in environmental stress and reproductive isolation. Accordingly, we set out to develop a broadly-useful panel of genetic markers and to construct a linkage map dense enough for quantitative trait locus detection in an interval mapping framework for T. californicus--a first for copepods.ResultsOne hundred and ninety Single Nucleotide Polymorphisms (SNPs) were used to genotype our mapping population of 250 F2 larvae. We were able to construct a linkage map with an average intermarker distance of 1.8 cM, and a maximum intermarker distance of 10.3 cM. All markers were assembled into linkage groups, and the 12 linkage groups corresponded to the 12 known chromosomes of T. californicus. We estimate a total genome size of 401.0 cM, and a total coverage of 73.7%. Seventy five percent of the mapped markers were detected in 9 additional populations of T. californicus. Of available model arthropod genomes, we were able to show more colocalized pairs of homologues between T. californicus and the honeybee Apis mellifera, than expected by chance, suggesting preserved macrosynteny between Hymenoptera and Copepoda.ConclusionsOur study provides an abundance of linked markers spanning all chromosomes. Many of these markers are also found in multiple populations of T. californicus, and in two other species in the genus. The genomic resource we have developed will enable mapping throughout the geographical range of this species and in closely related species. This linkage map will facilitate genome sequencing, mapping and assembly in an ecologically and taxonomically interesting group for which genomic resources are currently under development.

Hybrid Breakdown Weakens under Thermal Stress in Population Crosses of the Copepod Tigriopus californicus

The outcome of hybridization can be impacted by environmental conditions, which themselves can contribute to reproductive isolation between taxa. In crosses of genetically divergent populations, hybridization can have both negative and positive impacts on fitness, the balance between which might be tipped by changes in the environment. Genetically divergent populations of the intertidal copepod Tigriopus californicus have been shown to differ in thermal tolerance at high temperatures along a latitudinal gradient. In this study, a series of crosses were made between pairs of genetically divergent populations of T. californicus, and the thermal tolerance of these hybrids was tested. In most cases, the first-generation hybrids had relatively high thermal tolerance and the second-generation hybrids were not generally reduced below the less-tolerant parental population for high temperature tolerance. This pattern contrasts with previous studies from crosses of genetically divergent populations of this copepod, which often shows hybrid breakdown in these second-generation hybrids for other measures of fitness. These results suggest that high temperature stress could either increase the positive impacts of hybridization or decrease the negative impacts of hybridization resulting in lowered hybrid breakdown in these population crosses.

Complex Deleterious Interactions Associated with Malic Enzyme May Contribute to Reproductive Isolation in the Copepod Tigriopus californicus

Dobzhansky-Muller incompatibilities can result from the interactions of more than a single pair of interacting genes and there are several different models of how such complex interactions can be structured. Previous empirical work has identified complex conspecific epistasis as a form of complex interaction that has contributed to postzygotic reproductive isolation between taxa, but other forms of complexity are also possible. Here, I probe the genetic basis of reproductive isolation in crosses of the intertidal copepod Tigriopus californicus by looking at the impact of markers in genes encoding metabolic enzymes in F2 hybrids. The region of the genome associated with the locus ME2 is shown to have strong, repeatable impacts on the fitness of hybrids in crosses and epistatic interactions with another chromosomal region marked by the GOT2 locus in one set of crosses. In a cross between one of these populations and a third population, these two regions do not appear to interact despite the continuation of a large effect of the ME2 region itself in both crosses. The combined results suggest that the ME2 chromosomal region is involved in incompatibilities with several unique partners. If these deleterious interactions all stem from the same factor in this region, that would suggest a different form of complexity from complex conspecific epistasis, namely, multiple independent deleterious interactions stemming from the same factor. Confirmation of this idea will require more fine-scale mapping of the interactions of the ME2 region of the genome.

Limited potential for adaptation to climate change in a broadly distributed marine crustacean.

The extent to which acclimation and genetic adaptation might buffer natural populations against climate change is largely unknown. Most models predicting biological responses to environmental change assume that species' climatic envelopes are homogeneous both in space and time. Although recent discussions have questioned this assumption, few empirical studies have characterized intraspecific patterns of genetic variation in traits directly related to environmental tolerance limits. We test the extent of such variation in the broadly distributed tidepool copepod Tigriopus californicus using laboratory rearing and selection experiments to quantify thermal tolerance and scope for adaptation in eight populations spanning more than 17° of latitude. Tigriopus californicus exhibit striking local adaptation to temperature, with less than 1 per cent of the total quantitative variance for thermal tolerance partitioned within populations. Moreover, heat-tolerant phenotypes observed in low-latitude populations cannot be achieved in high-latitude populations, either through acclimation or 10 generations of strong selection. Finally, in four populations there was no increase in thermal tolerance between generations 5 and 10 of selection, suggesting that standing variation had already been depleted. Thus, plasticity and adaptation appear to have limited capacity to buffer these isolated populations against further increases in temperature. Our results suggest that models assuming a uniform climatic envelope may greatly underestimate extinction risk in species with strong local adaptation.

Long-term experimental hybrid swarms between moderately incompatible Tigriopus californicus populations: hybrid inferiority in early generations yields to hybrid superiority in later generations

The deleterious effects of hybridization are a serious concern for the conservation and management of species, particularly when populations mix as a result of human activity. Outbreeding depression is the typical result observed in early-generation interpopulation hybrids of Tigriopus californicus. We examined both controlled crosses and long-term, freely-mating, experimental hybrid populations composed of southern California populations Royal Palms (RP) and San Diego (SD). Controlled crosses included parentals plus all reciprocal F1, F2, F3 and backcross cohorts, and only F2 cohorts showed significant declines in fitness compared to midparent values, indicating recovery in the F3. For long-term studies, four treatment groups were initiated: 100% RP, 100% SD, 50% RP: 50% SD, and 80% RP: 20% SD. Replicates were surveyed at 3-month intervals for morphometric, census and fitness measures. Fitness of hybrid treatments showed declines relative to midparent values followed by rapid recovery, with two hybrid replicates ultimately showing higher fitness than parentals at the final 15-month time-point (up to 20 generations). In contrast, both males and females in hybrid treatments were larger than the midparent for several morphometric characters at the first time-point, and smaller than the midparent at the final time-point, indicating a possible tradeoff between fitness and body size. Microsatellites for a subset of samples revealed extensive introgression in hybrid treatments. This adds to previous evidence that hybrid breakdown in early generations may be a temporary phenomenon followed by the persistence of highly fit recombinant genotypes.

Interpopulation hybridization results in widespread viability selection across the genome in Tigriopus californicus

BackgroundGenetic interactions within hybrids influence their overall fitness. Understanding the details of these interactions can improve our understanding of speciation. One experimental approach is to investigate deviations from Mendelian expectations (segregation distortion) in the inheritance of mapped genetic markers. In this study, we used the copepod Tigriopus californicus, a species which exhibits high genetic divergence between populations and a general pattern of reduced fitness in F2 interpopulation hybrids. Previous studies have implicated both nuclear-cytoplasmic and nuclear-nuclear interactions in causing this fitness reduction. We identified and mapped population-diagnostic single nucleotide polymorphisms (SNPs) and used these to examine segregation distortion across the genome within F2 hybrids.ResultsWe generated a linkage map which included 45 newly elucidated SNPs and 8 population-diagnostic microsatellites used in previous studies. The map, the first available for the Copepoda, was estimated to cover 75% of the genome and included markers on all 12 T. californicus chromosomes. We observed little segregation distortion in newly hatched F2 hybrid larvae (fewer than 10% of markers at p < 0.05), but strikingly higher distortion in F2 hybrid adult males (45% of markers at p < 0.05). Hence, segregation distortion was primarily caused by selection against particular genetic combinations which acted between hatching and maturity. Distorted markers were not distributed randomly across the genome but clustered on particular chromosomes. In contrast to other studies in this species we found little evidence for cytonuclear coadaptation. Instead, different linkage groups exhibited markedly different patterns of distortion, which appear to have been influenced by nuclear-nuclear epistatic interactions and may also reflect genetic load carried within the parental lines.ConclusionAdult male F2 hybrids between two populations of T. californius exhibit dramatic segregation distortion across the genome. Distorted loci are clustered within specific linkage groups, and the direction of distortion differs between chromosomes. This segregation distortion is due to selection acting between hatching and adulthood.

Complex Deleterious Interactions Associated with Malic Enzyme May Contribute to Reproductive Isolation in the Copepod Tigriopus californicus

Complex Deleterious Interactions Associated with Malic Enzyme May Contribute to Reproductive Isolation in the Copepod Tigriopus californicus

Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus

The accumulation of genetic incompatibilities between isolated populations is thought to lead to the evolution of intrinsic postzygotic isolation. The molecular basis for these mechanisms, however, remains poorly understood. The intertidal copepod Tigriopus californicus provides unique opportunities for addressing mechanistic questions regarding the early stages of speciation; hybrids between highly divergent populations are fertile and viable, but exhibit reduced fitness at the F(2) or later generations. Given the current scarcity of genomic information in taxa at incipient stages of reproductive isolation, we utilize high-throughout 454 pyrosequencing to characterize a substantial fraction of protein-coding regions (the transcriptome) of T. californicus. Our sequencing effort was divided equally between two divergent populations in order to estimate levels of divergence and to reveal patterns of selection across the transcriptome. Assembly of sequences generated over 40,000 putatively unique transcripts (unigenes) for each population, 19,622 of which were orthologous between populations. BLAST searches of public databases determined protein identity and functional features for 15,402 and 12,670 unigenes, respectively. Based on rates of nonsynonymous and synonymous substitutions in 5897 interpopulation orthologs (those >150 bp and with at least 2X coverage), we identified 229 potential targets of positive selection. Many of these genes are predicted to be involved in several metabolic processes, and to function in hydrolase, peptidase and binding activities. The library of T. californicus coding regions, annotated with their predicted functions and level of divergence, will serve as an invaluable resource for elucidating molecular mechanisms underlying the early stages of speciation.

The nature of interactions that contribute to postzygotic reproductive isolation in hybrid copepods

Deleterious interactions within the genome of hybrids can lower fitness and result in postzygotic reproductive isolation. Understanding the genetic basis of these deleterious interactions, known as Dobzhansky-Muller incompatibilities, is the subject of intense current study that seeks to elucidate the nature of these deleterious interactions. Hybrids from crosses of individuals from genetically divergent populations of the intertidal copepod Tigriopus californicus provide a useful model in which to study Dobzhansky-Muller incompatibilities. Studies of the basis of postzygotic reproductive isolation in this species have revealed a number of patterns. First, there is evidence for a breakdown in genomic coadaptation between mtDNA-encoded and nuclear-encoded proteins that can result in a reduction in hybrid fitness in some crosses. It appears from studies of the individual genes involved in these interactions that although this coadaptation could lead to asymmetries between crosses, patterns of genotypic viabilities are not often consistent with simple models of genomic coadaptation. Second, there is a large impact of environmental factors on these deleterious interactions suggesting that they are not strictly intrinsic in nature. Temperature in particular appears to play an important role in determining the nature of these interactions. Finally, deleterious interactions in these hybrid copepods appear to be complex in terms of the number of genetic factors that interact to lead to reductions in hybrid fitness. This complexity may stem from three or more factors that all interact to cause a single incompatibility or the same factor interacting with multiple other factors independently leading to multiple incompatibilities.