Genomic and Pathologic Profiling of Very Well-Differentiated Gastric Adenocarcinoma of Intestinal Type: A Study With Emphasis on Diffuse-Type Transformation.

Very well-differentiated adenocarcinoma of intestinal type is a distinct subtype of gastric cancer characterized by anastomosing glands with a hand-in-hand pattern and low-grade cytologic atypia resembling intestinal metaplasia. This is a slow-growing neoplasm with an indolent clinical course; however, a subset demonstrates transformation into adenocarcinoma with higher-grade histology, typically diffuse-type carcinoma, and behaves aggressively. This study aimed to better characterize the genomic and pathologic features, with a focus on factors associated with diffuse-type transformation. A total of 58 cases with (n=31) and without (n=27) diffuse-type transformation were analyzed for molecular and pathologic features. First, comprehensive deep DNA sequencing was conducted in 18 cases (discovery cohort), followed by a digital droplet polymerase chain reaction of hot spot RHOA mutations in 40 cases (validation cohort). In total, RHOA mutations were the most common alteration (34%), followed by loss of ARID1A (12%), p53 alterations (10%), and CLDN18 :: ARHGAP26/6 fusions (3.4%). FGFR2 amplification was identified in an advanced case with a p53 alteration. Altered p53 expression was recognized only in higher-grade components and was significantly associated with advanced disease ( P =0.0015) and diffuse-type transformation ( P =0.026). A mixed mucin phenotype was also strongly correlated with advanced disease ( P <0.001) and diffuse-type transformation ( P <0.001). Decreased E-cadherin expression was frequently observed (74%) in poorly cohesive components. This study demonstrated that a subset of RHOA -mutant diffuse-type gastric cancers develops through the transformation of very well-differentiated adenocarcinoma of intestinal type. Our observations suggest a mixed mucin phenotype as a risk factor and alterations in p53 and E-cadherin as drivers of diffuse-type transformation.

Genomic and epigenomic integrative subtypes of renal cell carcinoma in a Japanese cohort.

Renal cell carcinoma (RCC) comprises several histological types characterised by different genomic and epigenomic aberrations; however, the molecular pathogenesis of each type still requires further exploration. We perform whole-genome sequencing of 128 Japanese RCC cases of different histology to elucidate the significant somatic alterations and mutagenesis processes. We also perform transcriptomic and epigenomic sequencing to identify distinguishing features, including assay for transposase-accessible chromatin sequencing (ATAC-seq) and methyl sequencing. Genomic analysis reveals that the mutational signature differs among the histological types, suggesting that different carcinogenic factors drive each histology. From the ATAC-seq results, master transcription factors are identified for each histology. Furthermore, clear cell RCC is classified into three epi-subtypes, one of which expresses highly immune checkpoint molecules with frequent loss of chromosome 14q. These genomic and epigenomic features may lead to the development of effective therapeutic strategies for RCC.

Sleeping Beauty transposon mutagenesis identified genes and pathways involved in inflammation-associated colon tumor development.

Chronic inflammation promotes development and progression of colorectal cancer (CRC). To comprehensively understand the molecular mechanisms underlying the development and progression of inflamed CRC, we perform in vivo screening and identify 142 genes that are frequently mutated in inflammation-associated colon tumors. These genes include senescence and TGFß-activin signaling genes. We find that TNFα can induce stemness and activate senescence signaling by enhancing cell plasticity in colonic epithelial cells, which could act as a selective pressure to mutate senescence-related genes in inflammation-associated colonic tumors. Furthermore, we show the efficacy of the Cdk4/6 inhibitor in vivo for inflammation-associated colonic tumors. Finally, we functionally validate that Arhgap5 and Mecom are tumor suppressor genes, providing possible therapeutic targets for CRC. Thus, we demonstrate the importance of the inactivation of senescence pathways in CRC development and progression in an inflammatory microenvironment, which can help progress toward precision medicine.

Oncogenic structural aberration landscape in gastric cancer genomes.

Structural variants (SVs) are responsible for driver events in gastric cancer (GC); however, their patterns and processes remain poorly understood. Here, we examine 170 GC whole genomes to unravel the oncogenic structural aberration landscape in GC genomes and identify six rearrangement signatures (RSs). Non-random combinations of RSs elucidate distinctive GC subtypes comprising one or a few dominant RS that are associated with specific driver events (BRCA1/2 defects, mismatch repair deficiency, and TP53 mutation) and epidemiological backgrounds. Twenty-seven SV hotspots are identified as GC driver candidates. SV hotspots frequently constitute complexly clustered SVs involved in driver gene amplification, such as ERBB2, CCNE1, and FGFR2. Further deconstruction of the locally clustered SVs uncovers amplicon-generating profiles characterized by super-large SVs and intensive segmental amplifications, contributing to the extensive amplification of GC oncogenes. Comprehensive analyses using adjusted SV allele frequencies indicate the significant involvement of extra-chromosomal DNA in processes linked to specific RSs.

Multiancestry genomic and transcriptomic analysis of gastric cancer.

Gastric cancer is among the most common malignancies worldwide, characterized by geographical, epidemiological and histological heterogeneity. Here, we report an extensive, multiancestral landscape of driver events in gastric cancer, involving 1,335 cases. Seventy-seven significantly mutated genes (SMGs) were identified, including ARHGAP5 and TRIM49C. We also identified subtype-specific drivers, including PIGR and SOX9, which were enriched in the diffuse subtype of the disease. SMGs also varied according to Epstein-Barr virus infection status and ancestry. Non-protein-truncating CDH1 mutations, which are characterized by in-frame splicing alterations, targeted localized extracellular domains and uniquely occurred in sporadic diffuse-type cases. In patients with gastric cancer with East Asian ancestry, our data suggested a link between alcohol consumption or metabolism and the development of RHOA mutations. Moreover, mutations with potential roles in immune evasion were identified. Overall, these data provide comprehensive insights into the molecular landscape of gastric cancer across various subtypes and ancestries.

Genomic landscape of chemical-induced lung tumors under Nrf2 different expression levels.

The transcription factor Nrf2 plays a crucial role in the anti-oxidative stress response, protection of DNA from injury and DNA repair mechanisms. Nrf2 activity reduces cancer initiation, but how Nrf2 affects whole-genome alterations upon carcinogenic stimulus remains unexplored. Although recent genome-wide analysis using next-generation sequencing revealed landscapes of nucleotide mutations and copy number alterations in various human cancers, genomic changes in murine cancer models have not been thoroughly examined. We elucidated the relationship between Nrf2 expression levels and whole exon mutation patterns using an ethyl-carbamate (urethane)-induced lung carcinogenesis model employing Nrf2-deficient and Keap1-kd mice, the latter of which express high levels of Nrf2. Exome analysis demonstrated that single nucleotide and trinucleotide mutation patterns and the Kras mutational signature differed significantly and were dependent on the expression level of Nrf2. The Nrf2-deficient tumors exhibited fewer copy number alterations relative to the Nrf2-wt and Keap1-kd tumors. The observed trend in genomic alterations likely prevented the Nrf2-deficient tumors from progressing into malignancy. For the first time, we present whole-exome sequencing results for chemically-induced lung tumors in the Nrf2 gain or loss of function mouse models. Our results demonstrate that different Nrf2 expression levels lead to distinct gene mutation patterns that underly different oncogenic mechanisms in each tumor genotype.

Diffusely infiltrating glioma with CREBBP-BCORL1 fusion showing overexpression of not only BCORL1 but BCOR: A case report.

BCORL1 encodes a transcriptional corepressor homolog to BCOR. BCORL1 rearrangements have been previously described as rare events, and among them, CREBBP-BCORL1 has been reported only in 2 cases of ossifying fibromyxoid tumors. Herein, we present the first case of diffusely infiltrating glioma with CREBBP-BCORL1 involving a 17-year-old female patient. Histologically, the tumor was composed of a diffusely infiltrative proliferation of small tumor cells with moderate cellularity showing prominent microcystic formation. DNA methylation analysis revealed that the current case and a previously reported anaplastic ependymoma with EP300-BCORL1 were clustered together in close proximity to but distinct from methylation class high-grade neuroepithelial tumor with BCOR alteration. RNA sequencing demonstrated high mRNA expression of not only BCORL1 but BCOR, and the latter was compatible with diffuse nuclear expression of BCOR detected by immunohistochemistry. Our findings suggest that central nervous system tumors with CREBBP/EP300-BCORL1 may exhibit diverse morphologies but form a distinct DNA methylation group and that BCORL1 fusion genes may lead to upregulation of both BCOR and BCORL1.

Transcriptome and methylome analysis of CNS germ cell tumor finds its cell-of-origin in embryogenesis and reveals shared similarities with testicular counterparts.

BACKGROUND: CNS germ cell tumors (GCTs) predominantly develop in pediatric and young adult patients with variable responses to surgery, radiation, and chemotherapy. This study aimed to examine the complex and largely unknown pathogenesis of CNS GCTs. METHODS: We used a combined transcriptomic and methylomic approach in 84 cases and conducted an integrative analysis of the normal cells undergoing embryogenesis and testicular GCTs. RESULTS: Genome-wide transcriptome analysis in CNS GCTs indicated that germinoma had a transcriptomic profile representative of primitive cells during early embryogenesis with high meiosis/mitosis potentials, while nongerminomatous GCTs (NGGCTs) had differentiated phenotypes oriented toward tissue formation and organogenesis. Co-analysis with the transcriptome of human embryonic cells revealed that germinomas had expression profiles similar to those of primordial germ cells, while the expression profiles of NGGCTs were similar to those of embryonic stem cells. Some germinoma cases were characterized by extensive immune-cell infiltration and high expression of cancer-testis antigens. NGGCTs had significantly higher immune-cell infiltration, characterized by immune-suppression phenotype. CNS and testicular GCTs (TGCTs) had similar mutational profiles; TGCTs showed enhanced copy number alterations. Methylation analysis clustered germinoma/seminoma and nongerminoma/nonseminoma separately. Germinoma and seminoma were co-categorized based on the degree of the tumor microenvironment balance. CONCLUSIONS: These results suggested that the pathophysiology of GCTs was less dependent on their site of origin and more dependent on the state of differentiation as well as on the tumor microenvironment balance. This study revealed distinct biological properties of GCTs, which will hopefully lead to future treatment development.

Identification of novel SSX1 fusions in synovial sarcoma.

Synovial sarcoma is characterized by variable epithelial differentiation and specific SS18-SSX gene fusions. The diagnosis is primarily based on phenotype, but fusion gene detection is increasingly being considered indispensable, with SS18 break-apart fluorescence in situ hybridization (FISH) being favored in many laboratories. However, SS18 FISH assay produces negative or atypical results in a minority of cases, leaving uncertainties in diagnosis and management. Here, we analyzed this challenging subset of SS18 FISH-negative/atypical synovial sarcoma using RNA sequencing and monoclonal antibodies that recognize SS18-SSX and the SSX C-terminus. Among 99 synovial sarcoma cases that were previously subjected to SS18 break-apart FISH, eight cases were reported as negative and three cases were indeterminate, owing to atypical signal patterns. Three of these 11 tumors (two monophasic and one biphasic) harbored novel EWSR1-SSX1 fusions, were negative for SS18-SSX staining, and were positive for SSX C-terminus staining. One monophasic tumor harbored a novel MN1-SSX1 fusion, and showed negative SS18-SSX expression and positive SSX C-terminus staining. Another monophasic tumor carried an SS18L1-SSX1 fusion, and was weakly positive for SS18-SSX, while SMARCB1 expression was reduced. The presence of these novel and/or rare fusions was confirmed using RT-PCR and Sanger sequencing. EWSR1-SSX1 was further validated by EWSR1 FISH assay. The remaining six tumors (five monophasic and one biphasic) showed strong SS18-SSX expression, and RNA sequencing successfully performed in three cases identified canonical SS18-SSX2 fusions. Based on a DNA methylation-based unsupervised clustering, the tumors with EWSR1-SSX1 and SS18L1-SSX1 clustered with synovial sarcoma, while the MN1-SSX1-positive tumor was not co-clustered despite classic histology and immunoprofile. In summary, we discovered novel and rare SSX1 fusions to non-SS18 genes in synovial sarcoma. The expanded genetic landscape carries significant diagnostic implications and advances our understanding of the oncogenic mechanism.

Comprehensive Genomic Profiling of Neuroendocrine Carcinomas of the Gastrointestinal System.

The neuroendocrine carcinoma of the gastrointestinal system (GIS-NEC) is a rare but highly malignant neoplasm. We analyzed 115 cases using whole-genome/exome sequencing, transcriptome sequencing, DNA methylation assays, and/or ATAC-seq and found GIS-NECs to be genetically distinct from neuroendocrine tumors (GIS-NET) in the same location. Clear genomic differences were also evident between pancreatic NECs (Panc-NEC) and nonpancreatic GIS-NECs (Nonpanc-NEC). Panc-NECs could be classified into two subgroups (i.e., "ductal-type" and "acinar-type") based on genomic features. Alterations in TP53 and RB1 proved common in GIS-NECs, and most Nonpanc-NECs with intact RB1 demonstrated mutually exclusive amplification of CCNE1 or MYC. Alterations of the Notch gene family were characteristic of Nonpanc-NECs. Transcription factors for neuroendocrine differentiation, especially the SOX2 gene, appeared overexpressed in most GIS-NECs due to hypermethylation of the promoter region. This first comprehensive study of genomic alterations in GIS-NECs uncovered several key biological processes underlying genesis of this very lethal form of cancer. SIGNIFICANCE: GIS-NECs are genetically distinct from GIS-NETs. GIS-NECs arising in different organs show similar histopathologic features and share some genomic features, but considerable differences exist between Panc-NECs and Nonpanc-NECs. In addition, Panc-NECs could be classified into two subgroups (i.e., "ductal-type" and "acinar-type") based on genomic and epigenomic features. This article is highlighted in the In This Issue feature, p. 587.

Ependymoma-like tumor with mesenchymal differentiation harboring C11orf95-NCOA1/2 or -RELA fusion: A hitherto unclassified tumor related to ependymoma.

Recurrent fusion genes involving C11orf95, C11orf95-RELA, have been identified only in supratentorial ependymomas among primary CNS tumors. Here, we report hitherto histopathologically unclassifiable high-grade tumors, under the tentative label of "ependymoma-like tumors with mesenchymal differentiation (ELTMDs)," harboring C11orf95-NCOA1/2 or -RELA fusion. We examined the clinicopathological and molecular features in five cases of ELTMDs. Except for one adult case (50 years old), all cases were in children ranging from 1 to 2.5 years old. All patients presented with a mass lesion in the cerebral hemisphere. Histologically, all cases demonstrated a similar histology with a mixture of components. The major components were embryonal-appearing components forming well-delineated tumor cell nests composed of small uniform cells with high proliferative activity, and spindle-cell mesenchymal components with a low- to high-grade sarcoma-like appearance. The embryonal-appearing components exhibited minimal ependymal differentiation including a characteristic EMA positivity and tubular structures, but histologically did not fit with ependymoma because they lacked perivascular pseudorosettes, a histological hallmark of ependymoma, formed well-delineated nests, and had diffuse and strong staining for CAM5.2. Molecular analysis identified C11orf95-NCOA1, -NCOA2, and -RELA in two, one, and two cases, respectively. t-distributed stochastic neighbor embedding analysis of DNA methylation data from two cases with C11orf95-NCOA1 or -NCOA2 and a reference set of 380 CNS tumors revealed that these two cases were clustered together and were distinct from all subgroups of ependymomas. In conclusion, although ELTMDs exhibited morphological and genetic associations with supratentorial ependymoma with C11orf95-RELA, they cannot be regarded as ependymoma. Further analyses of more cases are needed to clarify their differences and similarities.

Molecular detection and clinicopathological characteristics of advanced/recurrent biliary tract carcinomas harboring the FGFR2 rearrangements: a prospective observational study (PRELUDE Study).

BACKGROUND: Fibroblast growth factor receptor 2 (FGFR2) rearrangement is expected to be a novel therapeutic target in advanced/recurrent biliary tract cancer (BTC). However, efficient detection and the exact frequency of FGFR2 rearrangements among patients with advanced/recurrent BTC have not been determined, and the clinical characteristics of FGFR2 rearrangement-positive patients have not been fully elucidated. We aimed to determine the frequency of FGFR2 rearrangement-positive patients among those with advanced/recurrent BTC and elucidate their clinicopathological characteristics. METHODS: Paraffin-embedded tumor samples from formalin-fixed surgical or biopsy specimens of patients with advanced/recurrent BTC were analyzed for positivity of FGFR2 rearrangement by fluorescent in situ hybridization (FISH). RNA sequencing was performed on samples from all FISH-positive and part of FISH-negative patients. RESULTS: A total of 445 patients were enrolled. FISH was performed on 423 patients (272 patients with intrahepatic cholangiocarcinoma (ICC), 83 patients with perihilar cholangiocarcinoma (PCC), and 68 patients with other BTC). Twenty-one patients with ICC and four patients with PCC were diagnosed as FGFR2-FISH positive. Twenty-three of the 25 FISH-positive patients (20 ICC and 3 PCC) were recognized as FGFR2 rearrangement positive by targeted RNA sequencing. Younger age (≤ 65 years; p = 0.018) and HCV Ab- and/or HBs Ag-positivity (p = 0.037) were significantly associated with the presence of FGFR2 rearrangement (logistic regression). CONCLUSIONS: FGFR2 rearrangement was identified in ICC and PCC patients, and was associated with younger age and history of hepatitis viral infection.

Ependymoma with C11orf95-MAML2 fusion: presenting with granular cell and ganglion cell features.

C11orf95-RELA fusion or, less frequently, YAP1 fusion is recurrently detected in most cases of supratentorial ependymoma. Other fusions have rarely been reported in some cases of supratentorial ependymoma, and little is known about their pathological or clinical features. Here, we present a case of supratentorial ependymoma with unusual pathological findings and C11orf95-MAML2 fusion. A 23-year-old man was admitted to our hospital because of headache and vomiting. Magnetic resonance imaging revealed a cystic lesion in the right frontal lobe, and gross total resection of the tumor was performed. Pathologically, the tumor was mainly composed of typical ependymal lesions with perivascular pseudorosettes and contained some atypical lesions, with granular and ganglion cell features. The tumor was diagnosed as anaplastic ependymoma, which was classified as grade III on the World Health Organization scale, and found to be RELA fusion-positive in the DNA methylation analysis. However, the tumor was negative for C11orf95-RELA fusion, and RNA sequencing detected C11orf95-MAML2 fusion. The patient has not received adjuvant therapy and has remained alive without any evidence of disease for 30 months, suggesting that the prognosis might be better than that of typical C11orf95-RELA fusion-positive ependymoma.

Primary spinal intramedullary Ewing-like sarcoma harboring CIC-DUX4 translocation: a similar cytological appearance as its soft tissue counterpart but no lobulation in association with desmoplastic stroma.

The CIC-DUX4 translocation is the most common genetic alteration of small round cell sarcomas without EWSR1 rearrangement. These "Ewing-like sarcomas" usually occur in peripheral soft tissues, and rare primary central nervous system (CNS) tumors have been described. We report a rare case of primary spinal intramedullary Ewing-like sarcoma harboring CIC-DUX4 translocation. A 23-year-old man presented with weakness in the extremities. Magnetic resonance imaging revealed a large intramedullary tumor spanning C3-C5 with heterogeneous enhancement following gadolinium administration. Histologically, most of the tumor displayed dense myeloid proliferation composed of medium- to slightly small-sized primitive cells. Postoperatively, he received local adjuvant radiation therapy without tumor progression for 10 months. Target RNA sequencing analysis revealed the CIC-DUX4 fusion gene. Methylation array analysis resulted in a diagnosis of "methylation class CNS Ewing sarcoma family tumor with CIC alteration". Although this tumor lacked characteristic histological features such as lobular structures in association with desmoplastic stroma, relatively uniform nuclei with prominent nucleoli and eosinophilic cytoplasm, which are often found in CIC-rearranged sarcomas of soft tissue, were identified. Recently, many CNS and soft tissue tumors require genetic analysis for precise diagnosis. To consider certain molecular testing, careful histological examination is essential.

Expanding the clinicopathologic and molecular spectrum of BCOR-associated sarcomas in adults.

AIMS: BCOR gene alteration is a genetic signature of rare subsets of sarcomas. Most BCOR-associated sarcomas thus far reported are in the pediatric population, except for uterine sarcomas. We studied seven cases of BCOR-associated non-uterine sarcomas in adult patients. METHODS AND RESULTS: The patients were four men and three women ranging from 26 to 71 years in age. Three tumors, two of which primarily affected the kidney, showed BCOR-CCNB3. One tumor with a ZC3H7B-BCOR occurred in the chest wall, and a tumor with a novel CIITA-BCOR was found in the sinonasal tract. Two tumors in the lung and breast harbored exon 15 internal tandem duplications of BCOR, a highly unexpected observation in this age group. All seven sarcomas consisted of dense proliferations of uniform round to spindle cells with fine chromatin within vascular stroma. BCOR-CCNB3 sarcomas showed swirling fascicular growth. The tumor with the ZC3H7B-BCOR fusion showed a multinodular growth of spindle cells, and the tumors with the CIITA-BCOR fusion showed palisading of oval cells. Both tumors with BCOR internal tandem duplication demonstrated nested to palisading growth of round cells within sclerotic non-myxoid stroma. All seven sarcomas diffusely expressed BCOR and SATB2 immunohistochemically, with all three BCOR-CCNB3 sarcomas being immunopositive for CCNB3. BCOR alterations were confirmed by RNA sequencing, polymerase chain reaction, Sanger sequencing, and/or fluorescence in situ hybridization. CONCLUSIONS: This study expands the clinicopathologic and molecular spectrum of BCOR-associated sarcomas, and emphasizes the importance of being aware of this entity in the differential diagnosis of adult non-uterine sarcomas.

KMT2A (MLL) fusions in aggressive sarcomas in young adults.

AIM: To characterise unclassifiable sarcomas by use of a combined histological and molecular approach. METHODS AND RESULTS: Using RNA sequencing, we identified in-frame fusions involving KMT2A (MLL) in two cases. Case 1 was a 20-year-old woman with a deep soft tissue mass in the thigh. The tumour consisted of monomorphic round, epithelioid and spindle cells in a highly sclerotic background that were focally immunopositive for CD34, CD31, and ERG. Case 2 was a 30-year-old woman with a tumour that affected the femur and surrounding soft tissue. The tumour consisted of monomorphic round to spindle cells that were immunopositive for BCOR, Wilms tumour 1, and NKX2-2. Both tumours were aggressive and had metastasised to the lung; both patients died within a few years. RNA sequencing identified a YAP1 (exon 5)-KMT2A (exon 4) fusion in case 1 and a VIM (exon 4)-KMT2A (exon 2) fusion in case 2, both of which were confirmed by reverse transcription polymerase chain reaction, Sanger sequencing, and fluorescence in-situ hybridisation. The fusion protein structure was different from that in acute leukaemia, suggesting a novel oncogenic mechanism. CONCLUSIONS: KMT2A fusions account for a subset of aggressive unclassifiable sarcomas in young adults. Although it is presently unclear whether these sarcomas belong to a single group, the well-established role of KMT2A fusions as drivers of acute leukaemia and a recent publication regarding identification of YAP1-KMT2A in one unclassifiable sarcoma support the significance of these fusions. Further studies on additional cases are necessary to fully understand the clinicopathological and molecular aspects of KMT2A-rearranged sarcomas.

Integrated genetic and epigenetic analysis of myxofibrosarcoma.

Myxofibrosarcoma (MFS) is a common adult soft tissue sarcoma characterized by an infiltrative growth pattern and a high local recurrence rate. Here we report the genetic and epigenetic landscape of MFS based on the results of whole-exome sequencing (N = 41), RNA sequencing (N = 29), and methylation analysis (N = 41), using 41 MFSs as a discovery set, and subsequent targeted sequencing of 140 genes in the entire cohort of 99 MFSs and 17 MFSs' data from TCGA. Fourteen driver genes are identified, including potentially actionable therapeutic targets seen in 37% of cases. There are frequent alterations in p53 signaling (51%) and cell cycle checkpoint genes (43%). Other conceivably actionable driver genes including ATRX, JAK1, NF1, NTRK1, and novel oncogenic BRAF fusion gene are identified. Methylation patterns cluster into three subtypes associated with unique combinations of driver mutations, clinical outcomes, and immune cell compositions. Our results provide a valuable genomic resource to enable the design of precision medicine for MFS.

Epigenome mapping of human normal purified hepatocytes: personal epigenome variation and genome-epigenome correlation.

AIM: The aim of this study was to reveal the epigenome landscape of human normal hepatocytes. MATERIALS & METHODS: Cells purified from partial hepatectomy specimens of Japanese patients were subjected to whole-genome bisulfite sequencing using postbisulfite adaptor tagging, chromatin immunoprecipitation sequencing, RNA sequencing and whole-genome sequencing. RESULTS: CHG and CHH methylations were inversely associated with gene expression. Histone modification profiles of personal differentially methylated regions (pDMRs) differed considerably among samples. pDMRs were observed around the transcription start sites of genes whose expression is reportedly regulated by CpG methylation. pDMRs were frequently observed in the vicinity of single-nucleotide variations and insertions/deletions. CONCLUSION: Genetic variations may induce epigenetic variations, generating individual differences in the phenotypes of normal hepatocytes through variations in expression.