Search

The Antarctic toothfish, Dissostichus mawsoni, belongs to the Nototheniidae family and is distributed in sub-zero temperatures below S60° latitude in the Southern Ocean. Therefore, it is an attractive model species to study the stenothermal cold-adapted character state. In this study, we successfully generated highly contiguous genome sequences of D. mawsoni, which contained 1 062 scaffolds with a N50 length of 36.98 Mb and longest scaffold length of 46.82 Mb. Repetitive elements accounted for 40.87% of the genome. We also inferred 32 914 protein-coding genes using in silico gene prediction and transcriptome sequencing and detected splicing variants using Isoform-Sequencing (Iso-Seq), which will be invaluable resource for further exploration of the adaptation mechanisms of Antarctic toothfish. This new high-quality reference genome of D. mawsoni provides a fundamental resource for a deeper understanding of cold adaptation and conservation of species.

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome-level reference genome for T. obscurus was constructed using long-read Pacific Biosciences (PacBio) Sequel sequencing and a Hi-C-based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single-molecule real-time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome-scale scaffolds using a Hi-C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one-to-one orthologous relationship was observed between the two species, supporting the chromosome-level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu-specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.

The genus Pogonophryne is a speciose group that includes 28 species inhabiting the coastal or deep waters of the Antarctic Southern Ocean. The genus has been divided into five species groups, among which the P. albipinna group is the most deep-living group and is characterized by a lack of spots on the top of the head. Here, we carried out genome survey sequencing of P. albipinna using the Illumina HiSeq platform to estimate the genomic characteristics and identify genome-wide microsatellite motifs. The genome size was predicted to be ∼883.8 Mb by K-mer analysis (K = 25), and the heterozygosity and repeat ratio were 0.289 and 39.03%, respectively. The genome sequences were assembled into 571624 contigs, covering a total length of ∼819.3 Mb with an N50 of 2867 bp. A total of 2217422 simple sequence repeat (SSR) motifs were identified from the assembly data, and the number of repeats decreased as the length and number of repeats increased. These data will provide a useful foundation for the development of new molecular markers for the P. albipinna group as well as for further whole-genome sequencing of P. albipinna.

The southern elephant seal Mirounga leonina is the largest member of the phocid seals and is a highly sexually dimorphic predator at circumpolar regions. In this study, we generated the first high-quality de novo assembly of the blood transcriptome of M. leonina using the Illumina MiSeq platform. A total of 40.5 million raw reads were obtained and assembled using the Oases assembly platform, resulting in 46,445 contigs with an N50 value of 634 bp. We performed functional gene annotations through pathway analyses of the Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Ortholog comparison revealed that a high proportion of the M. leonina blood transcriptome had significant sequence homology within pinnipeds. M. leonine, a deep and long-diving seal, is routinely exposed to progressive hypoxia during dives. KEGG analyses detected an intact hypoxia-inducible factor-1 (HIF-1) signaling pathway, which is a key metabolism pathway in the adaptive homeostatic responses to hypoxia. In the blood transcriptome, many essential genes involved in innate immunity were detected, which suggests that the blood serves as a host defense modulator against exogenous infections in M. leonina. This genomic resource will be useful for understanding adaptive metabolism upon repeated breath-hold dives and determining the health status of southern elephant seals without the need to sacrifice them in experiments.

The southern elephant seal Mirounga leonina is the largest phocid seal and one of the two species of elephant seals. They are listed as &lsquo, least concern&rsquo, by the International Union for Conservation of Nature (IUCN) Red List of Threatened Species 2015. Here, we have assembled the reference genome for M. leonina using the 10&times, chromium sequencing platform. The final genome assembly of M. leonina was 2.42 Gb long, with a contig N50 length of 54 Mb and a maximum length of 111.6 Mb. The M. leonina genome contained 20,457 predicted protein-coding genes and possessed 41.51% repeated sequences. The completeness of the M. leonina genome was evaluated using benchmarking universal single-copy orthologous genes (BUSCOs): the assembly was highly complete, containing 95.6% of the core set of mammalian genes. The high-quality genomic information on M. leonina will be essential for further understanding of adaptive metabolism upon repeated breath-hold dives and the exploration of molecular mechanisms contributing to its unique biochemical and physiological characteristics. The southern elephant seal genome project was deposited at NCBI (National Center for Biotechnology Information) under BioProject number PRJNA587380.

The complete mitochondrial genome of Trematomus loennbergii was studied using NGS technology with PacBio platform. The mitochondrial genome size was 19,374bp and it had 13 protein-coding genes, 22 tRNAs and 2 rRNAs. There were 4 types of stop codons which were TAA, TAG, AGG and T(AA) but start codon type was only one (ATG). The contents of GC were 44.09% and AT contents were 55.91%. To conduct phylogenetic analysis, 12 species in 3 families were used. The result suggested that T. loennbergii was close to Pagothenia borchgrevinki in Nototheniidae. This study would provide a fundamental data for molecular evolution of T. loennbergii.

The complete mitochondrial genome of Eaton’s skate Bathyraja eatonii was studied using the long-read technology, PacBio Sequel System. The complete mitochondrial genome form of B. eatonii was 16,69...

The brine shrimp Artemia has a ZW sex determination system with ZW chromosomes in females and ZZ chromosomes in males. Artemia has been considered a promising model organism for ZW sex-determining systems, but the genes involved in sex determination and differentiation of Artemia have not yet been identified. Here, we conducted transcriptome sequencing of female and male A. franciscana using PacBio Iso-Seq and Illumina RNA-Seq techniques to identify candidate sex determination genes. Among the 42,566 transcripts obtained from Iso-Seq, 23,514 were analyzed. Of these, 2065 (8.8%) were female specific, 2513 (10.7%) were male specific, and 18,936 (80.5%) were co-expressed in females and males. Based on GO enrichment analysis and expression values, we found 10 female-biased and 29 male-biased expressed genes, including DMRT1 and Sad genes showing male-biased expression. Our results showed that DMRT1 has three isoforms with five exons, while Sad has seven isoforms with 2–11 exons. The Sad gene is involved in ecdysteroid signaling related to molting and metamorphosis in arthropods. Further studies on ecdysteroid biosynthetic genes are needed to improve our understanding of Artemia sex determination. This study will provide a valuable resource for sex determination and differentiation studies on Artemia and other crustaceans with ZW systems.

The full-length mitochondrial genome of Muraenolepis orangiensis (Vaillant, 1888) was studied using PacBio platform and it is first report in a Muraenolepididae family. The circular form of mitochondria genome is 16,833 bp including 13 protein-coding genes, two rRNA, and 22 tRNA. Start codon of 13 protein-coding genes was only ATG but three types of stop codons (TAA, T(AA), and TAG) were detected. To evaluate evolutionary position of M. orangiensis, the phylogenetic tree with other 13 Antarctic fishes belonged to five families were showed that M. orangiensis is unique cluster as a Muraenolepididae family and this study would provide fundamental data to understand the evolutionary relationship of fishes founded in Antarctic area.

The complete mitochondrial genome of Notothenia rossii was obtained using PacBio Sequel long-read sequencing platform. The mitogenome of N. rossii was circular form and 18,274 bp long, which consists of 13 protein-coding genes, 24 tRNAs, 2 rRNAs, and non-coding control region. Particularly, we found duplicated tRNAThr and tRNAPro in addition to the typical 22 tRNAs. The phylogenetic tree revealed that N. rossii was most closely related to N. coriiceps among species in the Nototheniidae clade within the suborder Notothenioidei.

The complete mitochondrial genome of Macrourus witsoni was determined in this study by the Long-read Technology, such as PacBio Sequel System. The Long-read Technology, which can sequence continuously the whole vertebrate mitochondrial genome, allows more accurate genomes to be completed. The circular form of its mitochondrial genome was 16,714bp, which contained 13 protein-coding genes, 22 tRNA, and 2 rRNA. The gene orders of M.witsoni was identical to that of the other species of Macrouridae family. Phylogenetic analysis indicated M. witsoni was mostly close to C.kishinouyei in the Macrouridae family.

Ciliates are a diverse species group of the Protozoa, and nuclear and mitochondrial genes have been utilized to discover new species and discriminate closely related species. The mitochondrial cytochrome c oxidase subunit 1 (CO1) gene has been used to discriminate metazoan species and has also been applied for some groups in the phylum Ciliophora. However, it is difficult to produce a universal primer as a standard barcode, because unlike metazoans, mitochondrial DNA sequences of ciliates are long and highly variable. Therefore, to design the new primer set, we sequenced the mitochondrial genomes of two pseudokeronopsids in the class Spirotrichea using next-generation sequencing technology (HiSeq™ 2000). Based on putative CO1 gene fragments of the pseudokeronopsids, we designed the new primer set and successfully sequenced the CO1 of 69 populations representing 47 species (five orders, 14 families, and 27 genera). We found that CO1 showed higher resolution for separating congeneric species than did nuclear SSU rRNA gene sequences, and we identified some putative cryptic species.

The chinstrap (Pygoscelis antarcticus) and gentoo (P. papua) penguins are distributed throughout Antarctica and the sub-Antarctic islands. In this study, high-quality de novo assemblies of blood transcriptomes from these penguins were generated using the Illumina MiSeq platform. A total of 22.2 and 21.8 raw reads were obtained from chinstrap and gentoo penguins, respectively. These reads were assembled using the Oases assembly platform and resulted in 26,036 and 21,854 contigs with N50 values of 929 and 933 base pairs, respectively. Functional gene annotations through pathway analyses of the Gene Ontology, EuKaryotic Orthologous Groups, and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were performed for each blood transcriptome, resulting in a similar compositional order between the two transcriptomes. Ortholog comparisons with previously published transcriptomes from the Adélie (P. adeliae) and emperor (Aptenodytes forsteri) penguins revealed that a high proportion of the four penguins’ transcriptomes had significant sequence homology. Because blood and tissues of penguins have been used to monitor pollution in Antarctica, immune parameters in blood could be important indicators for understanding the health status of penguins and other Antarctic animals. In the blood transcriptomes, KEGG analyses detected many essential genes involved in the major innate immunity pathways, which are key metabolic pathways for maintaining homeostasis against exogenous infections or toxins. Blood transcriptome studies such as this may be useful for checking the immune and health status of penguins without sacrifice.

In this study, we determined the complete chloroplast (cp) genomes from two varieties of Cannabis sativa. The genome sizes were 153,848 bp (the Korean non-drug variety, Cheungsam) and 153,854 bp (the African variety, Yoruba Nigeria). The genome structures were identical with 131 individual genes [86 protein-coding genes (PCGs), eight rRNA, and 37 tRNA genes]. Further, except for the presence of an intron in the rps3 genes of two C. sativa varieties, the cp genomes of C. sativa had conservative features similar to that of all known species in the order Rosales. To verify the position of C. sativa within the order Rosales, we conducted phylogenetic analysis by using concatenated sequences of all PCGs from 17 complete cp genomes. The resulting tree strongly supported monophyly of Rosales. Further, the family Cannabaceae, represented by C. sativa, showed close relationship with the family Moraceae. The phylogenetic relationship outlined in our study is well congruent with those previously shown for the order Rosales.