Population genetic structure and comparative diversity of smallmouth bass Micropterus dolomieu: congruent patterns from two genomes

Genetic diversity and divergence patterns of smallmouth bass Micropterus dolomieu spawning groups are analysed across its northern native range with mtDNA cytochrome b gene sequences and eight unlinked nuclear DNA microsatellite loci. Results reveal high levels of genetic variability and significant differences in allelic representation among populations (mtDNA: mean ± s.e ., H_D = 0·50 ± 0·06, mean ± s.e ., θ_ST = 0·41 ± 0·02 and microsatellites: mean ± s.e . H_O = 0·46 ± 0·03, mean ± s.e . θ_ST = 0·25 ± 0·01). The distributions of 28 variant mtDNA haplotypes, which differ by an average of 3·94 nucleotides (range = 1–8), denote divergent representation among geographic areas. Microsatellite data support nine primary population groups, whose high self‐assignment probabilities likewise display marked divergence. Genetic patterns demonstrate: (1) high genetic diversity in both genomes, (2) significant divergence among populations, probably resulting from natal site homing and low lifetime migration, (3) support for three post‐glacial refugia that variously contributed to the current northern populations, which remain evident today despite waterway connectivity and (4) a weak yet significant genetic isolation by geographic distance pattern, indicating that other processes affect the differences among populations, such as territoriality and site fidelity.     

No genetic differentiation between geographically isolated populations of Clarias macrocephalus Günther in Malaysia revealed by sequences of mtDNA Cytochrome b and D‐loop gene regions

In the present study, we assessed the genetic variation of three Clarias macrocephalus Günther populations collected from Kedah, Perlis and Kelantan (Peninsular Malaysia) using sequences of partial mitochondrial cytochrome b (Cyt b) and D‐loop genes. A total of 57 individuals were sequenced and 1470 bp were obtained (1053 bp Cyt‐b; 417 bp D‐loop). The analysis revealed 21 haplotypes based on 81 polymorphic sites. Nucleotide diversity (π) was 0.003 in all populations while haplotype diversity ranged from 0.657 to 0.765. No significant genetic differentiation among the three populations was observed. Nevertheless, a number of private haplotypes was discovered, providing valuable information for selective breeding programs.     

Low mtDNA variation and shallow population structure of the Chinese pomfret Pampus chinensis along the China coast

In the present study, population genetic structure and genetic diversity of the Chinese pomfret Pampus chinensis, along the China coast were investigated and compared with that from Indonesia using mitochondrial DNA cytochrome b gene sequences. A total of 28 variable sites (including 18 singleton sites and 10 parsimony information sites) were observed and 23 haplotypes were defined in 330 individuals from 11 localities. The haplotype diversity (H_D) of the populations ranged from 0·540 to 0·828, the nucleotide diversity (π) ranged from 0·081 to 0·295%. Pairwise F_ST statistics showed that significant genetic divergence occurred among populations from different geographical regions. The high dispersal capabilities, geographic segregation and ocean currents may be responsible for the present population genetic structure in this species. In addition, a population expansion event during the late Pleistocene period was inferred. The time of population expansion was estimated to occur about 117 000–169 000 years ago.     

Genetic bottleneck in invasive species: the potato tuber moth adds to the list

The level of genetic diversity within populations of introduced species has received increasing attention as an important factor influencing their survival and adaptive potential. We examined this issue with the Guatemalan potato tuber moth Tecia solanivora, an agricultural pest which has successfully invaded South America and the Canary Islands within the last 20 years. To analyse changes in T. solanivora genetic diversity, the mitochondrial marker cytochrome b was sequenced from individuals collected across its known distribution area. High haplotypic diversity was observed in Guatemala, whereas only three haplotypes have been found in Venezuela and a single one in the remaining invaded South American countries and the Canary Islands. Invasive haplotypes were not observed in our samples from Guatemala but are closely related to Guatemalan haplotypes. These results are consistent with the hypotheses that (i) either a few individuals were introduced into Venezuela leading to a strong initial genetic bottleneck, or selection pressure may have lead to the disparition of all but a reduced number of introduced haplotypes, (ii) a second bottleneck occured between Venezuela and Colombia, and (iii) the invasion of the Canary Islands originated from South America. We further reviewed the recent literature to compare this change in genetic diversity with those reported for other invasive species. We quantified the changes in genetic diversity between native and introduced ranges for 57 biological invasions. We found that the genetic homogenization in T. solanivora was among the strongest reported and discuss factors that can explain the success of invasive populations with low genetic diversity.     

Mitochondrial DNA variation in natural populations of endangered Indian Feather-Back Fish, Chitala chitala

Genetic variation at mitochondrial cytochrome b (cyt b) and D-loop region reveals the evidence of population sub-structuring in Indian populations of highly endangered primitive feather-back fish Chitala chitala. Samples collected through commercial catches from eight riverine populations from different geographical locations of India were analyzed for cyt b region (307 bp) and D-loop region (636–716 bp). The sequences of the both the mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The patterns of genetic diversity, haplotypes networks clearly indicated two distinct mitochondrial lineages and mismatch distribution strongly suggest a historical influence on the genetic structure of C. chitala populations. The baseline information on genetic variation and the evidence of population sub-structuring generated from this study would be useful for planning effective strategies for conservation and rehabilitation of this highly endangered species.     

Conservation genetics in hypersaline inland waters: mitochondrial diversity and phylogeography of an endangered Iberian beetle (Coleoptera: Hydraenidae)

Saline inland waters are globally threatened habitats harbouring many specialised endemic species, which often have restricted geographic ranges, and occur as highly isolated populations. We studied the genetic variation and phylogeography of Ochthebius glaber Montes and Soler, a rare and endangered water beetle endemic to hypersaline streams in the South and Southeast of the Iberian Peninsula. We used a 633 bp fragment of cytochrome oxidase subunit 1 gene to determine the genetic diversity and phylogeographic structure within this species, and interpret this in the light of the speciesȁ9 conservation requirements. Thirteen populations were sampled across the speciesȁ9 geographic range, and genetic diversity found to be very high, with 37 haplotypes across the 71 specimens examined (p-distance 0.2–7.3%, average 3.1±0.4). Phylogeographic analyses revealed a surprisingly high degree of geographical structure, detectable among populations separated by relatively short geographical distances, with three main groups of haplotypes which have apparently been isolated for significant periods of time. Past fragmentation and contiguous range expansion events were inferred as the main causes of the detected geographical associations of haplotypes. The establishment of independent evolutionary lineages as conservation units is particularly important for species inhabiting saline habitats such as O. glaber, which is endangered by habitat loss across most of its distribution. However, given the natural instability of hypersaline environments, the conservation of a network of populations and potential habitats would be necessary to enable the preservation of the process generating and maintaining the diversity of the species.     

Phylogeography of Podocarpus matudae (Podocarpaceae): pre‐Quaternary relicts in northern Mesoamerican cloud forests

Aim Cloud forests of northern Mesoamerica represent the northern and southern limit of the contact zone between species otherwise characteristic of North or South America. Several phylogeographic studies featuring temperate conifer species have improved our understanding of species responses to environmental changes. In contrast, conifer species that presumably colonized northern Mesoamerica from South America are far less studied. A phylogeographic study of Podocarpus matudae (Podocarpaceae) was conducted to identify any major evolutionary divergences or disjunctions across its range and to determine if its current distribution is associated with pre‐Quaternary climatic and/or long‐distance dispersal events. Location Northern Mesoamerica (Mexico and Guatemala). Methods Sixteen populations (157 individuals) of P. matudae were screened for variation at two plastid DNA markers. The intra‐specific phylogenetic relationships among haplotypes were reconstructed using Bayesian inference. Population genetic analyses were undertaken to gain insight into the evolutionary history of these populations. To test whether genetic divergence among populations occurred at different time‐scales plastid DNA sequence data and fossil‐ and coalescent‐based calibrations were integrated. Results The combination of plastid markers yielded 11 haplotypes. Differentiation among populations based on DNA variation (G_ST) (0.707, SE 0.0807) indicated a clear population structure in P. matudae. Differentiation for ordered alleles (N_ST) (0.811, SE 0.0732) was higher than that for G_ST, indicating phylogeographical structure in P. matudae. Most of the total variation (81.3%, P < 0.0001) was explained by differences among populations. The estimated divergence time between the unique haplotypes from a Guatemalan population and the two most common haplotypes from the Sierra Madre Oriental in Mexico was between 10 and 20 Ma, and further haplotype divergence in the poorly resolved clade of the Sierra Madre Oriental occurred between 3 and 0.5 Ma. Main conclusions Divergence estimations support the hypothesis that extant Podocarpus matudae populations are pre‐Quaternary relicts. This finding is consistent with fossil and pollen data that support a Miocene age for temperate floristic elements in Mesoamerican cloud forests, whereas further haplotype divergence within the Sierra Madre Oriental, Chiapas and Guatemala occurred more recently, coinciding with Pleistocene cloud forest refugia.     

Genetic diversity shaped by historical and recent factors in the live‐bearing twoline skiffia Neotoca bilineata

The endangered twoline skiffia Neotoca bilineata, a viviparous fish of the subfamily Goodeinae, endemic to central Mexico (inhabiting two basins, Cuitzeo and Lerma‐Santiago) was evaluated using genetic and habitat information. The genetic variation of all remaining populations of the species was analysed using both mitochondrial and microsatellite markers and their habitat conditions were assessed using a water quality index (I_WQ). An 80% local extinction was found across the distribution of N. bilineata. The species was found in three of the 16 historical localities plus one previously unreported site. Most areas inhabited by the remaining populations had I_WQ scores unsuitable for the conservation of freshwater biodiversity. Populations showed low but significant genetic differentiation with both markers (mtDNA φ_ST = 0·076, P < 0·001; microsatellite F_ST = 0·314, P < 0·001). Borbollon, in the Cuitzeo Basin, showed the highest level of differentiation and was identified as a single genetic unit by Bayesian assignment methods. Rio Grande de Morelia and Salamanca populations showed the highest genetic diversity and also a high migration rate facilitated by an artificial channel that connected the two basins. Overall, high genetic diversity values were observed compared with other freshwater fishes (average N_a = 16 alleles and loci and mean ±s.d . H_o = 0·63 ± 0·10 and nucleotide diversity π = 0·006). This suggests that the observed genetic diversity has not diminished as rapidly as the species' habitat destruction. No evidence of correlation between habitat conditions and genetic diversity was found. The current pattern of genetic diversity may be the result of both historical factors and recent modifications of the hydrological system. The main threat to the species may be the rapid habitat deterioration and associated demographic stochasticity rather than genetic factors.     

Assessing the use of mitochondrial cox1 gene and cox2-3 spacer for genetic diversity study of Malaysian Gracilaria changii (Gracilariaceae, Rhodophyta) from Peninsular Malaysia

Advances in DNA-based genetic markers provide the essential tools in measurement of genetic diversity relating to the evolution, biogeography, and systematics of red algae by exploiting genetic variation in the entire genome of organisms. The understanding of genetic diversity in Gracilaria changii (Gracilariaceae, Rhodophyta) will provide valuable information for conservation, plant breeding management, and strain selection for cultivation. However, information of intraspecific genetic variation is still rudimentary. In this study, two mitochondrial encoded markers, cytochrome oxidase subunit 1 (cox1) and intergenic spacer between the cytochrome oxidase subunits 2 and 3 (cox2-3 spacer) were used to investigate genetic diversity in 40 individuals of G. changii collected from 11 different geographically distinct populations from Peninsular Malaysia. Seven distinct mitochondrial haplotypes were identified with the cox1 gene and three mitochondrial haplotypes with the cox2-3 spacer. Intraspecific nucleotide differences ranged from 0 to 6 bp for the cox1 and 0–4 bp for the cox2-3 spacer, respectively. This is the first report comparing the suitability of mitochondrial markers of the cox1 gene and the cox2-3 spacer for genetic diversity studies on G. changii. The present study showed that the cox1 gene is a potential molecular marker to infer intraspecific genetic variation in red macroalgae. The cox1 marker is more variable compared to the cox2-3 spacer and revealed genetic variation and phylogeographic structure for this ecologically and economically important species.     

mtDNA evidence for a local northern latitude Pleistocene refugium for the root vole (Microtus oeconomus,Arvicolinae, Rodentia) from Eastern Poland

We analysed the genetic structure of 33 populations of the root vole (tundra vole, Microtus oeconomus, Pallas, 1776) inhabiting their typical habitats, located at different distances from the southern boundary of the species’ range (52°14′–53°56′ N) in eastern Poland. We determined its phylogeographic pattern as well as the possible occurrence of a small, local high‐latitude refugium of this species in southern Poland, previously suggested in palaeontological studies. 908 bp of cytochrome b sequences were analysed from 439 root voles, and 21 mtDNA cytb haplotypes belonging to the Central European (CE) phylogroup were found. Haplotype diversity in the examined populations varied between 0 and 0.872 (mean: 0.425 ± 0.332), while nucleotide diversity ranged between 0 and 0.62% (mean: 0.235% ± 0.217). Within the CE phylogroup of M. oeconomus, we identified with high bootstrap support a newly separated group of M. oeconomus that evolved from CE, denoted CE‐PL S. This group is located in the southern and central part of eastern Poland and most likely diverged from phylogroup CE in a small, cryptic refugium situated in southern Poland, in the Kraków‐Cz?stochowa Upland and/or the Holy Cross Mountains during the LGM and Younger Dryas.     

Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south‐eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non‐dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61–5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America–The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.     

Population structure of the olive flounder (Paralichthys olivaceus) in Korea inferred from microsatellite marker analysis

The population structure of olive flounder Paralichthys olivaceus was estimated using nine polymorphic microsatellite (MS) loci in 459 individuals collected from eight populations, including five wild and three hatchery populations in Korea. Genetic variation in hatchery (mean number of alleles per locus, A = 10·2–12·1; allelic richness, A_R = 9·3–10·1; observed heterozygosity, H_O = 0·766–0·805) and wild (mean number of alleles per locus, A = 11·8–19·6; allelic richness, A_R = 10·9–16·1; observed heterozygosity, H_O = 0·820–0·888) samples did not differ significantly, suggesting a sufficient level of genetic variation in these well‐managed hatchery populations, which have not lost a substantial amount of genetic diversity. Neighbour‐joining tree and principal component analyses showed that genetic separation between eastern and pooled western and southern wild populations in Korea was probably influenced by restricted gene flow between regional populations due to the barrier effects of sea currents. The pooled western and southern populations are genetically close, perhaps because larval dispersal may depend on warm currents. One wild population (sample from Wando) was genetically divergent from the main distribution, but it was genetically close to hatchery populations, indicating that the genetic composition of the studied populations may be affected by hydrographic conditions and the release of fish stocks. The estimated genetic population structure and potential applications of MS markers may aid in the proper management of P. olivaceus populations.     

Population diversity of Cistopus indicus inferred from mitochondrial DNA (Cytochrome b) variation from China and Vietnam

ABSTRACT

The octopus Cistopus indicus is an important target of cephalopod fisheries in China. It is widely distributed in the South Pacific and tropical Indian Ocean, from the South China Sea, the Philippines, Malaysia, to Indian and Pakistan seas. We collected specimens from five sites in China and Vietnam (Zhoushan, Wenzhou, Shacheng, Zhanjiang and Mangjie). A fragment of 675bp of cytochrome b (Cytb) was amplified from 95 individuals. A total of 27 haplotypes and 78 variable nucleotide sites was observed. High haplotype diversity and low nucleotide diversity were observed in all populations. The phylogenetic analysis separated these populations into two clades; one was composed of three populations (Zhoushan, Wenzhou and Shacheng), the other of two (Zhanjiang, Mangjie). AMOVA analysis detected that 4.67% of the genetic variation occurred within populations and 95.33% occurred among populations. FST values ranged from 0.014 to 0.993, highlighting the high genetic variation among the populations. Assuming a molecular clock with a rate of 2.15–2.6%/Ma for the Cytb gene, the two clades may have diverged 2.88–3.49 million years ago (Pliocene). Neutral evolution tests and mismatch distribution analysis suggested recent population expansion. The present results revealed valuable information for genetic assessment, management and conservation of this species.     

Population genetic structure and larval dispersal potential of spottedtail goby Synechogobius ommaturus in the north‐west Pacific

Larval dispersal may have an important effect on genetic structure of benthic fishes. To examine the population genetic structure of spottedtail goby Synechogobius ommaturus, a 478 base pair (bp) fragment of the hypervariable portion of the mtDNA control region was sequenced and used to interpret life‐history characteristics and larval dispersal strategy. Individuals (n = 186) from 10 locations on the coasts of China and Korea were analysed and 44 haplotypes were obtained. The levels of haplotype and nucleotide diversity were higher in East China Sea populations than in other populations. Both the phylogenetic tree and the minimum spanning tree showed that no significant genealogical structures corresponding to sampling locations existed. AMOVA and pair‐wise F_ST revealed significant genetic differentiation between populations from Korea and China. A significant isolation by distance pattern was observed in this species (r = 0·53, P < 0·001). Both mismatch distribution analysis and neutrality tests showed S. ommaturus to have experienced a recent population expansion. These results suggest that the Pleistocene ice ages had a major effect on the phylogeographic pattern of S. ommaturus, that larvae might avoid offshore dispersal and that dispersal of larvae may maintain a migration–drift equilibrium.     

Genetic diversity and maternal origin of Northeast African and South American donkey populations

To investigate the mtDNA variation and origin of maternal lineages in South American donkeys and to reassess the domestication of donkeys in northeast Africa, we analyzed sequences (489 bp of the D‐loop) from 323 domestic donkeys sampled from Peru, Brazil, Ethiopia and Egypt. Altogether, the 323 sequences displayed 53 different haplotypes (45 in Ethiopia, 14 in Egypt, eight in Peru and six in Brazil). Among the four populations, Egyptian donkeys possessed the highest haplotype diversity (0.910 ± 0.032), followed by Brazilian donkeys (0.879 ± 0.060). The Clade I haplotypes dominated in Peruvian donkeys (65%), whereas Clade II haplotypes dominated in Brazilian donkeys (67%). Estimates of F_ST values showed a high genetic differentiation between Peruvian and Brazilian donkey populations (F_ST = 0.4066), which could be explained by the complex introduction history of South American donkeys. Phylogeographic analysis indicates that northeast Africa could be the most probable domestication center for Clade I donkeys. Analysis of molecular variance confirmed a weak genetic structure in domestic donkey populations among four continents (Europe, Asia, Africa and South America).     

Phylogeographic differentiation in the mitochondrial control region in the koala, Phascolarctos cinereus (Goldfuss 1817)

The koala, Phascolarctos cinereus, is a geographically widespread species endemic to Australia, with three currently recognized subspecies: P.c. adustus, P.c. cinereus, and P.c. victor. Intraspecific variation in the mitochondrial DNA (mtDNA) control region was examined in over 200 animals from 16 representative populations throughout the species’ range. Eighteen different haplotypes were defined in the ≈ 860 bp mtDNA control region, as determined by heteroduplex analysis/temperature gradient gel electrophoresis (HDA/TGGE). Any single population typically possessed only one or two haplotypes yielding an average within-population haplotypic diversity of 0.180 ± 0.003, and nucleotide diversity of 0.16%. Overall, mtDNA control region sequence diversity between populations averaged 0.67%, and ranged from 0% to 1.56%. Nucleotide divergence between populations averaged 0.51%, and ranged from 0% to 1.53%. Neighbour-joining methods revealed limited phylogenetic distinction between geographically distant populations of koalas, and tentative support for a single evolutionarily significant unit (ESU). This is consistent with previous suggestions that the morphological differences formalized by subspecific taxonomy may be interpreted as clinal variation. Significant differentiation in mtDNA-haplotype frequencies between localities suggested that little gene flow currently exists among populations. When combined with microsatellite analysis, which has revealed substantial differentiation among koala populations, we conclude that the appropriate short-term management unit (MU) for koalas is the local population.     

Analysis of Amblyomma sculptum haplotypes in an area endemic for Brazilian spotted fever

Amblyomma sculptum (Ixodida: Ixodidae) Berlese, 1888, a member of the Amblyomma cajennense complex, is the major vector of Brazilian spotted fever (BSF) in southeastern Brazil. In this study, the genetic diversity of A. sculptum populations in the state of Rio de Janeiro (RJ), Brazil, was investigated because genetic variability in tick populations may be related to vector competence. Samples of A. sculptum from 19 municipalities in 7 regions of RJ were subjected to DNA extraction, amplification and sequencing of D‐loop, cytochrome oxidase II and 12S rDNA mitochondrial genes. These sequences were used to map the genetic diversity of this tick. Amblyomma sculptum populations are genetically diverse in RJ, especially in the South Centre and Highland regions. Few unique haplotypes were observed in all populations, and the majority of genetic variation found was among ticks within each population. Phylogenetic reconstruction reinforced the assumption that all the haplotypes identified in RJ belong to A. sculptum. However, some RJ haplotypes are closer to A. sculptum from Argentina than to A. sculptum from elsewhere in Brazil. In RJ, A. sculptum has high genetic diversity, although little genetic differentiation. Observations also indicated a high level of gene flow among the studied populations and no evidence of population structure according to region in RJ.     

Mitochondrial population structure and post‐glacial dispersal of longnose sucker Catostomus catostomus in Labrador,Canada: evidence for multiple refugial origins and limited ongoing gene flow

Two hundred and eighty‐seven longnose sucker Catostomus catostomus were collected from 14 lakes in Labrador, 52 from three lakes in Ontario, 43 from two lakes in British Columbia and 32 from a lake in Yukon; a total of 414 in all. The resulting 34 haplotypes (20 in Labrador) contained moderate haplotypic diversity (h = 0·657) and relatively low nucleotide diversity (π = 3·730 × 10^?3. Mean ?_ST (0·453, P < 0·05) over all populations revealed distinct genetic structuring among C. catostomus populations across Canada, based on province, which was validated by the analysis and spatial analysis of molecular variance (c. 80% variation between provinces). These results probably reflect the historical imprint of recolonization from different refugia and possibly indicate limited ongoing gene flow within provinces. A haplotype network revealed one major and two minor clades within Labrador that were assigned to the Atlantic, Beringian and Mississippian refugia, respectively, with tests of neutrality and mismatch distribution indicative of a recent population expansion in Labrador, dated between c. 3500 and 8300 years ago.     

High gene flow in Girella punctata (Perciformes,Kyphosidae) among the Japanese Islands inferred from partial sequence of the control region in mitochondrial DNA

Sequencing analysis of the partial control region (c. 450 bp) detected 88 haplotypes from 249 individuals of Girella punctata collected from coastal waters of nine locations in the Japanese Islands. A single haplotype was the most numerous at all sampling locations, and no significant genetic difference was found among G. punctata samples collected from various locations (F_ST=?0·0274 to 0·0247) with high haplotype diversity and low nucleotide diversity, suggesting that G. punctata population had experienced a bottleneck followed by rapid population growth around the Japanese Islands and might be affected by larval dispersal in association with warm currents.