Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

A time-course transcriptome analysis revealing the potential molecular mechanism of early gonadal differentiation in the Chinese giant salamander

The Chinese giant salamander (CGS) Andrias davidianus is the largest extant amphibian and has recently become an important species for aquaculture with high economic value. Meanwhile, its wild populations and diversity are in urgent need of protection. Exploring the mechanism of its early gonadal differentiation will contribute to the development of CGS aquaculture and the recovery of its wild population. In this study, transcriptomic and phenotypic research was conducted on the critical time points of early gonadal differentiation of CGS. The results indicate that around 210 days post-hatching (dph) is the critical window for female CGS's gonadal differentiation, while 270 dph is that of male CGS. Besides, the TRPM1 gene may be the crucial gene among many candidates determining the sex of CGS. More importantly, in our study, key genes involved in CGS's gonadal differentiation and development are identified and their potential pathways and regulatory models at early stage are outlined. This is an initial exploration of the molecular mechanisms of CGS's early gonadal differentiation at multiple time points, providing essential theoretical foundations for its captive breeding and offering unique insights into the conservation of genetic diversity in wild populations from the perspective of sex development.

Identification and Expression Analysis of Wnt2 Gene in the Sex Differentiation of the Chinese Soft-Shelled Turtle (Pelodiscus sinensis)

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important freshwater aquaculture animal in China. The Wnt gene family plays important regulatory roles in the development and growth of mammals. However, the precise function of these family genes has not been well understood in the sex differentiation of Chinese soft-shelled turtles. Here, we cloned a member of the Wnt family, Wnt2, which obtained a 1077 bp open reading frame that encoded a 358-aa protein. The putative amino acid sequences of proteins are exceeded 80% identical to other turtles. The expression level of Wnt2 peaked at the 14th stage both in female and male embryos during the early gonadal differentiation period of Chinese soft-shelled turtles, which occurred before gonadal differentiation. Wnt2 mRNA was expressed at higher levels in the brains and gonads of mature P. sinensis females compared with those in mature males. Wnt agonists significantly affected the expression level of Wnt2 during the gonadal differentiation period. After Wnt agonists (1.0 μg/μL, 2.5 μg/μL, 5.0 μg/μL) treatment, the expression level of the Wnt2 generally appeared to have an inverted-V trend over time in female embryonic gonads. The results suggested that Wnt2 may participate in the regulation of gonad development in P. sinensis during the early embryonic stages. These results could provide a theoretical basis for the reproduction process of the Chinese soft-shelled turtle.

Potential antagonistic relationship of fgf9 and rspo1 genes in WNT4 pathway to regulate the sex differentiation in Chinese giant salamander (Andrias davidianus)

Farmed chinese giant salamander (Andrias davidianus) was an important distinctive economically amphibian that exhibited male-biased sexual size dimorphism. Fgf9 and rspo1 genes antagonize each other in Wnt4 signal pathway to regulate mammalian gonadal differentiation has been demonstrated. However, their expression profile and function in A. davidianus are unclear. In this study, we firstly characterized fgf9 and rspo1 genes expression in developing gonad. Results showed that fgf9 expression level was higher in testes than in ovaries and increased from 1 to 6 years while rspo1 expression was higher in ovaries than in testes. In situ hybridization assay showed that both fgf9 and rspo1 genes expressed at 62 dpf in undifferentiated gonad, and fgf9 gene was mainly expressed in spermatogonia and sertoli cells in testis while strong positive signal of rspo1 was detected in granular cell in ovary. During sex-reversal, fgf9 expression was significantly higher in reversed testes and normal testes than in ovaries, and opposite expression pattern was detected for rspo1. When FH535 was used to inhibit Wnt/β-catenin pathway, expression of rspo1, wnt4 and β-catenin was down-regulated. Conversely, expression of fgf9, dmrt1, ftz-f1 and cyp17 were up-regulated. Furthermore, when rspo1 and fgf9 were knocked down using RNAi technology, respectively. We observed that female biased genes were down regulated in ovary primordial cells after rspo1 was knocked down, while the opposite expression profile was observed in testis primordial cells after fgf9 was knocked down. These results suggested that fgf9 and rspo1 played an antagonistic role to regulate sex differentiation in the process of the gonadal development and provided a foundation for further functional characterizations. The data also provided basic information for genome editing breeding to improve the Chinese giant salamander farming industry.

Comprehensive Transcriptome Analysis of Gonadal and Somatic Tissues for Identification of Sex-Related Genes in the Largemouth Bass Micropterus salmoides

Largemouth bass (Micropterus salmoides) is an economically important fish. It can spawn many times during a breeding season, and there are no obvious morphological characteristics to distinguish male and female juvenile fish. So far, little is known about the genes regulating their sexual development in this species. Here, we performed RNA sequencing (RNA-Seq) analysis of the testis, ovary, and somatic tissue to identify sex-related genes in the largemouth bass. A total of 51,672 unigenes were obtained via the transcriptome analysis, and 5900 differential expression genes (DEGs), including 3028 up-regulated and 2872 down-regulated DEGs, were obtained in the somatic tissue, testis, and ovary. DEGs were retrieved by making comparisons: somatic tissue vs testis (1733-up and 1382-down), testis vs ovary (841-up and 807-down), and ovary vs somatic tissue (454-up and 683-down). Finally, functional annotation identified 22 key sex-related DEGs, including 13 testis-biased DEGs (dmrt1, cyp11b1, sox9, spata4, spata22, spata17, fshr, fem-1a, wt1, daz1, amh, vasa, and piwi1) and 9 ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, bmp15, fem-1b, fig. la, and piwi2). This result was further confirmed by the tissue expression detection via RT-PCR and RT-qPCR. Protein-protein interacting (PPI) network analysis revealed that the testis-specific dmrt1 interacts directly with the testis-biased DEGs (cyp11b1 and spata4) and the ovary-biased DEGs (foxl2, gdf9, zp3, sox3, cyp19a, and bmp15), suggesting that the dmrt1 as a sex-determining gene can play a dual role through inducing the testis-biased DEGs and inhibiting the ovary-biased DEGs during the testicular development. Our present results provide useful molecular data for a better understanding of sexual development in the largemouth bass.

Molecular Identification of Nocardia seriolae and Comparative Analysis of Spleen Transcriptomes of Hybrid Snakehead (Channa maculata Female × Channa argus Male) With Nocardiosis Disease

Hybrid snakehead (Channa maculata female × Channa argus male) is a new freshwater aquaculture fish species in southern China. During intensive aquaculture, hybrid snakeheads are often infected by Nocardia seriolae. In this study, hybrid snakehead infected suspiciously by N. seriolae in an artificial breeding pond were examined. Diseased hybrid snakeheads swam slowly without food intake, and the clinical symptoms included skin wound, anal swelling and ascites, and white granulomatous in liver, spleen, and kidney of fish. Through bacterial isolation, 16S rDNA sequencing, fluorescence in situ hybridization (FISH) and artificial infection experiment, the pathogen was identified as N. seriolae. Furthermore, the spleen samples from diseased and healthy male hybrid snakeheads in the same pond were used for RNA-Seq analysis. A total of 3,512 unique transcripts (unigenes) were identified as differentially expressed genes (DEGs), and 1,886 of them were up-regulated in diseased fish. The expression patterns of 20 DEGs were verified by quantitative polymerase chain reaction (qPCR). Several immune-related pathways and many immune-related genes were identified. qPCR results showed that the expression patterns of immune-related genes in the liver and kidney of diseased fish were comparable to that in the spleen. This study provides deep-sequencing data of hybrid snakehead spleen and will help understand the immune response of hybrid snakehead to N. seriolae. It is also helpful for the biomarker screening of fish-borne Nocardia spp. and the breeding of nocardiosis-resistant fish species.

Gene Expression of Takifugu rubripes Gonads During AI- or MT-induced Masculinization and E2-induced Feminization

Elucidating the global molecular changes that occur during aromatase inhibitor (AI)- or 17α-methyltestosterone (MT)-induced masculinization and estradiol-17β (E2)-induced feminization is critical to understanding the roles that endocrine and genetic factors play in regulating the process of sex differentiation in fish. Here, fugu larvae were treated with AI (letrozole), MT, or E2 from 25 to 80 days after hatching (dah), and gonadal transcriptomic analysis at 80 dah was performed. The expression of dmrt1, gsdf, foxl2, and other key genes (star, hsd3b1, cyp11c1, cyp19a1a, etc.) involved in the steroid hormone biosynthesis pathway were found be altered. The expression of dmrt1, gsdf, cyp19a1a, and foxl2 was further verified by quantitative polymerase chain reaction. In the control group, the expression of dmrt1 and gsdf was significantly higher in XY larvae than in XX larvae, while the expression of foxl2 and cyp19a1a was significantly higher in XX larvae than in XY larvae (P < .05). AI treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 and gsdf in XX larvae. MT treatment suppressed the expression of foxl2, cyp19a1a, dmrt1, and gsdf in XX larvae. E2 treatment suppressed the expression of dmrt1 and gsdf, but did not restore the expression of foxl2 and cyp19a1a in XY larvae. The shared response following AI, MT, and E2 treatment suggested that these genes are essential for sex differentiation. This finding offers some insight into AI or MT-induced masculinization, and E2-induced femininization in fugu.

Semen and oocyte collection, sperm cryopreservation and IVF with the threatened North American giant salamander Cryptobranchus alleganiensis

Semen of high to moderate quality was collected following the hormonal induction of North American giant salamanders Cryptobranchus alleganiensis. Oocytes from one female yielded the first C. alleganiensis produced while maintained in aquaria under human care and the first externally fertilising salamander produced with cryopreserved spermatozoa and IVF. Further research is needed with North American giant salamanders to establish reliable techniques to produce large numbers of viable offspring, along with the application of cryopreserved spermatozoa.

Full-length transcriptome assembly of Andrias davidianus (Amphibia: Caudata) skin via hybrid sequencing

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2039 bp and N50 of 2172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.

A high-density genetic map construction and sex-related loci identification in Chinese Giant salamander

BACKGROUND

The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci.

RESULTS

In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively.

CONCLUSION

We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

Environmental Temperatures Affect the Gastrointestinal Microbes of the Chinese Giant Salamander

An increasing number of studies have shown that warming also influences the animal gut microbiome (altering the community structure and decreasing its diversity), which might further impact host fitness. Here, based on an analysis of the stomach and gut (the entire intestine: from the anterior intestine to the cloaca) microbiome in laboratory larva of giant salamanders (Andrias davidianus) under different living water temperatures (5, 15, and 25°C) at two sample time points (80 and 330 days after the acclimation), we investigated the potential effect of temperature on the gastrointestinal microbiome community. We found the significant Interaction between sampling time and temperature, or type (stomach and gut) on Shannon index in the gastrointestinal microbiome of the giant salamanders. We also found the significant difference in Shannon index among temperature groups within the same sample type (stomach or gut) at each sample time. 10% of variation in microbiome community could be explained by temperature alone in the total samples. Both the stomach and gut microbiomes displayed the highest similarity in the microbiome community (significantly lowest pairwise unweighted Unifrac distance) in the 25-degree group between the two sampling times compared to those in the 5-degree and 15-degree groups. Moreover, the salamanders in the 25°C treatment showed the highest food intake and body mess compared to that of other temperature treatments. A significant increase in the abundance of Firmicutes in the gastrointestinal microbiome on day 330 with increasing temperatures might be caused by increased host metabolism and food consumption. Therefore, we speculate that the high environmental temperature might indirectly affect both alpha and beta diversity of the gastrointestinal microbiome.

Genome-Wide DNA Methylation and RNA Analysis Reveal Potential Mechanism of Resistance to Streptococcus agalactiae in GIFT Strain of Nile Tilapia (Oreochromis niloticus )

Streptococcus agalactiae is an important pathogenic bacterium causing great economic loss in Nile tilapia (Oreochromis niloticus) culture. Resistant and susceptible groups sharing the same genome showed significantly different resistance to S. agalactiae in the genetically improved farmed tilapia strain of Nile tilapia. The resistance mechanism is unclear. We determined genome-wide DNA methylation profiles in spleen of resistant and susceptible O. niloticus at 5 h postinfection with S. agalactiae using whole-genome bisulfite sequencing. The methylation status was higher in the spleen samples from resistant fish than in the susceptible group. A total of 10,177 differentially methylated regions were identified in the two groups, including 3725 differentially methylated genes (DMGs) (3129 hyper-DMGs and 596 hypo-DMGs). The RNA sequencing showed 2374 differentially expressed genes (DEGs), including 1483 upregulated and 891 downregulated. Integrated analysis showed 337 overlapping DEGs and DMGs and 82 overlapping DEGs and differentially methylated region promoters. By integrating promoter DNA methylation with gene expression, we revealed four immune-related genes (Arnt2, Nhr38, Pcdh10, and Ccdc158) as key factors in epigenetic mechanisms contributing to pathogen resistance. Our study provided systematic methylome maps to explore the epigenetic mechanism and reveal the methylation loci of pathogen resistance and identified methylation-regulated genes that are potentially involved in defense against pathogens.

RNA-seq analysis provides insight into molecular adaptations of Andrias davidianus

The Chinese giant salamander Andrias davidianus is regarded as an ideal model for studying local adaptations, such as longevity, tolerance to starvation, and cutaneous respiration. Transcriptome analysis is useful for studying the large and complex genomes of amphibians. Based on the coding gene set of adult A. davidianus, dozens of A. davidianus-specific genes were identified and three signaling pathway (JAK-STAT, HIF-1, and FoxO) genes were expanded as compared with other amphibians. The results of the pathway analysis of A. davidianus-specific genes indicated that the molecular adaptation of A. davidianus may have required a more rapid evolution of the immune system. Additionally, for the first time, the gene expressions in different parts of the skin tissue were compared. The results of the comparison analysis demonstrated that lateral skin could be more focused on mucus secretion, dorsal skin on immunity and melanogenesis, and abdominal skin on water and salt metabolism. This study provides the first insight into studying longevity and starvation tolerance in A. davidianus, and offers a basis for further investigation of the molecular mechanisms of adaptations in amphibians.

Genome-wide RAD sequencing to identify a sex-specific marker in Chinese giant salamander Andrias davidianus

Background

Chinese giant salamander Andrias davidianus is an endangered species. The success of artificial breeding provides a useful way to protect this species. However, the method to identify the sex and mechanism of sex determination were unclear which hinder the improvement of the artificial breeding. Detection of a sex specific marker provides an effective approach to identify genetic sex and investigate the sex determination mechanism.

Results

We used restriction-site-associated DNA (RAD) sequencing to isolate a sex-specific genetic marker in A. davidianus to expand knowledge of the sex determination mechanism. Four male and four female specimens were subjected to RAD sequencing, which generated 934,072,989 reads containing approximately 134.4 Gb of sequences. The first round of comparison of the assembled sequence against the opposite sex raw reads revealed 19,097 female and 17,994 male unmatched sequences. Subsequently, 19,097 female sequences were subjected to a BLAST search against male genomic data, which revealed 308 sequences unmapped to the male genome. One hundred of these were randomly selected and validated by PCR in five male and five female specimens, and four putative sex-specific sequences were produced. Further validation was performed by PCR in another 24 females and 24 males, and all female individuals exhibited the expected specific bands, while the males did not. To apply the sex-specific marker, three specimens reversed from genetic female to physiological male were found in a group exposed to elevated temperature, and 13 individuals reversed from genetic male to physiological female were obtained in a 17β-estradiol exposed group.

Conclusion

This is the first report of a sex-specific marker in A. davidianus and may have potential for elucidation of its sex determination mechanism and, hence, its conservation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5771-5) contains supplementary material, which is available to authorized users.

The First Transcriptome Assembly of Yenyuan Stream Salamander (Batrachuperus yenyuanensis) Provides Novel Insights into Its Molecular Evolution

The Yenyuan stream salamander (Batrachuperus yenyuanensis) has been previously evaluated with regards to phylogeny, population genetics, and hematology, but genomic information is sparse due to the giant genome size of salamanders which contain highly repetitive sequences, thus resulting in the lack of a complete reference genome. This study evaluates the encoding genetic sequences and provides the first transcriptome assembly of Yenyuan stream salamander based on mixed samples from the liver, spermary, muscle and spleen tissues. Using this transcriptome assembly and available encoding sequences from other vertebrates, the gene families, phylogenetic status, and species divergence time were compared or estimated. A total of 13,750 encoding sequences were successfully obtained from the transcriptome assembly of Yenyuan stream salamander, estimated to contain 40.1% of the unigenes represented in tetrapod databases. A total of 88.79% of these genes could be annotated to a biological function by current databases. Through gene family clustering, we found multiple possible isoforms of the Scribble gene-whose function is related to regeneration-based on sequence similarity. Meanwhile, we constructed a robust phylogenetic tree based on 56 single-copy orthologues, which indicates that based on phylogenetic position, the Yenyuan stream salamander presents the closest relationship with the Chinese giant salamander (Andrias davidianus) of the investigated vertebrates. Based on the fossil-calibrated phylogeny, we estimated that the lineage divergence between the ancestral Yenyuan stream salamander and the Chinese giant salamander may have occurred during the Cretaceous period (~78.4 million years ago). In conclusion, this study not only provides a candidate gene that is valuable for exploring the remarkable capacity of regeneration in the future, but also gives an interesting insight into the understanding of Yenyuan stream salamander by this first transcriptome assembly.