Long-read sequencing and de novo genome assembly of marine medaka (Oryzias melastigma)

Genome-wide identification of 769 G protein-coupled receptor (GPCR) genes from the marine medaka Oryzias melastigma

The marine medaka Oryzias melastigma is a useful fish model for marine and estuarine ecotoxicology studies and can be applied to field-based population genomics because of its distribution in Asian estuaries and other coastal areas. We identified 769 full-length G protein-coupled receptor genes in the O. melastigma genome and classified them into five distinct classes. A phylogenetic comparison of GPCR genes in O. melastigma to humans and two other small fish species revealed a high-level orthological relationship. Purinergic and chemokine receptors were highly differentiated in humans whereas significant differentiation of chemosensory receptors was evident in fish species. Our results suggest that the GPCR gene families among the species used in this study exhibit evidence of sporadic evolutionary processes. These results may help improve our understanding of the advanced repertoires of GPCR and expand our knowledge of physiological mechanisms of fish in response to various environmental stimuli.

Methyltransferase-like 3 suppresses red spotted grouper nervous necrosis virus and viral hemorrhagic septicemia virus infection by enhancing type I interferon responses in sea perch

Methylation at the N6 position of adenosine (m6A) is the most abundant internal mRNA modification in eukaryotes, tightly associating with regulation of viral life circles and immune responses. Here, a methyltransferase-like 3 homolog gene from sea perch (Lateolabrax japonicus), designated LjMETTL3, was cloned and characterized, and its negative role in fish virus pathogenesis was uncovered. The cDNA of LjMETTL3 encoded a 601-amino acid protein with a MT-A70 domain, which shared the closest genetic relationship with Echeneis naucrates METTL3. Spatial expression analysis revealed that LjMETTL3 was more abundant in the immune tissues of sea perch post red spotted grouper nervous necrosis virus (RGNNV) or viral hemorrhagic septicemia virus (VHSV) infection. LjMETTL3 expression was significantly upregulated at 12 and 24 h post RGNNV and VHSV infection in vitro. In addition, ectopic expression of LjMETTL3 inhibited RGNNV and VHSV infection in LJB cells at 12 and 24 h post infection, whereas knockdown of LjMETTL3 led to opposite effects. Furthermore, we found that LjMETTL3 may participate in boosting the type I interferon responses by interacting with TANK-binding kinase. Taken together, these results disclosed the antiviral role of fish METTL3 against RGNNV and VHSV and provided evidence for understanding the potential mechanisms of fish METTL3 in antiviral innate immunity.

Hypoxia-induced epigenetic transgenerational miRNAs dysregulation involved in reproductive impairment of ovary

Hypoxia is a potent endocrine disruptor that is posing serious problems to the fish reproductive systems. Our previous studies reported that hypoxia could cause a transgenerational impairment of ovarian development and interfere hatching success in F2 offspring of marine medaka fish (Oryzias melastigma) through epigenetic regulation. As part of the epigenetic regulation, we investigated the involvement of microRNAs (miRNAs) in hypoxia-induced transgenerational reproductive impairments. In the present study, we used comparative small RNA sequencing to reveal that hypoxia caused miRNA dysregulation in ovaries of F0 hypoxia group and F2 transgenerational group. We found 4 common dysregulated miRNA in the F0 and F2 generations. Furthermore, integrated miRNA-mRNA analysis, followed by gene ontology enrichment analysis on the hypoxia-dysregulated miRNA-target genes further highlighted the importance of these dysregulated miRNAs in biological processes related to reproduction. More importantly, we identified 3 miRNA-mRNA pairs (novel miRNA-525-DIAPH2, novel miRNA-525-MYOCD, and novel miRNA-525-RAI14) that might play epigenetic roles in hypoxia-induced reproductive impairment. For the first time, our findings suggested the involvement of miRNA in hypoxia-induced reproductive impairments may be inherited via a transgenerational manner.

Genomic variations and epigenomic landscape of the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel

BACKGROUND

The teleost medaka (Oryzias latipes) is a well-established vertebrate model system, with a long history of genetic research, and multiple high-quality reference genomes available for several inbred strains. Medaka has a high tolerance to inbreeding from the wild, thus allowing one to establish inbred lines from wild founder individuals.

RESULTS

We exploit this feature to create an inbred panel resource: the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel. This panel of 80 near-isogenic inbred lines contains a large amount of genetic variation inherited from the original wild population. We use Oxford Nanopore Technologies (ONT) long read data to further investigate the genomic and epigenomic landscapes of a subset of the MIKK panel. Nanopore sequencing allows us to identify a large variety of high-quality structural variants, and we present results and methods using a pan-genome graph representation of 12 individual medaka lines. This graph-based reference MIKK panel genome reveals novel differences between the MIKK panel lines and standard linear reference genomes. We find additional MIKK panel-specific genomic content that would be missing from linear reference alignment approaches. We are also able to identify and quantify the presence of repeat elements in each of the lines. Finally, we investigate line-specific CpG methylation and performed differential DNA methylation analysis across these 12 lines.

CONCLUSIONS

We present a detailed analysis of the MIKK panel genomes using long and short read sequence technologies, creating a MIKK panel-specific pan genome reference dataset allowing for investigation of novel variation types that would be elusive using standard approaches.

RNA-seq analyses of Marine Medaka (Oryzias melastigma) reveals salinity responsive transcriptomes in the gills and livers

Increasing salinity levels in marine and estuarine ecosystems greatly influence developmental, physiological and molecular activities of inhabiting fauna. Marine medaka (Oryzias melastigma), a euryhaline research model, has extraordinary abilities to survive in a wide range of aquatic salinity. To elucidate how marine medaka copes with salinity differences, the responses of Oryzias melastigma after being transferred to different salt concentrations [0 practical salinity units (psu), 15 psu, 30 psu (control), 45 psu] were studied at developmental, histochemical and transcriptome levels in the gill and liver tissues. A greater number of gills differentially expressed genes (DEG) under 0 psu (609) than 15 psu (157) and 45 psu (312), indicating transcriptomic adjustments in gills were more sensitive to the extreme hypotonic environment. A greater number of livers DEGs were observed in 45 psu (1,664) than 0 psu (87) and L15 psu (512), suggesting that liver was more susceptible to hypertonic environment. Further functional analyses of DEGs showed that gills have a more immediate response, mainly in adjusting ion balance, immune and signal transduction. In contrast, DEGs in livers were involved in protein synthesis and processing. We also identified common DEGs in both gill and liver and found they were mostly involved in osmotic regulation of amino sugar and nucleotide sugar metabolism and steroid biosynthesis. Additionally, salinity stresses showed no significant effects on most developmental and histochemical parameters except increased heartbeat with increasing salinity and decreased glycogen after transferred from stable conditions (30 psu) to other salinity environments. These findings suggested that salinity-stress induced changes in gene expressions could reduce the effects on developmental and histochemical parameters. Overall, this study provides a useful resource for understanding the molecular mechanisms of fish responses to salinity stresses.

Genome Assembly of Salicaceae Populus deltoides (Eastern Cottonwood) I-69 Based on Nanopore Sequencing and Hi-C Technologies

Populus deltoides has important ecological and economic values, widely used in poplar breeding programs due to its superior characteristics such as rapid growth and resistance to disease. Although the genome sequence of P. deltoides WV94 is available, the assembly is fragmented. Here, we reported an improved chromosome-level assembly of the P. deltoides cultivar I-69 by combining Nanopore sequencing and chromosome conformation capture (Hi-C) technologies. The assembly was 429.3 Mb in size and contained 657 contigs with a contig N50 length of 2.62 Mb. Hi-C scaffolding of the contigs generated 19 chromosome-level sequences, which covered 97.4% (418 Mb) of the total assembly size. Moreover, repetitive sequences annotation showed that 39.28% of the P. deltoides genome was composed of interspersed elements, including retroelements (23.66%), DNA transposons (6.83%), and unclassified elements (8.79%). We also identified a total of 44 362 protein-coding genes in the current P. deltoides assembly. Compared with the previous genome assembly of P. deltoides WV94, the current assembly had some significantly improved qualities: the contig N50 increased 3.5-fold and the proportion of gaps decreased from 3.2% to 0.08%. This high-quality, well-annotated genome assembly provides a reliable genomic resource for identifying genome variants among individuals, mining candidate genes that control growth and wood quality traits, and facilitating further application of genomics-assisted breeding in populations related to P. deltoides.