Comparative phylogeography of codistributed aquatic insects in western North America: insights into dispersal and regional patterns of genetic structure

相似文献   

Invertebrate population genetics across Earth's largest habitat: The deep‐sea floor

Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non‐chemosynthetic ecosystems on the deep‐sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow‐water species. Generally, populations at similar depths were well connected over 100s–1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s–1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean‐wide (under 4%), and 48% were Atlantic‐focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, “ecosystem engineers” and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single‐locus mitochondrial genes revealing a common pattern of non‐neutrality, consistent with demographic instability or selective sweeps; similar to deep‐sea hydrothermal vent fauna. The absence of a clear difference between vent and non‐vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single‐locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta‐analyses where broad inferences about deep‐sea ecology could be made.     

Metapopulation patterns of additive and nonadditive genetic variance in the sea bass (Dicentrarchus labrax)

Describing and explaining the geographic within‐species variation in phenotypes (“phenogeography”) in the sea over a species distribution range is central to our understanding of a variety of eco‐evolutionary topics. However, phenogeographic studies that have a large potential to investigate adaptive variation are overcome by phylogeographic studies, still mainly focusing on neutral markers. How genotypic and phenotypic data could covary over large geographic scales remains poorly understood in marine species. We crossed 75 noninbred sires (five origins) and 26 dams (two origins; each side of a hybrid zone) in a factorial diallel cross in order to investigate geographic variation for early survival and sex ratio in the metapopulation of the European sea bass (Dicentrarchus labrax), a highly prized marine fish species. Full‐sib families (N = 1,950) were produced and reared in a common environment. Parentage assignment of 7,200 individuals was performed with seven microsatellite markers. Generalized linear models showed significant additive effects for both traits and pleiotropy between traits. A significant nonadditive genetic effect was detected. Different expression of traits and distinct relative performances were found for reciprocal crosses involving populations located on each side of the main hybrid zone located at the Almeria‐Oran front, illustrating asymmetric reproductive isolation. The poor fitness performance observed for the Western Mediterranean population of sea bass is discussed as it represents the main source of seed hatchery production, but also because it potentially illustrates nonadaptive introgression and maladaptation.     

Is reproductive strategy a key factor in understanding the evolutionary history of Southern Ocean Asteroidea (Echinodermata)?

Life traits such as reproductive strategy can be determining factors of species evolutionary history and explain the resulting diversity patterns. This can be investigated using phylogeographic analyses of genetic units. In this work, the genetic structure of five asteroid genera with contrasting reproductive strategies (brooding: Diplasterias, Notasterias and Lysasterias versus broadcasting: Psilaster and Bathybiaster) was investigated in the Southern Ocean. Over 1,400 mtDNA cytochrome C oxidase subunit I (COI) sequences were analysed using five species delineation methods (ABGD, ASAP, mPTP, sGMYC and mGMYC), two phylogenetic reconstructions (ML and BA), and molecular clock calibrations, in order to examine the weight of reproductive strategy in the observed differences among phylogeographic patterns. We hypothesised that brooding species would show higher levels of genetic diversity and species richness along with a clearer geographic structuring than broadcasting species. In contrast, genetic diversity and species richness were not found to be significantly different between brooders and broadcasters, but broadcasters are less spatially structured than brooders supporting our initial hypothesis and suggesting more complex evolutionary histories associated to this reproductive strategy. Broadcasters' phylogeography can be explained by different scenarios including deep‐sea colonisation routes, bipolarity or cosmopolitanism, and sub‐Antarctic emergence for the genus Bathybiaster; Antarctic‐ New Zealand faunal exchanges across the Polar Front for the genus Psilaster. Brooders' phylogeography could support the previously formulated hypothesis of a past trans‐Antarctic seaway established between the Ross and the Weddell seas during the Plio‐Pleistocene. Our results also show, for the first time, that the Weddell Sea is populated by a mixed asteroid fauna originating from both the East and West Antarctic.     

Recent phylogeographic structure in a widespread ‘weedy’ Neotropical tree species,Cordia alliodora (Boraginaceae)

Aim Although hundreds of tree species have broad geographic ranges in the Neotropics, little is known about how such widespread species attained disjunct distributions around mountain, ocean and xeric barriers. Here, we examine the phylogeographic structure of a widespread and economically important tree, Cordia alliodora, to: (1) test the roles of vicariance and dispersal in establishing major range disjunctions, (2) determine which geographic regions and/or habitats contain the highest levels of genetic diversity, and (3) infer the geographic origin of the species. Location Twenty‐five countries in Central and South America, and the West Indies. Methods Chloroplast simple sequence repeats (cpSSR; eight loci) were assayed in 67 populations (240 individuals) sampled from the full geographic range of C. alliodora. Chloroplast (trnH–psbA) and nuclear (internal transcribed spacer, ITS) DNA sequences were sampled from a geographically representative subset. Genetic structure was determined with samova , structure and haplotype networks. Analysis of molecular variance (AMOVA) and rarefaction analyses were used to compare regional haplotype diversity and differentiation. Results Although the ITS region was polymorphic it revealed limited phylogeographic structure, and trnH–psbA was monomorphic. However, structure analysis of cpSSR variation recovered three broad demes spanning Central America (Deme 1), the Greater Antilles and the Chocó (Deme 2), and the Lesser Antilles and cis‐Andean South America (Deme 3). samova showed two predominant demes (Deme 1 + 2 and Deme 3). The greatest haplotype diversity was detected east of the Andes, while significantly more genetic variation was partitioned among trans‐Andean populations. Populations experiencing high precipitation seasonality (dry ecotype) had greater levels of genetic variation. Main conclusions Cordia alliodora displayed weak cis‐ and trans‐Andean phylogeographic structure based on DNA sequence data, indicative of historical dispersal around this barrier and genetic exchange across its broad range. The cpSSR data revealed phylogeographic structure corresponding to three biogeographic zones. Patterns of genetic diversity are indicative of an origin in the seasonally dry habitats of South America. Therefore, C. alliodora fits the disperser hypothesis for widespread Neotropical species. Dispersal is evident in the West Indies and the northern Andean cordilleras. The dry ecotype harbours genetic variation that is likely to represent the source for the establishment of populations under future warmer and drier climatic scenarios.     

A rangewide herbarium‐derived dataset indicates high levels of gene flow in black cherry (Prunus serotina)

Isolation by Distance (IBD) is a genetic pattern in which populations geographically closer to one another are more genetically similar to each other than populations which are farther apart. Black cherry (Prunus serotina Ehrh.) (Rosaceae) is a forest tree species widespread in eastern North America, and found sporadically in the southwestern United States, Mexico, and Guatemala. IBD has been studied in relatively few North American plant taxa, and no study has rigorously sampled across the range of such a widespread species. In this study, IBD and overall genetic structure were assessed in eastern black cherry (P. serotina Ehrh. var. serotina), the widespread variety of eastern North America. Eastern North America. Prunus serotina Ehrh. var. serotina (Rosaceae). Dense sampling across the entire range of eastern black cherry was made possible by genotyping 15 microsatellite loci in 439 herbarium samples from all portions of the range. Mantel tests and STRUCTURE analyses were performed to evaluate the hypothesis of IBD and genetic structure. Mantel tests demonstrated significant but weak IBD, while STRUCTURE analyses revealed no clear geographic pattern of genetic groups. The modest geographic/genetic structure across the eastern black cherry range suggests widespread gene flow in this taxon. This is consistent with P. serotina's status as a disturbance‐associated species. Further studies should similarly evaluate IBD in species characteristic of low‐disturbance forests.     

Phylogeographic inference using Bayesian model comparison across a fragmented chorus frog species complex

Fragmented species complexes provide an interesting system for investigating biogeographic history and the present distribution of genetic variation. Recent advances in sequencing technology and statistical phylogeography enable the collection and rigorous analysis of large multilocus data sets, but designing studies that produce meaningful phylogeographic inferences remains challenging. We implemented a Bayesian model comparison approach to investigate previous biogeographic hypotheses while simultaneously inferring the presence of genetic structure in a chorus frog species complex. The Illinois chorus frog (Pseudacris illinoensis), originally described as a subspecies of the broadly distributed Strecker's chorus frog (Pseudacris streckeri), occurs in small, disjunct regions associated with scarce sand prairie habitats that have been impacted by human development. We used high‐throughput sequencing to develop and collect a multitiered genetic data set comprised of three different marker types (23 anonymous nuclear sequence loci, four mitochondrial genes and 14 microsatellite loci) designed to address questions across different evolutionary timescales. Phylogenetic analyses uncovered a deep divergence between populations in the Edwards Plateau of central Texas and all other P. streckeri/P. illinoensis populations, but suggest the disjunct distribution of P. illinoensis occurred more recently. Our best‐supported migration model is consistent with the hypothesis that central Texas represented a refugium from which populations expanded via multiple routes. This model also indicates that disjunct northern and southern regions of P. illinoensis should be considered genetically distinct management units. Our study provides an evolutionary context for future studies and conservation efforts in P. illinoensis and demonstrates the utility of model‐based approaches for phylogeographic inference.     

Biogeography and host‐related factors trump parasite life history: limited congruence among the genetic structures of specific ectoparasitic lice and their rodent hosts

Parasites and hosts interact across both micro‐ and macroevolutionary scales where congruence among their phylogeographic and phylogenetic structures may be observed. Within southern Africa, the four‐striped mouse genus, Rhabdomys, is parasitized by the ectoparasitic sucking louse, Polyplax arvicanthis. Molecular data recently suggested the presence of two cryptic species within P. arvicanthis that are sympatrically distributed across the distributions of four putative Rhabdomys species. We tested the hypotheses of phylogeographic congruence and cophylogeny among the two parasite lineages and the four host taxa, utilizing mitochondrial and nuclear sequence data. Despite the documented host‐specificity of P. arvicanthis, limited phylogeographic correspondence and nonsignificant cophylogeny was observed. Instead, the parasite–host evolutionary history is characterized by limited codivergence and several duplication, sorting and host‐switching events. Despite the elevated mutational rates found for P. arvicanthis, the spatial genetic structure was not more pronounced in the parasite lineages compared with the hosts. These findings may be partly attributed to larger effective population sizes of the parasite lineages, the vagility and social behaviour of Rhabdomys, and the lack of host‐specificity observed in areas of host sympatry. Further, the patterns of genetic divergence within parasite and host lineages may also be largely attributed to historical biogeographic changes (expansion‐contraction cycles). It is thus evident that the association between P. arvicanthis and Rhabdomys has been shaped by the synergistic effects of parasite traits, host‐related factors and biogeography over evolutionary time.     

Genetic diversity of Cedrela fissilis (Meliaceae) in the Brazilian Atlantic Forest reveals a complex phylogeographic history driven by Quaternary climatic fluctuations

Quaternary climatic fluctuations have shaped the geographic distribution of lineages, potentially affecting the demography, genetic structure, and patterns of genetic diversity of extant species. Different phylogeographic scenarios have been proposed for plants in neotropical cloud forests during the Last Glacial Maximum based on paleoecological data: the dry refugia hypothesis (DRH) and the moist forest hypothesis. We specifically focus on the Brazilian Atlantic Forest (BAF) range of Cedrela fissilis (Meliaceae), sampling 410 specimens from 50 localities. Our study combines analyses of the genetic diversity, phylogeographic patterns, and past geographic distributions with a particular focus on highland populations. We identified 283 alleles across the 11 microsatellite loci, ranging from 18 to 33 alleles per locus, distributed across five genetic groups. Most populations of C. fissilis from the BAF exhibited a diffuse genetic structure, reflected in low pairwise F_ST values, which could be the consequence of high gene flow. In addition, the plastid data showed a connection between the western, southern, and eastern populations in the North‐East of Brazil, but no association between genetic data and elevation was observed. Habitat suitability projections over the past 140 000 years showed less fragmentation relative to the present, indicating a higher connectivity and gene flow. Our results provide support for both the moist forest as well as the DRH, suggesting that most likely, a mixture of these processes has acted through space and time.     

Genotyping‐by‐sequencing and ecological niche modeling illuminate phylogeography,admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest‐ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome‐wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal–Cascades (cluster 1), east‐slope Cascades–Sierra Nevadas–Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the “stable‐edge hypothesis” that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited “trailing‐edge” dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the “inland dispersal hypothesis” predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific‐coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable‐edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific‐coastal genetic lineage of quaking aspen.     

Phylogeography and ecological niche modeling unravel the evolutionary history of the African green toad,Bufotes boulengeri boulengeri (Amphibia: Bufonidae), through the Quaternary

Recent integration of ecological niche models in phylogeographic studies is improving our understanding of the processes structuring genetic variation across landscapes. Previous studies on the amphibian Bufotes boulengeri boulengeri uncovered a surprisingly weak intraspecific differentiation across the Maghreb region. We widely sampled this species from Morocco to Egypt and sequenced one nuclear and three mitochondrial (mtDNA) genes to determine the level of genetic variability across its geographic range. We evaluated these data with ecological niche modeling to reveal its evolutionary history in response to climate change during the Quaternary. Our results highlight some mtDNA phylogeographic structure within this species, with one haplogroup endemic to coastal Morocco, and one haplogroup widely distributed throughout North Africa. No or little genetic differentiation is observed between isolated populations from the Hoggar Mountains, the Sabha district and the islands of Kerkennah and Lampedusa, compared to others populations. This can be explained by the expansion of the distribution range of B. b. boulengeri during glacial periods. This might have facilitated the species’ dispersal and subsequent gene flow between most North African localities.     

Three‐dimensional post‐glacial expansion and diversification of an exploited oceanic fish

Vertical divergence in marine organisms is being increasingly documented, yet much remains to be carried out to understand the role of depth in the context of phylogeographic reconstruction and the identification of management units. An ideal study system to address this issue is the beaked redfish, Sebastes mentella – one of four species of ‘redfish’ occurring in the North Atlantic – which is known for a widely distributed ‘shallow‐pelagic’ oceanic type inhabiting waters between 250 and 550 m, and a more localized ‘deep‐pelagic’ population dwelling between 550 and 800 m, in the oceanic habitat of the Irminger Sea. Here, we investigate the extent of population structure in relation to both depth and geographic spread of oceanic beaked redfish throughout most of its distribution range. By sequencing the mitochondrial control region of 261 redfish collected over a decadal interval, and combining 160 rhodopsin coding nuclear sequences and previously genotyped microsatellite data, we map the existence of two strongly divergent evolutionary lineages with significantly different distribution patterns and historical demography, and whose genetic variance is mostly explained by depth. Combined genetic data, analysed via independent approaches, are consistent with a Late Pleistocene lineage split, where segregation by depth probably resulted from the interplay of climatic and oceanographic processes with life history and behavioural traits. The ongoing process of diversification in North Atlantic S. mentella may serve as an ‘hourglass’ to understand speciation and adaptive radiation in Sebastes and in other marine taxa distributed across a depth gradient.     

Inter‐oceanic variation in patterns of host‐associated divergence in a seabird ectoparasite

Aim Parasites with global distributions and wide host spectra provide excellent models for exploring the factors that drive parasite diversification. Here, we tested the relative force of host and geography in shaping population structure of a widely distributed and common ectoparasite of colonial seabirds, the tick Ixodes uriae. Location Two natural geographic replicates of the system: numerous seabird colonies of the North Pacific and North Atlantic Ocean basins. Methods Using eight microsatellite markers and tick samples from a suite of multi‐specific seabird colonies, we examined tick population structure in the North Pacific and compare patterns of diversity and structure to those in the Atlantic basin. Analyses included population genetic estimations of diversity and population differentiation, exploratory multivariate analyses, and Bayesian clustering approaches. These different analyses explicitly took into account both the geographic distance among colonies and host use by the tick. Results Overall, little geographic structure was observed among Pacific tick populations. However, host‐related genetic differentiation was evident, but was variable among host types and lower than in the North Atlantic. Main conclusions Tick population structure is concordant with the genetic structure observed in seabird host species within each ocean basin, where seabird populations tend to be less structured in the North Pacific than in the North Atlantic. Reduced tick genetic structure in the North Pacific suggests that host movement among colonies, and thus tick dispersal, is higher in this region. In addition to information on parasite diversity and gene flow, our findings raise interesting questions about the subtle ways that host behaviour, distribution and phylogeographic history shape the genetics of associated parasites across geographic landscapes.     

The presence of a cryptic barrier in the West Pacific Ocean suggests the effect of glacial climate changes on a widespread sea‐dispersed plant,Vigna marina (Fabaceae)

Ocean currents are an important driver of evolution for sea‐dispersed plants, enabling them to maintain reciprocal gene flow via sea‐dispersed diaspores and obtain wide distribution ranges. Although geographic barriers are known to be the primary factors shaping present genetic structure of sea‐dispersed plants, cryptic barriers which form clear genetic structure within oceanic regions are poorly understood. To test the presence of a cryptic barrier, we conducted a phylogeographic study together with past demographic inference for a widespread sea‐dispersed plant, Vigna marina, using 308 individuals collected from the entire Indo‐West Pacific (IWP) region. Chloroplast DNA variation showed strong genetic structure that separated populations into three groups: North Pacific (NP), South Pacific (SP) and Indian Ocean (IN) (F′_CT among groups = 0.954–1.000). According to the Approximate Bayesian computation inference, splitting time between NP and SP was approximately 20,200 years (95%HPD, 4,530–95,400) before present. Moreover, a signal of recent population expansion was detected in the NP group. This study clearly showed the presence of a cryptic barrier in the West Pacific region of the distributional range of V. marina. The locations of the cryptic barrier observed in V. marina corresponded to the genetic breaks found in other plants, suggesting the presence of a common cryptic barrier for sea‐dispersed plants. Demographic inference suggested that genetic structure related to this cryptic barrier has been present since the last glacial maximum and may reflect patterns of past population expansion from refugia.     

Land masses and oceanic currents drive population structure of Heritiera littoralis,a widespread mangrove in the Indo‐West Pacific

Phylogeographic forces driving evolution of sea‐dispersed plants are often influenced by regional and species characteristics, although not yet deciphered at a large spatial scale for many taxa like the mangrove species Heritiera littoralis. This study aimed to assess geographic distribution of genetic variation of this widespread mangrove in the Indo‐West Pacific region and identify the phylogeographic factors influencing its present‐day distribution. Analysis of five chloroplast DNA fragments’ sequences from 37 populations revealed low genetic diversity at the population level and strong genetic structure of H. littoralis in this region. The estimated divergence times between the major genetic lineages indicated that glacial level changes during the Pleistocene epoch induced strong genetic differentiation across the Indian and Pacific Oceans. In comparison to the strong genetic break imposed by the Sunda Shelf toward splitting the lineages of the Indian and Pacific Oceans, the genetic differentiation between Indo‐Malesia and Australasia was not so prominent. Long‐distance dispersal ability of H. littoralis propagules helped the species to attain transoceanic distribution not only across South East Asia and Australia, but also across the Indian Ocean to East Africa. However, oceanic circulation pattern in the South China Sea was found to act as a barrier creating further intraoceanic genetic differentiation. Overall, phylogeographic analysis in this study revealed that glacial vicariance had profound influence on population differentiation in H. littoralis and caused low genetic diversity except for the refugia populations near the equator which might have persisted through glacial maxima. With increasing loss of suitable habitats due to anthropogenic activities, these findings therefore emphasize the urgent need for conservation actions for all populations throughout the distribution range of H. littoralis.     

Phylogenetic placement of Cibicidoides wuellerstorfi (Schwager, 1866) from methane seeps and non‐seep habitats on the Pacific margin

Benthic foraminifera are among the most abundant groups found in deep‐sea habitats, including methane seep environments. Unlike many groups, no endemic foraminiferal species have been reported from methane seeps, and to our knowledge, genetic data are currently sparse for Pacific deep‐sea foraminifera. In an effort to understand the relationships between seep and non‐seep populations of the deep‐sea foraminifera Cibicidoides wuellerstorfi, a common paleo‐indicator species, specimens from methane seeps in the Pacific were analyzed and compared to one another for genetic similarities of small subunit rDNA (SSU rDNA) sequences. Pacific Ocean C. wuellerstorfi were also compared to those collected from other localities around the world (based on 18S gene available on Genbank, e.g., Schweizer et al., 2009). Results from this study revealed that C. wuellerstorfi living in seeps near Costa Rica and Hydrate Ridge are genetically similar to one another at the species level. Individuals collected from the same location that display opposite coiling directions (dextral and sinstral) had no species level genetic differences. Comparisons of specimens with genetic information available from Genbank (SSU rDNA) showed that Pacific individuals, collected for this study, are genetically similar to those previously analyzed from the North Atlantic and Antarctic. These observations provide strong evidence for the true cosmopolitan nature of C. wuellerstorfi and highlight the importance of understanding how these microscopic organisms are able to maintain sufficient genetic exchange to remain within the same species between seep and non‐seep habitats and over global distances.     

Deep phylogeographic divisions and long‐term persistence of forest invertebrates (Hexapoda: Collembola) in the North‐Western Mediterranean basin

The North‐Western Mediterranean basin is well known for its high number of relictual endemic taxa, and has been indicated as one of the world’s major biodiversity hotspots at the species level. A possible contributing factor may be long‐term persistence of populations and their prolonged stability. This study was designed to investigate the phylogeographic structure of three common species of the genus Lepidocyrtus (Hexapoda: Collembola), soil‐dwelling arthropods characterized by limited dispersal capabilities and generally associated with forest habitats. We observed a remarkable geographic structure, with numerous deeply divergent genetic lineages occupying islands as well as mainland sites with no apparent gene flow among most sites, even across distances of only tens of kilometres. The reconstructed time frame for the evolution of these lineages suggests divergence between 5 and 15 Ma. This indicates a remarkably ancient origin and long‐term persistence of individual lineages over a fine geographic scale despite the occurrence of abrupt sea level and climatic fluctuations in the area. This further suggests that currently recognized morphological species might be a serious underestimation of the true springtail biodiversity within this region.     

Amplicon sequencing of 42 nuclear loci supports directional gene flow between South Pacific populations of a hydrothermal vent limpet

In the past few decades, population genetics and phylogeographic studies have improved our knowledge of connectivity and population demography in marine environments. Studies of deep‐sea hydrothermal vent populations have identified barriers to gene flow, hybrid zones, and demographic events, such as historical population expansions and contractions. These deep‐sea studies, however, used few loci, which limit the amount of information they provided for coalescent analysis and thus our ability to confidently test complex population dynamics scenarios. In this study, we investigated population structure, demographic history, and gene flow directionality among four Western Pacific hydrothermal vent populations of the vent limpet Lepetodrilus aff. schrolli. These vent sites are located in the Manus and Lau back‐arc basins, currently of great interest for deep‐sea mineral extraction. A total of 42 loci were sequenced from each individual using high‐throughput amplicon sequencing. Amplicon sequences were analyzed using both genetic variant clustering methods and evolutionary coalescent approaches. Like most previously investigated vent species in the South Pacific, L. aff. schrolli showed no genetic structure within basins but significant differentiation between basins. We inferred significant directional gene flow from Manus Basin to Lau Basin, with low to no gene flow in the opposite direction. This study is one of the very few marine population studies using >10 loci for coalescent analysis and serves as a guide for future marine population studies.     

Host hybridization as a potential mechanism of lateral symbiont transfer in deep‐sea vesicomyid clams

Deep‐sea vesicomyid clams live in mutualistic symbiosis with chemosynthetic bacteria that are inherited through the maternal germ line. On evolutionary timescales, strictly vertical transmission should lead to cospeciation of host mitochondrial and symbiont lineages; nonetheless, examples of incongruent phylogenies have been reported, suggesting that symbionts are occasionally horizontally transmitted between host species. The current paradigm for vesicomyid clams holds that direct transfers cause host shifts or mixtures of symbionts. An alternative hypothesis suggests that hybridization between host species might explain symbiont transfers. Two clam species, Archivesica gigas and Phreagena soyoae, frequently co‐occur at deep‐sea hydrocarbon seeps in the eastern Pacific Ocean. Although the two species typically host gammaproteobacterial symbiont lineages marked by divergent 16S rRNA phylotypes, we identified a number of clams with the A. gigas mitotype that hosted symbionts with the P. soyoae phylotype. Demographic inference models based on genome‐wide SNP data and three Sanger sequenced gene markers provided evidence that A. gigas and P. soyoae hybridized in the past, supporting the hypothesis that hybridization might be a viable mechanism of interspecific symbiont transfer. These findings provide new perspectives on the evolution of vertically transmitted symbionts and their hosts in deep‐sea chemosynthetic environments.