The complete mitochondrial genome of the Chinese giant salamander, Andrias davidianus (Amphibia: Caudata)

Discovery of ex situ individuals of Andrias sligoi, an extremely endangered species and one of the largest amphibians worldwide

The South China giant salamander, Andrias sligoi, is one of the largest extant amphibian species worldwide. It was recently distinguished from another Chinese species, the Chinese giant salamander, Andrias davidianus, which is considered Critically Endangered according to the International Union for Conservation of Nature (IUCN) Red List. It appears too late to save this extremely rare and large amphibian in situ. Another extant species of the same genus, Andrias japonicus, inhabits Japan. However, the introduction of Chinese giant salamanders into some areas of Japan has resulted in hybridization between the Japanese and Chinese species. During our genetic screening of giant salamanders in Japan, we unexpectedly discovered four individuals of the South China giant salamander: two were adult males in captivity, and one had recently died. The last individual was a preserved specimen. In this study, we report these extremely rare individuals of A. sligoi in Japan and discuss the taxonomic and conservational implications of these introduced individuals.

Molecular cloning, tissues distribution, and function analysis of thioredoxin-like protein-1 (TXNL1) in Chinese giant salamanders Andrias davidianus

Thioredoxin-like protein-1 (TXNL1) is the member of thioredoxin superfamily, a family of thiol oxidoreductases. TXNL1 plays an important role in scavenging ROS and the maintenance of cellular redox balance. However, its physiological functions in Andrias davidianus have not been well understood. In the present study, the full-length cDNA encoding thioredoxin-like protein-1 (AdTXNL1) of A. davidianus was cloned, the mRNA tissue distribution was analyzed, and the function was characterized. The Adtxnl1 cDNA contained an open reading frame (ORF) of 870 bp encoding a polypeptide of 289 amino acids with the N-terminal TRX domain, a Cys³⁴-Ala³⁵-Pro³⁶-Cys³⁷ (CAPC) motif, and the C-terminal proteasome-interacting thioredoxin domain (PITH). The mRNA of AdTXNL1 was expressed in a wide range of tissues, with the highest level in the liver. The transcript level of AdTXNL1 was significantly up-regulated post Aeromonas hydrophila challenge in liver tissue. Moreover, the recombinant AdTXNL1 protein was produced and purified, and used to investigate the antioxidant activity. In the insulin disulfide reduction assay, rAdTXNL1 exhibited strong antioxidant capability. Altogether, the thioredoxin-like protein-1 may be involved in reduction/oxidation (redox) balance and as an important immunological gene in A. davidianus.

Changes of intestinal microbiota in the giant salamander (Andrias davidianus) during growth based on high-throughput sequencing

Despite an increasing appreciation of the importance of host–microbe interaction in healthy growth, information on gut microbiota changes of the Chinese giant salamander (Andrias davidianus) during growth is still lacking. Moreover, it is interesting to identify gut microbial structure for further monitoring A. davidianus health. This study explored the composition and functional characteristics of gut bacteria in different growth periods, including tadpole stage (ADT), gills internalization stage (ADG), 1 year age (ADY), 2 year age (ADE), and 3 year age (ADS), using high-throughput sequencing. The results showed that significant differences were observed in microbial community composition and abundance among different growth groups. The diversity and abundance of intestinal flora gradually reduced from larvae to adult stages. Overall, the gut microbial communities were mainly composed of Fusobacteriota, Firmicutes, Bacteroidota, and Proteobacteria. More specifically, the Cetobacterium genus was the most dominant, followed by Lactobacillus and Candidatus Amphibiichlamydia. Interestingly, Candidatus Amphibiichlamydia, a special species related to amphibian diseases, could be a promising indicator for healthy monitoring during A. davidianus growth. These results could be an important reference for future research on the relationship between the host and microbiota and also provide basic data for the artificial feeding of A. davidianus.

Using the Ratio of Urine Testosterone to Estrone-3-Glucuronide to Identify the Sex of Chinese Giant Salamanders (Andrias davidianus)

Minimally invasive sampling was used to determine the sex of Chinese giant salamanders (Andrias davidianus). Urine samples (n = 25) were collected from 6 adults in the breeding season and from 19 individuals (7 adults and 12 juveniles) in the non-breeding season. The hormone testosterone (T) and estrone-3-glucuronide (E1G) in urine were collected from Chinese giant salamanders (CGSs), and the hormone extracts were analyzed by enzyme immunoassays (EIA). The data demonstrated that the urine T concentration of the male CGSs was significantly higher than that of the females during the breeding season (p < 0.05) and even more pronounced during the non-breeding season (p < 0.01). The urine E1G concentration of the males was less pronounced than that of the females during the breeding season (p < 0.01) and significantly lower during the non-breeding season (p < 0.05). The urine T/E1G values of all the male salamanders were significantly higher than those of the females (p < 0.01) during both the breeding season and the non-breeding season. An interesting pattern was found in this study: the value of urine log10(T/E1G) of the male CGSs was higher than 1, whereas the value for the females was lower than 1, during both the breeding and non-breeding seasons, and in the adult and sub-adult age groups of CGSs. There were 25 salamanders in this study and the accuracy rate reached 100% by using a log10(T/E1G) value of 1. The results of the log10(T/E1G) value provide new insight into the future development of the sex identification of CGSs and also lay the foundation for accurate sex identification in the preparation for artificial release. This is the first study to show that the T/E1G ratio in urinary hormones is reliable for the sex identification of CGSs. Additionally, urinary hormone T/E1G measures are promising sex identification tools for amphibian or monomorphic species and for those whose secondary sex characteristics are visible only during the breeding season.

Development of an eDNA metabarcoding tool for surveying the world’s largest amphibian

Due to the overexploitation of farming, as well as habitat destruction, the wild population of Chinese giant salamander (CGS) Andrias davidianus, a species with seven genetically distinct lineages, has decreased by over 80% in the past 70 years. Traditional survey methods have proven to be unsuitable for finding this rare and elusive species. We evaluated the efficacy of environmental DNA (eDNA) sampling to detect CGS indirectly from its aquatic environment. We developed several species-specific primer sets; validated their specificity and sensitivity; and assessed their utility in silico, in the laboratory, and at two field sites harboring released farm-bred CGS. We detected the presence of CGS DNA by using polymerase chain reaction and Sanger sequencing. We also sequenced an amplicon mixture of seven haplotype-represented samples using high-throughput sequencing. Our eDNA methods could detect the presence of CGS at moderate densities reported across its range, proving them as a cost-effective way to establish broad-scale patterns of occupancy for CGS. In addition, our primers enabled the detection of mitochondrial lineage mixture or introduced individuals from geographically isolated populations of CGS.

The Immune System and the Antiviral Responses in Chinese Giant Salamander, Andrias davidianus

The Chinese giant salamander, belonging to an ancient amphibian lineage, is the largest amphibian existing in the world, and is also an important animal for artificial cultivation in China. However, some aspects of the innate and adaptive immune system of the Chinese giant salamander are still unknown. The Chinese giant salamander iridovirus (GSIV), a member of the Ranavirus genus (family Iridoviridae), is a prominent pathogen causing high mortality and severe economic losses in Chinese giant salamander aquaculture. As a serious threat to amphibians worldwide, the etiology of ranaviruses has been mainly studied in model organisms, such as the Ambystoma tigrinum and Xenopus. Nevertheless, the immunity to ranavirus in Chinese giant salamander is distinct from other amphibians and less known. We review the unique immune system and antiviral responses of the Chinese giant salamander, in order to establish effective management of virus disease in Chinese giant salamander artificial cultivation.

Identification and functional analysis of two interferon regulatory factor 3 genes and their involvement in antiviral immune responses in the Chinese giant salamander Andrias davidianus

Interferon regulatory factor 3 (IRF3), a crucial member of interferon regulatory factor (IRF) family, plays an important role in innate immunity in vertebrates. However, there are no reports on the characterization and especially their respective functional analysis of two IRF3 genes in some species. In this study, two IRF3 genes as well as their roles in the immune response were identified and investigated in Chinese giant salamander, Andrias davidianus. The complete open reading frames of AdIRF3A and AdIRF3B were 1, 113 bp and 1, 380 bp in length, encoding 370 and 459 amino acids, respectively. Both AdIRF3A and AdIRF3B protein contain an IRF and an IRF3 domain. Phylogenetic analysis indicated that AdIRF3s clustered together with other IRF3 proteins. Tissue distribution analysis showed that AdIRF3s were expressed in all tissues tested, with highest expression levels in blood. Both AdIRF3s actively responded to Chinese giant salamander iridovirus (GSIV) and poly (I:C) challenge in A. davidianus. AdIRF3A/B silencing significantly suppressed the DNA virus and viral RNA analog-induced expression of IFN-inducible genes. Luciferase reporter assay further confirmed the regulatory role of AdIRF3s in IFN signaling. These results provide new insights into the origin or evolution of IRF3 in amphibians and even in vertebrates.

The first evidence of four transcripts from two Interleukin 18 genes in animal and their involvement in immune responses in the largest amphibian Andrias davidianus

Interleukin 18 (IL-18), a member of IL-1 cytokine superfamily, is an important proinflammatory cytokine with multiple functions in both innate immunity and acquired immunity. However, the characteristics and functional roles of IL-18 remain largely unknown in amphibians, which were classed as major group of vertebrates. In the present study, two IL-18 genes (AdIL-18A and AdIL-18B) and four transcripts (AdIL-18A1, AdIL-18A2, AdIL-18B1 and AdIL-18B2) were firstly identified and characterized from Chinese giant salamander (Andrias davidianus). To the best of our knowledge, this is the first report on the presence of more than one gene copy or two transcripts of IL-18 in one species. The complete open reading frames of AdIL-18A1, AdIL-18A2, AdIL-18B1 and AdIL-18B2 were 588 bp, 603 bp, 591 bp and 606 bp, respectively. The putative AdIL-18 proteins possessed the typical IL-1 domains and phylogenetic analysis indicated that AdIL-18s grouped together with other vertebrate IL-18 proteins. The expression profiles of AdIL-18s were investigated under the challenges of Aeromonas hydrophila, Staphylococcus ureae and Poly (I:C) respectively, and the results suggested that AdIL-18s were involved in the immune responses against both bacterial and viral infections. Moreover, the expression levels of two NF-κBs (P100 and P105) and four proinflammatory cytokines (IL-1β, IL-6, TNF-α and IFN-γ) were inhibited in AdIL-18A1/A2-silenced cells when treated with bacteria and viral RNA analog. Additionally, the transcription levels of these immune-related cytokine genes were markedly induced when the lymphocytes were treated with recombinant AdIL-18A1 or AdIL-18A2 proteins, implying the involvement of AdIL-18s in triggering NF-κB signaling and proinflammatory responses. These results might provide new insights into the origin or evolution of IL-18 in amphibians and even in vertebrates.

Draft genome sequencing and annotation of a low-virulence Morganella morganii strain CQ-M7, a multidrug-resistant isolate from the giant salamander in China

OBJECTIVES

A multidrug-resistant Morganella morganii strain (CQ-M7), isolated from the kidney of a diseased Chinese giant salamander in China, was examined with whole genome sequencing to better understand drug tolerance and its pathogenicity.

METHODS

The draft genome of the investigated strain was assembled using HGA assembler and annotated using Rapid Annotations Subsystems Technology (RAST) server. The contigs were annotated by the appropriate bioinformatics tools available on the National Center for Biotechnology Information (NCBI) website. Antibiotic resistance genes were detected by PCR. Pathogenicity of the isolate was performed on 30 healthy Chinese giant salamanders with different infection dosages.

RESULTS

The CQ-M7 strain showed resistance to multiple antimicrobials, especially to aminoglycoside and β-lactam antibiotics. Seventeen drug-resistance genes were detected, which were related to β-lactams, aminoglycosides, fluoroquinolones, tetracyclines, peptide antibiotic, and fosfomycin resistance. Sequence analysis showed the assembled genome size to be 4 966 326bp with 51.16% of GC content, containing 4587 protein-coding genes, 71 pseudogenes, five rRNAs, 80 tRNAs, and five noncoding RNAs. The genome sequence was deposited in GenBank under accession number RQIJ00000000. Artificial infection results indicated that the CQ-M7 strain was a low-virulence strain for the Chinese giant salamander.

CONCLUSION

It is believed that this is the first draft genome of Chinese giant salamander original Morganella morganii strain harbouring multiple antibiotic resistance genes in China. The reported genome sequence could provide insights into antibiotic resistance mechanisms and control strategies of Morganella morganii.

Involvement of a short-type peptidoglycan recognition protein (PGRP) from Chinese giant salamanders Andrias davidianus in the immune response against bacterial infection

PGRPs (Peptidoglycan recognition proteins) could recognize peptidoglycan and play vital roles in innate immunity among different animals. Till present, the functions of PGRP have been studied in various animals, but few reports have studied the amphibian PGRPs. In the current research, a short type PGRP was identified from Chinese giant salamander and its involvement in the innate immunity was studied. The ORF of AdPGRP-SC2 cDNA was 573 bp, which encoded 190 amino acids, and contained a PGRP and an amidase_2 domain. The qPCR analysis revealed that AdPGRP-SC2 mRNA transcripts expressed in different tissues, with the highest expression level in muscle, intestine and spleen. Results of immune challenges with peptidoglycan (PGN) demonstrated that expression patterns of AdPGRP-SC2 were significantly up-regulated in erythrocyte and spleen at the early injection stage. The recombinant AdPGRP-SC2 protein was successfully produced and purified, and it could show binding affinity to different bacteria. In the presence of Zn²⁺, the rAdPGRP-SC2 could exhibit a broad PAMPs binding activities, strongly agglutinate bacteria and exhibit amidase enzyme activity. Collectively, these data indicate AdPGRP-SC2 could act as PRR to recognize the invading microorganisms and as the antimicrobial effectors during the innate immune response of A. davidianus.

Four complete mitochondrial genomes of living wild-type chinese giant salamander Andrias davidianus (Amphibia: Cryptobranchidae)

The Chinese giant salamander (CGS), Andrias davidianus (Amphibian, Caudata, Cryptobranchidae), is endemic to China. After overhunting in the 1990's, it is very difficult to find the CGS in the wild. Due to mating disorder, the captive breeding population is genetically confounded. The genetic backgrounds of all wild-release individuals in China are not explicit. Herein, we reported four living wild-type complete mitochondrial genomes of this species. The gene order and contents are identical to those found in typical vertebrates. Thirteen protein-coding genes (PCGs) of 7 A. davidianus (4 from this study, 3 retrieved from GenBank) and 11 other closely species retrieved from GenBank were used to reconstruct phylogenetic tree. The Maximum likelihood (ML) topology shown that the clade of CGS has two subclades with a high support (100%). This study provides partial fundamental information for further exploring the true genetic background of whole population of A. davidianus.

Rag1 and rag2 gene expressions identify lymphopoietic tissues in juvenile and adult Chinese giant salamander (Andrias davidianus)

Rag1 and rag2 are two closely linked recombination activating genes required for V(D)J recombination of antigen receptors in immature lymphocytes, whose expression can serve as marker to identify the lymphopoietic tissues. To study the development of lymphopoietic tissues in Chinese giant salamander (Andrias davidianus), the Chinese giant salamander rag1 and rag2 coding sequences were cloned and determined. High transcript levels of rag1 and rag2 were co-detected in the thymus before 14 months of age, whereas levels were lower in spleen, liver and kidney at all stage of development. The spatial expression patterns of rag1 and rag2 were studied in combination with igY and tcrβ gene expression using in situ hybridization. Significant transcript signals for rag1, rag2, tcrβ and igY were detected not only in the thymus and spleen but also the liver and kidney of juvenile and adult Chinese giant salamanders, which suggests that cells of lymphocyte lineage are present in multiple tissues of the Chinese giant salamander. This implies that lymphopoiesis may take place in these tissues. The tissue morphology of thymus suggested that the branched thymic primordium developed into mature organ with the development of thymocyte from juvenile to adult. These results not only confirm that as expected the thymus and spleen are primordial lymphopoietic tissues but also suggest that the liver and kidney provide site of lymphocyte differentiation in Chinese giant salamander.

Characterization and function of the T-box 1 gene in Chinese giant salamander Andrias davidianus

We characterized the Andrias davidianus T-box 1 (Tbx1) gene. Tbx1 expression was high in testis and low in other examined tissues. Immunohistochemistry detected tbx1 expression in somatic and germ cells 62 days post-hatching (dph), prior to gonad differentiation. At 210 dph, after gonad differentiation, tbx1 was expressed in spermatogonia and testis somatic cells and in granulosa cells in ovary. Tbx1 expression was up-regulated in ovary after high temperature treatment. In the neomale, tbx1 expression showed a similar profile to normal males, and vice-versa for genetic male. Over-expression of tbx1 in females after injection of TBX1 protein down-regulated the female-biased genes cyp19a and foxl2 and up-regulated the male-biased amh gene. When tbx1 was knocked down by tbx1/siRNA, cyp19a and foxl2 expression was up-regulated, and expression of amh, cyp26a, dmrt1, and wt1 was down-regulated. Results suggest that tbx1 influenced sex-related gene expression and participates in regulation of A. davidianus testis development.

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).

Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders

Background

Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication.

Results

We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years.

Conclusions

The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.

Identification of the first cathelicidin gene from skin of Chinese giant salamanders Andrias davidianus with its potent antimicrobial activity

Cathelicidins, as effector molecules, play important roles against infections and represent a crucial component of the innate immune system in vertebrates. They are widely studied in mammals, but little is known in amphibians. In the present study, we report the identification and characterization of a novel cathelicidin from Chinese giant salamander Andrias davidianus, which is the first study in Caudata amphibian. The cDNA sequence encodes a predicted 148-amino-acid polypeptide, which composed of a 20-residue signal peptide, a 94-residue conserved cathelin domain and a 34-residue mature peptide. From the multiple sequence alignments and phylogenetic analysis, AdCath shared conserved structure with other orthologs and clustered with other amphibian peptides. The tissue expression profiles revealed AdCath was highly expressed in skin. To study the function of AdCath gene, the AdCath precursor protein and mature peptide were recombinantly expressed and chemical synthesized respectively. The rAdCath protein could bind to LPS in a dose-dependent manner. When the concentrations of rAdCath protein and mature peptide were up to 22 μg/mL, they showed significantly cytotoxicity to human 293T cell lines. The rAdCath protein and synthetic peptide could exhibit antibacterial activities detected by the minimum inhibitory concentrations assay. From the SEM assay, the synthetic mature peptide could destroy the membrane of bacteria and cause loss of membrane integrity. Collectively, these findings characterized the first cathelicidin from A. davidianus, and highlighted its potential antimicrobial activities, indicating its important roles in the skin immune response against different bacteria.

Molecular cloning, characterization and evolutionary analysis of leptin gene in Chinese giant salamander, Andrias davidianus

Leptin is an important hormone possessing diverse physiological roles in mammals and teleosts. However, it has been characterized only in a few amphibian species, and its evolutions are still under debate. Here, the full length of the leptin (Adlep) cDNA of Chinese giant salamander (Andrias davidianus), an early diverging amphibian species, is characterized and according to the results of the primary sequence analysis, tertiary structure reconstruction and phylogenetic analysis is confirmed to be an ortholog of mammalian leptin. An intron was identified between the coding exons of A. davidianus leptin, which indicated that the leptin is present in the salamander genome and contains a conserved gene structure in vertebrates. Adlep is widely distributed but expression levels vary among different tissues, with highest expression levels in the muscle. Additionally, the leptin receptor and other genes were mapped to three known leptin signaling pathways, suggesting that the leptin signaling pathways are present in A. davidianus. Phylogenetic topology of leptins are consistent with the generally accepted evolutionary relationships of vertebrates, and multiple leptin members found in teleosts seem to be obtained through a Cluopeocephala-specific gene duplication event. Our results will lay a foundation for further investigations into the physiological roles of leptin in A. davidianus.

Transcriptome analysis of the Chinese giant salamander (Andrias davidianus) using RNA-sequencing

The Chinese giant salamander (Andrias davidianus) is an economically important animal on academic value. However, the genomic information of this species has been less studied. In our study, the transcripts of A. davidianus were obtained by RNA-seq to conduct a transcriptomic analysis. In total 132,912 unigenes were generated with an average length of 690 bp and N50 of 1263 bp by de novo assembly using Trinity software. Using a sequence similarity search against the nine public databases (CDD, KOG, NR, NT, PFAM, Swiss-prot, TrEMBL, GO and KEGG databases), a total of 24,049, 18,406, 36,711, 15,858, 20,500, 27,515, 36,705, 28,879 and 10,958 unigenes were annotated in databases, respectively. Of these, 6323 unigenes were annotated in all database and 39,672 unigenes were annotated in at least one database. Blasted with KEGG pathway, 10,958 unigenes were annotated, and it was divided into 343 categories according to different pathways. In addition, we also identified 29,790 SSRs. This study provided a valuable resource for understanding transcriptomic information of A. davidianus and laid a foundation for further research on functional gene cloning, genomics, genetic diversity analysis and molecular marker exploitation in A. davidianus.

Analysis on the expression and function of a chicken-type and goose-type lysozymes in Chinese giant salamanders Andrias davidianus

Lysozymes as an important immune factor, play vital roles in innate immune response against pathogen infection. In the present study, one c-type and g-type lysozymes were identified from Chinese giant salamander (Andrias davidianus). They shared highly conserved structural features with lysozymes from other species. Spatial expression analysis revealed that AdlysC transcript was most abundant in liver and stomach, and least in muscle and brain. In contrast, the expression level of AdlysG was most abundant in liver and least in muscle and skin. The transcription level of c-type and g-type lysozymes were up-regulated after Aeromonas hydrophila infection in liver and spleen, indicating their participations in the immune response. Moreover, the recombinant AdlysC and AdlysG protein were produced and purified, and were used to investigate the lysozyme activity at different pH and temperatures. The optimal lytic activity was determined at pH 6.0 and at a temperature of 30 °C. Through the minimal inhibitory concentration test, the rAdlysC and rAdlysG exhibited apparent antibacterial activity against both Gram-positive and Gram-negative bacteria with a variable concentration. In conclusion, it is the first report of lysozymes in A. davidianus, and c-type and g-type lysozymes should be involved in the innate immune response of A. davidianus.

The complete mitochondrial genome sequence of the Sichuan Digging Frog, Kaloula rugifera (Anura: Microhylidae) and its phylogenetic implications

The Sichuan Digging Frog (Kaloula rugifera) belongs to the family Dicroglossidae, which is endemic to northeastern Sichuan and southernmost Gansu provinces, in southwestern China. In this study, the complete mitochondrial genome of K. rugifera was sequenced. The mitogenome was 17,074bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding control region. As in other vertebrates, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes which are encoded on the light strand. The overall base composition of the K. rugifera is 30.32% A, 25.76% C, 29.72% T, and 14.20% G, which is consistent with the lowest frequency for G content in typical amphibian animals' mitochondrial genomes. The alignment of the Kaloula species control regions exhibited high genetic variability and rich A+T content. Besides, 3 types of tandem repeat units were also identified in the control region. Phylogenetic tree demonstrated that K. rugifera was clustered together with K. borealis and K. verrucosa and they had a close relationship with each other. The complete mitogenome of K. rugifera can provide an important data for the studies on phylogenetic relationship to further explore the taxonomic status of Kaloula species.

Identification and characterization of known and novel microRNAs in three tissues of Chinese giant salamander base on deep sequencing approach

MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in the regulation of gene expression. Although animal miRNAs have been extensively studied in model systems, less is known in other animal with limited genome sequence data, including Chinese giant salamander (Andrias davidianus). The identification of the full repertoire of miRNAs expressed in the liver, spleen and muscle of A. davidianus would significantly increase our understanding for physiological function of A. davidianus, in this ancient and endangered urodele amphibian. In this study, three independent small RNA libraries were constructed from the liver, spleen and muscle of A. davidianus. The libraries were subjected to high-throughput sequencing by using the Illumina deep sequencing. As a result, a total of 12,831,239, 13,592,195 and 9,887,531 raw reads representing 2,240,771, 1,363,266 and 1,964,252 clean reads per library were obtained separately. Through bioinformatics analysis, we identified total of 553 known miRNAs and 44 putative novel miRNAs in our small RNA dataset from liver, spleen and muscle tissues. Five known miRNAs (gga-miR-10a-5p, pma-miR-29d-5p, aca-miR-338-3p, hsa-miR-455-3p and ssa-miR-2184-5p_R-1) and three novel miRNAs (PC-5p-891_1763, PC-5p-32538_50 and PC-3p-33645_48) showed different expression in eight different tissues as revealed by stem-loop qPCR analysis. This study characterized the miRNA of A. davidianus for the first time, which provides an opportunity for further understanding of miRNA regulation function in A. davidianus ranked as living fossils.

Characterization of galectin-1 from Chinese giant salamanders Andrias davidianus and its involvements during immune response

Galectins are considered as a multifunctional protein which play essential roles in cell adhesion and apoptosis, inflammation, tumor progression and immune response. In spite of extensive studies of galectin importance in immune system among different animals, few studies have been devoted to their functions in amphibian. In the present study, we characterized one proto type of galectin (named AdGal1) from Chinese giant salamander Andrias davidianus and studied its function in immune response. AdGal1 cDNA possesses an open reading frame of 598 bp, which encodes a putative galectin of 134 amino acids containing one carbohydrate recognition domains (CRDs). The constitutive expression of mRNA transcripts was detected in a wide range of tissues, with the highest expression in kidney. Immune challenges with Aeromonas hydrophila and Chinese giant salamander iridovirus (GSIV), the transcript level of AdGal1 in kidney was significantly upregulated. The mature protein of AdGal1 was successfully expressed and purified in Escherichia coli BL21 (DE3). The recombinant AdGal1 (rAdGal1) could show bind activity to different Gram negative and Gram positive bacteria. It could also strongly agglutinate different kinds of bacteria at different concentrations. Collectively, these data from the present study indicate that AdGal1 is a vital pattern recognition receptor to recognize different microbes in the innate immune system of Andrias davidianus.

Sequencing and de novo transcriptome assembly of the Chinese giant salamander (Andrias davidianus)

Next-generation technologies for determination of genomics and transcriptomics composition have a wide range of applications. Andrias davidianus, has become an endangered amphibian species of salamander endemic in China. However, there is a lack of the molecular information. In this study, we obtained the RNA-Seq data from a pool of A. davidianus tissue including spleen, liver, muscle, kidney, skin, testis, gut and heart using Illumina HiSeq 2500 platform. A total of 15,398,997,600 bp were obtained, corresponding to 102,659,984 raw reads. A total of 102,659,984 reads were filtered after removing low-quality reads and trimming the adapter sequences. The Trinity program was used to de novo assemble 132,912 unigenes with an average length of 690 bp and N50 of 1263 bp. Unigenes were annotated through number of databases. These transcriptomic data of A. davidianus should open the door to molecular evolution studies based on the entire transcriptome or targeted genes of interest to sequence. The raw data in this study can be available in NCBI SRA database with accession number of SRP099564.

The complete mitochondrial genome of the Neophocaena asiaeorientails sunamrei (Phocaenidae: Neophocaena)

The complete sequence of the mitochondrial DNA (mtDNA) of the East Asian finless porpoise (Neophocaena asiaeorientails sunamrei) has been completed and annotated in this study. The circular genome is 16 385 bp in length, consisting of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and 2 non-coding areas. Ten protein-coding genes use ATG as the initiation codon, except that ND2 started with ATT, and ND3 and ND5 started with ATA. Canonical TAA and TAG termination codons are found in seven protein-coding genes, the remaining four (ND4, ND2, ND3, and COIII) have incomplete termination codons (T or TA), in addition, the COI ends with AGG, and the Cyt b ends with AGA. On the neighbor-joining phylogenetic tree, these results show that the relationship of Phocoenidae is closer with Delphinidae.

The complete mitochondrial genome of White-cheeked macaque (Macaca leucogenys)

White-cheeked macaque is a newly described species in genus Macaca. Here, the complete mitochondrial genome was firstly determined, which was deposited in Genbank with accession number KU564271. The length of the mitogenome is 16,494 base pairs (bp). It includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a control region. We used 12 mitochondrion genes of 15 species to constructed phylogenetic tree with three methods. The study will provide useful data for further studies on phylogenetic relationships as well as population structure, biodiversity conservation and conservation genetics.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

The complete mitochondrial genome of bearded pig, Sus barbatus, and comparative mitochondrial genomics of Cetartiodactyla

In this study, the complete mitochondrial genome sequence of bearded pig, Sus barbatus, with the total length of 16,480 bp, is determined for the first time. This mitogenome harbors 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region (D-loop). The overall base composition is A (34.80%), C (26.07%), G (13.12%), and T (26.01%), so the slight A-T bias (60.81%) was detected. Most of the genes are distributed on the H-strand, except for the ND6 subunit gene and eight tRNA genes. To obtain the phylogenetic relationship of the Cetartiodactyla, 11 mitochondrial genomes were used for phylogenetic analysis. The mitochondrial genome of S. barbatus presented here will contribute to a better understanding of the population genetics.

The complete mitochondrial genome sequence of Tylototriton taliangensis (Amphibia: Caudata)

Tylototriton taliangensis was listed as a Near Threatened amphibian in IUCN red list. In this study, we sequenced the complete mitochondrial (mt) genome of this species (GenBank: KP979646) and found it contains 16,265 base pairs, which encode 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA) and 1 control region (CR). We also found that almost all PCGs and tRNA genes are located on the H-strand, except for ND6 subunit gene and eight tRNA genes, which were distributed on the L-strand. The PCGs used "ATG" and "GTG" as the start codon, while used four types of stop codons. Almost all tRNA genes were folded into typical cloverleaf secondary structures. The L-strand replication origin (OL) and a non-coding region were also found. The new mitogenomic phylogenetic tree confirms the reciprocally monophyly of the genus Tylototriton, Echinotriton and Pleurodeles with high bootstrap value. The present study will provide information for future studies on the conservation genetics and phylogeny of this species and its relatives.

Cloning, sequence analysis and expression profiles of Toll-like receptor 7 from Chinese giant salamander Andrias davidianus

The Chinese giant salamander, Andrias davidianus, is the largest extant amphibian species in the world, which is of significance due to its specific position in the evolutionary history of vertebrates. Currently, limited information about the innate immune system of this animal is known. In this study, the toll-like receptor 7 (TLR7), designated CgsTLR7, was cloned from Chinese giant salamander, A. davidianus. The full-length cDNA of CgsTLR7 is 3747 bp, with an open reading frame of 3150 bp, encoding 1049 amino acids. The TLR family motifs, including the leucine-rich repeat (LRR) and Toll/interleukin (IL)-1 receptor (TIR) domain are conserved in CgsTLR7, which includes 19 LRRs and a TIR domain. The predicted amino acid sequence of CgsTLR7 has 71%, 65%, 63% and 55% identity with turtle, chicken, human and fugu TLR7 homologues, respectively. Phylogenetic analysis showed that CgsTLR7 is closest to that of frog TLR7 among the examined species. Quantitative real-time PCR analysis revealed broad expression of CgsTLR7 in tissues from apparently healthy Chinese giant salamanders with the highest expression in the liver and the lowest expression in the intestine. The mRNA expression was up-regulated and reached a peak level in the kidney, liver and spleen at 12 h, 24 h and 48 h after infecting the animals with the giant salamander iridovirus (GSIV), respectively. These results suggest that CgsTLR7 has a conserved gene structure and might play an important role in immune regulation against viral infections in the Chinese giant salamander.

The complete mitochondrial genome of Stump-tailed Macaques (Macaca arctoides)

Stump-tailed macaque (Macaca arctoides) has been enlisted as the Near Threatened species in the IUCN Red List. In this study, the complete mitochondial genome of M. arctoides was determined. The mitogenome was 16,559 bp in length with an A + T content of 56.8%, containing 2 ribosomal RNA genes (rRNA), 22 transfer RNA genes (tRNA), 13 protein-coding genes (PCGs) and a control region (D-loop). Ten protein-coding genes (ND1, COX1, COX2, ND4, COX3, ND4L, ATP8, ATP6, CYTB, ND6) started with ATG codons, while ND2, ND3 and ND5 initiated with ATT, ATC and ATA, respectively. Eight PCGs ended with complete termination codons except for COX3, ND3, ND4, and CYTB terminated with incompleted codon T.

The complete mitochondrial genome of Assamese Macaques (Macaca assamensis)

The complete mitochondrial genome of Assamese macaques (Macaca assamensis) was sequenced in this study. The genome is 16,542 bp long, consisting of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and 2 non-coding areas. Eight PCGs (COI, COII, ATP6, COIII, ND4L, ND4, ND6, CYTB) initiate with the start codon ATG and another two genes (ND1, ATP8) use GTG, while ND2, ND3 and ND5 start with ATT, ATC and ATA, respectively. Five genes (COII, ATP8, ATP6, ND4L and ND5) use the complete stop codon TAA, whereas four genes have incomplete stop codons, TA- (COIII) and T- - (ND3, ND4, CYTB), while others use standard canonical TAA as their termination codons. The largest non-coding control region with the length of 1091 bp is located between the tRNA-Pro and tRNA-Phe genes.

Extensive diversification of MHC in Chinese giant salamanders Andrias davidianus (Anda-MHC) reveals novel splice variants

A series of MHC alleles (including 26 class IA, 27 class IIA, and 17 class IIB) were identified from Chinese giant salamander Andrias davidianus (Anda-MHC). These genes are similar to classical MHC molecules in terms of characteristic domains, functional residues, deduced tertiary structures and genetic diversity. The majority of variation between alleles is found in the putative peptide-binding region (PBR), which is driven by positive Darwinian selection. The coexistence of two isoforms in MHC IA, IIA, and IIB alleles are shown: one full-length transcript and one novel splice variant. Despite lake of the external domains, these variants exhibit similar subcellular localization with the full-length transcripts. Moreover, the expression of MHC isoforms are up-regulated upon in vivo and in vitro stimulation with Andrias davidianus ranavirus (ADRV), suggesting their potential roles in the immune response. The results provide insights into understanding MHC variation and function in this ancient and endangered urodele amphibian.

The complete mitochondrial genome sequence analysis of Tibetan argali (Ovis ammon hodgsoni): implications of Tibetan argali and Gansu argali as the same subspecies

The genus Ovis (Bovidae, Artiodactyla) includes six species, i.e. Ovis ammon, Ovis aries, Ovis canadensis, Ovis dalli, Ovis nivicola and Ovis vignei. Based on morphology, geographical location, habitat, etc., the species O. ammon is divided into nine subspecies. The near threatened Tibetan argali is distributed across the Tibetan Plateau and its peripheral mountains, and believed to be one of the O. ammon subspecies (O. a. hodgsoni). However, considering its morphological features and distributions, a question has been proposed by some researchers about the subspecies status of Tibetan argali. In this study, we employed complete mitochondrial DNA (mtDNA) to explore the phylogenetic relationship and population genetic structure of Tibetan argali. The results revealed that the nucleotide composition, gene arrangement and codon usage pattern of the mitochondrial genome of Tibetan argali are similar to those of other caprines. Phylogenetic analyses showed that Tibetan argali was clustered with O. ammon. Interestingly, five Tibetan argali individuals and one of the three Gansu argali (O. a. dalailamae) individuals were clustered in the same branch, which is a sister group to other two Gansu argali individuals. Together with morphological characteristics, our results suggested that Tibetan argali and Gansu argali may belong to the same subspecies (O. a. hodgsoni) of O. ammon, rather than two different subspecies.

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.

Complete mitochondrial genome of a new vole Proedromys liangshanensis (Rodentia: Cricetidae) and phylogenetic analysis with related species: are there implications for the validity of the genus Proedromys?

AIM. The complete mitochondrial genome sequence of a newly discovered vole, Proedromys liangshanensis (Rodentia: Cricetidae: Arvicolinae), was determined. RESULTS. The mitogenome of P. liangshanensis is 16,296 bp in length. As with most other mammals, it contains the same gene order and an identical number of genes or regions, including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and one putative control region. The A+T content of the control region is 58.2%, the lowest proportion detected so far in Myomorpha. To confirm the phylogenetic position of P. liangshanensis, we carried out phylogenetic analyses based on complete mitochondrial genomic data using Bayesian, maximum parsimony, and maximum likelihood methods. CONCLUSION. All results revealed that P. liangshanenis is sister to Microtus. Although the results do not bear light on the validity of the genus Proedromys, based on the morphological characters, we suggest that Proedromys is an independent genus of equal rank to the genus Microtus.

The complete mitochondrial genome of dhole Cuon alpinus: phylogenetic analysis and dating evolutionary divergence within Canidae

The dhole (Cuon alpinus) is the only existent species in the genus Cuon (Carnivora: Canidae). In the present study, the complete mitochondrial genome of the dhole was sequenced. The total length is 16672 base pairs which is the shortest in Canidae. Sequence analysis revealed that most mitochondrial genomic functional regions were highly consistent among canid animals except the CSB domain of the control region. The difference in length among the Canidae mitochondrial genome sequences is mainly due to the number of short segments of tandem repeated in the CSB domain. Phylogenetic analysis was progressed based on the concatenated data set of 14 mitochondrial genes of 8 canid animals by using maximum parsimony (MP), maximum likelihood (ML) and Bayesian (BI) inference methods. The genera Vulpes and Nyctereutes formed a sister group and split first within Canidae, followed by that in the Cuon. The divergence in the genus Canis was the latest. The divarication of domestic dogs after that of the Canis lupus laniger is completely supported by all the three topologies. Pairwise sequence divergence data of different mitochondrial genes among canid animals were also determined. Except for the synonymous substitutions in protein-coding genes, the control region exhibits the highest sequence divergences. The synonymous rates are approximately two to six times higher than those of the non-synonymous sites except for a slightly higher rate in the non-synonymous substitution between Cuon alpinus and Vulpes vulpes. 16S rRNA genes have a slightly faster sequence divergence than 12S rRNA and tRNA genes. Based on nucleotide substitutions of tRNA genes and rRNA genes, the times since divergence between dhole and other canid animals, and between domestic dogs and three subspecies of wolves were evaluated. The result indicates that Vulpes and Nyctereutes have a close phylogenetic relationship and the divergence of Nyctereutes is a little earlier. The Tibetan wolf may be an archaic pedigree within wolf subspecies. The genetic distance between wolves and domestic dogs is less than that among different subspecies of wolves. The domestication of dogs was about 1.56-1.92 million years ago or even earlier.

The complete mitochondrial genome of the large-headed frog, Limnonectes bannaensis (Amphibia: Anura), and a novel gene organization in the vertebrate mtDNA

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of the large-headed frog, Limnonectes bannaensis (Amphibia, Anura) by using polymerase chain reaction (PCR). The entire mtDNA sequence is 16,867 bp in length with a novel case of tRNAs in vertebrates. This mt genome is characterized by three distinctive features: (1) a tandem duplication of tRNA(Met) gene is observed, (2) the tRNA(Ala), tRNA(Asn), tRNA(Cys) and tRNA(Glu) genes coded on the L-strand are absent from the L. bannaensis mtDNA, the tRNA(Cys) and tRNA(Glu) genes change into tRNA pseudogenes by reason of degenerative anticodon, and a noncoding sequence of 206 nt long (NC1) has replaced the original position of other two tRNAs, (3) besides NC1, another three noncoding spacers (NC2-4) longer than 50 bp are found in the broken WANCY region and the region NC3-ND5-NC4-ND6-PsiE-Cytb-CR of the new sequence. These features could be explained by a model of gene duplication and deletion. The new sequence data was used to assess the phylogenetic relationships among 25 species of Anura using neighbor-joining, Bayesian, and maximum likelihood methods, and the phylogenetic tree shows the rice frog Fejervarya limnocharis is closest to L. bannaensis in the study.