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.

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.

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.

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.

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.

Genome-wide identification of nuclear receptor (NR) genes and the evolutionary significance of the NR1O subfamily in the monogonont rotifer Brachionus spp.

Nuclear receptors (NRs) are a large family of transcription factors that are involved in many fundamental biological processes. NRs are considered to have originated from a common ancestor, and are highly conserved throughout the whole animal taxa. Therefore, the genome-wide identification of NR genes in an animal taxon can provide insight into the evolutionary tendencies of NRs. Here, we identified all the NR genes in the monogonont rotifer Brachionus spp., which are considered an ecologically key species due to their abundance and world-wide distribution. The NR family was composed of 40, 32, 29, and 32 genes in the genomes of the rotifers B. calyciflorus, B. koreanus, B. plicatilis, and B. rotundiformis, respectively, which were classified into seven distinct subfamilies. The composition of each subfamily was highly conserved between species, except for NR1O genes, suggesting that they have undergone sporadic evolutionary processes for adaptation to their different environmental pressures. In addition, despite the dynamics of NR evolution, the significance of the conserved endocrine system, particularly for estrogen receptor (ER)-signaling, in rotifers was discussed on the basis of phylogenetic analyses. The results of this study may help provide a better understanding the evolution of NRs, and expand our knowledge of rotifer endocrine systems.

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.

Identification and conservation of gene loss events of Hox gene clusters in the marine medaka (Oryzias melastigma).

In this study, the identification of the whole Hox gene clusters (46 Hox genes) in the marine medaka Oryzias melastigma was investigated using genome assembly and RNA-seq information. Moreover, the gene loss events of Hox gene clusters, which may occur during fish evolution, were examined for a better understanding of the evolutionary status of the gene lost events of the Hox gene cluster across fish species, particularly in the genus Oryzias.

Identification of Hox genes and rearrangements within the single homeobox (Hox) cluster (192.8 kb) of the cyclopoid copepod (Paracyclopina nana).

We report the first identification of the entire complement of the eight typical homeobox (hox) genes (lab, pb, Dfd, scr, antp, ubx, Abd-A, and Abd-B) and the ftz gene in a 192.8 kb region in the cyclopoid copepod Paracyclopina nana. A Hox3 gene ortholog was not present in the P. nana hox gene cluster, while the P. nana Dfd gene was transcribed in the opposite direction to the Daphnia pulex Dfd gene, but in the same direction as the Dfd genes of the fruit fly Drosophila melanogaster and red flour beetle Tribolium castaneum. The location of the lab and pb genes was switched in the P. nana hox cluster, while the order of the remaining hox genes was generally conserved with those of other arthropods. J. Exp. Zool. (Mol. Dev. Evol.) 9999B:XX-XX, 2016. © 2016 Wiley Periodicals, Inc.

Adverse Effects, Expression of the Bk-CYP3045C1 Gene, and Activation of the ERK Signaling Pathway in the Water Accommodated Fraction-Exposed Rotifer.

To examine the deleterious effects of the water accommodated fraction (WAF) of crude oil, the growth curve, fecundity, and lifespan of the monogonont rotifer (Brachionus koreanus) were measured for 24 h in response to three different doses (0.2×, 0.4×, and 0.8×) of WAFs. A higher dose of WAFs significantly reduced the fecundity and lifespan. A rotifer 32K microarray chip showed that the Bk-CYP3045C1 gene had the highest expression. Of the 25 entire CYP genes, the Bk-CYP3045C1 gene showed a significant expression for different doses and times in response to WAFs and chemical components of WAFs (naphthalene and phenanthrene); also, glutathione S-transferase genes, ABC transporter, and other genes showed dose responses upon exposure to 80% WAF over time. Different doses of WAFs increased the oxidative stress with an induction of reactive oxygen species (ROS) and a depletion of glutathione (GSH). Exposure to WAFs did not show toxic effects on survivability in B. koreanus; however, toxicity to WAFs was shown when piperonyl butoxide, a potent inhibitor of cytochrome P450 (CYP) enzymes, was added. This toxicity was dose-dependent. After WAFs exposure, p-ERK was activated over time in response to WAFs, which suggests that WAFs can be activated by the p-ERK signaling pathway.

Identification of the Full 46 Cytochrome P450 (CYP) Complement and Modulation of CYP Expression in Response to Water-Accommodated Fractions of Crude Oil in the Cyclopoid Copepod Paracyclopina nana.

The 46 cytochrome P450 (CYP) gene superfamily was identified in the marine copepod Paracyclopina nana after searching an RNA-seq database and comparing it with other copepod CYP gene families. To annotate the 46 Pn-CYP genes, a phylogenetic analysis of CYP genes was performed using a Bayesian method. Pn-CYP genes were separated into five different clans: CYP2, CYP3, CYP20, CYP26, and mitochondrial. Among these, the principal Pn-CYP genes involved in detoxification were identified by comparing them with those of the copepod Tigriopus japonicus and were examined with respect to their responses to exposure to a water-accommodated fraction (WAF) of crude oil and to the alkylated forms of two polycyclic aromatic hydrocarbons (PAHs; phenanthrene and fluorene). The expression of two Pn-CYP3027 genes (CYP3027F1 and CYP3027F2) was increased in response to WAF exposure and also was upregulated in response to the two alkylated PAHs. In particular, Pn-CYP3027F2 showed the most notable increase in response to 80% WAF exposure. These two responsive CYP genes (Pn-CYP3027F1 and CYP3027F2) were also phylogenetically clustered into the same clade of the WAF- and alkylated PAH-specific CYP genes of the copepod T. japonicus, suggesting that these CYP genes would be those chiefly involved in detoxification in response to WAF exposure in copepods. In this paper, we provide information on the copepod P. nana CYP gene superfamily and also speculate on its potential role in the detoxification of PAHs in marine copepods. Despite the nonlethality of WAF, Pn-CYP3027F2 was rapidly and significantly upregulated in response to WAF that may serve as a useful biomarker of 40% or higher concentration of WAF exposure. This paper will be helpful to better understand the molecular mechanistic events underlying the metabolism of environmental toxicants in copepods.

Genome-wide identification of nuclear receptor (NR) superfamily genes in the copepod Tigriopus japonicus.

BACKGROUND: Nuclear receptors (NRs) are a large superfamily of proteins defined by a DNA-binding domain (DBD) and a ligand-binding domain (LBD). They function as transcriptional regulators to control expression of genes involved in development, homeostasis, and metabolism. The number of NRs differs from species to species, because of gene duplications and/or lineage-specific gene losses during metazoan evolution. Many NRs in arthropods interact with the ecdysteroid hormone and are involved in ecdysone-mediated signaling in arthropods. The nuclear receptor superfamily complement has been reported in several arthropods, including crustaceans, but not in copepods. We identified the entire NR repertoire of the copepod Tigriopus japonicus, which is an important marine model species for ecotoxicology and environmental genomics. RESULTS: Using whole genome and transcriptome sequences, we identified a total of 31 nuclear receptors in the genome of T. japonicus. Nomenclature of the nuclear receptors was determined based on the sequence similarities of the DNA-binding domain (DBD) and ligand-binding domain (LBD). The 7 subfamilies of NRs separate into five major clades (subfamilies NR1, NR2, NR3, NR4, and NR5/6). Although the repertoire of NR members in, T. japonicus was similar to that reported for other arthropods, there was an expansion of the NR1 subfamily in Tigriopus japonicus. The twelve unique nuclear receptors identified in T. japonicus are members of NR1L. This expansion may be a unique lineage-specific feature of crustaceans. Interestingly, E78 and HR83, which are present in other arthropods, were absent from the genomes of T. japonicus and two congeneric copepod species (T. japonicus and Tigriopus californicus), suggesting copepod lineage-specific gene loss. CONCLUSIONS: We identified all NR receptors present in the copepod, T. japonicus. Knowledge of the copepod nuclear receptor repertoire will contribute to a better understanding of copepod- and crustacean-specific NR evolution.

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.

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.

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.

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.

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.

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.

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.