Molecular characterization, spatio-temporal expression patterns of crtc2 gene and its immune roles in yellow catfish (Pelteobagrus fulvidraco).

cAMP response element binding (CREB) protein 2 (CRTC2) is a transcriptional coactivator of CREB and plays an important role in the immune system. Thus far, the physiological roles of Crtc2 in teleost are still poorly understood. In this study, the crtc2 gene was identified and characterized from yellow catfish (Pelteobagrus fulvidraco; therefore, the gene is termed as pfcrtc2), and its evolutionary and molecular characteristics as well as potential immunity-related roles were investigated. Our results showed that the open reading frame of pfcrtc2 was 2346 bp in length, encoding a protein with 781 amino acids. Gene structure analysis revealed its existence of 14 exons and 13 introns. A phylogenetic analysis proved that the tree of crtc2 was clustered into five groups, exhibiting a similar evolutionary topology with species evolution. Multiple protein sequences alignment demonstrated high conservation of the crtc2 in various vertebrates with similar structure. Syntenic and gene structural comparisons further established that crtc2 was highly conserved, implying its similar roles in diverse vertebrates. Tissue distribution pattern detected by quantitative real-time PCR showed that the pfcrtc2 gene was almost expressed in all detected tissues except for eyes, with the highest expression levels in the gonad, indicating that Crtc2 may play important roles in various tissues. In addition, pfcrtc2 was transcribed at all developmental stages in yellow catfish, showing the highest expression levels at 12 h after fertilization. Finally, the transcriptional profiles of crtc2 were significantly increased in yellow catfishes injected with Aeromonas hydrophila or Poly I:C, which shared a consistent change pattern with four immune-related genes including IL-17A, IL-10, MAPKp38, and NF-κBp65, suggesting pfCrtc2 may play critical roles in preventing both exogenous bacteria and virus invasion. In summary, our findings lay a solid foundation for further studies on the functions of pfcrtc2, and provide novel genetic loci for developing new strategies to control disease outbreak in teleost.

Chromosome-level genome assemblies of vulnerable male and female elongate loach (Leptobotia elongata).

Endemic to the upper and middle reaches of the Yangtze River in China, elongate loach (Leptobotia elongata) has become a vulnerable species mainly due to overfishing and habitat destruction. Thus far, no genome data of this species are reported. As a result, lacking of such genomic information has restricted practical conservation and utilization of this economic fish. Here, we constructed chromosome-level genome assemblies for both male and female elongate loach by integration of MGI, PacBio HiFi and Hi-C sequencing technologies. Two primary genome assemblies (586-Mb and 589-Mb) were obtained for female and male fishes, respectively. Indeed, 98.22% and 98.61% of the contig sequences were anchored onto 25 chromosomes, with identification of 26.22% and 25.92% repeat contents in both assembled genomes. Meanwhile, a total of 25,215 and 25,253 protein-coding genes were annotated, of which 97.41% and 98.8% could be predicted with functions. Taken together, our genome data presented here provide a valuable genomic resource for in-depth evolutionary and functional research, as well as molecular breeding and conservation of this economic fish species.

Molecular Identification and Functional Characterization of LC-PUFA Biosynthesis Elongase (elovl2) Gene in Chinese Sturgeon (Acipenser sinensis).

Elongases of very-long-chain fatty acids (Elovls) are critical rate-limiting enzymes that are involved in LC-PUFA biosynthesis through catalyzing the two-carbon elongation of a pre-existing fatty acyl chain. Thus far, several Elovls have been extensively studied in teleost. However, the functional and physiological roles of Elovls in chondrichthyans have rarely been reported. In this study, we identified and characterized elovl2 from the endangered Chinese sturgeon (Acipenser sinensis) by whole genome scanning. The results show that the coding sequence of elovl2 was 894 bp in length, for a putative protein of 297 amnio acids. Comparative genomic analyses indicated that Chinese sturgeon elovl2 was evolutionarily conserved. Functional characterization in yeast demonstrated that the Chinese sturgeon Elovl2 could efficiently elongate C20 (ARA and EPA) and C22 (22:4n-6 and 22:5n-3) substrates, confirming its critical roles in LC-PUFA biosynthesis. Spatial and temporal expression analyses showed high elovl2 mRNA levels were detected in the liver and brain and showed an increase trend both in embryonic and post-hatching stages. Interestingly, diets with vegetable oils as lipid sources could significantly induce the high expression of elovl2 in Chinese sturgeon, implying that the endogenous LC-PUFA biosynthesis pathway was stimulated by lack of LC-PUFA in their diets. Our findings will enhance our understanding about the evolutionary and functional roles of elovl2 and provide novel insights into the LC-PUFA biosynthesis mechanism in vertebrates.

5 mC modification of steroid hormone biosynthesis-related genes orchestrates feminization of channel catfish induced by high-temperature.

To elucidate the mechanism behind channel catfish feminization induced by high temperature, gonad samples were collected from XY pseudo-females and wild-type females and subjected to high-throughput sequencing for Whole-Genome-Bisulfite-Seq (WGBS) and transcriptome sequencing (RNA-Seq). The analysis revealed 50 differentially methylated genes between wild-type females and XY pseudo-females, identified through the analysis of KEGG pathways and GO enrichment in the promoter of the genome and differentially methylated regions (DMRs). Among these genes, multiple differential methylation sites observed within the srd5a2 gene. Repeatability tests confirmed 7 differential methylation sites in the srd5a2 gene in XY pseudo-females compared to normal males, with 1 specific differential methylation site (16608174) distinguishing XY pseudo-females from normal females. Interestingly, the expression of these genes in the transcriptome showed no difference between wild-type females and XY pseudo-females. Our study concluded that methylation of the srd5a2 gene sequence leads to decreased expression, which inhibits testosterone synthesis while promoting the synthesis of 17ß-estradiol from testosterone. This underscores the significance of the srd5a2 gene in the sexual differentiation of channel catfish, as indicated by the ipu00140 KEGG pathway analysis.

Characterization of three lamp genes from largemouth bass (Micropterus salmoides): molecular cloning, expression patterns, and their transcriptional levels in response to fast and refeeding strategy.

Lysosomes-associated membrane proteins (LAMPs), a family of glycosylated proteins and major constituents of the lysosomal membranes, play a dominant role in various cellular processes, including phagocytosis, autophagy and immunity in mammals. However, their roles in aquatic species remain poorly known. In the present study, three lamp genes were cloned and characterized from Micropterus salmoides. Subsequently, their transcriptional levels in response to different nutritional status were investigated. The full-length coding sequences of lamp1, lamp2 and lamp3 were 1251bp, 1224bp and 771bp, encoding 416, 407 and 256 amino acids, respectively. Multiple sequence alignment showed that LAMP1-3 were highly conserved among the different fish species, respectively. 3-D structure prediction, genomic survey, and phylogenetic analysis were further confirmed that these genes are widely existed in vertebrates. The mRNA expression of the three genes was ubiquitously expressed in all selected tissues, including liver, brain, gill, heart, muscle, spleen, kidney, stomach, adipose and intestine, lamp1 shows highly transcript levels in brain and muscle, lamp2 displays highly expression level in heart, muscle and spleen, but lamp3 shows highly transcript level in spleen, liver and kidney. To analyze the function of the three genes under starvation stress in largemouth bass, three experimental treatment groups (fasted group and refeeding group, control group) were established in the current study. The results indicated that the expression of lamp1 was significant induced after starvation, and then returned to normal levels after refeeding in the liver. The expression of lamp2 and lamp3 exhibited the same trend in the liver. In addition, in the spleen and the kidney, the transcript level of lamp1 and lamp2 was remarkably increased in the fasted treatment group and slightly decreased in the refed treatment group, respectively. Collectively, our findings suggest that three lamp genes may have differential function in the immune and energetic organism in largemouth bass, which is helpful in understanding roles of lamps in aquatic species.

Characterization of TLR1 and expression profiling of TLR signaling pathway related genes in response to Aeromonas hydrophila challenge in hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂).

Toll-like receptor 1 (TLR1) mediates the innate immune response to a variety of microbes through recognizing cell wall components (such as bacterial lipoproteins) in mammals. However, the detailed molecular mechanism of TLR1 involved in pathogen immunity in the representative hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂) has not been well studied. In the present study, we identified the TLR1 gene from the hybrid yellow catfish, and further comparative synteny data from multiple species confirmed that the TLR1 gene is highly conserved in teleosts. Phylogenetic analysis revealed distinguishable TLR1s in diverse taxa, suggesting consistence in evolution of the TLR1 proteins with various species. Structural prediction indicated that the three-dimensional structures of TLR1 proteins are relatively conserved among different taxa. Positive selection analysis showed that purifying selection dominated the evolutionary process of TLR1s and TLR1-TIR domain in both vertebrates and invertebrates. Expression pattern analysis based on the tissue distribution showed that TLR1 mainly transcribed in the gonad, gallbladder and kidney, and the mRNA levels of TLR1 in kidney were remarkably up-regulated after Aeromonas hydrophila stimulation, indicating that TLR1 participates in the inflammatory responses to exogenous pathogen infection in hybrid yellow catfish. Homologous sequence alignment and chromosomal location indicated that the TLR signaling pathway is very conserved in the hybrid yellow catfish. The expression patterns of TLR signaling pathway related genes (TLR1- TLR2 - MyD88 - FADD - Caspase 8) were consistent after pathogen stimulation, revealing that the TLR signaling pathway is triggered and activated after A. hydrophila infection. Our findings will lay a solid foundation for better understanding the immune roles of TLR1 in teleosts, as well as provide basic data for developing strategies to control disease outbreak in hybrid yellow catfish.

Molecular characterization of a teleost-specific toll-like receptor 22 (tlr22) gene from yellow catfish (Pelteobagrus fulvidraco) and its transcriptional change in response to poly I:C and Aeromonas hydrophila stimuli.

Toll-like receptors (TLRs) are a class of pattern recognition receptors (PRRs) that can recognize pathogen-associated molecular patterns (PMPs) and play important roles in the innate immune system in vertebrates. In this study, we identified a teleost-specific tlr22 gene from yellow catfish (Pelteobagrus fulvidraco) and its immune roles in response to different pathogens were also determined. The open reading frame (ORF) of the tlr22 was 2892 bp in length, encoding a protein of 963 amino acids. Multiple protein sequences alignment, secondary and three-dimensional structure analyses revealed that TLR22 is highly conserved among different fish species. Phylogenetic analysis showed that the phylogenetic topology was divided into six families of TLR1, TLR3, TLR4, TLR5, TLR7 and TLR11, and TLR22 subfamily was clustered into TLR11 family. Meanwhile, synteny and gene structure comparisons revealed functional and evolutionary conservation of the tlr22 gene in teleosts. Furthermore, tlr22 gene was shown to be widely expressed in detected tissues except barbel and eye, with highest expression level in liver. The transcription of tlr22 was significantly increased in spleen, kidney, liver and gill tissues at different timepoints after Poly I:C infection, suggesting TLR22 plays critical roles in defensing virus invasion. Similarly, the transcription of tlr22 was also dramatically up-regulated in spleen, kidney and gill tissues with different patterns after Aeromonas hydrophila infection, indicating that TLR22 is also involved in resisting bacteria invasion. Our findings will provide a solid basis for the investigation the immune functions of tlr22 gene in teleosts, as well as provide useful information for disease control and treatment for yellow catfish.

Toll-like Receptor 3 in the Hybrid Yellow Catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂): Protein Structure, Evolution and Immune Response to Exogenous Aeromonas hydrophila and Poly (I:C) Stimuli.

As a major mediator of cellular response to viral infection in mammals, Toll-like receptor 3 (TLR3) was proved to respond to double-stranded RNA (dsRNA). However, the molecular mechanism by which TLR3 functions in the viral infection response in teleosts remains to be investigated. In this study, the Toll-like receptor 3 gene of the hybrid yellow catfish was identified and characterized by comparative genomics. Furthermore, multiple sequence alignment, genomic synteny and phylogenetic analysis suggested that the homologous TLR3 genes were unique to teleosts. Gene structure analysis showed that five exons and four introns were common components of TLR3s in the 12 examined species, and interestingly the third exon in teleosts was the same length of 194 bp. Genomic synteny analysis indicated that TLR3s were highly conserved in various teleosts, with similar organizations of gene arrangement. De novo predictions showed that TLR3s were horseshoe-shaped in multiple taxa except for avian (with a round-shaped structure). Phylogenetic topology showed that the evolution of TLR3 was consistent with the evolution of the studied species. Selection analysis showed that the evolution rates of TLR3 proteins were usually higher than those of TLR3-TIR domains, indicating that the latter were more conserved. Tissue distribution analysis showed that TLR3s were widely distributed in the 12 tested tissues, with the highest transcriptions in liver and intestine. In addition, the transcription levels of TLR3 were significantly increased in immune-related tissues after infection of exogenous Aeromonas hydrophila and poly (I:C). Molecular docking showed that TLR3 in teleosts forms a complex with poly (I:C). In summary, our present results suggest that TLR3 is a pattern recognition receptor (PRR) gene in the immune response to pathogen infections in hybrid yellow catfish.

Cadmium induced oxidative stress, endoplasmic reticulum (ER) stress and apoptosis with compensative responses towards the up-regulation of ribosome, protein processing in the ER, and protein export pathways in the liver of zebrafish.

The present study identified that exposure to 5, 10, and 20 µg/L Cd for 48 days reduced growth, increased Cd accumulation and levels of reactive oxygen species (ROS) and lipid peroxidation, and induced ER stress and cellular apoptosis in the liver in a dose-dependent manner. However, the survival rate was not affected by Cd. The increased production of ROS might result from reduced catalase (CAT) and copper/zinc-superoxide dismutase (Cu/Zn-SOD) activities, which might trigger ER stress pathways and subsequently induce apoptotic responses, ultimately leading to growth inhibition. Transcriptomic analyses indicated that the differentially expressed genes (DEGs) involved in metabolic pathways were significantly enriched and dysregulated by Cd, suggesting that metabolic disturbances may contribute to Cd toxicity. However, there were increases in glutathione peroxidase (GPX) activity, protein levels of metallothioneins (MTs) and heat shock protein 70 (HSP70), and mRNA levels of sod1, cat, gpx, mt2, and hsp70. Furthermore, DEGs related to ribosome, protein processing in the ER, and protein export pathways were significantly enriched and up-regulated by Cd. These increases may be compensatory responses following oxidative stress, ER stress, and apoptosis to resist negative effects. Taken together, we demonstrated that environmentally relevant levels of Cd induced adaptive responses with compensatory mechanisms in fish, which may help to maintain fish survival at the cost of growth.

Phylogenetic Analysis of Core Melanin Synthesis Genes Provides Novel Insights Into the Molecular Basis of Albinism in Fish.

Melanin is the most prevalent pigment in animals. Its synthesis involves a series of functional genes. Particularly, teleosts have more copies of these genes related to the melanin synthesis than tetrapods. Despite the increasing number of available vertebrate genomes, a few systematically genomic studies were reported to identify and compare these core genes for the melanin synthesis. Here, we performed a comparative genomic analysis on several core genes, including tyrosinase genes (tyr, tyrp1, and tyrp2), premelanosome protein (pmel), microphthalmia-associated transcription factor (mitf), and solute carrier family 24 member 5 (slc24a5), based on 90 representative vertebrate genomes. Gene number and mutation identification suggest that loss-of-function mutations in these core genes may interact to generate an albinism phenotype. We found nonsense mutations in tyrp1a and pmelb of an albino golden-line barbel fish, in pmelb of an albino deep-sea snailfish (Pseudoliparis swirei), in slc24a5 of cave-restricted Mexican tetra (Astyanax mexicanus, cavefish population), and in mitf of a transparent icefish (Protosalanx hyalocranius). Convergent evolution may explain this phenomenon since nonsense mutations in these core genes for melanin synthesis have been identified across diverse albino fishes. These newly identified nonsense mutations and gene loss will provide molecular guidance for ornamental fish breeding, further enhancing our in-depth understanding of human skin coloration.

The complete mitochondrial genome of the intertidal spider (Desis jiaxiangi) provides novel insights into the adaptive evolution of the mitogenome and the evolution of spiders.

BACKGROUND: Although almost all extant spider species live in terrestrial environments, a few species live fully submerged in freshwater or seawater. The intertidal spiders (genus Desis) built silk nests within coral crevices can survive submerged in high tides. The diving bell spider, Argyroneta aquatica, resides in a similar dynamic environment but exclusively in freshwater. Given the pivotal role played by mitochondria in supplying most energy for physiological activity via oxidative phosphorylation and the environment, herein we sequenced the complete mitogenome of Desis jiaxiangi to investigate the adaptive evolution of the aquatic spider mitogenomes and the evolution of spiders. RESULTS: We assembled a complete mitogenome of the intertidal spider Desis jiaxiangi and performed comparative mitochondrial analyses of data set comprising of Desis jiaxiangi and other 45 previously published spider mitogenome sequences, including that of Argyroneta aquatica. We found a unique transposition of trnL2 and trnN genes in Desis jiaxiangi. Our robust phylogenetic topology clearly deciphered the evolutionary relationships between Desis jiaxiangi and Argyroneta aquatica as well as other spiders. We dated the divergence of Desis jiaxiangi and Argyroneta aquatica to the late Cretaceous at ~ 98 Ma. Our selection analyses detected a positive selection signal in the nd4 gene of the aquatic branch comprising both Desis jiaxiangi and Argyroneta aquatica. Surprisingly, Pirata subpiraticus, Hypochilus thorelli, and Argyroneta aquatica each had a higher Ka/Ks value in the 13 PCGs dataset among 46 taxa with complete mitogenomes, and these three species also showed positive selection signal in the nd6 gene. CONCLUSIONS: Our finding of the unique transposition of trnL2 and trnN genes indicates that these genes may have experienced rearrangements in the history of intertidal spider evolution. The positive selection signals in the nd4 and nd6 genes might enable a better understanding of the spider metabolic adaptations in relation to different environments. Our construction of a novel mitogenome for the intertidal spider thus sheds light on the evolutionary history of spiders and their mitogenomes.

MRAP2 Interaction with Melanocortin-4 Receptor in SnakeHead (Channa argus).

The melanocortin-4 receptor (MC4R) plays an important role in the regulation of food intake and energy expenditure. Melanocortin-2 receptor accessory protein 2 (MRAP2) modulates trafficking, ligand binding, and signaling of MC4R. The Northern snakehead (Channa argus) is an economically important freshwater fish native to East Asia. To explore potential interaction between snakehead MC4R and MRAP2, herein we cloned snakehead mc4r and mrap2. The snakehead mc4r consisted of a 984 bp open reading frame encoding a protein of 327 amino acids, while snakehead mrap2 contained a 693 bp open reading frame encoding a protein of 230 amino acids. Synteny analysis indicated that mc4r was highly conserved with similar gene arrangement, while mrap2 contained two isoforms in teleost with different gene orders. Snakehead mc4r was primarily expressed in the brain, whereas mrap2 was expressed in the brain and intestine. Snakehead mc4r and mrap2 expression was modulated by fasting and refeeding. Further pharmacological experiments showed that the cloned snakehead MC4R was functional, capable of binding to peptide agonists and increasing intracellular cAMP production in a dose-dependent manner. Snakehead MC4R exhibited high constitutive activity. MRAP2 significantly decreased basal and agonist-stimulated cAMP signaling. These findings suggest that snakehead MC4R might be involved in energy balance regulation by interacting with MRAP2. Further studies are needed to elucidate MC4R in regulating diverse physiological processes in snakehead.

Complete mitochondrial genome of broadbanded cardinalfish (Ostorhinchus fasciatus) and phylogenetic analysis.

In this study, the complete mitochondrial genome (mitogenome) of Ostorhinchus fasciatus was first determined and its phylogenetic position was investigated. The mitogenome was 16568 bp long and showed a typical teleost orders, containing 13 protein-coding genes (PCGs), 2 ribosome RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and a D-loop region. The overall nucleotide composition included A, 25.89%; C, 30.40%; G, 17.46%; and T, 26.26%. Except for nad6 was located on the light strand, the other PCGs were encoded on the heavy strand. Phylogenetic analysis suggested that O. fasciatus shared a close relationship with Sphaeramia orbicularis and Pterapogon kauderni.

Characterization of two kcnk3 genes in rabbitfish (Siganus canaliculatus): Molecular cloning, distribution patterns and their potential roles in fatty acids metabolism and osmoregulation.

KCNK3 is a two-pore-domain (K2P) potassium channel involved in maintaining ion homeostasis, mediating thermogenesis, controlling breath and modulating electrical membrane potential. Although the functions of this channel have been widely described in mammals, its roles in fishes are still rarely known. Here, we identified two kcnk3 genes from the euryhaline rabbitfish (Siganus canaliculatus), and their roles related to fatty acids metabolism and osmoregulation were investigated. The open reading frames of kcnk3a and kcnk3b were 1203 and 1176 bp in length, encoding 400 and 391 amino acids respectively. Multiple sequences alignment and phylogenetic analysis revealed that the two isotypes of kcnk3 were extensively presented in fishes. Quantitative real-time PCRs indicated that both genes were widely distributed in examined tissues but showed different patterns. kcnk3a primary distributed in adipose, eye, heart, and spleen tissues, while kcnk3b was mainly detectable in heart, kidney, muscle and spleen tissues. In vivo experiments showed that fish fed diets with fish oil as dietary lipid (rich in long chain polyunsaturated fatty acids, LC-PUFA) induced higher mRNA expression levels of kcnk3 genes in comparison with fish fed with plant oil diet at two different salinity environments (32 and 15‰). Meanwhile, the expression levels of kcnk3 genes were higher in seawater (32‰) than that in brackish water (15‰) when fishes were fed with both types of feeds. In vitro experiments with rabbitfish hepatocytes showed that LC-PUFA significantly improved hepatic kcnk3a expression level compared with treatment of linolenic acid. These results suggest that two kcnk3 genes are widely existed and they might be functionally related to fatty acids metabolism and osmoregulation in the rabbitfish.

Molecular characterization of two leptin genes and their transcriptional changes in response to fasting and refeeding in Northern snakehead (Channa argus).

Leptin has been proved to play critical roles in energy metabolism, body weight regulation, food intake, reproduction and immunity in mammals. However, its roles are still largely unclear in fish. Here, we report two leptin genes (lepA and lepB) from the Northern snakehead (Channa argus) and their transcriptions in response to different feeding status. The snakehead lepA is 781 bp in length and contains a 480 bp open reading frame (ORF) encoding a 159-aa protein, while the snakehead lepB is 553 bp in length and contains a 477 bp ORF encoding a 158-aa protein. Multi-sequences alignment, three-dimensional (3D) model prediction, syntenic and genomic comparison, and phylogenetic analysis confirm two leptin genes are widely existing in teleost. Tissue distribution revealed that the two leptin genes exhibit different patterns. In a post-prandial experiment, the hepatic lepA and brain lepB showed a similar transcription pattern. In a long-term (2-week) fasting and refeeding experiment, the hepatic lepA and brain lepB showed a similar transcription change pattern induced by food deprivation stimulation but differential changes after refeeding. These findings suggest snakehead lepA and lepB are differential both in tissue distribution and molecular functions, and they might play as an important regulator in energy metabolism and food intake in fish, respectively.

Characterization of two kcnk3 genes in Nile tilapia (Oreochromis niloticus): Molecular cloning, tissue distribution, and transcriptional changes in various salinity of seawater.

As one important member of the two-pore-domain potassium channel (K2P) family, potassium channel subfamily K member 3 (KCNK3) has been reported for thermogenesis regulation, energy homeostasis, membrane potential conduction, and pulmonary hypertension in mammals. However, its roles in fishes are far less examined and published. In the present study, we identified two kcnk3 genes (kcnk3a and kcnk3b) in an euryhaline fish, Nile tilapia (Oreochromis niloticus), by molecular cloning, genomic survey and laboratory experiments to investigate their potential roles for osmoregulation. We obtained full-length coding sequences of the kcnk3a and kcnk3b genes (1209 and 1173 bp), which encode 402 and 390 amino acids, respectively. Subsequent multiple sequence alignments, putative 3D-structure model prediction, genomic survey and phylogenetic analysis confirmed that two kcnk3 paralogs are widely presented in fish genomes. Interestingly, a DNA fragment inversion of a kcnk3a cluster was found in Cypriniforme in comparison with other fishes. Quantitative real-time PCRs demonstrated that both the tilapia kcnk3 genes were detected in all the examined tissues with a similar distribution pattern, and the highest transcriptions were observed in the heart. Meanwhile, both kcnk3 genes in the gill were proved to have a similar transcriptional change pattern in response to various salinity of seawater, implying that they might be involved in osmoregulation. Furthermore, three predicted transcription factors (arid3a, arid3b, and arid5a) of both kcnk3 genes also showed a similar pattern as their target genes in response to the various salinity, suggesting their potential positive regulatory roles. In summary, we for the first time characterized the two kcnk3 genes in Nile tilapia, and demonstrated their potential involvement in osmoregulation for this economically important fish.

Genome Sequencing of the Japanese Eel (Anguilla japonica) for Comparative Genomic Studies on tbx4 and a tbx4 Gene Cluster in Teleost Fishes.

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

Molecular cloning of two kcnk3 genes from the Northern snakehead (Channa argus) for quantification of their transcriptions in response to fasting and refeeding.

Potassium channel subfamily K member 3 (KCNK3) has been reported to play important roles in membrane potential conduction, pulmonary hypertension and thermogenesis regulation in mammals. However, its roles remain largely unknown and scarce reports were seen in fish. In the present study, we for the first time identified two kcnk3 genes (kcnk3a and kcnk3b) from the carnivorous Northern snakehead (Channa argus) by molecular cloning and a genomic survey. Subsequently, their transcription changes in response to different feeding status were investigated. Full-length coding sequences of the kcnk3a and kcnk3b genes are 1203 and 1176 bp, encoding 400 and 391 amino acids, respectively. Multiple alignments, 3D-structure prediction and phylogenetic analysis further suggested that these kcnk3 genes may be highly conserved in vertebrates. Tissue distribution analysis by real-time PCR demonstrated that both the snakehead kcnk3s were widely transcribed in majority of the examined tissues but with different distribution patterns. In a short-term (24-h) fasting experiment, we observed that brain kcnk3a and kcnk3b genes showed totally opposite transcription patterns. In a long-term (2-week) fasting and refeeding experiment, we also observed differential change patterns for the brain kcnk3 genes. In summary, our findings suggest that the two kcnk3 genes are close while present different transcription responses to fasting and refeeding. They therefore can be potentially selected as novel target genes for improvement of production and quality of this economically important fish.

Complete mitochondrial genome of a kind of snakehead fish Channa siamensis and its phylogenetic consideration.

The snakehead fish, Channa siamensis, belongs to the genus of Channa (perciformes: Channidae) and was first reported by Günther in 1861. Despite it has been described approximately for 15 decades, the genetic information is limited and the taxon status of this kind of fish is still unclear. The primary objective of this study is to get more genomic data and calculate the taxon location of this kind of fish. The next generation sequencing method was used to obtain the whole mitochondrial DNA information, and bioinformatic analysis was performed to investigate the evolutionary status and taxon location of C. siamensis. The circular mitochondrial DNA was 16,570 bp in length, and which showed typical piscine structure and arrangement. The overall nucleotide composition was 29.28% A, 24.72% T, 30.71% C, 15.29% G, with 54.1% AT, respectively. Phylogenetic analyses using concatenated amino acid and nucleotide sequences of the 13 protein-coding genes with two different methods (Maximum likelihood and Bayesian analysis) both highly supported C. siamensis belongs to the genus Channa and shows a close relationship with C. micropeltes. These data will provide more useful information for a better understanding of the mitochondrial genomic diversities and evolution in fish as well as novel genetic markers for studying population genetics and species identification.

Complete mitochondrial genome of the snakehead (Channa gachua) and its phylogeny.

In present study, the mitochondrial genome (mitogenome) of Channa gachua was determined and the phylogenetic relationship of Channidae fish was reconsidered. The mitogenome of the C. gachua is 16547 bp in length, containing 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosome RNA genes (rRNAs), a control region (D-loop) and an origin region of replication on the light-strand (OL). The overall nucleotide composition is 28.32% A, 26.58% T, 29.41% C, 15.69% G, with 54.90% AT, respectively. Phylogenetic analyses revealed that C. gachua belongs to the genus Channa and shares a close relationship with C. marulius and C. striata.