New insights into the origin and the genetic status of the Balkan donkey from Serbia

The Balkan donkey (Equus asinus L.) is commonly regarded as a large‐sized, unselected, unstructured and traditionally managed donkey breed. We assessed the current genetic status of the three largest E. asinus populations in the central Balkans (Serbia) by analysing the variability of nuclear microsatellites and the mitochondrial (mtDNA) control region of 77 and 49 individuals respectively. We further analysed our mtDNA dataset along with 209 published mtDNA sequences of ancient and modern individuals from 19 European and African populations to provide new insights into the origin and the history of the Balkan donkey. Serbian donkey populations are highly genetically diverse at both the nuclear and mtDNA levels despite severe population decline. Traditional Balkan donkeys in Serbia are rather heterogeneous; we found two groups of individuals with similar phenotypic features, somewhat distinct nuclear backgrounds and different proportions of mtDNA haplotypes belonging to matrilineal Clades 1 and 2. Another group, characterized by larger body size, different coat colour, distinct nuclear gene pool and predominantly Clade 2 haplotypes, was delineated as the Banat donkey breed. The maternal landscape of the large Balkan donkey population is highly heterogeneous and more complex than previously thought. Given the two independent domestication events in donkeys, multiple waves of introductions into the Balkans from Greece are hypothesized. Clade 2 donkeys probably appeared in Greece prior to those belonging to Clade 1, whereas expansion and diversification of Clade 1 donkeys within the Balkans predated that of Clade 2 donkeys.     

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).     

Mitochondrial DNA variation of Nigerian domestic helmeted guinea fowl

We analyzed genetic diversity of 215 mitochondrial DNA (mtDNA) D‐loop sequences from seven populations of domesticated helmeted guinea fowl (Numida meleagris) in Nigeria and compared that with results of samples collected in Kenya (n = 4) and China (n = 22). In total, 241 sequences were assigned to 22 distinct haplotypes. Haplotype diversity in Nigeria was 0.693 ± 0.022. The network grouped most matrilines into two main haplogroups: A and B. There was an absence of a geographic signal, and two haplotypes dominated across all locations with the exception of the Kebbi population in the northwest of the country; AMOVA also confirmed this observation (F_ST = 0.035). The low genetic diversity may be a result of recent domestication, whereas the lack of maternal genetic structure likely suggests the extensive genetic intermixing within the country. Additionally, the differentiation of the Kebbi population may be due to a certain demographic history and/or artificial selection that shaped its haplotype profile. The current data do not permit us to make further conclusions; therefore, more research evidence from genetics and archaeology is still required.     

Genetic variation of the East Balkan Swine (Sus scrofa) in Bulgaria,revealed by mitochondrial DNA and Y chromosomal DNA

East Balkan Swine (EBS) Sus scrofa is the only aboriginal domesticated pig breed in Bulgaria and is distributed on the western coast of the Black Sea in Bulgaria. To reveal the breed's genetic characteristics, we analysed mitochondrial DNA (mtDNA) and Y chromosomal DNA sequences of EBS in Bulgaria. Nucleotide diversity (π_n) of the mtDNA control region, including two newly found haplotypes, in 54 EBS was higher (0.014 ± 0.007) compared with that of European (0.005 ± 0.003) and Asian (0.006 ± 0.003) domestic pigs and wild boar. The median‐joining network based on the mtDNA control region showed that the EBS and wild boar in Bulgaria comprised mainly two major mtDNA clades, European clade E1 (61.3%) and Asian clade A (38.7%). The coexistence of two mtDNA clades in EBS in Bulgaria may be the relict of historical pig translocation. Among the Bulgarian EBS colonies, the geographical differences in distribution of two mtDNA clades (E1 and A) could be attributed to the source pig populations and/or historical crossbreeding with imported pigs. In addition, analysis of the Y chromosomal DNA sequences for the EBS revealed that all of the EBS had haplotype HY1, which is dominant in European domestic pigs.     

Multisource noninvasive genetics of brown bears (Ursus arctos) in Greece reveals a highly structured population and a new matrilineal contact zone in southern Europe

In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (N _e) was estimated to be 199 individuals, and total combined population size (N _C) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe''s low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Maternal and paternal genetic diversity of ancient sheep in Estonia from the Late Bronze Age to the post‐medieval period and comparison with other regions in Eurasia

Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D‐loop of mitochondrial DNA (mtDNA; 213 bp) and the Y‐chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G‐oY1 of the Y‐chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.     

Limited mitochondrial DNA variation within South Africa's black rhino (Diceros bicornis minor) population and implications for management

The taxonomy of African black rhinoceros (Diceros bicornis) remains unresolved. Maintaining levels of genetic diversity and species rescue by reintroduction and restocking requires its resolution. We compared the sequences of the mitochondrial DNA (mtDNA) control region for a total of 101 D. bicornis from three subspecies: D. b. minor, D. b. michaeli and D. b. bicornis. A single unique haplotype was found within the 65 D. b. minor samples from KwaZulu‐Natal (KZN) Province, South Africa, 55 of which came from Hluhluwe‐iMfolozi Game Park (HiP) and Mkuzi Game Reserve (MGR) source populations. However, six different haplotypes were represented in eleven D. b. minor samples from Zimbabwe. Similarly, published autosomal microsatellite data indicate low levels of diversity within the KZN D. b. minor populations. The low levels of mtDNA diversity within the KZN metapopulation point to the possible need for genetic supplementation. However, there is a need to determine whether the low levels of genetic variation within KZN D. b. minor are a result of the recent bottleneck or whether KZN historically always had low diversity.     

Genetic diversity loss in a biodiversity hotspot: ancient DNA quantifies genetic decline and former connectivity in a critically endangered marsupial

The extent of genetic diversity loss and former connectivity between fragmented populations are often unknown factors when studying endangered species. While genetic techniques are commonly applied in extant populations to assess temporal and spatial demographic changes, it is no substitute for directly measuring past diversity using ancient DNA (aDNA). We analysed both mitochondrial DNA (mtDNA) and nuclear microsatellite loci from 64 historical fossil and skin samples of the critically endangered Western Australian woylie (Bettongia penicillata ogilbyi), and compared them with 231 (n = 152 for mtDNA) modern samples. In modern woylie populations 15 mitochondrial control region (CR) haplotypes were identified. Interestingly, mtDNA CR data from only 29 historical samples demonstrated 15 previously unknown haplotypes and detected an extinct divergent clade. Through modelling, we estimated the loss of CR mtDNA diversity to be between 46% and 91% and estimated this to have occurred in the past 2000–4000 years in association with a dramatic population decline. In addition, we obtained near‐complete 11‐loci microsatellite profiles from 21 historical samples. In agreement with the mtDNA data, a number of ‘new’ microsatellite alleles was only detected in the historical populations despite extensive modern sampling, indicating a nuclear genetic diversity loss >20%. Calculations of genetic diversity (heterozygosity and allelic rarefaction) showed that these were significantly higher in the past and that there was a high degree of gene flow across the woylie's historical range. These findings have an immediate impact on how the extant populations are managed and we recommend the implementation of an assisted migration programme to prevent further loss of genetic diversity. Our study demonstrates the value of integrating aDNA data into current‐day conservation strategies.     

Mitochondrial DNA D‐loop sequence variation in maternal lineages of Iranian native horses

To understand the origin and genetic diversity of Iranian native horses, mitochondrial DNA (mtDNA) D‐loop sequences were generated for 95 horses from five breeds sampled in eight geographical locations in Iran. Sequence analysis of a 247‐bp segment revealed a total of 27 haplotypes with 38 polymorphic sites. Twelve of 19 mtDNA haplogroups were identified in the samples. The most common haplotypes were found within haplogroup X2. Within‐population haplotype and nucleotide diversities of the five breeds ranged from 0.838 ± 0.056 to 0.974 ± 0.022 and 0.011 ± 0.002 to 0.021 ± 0.001 respectively, indicating a relatively high genetic diversity in Iranian horses. The identification of several ancient sequences common between the breeds suggests that the lineage of the majority of Iranian horse breeds is old and obviously originated from a vast number of mares. We found in all native Iranian horse breeds lineages of the haplogroups D and K, which is concordant with the previous findings of Asian origins of these haplogroups. The presence of haplotypes E and K in our study also is consistent with a geographical west–east direction of increasing frequency of these haplotypes and a genetic fusion in Iranian horse breeds.     

Phylogeography of roe deer (Capreolus capreolus) populations: the effects of historical genetic subdivisions and recent nonequilibrium dynamics

We sequenced 704 mitochondrial DNA (mtDNA) control-region nucleotides and genotyped 11 autosomal microsatellites (STR) in 617 European roe deer (Capreolus capreolus) samples, aiming to infer the species' phylogeographical structure. The mtDNA sequences were split in three distinct haplogroups, respectively, named: Clade West, sampled mainly in Iberia; Clade East, sampled mainly in Greece and in the Balkans; and Clade Central, which was widespread throughout Europe, including the eastern countries and Iberia, but not Greece. These clades might have originated in distinct Iberian and Balkanic refuges during the penultimate or the last glaciations. Clades East and West contributed little to the current postglacial mtDNA diversity in central Europe, which apparently was recolonized mainly by haplotypes belonging to Clade Central. A unique subclade within Clade Central grouped all the haplotypes sampled from populations of the Italian subspecies C. c. italicus. In contrast, haplotypes sampled in central and southern Spain joined both Clade Central and Clade West, suggesting that subspecies C. c. garganta has admixed origin. STR data support a genetic distinction of peripheral populations in north Iberia and southern Italy, and show the effects of anthropogenic disturbance in fragmented populations, which were recently reintroduced or restocked and not may be in mutation-drift equilibrium. Roe deer in central Europe are mainly admixed, while peripheral populations in north Portugal, the southern Italian Apennines and Greece represent the remains of refugial populations and should be managed accordingly.     

Low Genetic Diversity and Strong Geographical Structure of the Critically Endangered White-Headed Langur (Trachypithecus leucocephalus) Inferred from Mitochondrial DNA Control Region Sequences

Many Asian colobine monkey species are suffering from habitat destruction and population size decline. There is a great need to understand their genetic diversity, population structure and demographic history for effective species conservation. The white-headed langur (Trachypithecus leucocephalus) is a Critically Endangered colobine species endemic to the limestone karst forests in southwestern China. We analyzed the mitochondrial DNA (mtDNA) control region sequences of 390 fecal samples from 40 social groups across the main distribution areas, which represented one-third of the total extant population. Only nine haplotypes and 10 polymorphic sites were identified, indicating remarkably low genetic diversity in the species. Using a subset of 77 samples from different individuals, we evaluated genetic variation, population structure, and population demographic history. We found very low values of haplotype diversity (h = 0.570 ± 0.056) and nucleotide diversity (π = 0.00323 ± 0.00044) in the hypervariable region I (HVRI) of the mtDNA control region. Distribution of haplotypes displayed marked geographical pattern, with one population (Chongzuo, CZ) showing a complete lack of genetic diversity (having only one haplotype), whereas the other population (Fusui, FS) having all nine haplotypes. We detected strong population genetic structure among habit patches (Φ _ST = 0.375, P < 0.001). In addition, the Mantel test showed a significant correlation between the pairwise genetic distances and geographical distances among social groups in FS (correlation coefficient = 0.267, P = 0.003), indicting isolation-by-distance pattern of genetic divergence in the mtDNA sequences. Analyses of demographic history suggested an overall stable historical population size and modest population expansion in the last 2,000 years. Our results indicate different genetic diversity and possibly distinct population history for different local populations, and suggest that CZ and FS should be considered as one evolutionarily significant unit (ESU) and two management units (MUs) pending further investigation using nuclear markers.     

中国甜菜夜蛾地理种群的遗传分化与基因流

【目的】为了明确我国甜菜夜蛾Spodoptera exigua地理种群间的遗传分化及基因流,阐明该种害虫在我国的种群历史动态。【方法】本研究对采自我国20个地理种群的529头甜菜夜蛾样品进行线粒体COI基因序列测定与分析,利用DnaSP 5.0和Arlequin 3.11软件分析种群间遗传多样性、遗传分化、基因流水平及分子变异,构建了单倍型系统发育树与单倍型网络图。【结果】在所分析的所有529个序列样本中,共检测出10个单倍型,其中Hap_1为所有种群所共享。总群体遗传多样性较低（Hd=0.257±0.025,Pi=0.0007±0.0001,Kxy=0.323）,群体间遗传分化较小（F_ST=0.211）,基因流水平较高（N_m=1.870）。AMOVA分析表明,甜菜夜蛾遗传变异主要来自种群内,种群间变异水平较低。各种群间遗传分化程度与地理距离无显著相关性（R²=0.005,P>0.05）。各单倍型相互散布在不同种群中,未形成明显系统地理结构。中性检验（Tajima’s D=-2.177, P<0.05; Fu’s F_S=-8.629, P<0.05)与错配分布分析表明,我国甜菜夜蛾种群曾经历种群近期扩张。【结论】研究结果揭示,甜菜夜蛾各种群间的基因交流并未受到地理距离的影响,验证了甜菜夜蛾具有高度的迁飞能力。     

Population genetic structure and intraspecific genetic distance of Periplaneta americana (Blattodea: Blattidae) based on mitochondrial and nuclear DNA markers

The American cockroach (Periplaneta americana) is a globally invasive pest that can cause significant economic loss and threaten human health. Although it is abundant and lives in close proximity to humans, few studies have investigated the genetic diversity of P. americana. Our study analyzed 1,053 P. americana and other Periplaneta species' samples from different locations in China and the United States. A traditional tree‐based method using 17 unique mitochondrial COI haplotypes of P. americana and 20 haplotypes of the other Periplaneta species accurately identified P. americana with a barcoding threshold of 5.1%. To identify the population genetic structure of P. americana, we investigated wingless gene and pooled them with obtained mtDNA data for a combined analysis. Although the genetic diversity of the USA group was relatively higher than the China group, the number of haplotypes and alleles of both groups was small. The analysis of molecular variance (AMOVA), intraspecific phylogeny, and haplotype networks indicated that P. americana had very little global genetic differentiation. The weak geographic genetic structure might reflect the human‐mediated dispersal of P. americana. Despite no apparent phylogeographic assignment of mtDNA and nuclear lineages was observed in both BI trees, the integrated COI sequence data identified four distinct P. americana haplotype groups, showing four ancient maternal lineages of P. americana in China and the United States.     

Analysis of mitochondrial DNA in Bolivian llama,alpaca and vicuna populations: a contribution to the phylogeny of the South American camelids

The objectives of this work were to assess the mtDNA diversity of Bolivian South American camelid (SAC) populations and to shed light on the evolutionary relationships between the Bolivian camelids and other populations of SACs. We have analysed two different mtDNA regions: the complete coding region of the MT‐CYB gene and 513 bp of the D‐loop region. The populations sampled included Bolivian llamas, alpacas and vicunas, and Chilean guanacos. High levels of genetic diversity were observed in the studied populations. In general, MT‐CYB was more variable than D‐loop. On a species level, the vicunas showed the lowest genetic variability, followed by the guanacos, alpacas and llamas. Phylogenetic analyses performed by including additional available mtDNA sequences from the studied species confirmed the existence of the two monophyletic clades previously described by other authors for guanacos (G) and vicunas (V). Significant levels of mtDNA hybridization were found in the domestic species. Our sequence analyses revealed significant sequence divergence within clade G, and some of the Bolivian llamas grouped with the majority of the southern guanacos. This finding supports the existence of more than the one llama domestication centre in South America previously suggested on the basis of archaeozoological evidence. Additionally, analysis of D‐loop sequences revealed two new matrilineal lineages that are distinct from the previously reported G and V clades. The results presented here represent the first report on the population structure and genetic variability of Bolivian camelids and may help to elucidate the complex and dynamic domestication process of SAC populations.     

Ancient vs. recent processes as factors shaping the genetic variation of the European wild boar: are the effects of the last glaciation still detectable?

The European wild boar is an important game species, subjected to local extinctions and translocations in the past, and currently enormously and worryingly expanding in some areas where management is urgently required. Understanding the relative roles of ancient and recent events in shaping the genetic structure of this species is therefore not only an interesting scientific issue, but it represents also the basis for addressing future management strategies. In addition, several pig breeds descend from the European wild boar, but the geographical location of the domestication area(s) and the possible introgression of pig genomes into wild populations are still open questions. Here, we analysed the genetic variation in different wild boar populations in Europe. Ten polymorphic microsatellites were typed in 252 wild boars and the mtDNA control region was sequenced in a subset of 145 individuals. Some samples from different pig breeds were also analysed. Our results, which were obtained considering also 612 published mtDNA sequences, suggest that (i) most populations are similarly differentiated, but the major discontinuity is found along the Alps; (ii) except for the Italian populations, European wild boars show the signature of a postglacial demographic expansion; (iii) Italian populations seem to preserve a high proportion of preglaciation diversity; (iv) the demographic decline which occurred in some areas in the last few centuries did not produce a noticeable reduction of genetic variation; (v) signs of human-mediated gene flow among populations are weak, although in some regions the effects of translocations are detectable and a low degree of pig introgression can be identified; (vi) the hypothesis of an independent domestication centre in Italy is not supported by our data, which in turn confirm that Central European wild boar might have represented an important source for domestic breeds. We can therefore conclude that recent human activities had a limited effect on the wild boar genetic structure. It follows that areas with high variation and differentiation represent natural reservoirs of genetic diversity to be protected avoiding translocations. In this context controlling some populations by hunting is not expected to affect significantly genetic variation in this species.     

Genetic connectivity across marginal habitats: the elephants of the Namib Desert

Locally isolated populations in marginal habitats may be genetically distinctive and of heightened conservation concern. Elephants inhabiting the Namib Desert have been reported to show distinctive behavioral and phenotypic adaptations in that severely arid environment. The genetic distinctiveness of Namibian desert elephants relative to other African savanna elephant (Loxodonta africana) populations has not been established. To investigate the genetic structure of elephants in Namibia, we determined the mitochondrial (mt) DNA control region sequences and genotyped 17 microsatellite loci in desert elephants (n = 8) from the Hoanib River catchment and the Hoarusib River catchment. We compared these to the genotypes of elephants (n = 77) from other localities in Namibia. The mtDNA haplotype sequences and frequencies among desert elephants were similar to those of elephants in Etosha National Park, the Huab River catchment, the Ugab River catchment, and central Kunene, although the geographically distant Caprivi Strip had different mtDNA haplotypes. Likewise, analysis of the microsatellite genotypes of desert‐dwelling elephants revealed that they were not genetically distinctive from Etosha elephants, and there was no evidence for isolation by distance across the Etosha region. These results, and a review of the historical record, suggest that a high learning capacity and long‐distance migrations allowed Namibian elephants to regularly shift their ranges to survive in the face of high variability in climate and in hunting pressure.     

Population structure in two sympatric species of the Lake Tanganyika cichlid tribe Eretmodini: evidence for introgression

Patterns of genetic differentiation were analysed and compared in two sympatric species of the endemic Lake Tanganyika cichlid tribe Eretmodini by means of mitochondrial DNA (mtDNA) sequences of the control region and six microsatellite DNA loci. The sample area covers a total of 138 km of mostly uninterrupted rocky shoreline in the Democratic Republic of Congo and includes the entire distribution range of Tanganicodus cf. irsacae that stretches over a distance of 35 km. Both markers detected significant genetic differentiation within and between the two species. T. cf. irsacae contained lower overall genetic variation than Eretmoduscyanostictus, possibly due to its more restricted range of distribution and its smaller effective population sizes. Complete fixation of Tanganicodus mtDNA haplotypes was observed in Eretmodus at two localities, while at two other localities some Tanganicodus individuals possessed Eretmodus mtDNA haplotypes. Taking into account the relatively large average sequence divergence of 6.2% between the two species, as well as the geographical distribution of mtDNA haplotypes in the lake, the observed pattern is more likely to be a consequence of asymmetric introgression than of shared ancestral polymorphism. As there is significant population differentiation between sympatric Tanganicodus and Eretmodus populations, the events of introgressions may have happened after secondary contact, but our data provide no evidence for ongoing gene flow and suggest that both species are reproductively isolated at present time.     

Lack of genetic structure in endangered large yellow croaker Larimichthys crocea from China inferred from mitochondrial control region sequence data

The large yellow croaker Larimichthys crocea is an important commercial marine fish species in China. However, information about the population structure of this species is limited. In the present study, mitochondrial DNA (mtDNA) control region was sequenced from four populations of the yellow croaker in the southern Yellow Sea, East China Sea and South China Sea to investigate the genetic structure of this species. A total of 54 haplotypes were detected among 62 individuals of large yellow croaker. High levels of population genetic diversity were observed. Among the four populations, the haplotype diversity was between 0.9895 ± 0.0193 (Xiamen) and 1.0000 ± 0.0524 (Lvsi, Zhoushan). The nucleotide diversity ranged from 0.0208 ± 0.0108 (Xuwen) to 0.0246 ± 0.0138 (Lvsi). The results of AMOVA detected no significant differences among populations. The conventional F_ST statistics were negative and insignificant values. These indicated lack of population genetic structure throughout the Yellow Sea, East China Sea and South China Sea, and random mixing of individuals among the samples. Biological characters of large yellow croaker and lack of physical barrier in the studied area could be the reasons for lack of genetic structure in this species.