Transcriptome analysis of Edwardsiella piscicida during intracellular infection reveals excludons are involved with the activation of a mitochondrion-like energy generation program

Phylogenetic evidence suggests a shared ancestry between mitochondria and modern Proteobacteria, a phylum including several genera of intracellular pathogens. Studying these diverse pathogens, particularly during intracellular infection of their hosts, can reveal characteristics potentially representative of the mitochondrial-Proteobacterial ancestor by identifying traits shared with mitochondria. While transcriptomic approaches can provide global insights into intracellular acclimatization by pathogens, they are often limited by excess host RNAs in extracts. Here, we developed a method employing magnetic nanoparticles to enrich RNA from an intracellular Gammaproteobacterium, Edwardsiella piscicida, within zebrafish, Danio rerio, fin fibroblasts, enabling comprehensive exploration of the bacterial transcriptome. Our findings revealed that the intracellular E. piscicida transcriptome reflects a mitochondrion-like energy generation program characterized by the suppression of glycolysis and sugar transport, coupled with upregulation of the tricarboxylic acid (TCA) cycle and alternative import of simple organic acids that directly flux into TCA cycle intermediates or electron transport chain donors. Additionally, genes predicted to be members of excludons, loci of gene pairs antagonistically co-regulated by overlapping antisense transcription, are significantly enriched in the set of all genes with perturbed sense and antisense transcription, suggesting a general but important involvement of excludons with intracellular acclimatization. Notably, genes involved with the activation of the mitochondrion-like energy generation program, specifically with metabolite import and glycolysis, are also members of predicted excludons. Other intracellular Proteobacterial pathogens appear to employ a similar mitochondrion-like energy generation program, suggesting a potentially conserved mechanism for optimized energy acquisition from hosts centered around the TCA cycle.IMPORTANCEPhylogenetic evidence suggests that mitochondria and Proteobacteria, a phylum encompassing various intracellular pathogens, share a common ancestral lineage. In this study, we developed a novel method employing magnetic nanoparticles to explore the transcriptome of an aquatic Gammaproteobacterium, Edwardsiella piscicida, during intracellular infection of host cells. We show that the strategy E. piscicida uses to generate energy strikingly mirrors the function of mitochondria-energy generators devoid of glycolytic processes. Notably, several implicated genes are members of excludons-gene pairs antagonistically co-regulated by overlapping antisense transcription. Other intracellular Proteobacterial pathogens appear to adopt a similar mitochondrion-like energy generation program, indicating a possibly conserved strategy for optimized energy acquisition from hosts centered around the tricarboxylic acid cycle.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

The highly developed symbiotic system between the solar-powered nudibranch Pteraeolidia semperi and Symbiodiniacean algae

The intricate coexistence of Symbiodiniacean algae with a diverse range of marine invertebrates underpins the flourishing biodiversity observed within coral reef ecosystems. However, the breakdown of Symbiodiniaceae-host symbiosis endangers these ecosystems, necessitating urgent study of the symbiotic mechanisms. The symbiosis between nudibranchs and Symbiodiniaceae has been identified as an efficacious model for examining these mechanisms, yet a comprehensive understanding of their histological structures and cellular processes remains elusive. A meticulous histological exploration of the nudibranch Pteraeolidia semperi, employing optical, fluorescence, and electron microscopy, has revealed fine tubules extending to the body surface, with associated epithelial cells having been shown to adeptly encapsulate Symbiodiniaceae intracellularly. By tracing the stages of the "bleaching" in nudibranchs, it was inferred that algal cells, translocated via the digestive gland, are directly phagocytosed and expelled by these epithelial cells. Collectively, these insights contribute substantially to the scholarly discourse on critical marine symbiotic associations.

The complete mitochondrial genome of Spurilla braziliana MacFarland 1909 (Nudibranchia, Aeolidiidae)

Spurilla braziliana MacFarland 1909 is a morphologically diverse nudibranch found in the Pacific and Western Atlantic. The complete mitochondrial genome of S. braziliana has been constructed using next-generation sequencing technology. The mitochondrial genome is 14,291 bp and contains 13 protein-coding genes, 2 rRNA genes, and 23 tRNA genes. Molecular phylogenetic analysis using the maximum likelihood method revealed that S. braziliana is included in the superfamily Aeolidioidea and forms a monophyletic group with Berghia stephanieae, a nudibranch of the family Aeolidiidae. This study reinforces existing taxonomic insights and provides a basis for further molecular phylogenetic analysis.

Draft genome sequencing and assembly of Risso's dolphin, Grampus griseus

The Risso's dolphin (Grampus griseus) is one of the migratory marine mammals and they have commonly dispersed in tropical and temperate seas. It is a least concerned species in the IUCN red list of threatened species. However, their population size and factors affecting their population structure are unknown. Due to the wide distribution of this species, their populations might be genetically stable and less structured. To support genetic studies on dolphins and other marine mammals, we assembled the draft genome of Risso's dolphin that was found in Japan. The tissue samples were used to extract high molecular DNA and subjected to sequencing by Illumina HiSeq X, Oxford Nanopore MinION, and Bionano Saphyr. The assembled hybrid genome was 75.9% of complete eukaryotic BUSCOs and the genome size was 2.256 Gb with 2.042 Mb of scaffold N50. De novo assembly of this genome by Bionano Saphyr recovered 2.036 Gb total genome map length and structural variations. The gene structures of this draft genome were identified by BRAKER2, and 9947 genes were recovered. The data will be useful for future studies of cetaceans.

Evolution of nacre- and prisms-related shell matrix proteins in the pen shell, Atrina pectinata

The molluscan shell is a good model for understanding the mechanisms underlying biomineralization. It is composed of calcium carbonate crystals and many types of organic molecules, such as the matrix proteins, polysaccharides, and lipids. The pen shell Atrina pectinata (Pterioida, Pinnidae) has two shell microstructures: an outer prismatic layer and an inner nacreous layer. Similar microstructures are well known in pearl oysters (Pteriidae), such as Pinctada fucata, and many kinds of shell matrix proteins (SMPs) have been identified from their shells. However, the members of SMPs that consist of the nacreous and prismatic layers of Pinnidae bivalves remain unclear. In this study, we identified 114 SMPs in the nacreous and prismatic layers of A. pectinata, of which only seven were found in both microstructures. 54 of them were found to bind calcium carbonate. Comparative analysis of nine molluscan shell proteomes showed that 69 of 114 SMPs of A. pectinata were found to have sequential similarity with at least one or more SMPs of other molluscan species. For instance, nacrein, tyrosinase, Pif/BMSP-like, chitinase (CN), chitin-binding proteins, CD109, and Kunitz-type serine proteinase inhibitors are widely shared among bivalves and gastropods. Our results provide new insights for understanding the complex evolution of SMPs related to nacreous and prismatic layer formation in the pteriomorph bivalves.

Age-Associated Different Transcriptome Profiling in Zebrafish and Rats: an Insight into the Diversity of Vertebrate Aging

Most mammals, including humans, show obvious aging phenotypes, for example, loss of tissue plasticity and sarcopenia. In this regard, fish can be attractive models to study senescence because of their unique aging characteristics. The lifespan of fish varies widely, and several species can live for over 200 years. Moreover, some fish show anti-aging features and indeterminate growth throughout their life. Therefore, exploring the aging mechanism in fish could provide new insights into vertebrate aging. To this end, we conducted RNA sequencing (RNA-seq) assays for various organs and growth stages of zebrafish and compared the data with previously published RNA-seq data of rats. Age-associated differentially expressed genes (DEGs) for all zebrafish tissue samples reveal the upregulation of circadian genes and downregulation of hmgb3a. On one hand, a comparative analysis of DEG profiles associated with aging between zebrafish and rats identifies upregulation of circadian genes and downregulation of collagen genes as conserved transcriptome changes. On the other hand, in zebrafish, upregulation of autophagy-related genes in muscles and AP-1 transcription factor genes in various tissues is observed, which may imply fish-specific anti-aging characteristics. Consistent with our knowledge of mammalian aging, DEG profiles related to tissue senescence are observed in rats. We also detect age-associated downregulation of muscle homeostasis and differentiation-related genes in zebrafish gills, indicating a fish-specific senescence phenotype. Our results indicate both common and different aging profiles between fish and mammals, which could be used for future translational research.

Zebrafish Danio rerio myotomal muscle structure and growth from a spatial transcriptomics perspective

Fish exhibit different muscle structures and growth characteristics compared with mammals. We used a spatial transcriptomics approach and examined myotomal muscle sections from zebrafish. Adult muscles were divided into eight regions according to spatial gene expression characteristics. Slow muscle was located in the wedge-shaped region near the lateral line and at the base of the dorsal fin, intermediate muscle was located in a ribbon-shaped region adjacent to slow muscle, and fast muscle was located in the deep region of the trunk, surrounded by intermediate muscle; the interior of fast muscle was further divided into 6 parts by their transcriptomic features. Combined analysis of adult and larval data revealed that adult muscles contain specific regions similar to larval muscles. These regions showed active myogenesis and a high expression of genes associated with muscle hyperplasia. This is the first study to apply spatial transcriptomics to fish myotomal muscle structure and growth.

Induction of endoplasmic reticulum stress markers in an acromegaly model

Acromegaly is a growth hormone (GH) excess pathological condition in humans. Acromegaly is associated with somatic disfigurement and a wide range of systemic manifestations such as arthritis, neuropathy, carpal tunnel syndrome, reproductive disorders, metabolic disorders, and gastrointestinal complications. The influence of excess GH on the cellular level could aid in understanding the root causes of acromegaly-related health complications. Previously, we found that GH excess induces DNA damage to somatic cells and reduces the stem cells number and causes premature aging. In this study, an in-depth analysis of the acromegaly RNAseq data revealed the disruption of important biological cellular processes. Gene set enrichment analysis, heatmap, and enrichment analysis of acromegaly RNAseq data revealed induction of endoplasmic reticulum (ER) stress markers in various organs. Interestingly, the induction of ER stress was even more apparent than in aged zebrafish. Splicing of box-binding protein-1 (XBP1) mRNA is a hallmark of ER stress. Therefore, we quantified spliced XBP1 mRNA in different organs of our acromegaly model. Thus, our study emphasizes the importance of ER stress in GH oversecretion, which is important for understanding the health complications of acromegaly.

Exosome-derived small non-coding RNAs reveal immune response upon grafting transplantation in Pinctada fucata (Mollusca)

Exosomes, a subset of small extracellular vesicles, carry various nucleic acids, proteins, lipids, amino acids and metabolites. They function as a mode of intercellular communication and molecular transfer. Exosome cargo molecules, including small non-coding RNAs (sncRNAs), are involved in the immune response in various organisms. However, the role of exosome-derived sncRNAs in immune responses in molluscs remains unclear. Here, we aimed to reveal the sncRNAs involved in the immune response during grafting transplantation by the pearl oyster Pinctada fucata. Exosomes were successfully extracted from the P. fucata haemolymph during graft transplantation. Abundant microRNAs (miRNAs) and PIWI-interacting RNAs (piRNAs) were simultaneously discovered in P. fucata exosomes by small RNA sequencing. The expression patterns of the miRNAs and piRNAs at the grafting and initial stages were not substantially different, but varied significantly between the initial and later stages. Target prediction and functional analysis indicate that these miRNAs and piRNAs are related to immune response upon grafting transplantation, whereas piRNAs may also be associated with transposon silencing by targeting with genome transposon elements. This work provides the basis for a functional understanding of exosome-derived sncRNAs and helps to gain further insight into the PIWI/piRNA pathway function outside of germline cells in molluscs.

Construction of a chromosome-level Japanese stickleback species genome using ultra-dense linkage analysis with single-cell sperm sequencing

It is still difficult to construct the genomes of higher organisms as their genome sequences must be extended to the length of the chromosome by linkage analysis. In this study, we attempted to provide an innovative alternative to conventional linkage analysis by devising a method to genotype sperm using 10× Genomics single-cell genome sequencing libraries to generate a linkage map without interbreeding individuals. A genome was assembled using sperm from the Japanese stickleback Gasterosteus nipponicus, with single-cell genotyping yielding 1 864 430 very dense hetero-SNPs and an average coverage per sperm cell of 0.13×. In total, 1665 sperm were used, which is an order of magnitude higher than the number of recombinations used for conventional linkage analysis. We then improved the linkage analysis tool scaffold extender with low depth linkage analysis (SELDLA) to analyze the data according to the characteristics of the single-cell genotyping data. Finally, we were able to determine the chromosomal location (97.1%) and orientation (64.4%) of the contigs in the 456 Mb genome of G. nipponicus, sequenced using nanopores. This method promises to be a useful tool for determining the genomes of non-model organisms for which breeding systems have not yet been established by linkage analysis.

Selection of a reference gene for studies on lipid-related aquatic adaptations of toothed whales (Grampus griseus)

Toothed whales are one group of marine mammals that has developed special adaptations, such as echolocation for predation, to successfully live in a dynamic aquatic environment. Their fat metabolism may differ from that of other mammals because toothed whales have acoustic fats. Gene expression in the metabolic pathways of animals can change with respect to their evolution and environment. A real-time quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for studying the relative expressions of genes. However, since the accuracy of RT-qPCR data is totally dependent on the reference gene, the selection of the reference gene is an essential step. In this study, 10 candidate reference genes (ZC3H10, FTL, LGALS1, RPL27, GAPDH, FTH1, DCN, TCTP, NDUS5, and UBIM) were initially tested for amplification efficiency using RT-qPCR. After excluding DCN, the remaining nine genes, which are nearly 100% efficient, were selected for the gene stability analysis. Stable reference genes across eight different fat tissue, liver, and muscle samples from Grampus griseus were identified by four algorithms, which were provided in Genorm, NormFinder, BestKeeper, and Delta CT. Finally, a RefFinder comprehensive ranking was performed based on the stability values, and the nine genes were ranked as follows: LGALS1 > FTL > GAPDH > ZC3H10 > FTH1 > NDUS5 > TCTP > RPL27 > UBIM. The LGALS1 and FTL genes were identified as the most stable novel reference genes. The third-ranked gene, GAPDH, is a well-known housekeeping gene for mammals. Ultimately, we suggest the use of LGALS1 as a reliable novel reference gene for genomics studies on the lipid-related aquatic adaptations of toothed whales.

Metagenomic analysis provides functional insights into seasonal change of a non-cyanobacterial prokaryotic community in temperate coastal waters

The taxonomic compositions of marine prokaryotic communities are known to follow seasonal cycles, but functional metagenomic insights into this seasonality is still limited. We analyzed a total of 22 metagenomes collected at 11 time points over a 14-month period from two sites in Sendai Bay, Japan to obtain seasonal snapshots of predicted functional profiles of the non-cyanobacterial prokaryotic community. Along with taxonomic composition, functional gene composition varied seasonally and was related to chlorophyll a concentration, water temperature, and salinity. Spring phytoplankton bloom stimulated increased abundances of putative genes that encode enzymes in amino acid metabolism pathways. Several groups of functional genes, including those related to signal transduction and cellular communication, increased in abundance during the mid- to post-bloom period, which seemed to be associated with a particle-attached lifestyle. Alternatively, genes in carbon metabolism pathways were generally more abundant in the low chlorophyll a period than the bloom period. These results indicate that changes in trophic condition associated with seasonal phytoplankton succession altered the community function of prokaryotes. Our findings on seasonal changes of predicted function provide fundamental information for future research on the mechanisms that shape marine microbial communities.

Regulation of the Expression of the Myosin Heavy Chain (MYH) Gene myh14 in Zebrafish Development

The human sarcomeric myosin heavy chain gene MYH14 contains an intronic microRNA, miR-499. Our previous studies demonstrated divergent genomic organization and expression patterns of myh14/miR-499 among teleosts; however, the regulatory mechanism is partly known. In this study, we report the regulation of myh14 expression in zebrafish, Danio rerio. Zebrafish myh14 has three paralogs, myh14-1, myh14-2, and myh14-3. Detailed promoter analysis suggested that a 5710-bp 5'-flanking region of myh14-1 and a 5641-bp region of myh14-3 contain a necessary regulatory region to recapitulate specific expression during embryonic development. The 5'-flanking region of zebrafish myh14-1 and its torafugu ortholog shared two distal and a single proximal conserved region. The two distal conserved regions had no effect on zebrafish myh14-1 expression, in contrast to torafugu expression, suggesting an alternative regulatory mechanism among the myh14 orthologs. Comparison among the 5'-flanking regions of the myh14 paralogs revealed two conserved regions. Deletion of these conserved regions significantly reduced the promoter activity of myh14-3 but had no effect on myh14-1, indicating different cis-regulatory mechanisms of myh14 paralogs. Loss of function of miR-499 resulted in a marked reduction in slow muscle fibers in embryonic development. Our study identified different cis-regulatory mechanisms controlling the expression of myh14/miR-499 and an indispensable role of miR-499 in muscle fiber-type specification in zebrafish.

Another polymorphic mitochondrial genome of Grampus griseus and phylogeny of family Delphinidae

Risso’s dolphin (Grampus griseus Cuvier, 1812) is the only species of genus Grampus and a cosmopolitan marine inhabitant. Here, we report a polymorphic complete mitochondrial genome of G. griseus. The size of the total mitochondrial genome was 16,386 bp in length and contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region. 37 single nucleotide polymorphic sites (SNPs) were identified compared to the references. Based on the available total mitochondrial dolphin genomes’ phylogenetics, G. griseus has formed a clade with 0.1415 distance, sister to the following species of the subfamily Globicephalinae and the taxonomy of Orcinus orca still needs further investigations.

Discovery and functional understanding of MiRNAs in molluscs: a genome-wide profiling approach

Small non-coding RNAs play a pivotal role in gene regulation, repression of transposable element and viral activity in various organisms. Among the various categories of these small non-coding RNAs, microRNAs (miRNAs) guide post-translational gene regulation in cellular development, proliferation, apoptosis, oncogenesis, and differentiation. Here, we performed a genome-wide computational prediction of miRNAs to improve the understanding of miRNA observation and function in molluscs. As an initial step, hundreds of conserved miRNAs were predicted in 35 species of molluscs through genome scanning. Afterwards, the miRNAs' population, isoforms, organization, and function were characterized in detail. Furthermore, the key miRNA biogenesis factors, including AGO2, DGCR8, DICER, DROSHA, TRABP2, RAN, and XPO5, were elucidated based on homologue sequence searching. We also summarized the miRNAs' function in biomineralization, immune and stress response, as well as growth and development in molluscs. Because miRNAs play a vital role in various lifeforms, this study will provide insight into miRNA biogenesis and function in molluscs, as well as other invertebrates.

An Acromegaly Disease Zebrafish Model Reveals Decline in Body Stem Cell Number along with Signs of Premature Aging

In our previous publication, it was shown that growth hormone (GH) excess in acromegaly affects the cell integrity of somatic cells through increased DNA damage throughout the body and impaired DNA repair pathways. Acromegaly is a hormone disorder pathological condition that develops as a result of growth hormone over-secretion from the pituitary gland. We produced a zebrafish acromegaly model to gain a better understanding of the excess GH effects at the cellular level. Here we show that the acromegaly zebrafish model progressively reduced the number of stem cells in different organs and increased oxidative stress in stem cells. Importantly, the decline in the stem cells was even more apparent than in aged fish. The controversy and debate over the use of GH as an anti-aging therapy have been going on for several years. In this study, excess GH induced aging signs such as increased senescence-associated (SA)-β-galactosidase staining of abdominal skin and similarity of the pattern of gene expression between aged and acromegaly zebrafish. Thus, this study highlights the role of excess GH in acromegaly stem cells.

Trivalent Iron Is Responsible for the Yellow Color Development in the Nacre of Akoya Pearl Oyster Shells

The gold and cream colors of cultured Akoya pearls, as well as natural yellow nacre of pearl oyster shells, are thought to arise from intrinsic yellow pigments. While the isolation of the yellow pigments has been attempted using a large amount of gold pearls, the substance concerned is still unknown. We report here on the purification and characterization of yellow pigments from the nacre of Akoya pearl oyster shells. Two yellow components, YC1 and YC2, were isolated from the HCl-methanol (HCl-MeOH) extract from nacreous organic matrices obtained by decalcification of the shells with ethylenediaminetetraacetic acid (EDTA). Energy-dispersive X-ray and infrared spectroscopy analyses suggested that YC1 and YC2 precipitated under basic conditions are composed of Fe-containing inorganic and polyamide-containing organic compounds, respectively. YC1 solubilized under acidic conditions exhibited positive reactions to KSCN and K₄[Fe(CN)₆] reagents, showing the same ultraviolet-visible absorption spectrum as those of Fe(III)-containing compounds. In addition, X-ray absorption fine structure analysis supported the compound in the form of Fe(III). The total amount of Fe was approximately 2.6 times higher in the yellow than white nacre, and most Fe was fractionated into the EDTA-decalcifying and HCl-MeOH extracts. These results suggest that Fe(III) coordinated to EDTA-soluble and insoluble matrix compounds are mainly associated with yellow color development not only in the Akoya pearl oyster shells but also in the cultured Akoya pearls.

Identification and Characterization of microRNAs and Their Predicted Functions in Biomineralization in the Pearl Oyster (Pinctada fucata)

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.

Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

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

Piwi-interacting RNA (piRNA) expression patterns in pearl oyster (Pinctada fucata) somatic tissues

Piwi-interacting RNAs (piRNAs) belong to a recently discovered class of small non-coding RNAs whose best-understood function is repressing transposable element activity. Most piRNA studies have been conducted on model organisms and little is known about piRNA expression and function in mollusks. We performed high-throughput sequencing of small RNAs extracted from the mantle, adductor muscle, gill, and ovary tissues of the pearl oyster, Pinctada fucata. RNA species with sequences of approximately 30 nt were widely expressed in all tissues. Uridine at the 5' terminal and protection from β-elimination at the 3' terminal suggested that these were putative piRNAs. A total of 18.0 million putative piRNAs were assigned to 2.8 million unique piRNAs, and 35,848 piRNA clusters were identified. Mapping to the reference genome showed that 25% of the unique piRNAs mapped to multiple tandem loci on the scaffold. Expression patterns of the piRNA clusters were similar within the somatic tissues, but differed significantly between the somatic and gonadal tissues. These findings suggest that in pearl oysters piRNAs have important and novel functions beyond those in the germ line.

Effects of Nutrition Therapy in Older Stroke Patients Undergoing Rehabilitation: A Systematic Review and Meta-Analysis

OBJECTIVE

To systematically review evidence on the effects of nutrition therapy in older stroke patients undergoing rehabilitation and identify its effectiveness using meta-analysis.

METHODS

PubMed (MEDLINE), EMBASE (via Dialog), Cochrane Central Register of Controlled Trial, World Health Organization International Clinical Trials Registry Platform and Ichu-shi Web were searched for relevant articles. Randomized controlled trials investigating the effects of nutrition therapy compared to control interventions in older stroke patients undergoing rehabilitation were considered eligible. The primary outcome was activities of daily living (ADL), and secondary outcomes were all-cause mortality, infections, pneumonia incidence, disability level, walking ability, fall, stroke recurrence, and quality of life. The risk of bias of each trial was assessed using the Cochrane Collaboration Tool, and the quality of the body of evidence was assessed using the Grading of Recommendations Assessment, Development and Evaluation approach.

RESULTS

Eight randomized controlled trials with a total of 5484 participants were included in the meta-analysis. The meta-analysis for ADL showed no significant effects (mean difference, 4.16; 95% confidence interval [CI], -0.88 to 9.20; I2=53%, low-quality evidence). The meta-analyses for secondary outcomes revealed a significant effect of reduced infections (risk ratio, 0.65; 95% CI, 0.51 to 0.84; I2=0%; low-quality evidence), with no significant effects on the other outcomes.

CONCLUSION

Nutrition therapy had no statistically significant effect on ADL. However, it reduced the incidence of infections. More high-quality trials are warranted to clarify the effects of nutrition therapy in older stroke patients undergoing rehabilitation.

Transcriptomic analysis of immunoglobulin novel antigen receptor (IgNAR) heavy chain constant domains of brownbanded bamboo shark (Chiloscyllium punctatum)

Cartilaginous fish are the evolutionarily oldest group of animals which possess antibodies, T cell receptors and major histocompatibility complex (MHC). The immunoglobulin novel antigen receptor (IgNAR) found in cartilaginous fish is a heavy chain homodimer which lacks light chain. The presence of non-canonical cysteine molecules and lack of CDR2 region make it more significant. To synthesize active binding domains based on variable region of IgNAR (VNAR), knowledge on the constant region dynamics play a significant role. The IgNAR exhibit species variations in its primary sequence features; hence, this study was conducted to determine the IgNAR heavy chain constant domain of the brownbanded bamboo shark (Chiloscyllium punctatum). Peripheral blood leukocytes (PBL) isolated from adult bamboo sharks were used to synthesize a cDNA library. A total of four billion residues of two million sequences (average length 218.41 bp) were obtained. Assembled sequences were aligned with published cartilaginous fish IgNAR constant region sequences. Transcriptome analysis revealed two distinct types of IgNAR in the brownbanded bamboo shark. Also, constant-1 domain sequences displayed 13 unique sequences which may reflect the least number of IgNAR gene clusters. The phylogenetic analysis revealed the closest relationship with the nurse shark (Ginglymostoma cirratum) followed by the wobbegong shark (Orectolobus maculatus) which belong to the same order Orectolobiformes. Analysis of the constant domains of the brownbanded bamboo shark IgNAR revealed an evolutionarily conserved nature and this knowledge can be used to design primers for VNAR cloning. Furthermore, knowledge on the structural features in IgNAR constant domains that increase the stability could be useful in the process of stabilizing human immunoglobulins.

Comparison of Two Pearl Sacs Formed in the Same Recipient Oyster with Different Genetic Background Involved in Yellow Pigmentation in Pinctada fucata

Color is one of the most important factors determining the commercial value of pearls. Pinctada fucata is a well-known pearl oyster producing high-quality Akoya pearls. Phenotypic variation in amount of yellow pigmentation produces white and yellowish pearls. It has been reported that polymorphism of yellow pigmentation of Akoya pearls is genetically regulated, but the responsible gene(s) has remained unknown. Here, we prepared pearl sac pairs formed in the same recipient oyster but coming from donor oysters that differ in their color. These two pearl sacs produced pearls with different yellowness even in the same recipient oyster. Yellow tone of produced pearls was consistent with shell nacre color of donor oysters from which mantle grafts were prepared, indicating that donor oysters strongly contribute to the yellow coloration of Akoya pearls. We also conducted comparative RNA-seq analysis and retrieved several candidate genes involved in the pearl coloration. Whole gene expression patterns of pair sacs were not grouped by pearl color they produced, but grouped by recipient oysters in which they were grown, suggesting that the number of genes involved in the yellow coloration is quite small, and that recipient oyster affects gene expression of the majority of genes in the pearl sac.

Promoter analysis of the fish gene of slow/cardiac-type myosin heavy chain implicated in specification of muscle fiber types

Vertebrate skeletal muscles consist of heterogeneous tissues containing various types of muscle fibers, where specification of the fiber type is crucial for muscle development. Fish are an attractive experimental model to study the mechanisms of such fiber type specification because of the separated localization of slow and fast muscles in the trunk myotome. We examined regulation of expression of the torafugu gene of slow/cardiac-type myosin heavy chain, MYH _M5 , and isolated an operational promoter in order to force its tissue-specific expression across different fish species via the transgenic approach in zebrafish and medaka. This promoter activity was observed in adaxial cell-derived superficial slow muscle fibers under the control of a hedgehog signal. We also uncovered coordinated expression of MYH _M5 and Sox6b, which is an important transcriptional repressor for specification of muscle fiber types and participates in hedgehog signaling. Sequence comparison in the 5'-flanking region identified three conserved regions, CSR1-CSR3, between torafugu MYH _M5 and its zebrafish ortholog. Analysis of deletion mutants showed that CSR1 significantly stimulates gene expression in slow muscle fibers. In contrast, deletion of CSR3 resulted in ectopic expression of a reporter gene in fast muscle fibers. CSR3 was found to contain a putative Sox family protein-binding site. These results indicate that the dual mechanism causing inhibition in fast muscle fibers and activation in slow muscle fibers is essential for slow muscle fiber-specific gene expression in fish.

Novel Isoforms of N16 and N19 Families Implicated for the Nacreous Layer Formation in the Pearl Oyster Pinctada fucata

Although a wide variety of proteins and genes possibly related to the shell formation in bivalve have been identified, their functions have been only partially approved. We have recently performed deep sequencing of expressed sequence tags (ESTs) from the pearl oyster Pinctada fucata using a next-generation sequencer, identifying a dozen of novel gene candidates which are possibly associated with the nacreous layer formation. Among the ESTs, we focused on three novel isoforms (N16-6, N16-7, and N19-2) of N16 and N19 families with reference to five known genes in the families and determined the full-length cDNA sequences of these isoforms. Reverse transcription-polymerase chain reaction revealed that N16-6 was expressed in gill, gonad, adductor muscle, and mantle, whereas N16-7 exclusively in mantle. N19-2 was expressed in all tissues examined. In situ hybridization demonstrated their regional expression in mantle and pearl sac, which well corresponded to those shown by EST analysis previously reported. Shells in the pearl oyster injected with dsRNAs of N16-7 and N19-2 showed abnormal surface appearance in the nacreous layer. Taken together, novel isoforms in N16 and N19 families shown in this study are essential to form the nacreous layer.

Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes)

Next-generation sequencing enables genome-wide genotyping of a large population and further facilitates the construction of a genetic linkage map. Low-coverage whole-genome sequencing has been employed for genetic linkage map construction in several species. However, this strategy generally requires available high-quality reference genomes and/or designed inbred pedigree lines, which restrict the scope of application for non-model and unsequenced species. Here, using torafugu (Takifugu rubripes) as a test model, we propose a new strategy for ultrahigh-density genetic linkage map construction using low-coverage whole-genome sequencing of a haploid/doubled haploid (H/DH) population without above requirements. Low-coverage (≈1×) whole-genome sequencing data of 165 DH individuals were used for de novo assembly and further performed single nucleotide polymorphisms (SNPs) calling, resulting in the identification of 1,070,601 SNPs. Based on SNP genotypes and de novo assembly, genotypes were associated with short DNA segments and an ultrahigh-density linkage map was constructed containing information of 802,277 SNPs in 3090 unique positions. Comparative analyses showed near-perfect concordance between the present linkage map and the latest published torafugu genome (FUGU5). This strategy would facilitate ultrahigh-density linkage map construction in various sexually reproducing organisms for which H/DH populations can be generated.

Mucocutaneous inflammation in the Ikaros Family Zinc Finger 1-keratin 5-specific transgenic mice

BACKGROUND

Our genomewide association study documented an association between cold medicine-related Stevens-Johnson syndrome/toxic epidermal necrolysis (CM-SJS/TEN) and Ikaros Family Zinc Finger 1 (IKZF1). Few studies examined biological and pathological functions of IKZF1 in mucosal immunity. We hypothesized that IKZF1 contributes to the mucocutaneous inflammation.

METHODS

Human skin and conjunctival tissues were obtained for immunohistological studies. Primary human conjunctival epithelial cells (PHCjECs) and adult human epidermal keratinocytes (HEKa) also used for gene expression analysis. We also generated K5-Ikzf1-EGFP transgenic mice (Ikzf1 Tg) by introducing the Ik1 isoform into cells expressing keratin 5, which is expressed in epithelial tissues such as the epidermis and conjunctiva, and then examined them histologically and investigated gene expression of the epidermis. Moreover, Ikzf1 Tg were induced allergic contact dermatitis.

RESULTS

We found that human epidermis and conjunctival epithelium expressed IKZF1, and in PHCjECs and HEKa, the expression of IKZF1 mRNA was upregulated by stimulation with polyI:C, a TLR3 ligand. In Ikzf1 Tg, we observed dermatitis and mucosal inflammation including the ocular surface. In contact dermatitis model, inflammatory infiltrates in the skin of Ikzf1 Tg were significantly increased compared with wild type. Microarray analysis showed that Lcn2, Adh7, Epgn, Ifi202b, Cdo1, Gpr37, Duoxa1, Tnfrsf4, and Enpp5 genes were significantly upregulated in the epidermis of Ikzf1 Tg compared with wild type.

CONCLUSION

Our findings support the hypothesis that Ikaros might participate in mucocutaneous inflammation.

Evolution and Distribution of Teleost myomiRNAs: Functionally Diversified myomiRs in Teleosts

Myosin heavy chain (MYH) genes belong to a multigene family, and the regulated expression of each member determines the physiological and contractile muscle properties. Among these, MYH6, MYH7, and MYH14 occupy unique positions in the mammalian MYH gene family because of their specific expression in slow/cardiac muscles and the existence of intronic micro(mi) RNAs. MYH6, MYH7, and MYH14 encode miR-208a, miR-208b, and miR-499, respectively. These MYH encoded miRNAs are designated as myomiRs because of their muscle-specific expression and functions. In mammals, myomiRs and host MYHs form a transcription network involved in muscle fiber-type specification; thus, genomic positions and expression patterns of them are well conserved. However, our previous studies revealed divergent distribution and expression of MYH14/miR-499 among teleosts, suggesting the unique evolution of myomiRs and host MYHs in teleosts. Here, we examined distribution and expression of myomiRs and host MYHs in various teleost species. The major cardiac MYH isoforms in teleosts are an intronless gene, atrial myosin heavy chain (amhc), and ventricular myosin heavy chain (vmhc) gene that encodes an intronic miRNA, miR-736. Phylogenetic analysis revealed that vmhc/miR-736 is a teleost-specific myomiR that differed from tetrapoda MYH6/MYH7/miR-208s. Teleost genomes also contain species-specific orthologs in addition to vmhc and amhc, indicating complex gene duplication and gene loss events during teleost evolution. In medaka and torafugu, miR-499 was highly expressed in slow/cardiac muscles whereas the expression of miR-736 was quite low and not muscle specific. These results suggest functional diversification of myomiRs in teleost with the diversification of host MYHs.

The FGF23/KLOTHO Regulatory Network and Its Roles in Human Disorders

The functions of Klotho (KL) are multifaceted and include the regulation of aging and mineral metabolism. It was originally identified as the gene responsible for premature aging-like symptoms in mice and was subsequently shown to function as a coreceptor in the fibroblast growth factor (FGF) 23 signaling pathway. The discovery of KL as a partner for FGF23 led to significant advances in understanding of the molecular mechanisms underlying phosphate and vitamin D metabolism, and simultaneously clarified the pathogenic roles of the FGF23 signaling pathway in human diseases. These novel insights led to the development of new strategies to combat disorders associated with the dysregulated metabolism of phosphate and vitamin D, and clinical trials on the blockade of FGF23 signaling in X-linked hypophosphatemic rickets are ongoing. Molecular and functional insights on KL and FGF23 have been discussed in this review and were extended to how dysregulation of the FGF23/KL axis causes human disorders associated with abnormal mineral metabolism.

Bivalve-specific gene expansion in the pearl oyster genome: implications of adaptation to a sessile lifestyle

INTRODUCTION

Bivalve molluscs have flourished in marine environments, and many species constitute important aquatic resources. Recently, whole genome sequences from two bivalves, the pearl oyster, Pinctada fucata, and the Pacific oyster, Crassostrea gigas, have been decoded, making it possible to compare genomic sequences among molluscs, and to explore general and lineage-specific genetic features and trends in bivalves. In order to improve the quality of sequence data for these purposes, we have updated the entire P. fucata genome assembly.

RESULTS

We present a new genome assembly of the pearl oyster, Pinctada fucata (version 2.0). To update the assembly, we conducted additional sequencing, obtaining accumulated sequence data amounting to 193× the P. fucata genome. Sequence redundancy in contigs that was caused by heterozygosity was removed in silico, which significantly improved subsequent scaffolding. Gene model version 2.0 was generated with the aid of manual gene annotations supplied by the P. fucata research community. Comparison of mollusc and other bilaterian genomes shows that gene arrangements of Hox, ParaHox, and Wnt clusters in the P. fucata genome are similar to those of other molluscs. Like the Pacific oyster, P. fucata possesses many genes involved in environmental responses and in immune defense. Phylogenetic analyses of heat shock protein70 and C1q domain-containing protein families indicate that extensive expansion of genes occurred independently in each lineage. Several gene duplication events prior to the split between the pearl oyster and the Pacific oyster are also evident. In addition, a number of tandem duplications of genes that encode shell matrix proteins are also well characterized in the P. fucata genome.

CONCLUSIONS

Both the Pinctada and Crassostrea lineages have expanded specific gene families in a lineage-specific manner. Frequent duplication of genes responsible for shell formation in the P. fucata genome explains the diversity of mollusc shell structures. These duplications reveal dynamic genome evolution to forge the complex physiology that enables bivalves to employ a sessile lifestyle in the intertidal zone.

Three percent diquafosol ophthalmic solution as an additional therapy to existing artificial tears with steroids for dry-eye patients with Sjögren's syndrome

PURPOSE

To investigate the long-term results of 3% diquafosol ophthalmic solution as an alternative therapy to existing ophthalmic solutions, including topical immunosuppression, for the treatment of dry eye in patients with Sjögren's syndrome.

METHODS

This study involved 14 female dry-eye patients (mean age: 62.4 years) with Sjögren's syndrome who insufficiently responded to their current therapy. In all patients, 3% diquafosol ophthalmic solution was administered six times daily for 12 months in substitution for artificial tears and sodium hyaluronate ophthalmic solution. Their use of corticosteroid eye drops remained unchanged from that prior to the treatment with diquafosol sodium. The subjective symptoms assessed, and ocular signs including tear meniscus radius and the tear film breakup time, and ocular-surface epithelial damage score were examined at 1, 2, 3, 4, 5, 6, 9, and 12 months after initiating treatment.

RESULTS

Among the subjective symptoms, significant improvement was obtained in dryness at 2 months post treatment, in eye fatigue at 1, 2, 3, 4, and 12 months post treatment, and in pain at 1, 2, 6, and 12 months post treatment. Difficulty in opening the eye, foreign body sensation, and redness were also significantly ameliorated at various time-points. The tear meniscus radius and the tear film breakup time were significantly improved throughout the observation period, and the corneal epithelial staining scores were significantly decreased at 3 months post treatment.

CONCLUSIONS

In dry-eye patients with Sjögren's syndrome, treatment with 3% diquafosol ophthalmic solution improved both symptoms and signs, and that effectiveness was maintained for 12 months.

Deep sequencing, profiling and detailed annotation of microRNAs in Takifugu rubripes

Background

microRNAs (miRNAs) in fish have not been as extensively studied as those in mammals. The fish species Takifugu rubripes is an intensively studied model organism whose genome has been sequenced. The T. rubripes genome is approximately eight times smaller than the human genome, but has a similar repertoire of protein-coding genes. Therefore, it is useful for identifying non-coding genes, including miRNA genes. To identify miRNA expression patterns in different organs of T. rubripes and give fundamental information to aid understanding of miRNA populations in this species, we extracted small RNAs from tissues and performed deep sequencing analysis to profile T. rubripes miRNAs. These data will be of assistance in functional studies of miRNAs in T. rubripes.

Results

After analyzing a total of 139 million reads, we found miRNA species in nine tissues (fast and slow muscles, heart, eye, brain, intestine, liver, ovaries, and testes). We identified 1420 known miRNAs, many of which were strongly expressed in certain tissues with expression patterns similar to those described for other animals in previous reports. Most miRNAs were expressed in tissues other than the ovaries or testes. However, some miRNA families were highly abundant in the gonads, but expressed only at low levels in somatic tissue, suggesting specific function in germ cells. The most abundant isomiRs (miRNA variants) of many miRNAs had identical sequences in the 5′ region. However, isomiRs of some miRNAs, including fru-miR-462-5p, varied in the 5′ region in some tissues, suggesting that they may target different mRNA transcripts. Longer small RNAs (26–31 nt), which were abundant in the gonads, may be putative piRNAs because of their length and their origin from repetitive elements. Additionally, our data include possible novel classes of small RNAs.

Conclusions

We elucidated miRNA expression patterns in various organs of T. rubripes. Most miRNA sequences are conserved in vertebrates, indicating that the basic functions of vertebrate miRNAs share a common evolution. Some miRNA species exhibit different distributions of isomiRs between tissues, suggesting that they have a broad range of functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1622-1) contains supplementary material, which is available to authorized users.

Dramatic improvement in genome assembly achieved using doubled-haploid genomes

Improvement in de novo assembly of large genomes is still to be desired. Here, we improved draft genome sequence quality by employing doubled-haploid individuals. We sequenced wildtype and doubled-haploid Takifugu rubripes genomes, under the same conditions, using the Illumina platform and assembled contigs with SOAPdenovo2. We observed 5.4-fold and 2.6-fold improvement in the sizes of the N50 contig and scaffold of doubled-haploid individuals, respectively, compared to the wildtype, indicating that the use of a doubled-haploid genome aids in accurate genome analysis.

Characterization and antimicrobial susceptibility of motile aeromonads isolated from freshwater ornamental fish showing signs of septicaemia

A total of 74 phenotypically identified presumptive motile Aeromonas isolates recovered from septicaemic freshwater ornamental fish in Sri Lanka were genetically characterized by sequencing of rpoD and gyrB genes. rpoD/gyrB phylogeny confirmed only 53 isolates as Aeromonas, among which A. veronii was the predominant species (79.2%), followed by A. hydrophila (7.5%), A. caviae (5.7%), A. jandaei (1.9%), A. dhakensis (3.8%) and A. entero pelogenes (1.9%). The aeromonads confirmed by sequencing were further subjected to 16S rDNA PCR-RFLP which substantiated sequencing results for 83% of isolates. Fingerprinting of A. enteropelogenes (n = 42) using ERIC-PCR revealed no dominant clones, and the majority were genetically distinct. All isolates were screened by PCR for 7 virulence determinant genes (aer, act, ast, alt, fla, ser, exu) and 2 integrase encoding genes (intI1, intI2). Each isolate contained ≥3 of the virulence genes tested for, with a heterogeneous distribution. Of the isolates, 77% harboured the intI1 gene, while none had intI2. In vitro antimicrobial susceptibility testing showed highest resistances towards tetracycline (58.5%) and erythromycin (54.7%). Our results indicate the diverse range of aeromonads that could potentially be associated with motile aeromonad septicaemia in ornamental fish. This is the first isolation of A. dhakensis from a septicaemic ornamental fish since its original description from the same host.

DPP and DSP are Necessary for Maintaining TGF-β1 Activity in Dentin

Porcine dentin sialophosphoprotein (DSPP) is the most abundant non-collagenous protein in dentin. It is processed by proteases into 3 independent proteins: dentin sialoprotein (DSP), dentin glycoprotein (DGP), and dentin phosphoprotein (DPP). We fractionated DPP and DSP along with TGF-β activity by ion exchange (IE) chromatography from developing pig molars and measured their alkaline phosphatase (ALP)-stimulating activity in human periodontal (HPDL) cells with or without TGF-β receptor inhibitor. We then purified TGF-β-unbound or -bound DPP and DSP by reverse-phase high-performance liquid chromatography (RP-HPLC) using the ALP-HPDL system. The TGF-β isoform bound to DPP and DSP was identified as being TGF-β1 by both ELISA and LC-MS/MS analysis. We incubated carrier-free human recombinant TGF-β1 (CF-hTGF-β1) with TGF-β-unbound DPP or DSP and characterized the binding on IE-HPLC using the ALP-HPDL system. When only CF-hTGF-β1 was incubated, approximately 3.6% of the ALP-stimulating activity remained. DPP and DSP rescued the loss of TGF-β1 activity. Approximately 19% and 10% of the ALP stimulating activities were retained by the binding of TGF-β to DPP and DSP, respectively. The type I collagen infrequently bound to CF-hTGF-β1. We conclude that both DPP and DSP help retain TGF-β1 activity in porcine dentin.

Novel genes participating in the formation of prismatic and nacreous layers in the pearl oyster as revealed by their tissue distribution and RNA interference knockdown

In our previous publication, we identified novel gene candidates involved in shell formation by EST analyses of the nacreous and prismatic layer-forming tissues in the pearl oyster Pinctada fucata. In the present study, 14 of those genes, including two known genes, were selected and further examined for their involvement in shell formation using the RNA interference. Molecular characterization based on the deduced amino acid sequences showed that seven of the novel genes encode secretory proteins. The tissue distribution of the transcripts of the genes, as analyzed by RT-PCR and in situ hybridization, was mostly consistent with those obtained by the EST analysis reported previously. Shells in the pearl oysters injected with dsRNAs targeting genes 000027, 000058, 000081, 000096, 000113 (nacrein), 000118, 000133 and 000411 (MSI60), which showed expression specific to the nacreous layer forming tissues, showed abnormal surface appearance in this layer. Individuals injected with dsRNAs targeting genes 000027, 000113 and 000133 also exhibited abnormal prismatic layers. Individuals injected with dsRNAs targeting genes 000031, 000066, 000098, 000145, 000194 and 000200, which showed expression specific to prismatic layer forming tissues, displayed an abnormal surface appearance in both the nacreous and prismatic layers. Taken together, the results suggest that the genes involved in prismatic layer formation might also be involved in the formation of the nacreous layers.