The most critical step in the pearl formation during aquaculture is issued to the proliferation and differentiation of outer epithelial cells of mantle graft into pearl sac. This pearl sac secretes various matrix proteins to produce pearls by a complex physiological process which has not been well-understood yet. Here, we aimed to unravel the genes involved in the development of pearl sac and pearl, and the sequential expression patterns of different shell matrix proteins secreted from the pearl sac during pearl formation by pearl oyster Pinctada fucata using high-throughput transcriptome profiling. Principal component analysis (PCA) showed clearly different gene expression profiles between earlier (before 1 week) and later stages (1 week to 3 months) of grafting. Immune-related genes were highly expressed between 0 h – 24 h (donor dependent) and 48 h – 1 w (host dependent), and in the course of wound healing process pearl sac was developed by two weeks of graft transplantation. Moreover, for the first time, we identified some stem cell marker genes including ABCG2, SOX2, MEF2A, HES1, MET, NRP1, ESR1, STAT6, PAX2, FZD1 and PROM1 that were expressed differentially during the formation of pearl sac. The expression profiling of 192 biomineralization-related genes demonstrated that most of the shell matrix proteins (SMPs) involved in prismatic layer formation were first up-regulated and then gradually down-regulated indicating their involvement in the development of pearl sac and the onset of pearl mineralization. Most of the nacreous layer forming SMPs were up-regulated at 2 weeks after the maturation of pearl sac. Nacrein, MSI7 and shematrin involved in both layer formation were highly expressed during 0 h – 24 h, down-regulated up to 1 week and then up-regulated again after accomplishment of pearl sac formation. Using an RNA-seq approach we unraveled the expression pattern of the key genes involved in the development of pearl sac and pearl as a result of host immune response after grafting. These findings provide valuable information in understanding the molecular mechanism of pearl formation and immune response in P. fucata.

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs (sncRNAs) that perform crucial biological functions in metazoans and defend against transposable elements (TEs) in germ lines. Recently, ubiquitously expressed piRNAs were discovered in soma and germ lines using small RNA sequencing (sRNA-seq) in humans and animals, providing new insights into the diverse functions of piRNAs. However, the role of piRNAs has not yet been fully elucidated, and sRNA-seq studies continue to reveal different piRNA activities in the genome. In this review, we summarize a set of simplified processes for piRNA analysis in order to provide a useful guide for researchers to perform piRNA research suitable for their study objectives. These processes can help expand the functional research on piRNAs from previously reported sRNA-seq results in metazoans. Ubiquitously expressed piRNAs have been discovered in the soma and germ lines in Annelida, Cnidaria, Echinodermata, Crustacea, Arthropoda, and Mollusca, but they are limited to germ lines in Chordata. The roles of piRNAs in TE silencing, gene expression regulation, epigenetic regulation, embryonic development, immune response, and associated diseases will continue to be discovered via sRNA-seq.

Many animals change their body color for visual signaling and environmental adaptation. Some dragonflies show wax-based color change and ultraviolet (UV) reflection, but the biochemical properties underlying the phenomena are totally unknown. Here we investigated the UV-reflective abdominal wax of dragonflies, thereby identifying very long-chain methyl ketones and aldehydes as unique and major wax components. Little wax was detected on young adults, but dense wax secretion was found mainly on the dorsal abdomen of mature males of Orthetrum albistylum and O. melania, and pruinose wax secretion was identified on the ventral abdomen of mature females of O. albistylum and Sympetrum darwinianum. Comparative transcriptomics demonstrated drastic upregulation of the ELOVL17 gene, a member of the fatty acid elongase gene family, whose expression reflected the distribution of very long-chain methyl ketones. Synthetic 2-pentacosanone, the major component of dragonfly’s wax, spontaneously formed light-scattering scale-like fine structures with strong UV reflection, suggesting its potential utility for biomimetics.

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.

The biological process of pearl formation is an ongoing research topic, and a number of genes associated with this process have been identified. However, the involvement of microRNAs (miRNAs) in biomineralization in the pearl oyster, Pinctada fucata, is not well understood. In order to investigate the divergence and function of miRNAs in P. fucata, we performed a transcriptome analysis of small RNA libraries prepared from adductor muscle, gill, ovary, and mantle tissues. We identified 186 known and 42 novel miRNAs in these tissues. Clustering analysis showed that the expression patterns of miRNAs were similar among the somatic tissues, but they differed significantly between the somatic and ovary tissues. To validate the existence of the identified miRNAs, nine known and three novel miRNAs were verified by stem-loop qRT-PCR using U6 snRNA as an internal reference. The expression abundance and target prediction between miRNAs and biomineralization-related genes indicated that miR-1990c-3p, miR-876, miR-9a-3p, and novel-3 may be key factors in the regulatory network that act by controlling the formation of matrix proteins or the differentiation of mineralogenic cells during shell formation in mantle tissue. Our findings serve to further clarify the processes underlying biomineralization in P. fucata.

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.