Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

A species delimitation approach in the Trochulus sericeus/hispidus complex reveals two cryptic species within a sharp contact zone

Background  

Mitochondrial DNA sequencing increasingly results in the recognition of genetically divergent, but morphologically cryptic lineages. Species delimitation approaches that rely on multiple lines of evidence in areas of co-occurrence are particularly powerful to infer their specific status. We investigated the species boundaries of two cryptic lineages of the land snail genus Trochulus in a contact zone, using mitochondrial and nuclear DNA marker as well as shell morphometrics.     

Nuclear genes (but not mitochondrial DNA barcodes) reveal real species: Evidence from the Brenthis fritillary butterflies (Lepidoptera,Nymphalidae)

Discordance between entities revealed by nuclear versus mitochondrial genes is a common phenomenon in evolutionary and taxonomic studies. However, little attention has been paid to analysis of how such discordant entities correspond to traditional species detected through investigation of their morphology, ecology, and distribution. Here, we used one mitochondrial (COI, DNA barcode fragment) and four nuclear (CAD, Ca‐ATPase, arginine kinase, wg) genes to analyze the genetic structure of the taxonomically well‐studied butterfly genus Brenthis (Lepidoptera, Nymphalidae). Analysis of COI revealed multiple diverged allopatric and sympatric mitochondrial lineages within the known Brenthis species hinting at possible presence of unrecognized cryptic species. However, these multiple‐species hypotheses were not supported by further studies of nuclear genes and phenotypic traits. The discovered mitochondrial lineages did not correspond to the clusters revealed by nuclear genes. Simultaneously, we found a complete congruence between (a) traditional species boundaries, (b) clusters recognized by nuclear genes, and (c) clusters identified via cladistic analysis of phenotypic traits (genitalia and wing pattern characters, ecological preferences, and chromosome numbers). We conclude that in case of the genus Brenthis, nuclear genes rather than mtDNA barcodes reveal real species boundaries. Additionally, we suggest to support each DNA barcode‐based taxonomic conclusion by analysis of phased alleles of nuclear genes, avoiding widely used practice of nuclear and mitochondrial genes concatenation without any examination of interaction of these different types of data.     

Unexpected mosaic distribution of two hybridizing sibling lineages in the teleplanically dispersing snail Stramonita haemastoma suggests unusual postglacial redistribution or cryptic invasion

Molecular approaches have proven efficient to identify cryptic lineages within single taxonomic entities. Sometimes these cryptic lineages maybe previously unreported or unknown invasive taxa. The genetic structure of the marine gastropod Stramonita haemastoma has been examined in the Western Mediterranean and North‐Eastern Atlantic populations with mtDNA COI sequences and three newly developed microsatellite markers. We identified two cryptic lineages, differentially fixed for alternative mtDNA COI haplogroups and significantly differentiated at microsatellite loci. The mosaic distribution of the two lineages is unusual for a warm‐temperate marine invertebrate with a teleplanic larval stage. The Atlantic lineage was unexpectedly observed as a patch enclosed in the north of the Western Mediterranean Sea between eastern Spain and the French Riviera, and the Mediterranean lineage was found in Macronesian Islands. Although cyto‐nuclear disequilibrium is globally maintained, asymmetric introgression occurs in the Spanish region where the two lineages co‐occur in a hybrid zone. A first interpretation of our results is mito‐nuclear discordance in a stable postglacial hybrid zone. Under this hypothesis, though, the location of genetic discontinuities would be unusual among planktonic dispersers. An alternative interpretation is that the Atlantic lineage, also found in Senegal and Venezuela, has been introduced by human activities in the Mediterranean area and is introgressing Mediterranean genes during its propagation, as theoretically expected. This second hypothesis would add an additional example to the growing list of cryptic marine invasions revealed by molecular studies.     

Variation in Phenotype,Parasite Load and Male Competitive Ability across a Cryptic Hybrid Zone

Background

Molecular genetic studies are revealing an increasing number of cryptic lineages or species, which are highly genetically divergent but apparently cannot be distinguished morphologically. This observation gives rise to three important questions: 1) have these cryptic lineages diverged in phenotypic traits that may not be obvious to humans; 2) when cryptic lineages come into secondary contact, what are the evolutionary consequences: stable co-existence, replacement, admixture or differentiation and 3) what processes influence the evolutionary dynamics of these secondary contact zones?

Methodology/Principal Findings

To address these questions, we first tested whether males of the Iberian lizard Lacerta schreiberi from two highly genetically divergent, yet morphologically cryptic lineages on either side of an east-west secondary contact could be differentiated based on detailed analysis of morphology, coloration and parasite load. Next, we tested whether these differences could be driven by pre-copulatory intra-sexual selection (male-male competition). Compared to eastern males, western males had fewer parasites, were in better body condition and were more intensely coloured. Although subtle environmental variation across the hybrid zone could explain the differences in parasite load and body condition, these were uncorrelated with colour expression, suggesting that the differences in coloration reflect heritable divergence. The lineages did not differ in their aggressive behaviour or competitive ability. However, body size, which predicted male aggressiveness, was positively correlated with the colour traits that differed between genetic backgrounds.

Conclusions/Significance

Our study confirms that these cryptic lineages differ in several aspects that are likely to influence fitness. Although there were no clear differences in male competitive ability, our results suggest a potential indirect role for intra-sexual selection. Specifically, if lizards use the colour traits that differ between genetic backgrounds to assess the size of potential rivals or mates, the resulting fitness differential favouring western males could result in net male-mediated gene flow from west to east across the current hybrid zone.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

Endemic or introduced? Phylogeography of Asparagopsis (Florideophyceae) in Australia reveals multiple introductions and a new mitochondrial lineage

The red seaweed Asparagopsis taxiformis embodies five cryptic mitochondrial lineages (lineage 1–5) introduced worldwide as a consequence of human mediated transport and climate change. We compared globally collected mitochondrial cox2‐3 intergenic spacer sequences with sequences produced from multiple Australian locations and South Korea to identify Asparagopsis lineages and to reveal cryptic introductions. We report A. taxiformis lineage 4 from Cocos (Keeling) Islands, Australia, and the highly invasive Indo‐Pacific Mediterranean lineage 2 from South Korea and Lord Howe Island, Australia. Phylogeographic analysis showed a clear haplotype and geographic separation between western Australian and Great Barrier Reef (GBR) isolates belonging to the recently described lineage 5. The same lineage, however, was characterized by a substantial genetic and geographic break between the majority of Australian specimens and Asparagopsis collections from South Solitary Island, Southern GBR, Lord Howe Island, Kermadec Islands, Norfolk Island, New Caledonia and French Polynesia. The disjunct geographic distribution and sequence divergence between these two groups supports the recognition of a sixth cryptic A. taxiformis mitochondrial lineage. As climatic changes accelerate the relocation of biota and offer novel niches for colonization, periodic surveys for early detection of cryptic invasive seaweeds will be critical in determining whether eradication or effective containment of the aliens are feasible.     

Multiple cryptic species of sympatric generalists within the avian blood parasite Haemoproteus majoris

The avian haemosporidian parasite Haemoproteus majoris has been reported to infect a wide range of passerine birds throughout the Holarctic ecozone. Five cytochrome b (cyt b) lineages have been described as belonging to the morphological species H. majoris, and these form a tight phylogenetic cluster together with 13 undescribed lineages that differ from each other by < 1.2% in sequence divergence. Records in a database (MalAvi) that contains global findings of haemosporidian lineages generated by universal primers suggest that these lineages vary substantially in host distribution. We confirm this pattern in a data set collected at Lake Kvismaren, Sweden, where three of the generalist lineages have local transmission. However, whether these lineages represent intraspecific mitochondrial diversity or clusters of cryptic species has previously not been examined. In this study, we developed novel molecular markers to amplify the partial segments of four nuclear genes to determine the level of genetic diversity and gene phylogenies among the five morphologically described cyt b lineages of H. majoris. All five cyt b lineages were strongly associated with unique nuclear alleles at all four nuclear loci, indicating that each mitochondrial lineage represents a distinct biological species. Within lineages, there was no apparent association between nuclear alleles and host species, indicating that they form genetically unstructured populations across multiple host species.     

Phylogenetic analysis of Allium subg. Melanocrommyum infers cryptic species and demands a new sectional classification

Allium subgenus Melanocrommyum (Alliaceae) from Eurasia comprises about 150 mostly diploid species adapted to arid conditions. The group is taxonomically complicated with different and contradictory taxonomic treatments, and was thought to include a considerable number of hybrid species, as the taxa show an admixture of assumed morphological key characters. We studied the phylogeny of the subgenus, covering all existing taxonomic groups and their entire geographic distribution. We analyzed sequences of the nuclear rDNA internal transcribed spacer region (ITS) for multiple individuals of more than 100 species. Phylogenetic analyses of cloned and directly sequenced PCR products confirmed the monophyly of the subgenus, while most sections were either para- or polyphyletic. The splits of the large sections are supported by differences in the anatomy of flower nectaries. ITS data (i) demand a new treatment at sectional level, (ii) do not support the hypotheses of frequent gene flow among species, (iii) indicate that multiple rapid radiations occurred within different monophyletic groups of the subgenus, and (iv) detected separately evolving lineages within three morphologically clearly defined species (cryptic species). In two cases these lineages were close relatives, while in Allium darwasicum they fall in quite different clades in the phylogenetic tree. Fingerprint markers show that this result is not due to ongoing introgression of rDNA (ITS capture) but that genome-wide differences between both lineages exist. Thus, we report one of the rare cases in plants where morphologically indistinguishable diploid species occurring in mixed populations are non-sister cryptic species.     

Strong selection against hybrids maintains a narrow contact zone between morphologically cryptic lineages in a rainforest lizard

Phenotypically cryptic lineages comprise an important yet understudied part of biodiversity; in particular, we have much to learn about how these lineages are formed and maintained. To better understand the evolutionary significance of such lineages, we studied a hybrid zone between two morphologically cryptic phylogeographic lineages in the rainforest lizard, Lampropholis coggeri. Analyzing a multilocus genetic dataset through cline inference, individual-based methods and population measures of disequilibrium and using simulations to explore our genetic results in context of theoretical expectations, we inferred the processes maintaining this hybrid zone. We find that these lineages meet in a hybrid zone that is narrow (≈400 m) relative to inferred dispersal rate. Further, the hybrid zone exhibits substantial genetic disequilibrium and sharply coincident and largely concordant clines. Based on our knowledge about the region's biogeography, the species' natural history, and our simulation results, we suggest that strong selection against hybrids structures this system. As all clines show a relatively narrow range of introgression, we posit that this hybrid zone might not yet be in equilibrium. Nonetheless, our results clearly show that phylogeographic lineages can evolve substantial reproductive isolation without concomitant morphological diversification, suggesting that such lineages can constitute a significant component of evolutionary diversity.     

Deceptive single‐locus taxonomy and phylogeography: Wolbachia‐associated divergence in mitochondrial DNA is not reflected in morphology and nuclear markers in a butterfly species

The satyrine butterfly Coenonympha tullia (Nymphalidae: Satyrinae) displays a deep split between two mitochondrial clades, one restricted to northern Alberta, Canada, and the other found throughout Alberta and across North America. We confirm this deep divide and test hypotheses explaining its phylogeographic structure. Neither genitalia morphology nor nuclear gene sequence supports cryptic species as an explanation, instead indicating differences between nuclear and mitochondrial genome histories. Sex‐biased dispersal is unlikely to cause such mito‐nuclear differences; however, selective sweeps by reproductive parasites could have led to this conflict. About half of the tested samples were infected by Wolbachia bacteria. Using multilocus strain typing for three Wolbachia genes, we show that the divergent mitochondrial clades are associated with two different Wolbachia strains, supporting the hypothesis that the mito‐nuclear differences resulted from selection on the mitochondrial genome due to selective sweeps by Wolbachia strains.     

Cryptic lineages of a common alpine mayfly show strong life‐history divergence

Understanding ecological divergence of morphologically similar but genetically distinct species – previously considered as a single morphospecies – is of key importance in evolutionary ecology and conservation biology. Despite their morphological similarity, cryptic species may have evolved distinct adaptations. If such ecological divergence is unaccounted for, any predictions about their responses to environmental change and biodiversity loss may be biased. We used spatio‐temporally replicated field surveys of larval cohort structure and population genetic analyses (using nuclear microsatellite markers) to test for life‐history divergence between two cryptic lineages of the alpine mayfly Baetis alpinus in the Swiss Alps. We found that the more widespread and abundant cryptic lineage represents a ‘generalist’ with at least two cohorts per year, whereas the less abundant lineage is restricted to higher elevations and represents a ‘specialist’ with a single cohort per year. Importantly, our results indicate partial temporal segregation in reproductive periods between these lineages, potentially facilitating local coexistence and reproductive isolation. Taken together, our findings emphasize the need for a taxonomic revision: widespread and apparently generalist morphospecies can hide cryptic lineages with much narrower ecological niches and distribution ranges.     

DNA barcoding of bark and ambrosia beetles reveals excessive NUMTs and consistent east‐west divergence across Palearctic forests

A comprehensive DNA barcoding library is very useful for rapid identification and detection of invasive pest species. We tested the performance of species identification in the economically most damaging group of wood‐boring insects – the bark and ambrosia beetles – with particular focus on broad geographical sampling across the boreal Palearctic forests. Neighbour‐joining and Bayesian analyses of cytochrome oxidase I (COI) sequences from 151 species in 40 genera revealed high congruence between morphology‐based identification and sequence clusters. Inconsistencies with morphological identifications included the discovery of a likely cryptic Nearctic species of Dryocoetes autographus, the possible hybrid origin of shared mitochondrial haplotypes in Pityophthorus micrographus and P. pityographus, and a possible paraphyletic Xyleborinus saxeseni. The first record of Orthotomicus suturalis in North America was confirmed by DNA barcoding. The mitochondrial data also revealed consistent divergence across the Palearctic or Holarctic, confirmed in part by data from the large ribosomal subunit (28S). Some populations had considerable variation in the mitochondrial barcoding marker, but were invariant in the nuclear ribosomal marker. These findings must be viewed in light of the high number of nuclear insertions of mitochondrial DNA (NUMTs) detected in eight bark beetle species, suggesting the possible presence of additional cryptic NUMTs. The occurrence of paralogous COI copies, hybridization or cryptic speciation demands a stronger focus on data quality assessment in the construction of DNA barcoding databases.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male‐biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal‐biased X‐chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.     

Phylogeny and evolutionary history of birch mice Sicista Griffith, 1827 (Sminthidae,Rodentia): Implications from a multigene study

Phylogeny of birch mice is estimated using sequences of ten nuclear genes and one mitochondrial gene. Based on the results of tree reconstructions and molecular dating, five major lineages are recognized: “tianschanica,” “concolor,” “caudata,” “betulina,” and “caucasica.” It is established that the three latter lineages constitute a clade and that the long‐tailed birch mouse Sicista caudata is the sister group of the “caucasica” lineage. The “tianschanica” lineage is placed as the sister branch to all other species, however, with insufficient support. The cytochrome b tree is generally concordant with the nuclear topology. The molecular clock results suggest that the radiation among the main lineages occurred in the Late Miocene–Early Pliocene (6.0–4.7 Mya). The correspondence between molecular dating and the fossil record is discussed. Based on nuclear data, a high level of divergence between cryptic species in the “tianschanica” lineage is confirmed. Mitochondrial and nuclear data suggest the existence of a potential cryptic species within Sicista strandi.     

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co‐occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next‐generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein‐encoding genes, with the Complex IV genes (coxI‐III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.     

When vicars meet: a narrow contact zone between morphologically cryptic phylogeographic lineages of the rainforest skink, Carlia rubrigularis

Phylogeographic analyses of the fauna of the Australian wet tropics rainforest have provided strong evidence for long-term isolation of populations among allopatric refugia, yet typically there is no corresponding divergence in morphology. This system provides an opportunity to examine the consequences of geographic isolation, independent of morphological divergence, and thus to assess the broader significance of historical subdivisions revealed through mitochondrial DNA phylogeography. We have located and characterized a zone of secondary contact between two long isolated (mtDNA divergence > 15%) lineages of the skink Carlia rubrigularis using one mitochondrial and eight nuclear (two intron, six microsatellite) markers. This revealed a remarkably narrow (width < 3 km) hybrid zone with substantial linkage disequilibrium and strong deficits of heterozygotes at two of three nuclear loci with diagnostic alleles. Cline centers were coincident across loci. Using a novel form of likelihood analysis, we were unable to distinguish between sigmoidal and stepped cline shapes except at one nuclear locus for which the latter was inferred. Given estimated dispersal rates of 90-133 m x gen(-1/2) and assuming equilibrium, the observed cline widths suggest effective selection against heterozygotes of at least 22-49% and possibly as high as 70%. These observations reveal substantial postmating isolation, although the absence of consistent deviations from Hardy-Weinberg equilibrium at diagnostic loci suggests that there is little accompanying premating isolation. The tight geographic correspondence between transitions in mtDNA and those for nuclear genes and corresponding evidence for selection against hybrids indicates that these morphologically cryptic phylogroups could be considered as incipient species. Nonetheless, we caution against the use of mtDNA phylogeography as a sole criterion for defining species boundaries.     

Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Recent speciation and hybridization in Icelandic deep-sea isopods: An integrative approach using genomics and proteomics

The crustacean marine isopod species Haploniscus bicuspis (Sars, 1877) shows circum-Icelandic distribution in a wide range of environmental conditions and along well-known geographic barriers, such as the Greenland-Iceland-Faroe (GIF) Ridge. We wanted to explore population genetics, phylogeography and cryptic speciation as well as investigate whether previously described, but unaccepted subspecies have any merit. Using the same set of specimens, we combined mitochondrial COI sequences, thousands of nuclear loci (ddRAD), and proteomic profiles, plus selected morphological characters using confocal laser scanning microscopy (CLSM). Five divergent genetic lineages were identified by COI and ddRAD, two south and three north of the GIF Ridge. Assignment of populations to the three northern lineages varied and detailed analyses revealed hybridization and gene flow between them, suggesting a single northern species with a complex phylogeographic history. No apparent hybridization was observed among lineages south of the GIF Ridge, inferring the existence of two more species. Differences in proteomic profiles between the three putative species were minimal, implying an ongoing or recent speciation process. Population differentiation was high, even among closely associated populations, and higher in mitochondrial COI than nuclear ddRAD loci. Gene flow is apparently male-biased, leading to hybrid zones and instances of complete exchange of the local nuclear genome through immigrating males. This study did not confirm the existence of subspecies defined by male characters, which probably instead refer to different male developmental stages.