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.

Using Pool-seq to Search for Genomic Regions Affected by Hybrid Inviability in the copepod T. californicus.

The formation of reproductive barriers between allopatric populations involves the accumulation of incompatibilities that lead to intrinsic postzygotic isolation. The evolution of these incompatibilities is usually explained by the Dobzhansky-Muller model, where epistatic interactions that arise within the diverging populations, lead to deleterious interactions when they come together in a hybrid genome. These incompatibilities can lead to hybrid inviability, killing individuals with certain genotypic combinations, and causing the population's allele frequency to deviate from Mendelian expectations. Traditionally, hybrid inviability loci have been detected by genotyping individuals at different loci across the genome. However, this method becomes time consuming and expensive as the number of markers or individuals increases. Here, we test if a Pool-seq method can be used to scan the genome of F2 hybrids to detect genomic regions that are affected by hybrid inviability. We survey the genome of hybrids between 2 populations of the copepod Tigriopus californicus, and show that this method has enough power to detect even small changes in allele frequency caused by hybrid inviability. We show that allele frequency estimates in Pool-seq can be affected by the sampling of alleles from the pool of DNA during the library preparation and sequencing steps and that special considerations must be taken when aligning hybrid reads to a reference when the populations/species are divergent.

Adaptive evolution of M3 lysin--a candidate gamete recognition protein in the Mytilus edulis species complex.

Marine invertebrate gamete recognition proteins (GRPs) are classic examples of rapid adaptive evolution of reproductive proteins, and hybridizing Mytilus blue mussels allow us to study the evolution of GRPs during speciation following secondary contact. Even with frequent hybridization, positive selection drives divergence of M7 lysin, one of the three Mytilus egg vitelline envelope (VE) lysins. Mytilus trossulus and M. edulis form a broad hybrid zone in the Canadian Maritimes and eastern Maine, isolated by strong (but partial) gamete incompatibility. M7 lysin, however, is an unlikely GRP controlling this gametic incompativility, as earlier studies showed either weak or no positive selection and extensive introgression between the two species. We used reverse transcriptase-polymerase chain reaction and cloned several alleles of M3 lysin, a potent VE lysin encoded by a nonhomologous gene whose evolution has not been studied. McDonald-Kreitman and HKA tests reveal strong positive selection, which PAML branch-site models detect in 19.7% of the codons. Protein structure predictions show that replacements map exclusively to one face of the carbohydrate recognition domain (CRD) of this C-type lectin, with codons under positive selection localizing to CRD regions known to control ligand specificity. Polymorphism/divergence analyses show that selective sweep has purged M. edulis but not M. trossulus of polymorphism, and unique to M3 is an absence of fixed substitutions and broad haplotype sharing between M. edulis and Mediterranean M. galloprovincialis. Taken together, these results suggest that different lysins serve as GRPs in different Mytilus hybrid zones, with M3 likely co-opted to play this role in the western Atlantic.

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.

Admixture in Africanized honey bees (Apis mellifera) from Panamá to San Diego, California (U.S.A.).

The Africanized honey bee (AHB) is a New World amalgamation of several subspecies of the western honey bee (Apis mellifera), a diverse taxon historically grouped into four major biogeographic lineages: A (African), M (Western European), C (Eastern European), and O (Middle Eastern). In 1956, accidental release of experimentally bred "Africanized" hybrids from a research apiary in Sao Paulo, Brazil initiated a hybrid species expansion that now extends from northern Argentina to northern California (U.S.A.). Here, we assess nuclear admixture and mitochondrial ancestry in 60 bees from four countries (Panamá; Costa Rica, Mexico; U.S.A) across this expansive range to assess ancestry of AHB several decades following initial introduction and test the prediction that African ancestry decreases with increasing latitude. We find that AHB nuclear genomes from Central America and Mexico have predominately African genomes (76%-89%) with smaller contributions from Western and Eastern European lineages. Similarly, nearly all honey bees from Central America and Mexico possess mitochondrial ancestry from the African lineage with few individuals having European mitochondria. In contrast, AHB from San Diego (CA) shows markedly lower African ancestry (38%) with substantial genomic contributions from all four major honey bee lineages and mitochondrial ancestry from all four clades as well. Genetic diversity measures from all New World populations equal or exceed those of ancestral populations. Interestingly, the feral honey bee population of San Diego emerges as a reservoir of diverse admixture and high genetic diversity, making it a potentially rich source of genetic material for honey bee breeding.

Phenotypic and transcriptional response of Daphnia pulicaria to the combined effects of temperature and predation.

Daphnia, an ecologically important zooplankton species in lakes, shows both genetic adaptation and phenotypic plasticity in response to temperature and fish predation, but little is known about the molecular basis of these responses and their potential interactions. We performed a factorial experiment exposing laboratory-propagated Daphnia pulicaria clones from two lakes in the Sierra Nevada mountains of California to normal or high temperature (15°C or 25°C) in the presence or absence of fish kairomones, then measured changes in life history and gene expression. Exposure to kairomones increased upper thermal tolerance limits for physiological activity in both clones. Cloned individuals matured at a younger age in response to higher temperature and kairomones, while size at maturity, fecundity and population intrinsic growth were only affected by temperature. At the molecular level, both clones expressed more genes differently in response to temperature than predation, but specific genes involved in metabolic, cellular, and genetic processes responded differently between the two clones. Although gene expression differed more between clones from different lakes than experimental treatments, similar phenotypic responses to predation risk and warming arose from these clone-specific patterns. Our results suggest that phenotypic plasticity responses to temperature and kairomones interact synergistically, with exposure to fish predators increasing the tolerance of Daphnia pulicaria to stressful temperatures, and that similar phenotypic responses to temperature and predator cues can be produced by divergent patterns of gene regulation.

Genomic scans reveal multiple mito-nuclear incompatibilities in population crosses of the copepod Tigriopus californicus.

The evolution of intrinsic postzygotic isolation can be explained by the accumulation of Dobzhansky-Muller incompatibilities (DMI). Asymmetries in the levels of hybrid inviability and hybrid sterility are commonly observed between reciprocal crosses, a pattern that can result from the involvement of uniparentally inherited factors. The mitochondrial genome is one such factor that appears to participate in DMI in some crosses but the frequency of its involvement versus biparentally inherited factors is unclear. Here we assess the relative importance of incompatibilities between nuclear factors (nuclear-nuclear) versus those between mitochondrial and nuclear factors (mito-nuclear) in a species that lacks sex chromosomes. We used a Pool-seq approach to survey three crosses among genetically divergent populations of the copepod, Tigriopus californicus, for regions of the genome that are affected by hybrid inviability. Results from reciprocal crosses suggest that mito-nuclear incompatibilities are more common than nuclear-nuclear incompatibilities overall. These results suggest that in the presence of very high levels of nucleotide divergence between mtDNA haplotypes, mito-nuclear incompatibilities can be important for the evolution of intrinsic postzygotic isolation. This is particularly interesting considering this species lacks sex chromosomes, which have been shown to harbor a particularly high number of nuclear-nuclear DMI in several other species.

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.

Locally adapted populations of a copepod can evolve different gene expression patterns under the same environmental pressures.

As populations diverge in allopatry, but under similar thermal conditions, do similar thermal performance phenotypes evolve by maintaining similar gene expression patterns, or does genetic divergence lead to divergent patterns of gene expression between these populations? We used genetically divergent populations of the copepod Tigriopus californicus, whose performance at different thermal conditions is well characterized, to investigate transcriptome-wide expression responses under two different thermal regimes: (1) a nonvariable temperature regime and (2) a regime with variable temperature. Our results show the expression profiles of the response to these regimes differed substantially among populations, even for populations that are geographically close. This pattern was accentuated when populations were raised in the variable temperature environment. Less heat-tolerant populations mounted strong but divergent responses to the different thermal regimes, with a large heat-shock response observed in one population, and an apparent reduction in the expression of genes involved in basic cellular processes in the other. Our results suggest that as populations diverge in allopatry, they may evolve starkly different responses to changes in temperature, at the gene expression level, while maintaining similar thermal performance phenotypes.

Chromosome-Wide Impacts on the Expression of Incompatibilities in Hybrids of Tigriopus californicus.

Chromosome rearrangements such as inversions have been recognized previously as contributing to reproductive isolation by maintaining alleles together that jointly contribute to deleterious genetic interactions and postzygotic reproductive isolation. In this study, an impact of potential incompatibilities merely residing on the same chromosome was found in crosses of populations of the copepod Tigriopus californicus When genetically divergent populations of this copepod are crossed, hybrids show reduced fitness, and deviations from expected genotypic ratios can be used to determine regions of the genome involved in deleterious interactions. In this study, a set of markers was genotyped for a cross of two populations of T. californicus, and these markers show widespread deviations from Mendelian expectations, with entire chromosomes showing marked skew. Despite the importance of mtDNA/nuclear interactions in incompatibilities in this system in previous studies, in these crosses the expected patterns stemming from these interactions are not widely apparent. Females lack recombination in this species, and a striking difference is observed between male and female backcrosses. This suggests that the maintenance of multiple loci on individual chromosomes can enable some incompatibilities, perhaps playing a similar role in the initial rounds of hybridization to chromosomal rearrangements in preserving sets of alleles together that contribute to incompatibilities. Finally, it was observed that candidate pairs of incompatibility regions are not consistently interacting across replicates or subsets of these crosses, despite the repeatability of the deviations at many of the single loci themselves, suggesting that more complicated models of Dobzhansky-Muller incompatibilities may need to be considered.

Higher levels of sex chromosome heteromorphism are associated with markedly stronger reproductive isolation.

The two 'rules of speciation', Haldane's rule and the large X-effect, describe the genetic basis of postzygotic isolation, and have led to the realization that sex chromosomes play an important role in this process. However, a range of sex determination mechanisms exists in nature, not always involving sex chromosomes. Based on these 'rules of speciation', I test the hypothesis that the presence of sex chromosomes will contribute to a faster evolution of intrinsic postzygotic isolation. I show that taxa that do not have sex chromosomes evolve lower levels of postzygotic isolation than taxa with sex chromosomes, at a similar amount of genetic divergence. Taxa with young homomorphic sex chromosomes show an intermediate pattern compared to taxa with heteromorphic sex chromosomes and taxa without sex chromosomes. These results are consistent with predictions from the two 'rules of speciation', and emphasize the importance of sex chromosomes for the evolution of intrinsic postzygotic isolation.

Evolutionary consequences of introgression at M7 lysin, a gamete recognition locus, following secondary contact between blue mussel species.

Hybridizing populations of blue mussels, Mytilus edulis and Mytilus trossulus, in Cobscook Bay (eastern Maine) have been used by our laboratory to study the evolution of gamete incompatibility and molecular evolution of the vitelline coat lysin proteins expressed in sperm. The M7 lysin locus has been the most studied of the three lysins, but evidence for positive selection necessary to help confirm its role in gamete recognition in western Atlantic hybrid zones is contradictory. We developed an alternative test, based on rates of introgression at M7 lysin. Contrary to expectations, introgression at this locus is much higher (instead of much lower) than is introgression at neutral markers. In this article, we present simulations, constructed using synthetic populations of combinations of admixed genotypes, representing various hybrid offspring categories. Simulations produced variation in introgression across loci, but did not generate the massive introgression at M7 lysin observed in natural populations in Cobscook Bay. We consider these results in the context of selection operating on gamete recognition loci, both within and between species, during the third stage of allopatric speciation in Mytilus.