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.

METTL8 mRNA Methyltransferase Enhances Cancer Cell Migration via Direct Binding to ARID1A.

The association of RNA modification in cancer has recently been highlighted. Methyltransferase like 8 (METTL8) is an enzyme and its role in mRNA m3C modification has barely been studied. In this study, we found that METTL8 expression was significantly up-regulated in canine mammary tumor and investigated its functional roles in the tumor process, including cancer cell proliferation and migration. METTL8 expression was up-regulated in most human breast cancer cell lines tested and decreased by Yin Yang 1 (YY1) transcription factor knockdown, suggesting that YY1 is a regulating transcription factor. The knockdown of METTL8 attenuated tumor cell growth and strongly blocked tumor cell migration. AT-rich interactive domain-containing protein 1A (ARID1A) was identified as a candidate mRNA by METTL8. ARID1A mRNA binds to METTL8 protein. ARID1A mRNA expression was not changed by METTL8 knockdown, but ARID1A protein level was significantly increased. Collectively, our study indicates that METTL8 up-regulated by YY1 in breast cancer plays an important role in cancer cell migration through the mRNA modification of ARID1A, resulting in the attenuation of its translation.

Exploring the key communicator role of exosomes in cancer microenvironment through proteomics.

There have been many attempts to fully understand the mechanism of cancer behavior. Yet, how cancers develop and metastasize still remain elusive. Emerging concepts of cancer biology in recent years have focused on the communication of cancer with its microenvironment, since cancer cannot grow and live alone. Cancer needs to communicate with other cells for survival, and thus they secrete various messengers, including exosomes that contain many proteins, miRNAs, mRNAs, etc., for construction of the tumor microenvironment. Moreover, these intercellular communications between cancer and its microenvironment, including stromal cells or distant cells, can promote tumor growth, metastasis, and escape from immune surveillance. In this review, we summarized the role of proteins in the exosome as communicators between cancer and its microenvironment. Consequently, we present cancer specific exosome proteins and their unique roles in the interaction between cancer and its microenvironment. Clinically, these exosomes might provide useful biomarkers for cancer diagnosis and therapeutic tools for cancer treatment.

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.

Real-Time Longitudinal Evaluation of Tumor Blood Vessels Using a Compact Preclinical Fluorescence Imaging System.

Tumor angiogenesis is enhanced in all types of tumors to supply oxygen and nutrients for their growth and metastasis. With the development of anti-angiogenic drugs, the importance of technology that closely monitors tumor angiogenesis has also been emerging. However, to date, the technology for observing blood vessels requires specialized skills with expensive equipment, thereby limiting its applicability only to the laboratory setting. Here, we used a preclinical optical imaging system for small animals and, for the first time, observed, in real time, the entire process of blood vessel development in tumor-bearing mice injected with indocyanine green. Time-lapse sequential imaging revealed blood vessel volume and blood flow dynamics on a microscopic scale. Upon analyzing fluorescence dynamics at each stage of tumor progression, vessel volume and blood flow were found to increase as the tumor developed. Conversely, these vascular parameters decreased when the mice were treated with angiogenesis inhibitors, which suggests that the effects of drugs targeting angiogenesis can be rapidly and easily screened. The results of this study may help evaluate the efficacy of angiogenesis-targeting drugs by facilitating the observation of tumor blood vessels easily in a laboratory unit without large and complex equipment.

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.

Analysis of opposing histone modifications H3K4me3 and H3K27me3 reveals candidate diagnostic biomarkers for TNBC and gene set prediction combination.

Breast cancer encompasses a major portion of human cancers and must be carefully monitored for appropriate diagnoses and treatments. Among the many types of breast cancers, triple negative breast cancer (TNBC) has the worst prognosis and the least cases reported. To gain a better understanding and a more decisive precursor for TNBC, two major histone modifications, an activating modification H3K4me3 and a repressive modification H3K27me3, were analyzed using data from normal breast cell lines against TNBC cell lines. The combination of these two histone markers on the gene promoter regions showed a great correlation with gene expression. A list of signature genes was defined as active (highly enriched H3K4me3), including NOVA1, NAT8L, and MMP16, and repressive genes (highly enriched H3K27me3), IRX2 and ADRB2, according to the distribution of these histone modifications on the promoter regions. To further enhance the investigation, potential candidates were also compared with other types of breast cancer to identify signs specific to TNBC. RNA-seq data was implemented to confirm and verify gene regulation governed by the histone modifications. Combinations of the biomarkers based on H3K4me3 and H3K27me3 showed the diagnostic value AUC 93.28% with P-value of 1.16e-226. The results of this study suggest that histone modification analysis of opposing histone modifications may be valuable toward developing biomarkers and targets for TNBC. [BMB Reports 2020; 53(5): 266-271].