Comparison of structural variant callers for massive whole-genome sequence data.

BACKGROUND: Detecting structural variations (SVs) at the population level using next-generation sequencing (NGS) requires substantial computational resources and processing time. Here, we compared the performances of 11 SV callers: Delly, Manta, GridSS, Wham, Sniffles, Lumpy, SvABA, Canvas, CNVnator, MELT, and INSurVeyor. These SV callers have been recently published and have been widely employed for processing massive whole-genome sequencing datasets. We evaluated the accuracy, sequence depth, running time, and memory usage of the SV callers. RESULTS: Notably, several callers exhibited better calling performance for deletions than for duplications, inversions, and insertions. Among the SV callers, Manta identified deletion SVs with better performance and efficient computing resources, and both Manta and MELT demonstrated relatively good precision regarding calling insertions. We confirmed that the copy number variation callers, Canvas and CNVnator, exhibited better performance in identifying long duplications as they employ the read-depth approach. Finally, we also verified the genotypes inferred from each SV caller using a phased long-read assembly dataset, and Manta showed the highest concordance in terms of the deletions and insertions. CONCLUSIONS: Our findings provide a comprehensive understanding of the accuracy and computational efficiency of SV callers, thereby facilitating integrative analysis of SV profiles in diverse large-scale genomic datasets.

Diagnosis of Crohn's disease and ulcerative colitis using the microbiome.

BACKGROUND: Inflammatory bowel disease (IBD) is a multifactorial chronic inflammatory disease resulting from dysregulation of the mucosal immune response and gut microbiota. Crohn's disease (CD) and ulcerative colitis (UC) are difficult to distinguish, and differential diagnosis is essential for establishing a long-term treatment plan for patients. Furthermore, the abundance of mucosal bacteria is associated with the severity of the disease. This study aimed to differentiate and diagnose these two diseases using the microbiome and identify specific biomarkers associated with disease activity. RESULTS: Differences in the abundance and composition of the microbiome between IBD patients and healthy controls (HC) were observed. Compared to HC, the diversity of the gut microbiome in patients with IBD decreased; the diversity of the gut microbiome in patients with CD was significantly lower. Sixty-eight microbiota members (28 for CD and 40 for UC) associated with these diseases were identified. Additionally, as the disease progressed through different stages, the diversity of the bacteria decreased. The abundances of Alistipes shahii and Pseudodesulfovibrio aespoeensis were negatively correlated with the severity of CD, whereas the abundance of Polynucleobacter wianus was positively correlated. The severity of UC was negatively correlated with the abundance of A. shahii, Porphyromonas asaccharolytica and Akkermansia muciniphilla, while it was positively correlated with the abundance of Pantoea candidatus pantoea carbekii. A regularized logistic regression model was used for the differential diagnosis of the two diseases. The area under the curve (AUC) was used to examine the performance of the model. The model discriminated UC and CD at an AUC of 0.873 (train set), 0.778 (test set), and 0.633 (validation set) and an area under the precision-recall curve (PRAUC) of 0.888 (train set), 0.806 (test set), and 0.474 (validation set). CONCLUSIONS: Based on fecal whole-metagenome shotgun (WMS) sequencing, CD and UC were diagnosed using a machine-learning predictive model. Microbiome biomarkers associated with disease activity (UC and CD) are also proposed.

Mechanism mediated by a noncoding RNA, nc886, in the cytotoxicity of a DNA-reactive compound.

DNA-reactive compounds are harnessed for cancer chemotherapy. Their genotoxic effects are considered to be the main mechanism for the cytotoxicity to date. Because this mechanism preferentially affects actively proliferating cells, it is postulated that the cytotoxicity is specific to cancer cells. Nonetheless, they do harm normal quiescent cells, suggesting that there are other cytotoxic mechanisms to be uncovered. By employing doxorubicin as a representative DNA-reactive compound, we have discovered a cytotoxic mechanism that involves a cellular noncoding RNA (ncRNA) nc886 and protein kinase R (PKR) that is a proapoptotic protein. nc886 is transcribed by RNA polymerase III (Pol III), binds to PKR, and prevents it from aberrant activation in most normal cells. We have shown here that doxorubicin evicts Pol III from DNA and, thereby, shuts down nc886 transcription. Consequently, the instantaneous depletion of nc886 provokes PKR and leads to apoptosis. In a short-pulse treatment of doxorubicin, these events are the main cause of cytotoxicity preceding the DNA damage response in a 3D culture system as well as the monolayer cultures. By identifying nc886 as a molecular signal for PKR to sense doxorubicin, we have provided an explanation for the conundrum why DNA-damaging drugs can be cytotoxic to quiescent cells that have the competent nc886/PKR pathway.

Nc886, a Novel Suppressor of the Type I Interferon Response Upon Pathogen Intrusion.

Interferons (IFNs) are a crucial component in the innate immune response. Especially the IFN-ß signaling operates in most cell types and plays a key role in the first line of defense upon pathogen intrusion. The induction of IFN-ß should be tightly controlled, because its hyperactivation can lead to tissue damage or autoimmune diseases. Activation of the IFN-ß promoter needs Interferon Regulatory Factor 3 (IRF3), together with Nuclear Factor kappa-light-chain-enhancer of activated B cells (NF-κB) and Activator Protein 1 (AP-1). Here we report that a human noncoding RNA, nc886, is a novel suppressor for the IFN-ß signaling and inflammation. Upon treatment with several pathogen-associated molecular patterns and viruses, nc886 suppresses the activation of IRF3 and also inhibits NF-κB and AP-1 via inhibiting Protein Kinase R (PKR). These events lead to decreased expression of IFN-ß and resultantly IFN-stimulated genes. nc886's role might be to restrict the IFN-ß signaling from hyperactivation. Since nc886 expression is regulated by epigenetic and environmental factors, nc886 might explain why innate immune responses to pathogens are variable depending on biological settings.

Methylation of the CDX2 promoter in Helicobacter pylori-infected gastric mucosa increases with age and its rapid demethylation in gastric tumors is associated with upregulated gene expression.

Pathological changes in the epigenetic landscape of chromatin are hallmarks of cancer. The caudal-type homeobox gene CDX2 is not expressed in normal gastric epithelia but rather in adult intestinal epithelia, and it is overexpressed in intestinal metaplasia (IM). However, it remains unclear how CDX2 transcription is suppressed in normal gastric epithelial cells and overexpressed in IM. Here, we demonstrate that methylation of the CDX2 promoter increases with age in Helicobacter pylori-positive, noncancerous gastric tissue, whereas the promoter is demethylated in paired gastric tumors in which CDX2 is upregulated. Moreover, we also found that the CDX2 promoter is demethylated in IM as well as gastric tumor. Immunohistochemistry revealed that CDX2 is present in foci of parts of the gastric mucosae but highly expressed in IM as well as in gastric tumors, suggesting that the elevated level of CDX2 in IM and gastric tumors may be attributable to promoter demethylation. Our data suggest that CDX2 repression may be associated with promoter methylation in noncancerous H. pylori-positive mucosa but its upregulation might be attributable to increased promoter activity mediated by chromatin remodeling during gastric carcinogenesis.

ONECUT2 upregulation is associated with CpG hypomethylation at promoter-proximal DNA in gastric cancer and triggers ACSL5.

Many studies have focused on global hypomethylation or hypermethylation of tumor suppressor genes, but less is known about the impact of promoter hypomethylation of oncogenes. We previously showed that promoter methylation may gradually increase or decrease during the transition from gastric mucosa (GM) to intestinal metaplasia (IM) to gastric cancer (GC). In our study, we focused on regional CpG hypomethylation of the promoter-proximal DNA of the transcription factor ONECUT2 (OC2) in IM and GC cells. We validated the hypomethylation of promoter-proximal DNA of OC2 in 160 primary GCs, in which methylation level correlated negatively with OC2 mRNA level. IM and GC cells stained positively for OC2, whereas GM cells did not. Stable transfection of OC2 in GC cells promoted colony formation, cell migration, invasion and proliferation. Moreover, OC2 knockdown with a short hairpin RNA suppressed tumorigenesis in nude mice. In addition, chromatin immunoprecipitation coupled with DNA sequencing and RNA-seq analyses revealed that OC2 triggered ACSL5, which is strongly expressed in IM of the stomach but not in GM, indicating that OC2 and ACSL5 are early-stage biomarkers for GC. We also observed a high correlation between the levels of OC2 and ACSL5 mRNAs in the GENT database These results suggest that epigenetic alteration of OC2 upregulates its expression, which then activates ACSL5; thus, OC2 is induced in IM by epigenetic alteration and triggers ACSL5 expression, and thus OC2 and ACSL5 may cooperatively promote intestinal differentiation and GC progression.

Identification of a molecular signature of prognostic subtypes in diffuse-type gastric cancer.

BACKGROUND: Although recent advances in high-throughput technology have provided many insights into gastric cancer (GC), few reliable biomarkers for diffuse-type GC have been identified. Here, we aim to identify a prognostic and predictive signature of diffuse-type GC heterogeneity. METHODS: We analyzed RNA-seq-based transcriptome data to identify a molecular signature in 150 gastric tissue samples including 107 diffuse-type GCs. The predictive value of the signature was verified using other diffuse-type GC samples in three independent cohorts (n = 466). Log-rank and Cox regression analyses were used to estimate the association between the signature and prognosis. The signature was also characterized by somatic variant analyses and tissue microarray analysis between diffuse-type GC subtypes. RESULTS: Transcriptomic profiling of RNA-seq data identified a signature which revealed distinct subtypes of diffuse-type GC: the intestinal-like (INT) and core diffuse-type (COD) subtypes. The signature showed high predictability and independent clinical utility in diffuse-type GC prognosis in other patient cohorts (HR 2.058, 95% CI 1.53-2.77, P = 1.76 × 10-6). Integrative mutational and gene expression analyses demonstrated that the COD subtype was responsive to chemotherapy, whereas the INT subtype was responsive to immunotherapy with an immune checkpoint inhibitor (ICI). Tissue microarray analysis showed the practical utility of IGF1 and NXPE2 for predicting diffuse-type GC heterogeneity. CONCLUSIONS: We present a molecular signature that can identify diffuse-type GC patients who display different clinical behaviors as well as responses to chemotherapy or ICI treatment.

Integrated genomic analysis of recurrence-associated small non-coding RNAs in oesophageal cancer.

OBJECTIVE: Oesophageal squamous cell carcinoma (ESCC) is a heterogeneous disease with variable outcomes that are challenging to predict. A better understanding of the biology of ESCC recurrence is needed to improve patient care. Our goal was to identify small non-coding RNAs (sncRNAs) that could predict the likelihood of recurrence after surgical resection and to uncover potential molecular mechanisms that dictate clinical heterogeneity. DESIGN: We developed a robust prediction model for recurrence based on the analysis of the expression profile data of sncRNAs from 108 fresh frozen ESCC specimens as a discovery set and assessment of the associations between sncRNAs and recurrence-free survival (RFS). We also evaluated the mechanistic and therapeutic implications of sncRNA obtained through integrated analysis from multiple datasets. RESULTS: We developed a risk assessment score (RAS) for recurrence with three sncRNAs (microRNA (miR)-223, miR-1269a and nc886) whose expression was significantly associated with RFS in the discovery cohort (n=108). RAS was validated in an independent cohort of 512 patients. In multivariable analysis, RAS was an independent predictor of recurrence (HR, 2.27; 95% CI, 1.26 to 4.09; p=0.007). This signature implies the expression of ΔNp63 and multiple alterations of driver genes like PIK3CA. We suggested therapeutic potentials of immune checkpoint inhibitors in low-risk patients, and Polo-like kinase inhibitors, mammalian target of rapamycin (mTOR) inhibitors, and histone deacetylase inhibitors in high-risk patients. CONCLUSION: We developed an easy-to-use prognostic model with three sncRNAs as robust prognostic markers for postoperative recurrence of ESCC. We anticipate that such a stratified and systematic, tumour-specific biological approach will potentially contribute to significant improvement in ESCC treatment.

Dynamic changes in DNA methylation and hydroxymethylation when hES cells undergo differentiation toward a neuronal lineage.

DNA methylation and hydroxymethylation have been implicated in normal development and differentiation, but our knowledge is limited about the genome-wide distribution of 5-methylcytosine (5 mC) and 5-hydroxymethylcytosine (5 hmC) during cellular differentiation. Using an in vitro model system of gradual differentiation of human embryonic stem (hES) cells into ventral midbrain-type neural precursor cells and terminally into dopamine neurons, we observed dramatic genome-wide changes in 5 mC and 5 hmC patterns during lineage commitment. The 5 hmC pattern was dynamic in promoters, exons and enhancers. DNA hydroxymethylation within the gene body was associated with gene activation. The neurogenesis-related genes NOTCH1, RGMA and AKT1 acquired 5 hmC in the gene body and were up-regulated during differentiation. DNA methylation in the promoter was associated with gene repression. The pluripotency-related genes POU5F1, ZFP42 and HMGA1 acquired 5 mC in their promoters and were down-regulated during differentiation. Promoter methylation also acted as a locking mechanism to maintain gene silencing. The mesoderm development-related genes NKX2-8, TNFSF11 and NFATC1 acquired promoter methylation during neural differentiation even though they were already silenced in hES cells. Our findings will help elucidate the molecular mechanisms underlying lineage-specific differentiation of pluripotent stem cells during human embryonic development.

Integrative genomics analysis reveals the multilevel dysregulation and oncogenic characteristics of TEAD4 in gastric cancer.

Tumorigenesis is a consequence of failures of multistep defense mechanisms against deleterious perturbations that occur at the genomic, epigenomic, transcriptomic and proteomic levels. To uncover previously unrecognized genes that undergo multilevel perturbations in gastric cancer (GC), we integrated epigenomic and transcriptomic approaches using two recently developed tools: MENT and GENT. This integrative analysis revealed that nine Hippo pathway-related genes, including components [FAT, JUB, LATS2, TEA domain family member 4 (TEAD4) and Yes-associated protein 1 (YAP1)] and targets (CRIM1, CYR61, CTGF and ITGB2), are concurrently hypomethylated at promoter CpG sites and overexpressed in GC tissues. In particular, TEAD4, a link between Hippo pathway components and targets, was significantly hypomethylated at CpG site cg21637033 (P = 3.8 × 10(-) (20)) and overexpressed (P = 5.2 × 10(-) (10)) in 108 Korean GC tissues compared with the normal counterparts. A reduced level of methylation at the TEAD4 promoter was significantly associated with poor outcomes, including large tumor size, high-grade tumors and low survival rates. Compared with normal tissues, the TEAD4 protein was more frequently found in the nuclei of tumor cells along with YAP1 in 53 GC patients, demonstrating the posttranslational activation of this protein. Moreover, the knockdown of TEAD4 resulted in the reduced growth of GC cells both in vitro and in vivo. Finally, chromatin immunoprecipitation-sequencing and microarray analysis revealed the oncogenic properties of TEAD4 and its novel targets (ADM, ANG, ARID5B, CALD1, EDN2, FSCN1 and OSR2), which are involved in cell proliferation and migration. In conclusion, the multilevel perturbations of TEAD4 at epigenetic, transcriptional and posttranslational levels may contribute to GC development.

Variable allelic expression of imprinted genes in human pluripotent stem cells during differentiation into specialized cell types in vitro.

Genomic imprinting is an epigenetic phenomenon by which a subset of genes is asymmetrically expressed in a parent-of-origin manner. However, little is known regarding the epigenetic behaviors of imprinted genes during human development. Here, we show dynamic epigenetic changes in imprinted genes in hESCs during in vitro differentiation into specialized cell types. Out of 9 imprinted genes with single nucleotide polymorphisms, mono-allelic expression for three imprinted genes (H19, KCNQ1OT1, and IPW), and bi- or partial-allelic expression for three imprinted genes (OSBPL5, PPP1R9A, and RTL1) were stably retained in H9-hESCs throughout differentiation, representing imprinting stability. Three imprinted genes (KCNK9, ATP10A, and SLC22A3) showed a loss and a gain of imprinting in a lineage-specific manner during differentiation. Changes in allelic expression of imprinted genes were observed in another hESC line during in vitro differentiation. These findings indicate that the allelic expression of imprinted genes may be vulnerable in a lineage-specific manner in human pluripotent stem cells during differentiation.

Microarray screening and qRT-PCR evaluation of microRNA markers for forensic body fluid identification.

MicroRNAs (miRNA) are a class of small (∼22 nucleotides) noncoding RNAs that regulate diverse biological processes at the post-transcriptional level. MiRNAs have great potential for forensic body fluid identification because they are expressed in a tissue specific manner and are less prone to degradation. Previous studies reported several miRNAs as body fluid specific, but there are few overlaps among them. Here, we used a genome-wide miRNA microarray containing over 1700 miRNAs to assay 20 body fluid samples and identify novel miRNAs useful for forensic body fluid identification. Based on Shannon Entropy and Q-statistics, 203 miRNAs specifically expressed in each body fluid were first selected. Eight miRNAs were then selected as novel forensically relevant miRNA markers: miR-484 and miR-182 for blood, miR-223 and miR-145 for saliva, miR-2392 and miR-3197 for semen, and miR-1260b and miR-654-5p for vaginal secretions. When the eight selected miRNAs were evaluated in 40 additional body fluid samples by qRT-PCR, they showed high sensitivity and specificity for the identification of the target body fluid. We suggest that the eight miRNAs may be candidates for developing an effective molecular assay for forensic body fluid identification.

Whole Genome Sequencing Analysis of Soybean Diversity Across Different Countries and Selection Signature of Korean Soybean Accession.

Soybean is an important agricultural crop known for its high protein and oil content, contributing to essential nutritional and health benefits for humans. Domesticated in China over 5000 years ago, soybean has since adapted to diverse environments and spread worldwide. This study aimed to investigate the genomic characteristics and population structures of 2,317 publicly available soybean whole-genome sequences from diverse geographical regions, including China, Korea, Japan, Europe, North America and South America. We used large-scale whole-genome sequencing data to perform high-resolution analyses to reveal the genetic characteristics of soybean accessions. Soybean accessions from China and Korea exhibited landrace characteristics, indicating higher genetic diversity and adaptation to local environments. On the other hand, soybean accessions from Japan, the European Union, and South America were found to have low genetic diversity due to artificial selection and breeding for agronomic traits. We also identified key variants and genes associated with the ability to adapt to different environments. In Korean soybean accessions, we observed strong selection signals for isoflavone synthesis, an adaptive trait critical for improving soybean adaptability, survival, and reproductive success by mitigating environmental stress. Identifying specific genomic regions showing unique patterns of selective sweeps for genes such as HIDH, CYP73A11, IFS1, and CYP81E11 associated with isoflavone synthesis provided valuable insights into potential adaptation mechanisms. Our research has significantly improved our understanding of soybean diversity at the genetic level. We have identified key genetic variants and genes influencing adaptability, laying the foundation for future advances in genomics-based breeding programs and crop improvement efforts.

Proteogenomic Characterization Reveals Estrogen Signaling as a Target for Never-Smoker Lung Adenocarcinoma Patients without EGFR or ALK Alterations.

Never-smoker lung adenocarcinoma (NSLA) is prevalent in Asian populations, particularly in women. EGFR mutations and anaplastic lymphoma kinase (ALK) fusions are major genetic alterations observed in NSLA, and NSLA with these alterations have been well studied and can be treated with targeted therapies. To provide insights into the molecular profile of NSLA without EGFR and ALK alterations (NENA), we selected 141 NSLA tissues and performed proteogenomic characterization, including whole genome sequencing (WGS), transcriptomic, methylation EPIC array, total proteomic, and phosphoproteomic analyses. Forty patients with NSLA harboring EGFR and ALK alterations and seven patients with NENA with microsatellite instability were excluded. Genome analysis revealed that TP53 (25%), KRAS (22%), and SETD2 (11%) mutations and ROS1 fusions (14%) were the most frequent genetic alterations in NENA patients. Proteogenomic impact analysis revealed that STK11 and ERBB2 somatic mutations had broad effects on cancer-associated genes in NENA. DNA copy number alteration analysis identified 22 prognostic proteins that influenced transcriptomic and proteomic changes. Gene set enrichment analysis revealed estrogen signaling as the key pathway activated in NENA. Increased estrogen signaling was associated with proteogenomic alterations, such as copy number deletions in chromosomes 14 and 21, STK11 mutation, and DNA hypomethylation of LLGL2 and ST14. Finally, saracatinib, an Src inhibitor, was identified as a potential drug for targeting activated estrogen signaling in NENA and was experimentally validated in vitro. Collectively, this study enhanced our understanding of NENA NSLA by elucidating the proteogenomic landscape and proposed saracatinib as a potential treatment for this patient population that lacks effective targeted therapies. SIGNIFICANCE: The proteogenomic landscape in never-smoker lung cancer without known driver mutations reveals prognostic proteins and enhanced estrogen signaling that can be targeted as a potential therapeutic strategy to improve patient outcomes.

CYB5R3 functions as a tumor suppressor by inducing ER stress-mediated apoptosis in lung cancer cells via the PERK-ATF4 and IRE1α-JNK pathways.

Cytochrome b5 reductase 3 (CYB5R3) is involved in various cellular metabolic processes, including fatty acid synthesis and drug metabolism. However, the role of CYB5R3 in cancer development remains poorly understood. Here, we show that CYB5R3 expression is downregulated in human lung cancer cell lines and tissues. Adenoviral overexpression of CYB5R3 suppresses lung cancer cell growth in vitro and in vivo. However, CYB5R3 deficiency promotes tumorigenesis and metastasis in mouse models. Transcriptome analysis revealed that apoptosis- and endoplasmic reticulum (ER) stress-related genes are upregulated in CYB5R3-overexpressing lung cancer cells. Metabolomic analysis revealed that CYB5R3 overexpression increased the production of nicotinamide adenine dinucleotide (NAD+) and oxidized glutathione (GSSG). Ectopic CYB5R3 is mainly localized in the ER, where CYB5R3-dependent ER stress signaling is induced via activation of protein kinase RNA-like ER kinase (PERK) and inositol-requiring enzyme 1 alpha (IRE1α). Moreover, NAD+ activates poly (ADP-ribose) polymerase16 (PARP16), an ER-resident protein, to promote ADP-ribosylation of PERK and IRE1α and induce ER stress. In addition, CYB5R3 induces the generation of reactive oxygen species and caspase-9-dependent intrinsic cell death. Our findings highlight the importance of CYB5R3 as a tumor suppressor for the development of CYB5R3-based therapeutics for lung cancer.

Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses.

The role of the serine/glycine metabolic pathway (SGP) has recently been demonstrated in tumors; however, the pathological relevance of the SGP in thyroid cancer remains unexplored. Here, we perform metabolomic profiling of 17 tumor-normal pairs; bulk transcriptomics of 263 normal thyroid, 348 papillary, and 21 undifferentiated thyroid cancer samples; and single-cell transcriptomes from 15 cases, showing the impact of mitochondrial one-carbon metabolism in thyroid tumors. High expression of serine hydroxymethyltransferase-2 (SHMT2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is associated with low thyroid differentiation scores and poor clinical features. A subpopulation of tumor cells with high mitochondrial one-carbon pathway activity is observed in the single-cell dataset. SHMT2 inhibition significantly compromises mitochondrial respiration and decreases cell proliferation and tumor size in vitro and in vivo. Collectively, our results highlight the importance of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer and suggest that SHMT2 is a potent therapeutic target.

Aberrant upregulation of ASCL2 by promoter demethylation promotes the growth and resistance to 5-fluorouracil of gastric cancer cells.

Achaete scute-like 2 (ASCL2), a basic helix-loop-helix transcription factor, plays an essential role in the maintenance of adult intestinal stem cells. However, the function of ASCL2 in gastric cancer (GC) is poorly understood. Therefore, we investigated the roles and regulatory transcription mechanisms of ASCL2 in GC. Gene expression and methylation data analysis showed that ASCL2 was upregulated and hypomethylated in GC tissues. Using real-time RT-PCR and pyrosequencing analysis, we confirmed that ASCL2 was overexpressed and hypomethylated in GC tissues compared to adjacent normal tissues. We then investigated the mechanisms underlying the aberrant expression of ASCL2 in GC and found that treatment with a methylation inhibitor induced ASCL2 expression in GC cell lines. MBD-sequencing assay also revealed hypermethylation of the promoter region of ASCL2 in GC cell lines, which barely expressed the ASCL2 gene. Furthermore, ASCL2 expression levels were inversely correlated with GC patient survival. Ectopic overexpression of ASCL2 showed that ASCL2 increased cell growth and promoted resistance to 5-fluorouracil in GC cells. These results suggest that ASCL2 might play an important role in gastric tumor growth and chemoresistance, and could be a useful prognostic marker for GC patients.

Generation and analysis of whole-genome sequencing data in human mammary epithelial cells.

Breast cancer is the most common cancer worldwide, and advanced breast cancer with metastases is incurable mainly with currently available therapies. Therefore, it is essential to understand molecular characteristics during the progression of breast carcinogenesis. Here, we report a dataset of whole genomes from the human mammary epithelial cell system derived from a reduction mammoplasty specimen. This system comprises pre-stasis 184D cells, considered normal, and seven cell lines along cancer progression series that are immortalized or additionally acquired anchorage-independent growth. Our analysis of the whole-genome sequencing (WGS) data indicates that those seven cancer progression series cells have somatic mutations whose number ranges from 8,393 to 39,564 (with an average of 30,591) compared to 184D cells. These WGS data and our mutation analysis will provide helpful information to identify driver mutations and elucidate molecular mechanisms for breast carcinogenesis.

DNA methylome profiling of blood to identify individuals in a pair of monozygotic twins.

BACKGROUND: Short tandem repeat (STR) markers cannot be used to distinguish between genetically identical monozygotic (MZ) twins, causing problems in a case with an MZ twin as a suspect. Many studies have shown that in older MZ twins, there are significant differences in overall content and genomic distribution of methylation. OBJECTIVE: In this study, we analyzed the DNA methylome profile of blood to identify recurrent differentially methylated CpG sites (DMCs) to discriminate between MZ twins. METHODS: Blood samples were collected from 47 paired MZ twins. We performed the DNA methylation profiling using the HumanMethylation EPIC BeadChip platform and identified recurrent DMCs between MZ twins. Then, Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and motif enrichment analyses were performed to reveal the biological functions of recurrent DMCs. We collected DNA methylome data from the Gene Expression Omnibus (GEO) public database to verify the recurrent DMCs between MZ twins. RESULTS: We identified recurrent DMCs between MZ twin samples and observed that they were enriched in immune-related genes. In addition, we verified our DMCs in a public dataset. CONCLUSION: Our results suggest that the methylation level at recurrent DMCs between MZ twins may serve as a valuable biomarker for identification of individuals in a pair of MZ twins.

Identification of NIFTP-Specific mRNA Markers for Reliable Molecular Diagnosis of Thyroid Tumors.

Non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP) is a low-risk thyroid tumor with a favorable prognosis. Nonetheless, differentiating NIFTP from other thyroid tumors remains challenging, necessitating reliable diagnostic markers. This study is aimed at discovering NIFTP-specific mRNA markers through RNA sequencing analysis of thyroid tumor tissues. We performed mRNA expression profiling for 74 fresh frozen thyroid tissue samples, including NIFTP and benign and malignant follicular-cell-derived tumors. NIFTP/malignant tumors showed 255 downregulated genes and 737 upregulated genes compared to benign tumors. Venn diagram analysis revealed 19 significantly upregulated and 7 downregulated mRNAs in NIFTP. Akaike information criterion analysis allowed us to select OCLN, ZNF423, LYG1, and AQP5 mRNA markers. We subsequently developed a predictive model based on logistic regression analysis using these four mRNAs, which we validated in independent samples (n = 90) using a qRT-PCR assay. This model demonstrated high accuracy in predicting NIFTP in discovery dataset (AUC (area under the receiver operating characteristic) = 0.960) and the validation dataset (AUC = 0.757). Our results suggest that OCLN, ZNF423, LYG1, and AQP5 mRNA markers might serve as reliable molecular markers for identifying NIFTP among other thyroid tumors, ultimately aiding in accurate diagnosis and management of NIFTP patients.