Mitochondrial DNA Sequence Variation in a Sea Star (Leptasteriasspp.) Species Complex

A 551-bp region of a PCR product containing the putative mitochondrial control region and flanking sequences was analyzed for sequence variation among 19 sea stars representing 10 previously described PCR–RFLP haplotypes within a cryptic species complex (Leptasteriasspp.). Most (97%) of the sequence variation was interhaplotypic rather than intrahaplotypic, which greatly reduced the utility of sequence polymorphisms in this mtDNA region as markers of intrahaplotypic population structure and gene flow. The estimated number of transition and transversion substitutions per nucleotide site, corrected for multiple hits, was 0.0364 and 0.0158, respectively. Most of the sequence variation occurred in the first half of the putative control region. Phylogenetic analysis (both maximum parsimony and maximum likelihood) revealed three well-supported clades, but the position of two PCR-RFLP haplotypes was not completely resolved. Low intraspecific mtDNA sequence divergence over large geographic distances may be a general pattern for echinoderm species.     

Dual African Origins of Global Aedes aegypti s.l. Populations Revealed by Mitochondrial DNA

Background

Aedes aegypti is the primary global vector to humans of yellow fever and dengue flaviviruses. Over the past 50 years, many population genetic studies have documented large genetic differences among global populations of this species. These studies initially used morphological polymorphisms, followed later by allozymes, and most recently various molecular genetic markers including microsatellites and mitochondrial markers. In particular, since 2000, fourteen publications and four unpublished datasets have used sequence data from the NADH dehydrogenase subunit 4 mitochondrial gene to compare Ae. aegypti collections and collectively 95 unique mtDNA haplotypes have been found. Phylogenetic analyses in these many studies consistently resolved two clades but no comprehensive study of mtDNA haplotypes have been made in Africa, the continent in which the species originated.

Methods and Findings

ND4 haplotypes were sequenced in 426 Ae. aegypti s.l. from Senegal, West Africa and Kenya, East Africa. In Senegal 15 and in Kenya 7 new haplotypes were discovered. When added to the 95 published haplotypes and including 6 African Aedes species as outgroups, phylogenetic analyses showed that all but one Senegal haplotype occurred in a basal clade while most East African haplotypes occurred in a second clade arising from the basal clade. Globally distributed haplotypes occurred in both clades demonstrating that populations outside Africa consist of mixtures of mosquitoes from both clades.

Conclusions

Populations of Ae. aegypti outside Africa consist of mosquitoes arising from one of two ancestral clades. One clade is basal and primarily associated with West Africa while the second arises from the first and contains primarily mosquitoes from East Africa     

Polyphyly and species delimitation of Picea brachytyla (Pinaceae) based on population genetic data

Accurate species delimitation is fundamental to biodiversity conservation. The endangered spruce Picea brachytyla (Franch.) E. Pritz. was suggested to be polyphyletic based on a limited number of samples in previous studies. To evaluate polyphyly of P. brachytyla, we sampled 139 individuals from 16 populations across most of its distributional range, plus representatives of two related species, Picea likiangensis (Franch.) E. Pritz. and Picea wilsonii Mast. We sequenced 13 nuclear loci and three chloroplast and two mitochondrial loci for the following species delimitation. Phylogenetic analyses of nuclear loci grouped all individuals of P. brachytyla from Sichuan and Chongqing into one distinct lineage and those from Yunnan and Tibet (southern distribution) nested within the P. likiangensis species complex. Structure analyses confirmed this result. Networks of chloroplast DNA haplotypes similarly showed that P. brachytyla from the southern distribution nested within the P. likiangensis species complex, whereas haplotypes for the northern distribution comprised a separate and well-supported lineage. These results suggest that P. brachytyla from the southern distribution is a part of the P. likiangensis species complex and should be removed from P. brachytyla. Our study highlights the utility of population genetic evidence in delimitating endangered species and understanding the conservation status of such species.     

Phylogenetic Relationships and Species Delimitation in Pinus Section Trifoliae Inferrred from Plastid DNA

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.     

Distribution of nuclear mitochondrial pseudogenes in three pollinator fig wasps associated with Ficus pumila

Nuclear mitochondrial pseudogenes (NUMTs) are nuclear sequences transferred from mitochondrial genomes. Although widespread, their distribution patterns among populations or closely related species are rarely documented. We amplified and sequenced the mitochondrial cytochrome b (Cytb) gene to check for NUMTs in three fig wasp species that pollinate Ficus pumila (Wiebesia sp. 1, 2 and 3) in Southeastern China using direct and cloned sequencing. Unambiguous sequences (332) of 487 bp in length belonging to 33 haplotypes were found by direct sequencing. Their distribution was highly concordant with those of cytochrome c oxidase subunit I (COI). Obvious signs of co-amplification of NUMTs were indicated by their uneven distribution. NUMTs were observed in all individuals of 12 populations of Wiebesia sp. 3, and 13 individuals of three northern populations of Wiebesia sp. 1. Sequencing clones of potential co-amplification products confirmed that they were NUMTs. These NUMTs either clustered as NUMT clades basal to mtDNA Cytb clades (basal NUMTs), or together with Cytb haplotypes. Basal NUMTs had either stop codons or frame-shifting mutations resulting from deletion of a 106 bp fragment. In addition, no third codon or synonymous substitutions were detected within each NUMT clade. The phylogenetic tree indicated that basal NUMTs had been inserted into nuclei before divergence of the three species. No significant pairwise differences were detected in their ratios of third codon substitutions, suggesting that these NUMTs originated from one transfer event, with duplication in the nuclear genome resulting in the coexistence of the 381 bp copy. No significant substitution differences were detected between Cytb haplotypes and NUMTs that clustered with Cytb haplotypes. However, these NUMTs coexisted with Cytb haplotypes in multiple populations, suggesting that these NUMT haplotypes were recently inserted into the nuclear genome. Both basal and recently inserted NUMTs were rare events, and were absent in most populations of Wiebesia sp. 1 and 2. Further studies are needed to distinguish between mechanisms potentially generating this rarity, such as purifying selection, genetic drift or amplification failure.     

Assessing diversity of King Crab Lithodes spp. in the south‐eastern pacific using phylogeny and molecular species delimitation methods

The purpose of this study was to test the hypothesis that the genetic diversity of commercially significant species of King Crabs (Lithodes spp.) along the south‐eastern Pacific (SEP) comprises different independent evolutionary units (IEUs) with spatially isolated distribution. Nine localities from inner and open waters along the SEP Chilean coast (39°S‐55°S) were sampled. We analyzed sequences from 173 individuals for the mitochondrial gene Cytochrome oxidase I (COX‐I), 151 individuals for the Internal Transcribed Spacer 1 (ITS) and 135 for the structural ribosomal RNA (28S). Genetic delimitation was performed through three analytical methods: ABGD, GMYC, and its Bayesian implementation, bGMYC. Bayesian phylogenetic analyses and haplotype networks were also performed. Divergence time between clades was assessed for the COX‐I marker and estimated from known evolutionary rates for this marker in other crustacean species and fossil calibration from other Anomuran species. Delimitation analyses, phylogenetic analyses, and mitochondrial haplotype networks suggested the presence of two deeply divergent mitochondrial lineages of Lithodes in the SEP, referred to as Clade1 and Clade 2. Nuclear markers showed low phylogenetic resolution and therefore were unsuitable for molecular species delimitation. Divergence time analysis of the mitochondrial lineages suggests a separation between Clades of approximately 2.3 Mya. The divergence time obtained suggested that Pliocene glaciations and deglaciations cycles could be involved in hybridization events between Lithodes IEUs at southern tip of South American coasts. The different frequencies of Lithodes haplotypes in inner and open water environments along SEP coasts could be explained by events such as the last glacial maximum or by differences in the adaptation of each clade to different environments. These findings support the necessity of evaluating the taxonomic status of Lithodes individuals found along SEP coasts under an integrative taxonomy approach or through markers with other evolution rates than those already used.     

Identifying two moth species (Lepidoptera: Ditrysia) from Saudi Arabia using mitochondrial 16S rRNA sequences

The mitochondrial genetic markers are considered useful tools for discrimination between more closely related lepidopteran taxa. Therefore, the present study aimed to investigate the role of mitochondrial (mt) 16 s rRNA gene in the determination of the taxonomic position for two moth species within Ditrysia clade. Maximum likelihood analysis has indicated a well-supported dendrogram based on the Tamura-Nei model for the recovered lepidopterans. The mt 16 s rRNA query sequences from 24 species within seven families were analyzed. This analysis and bootstrap confidence revealed two major clades representing Glossata suborder within Lepidoptera, with a close relationship of Noctuoidea + (Pyraloidea (Hesperioidea + Papilionoidea)). The subfamily Heliothinae forming a sister group with Risobinae (Noctinae + Hadeninae). In addition, there is a clear observation about the close relation between Phycitinae + Galleriinae within Pyraloidea and Cyrestinae + Limenitidinae within Papilionoidea. The present study supported that the Helicoverpa and Meroptera species are the first accounts of these genera inhabiting Saudi Arabia.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Integrative species delimitation reveals cryptic diversity in the southern Appalachian Antrodiaetus unicolor (Araneae: Antrodiaetidae) species complex

Although species delimitation can be highly contentious, the development of reliable methods to accurately ascertain species boundaries is an imperative step in cataloguing and describing Earth's quickly disappearing biodiversity. Spider species delimitation remains largely based on morphological characters; however, many mygalomorph spider populations are morphologically indistinguishable from each other yet have considerable molecular divergence. The focus of our study, the Antrodiaetus unicolor species complex containing two sympatric species, exhibits this pattern of relative morphological stasis with considerable genetic divergence across its distribution. A past study using two molecular markers, COI and 28S, revealed that A. unicolor is paraphyletic with respect to A. microunicolor. To better investigate species boundaries in the complex, we implement the cohesion species concept and use multiple lines of evidence for testing genetic exchangeability and ecological interchangeability. Our integrative approach includes extensively sampling homologous loci across the genome using a RADseq approach (3RAD), assessing population structure across their geographic range using multiple genetic clustering analyses that include structure , principal components analysis and a recently developed unsupervised machine learning approach (Variational Autoencoder). We evaluate ecological similarity by using large‐scale ecological data for niche‐based distribution modelling. Based on our analyses, we conclude that this complex has at least one additional species as well as confirm species delimitations based on previous less comprehensive approaches. Our study demonstrates the efficacy of genomic‐scale data for recognizing cryptic species, suggesting that species delimitation with one data type, whether one mitochondrial gene or morphology, may underestimate true species diversity in morphologically homogenous taxa with low vagility.     

Influence of management regime and population history on genetic diversity and population structure of brown hares (Lepus europaeus) in an Italian province

In many areas, the management of overexploited populations of brown hare (Lepus europaeus) is based on annual restocking. While in some cases exotic hares are introduced, in some others hares are captured locally within protected areas and subsequently released into hunting grounds. We evaluated the genetic effects of this management regime in an Italian province where the brown hare population has recovered in the last few decades, by sequencing the hypervariable domain 1 of the mitochondrial control region and by genotyping eight autosomal microsatellites in hares sampled in both hunting and non-hunting areas. Both nuclear (H _e?=?0.68 and H _o?=?0.65) and mitochondrial variability (h?=?0.853 and π?=?0.012) were in line with other European populations. When comparing our data with mitochondrial sequences retrieved from GenBank, out of the 21 detected haplotypes, 14 were private to our study area. While 4.6 % of the individuals were found to carry haplotypes attributable to past introductions, 41.5 % grouped within a well-supported lineage, previously identified with a presumed native Italian taxon, L. e. meridiei. Despite the detectable geographic partitioning of mitochondrial haplotypes across the province, no genetic structure resulted from microsatellites analysis, indicating that no reproductive barriers exist among hares carrying different mitochondrial lineages. In conclusion, the local management seems to have contributed to the recovery of the species and to a full admixture of nuclear genes in the province. However, neither the extensive translocations nor the possible introductions of exotic heads seem to have completely undermined the local mitochondrial lineages.     

Disintegrating over space and time: Paraphyly and species delimitation in the Wehrle's Salamander complex

Accurate species delimitation is critical for biodiversity studies. However, species complexes characterized by introgression, high levels of population structure and subtle phenotypic differentiation can be challenging to delimit. Here, we report on a molecular systematic investigation of the woodland salamanders Plethodon wehrlei and Plethodon punctatus, which traditionally have been placed in the Plethodon wehrlei species group. To quantify patterns of genetic variation, we collected genetic samples from throughout the range of both species, including 22 individuals from nine populations of P. punctatus, and 60 individuals from 26 populations of P. wehrlei. From these samples, we sequenced three mtDNA loci (5596 base pairs) and five nuclear loci (3377 base pairs). We inferred time‐calibrated gene trees and species trees using BEAST 2.4.6, and we delimited putative species using a Bayesian implementation of the general mixed Yule‐coalescent model (bGMYC) and STRUCTURE. Finally, we validated putative species using the multispecies coalescent as implemented in Bayesian Phylogenetics and Phylogeography (BPP). We found substantial phylogeographic diversity in P. wehrlei, including multiple geographically cohesive clades and an inferred mitochondrial common ancestor at 11.5 myr (95% HPD: 9.6–13.6 myr) that separated populations formerly assigned to P. dixi from all other populations. We also found that P. punctatus is deeply nested within P. wehrlei, rendering the latter paraphyletic. After discussing the challenges faced by modern species delimitation methods, we recommend retaining P. punctatus because it is ecologically and phenotypically distinct. We further recommend that P. dixi be recognized as a valid species.     

Species boundaries in the Himantura uarnak species complex (Myliobatiformes: Dasyatidae)

Samples of the ‘Himantura uarnak’ species complex (H. leoparda, H. uarnak, H. undulata under their current definitions), mostly from the Coral Triangle, were analyzed using nuclear markers and mitochondrial DNA sequences. Genotypes at five intron loci showed four reproductively isolated clusters of individuals. The COI sequences showed four major mitochondrial lineages, each diagnostic of a cluster as defined by nuclear markers. No mitochondrial introgression was detected. The average Kimura-2 parameter nucleotide distance separating clades was 0.061–0.120 (net: 0.055–0.114), while the distance separating individuals within a clade was 0.002–0.008. Additional, partial cytochrome-b gene sequences were used to link these samples with previously published sequences of reference specimens of the three nominal species. One of the clusters was identified as H. undulata and another one, as H. uarnak, while two cryptic species were uncovered within the recently-described H. leoparda, challenging the current morphology-based taxonomy of species within the H. uarnak species complex.     

Population structure and species boundary delimitation of cryptic Dioryctria moths: an integrative approach

Accurate delimitation of species boundaries is especially important in cryptic taxa where one or more character sources are uninformative or are in conflict. Rather than relying on a single marker to delimit species, integrative taxonomy uses multiple lines of evidence such as molecular, morphological, behavioural and geographic characters to test species limits. We examine the effectiveness of this approach by testing the delimitation of two cryptic Nearctic species of Dioryctria (Lepidoptera: Pyralidae) using three independent molecular markers [cytochrome c oxidase I (COI), second internal transcribed spacer unit (ITS2), and elongation factor 1alpha (EF1alpha)], forewing variation and larval host plant association. Although mitochondrial DNA (mtDNA) haplotypes do not form reciprocally monophyletic clades, restricted gene flow between COI haplotype groups, and concordance with ITS2 genotypes, forewing variation and host plant associations support delimitation of two Nearctic species: eastern Dioryctria reniculelloides and western Dioryctria pseudotsugella. Conversely, EF1alpha genotype variation was incongruent with the two previous markers. A case of discordance between COI and ITS2 was detected, suggesting either introgression due to hybridization or retained ancestral polymorphism due to incomplete coalescence. This study is consistent with other similar literature where molecular loci in closely related species progress from shared to fixed haplotypes/alleles, and from polyphyletic to reciprocally monophyletic relationships, although loci may vary in these characteristics despite maintenance of genomic integrity between distinct species. In particular, mtDNA in other studies generally showed a lower rate of fixation of differences than did X-linked or autosomal loci, reinforcing the need to use an integrative approach for delimiting species.     

Genetic subdivision and biogeography of the Danubian rheophilic barb Barbus petenyi inferred from phylogenetic analysis of mitochondrial DNA variation

The barb Barbus petenyi is a cyprinid widely distributed throughout the mountain regions in the Danube River basin. Phylogenetic analysis of the DNA sequence variation at the mitochondrial cytochrome b gene over much of this range yielded three deep-branching (5.9-9.4% average divergence), well-supported haplotype clades with mutually exclusive geographic distributions and divergence times estimated to be in the Tertiary. The clades did not form an altogether monophyletic group as the most divergent one coalesced more recently with haplotypes of phylogenetically close species than with the other B. petenyi haplotypes. This pattern was supported by bootstrap and log-likelihood Shimodaira-Hasegawa tests. The other two were sister clades, but their distinctiveness was supported by previous allozyme data. Hence, from a taxonomic point of view, the current recognition of B. petenyi is erroneous, as it does not represent a single evolutionary lineage, and we suggest that three species be recognized instead. Substantial phylogeographic differences were evident among the three putative species, the two more southerly ones displaying significant structure, which suggested that they each survived in several glacial refugia throughout the Pleistocene. The phylogeographic pattern of multiple populations of rheophilic barbs with a history of long-term persistence and separation within the Danube River basin is novel within fishes and provides a hypothesis against which phylogeographic patterns among other similarly distributed rheophilic species may be compared.     

Species boundaries in Malagasy snakes of the genus Madagascarophis (Serpentes: Colubridae sensu lato) assessed by nuclear and mitochondrial markers

We studied mitochondrial divergence in 27 individuals of colubrid snakes of the genus Madagascarophis Mertens from most of its distribution area in Madagascar. Combined analyses of 16S rRNA and cytochrome b sequences identified six major clades which only partly agreed with previously proposed classifications. Analysis of nuclear DNA sequences of the c-mos gene as well as of ISSR fingerprints revealed consistent differences only among three clades which we consider as distinct species: a widespread Madagascarophis colubrinus (Schlegel), with M. citrinus (Boettger) as a junior synonym, a southern M. meridionalis Domergue, and a presumably undescribed species from the extreme north of Madagascar. The species M. ocellatus Domergue was not available for our study. Within M. colubrinus there are two populations from the north-west, each showing two divergent haplotypes with pairwise divergences of up to 5.2% in the cytochrome b gene. Maximum divergence in this gene within M. colubrinus was 7.1%. These high values emphasise that caution needs to be applied before genetic distance values are used for species delimitation. Phylogeographically, most of the genetic variation in M. colubrinus is found in northern Madagascar, indicating that the species might have originated in this region. Later one haplotype clade colonised western and eastern Madagascar, with a putative secondary introgression into north-western populations.     

Two Different High Throughput Sequencing Approaches Identify Thousands of De Novo Genomic Markers for the Genetically Depleted Bornean Elephant

High throughput sequencing technologies are being applied to an increasing number of model species with a high-quality reference genome. The application and analyses of whole-genome sequence data in non-model species with no prior genomic information are currently under way. Recent sequencing technologies provide new opportunities for gathering genomic data in natural populations, laying the empirical foundation for future research in the field of conservation and population genomics. Here we present the case study of the Bornean elephant, which is the most endangered subspecies of Asian elephant and exhibits very low genetic diversity. We used two different sequencing platforms, the Roche 454 FLX (shotgun) and Illumina, GAIIx (Restriction site associated DNA, RAD) to evaluate the feasibility of the two methodologies for the discovery of de novo markers (single nucleotide polymorphism, SNPs and microsatellites) using low coverage data. Approximately, 6,683 (shotgun) and 14,724 (RAD) SNPs were detected within our elephant sequence dataset. Genotyping of a representative sample of 194 SNPs resulted in a SNP validation rate of ∼ 83 to 94% and 17% of the loci were polymorphic with a low diversity (H _o = 0.057). Different numbers of microsatellites were identified through shotgun (27,226) and RAD (868) techniques. Out of all di-, tri-, and tetra-microsatellite loci, 1,706 loci had sufficient flanking regions (shotgun) while only 7 were found with RAD. All microsatellites were monomorphic in the Bornean but polymorphic in another elephant subspecies. Despite using different sample sizes, and the well known differences in the two platforms used regarding sequence length and throughput, the two approaches showed high validation rate. The approaches used here for marker development in a threatened species demonstrate the utility of high throughput sequencing technologies as a starting point for the development of genomic tools in a non-model species and in particular for a species with low genetic diversity.     

Single-locus DNA barcoding and species delimitation of the sandfly subgenus Evandromyia (Aldamyia)

Sandfly specimens from the subgenus Evandromyia (Aldamyia) Galati, 2003 (Diptera: Psychodidae: Phlebotominae) were collected between 2012 and 2019 from nine localities in seven Brazilian states, morphologically-identified, and then DNA barcoded by sequencing the mitochondrial cytochrome c oxidase subunit I (coi) gene. Forty-four new barcode sequences generated from 10 morphospecies were combined with 49 previously published sequences from the same subgenus and analysed using sequence-similarity methods (best-match criteria) to assess their ability at specimen identification, while four different species delimitation methods (ABGD, GMYC, PTP and TCS) were used to infer molecular operational taxonomic units (MOTUs). Overall, seven of the 11 morphospecies analysed were congruent with both the well-supported clades identified by phylogenetic analysis and the MOTUs inferred by species delimitation, while the remaining four morphospecies – E. carmelinoi, E. evandroi, E. lenti and E. piperiformis – were merged into a single well-supported clade/MOTU. Although E. carmelinoi, E. evandroi and E. lenti were indistinguishable using coi DNA barcodes, E. piperiformis did form a distinct phylogenetic cluster and could be correctly identified using best-match criteria. Despite their apparent morphological differences, we propose on the basis of the molecular similarity of their DNA barcodes that these latter four morphospecies should be considered members of a recently-diverged species complex.     

Phylogenetic taxonomy of Hymenochaete and related genera (Hymenochaetales)

Species delimitation in Hymenochaete sensu lato was evaluated based on partial nucLSU rDNA sequences, isoenzyme analyses and morphological data. Analyses of LSU data revealed Hymenochaete sensu stricto and Pseudochaete as distinct monophyletic genera. Transfer of Hydnochaete olivacea into genus Pseudochaete and Dichochaete setosa into genus Hymenochaete (as H. resupinata) was supported by the results of molecular analyses. A new species Hymenochaete koeljalgii from Tanzania was described. The species of genus Hymenochaete sensu stricto were divided into four well-supported clades possessing no distinctive morphological characters.     

Fine-scale spatial partitioning of genetic variation and evolutionary contestability in the invasive estuarine mussel Xenostrobus securis

Initial studies of Australian populations of the widely invasive mussel Xenostrobus securis raised the possibility of non-random spatial partitioning within estuaries of mitochondrial DNA (mtDNA) clades that have dispersed intercontinentally and those that have not. A fine-scale phylogeographic investigation was made to examine this possibility using cytochrome c oxidase DNA sequences and data from five microsatellite loci developed here. X. securis in the study region was found to comprise multiple robustly supported mtDNA clades, many of which were genetically very distinct. Frequently, the clades comprised numerous haplotypes that have narrow ranges and may have evolved in situ within drainages or nearby. Microsatellite data reveal significant intra- and inter-drainage differentiation but at lower levels than mtDNA. Variation in mtDNA appears to reflect not only recolonization after infrequent local extinction within drainages but also the low probability of migratory haplotypes successfully establishing during evolutionary contests, with resident haplotypes having selective or demographic advantages derived from local occupation. The patterns of mtDNA haplotype distribution suggest that central New South Wales is the source of the internationally invasive lineages of X. securis.