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.     

Mitochondrial phylogeography of a Beringian relict: the endemic freshwater genus of blackfish Dallia (Esociformes)

Mitochondrial genetic variability among populations of the blackfish genus Dallia (Esociformes) across Beringia was examined. Levels of divergence and patterns of geographic distribution of mitochondrial DNA lineages were characterized using phylogenetic inference, median‐joining haplotype networks, Bayesian skyline plots, mismatch analysis and spatial analysis of molecular variance (SAMOVA) to infer genealogical relationships and to assess patterns of phylogeography among extant mitochondrial lineages in populations of species of Dallia. The observed variation includes extensive standing mitochondrial genetic diversity and patterns of distinct spatial segregation corresponding to historical and contemporary barriers with minimal or no mixing of mitochondrial haplotypes between geographic areas. Mitochondrial diversity is highest in the common delta formed by the Yukon and Kuskokwim Rivers where they meet the Bering Sea. Other regions sampled in this study host comparatively low levels of mitochondrial diversity. The observed levels of mitochondrial diversity and the spatial distribution of that diversity are consistent with persistence of mitochondrial lineages in multiple refugia through the last glacial maximum.     

Host genetic determinants of the gut microbiota of wild mice

Identifying a common set of genes that mediate host–microbial interactions across populations and species of mammals has broad relevance for human health and animal biology. However, the genetic basis of the gut microbial composition in natural populations remains largely unknown outside of humans. Here, we used wild house mouse populations as a model system to ask three major questions: (a) Does host genetic relatedness explain interindividual variation in gut microbial composition? (b) Do population differences in the microbiota persist in a common environment? (c) What are the host genes associated with microbial richness and the relative abundance of bacterial genera? We found that host genetic distance is a strong predictor of the gut microbial composition as characterized by 16S amplicon sequencing. Using a common garden approach, we then identified differences in microbial composition between populations that persisted in a shared laboratory environment. Finally, we used exome sequencing to associate host genetic variants with microbial diversity and relative abundance of microbial taxa in wild mice. We identified 20 genes that were associated with microbial diversity or abundance including a macrophage‐derived cytokine (IL12a) that contained three nonsynonymous mutations. Surprisingly, we found a significant overrepresentation of candidate genes that were previously associated with microbial measurements in humans. The homologous genes that overlapped between wild mice and humans included genes that have been associated with traits related to host immunity and obesity in humans. Gene–bacteria associations identified in both humans and wild mice suggest some commonality to the host genetic determinants of gut microbial composition across mammals.     

Molecular biogeography and host relations of a parasitoid fly

Successful geographic range expansion by parasites and parasitoids may also require host range expansion. Thus, the evolutionary advantages of host specialization may trade off against the ability to exploit new host species encountered in new geographic regions. Here, we use molecular techniques and confirmed host records to examine biogeography, population divergence, and host flexibility of the parasitoid fly, Ormia ochracea (Bigot). Gravid females of this fly find their cricket hosts acoustically by eavesdropping on male cricket calling songs; these songs vary greatly among the known host species of crickets. Using both nuclear and mitochondrial genetic markers, we (a) describe the geographical distribution and subdivision of genetic variation in O. ochracea from across the continental United States, the Mexican states of Sonora and Oaxaca, and populations introduced to Hawaii; (b) demonstrate that the distribution of genetic variation among fly populations is consistent with a single widespread species with regional host specialization, rather than locally differentiated cryptic species; (c) identify the more‐probable source populations for the flies introduced to the Hawaiian islands; (d) examine genetic variation and substructure within Hawaii; (e) show that among‐population geographic, genetic, and host song distances are all correlated; and (f) discuss specialization and lability in host‐finding behavior in light of the diversity of cricket songs serving as host cues in different geographically separate populations.     

Origins of clonal diversity in the hypervariable asexual ostracode Cypridopsis vidua

Prior allozyme studies have indicated that populations of the asexual ostracode, Cypridopsis vidua (Müller), show extraordinary clonal diversity. Based on a joint examination of allozyme variation and sequence divergence at the mitochondrial cytochrome c oxidase I (COI) gene, the present analysis provides new insights concerning the origins of this variation. The results establish that populations of C. vidua in one recently deglaciated region of North America are not only allozymically diverse, but also include several divergent mitochondrial DNA (mtDNA) lineages. The extent of sequence divergence among these lineages is so large as to suggest their diversification over the past 7–8 million years. The patterning of genetic divergence among co-occurring clones makes it apparent that much of the mtDNA and allozyme diversity in local populations owes its origins to recurrent colonization events. However, in situ mutational diversification also appears to explain some variation. The mechanisms enabling the sustained coexistence of such a large array of closely allied genotypes remain unclear, but there is an apparent difference in equilibrium diversity between benthic and planktonic asexual organisms.     

Delimiting shades of gray: phylogeography of the Northern Fulmar,Fulmarus glacialis

The Northern Fulmar (Fulmarus glacialis) is a common tube‐nosed seabird with a disjunct Holarctic range. Taxonomic divisions within the Northern Fulmar have historically been muddled by geographical variation notably including highly polymorphic plumage. Recent molecular analyses (i.e., DNA barcoding) have suggested that genetic divergence between Atlantic and Pacific populations could be on par with those typically observed between species. We employ a multigene phylogenetic analysis to better explore the level of genetic divergence between these populations and to test an old hypothesis on the origin of the modern distribution of color morphs across their range. Additionally, we test whether mutations in the melanocortin‐1 receptor gene (MC1R) are associated with dark plumage in the Northern Fulmar. We confirmed that mitochondrial lineages in the Atlantic and Pacific populations are highly divergent, but nuclear markers revealed incomplete lineage sorting. Genetic divergence between these populations is consistent with that observed between many species of Procellariiformes and we recommend elevating these two forms to separate species. We also find that MC1R variation is not associated with color morph but rather is better explained by geographical divergence.     

Mitochondrial DNA diversity in an apomictic Daphnia complex from the Canadian High Arctic

Cyclic parthenogenesis is the ancestral mode of reproduction in the cladoceran crustacean, Daphnia pulex, but some populations have made the transition to obligate parthenogenesis and this is the only mode of reproduction known to occur in arctic populations. Melanism and polyploidy are also common in arctic populations of this species. Prior allozyme studies of arctic D. pulex revealed substantial levels of clonal diversity on a regional scale. Clonal groupings based on cluster analysis of allozyme genotypes do not conform to groupings based on the presence/absence of melanin or on ploidy level. In order to further elucidate genetic relationships among arctic D. pulex clones, mitochondrial DNA (mtDNA) variation was examined in 31 populations from two Canadian high-arctic sites. The data were also compared to a previous study of mtDNA variation in populations from a Canadian low-arctic site. Cladistic analysis of restriction site variation of the entire mitochondrial genome and nucleotide sequence variation of the mitochondrial control region was used to construct genetic relationships among mitochondrial genotypes. Three distinct mitochondrial lineages were detected. One lineage was associated with diploid, nonmelanic clones and is the same as the lineage that is found in temperate populations of D. pulex. The other two lineages (A & B) were associated with polyploid, melanic clones. Sequence divergence between the A and B lineages was 2.4%. Sequence divergence between D. pulex and either of these two lineages exceeded 3%. It is suggested that the melanic, polyploid clones are hybrids between males of D. pulex (and/or a closely related congener, D. pulicaria) and females of either of two ancestral melanic species that have mitochondrial lineages A and B. Geographic patterns of mitochondrial diversity in ‘melanic’ lineage B support the hypothesis of an high-arctic refuge for the ancestral species during the last glacial period.     

Major Histocompatibility Complex class IIb polymorphism influences gut microbiota composition and diversity

Animals harbour diverse communities of symbiotic bacteria, which differ dramatically among host individuals. This heterogeneity poses an immunological challenge: distinguishing between mutualistic and pathogenic members of diverse and host‐specific microbial communities. We propose that Major Histocompatibility class II (MHC) genotypes contribute to recognition and regulation of gut microbes, and thus, MHC polymorphism contributes to microbial variation among hosts. Here, we show that MHC IIb polymorphism is associated with among‐individual variation in gut microbiota within a single wild vertebrate population of a small fish, the threespine stickleback. We sampled stickleback from Cedar Lake, on Vancouver Island, and used next‐generation sequencing to genotype the sticklebacks’ gut microbiota (16S sequencing) and their MHC class IIb exon 2 sequences. The presence of certain MHC motifs was associated with altered relative abundance (increase or decrease) of some microbial Families. The effect sizes are modest and entail a minority of microbial taxa, but these results represent the first indication that MHC genotype may affect gut microbiota composition in natural populations (MHC‐microbe associations have also been found in a few studies of lab mice). Surprisingly, these MHC effects were frequently sex‐dependent. Finally, hosts with more diverse MHC motifs had less diverse gut microbiota. One implication is that MHC might influence the efficacy of therapeutic strategies to treat dysbiosis‐associated disease, including the outcome of microbial transplants between healthy and diseased patients. We also speculate that macroparasite‐driven selection on MHC has the potential to indirectly alter the host gut microbiota, and vice versa.     

Clonal diversity driven by parasitism in a freshwater snail

One explanation for the widespread abundance of sexual reproduction is the advantage that genetically diverse sexual lineages have under strong pressure from virulent coevolving parasites. Such parasites are believed to track common asexual host genotypes, resulting in negative frequency‐dependent selection that counterbalances the population growth‐rate advantage of asexuals in comparison with sexuals. In the face of genetically diverse asexual lineages, this advantage of sexual reproduction might be eroded, and instead sexual populations would be replaced by diverse assemblages of clonal lineages. We investigated whether parasite‐mediated selection promotes clonal diversity in 22 natural populations of the freshwater snail Melanoides tuberculata. We found that infection prevalence explains the observed variation in the clonal diversity of M. tuberculata populations, whereas no such relationship was found between infection prevalence and male frequency. Clonal diversity and male frequency were independent of snail population density. Incorporating ecological factors such as presence/absence of fish, habitat geography and habitat type did not improve the predictive power of regression models. Approximately 11% of the clonal snail genotypes were shared among 2–4 populations, creating a web of 17 interconnected populations. Taken together, our study suggests that parasite‐mediated selection coupled with host dispersal ecology promotes clonal diversity. This, in return, may erode the advantage of sexual reproduction in M. tuberculata populations.     

Patterns of morphological variation among populations of the widespread annual killifish Nothobranchius orthonotus are independent of genetic divergence and biogeography

Populations of annual killifish of the genus Nothobranchius occur in patchily distributed temporary pools in the East African savannah. Their fragmented distribution and low dispersal ability result in highly structured genetic clustering of their populations. In this study, we examined body shape variation in a widely distributed species, Nothobranchius orthonotus with known phylogeographic structure. We tested whether genetic divergence of major mitochondrial lineages forming two candidate species is congruent with phenotypic diversification, using linear and geometric morphometry analyses of body shape in 23 wild populations. We also conducted a common‐garden experiment with two wild‐derived populations to control for the effect of local environmental conditions on body shape. We identified different allometric trajectories for different mitochondrial lineages and candidate species in both sexes. However, in a principal components analysis of population‐level body shape, the separation among mitochondrial lineages was incomplete. Higher similarity of mitochondrial lineages belonging to different candidate species than that of same candidate species prevented distinction of the two candidate species on the basis of body shape. Analysis at the individual level demonstrated that N. orthonotus express high intrapopulation variability, with major overlap among individuals from all populations. In conclusion, we suggest that N. orthonotus be considered as a single species with an extensive geographic range, strong population genetic structure and high morphological variability.     

Kinship,inbreeding and fine‐scale spatial structure influence gut microbiota in a hindgut‐fermenting tortoise

Herbivorous vertebrates rely on complex communities of mutualistic gut bacteria to facilitate the digestion of celluloses and hemicelluloses. Gut microbes are often convergent based on diet and gut morphology across a phylogenetically diverse group of mammals. However, little is known about microbial communities of herbivorous hindgut‐fermenting reptiles. Here, we investigate how factors at the individual level might constrain the composition of gut microbes in an obligate herbivorous reptile. Using multiplexed 16S rRNA gene sequencing, we characterized the faecal microbial community of a population of gopher tortoises (Gopherus polyphemus) and examined how age, genetic diversity, spatial structure and kinship influence differences among individuals. We recovered phylotypes associated with known cellulolytic function, including candidate phylum Termite Group 3, suggesting their importance for gopher tortoise digestion. Although host genetic structure did not explain variation in microbial composition and community structure, we found that fine‐scale spatial structure, inbreeding, degree of relatedness and possibly ontogeny shaped patterns of diversity in faecal microbiomes of gopher tortoises. Our findings corroborate widespread convergence of faecal‐associated microbes based on gut morphology and diet and demonstrate the role of spatial and demographic structure in driving differentiation of gut microbiota in natural populations.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Temporal and spatial mosaics: deep host association and shallow geographic drivers shape genetic structure in a widespread pinworm,Rauschtineria eutamii (Nematoda: Oxyuridae)

Climate and host demographic cycling often shape both parasite genetic diversity and host distributions, processes that transcend a history of strict host–parasite association. We explored host associations and histories based on an evaluation of mitochondrial and nuclear sequences to reveal the underlying history and genetic structure of a pinworm, Rauschtineria eutamii, infecting ten species of western North American chipmunks (Rodentia:Tamias, subgenus Neotamias). Rauschtineria eutamii contains divergent lineages influenced by the diversity of hosts and variation across the complex topography of western North America. We recovered six reciprocally monophyletic R. eutamii mitochondrial clades, largely supported by a multilocus concordance tree, exhibiting divergence levels comparable with intraspecific variation reported for other nematodes. Phylogenetic relationships among pinworm clades suggest that R. eutamii colonized an ancestral lineage of western chipmunks and lineages persisted during historical isolation in diverging Neotamias species or species groups. Pinworm diversification, however, is incongruent and asynchronous relative to host diversification. Secondarily, patterns of shallow divergence were shaped by geography through events of episodic colonization reflecting an interaction of taxon pulses and ecological fitting among assemblages in recurrent sympatry. Pinworms occasionally infect geographically proximal host species; however, host switching may be unstable or ephemeral, as there is no signal of host switching in the deeper history of R. eutamii.     

Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species

While Quaternary climatic changes are considered by some to have been a major factor promoting speciation within the neotropics, others suggest that much of the neotropical species diversity originated before the Pleistocene. Using mitochondrial and nuclear sequence data, we evaluate the relative importance of Pleistocene and pre-Pleistocene events within the evolutionary history of the Mexican iguana Ctenosaura pectinata , and related species. Results support the existence of cryptic lineages with strong mitochondrial divergence (> 4%) among them. Some of these lineages form zones of secondary contact, with one of them hybridizing with C. hemilopha . Evolutionary network analyses reveal the oldest populations of C. pectinata to be those of the northern and southern Mexican coastal regions. Inland and mid-latitudinal coastal populations are younger in age as a consequence of a history of local extinction within these regions followed by re-colonization. Estimated divergence times suggest that C. pectinata originated during the Pliocene, whereas geographically distinct mitochondrial DNA lineages first started to diverge during the Pliocene, with subsequent divergence continuing through the Pleistocene. Our results highlight the influence of both Pliocene and Pleistocene events in shaping the geographical distribution of genetic variation within neotropical lowland organisms. Areas of high genetic diversity in southern Mexico were detected, this finding plus the high levels of genetic diversity within C. pectinata , have implications for the conservation of this threatened species.     

Postglacial colonization history reflects in the genetic structure of natural populations of Festuca rubra in Europe

We conducted a large‐scale population genetic survey of genetic diversity of the host grass Festuca rubra s.l., which fitness can be highly dependent on its symbiotic fungus Epichloë festucae, to evaluate genetic variation and population structure across the European range. The 27 studied populations have previously been found to differ in frequencies of occurrence of the symbiotic fungus E. festucae and ploidy levels. As predicted, we found decreased genetic diversity in previously glaciated areas in comparison with nonglaciated regions and discovered three major maternal genetic groups: southern, northeastern, and northwestern Europe. Interestingly, host populations from Greenland were genetically similar to those from the Faroe Islands and Iceland, suggesting gene flow also between those areas. The level of variation among populations within regions is evidently highly dependent on the postglacial colonization history, in particular on the number of independent long‐distance seed colonization events. Yet, also anthropogenic effects may have affected the population structure in F. rubra. We did not observe higher fungal infection rates in grass populations with lower levels of genetic variability. In fact, the fungal infection rates of E. festucae in relation to genetic variability of the host populations varied widely among geographical areas, which indicate differences in population histories due to colonization events and possible costs of systemic fungi in harsh environmental conditions. We found that the plants of different ploidy levels are genetically closely related within geographic areas indicating independent formation of polyploids in different maternal lineages.     

Large sequence divergence of mitochondrial DNA genotypes of the control region within populations of the african antelope, kob (Kobus kob)

The genetic differentiation among kob populations (Kobus kob) representing two recognized subspecies was examined using mitochondrial control region sequences. Two distinct lineages (estimated sequence divergence of 9.8%) exhibited different geographical distributions and do not coincide with previously recognized ranges of subspecies. The presence of the two lineages was further supported with sequences of mitochondrial cytochrome b gene. One lineage was predominant in the west and southern ranges of the populations studied and the other was commonly found in a more northern distribution (Murchison populations) in Uganda. Murchison and the geographically intermediate Toro populations (Uganda) represented the area of overlap. The existence of the two lineages in the area of overlap is hypothesized to have resulted from a range expansion and secondary contact of the two lineages of kob that evolved in allopatry. The existence of the kob during the Pleistocene offers a plausible explanation for the observed biogeographic pattern. Our mitochondrial data reveal two examples of discordance between a gene tree and presumed species tree as: (i) the two lineages co-occur in the kob subspecies, Kobus kob thomasi (Uganda kob); and (ii) the puku, which was included in the analysis because of its controversial taxonomic status (currently recognized as a distinct species from the kob), is paraphyletic with respect to the kob. Significant degrees of heterogeneity were detected between populations. Relatively high genetic variation was observed in the populations, however, the inclusion of distinct lineages influences the population structure and nucleotide diversity of the kob populations.     

Genetic structure and host–parasite co‐divergence: evidence for trait‐specific local adaptation

Host–parasite co‐evolution can lead to genetic differentiation among isolated host–parasite populations and local adaptation between parasites and their hosts. However, tests of local adaptation rarely consider multiple fitness‐related traits although focus on a single component of fitness can be misleading. Here, we concomitantly examined genetic structure and co‐divergence patterns of the trematode Coitocaecum parvum and its crustacean host Paracalliope fluviatilis among isolated populations using the mitochondrial cytochrome oxidase I gene (COI). We then performed experimental cross‐infections between two genetically divergent host–parasite populations. Both hosts and parasites displayed genetic differentiation among populations, although genetic structure was less pronounced in the parasite. Data also supported a co‐divergence scenario between C. parvum and P. fluviatilis potentially related to local co‐adaptation. Results from cross‐infections indicated that some parasite lineages seemed to be locally adapted to their sympatric (home) hosts in which they achieved higher infection and survival rates than in allopatric (away) amphipods. However, local, intrinsic host and parasite characteristics (host behavioural or immunological resistance to infections, parasite infectivity or growth rate) also influenced patterns of host–parasite interactions. For example, overall host vulnerability to C. parvum varied between populations, regardless of parasite origin (local vs. foreign), potentially swamping apparent local co‐adaptation effects. Furthermore, local adaptation effects seemed trait specific; different components of parasite fitness (infection and survival rates, growth) responded differently to cross‐infections. Overall, data show that genetic differentiation is not inevitably coupled with local adaptation, and that the latter must be interpreted with caution in a multi‐trait context.     

Phylogeography of the parasitic fly Batrachomyia in the Wet Tropics of north-east Australia, and susceptibility of host frog lineages in a mosaic contact zone

Parasites could differentially impact intraspecific host lineages due to genetic, phenotypic, ecological, or behavioural differences between the lineages, or the development of reproductive isolation between them. Batrachomyia (Diptera: Chloropidae) are flies that exclusively parasitize Australian frogs, and in the Wet Tropics rainforest of north‐east Australia larvae are largely restricted to the green‐eyed tree frog Litoria genimaculata (Anura: Hylidae). This frog species consists of two highly divergent genetic lineages that overlap in two nearby, but independent, contact zones. At one contact zone there has been extensive phenotypic divergence and speciation between the lineages whereas, at the other contact relatively lower levels of phenotypic divergence and reproductive isolation suggest that speciation has not occurred. In the present study we tested: (1) whether the deep phylogeographic divergence between northern and southern host populations is mirrored by congruent genetic structuring in the parasite populations and (2) whether the host lineages are differentially impacted by parasitism. We found that the two divergent frog lineages are parasitized by a single lineage of Batrachomyia, which exhibits strikingly little phylogeographic structuring. We found a significant difference in Batrachomyia prevalence between the host lineages at mixed lineage sites in both contact zones, with the magnitude and direction of this effect being consistent in both. The pattern did not differ between the two contacts even though recent phenotypic divergence and speciation has occurred between the lineages at one contact but not the other. Taken together, this suggests a fundamental difference in susceptibility between the genetically divergent host lineages. Using weight relative to body length as a measure of body condition, we found no differential impact of parasitism on the body condition of each host lineage, and no evidence that parasitism impacts the body condition of the host in general. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 92 , 593–603.     

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs

Subdivided Pleistocene glacial refugia, best known as “refugia within refugia”, provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and “ghost” mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD‐sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply‐diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation‐by‐distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.     

Global genetic diversity,lineage distribution,and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well‐known example of a worldwide‐distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.