Transcriptome-based genome-wide analysis reveals hybridization dynamics and genetic structure of Japanese giant salamanders

The Japanese giant salamander ( Andrias japonicus ), an apex predator and a Special Natural Monument in Japan, is threatened by hybridization with introduced Chinese giant salamanders ( Andrias davidianus ). This hybridization has caused genetic introgression and expansion of hybrid populations, posing a serious conservation risk. Because morphological identification of hybrids is occasionally unreliable and current genetic methods rely on limited markers, a genome-wide approach is required. However, the extremely large genome (∼50 Gb) of giant salamanders has hindered whole-genome analyses.

In this study, we conducted transcriptome-based analyses of Japanese giant salamanders, Chinese giant salamanders, and their hybrids, generating RNA-seq data from 34 individuals. A total of over 419,000 SNP candidates were identified, from which 4,457 high-confidence SNPs in highly expressed genes were selected for analysis. Population structure analyses for Nabari colony revealed that hybrid individuals form two major groups, corresponding to different degrees of genetic contribution from Japanese and Chinese lineages. Most hybrids were inferred to be F2 or backcross individuals, while F1 hybrids were rare. Mitochondrial analysis indicated that all hybrids possessed Japanese-type mitochondrial genome, suggesting male-mediated introgression from Chinese salamanders.

Differential expression analysis revealed enhanced stress-response pathways in hybrids and stronger antiviral responses in Japanese individuals. Using the axolotl genome as a reference, we constructed a virtual chromosomal map, identifying large haplotype blocks and supporting recent hybridization with limited recombination. This study provides a genome-wide framework for understanding hybridization dynamics and supports future conservation and evolutionary studies.