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

The Structure of Digestive Tract Coordinating Digestion and Respiration in an Air-Breathing Weatherloach, Misgurnus anguillicaudatus

To clarify how the digestive tract of the weatherloach, Misgurnus anguillicaudatus, serves a dual function of digestion and respiration simultaneously, the histological structures of its digestive tract, the passage of digesta and air passing through its intestine and the rate of intestinal evacuation have been studied. The results indicate that the digestive tract is divided into five functional regions, i.e., esophagus, anterior intestine, middle intestine, posterior intestine and rectum. The diverse intestinal structures have the specialized function of coordinating digestion and respiration. An X-ray barium meal examination showed in the normal breathing state, the contents of the intestine are diffusely semifluid, and air is distributed as bubbles in the dorsal intestine 2 h after feeding. After 5 h, the contents accumulated in the mid and posterior intestine, and gas flowed above the contents as bundles. After 8 h, the intestinal food was basically evacuated. In the intestinal air-breathing restricted group, the contents of the intestine remained diffuse, and a large number of digesta entered and remained in the rectum after 5 h. After the inhibition was relieved, the contents of the rectum were rapidly discharged. Measurement of the intestinal evacuation rate in the intestine showed that the evacuation of the intestinal contents lagged behind that of the normal group in the air-breathing restricted group. Compared to the normal state and inhibited GAB (gastrointestinal air breathing), we could deduce that GAB could promote the movement of the intestine.

Gut-immunity modulation in Lepidocephalichthys guntea during Aeromonas hydrophila-infection and recovery assessed with transcriptome data

The fish immune system, which consists of innate and adaptive immunologic processes, defends against viruses, bacteria, fungi, and parasites. The gut immunity is an integral part of the host immune system that controls immunological homeostasis, hosts' interactions with their microbiomes, and provides defence against a number of intestinal infections. Lepidocephalichthys guntea, a facultative air-breathing fish, was experimentally infected with Aeromonas hydrophila using intraperitoneal injection followed by bath challenge, and transcriptome data were used to examine the gut immune responses during disease progression and recovery from the diseased state without the use of medication. For the control or uninfected fish (FGC) and the infected fish that were kept for seven days (FGE1) and fifteen days (FGE2), separate water tanks were set up. Coding DNA sequences (CDS) for FGC and FGE1, FGC and FGE2, and FGE1 and FGE2 were analyzed for differential gene expression (DGE). The presence and expression of genes involved in the T cell receptor (TCR) signalling pathway, natural killer (NK) cell-mediated cytotoxicity pathway, and complement-mediated pathway, along with a large number of other immune-related proteins, and heat shock protein (HSPs) under various experimental conditions and its relationship to immune modulation of the fish gut was the primary focus of this study. Significant up-and-down regulation of these pathways shows that, in FGE1, the fish's innate immune system was engaged, whereas in FGE2, the majority of innate immune mechanisms were repressed, and adaptive immunity was activated. Expression of genes related to the immune system and heat-shock proteins was induced during this host's immunological response, and this information was then used to build a thorough network relating to immunity and the heat-shock response. This is the first study to examine the relationship between pathogenic bacterial infection, disease reversal, and modification of innate and adaptive immunity as well as heat shock response.

The chromosome-level genome and key genes associated with mud-dwelling behavior and adaptations of hypoxia and noxious environments in loach (Misgurnus anguillicaudatus)

BACKGROUND

The loach (Misgurnus anguillicaudatus), the most widely distributed species of the family Cobitidae, displays a mud-dwelling behavior and intestinal air-breathing, inhabiting the muddy bottom of extensive freshwater habitats. However, lack of high-quality reference genome seriously limits the interpretation of the genetic basis of specialized adaptations of the loach to the adverse environments including but not limited to the extreme water temperature, hypoxic and noxious mud environment.

RESULTS

This study generated a 1.10-Gb high-quality, chromosome-anchored genome assembly, with a contig N50 of 3.83 Mb. Multiple comparative genomic analyses found that proto-oncogene c-Fos (fos), a regulator of bone development, is positively selected in loach. Knockout of fos (ID: Mis0086400.1) led to severe osteopetrosis and movement difficulties, combined with the comparison results of bone mineral density, supporting the hypothesis that fos is associated with loach mud-dwelling behavior. Based on genomic and transcriptomic analysis, we identified two key elements involved in the intestinal air-breathing of loach: a novel gene (ID: mis0158000.1) and heat shock protein beta-1 (hspb1). The flavin-containing monooxygenase 5 (fmo5) genes, central to xenobiotic metabolism, undergone expansion in loach and were identified as differentially expressed genes in a drug stress trial. A fmo5^-/- (ID: Mis0185930.1) loach displayed liver and intestine injury, indicating the importance of this gene to the adaptation of the loach to the noxious mud.

CONCLUSIONS

Our work provides valuable insights into the genetic basis of biological adaptation to adverse environments.

Morphological and Molecular Functional Evidence of the Pharyngeal Sac in the Digestive Tract of Silver Pomfret, Pampus argenteus

The pharyngeal sac is a comparatively rare organ in the digestive tract among teleost fishes. However, our understanding of this remarkable organ in the silver pomfret (Pampus argenteus) is limited. In the present study, we examined the various morphological and histological characteristics of the pharyngeal sac using histochemical techniques and electron microscopy. The pharyngeal sac showed unique characteristics such as well-developed muscular walls, weakly keratinized epithelium, numerous goblet cells, and needle-like processes on the papillae. The porous cavity of the papillae contained numerous adipocytes and was tightly enveloped by type I collagen fibers. These structures might provide mechanical protection and excellent biomechanical properties for grinding and shredding prey. A comparison of gene expression levels between the pharyngeal sac and esophagus using RNA-seq showed that phenotype-associated genes (epithelial genes and muscle genes) were upregulated, whereas genes related to nutrient digestion and absorption were downregulated in the pharyngeal sac. These results support the role of the pharyngeal sac in shredding and predigesting food. Overall, these findings provide a clearer understanding of the pharyngeal sac morphology and explain the morphological adaptations of the digestive tract for feeding on gelatinous prey. To our knowledge, this is the first report on pharyngeal sac gene expression in P. argenteus.

Depletion of LOXL2 improves respiratory capacity: From air-breathing fish to mammal under hypoxia

Air-breathing fish are fascinating because of their ability to survive under hypoxia for a long time by using air-breathing organs (ABOs). Fish ABOs are thought to resemble the mammal lung all along. However, the link between the two has not been studied in depth. Here, we reported a markedly improved respiratory capacity in mice under hypoxia by inhibiting lysyl oxidase-like 2 (LOXL2), inspired from the intestinal air-breathing of loach (Misgurnus anguillicaudatus). Moreover, a posterior intestine (an ABO) transcriptome analysis revealed that the deletion of Loxl2b obviously inhibited PI3K-AKT and TGF-β signaling, meanwhile, induced VEGF signaling, which could cause vasodilation and angiogenesis to improve the air-breathing ability of loach. The same phenomenon was found in LOXL2-inhibition mice under hypoxia, which significantly prolonged their living period relative to wild-type (WT) mice. In addition, compared with WT loach, Loxl2b^-/- loach presented enhanced anaerobic metabolism, which could also make itself to better survive in hypoxic environment. This should be the magic of air-breathing fish! Supplied from air-breathing fish, this study provides a novel means of improving respiratory capacity in mammal under hypoxia.

Comparative Transcriptomic Analysis of Regenerated Skins Provides Insights into Cutaneous Air-Breathing Formation in Fish

Cutaneous air-breathing is one of the air-breathing patterns in bimodal respiration fishes, while little is known about its underlying formation mechanisms. Here, we first investigated the skin regeneration of loach (Misgurnus anguillicaudatus, a cutaneous air-breathing fish) and yellow catfish (Pelteobagrus fulvidraco, a water-breathing fish) through morphological and histological observations. Then, the original skins (OS: MOS, POS) and regenerated skins (RS: MRS, PRS) when their capillaries were the most abundant (the structural foundation of air-breathing in fish) during healing, of the two fish species were collected for high-throughput RNA-seq. A total of 56,054 unigenes and 53,731 unigenes were assembled in loach and yellow catfish, respectively. A total of 640 (460 up- and 180 down-regulated) and 4446 (2340 up- and 2106 down-regulated) differentially expressed genes (DEGs) were respectively observed in RS/OS of loach and yellow catfish. Subsequently, the two DEG datasets were clustered in GO, KOG and KEGG databases, and further analyzed by comparison and screening. Consequently, tens of genes and thirteen key pathways were targeted, indicating that these genes and pathways had strong ties to cutaneous skin air-breathing in loach. This study provides new insights into the formation mechanism of cutaneous air-breathing and also offers a substantial contribution to the gene expression profiles of skin regeneration in fish.

Endogenic upregulations of HIF/VEGF signaling pathway genes promote air breathing organ angiogenesis in bimodal respiration fish

Air-breathing has evolved independently serval times with a variety of air-breathing organs (ABOs) in fish. The physiology of the air-breathing in bimodal respiration fish has been well understood, while studies on molecular mechanisms of the character are very limited. In the present study, we first determined the gill indexes of 110 fish species including 25 and 85 kinds of bimodal respiration fishes and non-air-breathing fishes, respectively. Then combined with histological observations of gills and ABOs/non-ABOs in three bimodal respiration fishes and two non-air breathing fishes, we found that the bimodal respiration fish was always of a degeneration gill and a well-vascularized ABO. Meanwhile, a comparative transcriptome analysis of posterior intestines, namely a well vascularized ABO in Misgurnus anguillicaudatus and a non-ABO in Leptobotia elongata, was performed to expound molecular variations of the air-breathing character. A total of 5,003 orthologous genes were identified. Among them, 1,189 orthologous genes were differentially expressed, which were enriched in 14 KEGG pathways. More specially, the expressions of hemoglobin genes and various HIF/VEGF signaling pathway genes were obviously upregulated in the ABO of M. anguillicaudatus. Moreover, we found that HIF-1α, VEGFAa, and MAP2K1 were co-expressed dramatically higher in ABOs of bimodal respiration fishes than those of non-ABOs of non-air-breathing fishes. These results indicated that the HIF/VEGF pathway played an important role in ABO angiogenesis/formation to promote fish to do aerial respiration. This study will contribute to our understanding of molecular mechanisms of air-breathing in fish.

Fibronectin 1B Gene Plays an Important Role in Loach Barbel Air-Breathing

Loach (Misgurnus anguillicaudatus) is well known to perform air-breathing through the posterior intestine and skin. However, we find here for the first time a unique central vascular structure in the loach barbel, with a blood–gas diffusion distance as short as that of the posterior intestine. Under acute hypoxia, the distance of loach barbels became significantly shorter. Moreover, barbel removal significantly decreased air-breathing frequency of the loach. These findings imply that the barbel is another air-breathing organ of the loach. For further investigation of loach barbel air-breathing, a transcriptome analysis of barbels with air exposure treatment was performed. A total of 2546 differentially expressed genes (DEGs) between the T-XU (air exposure) and C-XU (control) group were identified, and 13 key DEGs related to barbel air-breathing were screened out. On this foundation, sequence, expression, and location analysis results indicated an important positive role of fibronectin 1b (fn1b) in loach barbel air-breathing. We further generated an fn1b-depletion loach (MT for short) using the CRISPR/Cas9 technique. It was indicated that depletion of fn1b could weaker barbel air-breathing ability. In conclusion, due to nonlethal and regenerative characteristics, the loach barbel, a newly discovered and fn1b-related fish air-breathing organ, can be a good model for fish air-breathing research.

Clock Gene Expression in the Eye Exhibits Circadian Oscillation and Light Responsiveness but is Not Necessary for Nocturnal Locomotor Activity of Japanese Loach, Misgurnus anguillicaudatus

There are few model fish that are both edible and suitable for use in the laboratory. The Japanese loach (Misgurnus anguillicaudatus) is a traditional food in Japan, but is highly neglected despite its great nutritional value. To understand its circadian system and photic input pathway for synchronization of physiological activities to environmental light-dark cycles, we measured locomotor activity under light-dark and constant dark (DD) conditions. Locomotor activity was found to be higher in the nighttime than daytime, and its rhythmicity was weakened under DD conditions. The nocturnal activity of the Japanese loach is mainly controlled by environmental light, rather than the circadian clock. We explored the circadian regulation and light-responsiveness of clock gene expression in the eyes of loaches. The daily expression profiles of its mRNA revealed that most of the examined Cry and Per genes were likely regulated by internal circadian and/or environmental light signals. Among the Opsin genes transcribed in the eye, we detected the retinal photopigment porphyropsin at the protein level, which was lower than in mice. This property of loach eyes prompted us to analyze the locomotor activities of eye-enucleated fish. As a result, they still showed nocturnal circadian activity. Thus, it is likely that extraocular photoreceptive tissue(s) also contribute to the photic input pathway, although loach eyes are a circadian photosensitive tissue. This suggests that the loach mainly uses not its vision but other stimuli, such as mechanical or chemical stimuli, detected by barbels, to coordinate its nocturnal behavior.

Isolation, identification and characterisation of an emerging fish pathogen, Acinetobacter pittii, from diseased loach (Misgurnus anguillicaudatus) in China

Although members of the genus Acinetobacter have emerged as important nosocomial pathogens causing severe human infections, there are few reports about their occurrence as fish pathogens. In this study, five bacterial strains were isolated from diseased loach (Misgurnus anguillicaudatus) cultured in a farm in China. The diseased loach displayed shedding of skin mucus and many petechial haemorrhages all over the body. Based on sequence analyses of 16S rRNA and rpoB genes, the isolates were identified as Acinetobacter pittii. An experimental infection assay confirmed their pathogenicity to loach. The results of artificial infection in zebrafish (Barchydanio rerio) and nematode (Caenorhabditis elegans) suggested that, as well as loach, these A. pittii isolates are pathogenic and highly virulent to these organisms. Multilocus sequence typing analysis revealed that all the isolates belong to sequence type (ST) 839, which may be the dominant clone causing fish disease and exhibits a close phylogenetic relationship with ST396 from human clinical samples in Korea or Taiwan China. This is the first report demonstrating that A. pittii is an emerging causal agent of mass mortality in loach and poses significant risks to fish culturing besides causing human clinical infection worldwide.

Production of a mutant of large-scale loach Paramisgurnus dabryanus with skin pigmentation loss by genome editing with CRISPR/Cas9 system

CRISPR/Cas9 system has been developed as a highly efficient genome editing technology to specifically induce mutations in a few aquaculture species. In this study, we described induction of targeted gene (namely tyrosinase, tyr) mutations in large-scale loach Paramisgurnus dabryanus, an important aquaculture fish species and a potential model organism for studies of intestinal air-breathing function, using the CRISPR/Cas9 system. Tyr gene in large-scale loach was firstly cloned and then its expressions were investigated. Two guide RNAs (gRNAs) were designed and separately transformed with Cas9 in the loach. 89.4% and 96.1% of injected loach juveniles respectively displayed a graded loss of pigmentation for the two gRNAs, in other words, for target 1 and target 2. We classified the injected loach juveniles into five groups according to their skin color phenotypes, including four albino groups and one wild-type-like group. And one of them was clear albino group, which was of high ornamental and commercial value. More than 50 clones for each albino transformant with a visible phenotype in each target were randomly selected and sequenced. Results obtained here showed that along with the increase of pigmentation, wild-type alleles appeared in the injected loach juveniles more often and insertion/deletion alleles less frequently. This study demonstrated that CRISPR/Cas9 system could be practically performed to modify large-scale loach tyr to produce an albino mutant of high ornamental and commercial value, and for the first time showed successful use of the CRISPR/Cas9 system for genome editing in a Cobitidae species.

Full-length transcriptome of Misgurnus anguillicaudatus provides insights into evolution of genus Misgurnus

Reconstruction and annotation of transcripts, particularly for a species without reference genome, plays a critical role in gene discovery, investigation of genomic signatures, and genome annotation in the pre-genomic era. This study generated 33,330 full-length transcripts of diploid M. anguillicaudatus using PacBio SMRT Sequencing. A total of 6,918 gene families were identified with two or more isoforms, and 26,683 complete ORFs with an average length of 1,497 bp were detected. Totally, 1,208 high-confidence lncRNAs were identified, and most of these appeared to be precursor transcripts of miRNAs or snoRNAs. Phylogenetic tree of the Misgurnus species was inferred based on the 1,905 single copy orthologous genes. The tetraploid and diploid M. anguillicaudatus grouped into a clade, and M. bipartitus showed a closer relationship with the M. anguillicaudatus. The overall evolutionary rates of tetraploid M. anguillicaudatus were significantly higher than those of other Misgurnus species. Meanwhile, 28 positively selected genes were identified in M. anguillicaudatus clade. These positively selected genes may play critical roles in the adaptation to various habitat environments for M. anguillicaudatus. This study could facilitate further exploration of the genomic signatures of M. anguillicaudatus and provide potential insights into unveiling the evolutionary history of tetraploid loach.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.

RNA-Seq Analysis of the Antioxidant Status and Immune Response of Portunus trituberculatus Following Aerial Exposure

Desiccation tolerance has been long considered as an important trait for the life survival under acute environmental stress. One of the biggest problems for modern commercial crab farming is desiccation during transportation; high mortality could occur following the aerial exposure. In this regard, here, we utilized RNA-seq-based transcriptome profiling to characterize the molecular responses of swimming crab in response to aerial exposure. In present study, following aerial exposure, the gill samples were sequenced at 0, 6, 12, and 18 h. And the sequenced reads were assembled into 274,594 contigs, with average length of 735.59 bp and N50 size of 1262 bp. After differential expression analysis, a total of 1572 genes were captured significantly differentially expressed, and were categorized into antioxidant/oxidative stress response, chaperones/heat shock proteins, immune alteration, cell proliferation/apoptosis, and cytoskeletal. Our analysis revealed the dramatic tissue oxidant stress and the alteration of the tissue epithelial integrity, especially many genes that have not been reported in crab species. With the limited functional information in crab, further studies are needed and underway in our lab to further characterize the key cellular actors governing the crab tolerance to aerial exposure. Taken together, our results provide molecular resources for further identification of key genes for desiccation tolerance, and to facilitate the molecular selection and breeding of desiccation tolerant strain and family.