The Japanese eel (Anguilla japonica), European eel (Anguilla anguilla), and American eel (Anguilla rostrata) are migratory, catadromous, temperate zone fish sharing several common life cycle features. The population genetics of panmixia in these eel species has already been investigated. Our extensive population genetics analysis was based on 1400 Gb of whole-genome sequence (WGS) data from 84 eels. It demonstrated that a Japanese eel group from the Kuma River differed from other populations of the same species. Even after removing the potential adapted/selected single nucleotide polymorphism (SNP) data, and with very small differences (fixation index [Fst] = 0.01), we obtained results consistently indicating that panmixia does not occur in Japanese eels. The life cycle of the Japanese eel is well-established and the Kuma River is in the center of its habitat. Nevertheless, simple reproductive isolation is not the probable cause of non-panmixia in this species. We propose that the combination of spawning area subdivision, philopatry, and habitat preference/avoidance accounts for the non-panmixia in the Japanese eel population. We named this hypothesis the “reproductive isolation like subset mapping” (RISM) model. This finding may be indicative of the initial stages of sympatric speciation in these eels.

Understanding the mechanism of high-temperature tolerance in cold-freshwater fish is crucial for predicting how certain species will cope with global warming. In this study, we investigated temperature tolerance in masu salmon (Oncorhynchus masou, known in Japan as ‘yamame’), an important aquaculture species. By selective breeding, we developed a group of yamame (F2) with high-temperature tolerance. This group was subjected to a high-temperature tolerance test and divided into two groups: High-temperature tolerant (HT) and non-high-temperature tolerant (NT). RNA was extracted from the gill and adipose fin tissues of each group, and the mRNA expression profiles were analyzed using RNA sequencing. A total of 2893 differentially expressed genes (DEGs) from the gill and 836 from the adipose fin were identified by comparing the HT and NT groups. Functional analyses were then performed to identify associated gene ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The HT group showed a high expression of heat shock protein 70 (HSP70) gene and enriched gene expression in the extracellular matrix (ECM), cell junction, and adhesion pathways in gill tissues compared to the NT group. The HT group also exhibited highly expressed genes in glycolysis and showed lower expression of the genes in the p53 signaling pathway in adipose fin tissues. Taken together, the difference of expression of some genes in the normal condition may be responsible for the difference in heat tolerance between the HT and NT yamame in the heat stress condition.

Metagenomic data have mainly been addressed by showing the composition of organisms based on a small part of a well-examined genomic sequence, such as ribosomal RNA genes and mitochondrial DNAs. On the contrary, whole metagenomic data obtained by the shotgun sequence method have not often been fully analyzed through a homology search because the genomic data in databases for living organisms on earth are insufficient. In order to complement the results obtained through homology-search-based methods with shotgun metagenomes data, we focused on the composition of protein domains deduced from the sequences of genomes and metagenomes, and we utilized them in characterizing genomes and metagenomes, respectively. First, we compared the relationships based on similarities in the protein domain composition with the relationships based on sequence similarities. We searched for protein domains of 325 bacterial species produced using the Pfam database. Next, the correlation coefficients of protein domain compositions between every pair of bacteria were examined. Every pairwise genetic distance was also calculated from 16S rRNA or DNA gyrase subunit B. We compared the results of these methods and found a moderate correlation between them. Essentially, the same results were obtained when we used partial random 100 bp DNA sequences of the bacterial genomes, which simulated raw sequence data obtained from short-read next-generation sequences. Then, we applied the method for analyzing the actual environmental data obtained by shotgun sequencing. We found that the transition of the microbial phase occurred because the seasonal change in water temperature was shown by the method. These results showed the usability of the method in characterizing metagenomic data based on protein domain compositions.