Polymeric immunoglobulin receptor in dojo loach (Misgurnus anguillicaudatus): Molecular characterization and expression analysis in response to bacterial and parasitic challenge

The polymeric immunoglobulin receptor (pIgR) is an essential component of the mucosal immune system in jawed vertebrates including teleost fish, which mediate transepithelial transport of secretory immunoglobulins (sIgs) to protect organisms against environmental pathogens. In this study, we firstly cloned and identified the pIgR from dojo loach (Misgurnus anguillicaudatus). The full-length cDNA of Ma-pIgR was of 1145 bp, containing an open reading frame (ORF) of 1101 bp encoded a predicted protein of 336 amino acids. The structure of Ma-pIgR is comprised of a signal peptide, a transmembrane region, an intracellular region and an extracellular region with two Ig-like domains (ILDs), which are similar to their counterparts described in other teleosts. Multiple sequence alignment and phylogenetic analysis showed the dojo loach is closely related to the fish family Cyprinidae. The transcriptional level of Ma-pIgR was detected by quantitative real-time PCR (qRT-PCR) in different tissues and high expression was found in liver, skin, kidney, eye, fin and gills. Two infection models of the loach with bacteria (Aeromonas hydrophila) and parasite (Ichthyophthirius multifiliis) were constructed for the first time. Histological studies showed the goblet cells in skin significantly increased and the ratio of gill length to width also significantly changed after challenged with A.hydrophila. Both challenge experiments resulted in the significant up-regulated expression of Ma-pIgR not only in kidney and spleen, but also in skin and gills. Our results suggest that pIgR may play an important role in skin and gill mucosal immunity to protect the loach against bacteria and parasite.

Fertility and ploidy of gametes of allodiploid and allotriploid loaches produced by diploid Misgurnus anguillicaudatus females and Paramisgurnus dabryanus males

Artificial and natural hybridization of dojo loach Misgurnus anguillicaudatus (2N = 50, DD for short) and large-scale loach Paramisgurnus dabryanus (2N = 48, PP for short) are well-grown. However, these hybrid loaches have not yet been examined for fertility and ploidy of gametes. Here, histological observations, artificial propagation, observations of embryonic development, larval morphology, and ploidy analyses were conducted to determine the fertility and ploidy of gametes of allodiploid (DP for short) and allotriploid (DDP for short) loaches, produced by DD females × PP males and induced from fertilized eggs of DD females × PP males by cold shock to prevent the second polar body release, respectively. The ovaries of DP and DDP included smaller number of eggs when compared with those of the control DD, while full-grown oocytes were observed. Testes of these two loaches were delayed-developed without spermatids or mature spermatozoa. Results obtained here showed that DP and DDP were fertility-weakened female and sterile male. Moreover, DP females and DDP females could, respectively, produce few viable haploid eggs and few viable haploid and diploid eggs. This study will provide valuable information for fish hybrid researches.

Growth performance, fatty-acid composition, lipid deposition and hepatic-lipid metabolism-related gene expression in juvenile pond loach Misgurnus anguillicaudatus fed diets with different dietary soybean oil levels

A 10 week feeding trial was conducted to evaluate the effects of different dietary soybean oil (SO) levels on growth performance, fatty-acid composition and lipid deposition in viscera, histology and histochemistry of liver, intestine and hepatic-lipid metabolism-related gene expressions in pond loach Misgurnus anguillicaudatus juveniles. Misgurnus anguillicaudatus (mean ± s.d. mass 0·40 ± 0·01 g) were fed five experimental diets containing SO at different concentrations: 0, 20, 32, 56 and 100% SO and a diet containing 100% fish oil (100% FO). The mass gains and specific growth rates of M. anguillicaudatus fed 20% SO and 100% FO diets were significantly higher than those of the other groups (P < 0·05). The lipid content of viscera and the amount of cytoplasmic vacuolation in the liver increased with incremental dietary SO level. Meanwhile, increasing dietary SO levels up-regulated the messenger (m)RNA levels of lipogenic genes (such as Δ6fad, scd, pparγ, fas and srebp-1) and down-regulated the mRNA levels of the lipolytic genes (such as pparα, cpt1, atgl and hsl) in the liver. The percentage of 20:4n-6 significantly (P < 0·05) increased with increasing dietary SO level, which might be correlated with the up-regulation of the mRNA level of Δ6fad. The highest levels of dietary SO, however, had a negative effect on growth performance, lipid deposition of viscera and histology and histochemstry of liver and intestine. The increased lipid accumulation induced by incremental dietary SO level probably occurred through different strategies for lipid metabolism as a result of competition between lipolysis and lipogenesis and between export and import of lipids in this species.

Molecular characterization and expression analysis of complement component 3 in dojo loach (Misgurnus anguillicaudatus)

The complement component 3 (C3) is a central component of complement system. All three pathways converge at formation of C3 convertases and share the terminal pathways of membrane attack complex (MAC) formation. In this study, three isoforms of C3 were discovered in Misgurnus anguillicaudatus, named "C3-1", "C3-2" and "C3-3", respectively. The full-length of C3-1 cDNA sequence was firstly identified and analyzed from dojo loach (Misgurnus anguillicaudatus). The Ma-C3-1 cDNA sequence comprised of 4509 bp encoding 1454 amino acids with a putative signal peptide of 20 amino acid residues. The deduced amino acid sequence showed that Ma-C3-1 has conserved residues and domain, which are known to be crucial for C3 function. Interestingly, an amino acid substitution of the highly conserved GCGEQ was discovered in Ma-C3-1. Phylogenetic analysis showed that Ma-C3-1 was closely related to Cyprinidae. The mRNA expression levels of three isoforms of C3 were detected in kidney, eye, spleen, gonad, heart, fin ray, gut, muscle, brain, gill, skin, blood and liver. The expression of Ma-C3-1 and Ma-C3-3 were mainly detected in liver, followed by spleen, gonad. However, the high expression of Ma-C3-2 was found in kidney, followed by blood and gonad. The morphological changes of gill and skin, and the expression pattern of these three isoforms C3 molecular following the infection with Aeromonas hydrophila were investigated. The mRNA expression levels of three C3 isoforms were up-regulated in the gill, skin, liver and spleen after infection with A.hydrophila. Similarly, challenge experiments resulted in significant up-regulated expression of other complement-relevant genes in gill, liver and skin, such as C4, C5, C8b, especially at 24 h and 36 h. These results suggest that complement system might play an important role not only in liver, but also in the mucosal tissues as gill and skin of teleost fish.

Hybrid asexuality as a primary postzygotic barrier between nascent species: On the interconnection between asexuality, hybridization and speciation

Although sexual reproduction is ubiquitous throughout nature, the molecular machinery behind it has been repeatedly disrupted during evolution, leading to the emergence of asexual lineages in all eukaryotic phyla. Despite intensive research, little is known about what causes the switch from sexual reproduction to asexuality. Interspecific hybridization is one of the candidate explanations, but the reasons for the apparent association between hybridization and asexuality remain unclear. In this study, we combined cross-breeding experiments with population genetic and phylogenomic approaches to reveal the history of speciation and asexuality evolution in European spined loaches (Cobitis). Contemporary species readily hybridize in hybrid zones, but produce infertile males and fertile but clonally reproducing females that cannot mediate introgressions. However, our analysis of exome data indicates that intensive gene flow between species has occurred in the past. Crossings among species with various genetic distances showed that, while distantly related species produced asexual females and sterile males, closely related species produce sexually reproducing hybrids of both sexes. Our results suggest that hybridization leads to sexual hybrids at the initial stages of speciation, but as the species diverge further, the gradual accumulation of reproductive incompatibilities between species could distort their gametogenesis towards asexuality. Interestingly, comparative analysis of published data revealed that hybrid asexuality generally evolves at lower genetic divergences than hybrid sterility or inviability. Given that hybrid asexuality effectively restricts gene flow, it may establish a primary reproductive barrier earlier during diversification than other "classical" forms of postzygotic incompatibilities. Hybrid asexuality may thus indirectly contribute to the speciation process.

Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus

Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC₃, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.

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.

Molecular cloning and mRNA expression pattern of Sox4 in Paramisgurnus dabryanus

Sox4 belonged to the SoxC subfamily of the Sox family, which play important roles in the development of the vertebrate gonad and nervous system. A Sox4 homologue was cloned from brain of Paramisgurnus dabryanus by using homologous cloning and rapid amplification of cDNA ends (RACE), designated as PdSox4. The full-length cDNA was 2163bp, containing the 759bp 5'-untranslated region, 267bp 3'-untranslated region and encoding a putative protein of 378 amino acids with a characteristic high mobility group box (HMG-box) DNA-binding domain of 79 amino acids with the specific motif (RPMNAFMVW). Alignment and phylogenetic analyses indicated that PdSox4 shares highly identical sequence with Sox4 homologues from different species. The signal peptide analysis predicted that PdSox4 is a non-secretory protein. The hydropathy profile of PdSox4 protein revealed that this protein is hydrophilic in nature. The expression profiles of PdSox4 in different developmental stages and various adult tissues of sexs were analyzed by quantitative real-time RT-PCR (qRT-PCR) and In situ hybridization (ISH). The results showed that PdSox4 was ubiquitously expressed during embryogenesis and various adult tissues, especially in central nervous system. Tissue distribution analyses revealed that PdSox4 was expression in developing germ cells. Taken together, these preliminary findings suggested that PdSox4 is highly conserved during vertebrate evolution and involved in a wide range of developmental processes including embryogenesis, neurogenesis and gonad development.

Genome-wide identification of genes probably relevant to the adaptation of schizothoracins (Teleostei: Cypriniformes) to the uplift of the Qinghai-Tibet Plateau

BACKGROUND

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated.

RESULTS

The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation.

CONCLUSIONS

By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.

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.

Gene expression profile changes induced by acute toxicity of [C16 mim]Cl in loach (Paramisgurnus dabryanus)

Ionic liquids (ILs) are widely used as reaction media in various commercial applications. Many reports have indicated that most ILs are poorly decomposed by microorganisms and are toxic to aquatic organisms. In this study, differential gene expression profiling was conducted using a suppression subtraction hybridization cDNA library from hepatic tissue of the loach (Paramisgurnus dabryanus) after exposure to 1-hexadecyl-3-methylimidazolium chloride ([C₁₆ mim]Cl), a representative IL. Two hundred and fifty-nine differentially expressed candidate genes, whose expression was altered by >2.0-fold by the [C₁₆ mim]Cl treatment, were identified, including 127 upregulated genes and 132 downregulated genes. A gene ontology analysis of the known genes isolated in this study showed that [C₁₆ mim]Cl-responsive genes were involved in cell cycle, stimulus response, defense response, DNA damage response, oxidative stress responses, and other biological responses. To identify candidate genes that may be involved in [C₁₆ mim]Cl-induced toxicity, 259 clones were examined by Southern blot macroarray hybridization, and 20 genes were further characterized using quantitative real-time polymerase chain reaction. Finally, six candidate genes were selected, including three DNA damage response genes, two toxic substance metabolic genes, and one stress protein gene. Our results indicate that these changes in gene expression are associated with [C₁₆ mim]Cl-induced toxicity, and that these six candidate genes can be promising biomarkers for detecting [C₁₆ mim]Cl-induced toxicity. Therefore, this study demonstrates the use of a powerful assay to identify genes potentially involved in [C₁₆ mim]Cl toxicity, and it provides a foundation for the further study of related genes and the molecular mechanism of [C₁₆ mim]Cl toxicity. © 2016 Wiley Periodicals, Inc. Environ Toxicol 32: 404-416, 2017.

A Time-Calibrated Mitogenome Phylogeny of Catfish (Teleostei: Siluriformes)

A very significant part of the world's freshwater ichthyofauna is represented by ancient, exceptionally diverse and cosmopolitan ray-finned teleosts of the order Siluriformes. Over the years, catfish have been established as an exemplary model for probing historical biogeography at various scales. Yet, several tantalizing gaps still exist in their phylogenetic history, timeline and mode of diversification. Here, we re-examine the phylogeny of catfish by assembling and analyzing almost all publicly available mitogenome data. We constructed an ingroup matrix of 62 full-length mitogenome sequences from 20 catfish families together with four cypriniform outgroups, spanning 15,557 positions in total. Partitioned maximum likelihood analyses and Bayesian relaxed clock dating using fossil age constraints provide some useful and novel insights into the evolutionary history of this group. Loricarioidei are recovered as the first siluriform group to diversify, rendering Neotropics the cradle of the order. The next deepest clade is the South American Diplomystoidei placed as a sister group to all the remaining Siluroidei. The two multifamilial clades of "Big Asia" and "Big Africa" are also recovered, albeit nodal support for the latter is poor. Within "Big Asia", Bagridae are clearly polyphyletic. Other interfamilial relationships, including Clariidae + Heteropneustidae, Doradidae + Auchenipteridae and Ictaluridae + Cranoglanididae are robustly resolved. Our chronogram shows that siluriforms have a Pangaean origin, at least as far back as the Early Cretaceous. The inferred timeline of the basal splits corroborates the "Out-of-South America" hypothesis and accords well with the fossil record. The divergence of Siluroidei most likely postdated the final separation of Africa and South America. An appealing case of phylogenetic affinity elaborated by biogeographic dispersal is exemplified by the Early Paleogene split between the Southeast Asian Cranoglanididae and Ictaluridae, with the latter radiating into North America's freshwater realm by Eocene. The end of Cretaceous probably concludes the major bout of diversification at the family level while with the dawn of the Cenozoic a prolific radiation is evident at the generic level.

Comprehensive Transcriptome Analysis Provides Evidence of Local Thermal Adaptation in Three Loaches (Genus: Misgurnus)

The geographic distribution of three Misgurnus species, M. anguillicaudatus, M. bipartitus, and M. mohoity, displays a specific pattern in China, coincident with temperature zones. In this study, we sequenced the transcriptomes of these three species and used the sequences to investigate the lineage-specific adaptations within the genus Misgurnus. In total, 51 orphan genes (19 in M. anguillicaudatus, 18 in M. bipartitus, and 14 in M. mohoity) that may contribute to the species-specific adaptations were identified. An analysis of 1392 one-to-one orthologous genes revealed significantly higher ratios of nonsynonymous-to-synonymous substitutions in the M. mohoity lineage than in M. anguillicaudatus. The genes displaying signatures of positive selection and rapid evolution in Misgurnus were involved in four function categories, (1) energy metabolism; (2) signal transduction; (3) membrane; and (4) cell proliferation or apoptosis, implying that these candidate genes play critical roles in the thermal adaptation of the fish to their living environments. We also detected more than five positively selected sites in cldn15lb and isca1, which function as important factors in paracellular Na⁺ transport and Fe/S cluster assembly, respectively. Overall, our study provides valuable insights into the adaptive evolution of loaches from different temperature zones in China and is a foundation for future studies to clarify the genetic basis of temperature adaptation in fishes.

Resolving Cypriniformes relationships using an anchored enrichment approach

BACKGROUND

Cypriniformes (minnows, carps, loaches, and suckers) is the largest group of freshwater fishes in the world (~4300 described species). Despite much attention, previous attempts to elucidate relationships using molecular and morphological characters have been incongruent. In this study we present the first phylogenomic analysis using anchored hybrid enrichment for 172 taxa to represent the order (plus three out-group taxa), which is the largest dataset for the order to date (219 loci, 315,288 bp, average locus length of 1011 bp).

RESULTS

Concatenation analysis establishes a robust tree with 97 % of nodes at 100 % bootstrap support. Species tree analysis was highly congruent with the concatenation analysis with only two major differences: monophyly of Cobitoidei and placement of Danionidae.

CONCLUSIONS

Most major clades obtained in prior molecular studies were validated as monophyletic, and we provide robust resolution for the relationships among these clades for the first time. These relationships can be used as a framework for addressing a variety of evolutionary questions (e.g. phylogeography, polyploidization, diversification, trait evolution, comparative genomics) for which Cypriniformes is ideally suited.

Oxidative stress, genotoxicity and cytotoxicity of 1-methyl-3-octylimidazolium chloride on Paramisgurnus dabryanus

This study evaluated the toxicity of 1-methyl-3-octylimidazolium chloride ([C₈mim]Cl) on Paramisgurnus dabryanus by enzyme analysis, comet assay, and apoptosis analysis. The study showed that [C₈mim]Cl had an obvious toxic effect inducing oxidative stress, genotoxicity, and cytotoxicity to fish liver cells. [C₈mim]Cl also induced changes in the activities of superoxide dismutase and catalase, and the glutathione content and malondialdehyde level in fish exposed at 20-80mgL^-1. With increased exposure concentration and time, the four antioxidant enzyme activities, three different comet parameters and apoptosis rates of tested cells were significantly increased, with significant differences (P<0.05 or P<0.01) observed between control group and each treatment group. This study shows that [C₈mim]Cl could be a threat to aquatic organism health when accidentally released into aquatic ecosystems.

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.

Transcriptomic Analysis of Compromise Between Air-Breathing and Nutrient Uptake of Posterior Intestine in Loach (Misgurnus anguillicaudatus), an Air-Breathing Fish

Dojo loach (Misgurnus anguillicaudatus) is an air-breathing fish species by using its posterior intestine to breathe on water surface. So far, the molecular mechanism about accessory air-breathing in fish is seldom addressed. Five cDNA libraries were constructed here for loach posterior intestines form T01 (the initial stage group), T02 (mid-stage of normal group), T03 (end stage of normal group), T04 (mid-stage of air-breathing inhibited group), and T05 (the end stage of air-breathing inhibited group) and subjected to perform RNA-seq to compare their transcriptomic profilings. A total of 92,962 unigenes were assembled, while 37,905 (40.77 %) unigenes were successfully annotated. 2298, 1091, and 3275 differentially expressed genes (fn1, ACE, EGFR, Pxdn, SDF, HIF, VEGF, SLC2A1, SLC5A8 etc.) were observed in T04/T02, T05/T03, and T05/T04, respectively. Expression levels of many genes associated with air-breathing and nutrient uptake varied significantly between normal and intestinal air-breathing inhibited group. Intraepithelial capillaries in posterior intestines of loaches from T05 were broken, while red blood cells were enriched at the surface of intestinal epithelial lining with 241 ± 39 cells per millimeter. There were periodic acid-schiff (PAS)-positive epithelial mucous cells in posterior intestines from both normal and air-breathing inhibited groups. Results obtained here suggested an overlap of air-breathing and nutrient uptake function of posterior intestine in loach. Intestinal air-breathing inhibition in loach would influence the posterior intestine's nutrient uptake ability and endothelial capillary structure stability. This study will contribute to our understanding on the molecular regulatory mechanisms of intestinal air-breathing in loach.

A Ploidy Difference Represents an Impassable Barrier for Hybridisation in Animals. Is There an Exception among Botiid Loaches (Teleostei: Botiidae)?

One of the most efficient mechanisms to keep animal lineages separate is a difference in ploidy level (number of whole genome copies), since hybrid offspring from parents with different ploidy level are functionally sterile. In the freshwater fish family Botiidae, ploidy difference has been held responsible for the separation of its two subfamilies, the evolutionary tetraploid Botiinae and the diploid Leptobotiinae. Diploid and tetraploid species coexist in the upper Yangtze, the Pearl River and the Red River basins in China. Interestingly, the species ‘Botia’ zebra from the Pearl River basin combines a number of morphological characters that otherwise are found in the diploid genus Leptobotia with morphological characters of the tetraploid genus Sinibotia, therefore the aim of the present study is to test weather ‘B.’ zebra is the result of a hybridisation event between species from different subfamilies with different ploidy level. A closer morphological examination indeed demonstrates a high similarity of ‘B.’ zebra to two co-occurring species, the diploid Leptobotia guilinensis and the tetraploid Sinibotia pulchra. These two species thus could have been the potential parental species in case of a hybrid origin of ‘B.’ zebra. The morphologic analysis further reveals that ‘B.’ zebra bears even the diagnostic characters of the genera Leptobotia (Leptobotiinae) and Sinibotia (Botiinae). In contrast, a comparison of six allozyme loci between ‘B.’ zebra, L. guilinensis and S. pulchra showed only similarities between ‘B.’ zebra and S. pulchra, not between ‘B.’ zebra and L. guilinensis. Six specimens of ‘B.’ zebra that were cytogenetically analysed were tetraploid with 4n = 100. The composition of the karyotype (18% metacentric, 18% submetacentric, 36% subtelocentric and 28% acrocentric chromosomes) differs from those of L. guilinensis (12%, 24%, 20% and 44%) and S. pulchra (20%, 26%, 28% and 26%), and cannot be obtained by any combination of genomes from L. guilinensis and S. pulchra. Phylogenetic reconstructions based on sequence data of the mitochondrial cytochrome b gene and the nuclear RAG-1 gene invariably places ‘Botia’ zebra as sister species to S. pulchra, while L. guilinensis is only distantly related. The presented combination of genetic data demonstrates that ‘B.’ zebra is not the result of a hybridisation, but a species of tetraploid genus Sinibotia with a striking morphological evolution towards an enormous similarity with a co-occurring, but not directly related species. The complete lack of knowledge of the ecology of these species, their main predators or their ecological interactions hampers any conclusion regarding the evolutionary advantage of such adaptation.

Toxic effects of imidacloprid on adult loach (Misgurnus anguillicaudatus)

The present investigation was aimed to assess the effects of imidacloprid on the survival, genetic materials, hepatic transaminase activity and histopathology of loach (Misgurnus anguillicaudatus). The values of LC50 (24, 48, 72 and 96h) of imidacloprid were 167.7, 158.6, 147.9 and 145.8mg/L, respectively, and the safety concentration was 42.55mg/L. The erythrocyte micronuclei assays and the comet assay results showed that imidacloprid had genetic toxic effect on the loach erythrocytes. To assess the physiological and biochemical damage caused by imidacloprid, the activities of hepatic glutamic-pyruvic transaminase (GPT) and glutamic-oxalacetic transaminase (GOT) were measured and their values declined in treatment groups. Histological examination of testis revealed that imidacloprid treatment resulted in disorganized lobules and cysts structures. In the present work, we also investigated the joint toxicity of pesticides commonly used in paddy fields (imidacloprid and lambda-cyhalothrin) on M. anguillicaudatus, and confirmed that a synergistic effect existing in the binary mixtures. The results of our study provide relevant and comparable toxicity information that are useful for safety application of pesticides.

Introgressive Hybridization and the Evolution of Lake-Adapted Catostomid Fishes

Hybridization has been identified as a significant factor in the evolution of plants as groups of interbreeding species retain their phenotypic integrity despite gene exchange among forms. Recent studies have identified similar interactions in animals; however, the role of hybridization in the evolution of animals has been contested. Here we examine patterns of gene flow among four species of catostomid fishes from the Klamath and Rogue rivers using molecular and morphological traits. Catostomus rimiculus from the Rogue and Klamath basins represent a monophyletic group for nuclear and morphological traits; however, the Klamath form shares mtDNA lineages with other Klamath Basin species (C. snyderi, Chasmistes brevirostris, Deltistes luxatus). Within other Klamath Basin taxa, D. luxatus was largely fixed for alternate nuclear alleles relative to C. rimiculus, while Ch. brevirostris and C. snyderi exhibited a mixture of these alleles. Deltistes luxatus was the only Klamath Basin species that exhibited consistent covariation of nuclear and mitochondrial traits and was the primary source of mismatched mtDNA in Ch. brevirostris and C. snyderi, suggesting asymmetrical introgression into the latter species. In Upper Klamath Lake, D. luxatus spawning was more likely to overlap spatially and temporally with C. snyderi and Ch. brevirostris than either of those two with each other. The latter two species could not be distinguished with any molecular markers but were morphologically diagnosable in Upper Klamath Lake, where they were largely spatially and temporally segregated during spawning. We examine parallel evolution and syngameon hypotheses and conclude that observed patterns are most easily explained by introgressive hybridization among Klamath Basin catostomids.

High-Throughput Sequencing Identifies MicroRNAs from Posterior Intestine of Loach (Misgurnus anguillicaudatus) and Their Response to Intestinal Air-Breathing Inhibition

MicroRNAs (miRNAs) exert important roles in animal growth, immunity, and development, and regulate gene expression at the post-transcriptional level. Knowledges about the diversities of miRNAs and their roles in accessory air-breathing organs (ABOs) of fish remain unknown. In this work, we used high-throughput sequencing to identify known and novel miRNAs from the posterior intestine, an important ABO, in loach (Misgurnus anguillicaudatus) under normal and intestinal air-breathing inhibited conditions. A total of 204 known and 84 novel miRNAs were identified, while 47 miRNAs were differentially expressed between the two small RNA libraries (i.e. between the normal and intestinal air-breathing inhibited group). Potential miRNA target genes were predicted by combining our transcriptome data of the posterior intestine of the loach under the same conditions, and then annotated using COG, GO, KEGG, Swissprot and Nr databases. The regulatory networks of miRNAs and their target genes were analyzed. The abundances of nine known miRNAs were validated by qRT-PCR. The relative expression profiles of six known miRNAs and their eight corresponding target genes, and two novel potential miRNAs were also detected. Histological characteristics of the posterior intestines in both normal and air-breathing inhibited group were further analyzed. This study contributes to our understanding on the functions and molecular regulatory mechanisms of miRNAs in accessory air-breathing organs of fish.

Transforming growth factor-β1 expression in endangered age-0 shortnose suckers (Chasmistes brevirostris) from Upper Klamath Lake, OR relative to histopathology, meristic, spatial, and temporal data

During July-September of 2008, 2009, and 2010 endangered age-0 juvenile shortnose suckers were sampled from Upper Klamath Lake, OR in a health evaluation that included the measurement of transforming growth factor - beta (TGF-β) expression in spleen in combination with a histopathology assessment. This analysis was performed to determine if the expression of this immuno-regulator could be used as a component of a larger health evaluation intended to identify potential risk-factors that may help to explain why very few of these fish survive to age-1. Potential associations between TGF-β1 expression, histopathological findings, meristic data as well as temporal and spatial data were evaluated using analysis-of-variance. In this analysis, the absence or presence of opercula deformity and hepatic cell necrosis were identified as significant factors in accounting for the variance in TGF-β1 expression observed in age-0 shortnose suckers (n = 122, squared multiple R = 0.989). Location of sample collection and the absence or presence of anchor worms (Lernaea spp.) were identified as significant cofactors. The actual mechanisms involved with these relationships have yet to be determined. The strength, however, of our findings support the concept of using TGF-β1 expression as part of a broader fish health assessment and suggests the potential for using additional immunologic measures in future studies. Specifically, our results indicate that the measure of TGF-β1 expression in age-0 shortnose sucker health assessments can facilitate the process of identifying disease risks that are associated with the documented lack of recruitment into the adult population.

Asexual Reproduction Does Not Apparently Increase the Rate of Chromosomal Evolution: Karyotype Stability in Diploid and Triploid Clonal Hybrid Fish (Cobitis, Cypriniformes, Teleostei)

Interspecific hybridization, polyploidization and transitions from sexuality to asexuality considerably affect organismal genomes. Especially the last mentioned process has been assumed to play a significant role in the initiation of chromosomal rearrangements, causing increased rates of karyotype evolution. We used cytogenetic analysis and molecular dating of cladogenetic events to compare the rate of changes of chromosome morphology and karyotype in asexually and sexually reproducing counterparts in European spined loach fish (Cobitis). We studied metaphases of three sexually reproducing species and their diploid and polyploid hybrid clones of different age of origin. The material includes artificial F1 hybrid strains, representatives of lineage originated in Holocene epoch, and also individuals of an oldest known age to date (roughly 0.37 MYA). Thereafter we applied GISH technique as a marker to differentiate parental chromosomal sets in hybrids. Although the sexual species accumulated remarkable chromosomal rearrangements after their speciation, we observed no differences in chromosome numbers and/or morphology among karyotypes of asexual hybrids. These hybrids possess chromosome sets originating from respective parental species with no cytogenetically detectable recombinations, suggesting their integrity even in a long term. The switch to asexual reproduction thus did not provoke any significant acceleration of the rate of chromosomal evolution in Cobitis. Asexual animals described in other case studies reproduce ameiotically, while Cobitis hybrids described here produce eggs likely through modified meiosis. Therefore, our findings indicate that the effect of asexuality on the rate of chromosomal change may be context-dependent rather than universal and related to particular type of asexual reproduction.

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.

Karyotype differentiation in 19 species of river loach fishes (Nemacheilidae, Teleostei): extensive variability associated with rDNA and heterochromatin distribution and its phylogenetic and ecological interpretation

BACKGROUND

Loaches of the family Nemacheilidae are one of the most speciose elements of Palearctic freshwater ichthyofauna and have undergone rapid ecological adaptations and colonizations. Their cytotaxonomy is largely unexplored; with the impact of cytogenetical changes on this evolutionary diversification still unknown. An extensive cytogenetical survey was performed in 19 nemacheilid species using both conventional (Giemsa staining, C- banding, Ag- and Chromomycin A3/DAPI stainings) and molecular (fluorescence in situ hybridization with 5S rDNA, 45S rDNA, and telomeric (TTAGGG)n probes) methods. A phylogenetic tree of the analysed specimens was constructed based on one mitochondrial (cytochrome b) and two nuclear (RAG1, IRBP) genes.

RESULTS

Seventeen species showed karyotypes composed of 2n = 50 chromosomes but differentiated by fundamental chromosome number (NF = 68-90). Nemachilichthys ruppelli (2n = 38) and Schistura notostigma (2n = 44-48) displayed reduced 2n with an elevated number of large metacentric chromosomes. Only Schistura fasciolata showed morphologically differentiated sex chromosomes with a multiple system of the XY1Y2 type. Chromomycin A3 (CMA3)- fluorescence revealed interspecific heterogeneity in the distribution of GC-rich heterochromatin including its otherwise very rare association with 5S rDNA sites. The 45S rDNA sites were mostly located on a single chromosome pair contrasting markedly with a pattern of two (Barbatula barbatula, Nemacheilus binotatus, N. ruppelli) to 20 sites (Physoschistura sp.) of 5S rDNA. The cytogenetic changes did not follow the phylogenetic relationships between the samples. A high number of 5S rDNA sites was present in species with small effective population sizes.

CONCLUSION

Despite a prevailing conservatism of 2n, Nemacheilidae exhibited a remarkable cytogenetic variability on microstructural level. We suggest an important role for pericentric inversions, tandem and centric fusions in nemacheilid karyotype differentiation. Short repetitive sequences, genetic drift, founder effect, as well as the involvement of transposable elements in the dispersion of ribosomal DNA sites, might also have played a role in evolutionary processes such as reproductive isolation. These remarkable dynamics of their genomes qualify river loaches as a model for the study of the cytogenetic background of major evolutionary processes such as radiation, endemism and colonization of a wide range of habitats.