Phylogeny and Evolution of Multiple Common Carp (Cyprinus carpio L.) Populations Clarified by Phylogenetic Analysis Based on Complete Mitochondrial Genomes

Genetic diversity of the common carp black strain population based on mtDNA (D-loop and cytb)

The common strain black carp (Cyprinus carpio var. baisenensis) is a culturally important carp strain that is raised and cultured in Guangxi Province, China. Its color reflects the interactions between the Burau people and their surrounding environment. The population of the common carp black strain was isolated and cultured in a rice-fish integration system. To explore the genetic diversity and protection of germplasm resources, we analyzed mitochondrial DNA (mtDNA) sequences, specifically the displacement loop (D-loop) and cytochrome b (Cytb), using single-nucleotide polymorphisms (SNP). We compared these sequences with those from four other local common carp populations. The study included a total of 136 adult common carps from five strain populations: the common black carp strain (HJ = 31), Jian (F = 30), Heilongjiang (H = 10), Songpu (S = 31), and Saijiang (SJ = 34). The results of the Cytb and D-loop analyses showed that the Heilongjiang carp (H) and Saijiang (SJ) populations had the highest levels of haplotype diversity (0.867 ± 0.034785) and nucleotide diversity (π = 0.0063 ± 0.000137 and 0.0093 ± 0.000411), respectively. On the other hand, the Common carp black strain population (HJ) exhibited the lowest haplotype diversity in both Cytb and D-loop, with haplotype 2 being the most commonly observed among the populations. Private haplotypes dominated the five common carp populations, which were significantly different at P<0.001. Furthermore, analyzing the coefficient of genetic differentiation (Fst), the highest genetic difference was observed between Saijiang (SJ) and Heilongjiang (H) (Fst = 0.963), whereas the lowest was observed between Songpu (S) and the Common carp black strain population (HJ) (Fst = 0.019) for the Cytb gene sequences. For the D-loop, the Common carp black strain population (HJ) and Songpu (S) (Fst = 0.7) had the highest values, and Heilongjiang (H) and Common black carp strain (HJ) had an Fst of 0.125. Additionally, the AMOVA analysis revealed a higher level of variance for the Cytb and D-loop genes, indicating lower genetic diversity within the local carp community. On the other hand, the phylogenetic tree analysis showed that the five carp populations were closely related and formed a distinct cluster. The distinct cluster of populations suggests a common ancestor or recent gene flow, possibly due to geographic proximity or migration history, and unique genetic characteristics, possibly due to adaptations or selective pressures. The results of this study provide valuable insights into the genetic diversity of the common strain black carp, which can have implications for conservation, breeding programs, evolutionary studies, and fisheries management.

Genetic sources and diversity of the paddy field carp in the Pearl River basin inferred from two mitochondrial loci

Paddy field carp (PF-carp) is an economically important fish cocultured with rice in traditional agricultural systems. Several distinctive strains of PF-carp have been formed through years of artificial and cross breeding. There is a concern about the status of germplasm resources among PF-carp, since little is known about the genetic sources, diversity, or differentiation. In this study we collected 17 PF-carp populations covering Daotian carp (DTL), Ru Yuan No. 1 (RY), Jinbian carp (JBL), Shaijiang carp (SJL), and Wu carp (WL) along the Pearl River basin to explore their genetic sources and diversity using concatenated sequences of the mitochondrial cytochrome b gene and the D-Loop region. According to the haplotype network analyses, 1, 9, and 57 haplotypes originated from Cyprinus carpio carpio, Cyprinus carpio haematopterus and Cyprinus carpio rubrofuscus, respectively, confirming that genetic introgression has occurred in Pearl River PF-carp populations and Cyprinus carpio carpio was the most common species for genetic origin. The results showed that RY exhibited the lowest level of nucleotide diversity (π = 0.0011) due to high-intensity breeding and was significantly differentiated from the other four strains. PF-carp strains in these remote traditional systems tended to experience artificial selection and a lack of farmer connection that gradually increased genetic differentiation among strains. Notably, three populations of JBL exhibited significant high-level differentiation, since they originated from mountainous areas hindering farmers from fry exchange. In contrast, no significant differentiation was uncovered in the WL populations, since this strain is the most popular cultured strain and has undergone artificial exchange of parents and fry in many cultured regions. This study helps us to understand the status of germplasm resources among PF-carp and to trace their genetic origin before being introduced for local cultivation.

The complete mitochondrial genome of the common strain black carp (Cyprinus carpio var. baisenensis)

The common strain black carp (Cyprinus carpio var. baisenensis), known for its black skin is cultured in the integrated rice-agriculture system and non-escape property under torrential floods. The total mitogenome length of Cyprinus carpio var. baisenensis obtained in this study was 16,478bp, consisting of 13 protein-coding genes, 22 tRNA genes, 2rRNA genes (large and small), a light strand origin of replication, and one major non-coding region. By providing the complete mitochondrial genomes of Cyprinus carpio var. baisenensis, we will further understand the phylogenetic relationships within genus.

Genetic analysis of whole mitochondrial genome of Lateolabrax maculatus (Perciformes: Moronidae) indicates the presence of two populations along the Chinese coast

The whole mitochondrial genome ofLateolabrax maculatus(Cuvier, 1828) was used to investigate the reasons for the observed patterns of genetic differentiation among 12 populations in northern and southern China. The haplotype diversity and nucleotide diversity ofL. maculatuswere 0.998 and 0.00169, respectively. Pairwise FSTvalues between populations ranged from 0.001 to 0.429, correlating positively with geographic distance. Genetic structure analysis and haplotype network analysis indicated that these populations were split into two groups, in agreement with geographic segregation and environment. Tajima’s D values, Fu’s Fs tests and Bayesian skyline plot (BSP) indicated that a demographic expansion event may have occurred in the history ofL. maculatus. Through selection pressure analysis, we found evidence of significant negative selection at the ATP6, ND3, Cytb, COX3, COX2 and COX1 genes. In our hypotheses, this study implied that demographic events and selection of local environmental conditions, including temperature, are responsible for population divergence. These findings are a step forward toward the understanding of the genetic basis of differentiation and adaptation, as well as conservation ofL. maculatus.

Morphological and molecular studies of Eudiplozoon nipponicum (Goto, 1891) and Eudiplozoon kamegaii sp. n. (Monogenea; Diplozoidae)

Eudiplozoon nipponicum (Goto, 1891) Khotenovsky, 1985 (Monogenea: Diplozoidae), is known to parasitise Cyprinus carpio Linnaeus and species of Carassius. In this study, we conducted a taxonomic re-examination of E. nipponicum using genetic analysis and morphological comparisons from different host species from a single water system. rDNA nucleotide sequences of the internal transcription spacer 2 (ITS-2) region (645 bp) showed interspecific-level genetic differences among diplozoids from species of Carassius and C. carpio (p-distance: 3.1-4.0%) but no difference among those from different species of Carassius (0-0.4%) or between those from C. carpio collected in Asia and Europe (0-1.1%). Large variation was observed among 346 bp cytochrome c oxidase subunit I (COI) sequences (0.3-16.0 %); the topology of the phylogenetic tree showed no relationship to host genera or geographical regions of origin. Morphological observation showed that average clamp size of diplozoids from C. carpio was larger than those from Carassius spp. The number of folds on the hindbody was 10-25 for diplozoids from C. carpio and 12-19 for those from Carassius spp. Thus, our ITS-2 sequence and morphological comparison results indicate that diplozoids from C. carpio and species of Carassius belong to different species. The scientific name E. nipponicum should be applied to the species infected to the type host, Carassius sp. of Nakabo (2013) (Japanese name ginbuna). The diplozoid infecting C. carpio (Eurasian type) should be established as a new species: Eudiplozoon kamegaii sp. n. A neotype of E. nipponicum is designated in this report because the original E. nipponicum specimens are thought to have been lost.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

Patterns of Geographical and Potential Adaptive Divergence in the Genome of the Common Carp (Cyprinus carpio)

The common carp, Cyprinus carpio, is a cyprinid fish species cultured in Europe and Asia. It accounts for >70% of freshwater aquaculture production worldwide. We conducted a population genomics analysis on C. carpio using high-throughput SNP genotyping of 2,198 individuals from 14 populations worldwide to determine the genetic architecture of common carp populations and the genetic bases for environmental adaptation. Structure analyses including phylogeny and principal component analysis were also conducted, showing distinct geographical patterns in European and Asian populations. The linkage disequilibrium block average lengths of the 14 populations ranged from 3.94 kb to 36.67 kb. Genes within selective sweep regions were identified by genome scanning among the different populations, including gdf6a, bmpr1b, and opsin5. Gene Ontology and KEGG enrichment analyses revealed potential trait-related loci and genes associated with body shape, scaling patterns, and skin color. This population genomics analysis may provide valuable clues for future genome-assisted breeding of C. carpio.

Tour around the globe: The case of invasive tapeworm Atractolytocestus huronensis (Cestoda: Caryophyllidea), a parasite of common carp

The monozoic tapeworm Atractolytocestus huronensis Anthony, 1958 (Cestoda: Caryophyllidea), an intestinal parasite of the common carp, is characterized by its invasive character and potential to colonize new territories. It was initially described from North America and has also been found in several European countries. The most recent findings of A. huronensis originated from China and South Africa; however, no data on genetic relationships of these populations were available. The current study provides the first molecular characterisation of A. huronensis from South Africa and China using a partial sequence of mitochondrial cytochrome c oxidase subunit 1 (cox1) and a complete ribosomal ITS2 spacer. Ribosomal and mitochondrial data were applied for phylogenetic analyses in order to assess the genetic interrelationships among global A. huronensis populations. Divergent intragenomic copies of ribosomal ITS2 were detected in all analysed specimens; the structure and frequency of the ITS2 variants of tapeworms from China and South Africa corresponded with the data on ITS2 paralogues observed previously in A. huronensis from Slovakia, the United States and the United Kingdom. The phylogenetic analysis of cox1 indicated that A. huronensis exist in two slightly differentiated clusters; one cluster was supported by all phylogenetic approaches (NJ, ML, BI) and was represented by samples from China, the USA and the UK. A second cluster was represented by tapeworms from continental Europe (Slovakia, Hungary, Romania, Croatia) and South Africa. Haplotype network analysis revealed that the highest population diversity occurs in China. The results provide useful pilot information about the interrelationships of A. huronensis on four continents and indicate that China, or the eastern Palaearctic, served as the original source population for the global expansion of this invasive tapeworm. Data on the origin and distribution of the common carp, the only specific host of A. huronensis, are also discussed.

Two mitochondrial genomes in Alcedinidae (Ceryle rudis/Halcyon pileata) and the phylogenetic placement of Coraciiformes

Coraciiformes comprises 209 species belonging to ten families with significant divergence on external morphologies and life styles. The phylogenetic placement of Coraciiformes was still in debate. Here, we determined the complete mitochondrial genomes (mitogenomes) of Crested Kingfisher (Ceryle rudis) and Black-capped Kingfisher (Halcyon pileata). The mitogenomes were 17,355 bp (C. rudis) and 17,612 bp (H. pileata) in length, and both of them contained 37 genes (two rRNA genes, 22 tRNA genes and 13 protein-coding genes) and one control region. The gene organizations and characters of two mitogenomes were similar with those of other mitogenomes in Coraciiformes, however the sizes and nucleotide composition of control regions in different mitogenomes were significantly different. Phylogenetic trees were constructed with both Bayesian and Maximum Likelihood methods based on mitogenome sequences from 11 families of six orders. The trees based on two different data sets supported the basal position of Psittacidae (Psittaciformes), the closest relationship between Cuculiformes (Cuculidae) and Trogoniformes (Trogonidae), and the close relationship between Coraciiformes and Piciformes. The phylogenetic placement of the clade including Cuculiformes and Trogoniformes has not been resolved in present study, which need further investigations with more molecular markers and species. The mitogenome sequences presented here provided valuable data for further taxonomic studies on Coraciiformes and other related groups.

Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China's coast, which are shaped by unknown mechanisms in a north-south divide.