The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates.

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.

A comparative transcriptomic study on developmental gonads provides novel insights into sex change in the protandrous black porgy (Acanthopagrus schlegelii).

Protandrous black porgy (Acanthopagrus schlegelii) is a popular and valuable commercial marine fish in China and East Asian countries. Controlling and managing its breeding has been an imperative step towards obtaining a sustainable supply of this fish in aquaculture production systems. Therefore, study on the molecular mechanisms of sex change in black porgy has both scientific and commercial importance. Previously, we identified some candidate genes related to sex determination and differentiation from a high-quality genome assembly of the black porgy. In the present study, transcriptome sequencing of developmental gonads (including testis, ovotestis and ovary) of black porgy was performed to further investigate the sex-change mechanisms. Our results showed that the highly expressed male-related genes (dmrt1, piwi1, piwi2, sox9, sox30 and amh) at the male phase were significantly down-regulated to a substantial degree at the intersexual stage, and the female-related genes (jnk1, vasa, wnt4, figla and foxl2) were distinctly up-regulated when the fish grows into a female adult, suggesting the potential roles of these genes in sex change of the black porgy. These data also support a previous hypothesis that the femaleness will be switched on when the testis is entering the degenerated stage through the diminished dmrt1 expression. Our transcriptome data provide a very useful genomic resource for future studies on sex change and practical aquaculture in the black porgy.

Genome and Transcriptome Sequencing of casper and roy Zebrafish Mutants Provides Novel Genetic Clues for Iridophore Loss.

casper has been a widely used transparent mutant of zebrafish. It possesses a combined loss of reflective iridophores and light-absorbing melanophores, which gives rise to its almost transparent trunk throughout larval and adult stages. Nevertheless, genomic causal mutations of this transparent phenotype are poorly defined. To identify the potential genetic basis of this fascinating morphological phenotype, we constructed genome maps by performing genome sequencing of 28 zebrafish individuals including wild-type AB strain, roy orbison (roy), and casper mutants. A total of 4.3 million high-quality and high-confidence homozygous single nucleotide polymorphisms (SNPs) were detected in the present study. We also identified a 6.0-Mb linkage disequilibrium block specifically in both roy and casper that was composed of 39 functional genes, of which the mpv17 gene was potentially involved in the regulation of iridophore formation and maintenance. This is the first report of high-confidence genomic mutations in the mpv17 gene of roy and casper that potentially leads to defective splicing as one major molecular clue for the iridophore loss. Additionally, comparative transcriptomic analyses of skin tissues from the AB, roy and casper groups revealed detailed transcriptional changes of several core genes that may be involved in melanophore and iridophore degeneration. In summary, our updated genome and transcriptome sequencing of the casper and roy mutants provides novel genetic clues for the iridophore loss. These new genomic variation maps will offer a solid genetic basis for expanding the zebrafish mutant database and in-depth investigation into pigmentation of animals.

Transcriptomic evidence of adaptive tolerance to high environmental ammonia in mudskippers.

Mudskippers are typical amphibious fishes and possess various strategies to ameliorate ammonia toxicity during exposure to environmental ammonia. The present study aimed to provide transcriptomic evidence through profiling the gill and liver transcriptomes of Boleophthalmus pectinirostris (BP) and Periophthalmus magnuspinnatus (PM), which were subjected to treatment with high environmental ammonia for up to 72 h. The results of gene function annotation showed that most of the differentially expressed genes were involved in metabolic pathways. After ammonia exposure, the protein and amino acid metabolism related genes in mudskippers were down-regulated, and PM had more down-regulated genes than BP. The expression levels of several representative genes involved in ammonia excretion in the gill were commonly increased. Interestingly, NH4+ transporting and H+ excreting related genes, including Na+/K+(NH4+)/2Cl- cotransporter (nkcc), Na+/K+(NH4+)-ATPase (nka), carbonic anhydrase 2 (ca2), H+-ATPase, Na+/H+ (NH4+)-exchanger (nhe), and carbonic anhydrase 15 (ca15), were up-regulated more significantly in BP than PM; however, the transcription levels of Rhesus glucoprotein b (Rhbg) and Rhesus glucoprotein c1 (Rhcg1), which constitute the NH3 transporting channels, were up-regulated more significantly in PM than BP. Furthermore, the present study provides molecular evidence for how mudskippers adopt partial amino acid catabolism to decrease the production of endogenous ammonia under high environmental ammonia loading.

The Sinocyclocheilus cavefish genome provides insights into cave adaptation.

BACKGROUND: An emerging cavefish model, the cyprinid genus Sinocyclocheilus, is endemic to the massive southwestern karst area adjacent to the Qinghai-Tibetan Plateau of China. In order to understand whether orogeny influenced the evolution of these species, and how genomes change under isolation, especially in subterranean habitats, we performed whole-genome sequencing and comparative analyses of three species in this genus, S. grahami, S. rhinocerous and S. anshuiensis. These species are surface-dwelling, semi-cave-dwelling and cave-restricted, respectively. RESULTS: The assembled genome sizes of S. grahami, S. rhinocerous and S. anshuiensis are 1.75 Gb, 1.73 Gb and 1.68 Gb, respectively. Divergence time and population history analyses of these species reveal that their speciation and population dynamics are correlated with the different stages of uplifting of the Qinghai-Tibetan Plateau. We carried out comparative analyses of these genomes and found that many genetic changes, such as gene loss (e.g. opsin genes), pseudogenes (e.g. crystallin genes), mutations (e.g. melanogenesis-related genes), deletions (e.g. scale-related genes) and down-regulation (e.g. circadian rhythm pathway genes), are possibly associated with the regressive features (such as eye degeneration, albinism, rudimentary scales and lack of circadian rhythms), and that some gene expansion (e.g. taste-related transcription factor gene) may point to the constructive features (such as enhanced taste buds) which evolved in these cave fishes. CONCLUSION: As the first report on cavefish genomes among distinct species in Sinocyclocheilus, our work provides not only insights into genetic mechanisms of cave adaptation, but also represents a fundamental resource for a better understanding of cavefish biology.

A Comparative Genomic and Transcriptomic Survey Provides Novel Insights into N-Acetylserotonin Methyltransferase (ASMT) in Fish.

Melatonin is a multifunctional bioactive molecule that plays comprehensive physiological roles in all living organisms. N-acetylserotonin methyltransferase (ASMT, also known as hydroxyindole O-methyltransferase or HIOMT) is the final enzyme for biosynthesis of melatonin. Here, we performed a comparative genomic and transcriptomic survey to explore the ASMT family in fish. Two ASMT isotypes (ASMT1 and ASMT2) and a new ASMT-like (ASMTL) are all extracted from teleost genomes on the basis of phylogenetic and synteny analyses. We confirmed that C-terminal of the ASMTL proteins (ASMTL-ASMT) is homology to the full length of ASMT1 and ASMT2. Our results also demonstrate that the two ASMT isotypes and their distribution in teleosts seem to be the result of combinations of whole-genome duplication (WGD) and gene loss. Differences were also observed in tissue distribution and relative transcript abundances of ASMT1, ASMT2 and ASMTL through transcriptomic analysis. Protein sequence alignment and 3D structure prediction of ASMTs and ASMTL suggest differential roles for these ASMT genes. In summary, our current work provides novel insights into the ASMT genes in fish by combination of genomic and transcriptomic data.

Prediction of Toxin Genes from Chinese Yellow Catfish Based on Transcriptomic and Proteomic Sequencing.

Fish venom remains a virtually untapped resource. There are so few fish toxin sequences for reference, which increases the difficulty to study toxins from venomous fish and to develop efficient and fast methods to dig out toxin genes or proteins. Here, we utilized Chinese yellow catfish (Pelteobagrus fulvidraco) as our research object, since it is a representative species in Siluriformes with its venom glands embedded in the pectoral and dorsal fins. In this study, we set up an in-house toxin database and a novel toxin-discovering protocol to dig out precise toxin genes by combination of transcriptomic and proteomic sequencing. Finally, we obtained 15 putative toxin proteins distributed in five groups, namely Veficolin, Ink toxin, Adamalysin, Za2G and CRISP toxin. It seems that we have developed a novel bioinformatics method, through which we could identify toxin proteins with high confidence. Meanwhile, these toxins can also be useful for comparative studies in other fish and development of potential drugs.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation.

INTRODUCTION: Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation. METHODS: Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation. RESULTS: We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation. CONCLUSION: These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Genome-Wide Identification, Sequence Alignment, and Transcription of Five Sex-Related Genes in Largemouth Bass (Micropterus Salmoides).

BACKGROUND: Largemouth bass (Micropterus Salmoides) is an economically important fish species in China. Most research has focused on its growth, disease resistance, and nutrition improvement. However, the sex-determining genes in largemouth bass are still unclear. The transforming growth factor-beta (TGF-ß) gene family, including amh, amhr2 and gsdf, plays an important role in the sex determination and differentiation of various fishes. These genes are potentially involved in sex determination in largemouth bass. METHODS: We performed a systematic analysis of 5 sex-related genes (amh, amhr2, gsdf, cyp19a1, foxl2) in largemouth bass using sequence alignment, collinearity analysis, transcriptome, and quantitative real-time polymerase chain reaction (qRT-PCR). This included a detailed assessment of their sequences, gene structures, evolutionary traits, and gene transcription patterns in various tissues including gonads, and at different developmental stages. RESULTS: Comparative genomics revealed that the 5 sex-related genes were highly conserved in various fish genomes. These genes did not replicate, mutate or lose in largemouth bass. However, some were duplicated (amh, amhr2 and gsdf), mutated (gsdf) or lost (amhr2) in other fishes. Some genes (e.g., gsdf) showed significant differences in genomic sequence between males and females, which may contribute to sex determination and sex differentiation in these fishes. qRT-PCR was applied to quantify transcription profiling of the 5 genes during gonadal development and in the adult largemouth bass. Interestingly, amh, amhr2 and gsdf were predominantly expressed in the testis, while cyp19a1 and foxl2 were mainly transcribed in the ovary. All 5 sex-related genes were differentially expressed in the testes and ovaries from the 56th day post-fertilization (dpf). We therefore speculate that male/female differentiation in the largemouth bass may begin at this critical time-point. Examination of the transcriptome data also allowed us to screen out several more sex-related candidate genes. CONCLUSIONS: Our results provide a valuable genetic resource for investigating the physiological functions of these 5 sex-related genes in sex determination and gonadal differentiation, as well as in the control of gonad stability in adult largemouth bass.

Effect of aroma-producing yeasts in high-salt liquid-state fermentation soy sauce and the biosynthesis pathways of the dominant esters.

Fermentation with excellent aroma-producing yeasts can enhance the flavour of soy sauce. In this work, Millerozyma farinosa CS2.23, Zygosaccharomyces rouxii CS2.42, and Candida parapsilosis CS2.53 were added to the high-salt liquid-state moromi to promote soy sauce fermentation. All three yeasts improved the TE of soy sauce, the highest of which reached 1.03 g/L with added CS2.42. Other quality indexes of soy sauce, including RS, TA, and AN, were not greatly affected. The volatile esters of soy sauce added to the three yeasts increased by 108.85%, 166.71%, and 113.61% compared with the control through GC-MS analysis. Obviously, CS2.42 had an excellent ability to produce esters. Studying the biosynthesis pathway of esters, CS2.42 has the best esterification ability, while CS2.53 has the advantage of alcoholysis ability. The exploration of the biosynthetic pathway of acetate and ethyl esters has laid a foundation for regulating esters in soy sauce fermentation.

The first Conus genome assembly reveals a primary genetic central dogma of conopeptides in C. betulinus.

Although there are various Conus species with publicly available transcriptome and proteome data, no genome assembly has been reported yet. Here, using Chinese tubular cone snail (C. betulinus) as a representative, we sequenced and assembled the first Conus genome with original identification of 133 genome-widely distributed conopeptide genes. After integration of our genomics, transcriptomics, and peptidomics data in the same species, we established a primary genetic central dogma of diverse conopeptides, assuming a rough number ratio of ~1:1:1:10s for the total genes: transcripts: proteins: post-translationally modified peptides. This ratio may be special for this worm-hunting Conus species, due to the high diversity of various Conus genomes and the big number ranges of conopeptide genes, transcripts, and peptides in previous reports of diverse Conus species. Only a fraction (45.9%) of the identified conotopeptide genes from our achieved genome assembly are transcribed with transcriptomic evidence, and few genes individually correspond to multiple transcripts possibly due to intraspecies or mutation-based variances. Variable peptide processing at the proteomic level, generating a big diversity of venom conopeptides with alternative cleavage sites, post-translational modifications, and N-/C-terminal truncations, may explain how the 133 genes and ~123 transcripts can generate thousands of conopeptides in the venom of individual C. betulinus. We also predicted many conopeptides with high stereostructural similarities to the putative analgesic ω-MVIIA, addiction therapy AuIB and insecticide ImI, suggesting that our current genome assembly for C. betulinus is a valuable genetic resource for high-throughput prediction and development of potential pharmaceuticals.

Comprehensive transcriptional changes in the liver of Kanglang white minnow (Anabarilius grahami) in response to the infection of parasite Ichthyophthirius multifiliis.

Abstract: The notorious parasite Ichthyophthirius multifiliis (Ich) has been recorded worldwide in fish species and causes white spot disease, posing major threats and resulting in severe losses to international fish production. Extensively effective strategies for treating Ich are not available yet, and genetic mechanisms of hosts in response to the parasite are still largely unknown. In this study, we selected Kanglang white minnow (KWM, Anabarilius grahami) to examine its liver transcriptional changes after Ich infection, as white spot disease is one bottleneck problem in exploring this economically important species. We divided the experimental fishes into three groups (control, early-infected, and late-infected) to examine differentially expressed genes (DEGs). A total of 831 DEGs were identified and classified into 128 significantly enriched GO (Gene Ontology) terms and 71 significantly enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Most of these terms or pathways were functionally enriched in immunity, inflammatory response, and apoptosis, such as nucleotide-binding oligomerization domain-like (NOD-like) receptor signaling, tumor necrosis factor (TNF) signaling, interleukin-17 (IL-17) signaling, and apoptosis pathways. We also identified 178 putative antimicrobial peptides (AMPs) and AMP precursors based on our previously reported genome assembly of KWM, and revealed that the expressional patterns varied according to different types. In summary, our work reported the first comprehensive transcriptional changes in KWM in response to the exogenous infection of Ich, which would lay a solid foundation for in-depth studies on disease defense or resistant strains selection in this valuable fish.

High-throughput identification of heavy metal binding proteins from the byssus of chinese green mussel (Perna viridis) by combination of transcriptome and proteome sequencing.

The Byssus, which is derived from the foot gland of mussels, has been proved to bind heavy metals effectively, but few studies have focused on the molecular mechanisms behind the accumulation of heavy metals by the byssus. In this study, we integrated high-throughput transcriptome and proteome sequencing to construct a comprehensive protein database for the byssus of Chinese green mussel (Perna viridis), aiming at providing novel insights into the molecular mechanisms by which the byssus binds to heavy metals. Illumina transcriptome sequencing generated a total of 55,670,668 reads. After filtration, we obtained 53,047,718 clean reads and subjected them to de novo assembly using Trinity software. Finally, we annotated 73,264 unigenes and predicted a total of 34,298 protein coding sequences. Moreover, byssal samples were analyzed by proteome sequencing, with the translated protein database from the foot transcriptome as the reference for further prediction of byssal proteins. We eventually determined 187 protein sequences in the byssus, of which 181 proteins are reported for the first time. Interestingly, we observed that many of these byssal proteins are rich in histidine or cysteine residues, which may contribute to the byssal accumulation of heavy metals. Finally, we picked one representative protein, Pvfp-5-1, for recombinant protein synthesis and experimental verification of its efficient binding to cadmium (Cd2+) ions.

The Complete Mitochondrial Genome of Glyptothorax macromaculatus Provides a Well-Resolved Molecular Phylogeny of the Chinese Sisorid Catfishes.

Previous phylogenetic analyses of the Chinese sisorid catfishes have either been poorly resolved or have not included all the 12 sisorid genera. Here, we successfully assembled the first complete mitochondrial genome of the sisorid fish Glyptothorax macromaculatus. Based on this novel mitochondrial genome and previously published mitochondrial genomes in the Sisoridae, we generated maximum likelihood and Bayesian phylogenies. We dated our preferred topology using fossil calibration points. We also tested the protein-coding genes in the mitochondrial genomes of the glyptosternoid fishes for signals of natural selection by comparing the nucleotide substitution rate along the branch ancestral to the glyptosternoid fishes to other branches in our topology. The mitochondrial sequence structure of G. macromaculatus was similar to those known from other vertebrates, with some slight differences. Our sisorid phylogenies were well-resolved and well-supported, with exact congruence between the different phylogenetic methods. This robust phylogeny clarified the relationships among the Chinese sisorid genera and strongly supported the division of the family into three main clades. Interestingly, the glyptosternoid divergence time predicted by our molecular dating analysis coincided with the uplift of the Tibetan Plateau, suggesting that geology may have influenced speciation in the Sisoridae. Among the mitochondrial protein-coding genes, atp8 may have most rapidly evolved, and atp6 may have been subjected to positive selection pressure to adapt to high elevations. In summary, this study provided novel insights into the phylogeny, evolution and high-altitude adaptions of the Chinese sisorid fishes.

A new species of bandy-bandy (Vermicella: Serpentes: Elapidae) from the Weipa region, Cape York, Australia.

Bandy-bandies (genus Vermicella) are small (50-100cm) black and white burrowing elapids with a highly specialised diet of blindsnakes (Typhlopidae). There are currently 5 recognized species in the genus, all located in Australia, with Vermicella annulata the most encountered species with the largest distribution. Morphological and mitochondrial analyses of specimens collected from the Weipa area, Cape York, Queensland reveal the existence of a new species, which we describe as Vermicella parscauda sp. nov. Mitochondrial DNA analysis (16S and ND4) and external morphological characteristics indicate that the closest relatives of the new species are not V. annulata, which also occurs on Cape York, but rather species from Western Australia and the Northern Territory (V. intermedia and V. multifasciata) which, like V. parscauda, occupy monsoon habitats. Internasal scales are present in V. parscauda sp. nov., similar to V. annulata, but V. intermedia and V. multifasciata do not have nasal scales. V. parscauda sp. nov. has 55-94 black dorsal bands and mottled or black ventral scales terminating approximately 2/3rds of the body into formed black rings, suggesting that hyper-banding is a characteristic of the tropical monsoon snakes (V. intermedia, V. multifasciata and V. parscauda). The confined locality, potential habitat disruption due to mining activities, and scarcity of specimens indicates an urgent conservation concern for this species.

Assessing the genetic diversity of the critically endangered Chinese sturgeon Acipenser sinensis using mitochondrial markers and genome-wide single-nucleotide polymorphisms from RAD-seq.

As a living fossil, the endangered Chinese sturgeon (Acipenser sinensis) has been considered a national treasure in China. Here, the famous Gezhouba Dam and Three Gorges Dam on the Yangtze River were built in 1988 and 2006, for economic purposes. The natural population of Chinese sturgeon has declined since then, as these dams block its migratory route to the original spawning grounds in the middle reaches of the Yangtze River. In 2013 and 2014, there was an absence of spawning where it typically happened near the Gezhouba Dam. Nevertheless, from April to June in 2015, over 1,000 larvae with different body lengths (10-35 cm) were detected along the Shanghai Yangtze Estuary; but only little is currently known about the population genetic structure of the Chinese sturgeon. Herein, we inferred population genetic parameters from 462 available Chinese sturgeon specimens based on a 421-bp fragment of the mitochondrial DNA (mtDNA) D-loop region and 1,481,620 SNPs (single-nucleotide polymorphisms) generated by restriction site-associated DNA sequencing (RAD-seq). For the D-loop dataset, 15 haplotypes were determined. Randomly picked 23 individuals, representing the 15 D-loop haplotype groups, were subsequently used for further RAD-seq validation. The average nucleotide diversity calculated from the mtDNA and RAD datasets was 0.0086 and 0.000478, respectively. The overall effective female population size was calculated to be 1,255 to 2,607, and the long-term effective population size was estimated to range from 11,950 to 119,500. We observed that the genetic variability and the effective female population size of the current population in the Yangtze River are severely low, which are similar to the data reported over 10 years ago. The deduced relatively small effective population of female fish, limiting the genetic connectivity among Chinese sturgeon, should be considered a serious threat to this endangered species.

De novo assembly and comparative transcriptome analysis of the foot from Chinese green mussel (Perna viridis) in response to cadmium stimulation.

The Chinese green mussel, Perna viridis, is a marine bivalve with important economic values as well as biomonitoring roles for aquatic pollution. Byssus, secreted by the foot gland, has been proved to bind heavy metals effectively. In this study, using the RNA sequencing technology, we performed comparative transcriptomic analysis on the mussel feet with or without inducing by cadmium (Cd). Our current work is aiming at providing insights into the molecular mechanisms of byssus binding to heavy metal ions. The transcriptome sequencing generated a total of 26.13-Gb raw data. After a careful assembly of clean data, we obtained a primary set of 105,127 unigenes, in which 32,268 unigenes were annotated. Based on the expression profiles, we identified 9,048 differentially expressed genes (DEGs) between Cd treatment (50 or 100 µg/L) at 48 h and the control, suggesting an extensive transcriptome response of the mussels during the Cd stimulation. Moreover, we observed that the expression levels of 54 byssus protein coding genes increased significantly after the 48-h Cd stimulation. In addition, 16 critical byssus protein coding genes were picked for profiling by quantitative real-time PCR (qRT-PCR). Finally, we reached a primary conclusion that high content of tyrosine (Tyr), cysteine (Cys), histidine (His) residues or the special motif plays an important role in the accumulation of heavy metals in byssus. We also proposed an interesting model for the confirmed byssal Cd accumulation, in which biosynthesis of byssus proteins may play simultaneously critical roles since their transcription levels were significantly elevated.

The complete mitochondrial genome of horn-nosed boxfish (Ostracion rhinorhynchos).

The horn-nosed boxfish, Ostracion rhinorhynchos (Tetraodontiformes: Ostraciidae) is a toxic marine species inhabiting tropical coral reefs. In this study, we first reported its whole mitochondrial genome sequence. The complete mitochondrial genome, 16 483 bp with an AT ratio of 56.8%, is composed of 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs and an 826-bp D-loop control region. The molecular-based phylogenetic tree indicated that O. rhinorhynchos has close affinities with fishes from family Ostraciidae as expected.

The complete mitochondrial genome of Indonesian snakehead, Channa micropeltes (Channiformes, Channidae).

Here we have characterized the complete mitochondrial genome of Indonesian snakehead, Channa micropeltes, and described its organization in this paper. The complete mitochondrial genome, 16,567 bp, is composed of 32.54% A, 14.05% G, 25.82% T and 27.59% C. It includes 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a 921-bp D-loop control region. Phylogenetic analysis with five more fish species demonstrated that the Indonesian snakehead is most closely related to the Great snakehead (Channa marulius). Our mitochondrial genomic data will also be valuable for the study on the mitochondrial evolution of fishes.

The complete mitochondrial genome sequence of the giant mudskipper, Periophthalmodon schlosseri (Perciformes: gobiidae).

The complete mitochondrial genome sequence of the giant mudskipper, Periophthalmodon schlosseri, was first reported in this study. The circle genome is 16 470 bp in length, with the base composition of 28.4% T, 15.4% G, 29.2% A and 27.0% C. The mitogenome includes 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and one D-loop region. Only the NADH dehydrogenase subunit 6 (NAD6) and eight tRNA genes are encoded on the light strand. The mitochondrial gene arrangement of P. schlosseri is similar to those of most other gobies. The phylogenetic analysis using the neighbor-joining method showed that the kinship between Periophthalmodon and Periophthalmus is closer than those between Periophthalmodon and other selected genera. Our complete mitogenome data are going to provide the basis for taxonomic and phylogenetic research of amphibious mudskippers.