Comparative genomic analysis of copepod humoral immunity genes with sex-biased expression in Labidocera rotunda.

Studies of innate immune system function in invertebrates have contributed significantly to our understanding of the mammalian innate immune system. However, in-depth research on innate immunity in marine invertebrates remains sparse. We generated the first de novo genome and transcriptome sequences of copepod Labidocera rotunda using Illumina paired-end data and conducted a comparative genome analysis including five crustaceans (four copepods and one branchiopod species). We cataloged the presence of Toll, Imd, JAK/STAT, and JNK pathway components among them and compared them with 17 previously reported diverse arthropod species representative of insects, myriapods, chelicerates, and malacostracans. Our results indicated that copepod Gram-negative binding proteins may function in direct digestion or pathogen killing. The phylogenetic analysis of arthropod TEP and copepod-specific GCGEQ motif patterns suggested that the evolutionary history of copepod TEPs may have diverged from that of other arthropods. We classified the copepod Toll-like receptors identified in our analysis as either vertebrate or protostome types based on their cysteine motifs and the tree built with their Toll/interleukin-1 receptor domains. LrotCrustin, the first copepod AMP, was identified based on the structure of its WAP domain and deep-learning AMP predictors. Gene expression level analysis of L. rotunda innate immunity-related transcripts in each sex showed higher Toll pathway-related expression in male L. rotunda than in females, which may reflect an inverse correlation between allocation of reproductive investment and elevated immune response in males. Taken together, the results of our study provide insight into copepod innate immunity-related gene families and illuminate the evolutionary potential of copepods relative to other crustaceans.

Do red tide events promote an increase in zooplankton biodiversity?

Red tides occurring off the southern coast of Korea impact the marine ecosystem and aquaculture industries. Zooplankton are crucial in the food web, connecting primary producers to higher predators and interact diversely with red tide organisms. This study explores dynamics of the zooplankton community over seven years including three red tide and four non-red tide years in Tongyeong using metabarcoding. In non-red tide years, zooplankton diversity showed typical seasonal patterns, increasing from June to early October. However, during red tide years, diversity remained high, with a shift in species composition-decreased Copepoda and increased Branchiopoda, Echinodermata, Malacostraca, and Annelida. Diversity indices were significantly higher in red tide years across all periods except for the richness in "after" that showed an insignificant higher value. The differences in zooplankton assemblages across periods were influenced by surface temperatures and the density of the red tide-causing alga Margalefidinium polykrikoides. Eight species emerged as indicator species and showed direct correlations with M. polykrikoides and among them, seven species were indicator species for red tide occurrence years. The ecological characteristics of M. polykrikoides blooms and their recurrent occurrences over several decades suggest that zooplankton may adapt to the toxins and use these blooms as spawning cues. Overall, this study provides comprehensive understanding on changes in zooplankton communities during red tide events, offering novel insights into their ecology.

Whole genome sequencing and genome characterization of Aichivirus isolated from Korean adults.

The whole-genome sequence (WGS) analysis of Aichivirus (AiV) identified in Korea was performed in this study. Using Sanger and Nanopore sequencing, the 8228-nucleotide-long genomic sequence of AiV (OQ121963) was determined and confirmed to belong to genotype A. The full-length genome of OQ121963 consisted of a 7296 nt open reading frame (ORF) that encodes a single polyprotein, and 5' UTR (676 nt) and 3' UTR (256 nt) at 5' and 3' ends, respectively. The ORF consisted of leader protein (L), structural protein P1 (VP0, VP1, and VP3), and nonstructural protein P2 (2A, 2B, and 2C) and P3 (3A, 3B, 3C, and 3D). The secondary structure analysis of the 5' UTR identified only stem-loop C (SL-C) and not SL-A and SL-B. The variable region of the AiV genome was analyzed by MegAlign Pro and reconfirmed by SimPlot analysis using 16 AiV whole genomes known to date. Among the entire regions, structural protein region P1 showed the lowest amino acid identity (96.07%) with reference sequence AB040749 (originated in Japan; genotype A), while the highest amino acid identity (98.26%) was confirmed in the 3D region among nonstructural protein region P2 and P3. Moreover, phylogenetic analysis of the WGS of OQ121963 showed the highest homology (96.96%) with JX564249 (originated in Taiwan; genotype A) and lowest homology (90.14%) with DQ028632 (originated in Brazil; genotype B). Therefore, the complete genome characterization of OQ121963 and phylogenetic analysis of the AiV conducted in this study provide useful information allowing to improve diagnostic tools and epidemiological studies of AiVs.

Highly contiguous genome assembly and gene annotation of the short-finned eel (Anguilla bicolor pacifica).

In East Asia, anguillid eels are commercially important. However, unlike other species, they have not been successfully cultivated throughout their lifecycle. Facing population decline due to overharvesting and environmental pressures, the industry is turning to alternatives, such as Anguilla bicolor pacifica (short-finned eel). However, genomic data for short-finned eels are unavailable. Here, we present in-depth whole-genome sequencing results for short-finned eel obtained using two sequencing platforms (PacBio Revio, and Illumina). In this study, we achieved a highly contiguous genome assembly of the short-finned eel, comprising 19 pseudochromosomes encompassing 99.76% of the 1.087 Gb genome sequence with an N50 of 16.88 and 61.07 Mb from contig and scaffold, respectively. Transcripts from four different tissues led to the annotation of 23,095 protein-coding genes in the eel genome, 98.66% of which were functionally annotated. This high-quality genome assembly, along with the annotation data, provides a foundation for future functional genomic studies of short-finned eels.

The occurrence of positive selection on BicA transporter of Microcystis aeruginosa.

The rapid growth of cyanobacteria, particularly Microcystis aeruginosa, poses a significant threat to global water security. The proliferation of toxic Microcystis aeruginosa raises concerns due to its potential harm to human health and socioeconomic impacts. Dense blooms contribute to spatiotemporal inorganic carbon depletion, promoting interest in the roles of carbon-concentrating mechanisms (CCMs) for competitive carbon uptake. Despite the importance of HCO3- transporters, genetic evaluations and functional predictions in M. aeruginosa remain insufficient. In this study, we explored the diversity of HCO3- transporters in the genomes of 46 strains of M. aeruginosa, assessing positive selection for each. Intriguingly, although the Microcystis BicA transporter became a partial gene in 23 out of 46 genomic strains, we observed significant positive sites. Structural analyses, including predicted 2D and 3D models, confirmed the structural conservation of the Microcystis BicA transporter. Our findings suggest that the Microcystis BicA transport likely plays a crucial role in competitive carbon uptake, emphasizing its ecological significance. The ecological function of the Microcystis BicA transport in competitive growth during cyanobacterial blooms raises important questions. Future studies require experimental confirmation to better understand the role of the Microcysits BicA transporter in cyanobacterial blooms dynamics.

The m6A writer RBM15 drives the growth of triple-negative breast cancer cells through the stimulation of serine and glycine metabolism.

N6-adenosine methylation (m6A) is critical for controlling cancer cell growth and tumorigenesis. However, the function and detailed mechanism of how m6A methyltransferases modulate m6A levels on specific targets remain unknown. In the current study, we identified significantly elevated levels of RBM15, an m6A writer, in basal-like breast cancer (BC) patients compared to nonbasal-like BC patients and linked this increase to worse clinical outcomes. Gene expression profiling revealed correlations between RBM15 and serine and glycine metabolic genes, including PHGDH, PSAT1, PSPH, and SHMT2. RBM15 influences m6A levels and, specifically, the m6A levels of serine and glycine metabolic genes via direct binding to target RNA. The effects of RBM15 on cell growth were largely dependent on serine and glycine metabolism. Thus, RBM15 coordinates cancer cell growth through altered serine and glycine metabolism, suggesting that RBM15 is a new therapeutic target in BC.

Genomic analysis uncovers that cold-inducible RNA binding protein is associated with estrogen receptor in breast cancer.

BACKGROUND: RNA-binding proteins (RBPs) perform various biological functions in humans and are associated with several diseases, including cancer. Therefore, RBPs have emerged as novel therapeutic targets. Although recent investigations have shown that RBPs have crucial functions in breast cancer (BC), detailed research is underway to determine the RBPs that are closely related to cancers. OBJECTIVE: To provide an insight into estrogen receptor (ER) regulation by cold-inducible RNA binding protein (CIRBP) as a novel therapeutic target. RESULTS: By analyzing the genomic data, we identified a potential RBP in BC. We found that CIRBP is highly correlated with ER function and influences clinical outcomes, such as patient survival and endocrine therapy responsiveness. In addition, CIRBP influences the proliferation of BC cells by directly binding to ER-RNA. CONCLUSION: Our results suggest that CIRBP is a novel upstream regulator of ER and that the interplay between CIRBP and ER may be associated with the clinical relevance of BC.

Transcriptome and functional analyses of phenotypic plasticity in sea grape Caulerpa okamurae.

Caulerpa is a marine green macroalga distinguished by a large single cell with multiple nuclei. It also exhibits remarkable morphological intraspecies variations, in response to diverse environmental types. However, the molecular mechanisms underlying this phenotypic plasticity remain poorly understood. In this work, we compare the transcriptomes of Caulerpa okamurae Weber Bosse, 1897 displaying altered phenotypes of cultivation and natural phenotypes and investigate significantly regulated genes and their biological functions using differential expression analyses. We observe light-harvesting complex upregulation and cellular framework stability downregulation in altered phenotypes compared to the natural phenotypes. Intertidal macrophytes reduce light capture to avoid photodamage and regulate their morphology to protect against wave damage. In contrast, the lower light conditions and the cultivation environment augment light capture and increase a morphology prioritizing light trapping. Moreover, the addition of simulated wave-sweeping stimuli induces a return to the natural morphology under high-light conditions, showing how mechanical stress affects morphological organization in C. okamurae. We provide detailed gene expression patterns in C. okamurae under varying light intensities and water conditions, suggesting a distinct influence on its morphological traits.

BRD2-specific inhibitor, BBC0403, inhibits the progression of osteoarthritis pathogenesis in osteoarthritis-induced C57BL/6 male mice.

BACKGROUND AND PURPOSE: The discovery of new bromo- and extra-terminal inhibitors presents new drugs to treat osteoarthritis (OA). EXPERIMENTAL APPROACH: The new drug, BBC0403, was identified in the DNA-encoded library screening system by searching for compounds that target BRD (bromodomain-containing) proteins. The binding force with BRD proteins was evaluated using time-resolved fluorescence energy transfer (TR-FRET) and binding kinetics assays. Subsequently, in vitro and ex vivo analyses demonstrated the effects of the BRD2 inhibitor, BBC0403, on OA. For animal experiments, medial meniscus destabilization was performed to create a 12-week-old male C57BL/6 mouse model, and intra-articular (i.a.) injections were administered. Histological and immunohistochemical analyses were then performed. The underlying mechanism was confirmed by gene set enrichment analysis (GSEA) using RNA-seq. KEY RESULTS: TR-FRET and binding kinetics assays revealed that BBC0403 exhibited higher binding specificity for BRD2 compared to BRD3 and BRD4. The anti-OA effects of BBC0403 were tested at concentrations of 5, 10 and 20 µM (no cell toxicity in the range tested). The expression of catabolic factors, prostaglandin E2 (PGE2) production and extracellular matrix (ECM) degradation was reduced. Additionally, the i.a. injection of BBC0403 prevented OA cartilage degradation in mice. Finally, BBC0403 was demonstrated to suppress NF-κB and MAPK signalling pathways. CONCLUSION AND IMPLICATIONS: This study demonstrated that BBC0403 is a novel BRD2-specific inhibitor and a potential i.a.-injectable therapeutic agent to treat OA.

Dietary exposure of the water flea Daphnia galeata to microcystin-LR.

Harmful substances like the cyanotoxin microcystin-leucine-arginine (MC-LR) are commonly found in eutrophic freshwater environments, posing risks to aquatic organisms. The water flea, Daphnia, is a well-established model organism for environmental toxicology research. Nevertheless, there is currently insufficient research on the genes that respond to MC-LR in Daphnia galeata. This study aimed to gain insights into the notable genes that react significantly to MC-LR. In this study, we generated an extensive RNA-Seq sequences isolated from the D. galeata HK strain, Han River in Korea. This strain was nourished with a diet of the green microalga Chlorella vulgaris and treated with pure MC-LR at a concentration of 36 ug/L. The transcriptome profile in response to the MC-LR treatment was obtained and 336 differentially expressed genes were subjected to Gene Ontology (GO) and euKaryotic Orthologous Groups of proteins analyses. GO enrichment analysis showed that chemical stimulus, amino sugar metabolic and catabolic process, oxidative stress, and detoxification were highly enriched, in reverse, proteolysis and fucosylation were underpresented. Detoxification process related genes such as peroxidase-like, chorion, and thyroid peroxidase-like were enriched for eliminating or neutralizing MC_LR from an organism's body. Furthermore, functional protein classification revealed an upregulation of lipid and inorganic ion transport processes, while amino acid and carbohydrate transport processes were found to be downregulated. These findings offer insights into how organisms respond to ecotoxic stimuli, providing valuable information for understanding adaptation or defense pathways.

Comparative pangenome analysis of Aspergillus flavus and Aspergillus oryzae reveals their phylogenetic, genomic, and metabolic homogeneity.

Aspergillus flavus and Aspergillus oryzae are closely related fungal species with contrasting roles in food safety and fermentation. To comprehensively investigate their phylogenetic, genomic, and metabolic characteristics, we conducted an extensive comparative pangenome analysis using complete, dereplicated genome sets for both species. Phylogenetic analyses, employing both the entirety of the identified single-copy orthologous genes and six housekeeping genes commonly used for fungal classification, did not reveal clear differentiation between A. flavus and A. oryzae genomes. Upon analyzing the aflatoxin biosynthesis gene clusters within the genomes, we observed that non-aflatoxin-producing strains were dispersed throughout the phylogenetic tree, encompassing both A. flavus and A. oryzae strains. This suggests that aflatoxin production is not a distinguishing trait between the two species. Furthermore, A. oryzae and A. flavus strains displayed remarkably similar genomic attributes, including genome sizes, gene contents, and G + C contents, as well as metabolic features and pathways. The profiles of CAZyme genes and secondary metabolite biosynthesis gene clusters within the genomes of both species further highlight their similarity. Collectively, these findings challenge the conventional differentiation of A. flavus and A. oryzae as distinct species and highlight their phylogenetic, genomic, and metabolic homogeneity, potentially indicating that they may indeed belong to the same species.

AMPK-HIF-1α signaling enhances glucose-derived de novo serine biosynthesis to promote glioblastoma growth.

BACKGROUND: Cancer cells undergo cellular adaptation through metabolic reprogramming to sustain survival and rapid growth under various stress conditions. However, how brain tumors modulate their metabolic flexibility in the naturally serine/glycine (S/G)-deficient brain microenvironment remain unknown. METHODS: We used a range of primary/stem-like and established glioblastoma (GBM) cell models in vitro and in vivo. To identify the regulatory mechanisms of S/G deprivation-induced metabolic flexibility, we employed high-throughput RNA-sequencing, transcriptomic analysis, metabolic flux analysis, metabolites analysis, chromatin immunoprecipitation (ChIP), luciferase reporter, nuclear fractionation, cycloheximide-chase, and glucose consumption. The clinical significances were analyzed in the genomic database (GSE4290) and in human GBM specimens. RESULTS: The high-throughput RNA-sequencing and transcriptomic analysis demonstrate that the de novo serine synthesis pathway (SSP) and glycolysis are highly activated in GBM cells under S/G deprivation conditions. Mechanistically, S/G deprivation rapidly induces reactive oxygen species (ROS)-mediated AMP-activated protein kinase (AMPK) activation and AMPK-dependent hypoxia-inducible factor (HIF)-1α stabilization and transactivation. Activated HIF-1α in turn promotes the expression of SSP enzymes phosphoglycerate dehydrogenase (PHGDH), phosphoserine aminotransferase 1 (PSAT1), and phosphoserine phosphatase (PSPH). In addition, the HIF-1α-induced expression of glycolytic genes (GLUT1, GLUT3, HK2, and PFKFB2) promotes glucose uptake, glycolysis, and glycolytic flux to fuel SSP, leading to elevated de novo serine and glycine biosynthesis, NADPH/NADP+ ratio, and the proliferation and survival of GBM cells. Analyses of human GBM specimens reveal that the levels of overexpressed PHGDH, PSAT1, and PSPH are positively correlated with levels of AMPK T172 phosphorylation and HIF-1α expression and the poor prognosis of GBM patients. CONCLUSION: Our findings reveal that metabolic stress-enhanced glucose-derived de novo serine biosynthesis is a critical metabolic feature of GBM cells, and highlight the potential to target SSP for treating human GBM.

Integrated genomics and phenotype microarray analysis of Saccharomyces cerevisiae industrial strains for rice wine fermentation and recombinant protein production.

The industrial potential of Saccharomyces cerevisiae has extended beyond its traditional use in fermentation to various applications, including recombinant protein production. Herein, comparative genomics was performed with three industrial S. cerevisiae strains and revealed a heterozygous diploid genome for the 98-5 and KSD-YC strains (exploited for rice wine fermentation) and a haploid genome for strain Y2805 (used for recombinant protein production). Phylogenomic analysis indicated that Y2805 was closely associated with the reference strain S288C, whereas KSD-YC and 98-5 were grouped with Asian and European wine strains, respectively. Particularly, a single nucleotide polymorphism (SNP) in FDC1, involved in the biosynthesis of 4-vinylguaiacol (4-VG, a phenolic compound with a clove-like aroma), was found in KSD-YC, consistent with its lack of 4-VG production. Phenotype microarray (PM) analysis showed that KSD-YC and 98-5 displayed broader substrate utilization than S288C and Y2805. The SNPs detected by genome comparison were mapped to the genes responsible for the observed phenotypic differences. In addition, detailed information on the structural organization of Y2805 selection markers was validated by Sanger sequencing. Integrated genomics and PM analysis elucidated the evolutionary history and genetic diversity of industrial S. cerevisiae strains, providing a platform to improve fermentation processes and genetic manipulation.

Novel molecule BBC0901 inhibits BRD4 and acts as a catabolic regulator in the pathogenesis of osteoarthritis.

Osteoarthritis (OA) is induced by matrix degradation and inflammation mediated by bromo-domain-containing protein 4 (BRD4)-dependent catabolic factors. BRD4 acts as both a transcriptional regulator and an epigenetic reader. BBC0901 was identified as an inhibitor of BRD4 using a DNA-encoded library screening system. We aimed to demonstrate the effects of BBC0901 on OA pathogenesis by in vitro, ex vivo, and in vivo analyses. BBC0901 inhibited the expression of catabolic factors that degrade cartilage without significantly affecting the viability of mouse articular chondrocytes. Additionally, ex vivo experiments under conditions mimicking OA showed that BBC0901 suppressed extracellular matrix degradation. RNA sequencing analysis of gene expression patterns showed that BBC0901 inhibited the expression of catabolic factors, such as matrix metalloproteinases (MMPs) and cyclooxygenase (COX)2, along with reactive oxygen species (ROS) production. Furthermore, intra-articular (IA) injection of BBC0901 into the knee joint blocked osteoarthritic cartilage destruction by inhibition of MMP3, MMP13, COX2, interleukin (IL)6, and ROS production, thereby obstructing the nuclear factor kappa-light-chain-enhancer of activated B cell and mitogen activated protein kinase signaling. In conclusion, BBC0901-mediated BRD4 inhibition prevented OA development by attenuating catabolic signaling and hence, can be considered a promising IA therapeutic for OA.

Blockade of Activin Receptor IIB Protects Arthritis Pathogenesis by Non-Amplification of Activin A-ACVR2B-NOX4 Axis Pathway.

Although activin receptor IIB (ACVR2B) is emerging as a novel pathogenic receptor, its ligand and assembled components (or assembly) are totally unknown in the context of osteoarthritis (OA) pathogenesis. The present results suggest that upregulation of ACVR2B and its assembly could affect osteoarthritic cartilage destruction. It is shown that the ACVR2B ligand, activin A, regulates catabolic factor expression through ACVR2B in OA development. Activin A Tg mice (Col2a1-Inhba) exhibit enhanced cartilage destruction, whereas heterozygous activin A KO mice (Inhba+/- ) show protection from cartilage destruction. In silico analysis suggests that the Activin A-ACVR2B axis is involved in Nox4-dependent ROS production. Activin A Tg:Nox4 KO (Col2a1-Inhba:Nox4-/- ) mice show inhibition of experimental OA pathogenesis. NOX4 directly binds to the C-terminal binding site on ACVR2B-ACVR1B and amplifies the pathogenic signal for cartilage destruction through SMAD2/3 signaling. Together, the findings reveal that the ACVR2B assembly, which comprises Activin A, ACVR2B, ACVR1B, Nox4, and AP-1-induced HIF-2α, accelerates OA development. Furthermore, it is shown that shRNA-mediated ACVR2B knockdown or trapping ligands of ACVR2B abrogate OA development by competitively disrupting the ACVR2B-Activin A interaction. These results suggest that the ACVR2B assembly is required to amplify osteoarthritic cartilage destruction and could be a potential therapeutic target in efforts to treat OA.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome.

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

Galectin-4 increases the ability of M2 macrophages to enhance antiviral CD4+ T-cell responses.

Galectin-4 (Gal-4) is a ß-galactoside-binding protein belonging to the galectin family. Although Gal-4 is known to be involved in several physiologic processes of the gastrointestinal tract, its immunomodulatory roles remain unclear. In this study, we investigated whether Gal-4 influences the function of M1 and M2 macrophages. Gal-4 treatment drove more robust changes in the gene expression of M2 macrophages compared to M1 macrophages. Antiviral immune response-related genes were significantly upregulated in Gal-4-treated M2 macrophages. Gal-4 significantly enhanced the immunostimulatory activity of M2 macrophages upon Toll-like receptor 7 stimulation or infection with lymphocytic choriomeningitis virus (LCMV). Moreover, the antibody production against LCMV infection and the antiviral CD4+ T-cell responses, but not the antiviral CD8+ T-cell responses, were greatly increased by Gal-4-treated M2 macrophages in vivo. The present results indicate that Gal-4 enhances the ability of M2 macrophages to promote antiviral CD4+ T-cell responses. Thus, Gal-4 could be used to boost antiviral immune responses.

Genomic and functional features of yeast species in Korean traditional fermented alcoholic beverage and soybean products.

In this review, we describe the genomic and physiological features of the yeast species predominantly isolated from Nuruk, a starter for traditional Korean rice wines, and Jang, a traditional Korean fermented soy product. Nuruk and Jang have several prevalent yeast species, including Saccharomycopsis fibuligera, Hyphopichia burtonii, and Debaryomyces hansenii complex, which belong to the CUG clade showing high osmotic tolerance. Comparative genomics revealed that the interspecies hybridization within yeast species for generating heterozygous diploid genomes occurs frequently as an evolutional strategy in the fermentation environment of Nuruk and Jang. Through gene inventory analysis based on the high-quality reference genome of S. fibuligera, new genes involved in cellulose degradation and volatile aroma biosynthesis and applicable to the production of novel valuable enzymes and chemicals can be discovered. The integrated genomic and transcriptomic analysis of Hyphopichia yeasts, which exhibit strong halotolerance, provides insights into the novel mechanisms of salt and osmo-stress tolerance for survival in fermentation environments with a low-water activity and high-concentration salts. In addition, Jang yeast isolates, such as D. hansenii, show probiotic potential for the industrial application of yeast species beyond fermentation starters to diverse human health sectors.

A New Branchipolynoe (Aphroditiformia: Polynoidae) Scale Worm from the Onnuri Deep-sea Hydrothermal Vent Field, Northern Central Indian Ridge.

Deep-sea hydrothermal vents are dynamic environments with exotic fauna, including bathymodiolin mussels and scale worm annelids that are often in close association. In this study, we found a new species of Branchipolynoe (Aphroditiformia: Polynoidae) living in the recently discovered mussel Gigantidas vrijenhoeki in deep-sea hydrothermal vents and methane seeps at 2,014-2,023 m depth. Based on the morphology and full mitochondrial genome sequences of specimens of Branchipolynoe from the Onnuri vent field (OVF) on the northern Central Indian Ridge, we describe them as a new species: Branchipolynoe onnuriensis sp. nov. This species resembles B. longqiensis and B. tjiasmantoi, but can be distinguished from these species by the shape of the notopodial acicular lobe and the tips of the subacicular neurochaetae. This identity is well-supported by genetic distance and phylogenetic analyses based on the mitochondrial c oxidase subunit I (COI) gene, with the new species being closest to the Western Pacific species B. tjiasmantoi. Phylogenetic analyses support close relationships between the Indian Ocean and Western Pacific hydrothermal polychaetes. Our data provide a foundation for exploring the evolutionary relationship between scale worms and bathymodiolin mussels.

The complete mitochondrial genome of Labidocera rotunda Mori, 1929 (Copepoda: Calanoida) from Jeju Island, Korea.

We sequenced the complete mitochondrial genome of the copepod Labidocera rotunda (family Pontellidae) collected from Ihotaewoo Beach in Jeju, Korea. The mitochondrial genome was 16,564 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNAs, and two ribosomal RNAs. The concatenated phylogenetic tree of L. rotunda was reconstructed using the maximum-likelihood method based on the eight PCGs obtained from eight species of copepods including L. rotunda. The results of the phylogeny analysis showed that L. rotunda was closely related to the family Temoridae among the three families. The complete mitochondrial genome of L. rotunda analyzed for the first time in this study provides insight into the phylogenetic and evolutionary relationship of Labidocera.