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.

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.

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.

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.

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.

Structural and functional analyses of a TIMP and MMP in the ligament of Pinctada fucata.

The bivalve hinge ligament is the hard tissue that functions to open and close shells. The ligament contains fibrous structures consisting of aragonite crystals surrounded by a dense organic matrix. This organic matrix may contribute to the formation of fibrous aragonite crystals, but the mechanism underlying this formation remains unclear. In this study, we identified a novel ligament-specific protein, Pinctada fucata tissue inhibitor of metalloproteinase (PfTIMP), from the fibrous organic matrix between aragonite crystals in the ligament using the amino acid sequence and cDNA cloning methods. PfTIMP consists of 143 amino acid residues and has a molecular weight of 13,580.4. To investigate the activity of PfTIMP, inhibition of matrix metalloproteinase (MMP) activity was measured. PfTIMP strongly inhibited human MMP13 and MMP9. Eight MMP homologs were identified from a P. fucata genomic database by BLAST search. To identify the specific MMP that may contribute to ligament formation, the expression level of each MMP was measured in the mantle isthmus, which secretes the ligament. The expression of MMP54089 increased after scratching of the ligament, while the expressions of other MMPs did not increase after doing the same operation. To identify the role of MMP54089 in forming the ligament structure, double stranded (ds) RNA targeting MMP54089 was injected into living P. fucata to suppress the function of MMP54089. Scanning electron microscopic images showed disordered growing surfaces of the ligament in individuals injected with MMP54089-specific dsRNA. These results suggest that PfTIMP and MMP54089 play important roles in the formation of the fibrous ligament structure.

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.

Stimulatory and inhibitory mechanisms of slow muscle-specific myosin heavy chain gene expression in fish: transient and transgenic analysis of torafugu MYH(M86-2) promoter in zebrafish embryos.

The myosin heavy chain gene, MYHM86-2, exhibited restricted expression in slow muscle fibers of torafugu embryos and larvae, suggesting its functional roles for embryonic and larval muscle development. However, the transcriptional mechanisms involved in its expression are still ambiguous. The present study is the first extensive analysis of slow muscle-specific MYHM86-2 promoter in fish for identifying the cis-elements that are crucial for its expression. Combining both transient transfection and transgenic approaches, we demonstrated that the 2614bp 5'-flanking sequences of MYHM86-2 contain a sufficient promoter activity to drive gene expression specific to superficial slow muscle fibers. By cyclopamine treatment, we also demonstrated that the differentiation of such superficial slow muscle fibers depends on hedgehog signaling activity. The deletion analyses defined an upstream fragment necessary for repressing ectopic MYHM86-2 expression in the fast muscle fibers. The transcriptional mechanism that prevents MYHM86-2 expression in the fast muscle fibers is mediated through Sox6 binding elements. We also demonstrated that Sox6 may function as a transcriptional repressor of MYHM86-2 expression. We further discovered that nuclear factor of activated T cells (NFAT) binding elements plays a key role and myocyte enhancer factor-2 (MEF2) binding elements participate in the transcriptional regulation of MYHM86-2 expression.

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.

Tissue Localization of Tetrodotoxin in the Flatworm Planocera multitentaculata (Platyhelminthes: Polycladida).

Tetrodotoxin (TTX), a pufferfish toxin, is a highly potent neurotoxin that has been found in a wide variety of animals. The TTX-bearing flatworm Planocera multitentaculata possesses a large amount of TTX and is considered responsible for the toxification of TTX-bearing animals such as pufferfish (Takifugu and Chelonodon) and the toxic goby Yongeichthys criniger. However, the mechanism underlying TTX accumulation in flatworms remains unclear. Previous studies have been limited to identifying the distribution of TTX in multiple organs, such as the digestive organs, genital parts, and the remaining tissues of flatworms. Here, we performed liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and immunohistochemical staining using a monoclonal anti-TTX antibody to elucidate the detailed localization of TTX in the tissues and organs of the flatworm P. multitentaculata. Immunohistochemical staining for P. multitentaculata showed that TTX-specific signals were detected not only in the ovaries and pharynx but also in many other tissues and organs, whereas no signal was detected in the brain, Lang's vesicle, and genitalia. In addition, combined with LC-MS/MS analysis, it was revealed for the first time that TTX accumulates in high concentrations in the basement membrane and epidermis. These findings robustly support the hypotheses of "TTX utilization protection from predators."

Japanese Planocerid Flatworms: Difference in Composition of Tetrodotoxin and Its Analogs and the Effects of Ingestion by Toxin-Bearing Fishes in the Ryukyu Islands, Japan.

Tetrodotoxin (TTX), known as pufferfish toxin, is a potent neurotoxin blocking sodium channels in muscle and nerve tissues. TTX has been detected in various taxa other than pufferfish, including marine polyclad flatworms, suggesting that pufferfish toxin accumulates in fish bodies via food webs. The composition of TTX and its analogs in the flatworm Planocera multitentaculata was identical to those in wild grass puffer Takifugu alboplumbeus. Previously, Planocera sp. from Okinawa Island, Japan, were reported to possess high level of TTX, but no information was available on TTX analogs in this species. Here we identified TTX and analogs in the planocerid flatworm using high-resolution liquid chromatography-mass spectrometry, and compared the composition of TTX and analogs with those of another toxic and non-toxic planocerid species. We show that the composition of TTX and several analogs, such as 5,6,11-trideoxyTTX, dideoxyTTXs, deoxyTTXs, and 11-norTTX-6(S)-ol, of Planocera sp. was identical to those of toxic species, but not to its non-toxic counterpart. The difference in the toxin composition was reflected in the phylogenetic relationship based on the mitochondrial genome sequence. A toxification experiment using predatory fish and egg plates of P. multitentaculata demonstrated that the composition of TTX and analogs in wild T. alboplumbeus juveniles was reproduced in artificially toxified pufferfish. Additionally, feeding on the flatworm egg plates enhanced the signal intensities of all TTX compounds in Chelonodon patoca and that of deoxyTTXs in Yongeichthys criniger.

Evolution of the myosin heavy chain gene MYH14 and its intronic microRNA miR-499: muscle-specific miR-499 expression persists in the absence of the ancestral host gene.

BACKGROUND: A novel sarcomeric myosin heavy chain gene, MYH14, was identified following the completion of the human genome project. MYH14 contains an intronic microRNA, miR-499, which is expressed in a slow/cardiac muscle specific manner along with its host gene; it plays a key role in muscle fiber-type specification in mammals. Interestingly, teleost fish genomes contain multiple MYH14 and miR-499 paralogs. However, the evolutionary history of MYH14 and miR-499 has not been studied in detail. In the present study, we identified MYH14/miR-499 loci on various teleost fish genomes and examined their evolutionary history by sequence and expression analyses. RESULTS: Synteny and phylogenetic analyses depict the evolutionary history of MYH14/miR-499 loci where teleost specific duplication and several subsequent rounds of species-specific gene loss events took place. Interestingly, miR-499 was not located in the MYH14 introns of certain teleost fish. An MYH14 paralog, lacking miR-499, exhibited an accelerated rate of evolution compared with those containing miR-499, suggesting a putative functional relationship between MYH14 and miR-499. In medaka, Oryzias latipes, miR-499 is present where MYH14 is completely absent in the genome. Furthermore, by using in situ hybridization and small RNA sequencing, miR-499 was expressed in the notochord at the medaka embryonic stage and slow/cardiac muscle at the larval and adult stages. Comparing the flanking sequences of MYH14/miR-499 loci between torafugu Takifugu rubripes, zebrafish Danio rerio, and medaka revealed some highly conserved regions, suggesting that cis-regulatory elements have been functionally conserved in medaka miR-499 despite the loss of its host gene. CONCLUSIONS: This study reveals the evolutionary history of the MYH14/miRNA-499 locus in teleost fish, indicating divergent distribution and expression of MYH14 and miR-499 genes in different teleost fish lineages. We also found that medaka miR-499 was even expressed in the absence of its host gene. To our knowledge, this is the first report that shows the conversion of intronic into non-intronic miRNA during the evolution of a teleost fish lineage.

Initiating the mollusk genomics annotation community: toward creating the complete curated gene-set of the Japanese Pearl Oyster, Pinctada fucata.

The genome sequence of the Japanese pearl oyster, the first draft genome from a mollusk, was published in February 2012. In order to curate the draft genome assemblies and annotate the predicted gene models, two annotation Jamborees were held in Okinawa and Tokyo. To date, 761 genes have been surveyed and curated. A preparatory meeting and a debriefing were held at the Misaki Marine Biological Station before and after the Jamborees. These four events, in conjunction with the sequence-decoding project, have facilitated the first series of gene annotations. Genome annotators among the Jamboree participants added 22 functional categories to the annotation system to date. Of these, 17 are included in Generic Gene Ontology. The other five categories are specific to molluskan biology, such as "Byssus Formation" and "Shell Formation", including Biomineralization and Acidic Proteins. A total of 731 genes from our latest version of gene models are annotated and classified into these 22 categories. The resulting data will serve as a useful reference for future genomic analyses of this species as well as comparative analyses among mollusks.

The diversity of shell matrix proteins: genome-wide investigation of the pearl oyster, Pinctada fucata.

In molluscs, shell matrix proteins are associated with biomineralization, a biologically controlled process that involves nucleation and growth of calcium carbonate crystals. Identification and characterization of shell matrix proteins are important for better understanding of the adaptive radiation of a large variety of molluscs. We searched the draft genome sequence of the pearl oyster Pinctada fucata and annotated 30 different kinds of shell matrix proteins. Of these, we could identified Perlucin, ependymin-related protein and SPARC as common genes shared by bivalves and gastropods; however, most gastropod shell matrix proteins were not found in the P. fucata genome. Glycinerich proteins were conserved in the genus Pinctada. Another important finding with regard to these annotated genes was that numerous shell matrix proteins are encoded by more than one gene; e.g., three ACCBP-like proteins, three CaLPs, five chitin synthase-like proteins, two N16 proteins (pearlins), 10 N19 proteins, two nacreins, four Pifs, nine shematrins, two prismalin-14 proteins, and 21 tyrosinases. This diversity of shell matrix proteins may be implicated in the morphological diversity of mollusc shells. The annotated genes reported here can be searched in P. fucata gene models version 1.1 and genome assembly version 1.0 ( http://marinegenomics.oist.jp/pinctada_fucata ). These genes should provide a useful resource for studies of the genetic basis of biomineralization and evaluation of the role of shell matrix proteins as an evolutionary toolkit among the molluscs.

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.

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.

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.

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.

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.