Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging, but endangered, freshwater mussel

Historical fragmentation and stepping-stone gene flow led to population genetic differentiation in a coastal seabird

Understanding the forces that shape population genetic structure is fundamental both for understanding evolutionary trajectories and for conservation. Many factors can influence the geographic distribution of genetic variation, and the extent to which local populations differ can be especially difficult to predict in highly mobile organisms. For example, many species of seabirds are essentially panmictic, but some show strong structure. Pigeon Guillemots (Cepphus columba; Charadriiformes: Alcidae) breed in small colonies scattered along the North Pacific coastline and feed in shallow nearshore waters year-round. Given their distribution, gene flow is potentially lower and population genetic structure is stronger than in most other high-latitude Northern Hemisphere seabirds. We screened variation in the mitochondrial control region, four microsatellite loci, and two nuclear introns in 202 Pigeon Guillemots representing three of five subspecies. Mitochondrial sequences and nuclear loci both showed significant population differences, although structure was weaker for the nuclear loci. Genetic differentiation was correlated with geographic distance between sampling locations for both the mitochondrial and nuclear loci. Mitochondrial gene trees and demographic modeling both provided strong evidence for two refugial populations during the Pleistocene glaciations: one in the Aleutian Islands and one farther east and south. We conclude that historical fragmentation combined with a stepping-stone model of gene flow led to the relatively strong population differentiation in Pigeon Guillemots compared to other high-latitude Northern Hemisphere seabird species. Our study adds to growing evidence that Pleistocene glaciation events affected population genetic structure not only in terrestrial species but also in coastal marine animals.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Mitochondrial sequence data reveal population structure within Pustulosa pustulosa

Unionid mussels are among the most imperiled group of organisms in North America, and Pustulosa pustulosa is a freshwater species with a relatively wide latitudinal distribution that extends from southern Ontario, Canada, to Texas, USA. Considerable morphological and geographic variation in the genus Pustulosa (formerly Cyclonaias) has led to uncertainty over species boundaries, and recent studies have suggested revisions to species-level classifications by synonymizing C. aurea, C. houstonensis, C. mortoni, and C. refulgens with C. pustulosa (currently P. pustulosa). Owing to its wide range and shallow phylogenetic differentiation, we analyzed individuals of P. pustulosa using mitochondrial DNA sequence data under a population genetics framework. We included 496 individuals, which were comprised of 166 samples collected during this study and 330 additional sequences retrieved from GenBank. Pairwise ΦST measures based on ND1 data suggested there may be up to five major geographic groups present within P. pustulosa. Genetic differentiation between regions within Texas was higher compared to populations from the Mississippi and Great Lakes populations, which may reflect differences in historical connectivity. Mitochondrial sequence data also revealed varying demographic histories for each major group suggesting each geographic region has also experienced differential population dynamics in the past. Future surveys should consider exploring variation within species after phylogeographic delimitation has been performed. In this study, we begin to address this need for freshwater mussels via the P. pustulosa system.

Cryptic Lineage and Genetic Structure of Acanthopagrus pacificus Populations in a Natural World Heritage Site Revealed by Population Genetic Analysis

Recent studies have revealed extensive genetic differentiation among some populations of marine taxa that were previously believed to be essentially homogeneous because larvae are widely dispersed in ocean currents. Acanthopagrus pacificus is a commercially and ecologically important teleost fish that is endemic to shallow coastal waters and estuaries of some tropical and sub-tropical areas in the West Pacific Ocean. Here, we examined genetic structure and the inferred demographic history of A. pacificus populations from mtDNA control region sequence data. A 677–678 base-pair fragment was sequenced from 159 individuals sampled at three localities across the West Pacific Ocean. Haplotype diversity was high, ranging from 0.915 to 0.989, while nucleotide diversity was low to medium, ranging from 0.8% to 2.60%. Analysis of molecular variance (AMOVA) showed significant genetic subdivision (FST = 0.155, p < 0.05) among sampled populations while pairwise FST estimates also revealed strong genetic differentiation among populations indicating that gene flow was restricted. Two distinct cryptic lineages were identified that were estimated to have diverged during the Pleistocene. In summary, contemporary factors including regional oceanic currents and self-recruitment are considered to have played significant roles in producing the population structure in this fish. In particular, the genetic information generated in the current study will allow appropriate fisheries management and conservation strategies to be developed for this important local fish in the waters around Iriomotejima Island, a World Heritage site.

Population structure and gene flow in the Sheepnose mussel ( Plethobasus cyphyus ) and their implications for conservation

North American freshwater mussel species have experienced substantial range fragmentation and population reductions. These impacts have the potential to reduce genetic connectivity among populations and increase the risk of losing genetic diversity. Thirteen microsatellite loci and an 883 bp fragment of the mitochondrial ND1 gene were used to assess genetic diversity, population structure, contemporary migration rates, and population size changes across the range of the Sheepnose mussel (Plethobasus cyphyus). Population structure analyses reveal five populations, three in the Upper Mississippi River Basin and two in the Ohio River Basin. Sampling locations exhibit a high degree of genetic diversity and contemporary migration estimates indicate that migration within river basins is occurring, although at low rates, but there is no migration is occurring between the Ohio and Mississippi river basins. No evidence of bottlenecks was detected, and almost all locations exhibited the signature of population expansion. Our results indicate that although anthropogenic activity has altered the landscape across the range of the Sheepnose, these activities have yet to be reflected in losses of genetic diversity. Efforts to conserve Sheepnose populations should focus on maintaining existing habitats and fostering genetic connectivity between extant demes to conserve remaining genetic diversity for future viable populations.

RAD-tag and mitochondrial DNA sequencing reveal the genetic structure of a widespread and regionally imperiled freshwater mussel, Obovaria olivaria (Bivalvia: Unionidae)

Obovaria olivaria is a species of freshwater mussel native to the Mississippi River and Laurentian Great Lakes-St. Lawrence River drainages of North America. This mussel has experienced population declines across large parts of its distribution and is imperiled in many jurisdictions. Obovaria olivaria uses the similarly imperiled Acipenser fulvescens (Lake Sturgeon) as a host for its glochidia. We employed mitochondrial DNA sequencing and restriction site-associated DNA sequencing (RAD-seq) to assess patterns of genetic diversity and population structure of O. olivaria from 19 collection locations including the St. Lawrence River drainage, the Great Lakes drainage, the Upper Mississippi River drainage, the Ohioan River drainage, and the Mississippi Embayment. Heterozygosity was highest in Upper Mississippi and Great Lakes populations, followed by a reduction in diversity and relative effective population size in the St. Lawrence populations. Pairwise F ST ranged from 0.00 to 0.20, and analyses of genetic structure revealed two major ancestral populations, one including all St. Lawrence River/Ottawa River sites and the other including remaining sites; however, significant admixture and isolation by river distance across the range were evident. The genetic diversity and structure of O. olivaria is consistent with the existing literature on Acipenser fulvescens and suggests that, although northern and southern O. olivaria populations are genetically distinct, genetic structure in O. olivaria is largely clinal rather than discrete across its range. Conservation and restoration efforts of O. olivaria should prioritize the maintenance and restoration of locations where O. olivaria remain, especially in northern rivers, and to ensure connectivity that will facilitate dispersal of Acipenser fulvescens and movement of encysted glochidia.

Genetic and morphological characterization of the freshwater mussel clubshell species complex (Pleurobema clava and Pleurobema oviforme) to inform conservation planning

The shell morphologies of the freshwater mussel species Pleurobema clava (federally endangered) and Pleurobema oviforme (species of concern) are similar, causing considerable taxonomic confusion between the two species over the last 100 years. While P. clava was historically widespread throughout the Ohio River basin and tributaries to the lower Laurentian Great Lakes, P. oviforme was confined to the Tennessee and the upper Cumberland River basins. We used two mitochondrial DNA (mtDNA) genes, 13 novel nuclear DNA microsatellite markers, and shell morphometrics to help resolve this taxonomic confusion. Evidence for a single species was apparent in phylogenetic analyses of each mtDNA gene, revealing monophyletic relationships with minimal differentiation and shared haplotypes. Analyses of microsatellites showed significant genetic structuring, with four main genetic clusters detected, respectively, in the upper Ohio River basin, the lower Ohio River and Great Lakes, and upper Tennessee River basin, and a fourth genetic cluster, which included geographically intermediate populations in the Ohio and Tennessee river basins. While principal components analysis (PCA) of morphometric variables (i.e., length, height, width, and weight) showed significant differences in shell shape, only 3% of the variance in shell shape was explained by nominal species. Using Linear Discriminant and Random Forest (RF) analyses, correct classification rates for the two species' shell forms were 65.5% and 83.2%, respectively. Random Forest classification rates for some populations were higher; for example, for North Fork Holston (HOLS), it was >90%. While nuclear DNA and shell morphology indicate that the HOLS population is strongly differentiated, perhaps indicative of cryptic biodiversity, we consider the presence of a single widespread species the most likely biological scenario for many of the investigated populations based on our mtDNA dataset. However, additional sampling of P. oviforme populations at nuclear loci is needed to corroborate this finding.

Song type variations of Louisiana Waterthrush (Parkesia motacilla) and their geographic distributions

Louisiana Waterthrush (Parkesia motacilla) is a familiar singer in the Western Hemisphere family Parulidae, yet apparent geographic variations in its song and potentially related causal mechanisms have not received detailed examination in previously published studies. Here, we analyzed song pattern variations of 651 Louisiana Waterthrush singers in audio spectrogram recordings obtained from our field work and publicly accessible bioacoustics archives. Visual and auditory assessment of the introductory note sequence of each song identified three distinct song types (A, B, and C) and most of the songs were assigned to one of these types. Linear Discriminant Analysis and Random Forest methods were used to verify the assignments and showed strong agreement for Type A with slightly less agreement on Types B and C. User error rates (proportion of the Linear Discriminant Analysis classifications that were incorrect) were low for Types A and B, and somewhat higher for Type C, while producer error rates (proportion of the song type for which the Linear Discriminant Analysis was incorrect) were somewhat higher for Types A and C than the minimal levels achieved for Type B. Our findings confirmed that most between-individual variation was in the number of notes and note sequence duration while most within-individual variation resulted from the percent of downstrokes. The location of each singer was plotted on a map of the breeding range and results suggested the song types have large-scale discrete geographic distributions that co-occur in some regions but not range-wide. Evaluation of the distributions provided tentative support for a hypothesis that two of the song types may independently exhibit congruence with the geographic extent of Pleistocene glacial boundaries and the third song type may be distinguished by a lack of congruence, but further investigation is needed to elucidate whether the song variations represent subpopulations with three separate evolutionary histories.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Genome-wide SNPs redefines species boundaries and conservation units in the freshwater mussel genus Cyprogenia of North America

Detailed information on species delineation and population genetic structure is a prerequisite for designing effective restoration and conservation strategies for imperiled organisms. Phylogenomic and population genomic analyses based on genome-wide double digest restriction-site associated DNA sequencing (ddRAD-Seq) data has identified three allopatric lineages in the North American freshwater mussel genus Cyprogenia. Cyprogenia stegaria is restricted to the Eastern Highlands and displays little genetic structuring within this region. However, two allopatric lineages of C. aberti in the Ozark and Ouachita highlands exhibit substantial levels (mean uncorrected F_ST = 0.368) of genetic differentiation and each warrants recognition as a distinct evolutionary lineage. Lineages of Cyprogenia in the Ouachita and Ozark highlands are further subdivided reflecting structuring at the level of river systems. Species tree inference and species delimitation in a Bayesian framework using single nucleotide polymorphisms (SNP) data supported results from phylogenetic analyses, and supports three species of Cyprogenia over the currently recognized two species. A comparison of SNPs generated from both destructively and non-destructively collected samples revealed no significant difference in the SNP error rate, quality and amount of ddRAD sequence reads, indicating that nondestructive or trace samples can be effectively utilized to generate SNP data for organisms for which destructive sampling is not permitted.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

Riverscape genetic variation, migration patterns, and morphological variation of the threatened Round Rocksnail, Leptoxis ampla

Within riverine systems, headwater populations are hypothesized to harbour higher amounts of genetic distinctiveness than populations in the main stem of a river and display increased genetic diversity in large, downstream habitats. However, these hypotheses were mostly developed with insects and fish, and they have not been tested on many invertebrate lineages. Pleuroceridae gastropods are of particular ecological importance to rivers of eastern North America, sometimes comprising over 90% of macroinvertebrate biomass. Yet, virtually nothing is known of pleurocerid landscape genetics, including whether genetic diversity follows predictions made by hypotheses developed on more mobile species. Moreover, the commonly repeated hypothesis that intraspecific morphological variation in gastropods results from ecophenotypic plasticity has not been well tested on pleurocerids. Using 2bRAD-seq to discover single nucleotide polymorphisms, we show that the threatened, Cahaba River endemic pleurocerid, Leptoxis ampla, has limited gene flow among populations and that migration is downstream-biased, conflicting with previous hypotheses. Both tributary and main stem populations harbour unique genomic profiles, and genetic diversity was highest in downstream populations. Furthermore, L. ampla shell morphology was more correlated with genetic differences among individuals and populations than habitat characteristics. We anticipate similar genetic and demographic patterns to be seen in other pleurocerids, and hypotheses about gene flow and population demographics that were based on more mobile taxa often, but not always, apply to freshwater gastropods. From a conservation standpoint, genetic structure of L. ampla populations suggests distinctive genetic diversity is lost with localized extirpation, a phenomenon common across the range of Pleuroceridae.

Phylogeography of Bellamya (Mollusca: Gastropoda: Viviparidae) snails on different continents: contrasting patterns of diversification in China and East Africa

Background

Species diversity is determined by both local environmental conditions that control differentiation and extinction and the outcome of large-scale processes that affect migration. The latter primarily comprises climatic change and dynamic landscape alteration. In the past few million years, both Southeast Asia and Eastern Africa experienced drastic climatic and geological oscillations: in Southeast Asia, especially in China, the Tibetan Plateau significantly rose up, and the flow of the Yangtze River was reversed. In East Africa, lakes and rivers experienced frequent range expansions and regressions due to the African mega-droughts. To test how such climatic and geological histories of both regions relate to their respective regional species and genetic diversity, a large scale comparative phylogeographic study is essential. Bellamya, a species rich freshwater snail genus that is widely distributed across China and East Africa, represents a suitable model system to address this question. We sequenced mitochondrial and nuclear DNA for members of the genus from China and used published sequences from Africa and some other locations in Asia to investigate their phylogeny and distribution of genetic diversity.

Results

Our phylogenetic analysis revealed two monophyletic groups, one in China and one in East Africa. Within the Chinese group, Bellamya species show little genetic differentiation. In contrast, we observe fairly deep divergence among the East African lakes with almost every lake possessing its unique clade. Our results show that strong divergence does not necessarily depend on intrinsic characteristics of a species, but rather is related to the landscape dynamics of a region.

Conclusion

Our phylogenetic results suggest that the Bellamya in China and East Africa are independent phylogenetic clades with different evolutionary trajectories. The different climate and geological histories likely contributed to the diverging evolutionary patterns. Repeated range expansions and regressions of lakes likely contributed to the great divergence of Bellamya in East Africa, while reversal of the river courses and intermingling of different lineages had an opposite effect on Bellamya diversification in China.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1397-0) contains supplementary material, which is available to authorized users.

Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa

The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.

Predicting the effects of climate change on population connectivity and genetic diversity of an imperiled freshwater mussel, Cumberlandia monodonta (Bivalvia: Margaritiferidae), in riverine systems

In the face of global climate change, organisms may respond to temperature increases by shifting their ranges poleward or to higher altitudes. However, the direction of range shifts in riverine systems is less clear. Because rivers are dendritic networks, there is only one dispersal route from any given location to another. Thus, range shifts are only possible if branches are connected by suitable habitat, and stream-dwelling organisms can disperse through these branches. We used Cumberlandia monodonta (Bivalvia: Unionoida: Margaritiferidae) as a model species to investigate the effects of climate change on population connectivity because a majority of contemporary populations are panmictic. We combined ecological niche models (ENMs) with population genetic simulations to investigate the effects of climate change on population connectivity and genetic diversity of C. monodonta. The ENMs were constructed using bioclimatic and landscape data to project shifts in suitable habitat under future climate scenarios. We then used forward-time simulations to project potential changes in genetic diversity and population connectivity based on these range shifts. ENM results under current conditions indicated long stretches of highly suitable habitat in rivers where C. monodonta persists; populations in the upper Mississippi River remain connected by suitable habitat that does not impede gene flow. Future climate scenarios projected northward and headwater-ward range contraction and drastic declines in habitat suitability for most extant populations throughout the Mississippi River Basin. Simulations indicated that climate change would greatly reduce genetic diversity and connectivity across populations. Results suggest that a single, large population of C. monodonta will become further fragmented into smaller populations, each of which will be isolated and begin to differentiate genetically. Because C. monodonta is a widely distributed species and purely aquatic, our results suggest that persistence and connectivity of stream-dwelling organisms will be significantly altered in response to future climate change.

The eastern migratory caribou: the role of genetic introgression in ecotype evolution

Understanding the evolutionary history of contemporary animal groups is essential for conservation and management of endangered species like caribou (Rangifer tarandus). In central Canada, the ranges of two caribou subspecies (barren-ground/woodland caribou) and two woodland caribou ecotypes (boreal/eastern migratory) overlap. Our objectives were to reconstruct the evolutionary history of the eastern migratory ecotype and to assess the potential role of introgression in ecotype evolution. STRUCTURE analyses identified five higher order groups (i.e. three boreal caribou populations, eastern migratory ecotype and barren-ground). The evolutionary history of the eastern migratory ecotype was best explained by an early genetic introgression from barren-ground into a woodland caribou lineage during the Late Pleistocene and subsequent divergence of the eastern migratory ecotype during the Holocene. These results are consistent with the retreat of the Laurentide ice sheet and the colonization of the Hudson Bay coastal areas subsequent to the establishment of forest tundra vegetation approximately 7000 years ago. This historical reconstruction of the eastern migratory ecotype further supports its current classification as a conservation unit, specifically a Designatable Unit, under Canada's Species at Risk Act. These findings have implications for other sub-specific contact zones for caribou and other North American species in conservation unit delineation.

Northern Bobwhite (Colinus virginianus) Mitochondrial Population Genomics Reveals Structure, Divergence, and Evidence for Heteroplasmy

Herein, we evaluated the concordance of population inferences and conclusions resulting from the analysis of short mitochondrial fragments (i.e., partial or complete D-Loop nucleotide sequences) versus complete mitogenome sequences for 53 bobwhites representing six ecoregions across TX and OK (USA). Median joining (MJ) haplotype networks demonstrated that analyses performed using small mitochondrial fragments were insufficient for estimating the true (i.e., complete) mitogenome haplotype structure, corresponding levels of divergence, and maternal population history of our samples. Notably, discordant demographic inferences were observed when mismatch distributions of partial (i.e., partial D-Loop) versus complete mitogenome sequences were compared, with the reduction in mitochondrial genomic information content observed to encourage spurious inferences in our samples. A probabilistic approach to variant prediction for the complete bobwhite mitogenomes revealed 344 segregating sites corresponding to 347 total mutations, including 49 putative nonsynonymous single nucleotide variants (SNVs) distributed across 12 protein coding genes. Evidence of gross heteroplasmy was observed for 13 bobwhites, with 10 of the 13 heteroplasmies involving one moderate to high frequency SNV. Haplotype network and phylogenetic analyses for the complete bobwhite mitogenome sequences revealed two divergent maternal lineages (d_XY = 0.00731; F_ST = 0.849; P < 0.05), thereby supporting the potential for two putative subspecies. However, the diverged lineage (n = 103 variants) almost exclusively involved bobwhites geographically classified as Colinus virginianus texanus, which is discordant with the expectations of previous geographic subspecies designations. Tests of adaptive evolution for functional divergence (MKT), frequency distribution tests (D, F_S) and phylogenetic analyses (RAxML) provide no evidence for positive selection or hybridization with the sympatric scaled quail (Callipepla squamata) as being explanatory factors for the two bobwhite maternal lineages observed. Instead, our analyses support the supposition that two diverged maternal lineages have survived from pre-expansion to post-expansion population(s), with the segregation of some slightly deleterious nonsynonymous mutations.

Taxonomy and distribution of freshwater pearl mussels (Unionoida: Margaritiferidae) of the Russian Far East

The freshwater pearl mussel family Margaritiferidae includes 13 extant species, which are all listed by IUCN as endangered or vulnerable taxa. In this study, an extensive spatial sampling of Margaritifera spp. across the Russian Far East (Amur Basin, Kamchatka Peninsula, Kurile Archipelago and Sakhalin Island) was conducted for a revision of their taxonomy and distribution ranges. Based on their DNA sequences, shell and soft tissue morphology, three valid species were identified: Margaritifera dahurica (Middendorff, 1850), M. laevis (Haas, 1910) and M. middendorffi (Rosén, 1926). M. dahurica ranges across the Amur basin and some of the nearest river systems. M. laevis is distributed in Japan, Sakhalin Island and the Kurile Archipelago. M. middendorffi was previously considered an endemic species of the Kamchatka. However, it is widespread in the rivers of Kamchatka, Sakhalin Island, the Kurile Islands (across the Bussol Strait, which is the most significant biogeographical boundary within the archipelago), and, likely, in Japan. The Japanese species M. togakushiensis Kondo & Kobayashi, 2005 seems to be conspecific with M. middendorffi because of similar morphological patterns, small shell size (<100 mm long) and overlapped ranges, but it is in need of a separate revision. Phylogenetic analysis reveals that two NW Pacific margaritiferid species, M. laevis and M. middendorffi, formed a monophyletic 18S rDNA clade together with the North American species M. marrianae and M. falcata. The patterns that were found in these Margaritifera spp. are similar to those of freshwater fishes, indicating multiple colonizations of Eastern Asia by different mitochondrial lineages, including an ancient Beringian exchange between freshwater faunas across the Pacific.