Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus)

Multiple genetic lineages of anadromous migratory Mekong catfish Pangasius krempfi revealed by mtDNA control region and cytochrome b

Population genetic structure of migratory fishes can reflect ecological and evolutionary processes. Pangasius krempfi is a critically important anadromous catfish in the Mekong River, and its migration pathways and genetic structure have attracted much interest. To investigate, we quantified the genetic diversity of this species using the control region (D-loop) and Cytochrome b (Cytb) of the mitochondrial genome. Fish were sampled (n = 91) along the Mekong tributaries from upstream to estuaries and coastal areas in the Mekong Delta and compared to three samples from Pakse (Laos). The D-loop haplotype (0.941 ± 0.014) and nucleotide diversity (0.0083 ± 0.0005) were high in all populations, but that of Cytb was low (0.331 ± 0.059 and 0.00063 ± 0.00011, respectively). No genetic difference was detected between populations, indicating strong gene flow and confirming a long migration distance for this species. Pangasius krempfi was not genetically structured according to geographical populations but was delineated into three haplogroups, suggesting multiple genetic lineages. The presence of haplogroups in each sampling location implies that migration downstream is random but parallel when the fish enter two river tributaries bifurcating from the main Mekong River. Individuals can also migrate along the coast, far from the estuaries, suggesting a longer migration path than previously reported, which is crucial for maintaining diverse genetic origin and migration pathways for P. krempfi.

Phylogenetic analyses of Eurasian lynx (Lynx lynx Linnaeus, 1758) including new mitochondrial DNA sequences from Iran

The Eurasian lynx (Lynx lynx) is one of the widespread felids in Eurasia; however, relatively little is known about the Asian subspecies, and especially the Iranian populations, which comprise the most southwestern part of its range. The current study aimed to assess the phylogenetic status of Iranian populations relative to other populations of Eurasia, by sequencing a 613 bp fragment of the mitochondrial control region. In total, 44 haplotypes were recorded from 83 sequences throughout Eurasia, two of which were found in Iran. The haplotype (H1) is dominant in all Iranian lynx populations and identical to specimens from SW Russia and central China. The second haplotype (H2) is unique and was recorded only from Ghazvin Province in the central Alborz Mountains. Both haplotypes occur in Ghazvin Province. The phylogenetic tree and a median-joining network identified four clades (i.e., East, West 1, West 2, and South). These results are congruent with previous studies and suggest that Eurasian lynx was restricted to the southern part of its range during the glacial maxima and expanded from there to East Asia and to Europe during several independent re-colonization events. The Caucasus region most like plays an important role as a refugium during glacial cycles.

Haplotypes traceability and genetic variability of the breeding population of pacu (Piaractus mesopotamicus) revealed by mitochondrial DNA

The main objective of this study was to estimate the genetic diversity levels and haplotype traceability in pacu Piaractus mesopotamicus from the breeding program located in Brazil by analyses of the mitochondrial DNA control region (mtDNA). Moreover, broodstocks from eight commercial fish farms were used for comparative evaluation, four from Brazil (Br1-Br4) and four from Argentina (Ar1-Ar4). The descriptive results revealed 47 polymorphic sites and 51 mutations, which evidenced 34 haplotypes. Ten haplotypes were shared among fish farms and 24 were exclusive. The nucleotide diversity (π) ranged from 0.00031 to 0.01462 and haplotype diversity (Hd) from 0.125 to 0.868. The analysis of molecular variance (AMOVA) indicated high structure present in the analyzed stocks (FST = 0.13356 and ФST = 0.52707). The genetic diversity was high in most of the commercial broodstocks, especially those from Brazil. We observed seven haplotypes in the genetic breeding population, of which four were exclusive and three shared among the commercial fish farms. The genetic diversity was moderate (π = 0.00265 and Hd = 0.424) and considered appropriated for this breeding population of pacu. Our results provide support for the genetic diversity maintenance and mtDNA traceability of pacu commercial broodstocks.

The complete mitochondrial genomes of three imperiled cyprinid fishes Bonytail (Gila elegans), Rio Grande Silvery Minnow (Hybognathus amarus) and Loach Minnow (Tiaroga cobitis)

Gila elegans, Hybognathus amarus, and Tiaroga cobitis (Family Cyprinidae, Order Cypriniformes) are endemic and endangered fishes in the southwestern United States. We present complete mitochondrial genomes for each species. Each mitochondrion consisted of 13 protein-coding genes, 2 ribosomal (rRNA) genes, 22 transfer RNA (tRNA) genes, and a single control region (D-loop), and gene order was consistent with other cyprinid fishes. Total genome lengths were 16,593 base pairs (bp) for G. elegans, 16,705 bp for H. amarus, and 16,802 for T. cobitis. The GC content in G. elegans and H. amarus was 44%, but higher in T. cobitis at 48%. Phylogenetic trees were generated to confirm relationships inferred via novel mitogenomes, and best-supported trees were consistent with previous research.

The spatio-temporal features of chicken mitochondrial ND2 gene heteroplasmy and the effects of nutrition factors on this gene

Mitochondrial heterogeneity is the presence of two or more types of mitochondrial (mt)DNA in the same individual/tissue/cell. It is closely related to animal health and disease. ND2 is a protein-coding gene in mtDNA, which participates in mitochondrial respiratory chain and oxidative phosphorylation. In previous studies, we observed that the mt.A5703T and mt.T5727G sites in the ND2 gene were the heteroplasmic variation sites. We used pyrophosphate sequencing technology to examine chicken mt.A5703T and mt.T5727G heteroplasmic sites in the ND2 gene, in different tissues and at different development stages in chickens. We also investigated whether nutritional factors could affect the mt.A5703T and mt.T5727G heteroplasmy. Our results showed that chicken mt.A5703T and mt.T5727G heteroplasmy had clear spatio-temporal specificities, which varied between tissues/development stages. The mtDNA heterogeneity was relatively stable upon nutrition intervention, 30% dietary energy restriction (from 18 to 48 days old) and different types of dietary fats (at 5% concentration, from 1 to 42 days old) did not change the breast muscle heteroplasmy of broilers at the mt.A5703T and mt.T5727G sites. In addition, multiple potential heteroplasmic sites were detected by clone sequencing in the ND2 region, which potentially reflected abundant heteroplasmy in the chicken mitochondrial genome. These results provide an important reference for further research on heteroplasmy in chicken mitochondria.

Simulations indicate that scores of lionfish (Pterois volitans) colonized the Atlantic Ocean

The invasion of the western Atlantic Ocean by the Indo-Pacific red lionfish (Pterois volitans) has had devastating consequences for marine ecosystems. Estimating the number of colonizing lionfish can be useful in identifying the introduction pathway and can inform policy decisions aimed at preventing similar invasions. It is well-established that at least ten lionfish were initially introduced. However, that estimate has not faced probabilistic scrutiny and is based solely on the number of haplotypes in the maternally-inherited mitochondrial control region. To rigorously estimate the number of lionfish that were introduced, we used a forward-time, Wright-Fisher, population genetic model in concert with a demographic, life-history model to simulate the invasion across a range of source population sizes and colonizing population fecundities. Assuming a balanced sex ratio and no Allee effects, the simulations indicate that the Atlantic population was founded by 118 (54–514, 95% HPD) lionfish from the Indo-Pacific, the Caribbean by 84 (22–328, 95% HPD) lionfish from the Atlantic, and the Gulf of Mexico by at least 114 (no upper bound on 95% HPD) lionfish from the Caribbean. Increasing the size, and therefore diversity, of the Indo-Pacific source population and fecundity of the founding population caused the number of colonists to decrease, but with rapidly diminishing returns. When the simulation was parameterized to minimize the number of colonists (high θ and relative fecundity), 96 (48–216, 95% HPD) colonists were most likely. In a more realistic scenario with Allee effects (e.g., 50% reduction in fecundity) plaguing the colonists, the most likely number of lionfish increased to 272 (106–950, 95% HPD). These results, in combination with other published data, support the hypothesis that lionfish were introduced to the Atlantic via the aquarium trade, rather than shipping. When building the model employed here, we made assumptions that minimize the number of colonists, such as the lionfish being introduced in a single event. While we conservatively modelled the introduction pathway as a single release of lionfish in one location, it is more likely that a combination of smaller and larger releases from a variety of aquarium trade stakeholders occurred near Miami, Florida, which could have led to even larger numbers of colonists than simulated here. Efforts to prevent future invasions via the aquarium trade should focus on the education of stakeholders and the prohibition of release, with adequate rewards for compliance and penalties for violations.

Genetic Diversity and Structure Analysis of Percocypris pingi (Cypriniformes: Cyprinidae): Implications for Conservation and Hatchery Release in the Yalong River

Percocypris pingi is a near threatened cyprinid species, which has suffered a dramatic decline due to anthropogenic factors. As one response to this decline, hatchery release for P. pingi has been conducted in the lower reaches of the Yalong River since 2012. To understand the conservation status of this species and the potential impact of the release of hatchery-reared fish, we studied the genetic diversity and population structure of wild and hatchery populations of P. pingi. Two hatchery populations (Jinping [JPH] and Ya’an [YAH]) and two wild populations (Muli [MLW] and Woluo [WLW]) of P. pingi were analyzed based on microsatellite markers and the mitochondrial DNA control region. The results showed that P. pingi possesses moderate levels of genetic diversity, with observed heterozygosities ranging from 0.657 to 0.770 and nucleotide diversities ranging from 0.00212 to 0.00491. Our results also suggested WLW harbors considerable proportions of genetic diversity in this species and serves as a refuge for P. pingi during anthropogenic disturbance, thus playing an important role for the conservation of P. pingi populations. Microsatellite and mitochondrial markers both indicated close genetic relationships between YAH and MLW, JPH and WLW, respectively. The results to some extent reflected the geographical provenances for original broodstocks of the two hatchery populations, which provide some practical guidance for hatchery release of P. pingi. The existence of remarkable genetic divergence distributed along limited geographical range (approximately 10 kilometers) suggests the two wild populations should be regarded at least as two distinct evolutionary significant units (ESUs) and management units (MUs). Considering reduced intra-population genetic variation in hatchery population for release and significant genetic compositions of the two hatchery populations, some appropriate breeding strategies were proposed to benefit conservation of P. pingi.

Control region sequences indicate that multiple externae represent multiple infections by Sacculina carcini (Cirripedia: Rhizocephala)

The rhizocephalan barnacle, Sacculina carcini, is a common parasite of the European shore crab, Carcinus maenas, in which it causes significant detrimental physical and behavioral modifications. In the vast majority of cases, the external portion of the parasite is present in the form of a single sac-like externa; in rare cases, double or even triple externae may occur on the same individual host. Here, we use a highly variable DNA marker, the mitochondrial control region (CR), to investigate whether multiple externae in S. carcini represent infection by multiple parasites or asexual cloning developed by a single parasite individual. Sequences for multiple externae from C. maenas hosts from the Danish inlet, Limfjorden, and from the mud flates at Roscoff, France, were compared. In almost all cases, double or triple externae from an individual host yielded different haplotypes. In the few cases where identical haplotypes were identified from externae on a multiple-infected host, this always represented the most commonly found haplotype in the population. This indicates that in Sacculina carcini, the presence of multiple externae on a single host reflects infection by different individual parasites. A haplotype network of CR sequences also suggests a degree of geographical partitioning, with no shared haplotypes between the Limfjorden and Roscoff. Our data represent the first complete CR sequences for a rhizocephalan, and a unique gene order was also revealed. Although the utility of CR sequences for population-level work must be investigated further, the CR has proved a simple to use and highly variable marker for studies of S. carcini and can easily be applied to a variety of studies in this important parasite.

Molecular analysis of phylogeographic subspecies in three Ponto-Caspian sturgeon species

Sturgeons (Order Acipenseriformes) represent an extremely valuable natural resource that is now facing depletion. In the current study we evaluate if the traditional classification in subspecies of Acipenser gueldenstaedtii, Acipenser stellatus and Huso huso, endemic to Ponto-Caspian region is sustained by molecular analysis and if these represent Evolutionary Significant Units (ESUs) that should be managed separately in conservation programs. To examine the classification of taxonomic entities we sequenced a fragment of the mitochondrial control region in case of three sturgeon species that inhabit the North-western of Black Sea and migrate for reproduction in the Lower Danube. Beside these sequences, we used previously published sequences from sturgeon individuals sampled in the Black Sea, Azov Sea and Caspian Sea. We determined the genetic diversity and genetic differentiation, conducted a Population Aggregation Analysis (PAA) and inferred an intraspecific molecular phylogeny and haplotype network. The results indicated a low level of genetic differentiation between the geographically designated subspecies and did not support a significant divergence or reciprocal monophyly between them. Our results confirm previous genetic studies with smaller samples sizes, but additional analyses including nuclear markers should be conducted for proper recommendations aiming at the development of conservation programs.

Species identification through mitochondrial rRNA genetic analysis

Inter-species and intraspecific variations in mitochondrial DNA (mtDNA) were observed in a bioinformatics analysis of the mitochondrial genomic sequences of 11 animal species. Some highly conserved regions were identified in the mitochondrial 12S and 16S ribosomal RNA (rRNA) genes of these species. To test whether these sequences are universally conserved, primers were designed to target the conserved regions of these two genes and were used to amplify DNA from 21 animal tissues, including two of unknown origin. By sequencing these PCR amplicons and aligning the sequences to a database of non-redundant nucleotide sequences, it was confirmed that these amplicons aligned specifically to mtDNA sequences from the expected species of origin. This molecular technique, when combined with bioinformatics, provides a reliable method for the taxonomic classification of animal tissues.

Shallow Population Genetic Structures of Thread-sail Filefish (Stephanolepis cirrhifer) Populations from Korean Coastal Waters

Genetic diversities, population genetic structures and demographic histories of the thread-sail filefish Stephanolepis cirrhifer were investigated by nucleotide sequencing of 336 base pairs of the mitochondrial DNA (mtDNA) control region in 111 individuals collected from six populations in Korean coastal waters. A total of 70 haplotypes were defined by 58 variable nucleotide sites. The neighbor-joining tree of the 70 haplotypes was shallow and did not provide evidence of geographical associations. Expansion of S. cirrhifer populations began approximate 51,000 to 102,000 years before present, correlating with the period of sea level rise since the late Pleistocene glacial maximum. High levels of haplotype diversities (0.974±0.029 to 1.000±0.076) and nucleotide diversities (0.014 to 0.019), and low levels of genetic differentiation among populations inferred from pairwise population F_ST values (−0.007 to 0.107), support an expansion of the S. cirrhifer population. Hierarchical analysis of molecular variance (AMOVA) revealed weak but significant genetic structures among three groups (F_CT = 0.028, p<0.05), and no genetic variation within groups (0.53%; F_SC = 0.005, p = 0.23). These results may help establish appropriate fishery management strategies for stocks of S. cirrhifer and related species.

Evidence for panmixia despite barriers to gene flow in the southern African endemic, Caffrogobius caffer (Teleostei: Gobiidae)

Background

Oceanography and life-history characteristics are known to influence the genetic structure of marine species, however the relative role that these factors play in shaping phylogeographic patterns remains unresolved. The population genetic structure of the endemic, rocky shore dwelling Caffrogobius caffer was investigated across a known major oceanographic barrier, Cape Agulhas, which has previously been shown to strongly influence genetic structuring of South African rocky shore and intertidal marine organisms. Given the variable and dynamic oceanographical features of the region, we further sought to test how the pattern of gene flow between C. caffer populations is affected by the dominant Agulhas and Benguela current systems of the southern oceans.

Results

The variable 5' region of the mtDNA control region was amplified for 242 individuals from ten localities spanning the distributional range of C. caffer. Fifty-five haplotypes were recovered and in stark contrast to previous phylogeographic studies of South African marine species, C. caffer showed no significant population genetic structuring along 1300 km of coastline. The parsimony haplotype network, AMOVA and SAMOVA analyses revealed panmixia. Coalescent analyses reveal that gene flow in C. caffer is strongly asymmetrical and predominantly affected by the Agulhas Current. Notably, there was no gene flow between the east coast and all other populations, although all other analyses detect no significant population structure, suggesting a recent divergence. The mismatch distribution suggests that C. caffer underwent a population expansion at least 14 500 years ago.

Conclusion

We propose several possible life-history adaptations that could have enabled C. caffer to maintain gene flow across its distributional range, including a long pelagic larval stage. We have shown that life-history characteristics can be an important contributing factor to the phylogeography of marine species and that the effects of oceanography do not necessarily suppress its influence on effective dispersal.