Exosome proteomes reveal glycolysis-related enzyme enrichment in primary canine mammary gland tumor compared to metastases.

OBJECTIVE: Numerous evidence has highlighted the differences between primary tumors and metastases. Nonetheless, the differences in exosomal proteins derived from primary tumor and metastases remain elusive. Here, we aimed to identify differentially expressed exosomal proteins from primary canine mammary gland tumor and metastases to understand how they shape their own tumor microenvironment. METHODS: We clearly distinguished primary canine mammary gland tumors (CHMp) from metastases (CHMm) and profiled the proteins within their secreted exosomes using LC-MS/MS. Moreover, the abundance of glycolysis enzymes (GPI, LDHA) in CHMp exosome was verified with Western blotting, To broaden the scope, we extended to human colorectal cancer-derived exosomes (SW480 vs. SW620) for comparison. RESULTS: We identified significant differences in 87 and 65 proteins derived from CHMp and CHMm, respectively. Notably, glycolysis enzymes (GPI, LDHA, LDHB, TPI1, and ALDOA) showed specific enrichment in exosomes from the primary tumor. CONCLUSION: We observed significant differences in the cellular proteome between primary tumors and metastases, and intriguingly, we identified a parallel heterogeneity the protein composition of exosomes. Specifically, we reported that glycolysis enzymes were significantly enriched in CHMp exosomes compared to CHMm exosomes. We further demonstrated that this quantitative difference in glycolysis enzymes persisted across primary and metastases, extending to human colorectal cancer-derived exosomes (SW480 vs. SW620). Our findings of the specific enrichment of glycolysis enzymes in primary tumor-derived exosomes contribute to a better understanding of tumor microenvironment modulation and heterogeneity between primary tumors and metastases.

ARID3C Acts as a Regulator of Monocyte-to-Macrophage Differentiation Interacting with NPM1.

ARID3C is a protein located on human chromosome 9 and expressed at low levels in various organs, yet its biological function has not been elucidated. In this study, we investigated both the cellular localization and function of ARID3C. Employing a combination of LC-MS/MS and deep learning techniques, we identified NPM1 as a binding partner for ARID3C's nuclear shuttling. ARID3C was found to predominantly localize with the nucleus, where it functioned as a transcription factor for genes STAT3, STAT1, and JUNB, thereby facilitating monocyte-to-macrophage differentiation. The precise binding sites between ARID3C and NPM1 were predicted by AlphaFold2. Mutating this binding site prevented ARID3C from interacting with NPM1, resulting in its retention in the cytoplasm instead of translocation to the nucleus. Consequently, ARID3C lost its ability to bind to the promoters of target genes, leading to a loss of monocyte-to-macrophage differentiation. Collectively, our findings indicate that ARID3C forms a complex with NPM1 to translocate to the nucleus, acting as a transcription factor that promotes the expression of the genes involved in monocyte-to-macrophage differentiation.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus.

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

KOREF_S1: phased, parental trio-binned Korean reference genome using long reads and Hi-C sequencing methods.

BACKGROUND: KOREF is the Korean reference genome, which was constructed with various sequencing technologies including long reads, short reads, and optical mapping methods. It is also the first East Asian multiomic reference genome accompanied by extensive clinical information, time-series and multiomic data, and parental sequencing data. However, it was still not a chromosome-scale reference. Here, we updated the previous KOREF assembly to a new chromosome-level haploid assembly of KOREF, KOREF_S1v2.1. Oxford Nanopore Technologies (ONT) PromethION, Pacific Biosciences HiFi-CCS, and Hi-C technology were used to build the most accurate East Asian reference assembled so far. RESULTS: We produced 705 Gb ONT reads and 114 Gb Pacific Biosciences HiFi reads, and corrected ONT reads by Pacific Biosciences reads. The corrected ultra-long reads reached higher accuracy of 1.4% base errors than the previous KOREF_S1v1.0, which was mainly built with short reads. KOREF has parental genome information, and we successfully phased it using a trio-binning method, acquiring a near-complete haploid-assembly. The final assembly resulted in total length of 2.9 Gb with an N50 of 150 Mb, and the longest scaffold covered 97.3% of GRCh38's chromosome 2. In addition, the final assembly showed high base accuracy, with <0.01% base errors. CONCLUSIONS: KOREF_S1v2.1 is the first chromosome-scale haploid assembly of the Korean reference genome with high contiguity and accuracy. Our study provides useful resources of the Korean reference genome and demonstrates a new strategy of hybrid assembly that combines ONT's PromethION and PacBio's HiFi-CCS.

Enhanced CO2 electroconversion of Rhodobacter sphaeroides by cobalt-phosphate complex assisted water oxidation.

CO2 can be a next generation feedstock for electricity-driven bioproduction due to its abundance and availability. Microbial electrosynthesis (MES), a promising technique for CO2 electroconversion, provides an attractive route for the production of valuable products from CO2, but issues surrounding efficiency and reasonable productivity should be resolved. Improving the anode performance for water oxidation under neutral pH is one of the most important aspects to advance current MES. Here, we introduce cobalt-phosphate (Co-Pi) assisted water oxidation at the counter electrode (i.e., anode) to upgrade the MES performance at pH 7.0. We show that CO2 can be converted by photochemoautotrophic bacterium, Rhodobacter sphaeroides into organic acids and carotenoids in the MES reactor. Planktonic cells of R. sphareroides in the Co-Pi anode equipped MES reactor was ca. 1.5-fold higher than in the control condition (w/o Co-Pi). The faradaic efficiency of the Co-Pi anode equipped MES reactor was remarkably higher (58.3%) than that of the bare anode (27.8%). While the system can improve the CO2 electroconversion nonetheless there are some further optimizations are necessary.

Regulation of Reactive Oxygen Species Promotes Growth and Carotenoid Production Under Autotrophic Conditions in Rhodobacter sphaeroides.

Industrial demand for capture and utilization using microorganisms to reduce CO2, a major cause of global warming, is significantly increasing. Rhodobacter sphaeroides is a suitable strain for the process of converting CO2 into high-value materials because it can accept CO2 and has various metabolic pathways. However, it has been mainly studied for heterotrophic growth that uses sugars and organic acids as carbon sources, not autotrophic growth. Here, we report that the regulation of reactive oxygen species is critical for growth when using CO2 as a sole carbon source in R. sphaeroides. In general, the growth rate is much slower under autotrophic conditions compared to heterotrophic conditions. To improve this, we performed random mutagenesis using N-methyl-N'-nitro-N-nitrosoguanidine (NTG). As a result, we selected the YR-1 strain with a maximum specific growth rate (µ) 1.44 day-1 in the early growth phase, which has a 110% faster growth rate compared to the wild-type. Based on the transcriptome analysis, it was confirmed that the growth was more sensitive to reactive oxygen species under autotrophic conditions. In the YR-1 mutant, the endogenous contents of H2O2 levels and oxidative damage were reduced by 33.3 and 42.7% in the cells, respectively. Furthermore, we measured that concentrations of carotenoids, which are important antioxidants. The total carotenoid is produced 9.63 g/L in the YR-1 mutant, suggesting that the production is 1.7-fold higher than wild-type. Taken together, our observations indicate that controlling ROS promotes cell growth and carotenoid production under autotrophic conditions.

LT1, an ONT long-read-based assembly scaffolded with Hi-C data and polished with short reads.

We present LT1, the first high-quality human reference genome from the Baltic States. LT1 is a female de novo human reference genome assembly, constructed using 57× nanopore long reads and polished using 47× short paired-end reads. We utilized 72 GB of Hi-C chromosomal mapping data for scaffolding, to maximize assembly contiguity and accuracy. The contig assembly of LT1 was 2.73 Gbp in length, comprising 4490 contigs with an NG50 value of 12.0 Mbp. After scaffolding with Hi-C data and manual curation, the final assembly has an NG50 value of 137 Mbp and 4699 scaffolds. Assessment of gene prediction quality using Benchmarking Universal Single-Copy Orthologs (BUSCO) identified 89.3% of the single-copy orthologous genes included in the benchmark. Detailed characterization of LT1 suggests it has 73,744 predicted transcripts, 4.2 million autosomal SNPs, 974,616 short indels, and 12,079 large structural variants. These data may be used as a benchmark for further in-depth genomic analyses of Baltic populations.

Comparative analysis of 7 short-read sequencing platforms using the Korean Reference Genome: MGI and Illumina sequencing benchmark for whole-genome sequencing.

BACKGROUND: DNBSEQ-T7 is a new whole-genome sequencer developed by Complete Genomics and MGI using DNA nanoball and combinatorial probe anchor synthesis technologies to generate short reads at a very large scale-up to 60 human genomes per day. However, it has not been objectively and systematically compared against Illumina short-read sequencers. FINDINGS: By using the same KOREF sample, the Korean Reference Genome, we have compared 7 sequencing platforms including BGISEQ-500, DNBSEQ-T7, HiSeq2000, HiSeq2500, HiSeq4000, HiSeqX10, and NovaSeq6000. We measured sequencing quality by comparing sequencing statistics (base quality, duplication rate, and random error rate), mapping statistics (mapping rate, depth distribution, and percent GC coverage), and variant statistics (transition/transversion ratio, dbSNP annotation rate, and concordance rate with single-nucleotide polymorphism [SNP] genotyping chip) across the 7 sequencing platforms. We found that MGI platforms showed a higher concordance rate for SNP genotyping than HiSeq2000 and HiSeq4000. The similarity matrix of variant calls confirmed that the 2 MGI platforms have the most similar characteristics to the HiSeq2500 platform. CONCLUSIONS: Overall, MGI and Illumina sequencing platforms showed comparable levels of sequencing quality, uniformity of coverage, percent GC coverage, and variant accuracy; thus we conclude that the MGI platforms can be used for a wide range of genomics research fields at a lower cost than the Illumina platforms.

Chromosome-scale assembly comparison of the Korean Reference Genome KOREF from PromethION and PacBio with Hi-C mapping information.

BACKGROUND: Long DNA reads produced by single-molecule and pore-based sequencers are more suitable for assembly and structural variation discovery than short-read DNA fragments. For de novo assembly, Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) are the favorite options. However, PacBio's SMRT sequencing is expensive for a full human genome assembly and costs more than $40,000 US for 30× coverage as of 2019. ONT PromethION sequencing, on the other hand, is 1/12 the price of PacBio for the same coverage. This study aimed to compare the cost-effectiveness of ONT PromethION and PacBio's SMRT sequencing in relation to the quality. FINDINGS: We performed whole-genome de novo assemblies and comparison to construct an improved version of KOREF, the Korean reference genome, using sequencing data produced by PromethION and PacBio. With PromethION, an assembly using sequenced reads with 64× coverage (193 Gb, 3 flowcell sequencing) resulted in 3,725 contigs with N50s of 16.7 Mb and a total genome length of 2.8 Gb. It was comparable to a KOREF assembly constructed using PacBio at 62× coverage (188 Gb, 2,695 contigs, and N50s of 17.9 Mb). When we applied Hi-C-derived long-range mapping data, an even higher quality assembly for the 64× coverage was achieved, resulting in 3,179 scaffolds with an N50 of 56.4 Mb. CONCLUSION: The pore-based PromethION approach provided a high-quality chromosome-scale human genome assembly at a low cost with long maximum contig and scaffold lengths and was more cost-effective than PacBio at comparable quality measurements.

The genome of the marine monogonont rotifer Brachionus plicatilis: Genome-wide expression profiles of 28 cytochrome P450 genes in response to chlorpyrifos and 2-ethyl-phenanthrene.

Brachionus spp. (Rotifera: Monogononta) are globally distributed in aquatic environments and play important roles in the aquatic ecosystem. The marine monogonont rotifer Brachionus plicatilis is considered a suitable model organism for ecology, evolution, and ecotoxicology. In this study, we assembled and characterized the B. plicatilis genome. The total length of the assembled genome was 106.9 Mb and the number of final scaffolds was 716 with an N50 value of 1.15 Mb and a GC content of 26.75%. A total of 20,154 genes were annotated after manual curation. To demonstrate the use of whole genome data, we targeted one of the main detoxifying enzyme of phase I detoxification system and identified in a total of 28 cytochrome P450 s (CYPs). Based on the phylogenetic analysis using the maximum likelihood, 28 B. plicatilis-CYPs were apparently separated into five different clans, namely, 2, 3, 4, mitochondrial (MT), and 46 clans. To better understand the CYPs-mediated xenobiotic detoxification, we measured the mRNA expression levels of 28 B. plicatilis CYPs in response to chlorpyrifos and 2-ethyl-phenanthrene. Most B. plicatilis CYPs were significantly modulated (P < 0.05) in response to chlorpyrifos and 2-ethyl-phenanthrene. In addition, xenobiotic-sensing nuclear receptor (XNR) response element sequences were identified in the 5 kb upstream of promoter regions of 28 CYPs from the genome of B. plicatilis, indicating that these XNR can be associated with detoxification of xenobiotics. Overall, the assembled B. plicatilis genome presented here will be a useful resource for a better understanding the molecular ecotoxicology in the view of molecular mechanisms underlying toxicological responses, particularly on xenobiotic detoxification in this species.

The Draft Genome of an Octocoral, Dendronephthya gigantea.

Coral reefs composed of stony corals are threatened by global marine environmental changes. However, soft coral communities of octocorallian species, appear more resilient. The genomes of several cnidarians species have been published, including from stony corals, sea anemones, and hydra. To fill the phylogenetic gap for octocoral species of cnidarians, we sequenced the octocoral, Dendronephthya gigantea, a nonsymbiotic soft coral, commonly known as the carnation coral. The D. gigantea genome size is ∼276 Mb. A high-quality genome assembly was constructed from PacBio long reads (29.85 Gb with 108× coverage) and Illumina short paired-end reads (35.54 Gb with 128× coverage) resulting in the highest N50 value (1.4 Mb) reported thus far among cnidarian genomes. About 12% of the genome is repetitive elements and contained 28,879 predicted protein-coding genes. This gene set is composed of 94% complete BUSCO ortholog benchmark genes, which is the second highest value among the cnidarians, indicating high quality. Based on molecular phylogenetic analysis, octocoral and hexacoral divergence times were estimated at 544 MYA. There is a clear difference in Hox gene composition between these species: unlike hexacorals, the Antp superclass Evx gene was absent in D. gigantea. Here, we present the first genome assembly of a nonsymbiotic octocoral, D. gigantea to aid in the comparative genomic analysis of cnidarians, including stony and soft corals, both symbiotic and nonsymbiotic. The D. gigantea genome may also provide clues to mechanisms of differential coping between the soft and stony corals in response to scenarios of global warming.

The complete mitochondrial genome of a Dokdo shrimp, Lebbeus groenlandicus.

Lebbeus groenlandicus is a shrimp species indigenous to the Dokdo islands in the East Sea of Korea. We report the 17,399 bp mitochondrial genome (mitogenome) of the species that consists of 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region (CR). A maximum-likelihood tree, constructed with 18 prawn and 45 shrimp mitogenomes, confirmed that L. groenlandicus occupies the most basal position within the Caridea infra-order and is closely related to Pandalidae shrimps.

The genome of the marine medaka Oryzias melastigma.

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field-based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole-genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.

The genome of the freshwater monogonont rotifer Brachionus calyciflorus.

Monogononta is the most speciose class of rotifers, with more than 2,000 species. The monogonont genus Brachionus is widely distributed at a global scale, and a few of its species are commonly used as ecological and evolutionary models to address questions related to aquatic ecology, cryptic speciation, evolutionary ecology, the evolution of sex and ecotoxicology. With the importance of Brachionus species in many areas of research, it is remarkable that the genome has not been characterized. This study aims to address this lacuna by presenting, for the first time, the whole-genome assembly of the freshwater species Brachionus calyciflorus. The total length of the assembled genome was 129.6 Mb, with 1,041 scaffolds. The N50 value was 786.6 kb, and the GC content was 24%. A total of 16,114 genes were annotated with repeat sequences, accounting for 21% of the assembled genome. This assembled genome may form a basis for future studies addressing key questions on the evolution of monogonont rotifers. It will also provide the necessary molecular resources to mechanistically investigate ecophysiological and ecotoxicological responses.

Identification and characterization of homeobox (Hox) genes and conservation of the single Hox cluster (324.6 kb) in the water flea Daphnia magna.

We report the complete sequence analysis of the entire complement of eight typical homeobox (Hox) genes (Lab, Pb, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) and two other genes (Hox3 and Ftz) in a 324.6-kb region in the water flea Daphnia magna. In the cluster of D. magna Hox genes, we found one long interspersed nuclear element (LINE)/R2-NeSL between Ubx and Abd-A that was not present in Daphnia pulex Hox genes. In basal expression of Hox genes at different developmental stages, biothorax complex genes (Ubx, Abd-A, and Abd-B) and some antennapedia complex genes (Lab, Scr, Antp) were moderately expressed, but the Hox3 gene was barely expressed. Three homeobox genes (Antp, Ubx, Abd-A) were highly expressed at 6-7 days after release from the brood chamber and/or in the adult stage. The structural array and transcribed orientation of Dm-Hox genes were identical to those of the sister species D. pulex (∼340 kb), indicating that the Hox gene structure in daphnids is highly conserved. However, Dm- and Dp-Hox3, -deformed (Dfd), and -fushi tarazu (Ftz) genes varied from orthologous genes in pancrustacean species.

Identification of 74 cytochrome P450 genes and co-localized cytochrome P450 genes of the CYP2K, CYP5A, and CYP46A subfamilies in the mangrove killifish Kryptolebias marmoratus.

BACKGROUND: The mangrove killifish Kryptolebias marmoratus is the only vertebrate that reproduces by self-fertilizing and is an important model species in genetics and marine ecotoxicology. Using whole-genome and transcriptome sequences, we identified all members of the cytochrome P450 (CYP) family in this model teleost and compared them with those of other teleosts. RESULTS: A total of 74 cytochrome P450 genes and one pseudogene were identified in K. marmoratus. Phylogenetic analysis indicated that the CYP genes in clan 2 were most expanded, while synteny analysis with other species showed orthologous relationships of CYP subfamilies among teleosts. In addition to the CYP2K expansions, five tandem duplicated gene copies of CYP5A were observed. These features were unique to K. marmoratus. CONCLUSIONS: These results shed a light on CYP gene evolution, particularly the co-localized CYP2K, CYP5A, and CYP46A subfamilies in fish. Future studies of CYP expression could identify specific endogenous and exogenous environmental factors that triggered the evolution of tandem CYP duplication in K. marmoratus.

Genome-wide identification of 31 cytochrome P450 (CYP) genes in the freshwater rotifer Brachionus calyciflorus and analysis of their benzo[α]pyrene-induced expression patterns.

While marine invertebrate cytochrome P450 (CYP) genes and their roles in detoxification mechanisms have been studied, little information is available regarding freshwater rotifer CYPs and their functions. Here, we used genomic sequences and RNA-seq databases to identify 31 CYP genes in the freshwater rotifer Brachionus calyciflorus. The 31 Bc-CYP genes with a few tandem duplications were clustered into CYP 2, 3, 4, mitochondrial, and 46 clans with two marine rotifers Brachionus plicatilis and Brachionus koreanus. To understand the molecular responses of these 31 Bc-CYP genes, we also examined their expression patterns in response to benzo[α]pyrene (B[α]P). Three Bc-CYP genes (Bc-CYP3044B3, Bc-CYP3049B4, Bc-CYP3049B6) were significantly upregulated (P<0.05) in response to B[α]P, suggesting that these CYP genes can be involved in detoxification in response to B[α]P exposure. These genes might be useful as biomarkers of B[α]P exposure in B. calyciflorus. Overall, our findings expand the repertoire of known CYPs and shed light on their potential roles in xenobiotic detoxification in rotifers.

Modification of the plasma complement protein profile by exogenous estrogens is indicative of a compromised immune competence in marine medaka (Oryzias melastigma).

Growing evidence suggests that the immune system of teleost is vulnerable to xenoestrogens, which are ubiquitous in the marine environment. This study detected and identified the major circulatory immune proteins deregulated by 17α-ethinylestradiol (EE2), which may be linked to fish susceptibility to pathogens in the marine medaka, Oryzias melastigma. Fish immune competence was determined using a host resistance assay to pathogenic bacteria Edwardsiella tarda. Females were consistently more susceptible to infection-induced mortality than males. Exposure to EE2 could narrow the sex gap of mortality by increasing infection-induced death in male fish. Proteomic analysis revealed that the major plasma immune proteins of adult fish were highly sexually dimorphic. EE2 induced pronounced sex-specific changes in the plasma proteome, with the male plasma composition clearly becoming "feminised". Male plasma was found to contain a higher level of fibrinogens, WAP63 and ependymin-2-like protein, which are involved in coagulation, inflammation and regeneration. For the first time, we demonstrated that expression of C1q subunit B (C1Q), an initiating factor of the classical complement pathway, was higher in males and was suppressed in both sexes in response to EE2 and bacterial challenge. Moreover, cleavage and post-translational modification of C3, the central component of the complement system, could be altered by EE2 treatment in males (C3dg down; C3g up). Multiple regression analysis indicated that C1Q is possibly an indicator of fish survival, which warrants further confirmation. The findings support the potential application of plasma immune proteins for prognosis/diagnosis of fish immune competence. Moreover, this study provides the first biochemical basis of the sex-differences in fish immunity and how these differences might be modified by xenoestrogens.

Effects of water accommodated fractions (WAFs) of crude oil in two congeneric copepods Tigriopus sp.

Oil pollution has deleterious effects on marine ecosystems. However, the toxicity of crude oil towards Antarctic marine organisms has not been well studied. We compared the deleterious effects of water accommodated fractions (WAFs) of crude oil on reproduction, intracellular reactive oxygen species (ROS) levels, and antioxidant enzymatic activity in Antarctic (Tigriopus kingsejongensis) and temperate (Tigriopus japonicus) copepods. Reproductive rates of T. kingsejongensis and T. japonicus were significantly reduced (P < 0.05) in response to WAFs. Furthermore, T. kingsejongensis showed elevated levels of ROS and higher antioxidant enzyme (glutathione peroxidase [GPx]) activity than T. japonicus in response to WAFs. CYP genes from congeneric copepods were identified and annotated to better understand molecular detoxification mechanisms. We observed significant up-regulation (P < 0.05) of Tk-CYP3024A3 and Tj-CYP3024A2 in response to WAFs, suggesting that CYP genes may contribute to the detoxification mechanism in response to WAF exposure. These finding also suggest that WAFs may induce oxidative stress, leading to reproductive impairment in copepods. Furthermore, Tk-CYP3024A3 and Tj-CYP3024A2 genes can be considered as potential biomarkers of WAF toxicity in the congeneric copepods T. kingsejongensis and T. japonicus. This study will be helpful for enhancing our knowledge on the harmful effects of WAFs in Antarctic and temperate copepods and provides insight into the underlying molecular mechanisms.

Genome-wide identification of 52 cytochrome P450 (CYP) genes in the copepod Tigriopus japonicus and their B[α]P-induced expression patterns.

Cytochrome P450s (CYPs) are enzymes with a heme-binding domain that are found in all living organisms. CYP enzymes have important roles associated with detoxification of xenobiotics and endogenous compounds (e.g. steroids, fatty acids, and hormones). Although CYP enzymes have been reported in several invertebrates, including insects, little is known about copepod CYPs. Here, we identified the entire repertoire of CYP genes (n=52) from whole genome and transcriptome sequences of the benthic copepod Tigriopus japonicus, including a tandem duplication (CYP3026A3, CYP3026A4, CYP3026A5), and examined patterns of gene expression over various developmental stages and in response to benzo[α]pyrene (B[α]P) exposure. Through phylogenetic analysis, the 52 T. japonicus CYP genes were assigned to five distinct clans: CYP2 (22 genes), CYP3 (19 genes), CYP4 (two genes), CYP20 (one gene), and mitochondrial (eight genes). Developmental stage and gender-specific expression patterns of the 52 T. japonicus CYPs were analyzed. CYP3022A1 was constitutively expressed during all developmental stages. CYP genes in clans 2 and 3 were induced in response to B[α]P, suggesting that these differentially modulated CYP transcripts are likely involved in defense against exposure to B[α]P and other pollutants. This study enhances our understanding of the repertoire of CYP genes in copepods and of their potential role in development and detoxification in copepods.