Genome survey of Misgurnus anguillicaudatus to identify genomic information, simple sequence repeat (SSR) markers, and mitochondrial genome

BACKGROUND

The dojo loach Misgurnus anguillicaudatus is an important economic species in Asia because of its nutritional value and broad environmental adaptability. Despite its economic importance, genomic data for M. anguillicaudatus is currently unavailable.

METHODS AND RESULTS

In the present study, we conducted a genome survey of M. anguillicaudatus using next-generation sequencing technology. Its genome size was estimated to be 1105.97 Mb by using K-mer analysis, and its heterozygosity ratio, repeat sequence content, GC content were 1.45%, 58.98%, and 38.03%, respectively. A total of 376,357 microsatellite motifs were identified, and mononucleotides, with a frequency of 42.57%, were the most frequently repeated motifs, followed by 40.83% dinucleotide, 7.49% trinucleotide, 8.09% tetranucleotide, and 0.91% pentanucleotide motifs. The AC/GT, AAT/ATT, and ACAG/CTGT repeats were the most abundant motifs among dinucleotide, trinucleotide, and tetranucleotide motifs, respectively. Besides, the complete mitochondrial genome was sequenced. Based on the Maximum Likelihood and Bayesian inference analyses, M. anguillicaudatus yingde in this study was the "introgressed" mitochondrial type. Seventy microsatellite loci were randomly selected from detected SSR loci to test polymorphic, of which, 20 microsatellite loci were assessed in 30 individuals from a wild population. The number of alleles (Na), observed heterozygosity (Ho), and expected heterozygosity (He) per locus ranged from 7 to 19, 0.400 to 0.933, and 0.752 to 0.938, respectively. All 20 loci were highly informative (PIC > 0.700). Eight loci deviated from Hardy-Weinberg equilibrium after Bonferroni correction (P < 0.05).

CONCLUSIONS

This is the first report of genome survey sequencing in M. anguillicaudatus, genome information, mitochondrial genome, and microsatellite markers will be valuable for further studies on population genetic analysis, natural resource conservation, and molecular marker-assisted selective breeding.

Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae

BACKGROUND

Loaches of Cobitinae, widely distributed in Eurasian continent, have high economic, ornamental and scientific value. However, the phylogeny of Cobitinae fishes within genera or family level remains complex and controversial. Up to now, about 60 Cobitinae mitogenomes had been deposited in GenBank, but their integrated characteristics were not elaborated.

RESULTS

In this study, we sequenced and analyzed the complete mitogenomes of a female Cobits macrostigma. Then we conducted a comparative mitogenome analysis and revealed the conserved and unique characteristics of 58 Cobitinae mitogenomes, including C. macrostigma. Cobitinae mitogenomes display highly conserved tRNA secondary structure, overlaps and non-coding intergenic spacers. In addition, distinct base compositions were observed among different genus and significantly negative linear correlation between AT% and AT-skew were found among Cobitinae, genus Cobitis and Pangio mitogenomes, respectively. A specific 3 bp insertion (GCA) in the atp8-atp6 overlap was identified as a unique feature of loaches, compared to other Cypriniformes fish. Additionally, all protein coding genes underwent a strong purifying selection. Phylogenetic analysis strongly supported the paraphyly of Cobitis and polyphyly of Misgurnus. The strict molecular clock predicted that Cobitinae might have split into northern and southern lineages in the late Eocene (42.11 Ma), furthermore, mtDNA introgression might occur (14.40 Ma) between ancestral species of Cobitis and ancestral species of Misgurnus.

CONCLUSIONS

The current study represents the first comparative mitogenomic and phylogenetic analyses within Cobitinae and provides new insights into the mitogenome features and evolution of fishes belonging to the cobitinae family.

Effect on amino acid and mineral content of the loach (Misgurnus anguillicaudatus) by adding Fe (II) chelating hairtail protein hydrolysates (Fe (II)-HPH) to the feed

To study the effect on amino acid and mineral content of the loach meat by adding Fe (II) chelating hairtail protein hydrolysates (Fe (II)-HPH) to the feed. A total of 100 healthy loaches were selected. After 1 week's adaptive feeding, they were randomly divided into five groups and fed with feeds containing of Fe (II)-HPH (0, 0.5, 1, 2, and 4 g/kg). On the 40th day, detection work of general nutrients (moisture, ash, crude protein, and crude fat), mineral elements (Fe, Mn, Cu, Zn, Na, K, and Ca), amino acid and amino acid score (AAS), Chemical Score (CS) and essential amino acid index (EAAI) indexes were done. The results show that crude protein has the highest content while crude fat has the lowest when amount of added Fe (II)-HPH in feed is 2 g/kg. The Fe content is significantly improved while amount of added is 1, 2, 4 g/kg. The Ca content is significantly improved and the Zn content is significantly improved while amount of added was 2 g/kg. Mn contents are significantly lower than control while amount of added is 4 g/kg. Based on analysis of amino acids in each group, the nutritional value of loach meat with 2 g/kg Fe (II)-HPH addition amount is relatively high, total amount of essential amino acids increases significantly, and EAA/TAA and EAA/NEAA improve significantly. In conclusion, adding 2 g/kg Fe (II)-HPH to feed could improve the nutritional values of loach meat.

Molecular and morphological analyses revealed a cryptic species of dojo loach Misgurnus anguillicaudatus (Cypriniformes: Cobitidae) in Japan

Although it has been reported that populations of the Japanese dojo loach Misgurnus anguillicaudatus (Cypriniformes: Cobitidae) belong to two distinct mitochondrial (mt)DNA (Type I and Type II), the taxonomic status of the species remains unresolved. To address this question, nuclear DNA and morphological analyses were performed on M. anguillicaudatus population in the Nakaikemi Wetland, where Type I and Type II lineages are sympatric. Results suggest the existence of a cryptic species (Type I) within the Japanese dojo loach.

Hepatic transcriptome analysis and identification of differentially expressed genes response to dietary oxidized fish oil in loach Misgurnus anguillicaudatus

RNA sequencing and short-read assembly were utilized to produce a transcriptome of livers from loaches (Misgurnus anguillicaudatus) fed with three different diets respectively containing fresh fish oil (FO group), medium oxidized fish oil (MO group) and high oxidized fish oil (HO group). A total of 60,663 unigenes were obtained in this study, with mean length 848.74 bp. 50,814, 49,584 and 49,814 unigenes were respectively obtained from FO, MO and HO groups. There were 2,343 differentially expressed genes between FO and MO, with 855 down- and 1,488 up-regulated genes in the MO group. 2,813 genes were differentially expressed between FO and HO, including 1,256 down- and 1,552 up-regulated genes in the HO group. 2,075 differentially expressed genes were found in the comparison of MO and HO, including 1,074 up- and 1,001 down-regulated genes in the MO group. Some differentially expressed genes, such as fatty acid transport protein (fatp), fatty acid binding protein (fabp), apolipoprotein (apo), peroxisome proliferator activated receptor-gamma (ppar-γ), acetyl-CoA synthetase (acs) and arachidonate 5-lipoxygenase (alox5), were involved in lipid metabolism, suggesting these genes in the loach were responsive to dietary oxidized fish oil. Results of transcriptome profilings here were validated using quantitative real time PCR in fourteen randomly selected unigenes. The present study provides insights into hepatic transcriptome profile of the loach, which is a valuable resource for studies of loach genomics. More importantly, this study identifies some important genes responsible for dietary oxidized fish oil, which will benefit researches of lipid metabolism in fish.

Molecular Evidence for Multiple Origins of the European Spined Loaches (Teleostei, Cobitidae)

We present a phylogenetic investigation of the Northern Clade, the major monophyletic clade within the freshwater fish family Cobitidae, one of the most prominent families of freshwater fishes found in Asian and European waters. Phylogenetic reconstructions based on the cytochrome b and RAG-1 genes show the genera Microcobitis, Sabanejewia, Koreocobitis and Kichulchoia as monophyletic groups. These reconstructions also show a Cobitis sensu lato and a Misgurnus sensu lato group. The Cobitis sensu lato group includes all species of Cobitis, Iksookimia, Niwaella and Kichulchoia, while the Misgurnus sensu lato group includes Misgurnus, Paramisgurnus and Koreocobitis. Although the monophyly of both the Cobitis sensu lato and Misgurnus sensu lato groups is supported, relationships within the groups are incongruent with current generic definitions. The absence of monophyly of most genera included in the Cobitis sensu lato group (Cobitis, Iksookimia and Niwaella) or their low genetic differentiation (Kichuchoia) supports their consideration as synonyms of Cobitis. Molecular phylogenies indicate that the Asian species of Misgurnus experienced a mitochondrial introgression from a lineage of Cobitis. We also find two nuclear haplotypes in the same Cobitis species from the Adriatic area that, in the absence of morphological differentiation, may indicate molecular introgression. Most lineages within the Northern Clade consist of species found in East Asia. However, some lineages also contain species from Europe and Asia Minor. The phylogenetic relationships presented here are consistent with previous studies suggesting an East Asian origin of the Northern Clade. According to the current distributions and phylogenetic relationships of the Misgurnus sensu lato and Cobitis clade lineages, particularly of M. fossilis and C. melanoleuca, the range expansion of East Asian species into Europe was most likely via Siberia into Northern and Central Europe. Phylogenetic analyses also show that the Cobitis sensu lato group consists of two clear subgroups (I and II), each presenting geographical differences. Subgroup I is distributed exclusively in East Asian drainages with an Eastern European offshoot (C. melanoleuca), whereas Subgroup II includes species widespread throughout Europe (including the Mediterranean), Asia Minor, the Black Sea and the Caucasus, with some lineages related to species restricted to East Asia.

Comparative analysis of mitochondrial genomes in distinct nuclear ploidy loach Misgurnus anguillicaudatus and its implications for polyploidy evolution

Misgurnus anguillicaudatus has several natural ploidy types. To investigate whether nuclear polyploidy have an impact on mitochondrial DNA (mtDNA), the complete mitochondrial genomes (mitogenomes) of five distinct ploidy M. anguillicaudatus (natural diploid, triploid, tetraploid, pentaploid and hexaploid), which were collected in central China, were sequenced and analyzed. The five mitogenomes share the same gene arrangement and have similar gene size, base composition and codon usage pattern. The most variable regions of the mitogenome were the protein-coding genes, especially the ND4L (5.39% mutation rate). Most variations occurred in tetraploids. The phylogenetic tree showed that the tetraploid M. anguillicaudatus separated early from other ploidy loaches. Meanwhile, the mitogenomes from pentaploids, and hexaploids have the closest phylogenetic relations, but far from that of tetraploids, implying that pentaploids and hexaploids could not be formed from tetraploids, possibly from the diploids and triploids. The genus Misgurnus species were divided into two divergent inter-genus clades, and the five ploidy M. anguillicaudatus were monophyletic, which support the hypotheses about the mitochondrial introgression in loach species.