Masaki Fukuyo

Evolutionary maintenance of selfish homing endonuclease genes in the absence of horizontal transfer

Homing endonuclease genes are “selfish” mobile genetic elements whose endonuclease promotes the spread of its own gene by creating a break at a specific target site and using the host machinery to repair the break by copying and inserting the gene at this site. Horizontal transfer across the boundary of a species or population within which mating takes place has been thought to be necessary for their evolutionary persistence. This is based on the assumption that they will become fixed in a host population, where opportunities of homing will disappear, and become susceptible to degeneration. To test this hypothesis, we modeled behavior of a homing endonuclease gene that moves during meiosis through double-strand break repair. We mathematically explored conditions for persistence of the homing endonuclease gene and elucidated their parameter dependence as phase diagrams. We found that, if the cost of the pseudogene is lower than that of the homing endonuclease gene, the 2 forms can persist in a population through autonomous periodic oscillation. If the cost of the pseudogene is higher, 2 types of dynamics appear that enable evolutionary persistence: bistability dependent on initial frequency or fixation irrespective of initial frequency. The prediction of long persistence in the absence of horizontal transfer was confirmed by stochastic simulations in finite populations. The average time to extinction of the endonuclease gene was found to be thousands of meiotic generations or more based on realistic parameter values. These results provide a solid theoretical basis for an understanding of these and other extremely selfish elements.

Genetic and epigenetic aberrations occurring in colorectal tumors associated with serrated pathway

To clarify molecular alterations in serrated pathway of colorectal cancer (CRC), we performed epigenetic and genetic analyses in sessile serrated adenoma/polyps (SSA/P), traditional serrated adenomas (TSAs) and high‐methylation CRC. The methylation levels of six Group‐1 and 14 Group‐2 markers, established in our previous studies, were analyzed quantitatively using pyrosequencing. Subsequently, we performed targeted exon sequencing analyses of 126 candidate driver genes and examined molecular alterations that are associated with cancer development. SSA/P showed high methylation levels of both Group‐1 and Group‐2 markers, frequent BRAF mutation and occurrence in proximal colon, which were features of high‐methylation CRC. But TSA showed low‐methylation levels of Group‐1 markers, less frequent BRAF mutation and occurrence at distal colon. SSA/P, but not TSA, is thus considered to be precursor of high‐methylation CRC. High‐methylation CRC had even higher methylation levels of some genes, e.g., MLH1, than SSA/P, and significant frequency of somatic mutations in nonsynonymous mutations (p < 0.0001) and insertion/deletions (p = 0.002). MLH1‐methylated SSA/P showed lower methylation level of MLH1 compared with high‐methylation CRC, and rarely accompanied silencing of MLH1 expression. The mutation frequencies were not different between MLH1‐methylated and MLH1‐unmethylated SSA/P, suggesting that MLH1 methylation might be insufficient in SSA/P to acquire a hypermutation phenotype. Mutations of mismatch repair genes, e.g., MSH3 and MSH6, and genes in PI3K, WNT, TGF‐β and BMP signaling (but not in TP53 signaling) were significantly involved in high‐methylation CRC compared with adenoma, suggesting importance of abrogation of these genes in serrated pathway.

TET2 functions as a resistance factor against DNA methylation acquisition during Epstein-Barr virus infection

Extensive DNA methylation is observed in gastric cancer with Epstein-Barr virus (EBV) infection, and EBV infection is the cause to induce this extensive hypermethylaton phenotype in gastric epithelial cells. However, some 5′ regions of genes do not undergo de novo methylation, despite the induction of methylation in surrounding regions, suggesting the existence of a resistance factor against DNA methylation acquisition. We conducted an RNA-seq analysis of gastric epithelial cells with and without EBV infection and found that TET family genes, especially TET2, were repressed by EBV infection at both mRNA and protein levels. TET2 was found to be downregulated by EBV transcripts, e.g. BARF0 and LMP2A, and also by seven human miRNAs targeting TET2, e.g., miR-93 and miR-29a, which were upregulated by EBV infection, and transfection of which into gastric cells repressed TET2. Hydroxymethylation target genes by TET2 were detected by hydroxymethylated DNA immunoprecipitation sequencing (hMeDIP-seq) with and without TET2 overexpression, and overlapped significantly with methylation target genes in EBV-infected cells. When TET2 was knocked down by shRNA, EBV infection induced de novo methylation more severely, including even higher methylation in methylation-acquired promoters or de novo methylation acquisition in methylation-protected promoters, leading to gene repression. TET2 knockdown alone without EBV infection did not induce de novo DNA methylation. These data suggested that TET2 functions as a resistance factor against DNA methylation in gastric epithelial cells and repression of TET2 contributes to DNA methylation acquisition during EBV infection.

Epstein–Barr virus infection induces genome-wide de novo DNA methylation in non-neoplastic gastric epithelial cells

Epstein–Barr virus (EBV)‐positive gastric cancer (GC) shows a higher DNA methylation epigenotype. EBV infection can causally induce genome‐wide aberrant DNA methylation, as previously demonstrated by in vitro infection experiments in the low‐methylation GC cell line MKN7. However, whether EBV exerts DNA methylation remodelling properties in non‐neoplastic epithelial cells remains unclear. Here we performed post‐infection time‐series DNA methylation analyses using the immortalized normal gastric epithelial cell line GES1. Genome‐wide analysis using Illuminas Infinium 450 k BeadArray demonstrated global de novo DNA methylation from post‐infection day 17, which was completed by 28 days in a manner similar to that observed in MKN7 cells. De novo methylation of all types of GC‐specific methylation marker genes was observed, indicating that EBV infection is sufficient for gastric epithelial cells to acquire an EBV‐positive GC epigenotype. Pyrosequencing demonstrated that methylation of the viral genome preceded that of the host cellular genome, suggesting the existence of well‐ordered mechanisms that induce methylation. Spatiotemporal representation with differential models revealed dynamic alterations of DNA methylation in promoter regions, occurring from lower‐CpG peripheral regions and extending to higher‐CpG core regions. In summary, EBV infection exerted powerful pressure to induce global de novo DNA methylation in non‐neoplastic cells within a month in a spatiotemporally well‐ordered manner. Copyright

Combined Secretomics and Transcriptomics Revealed Cancer-Derived GDF15 is Involved in Diffuse-Type Gastric Cancer Progression and Fibroblast Activation

Gastric cancer is classified into two subtypes, diffuse and intestinal. The diffuse-type gastric cancer (DGC) has poorer prognosis, and the molecular pathology is not yet fully understood. The purpose of this study was to identify functional secreted molecules involved in DGC progression. We integrated the secretomics of six gastric cancer cell lines and gene expression analysis of gastric cancer tissues with publicly available microarray data. Hierarchical clustering revealed characteristic gene expression differences between diffuse- and intestinal-types. GDF15 was selected as a functional secreted molecule owing to high expression only in fetal tissues. Protein expression of GDF15 was higher in DGC cell lines and tissues. Serum levels of GDF15 were significant higher in DGC patients as compared with healthy individuals and chronic gastritis patients, and positively correlated with wall invasion and lymph node metastasis. In addition, the stimulation of GDF15 on NIH3T3 fibroblast enhanced proliferation and up-regulated expression of extracellular matrix genes, which were similar to TGF-β stimulation. These results indicate that GDF15 contributes to fibroblast activation. In conclusion, this study revealed that GDF15 may be a novel functional secreted molecule for DGC progression, possibly having important roles for cancer progression via the affecting fibroblast function, as well as TGF-β.

Histone modification alteration coordinated with acquisition of promoter DNA methylation during Epstein-Barr virus infection

Aberrant DNA hypermethylation is a major epigenetic mechanism to inactivate tumor suppressor genes in cancer. Epstein-Barr virus positive gastric cancer is the most frequently hypermethylated tumor among human malignancies. Herein, we performed comprehensive analysis of epigenomic alteration during EBV infection, by Infinium HumanMethylation 450K BeadChip for DNA methylation and ChIP-sequencing for histone modification alteration during EBV infection into gastric cancer cell line MKN7. Among 7,775 genes with increased DNA methylation in promoter regions, roughly half were “DNA methylation-sensitive” genes, which acquired DNA methylation in the whole promoter regions and thus were repressed. These included anti-oncogenic genes, e.g. CDKN2A. The other half were “DNA methylation-resistant” genes, where DNA methylation is acquired in the surrounding of promoter regions, but unmethylated status is protected in the vicinity of transcription start site. These genes thereby retained gene expression, and included DNA repair genes. Histone modification was altered dynamically and coordinately with DNA methylation alteration. DNA methylation-sensitive genes significantly correlated with loss of H3K27me3 pre-marks or decrease of active histone marks, H3K4me3 and H3K27ac. Apoptosis-related genes were significantly enriched in these epigenetically repressed genes. Gain of active histone marks significantly correlated with DNA methylation-resistant genes. Genes related to mitotic cell cycle and DNA repair were significantly enriched in these epigenetically activated genes. Our data show that orchestrated epigenetic alterations are important in gene regulation during EBV infection, and histone modification status in promoter regions significantly associated with acquisition of de novo DNA methylation or protection of unmethylated status at transcription start site.

A Novel Approach to Helicobacter pylori Pan-Genome Analysis for Identification of Genomic Islands

Genomes of a given bacterial species can show great variation in gene content and thus systematic analysis of the entire gene repertoire, termed the pan-genome, is important for understanding bacterial intra-species diversity, population genetics, and evolution. Here, we analyzed the pan-genome from 30 completely sequenced strains of the human gastric pathogen Helicobacter pylori belonging to various phylogeographic groups, focusing on 991 accessory (not fully conserved) orthologous groups (OGs). We developed a method to evaluate the mobility of genes within a genome, using the gene order in the syntenically conserved regions as a reference, and classified the 991 accessory OGs into five classes: Core, Stable, Intermediate, Mobile, and Unique. Phylogenetic networks based on the gene content of Core and Stable classes are highly congruent with that created from the concatenated alignment of fully conserved core genes, in contrast to those of Intermediate and Mobile classes, which show quite different topologies. By clustering the accessory OGs on the basis of phylogenetic pattern similarity and chromosomal proximity, we identified 60 co-occurring gene clusters (CGCs). In addition to known genomic islands, including cag pathogenicity island, bacteriophages, and integrating conjugative elements, we identified some novel ones. One island encodes TerY-phosphorylation triad, which includes the eukaryote-type protein kinase/phosphatase gene pair, and components of type VII secretion system. Another one contains a reverse-transcriptase homolog, which may be involved in the defense against phage infection through altruistic suicide. Many of the CGCs contained restriction-modification (RM) genes. Different RM systems sometimes occupied the same (orthologous) locus in the strains. We anticipate that our method will facilitate pan-genome studies in general and help identify novel genomic islands in various bacterial species.

TP53 mutation at early stage of colorectal cancer progression from two types of laterally spreading tumors

Although most sporadic colorectal cancers (CRC) are thought to develop from protruded adenomas through the adenoma–carcinoma sequence, some CRC develop through flat lesions, so‐called laterally spreading tumors (LST). We previously analyzed epigenetic aberrations in LST and found that LST are clearly classified into two molecular subtypes: intermediate‐methylation with KRAS mutation and low‐methylation with absence of oncogene mutation. Intermediate‐methylation LST were mostly granular type LST (LST‐G) and low‐methylation LST were mostly non‐granular LST (LST‐NG). In the present study, we conducted a targeted exon sequencing study including 38 candidate CRC driver genes to gain insight into how these genes modulate the development of LST. We identified a mean of 11.5 suspected nonpolymorphic variants per sample, including indels and non‐synonymous mutations, although there was no significant difference in the frequency of total mutations between LST‐G and LST‐NG. Genes associated with RTK/RAS signaling pathway were mutated more frequently in LST‐G than LST‐NG (P = 0.004), especially KRAS mutation occurring at 70% (30/43) of LST‐G but 26% (13/50) of LST‐NG (P < 0.0001). Both LST showed high frequency of APC mutation, even at adenoma stage, suggesting its involvement in the initiation stage of LST, as it is involved at early stage of colorectal carcinogenesis via adenoma‐carcinoma sequence. TP53 mutation was never observed in adenomas, but was specifically detected in cancer samples. TP53 mutation occurred during development of intramucosal cancer in LST‐NG, but during development of cancer with submucosal invasion in LST‐G. It is suggested that TP53 mutation occurs in the early stages of cancer development from adenoma in both LST‐G and LST‐NG, but is involved at an earlier stage in LST‐NG.

DNA methylation epigenotype and clinical features of NRAS-mutation(+) colorectal cancer

Sporadic colorectal cancer (CRC) is classified into several molecular subtypes. We previously established two groups of DNA methylation markers through genome‐wide DNA methylation analysis to classify CRC into distinct subgroups: high‐, intermediate‐, and low‐methylation epigenotypes (HME, IME, and LME, respectively). HME CRC, also called CpG island methylator phenotype (CIMP)‐high CRC, shows methylation of both Group 1 markers (CIMP markers) and Group 2 markers, while IME/CIMP‐low CRC shows methylation of Group 2, but not of Group 1 markers, and LME CRC shows no methylation of either Group 1 or Group 2 markers. While BRAF‐ and KRAS‐mutation(+) CRC strongly correlated with HME and IME, respectively, clinicopathological features of NRAS‐mutation(+) CRC, including association with DNA methylation, remain unclear. To characterize NRAS‐mutation(+) CRC, the methylation levels of 19 methylation marker genes (6 Group 1 and 13 Group 2) were analyzed in 61 NRAS‐mutation(+) and 144 NRAS‐mutation(−) CRC cases by pyrosequencing, and their correlation with clinicopathological features was investigated. Different from KRAS‐mutation(+) CRC, NRAS‐mutation(+) CRC significantly correlated with LME. NRAS‐mutation(+) CRC showed significantly better prognosis than KRAS‐mutation(+) CRC (P = 3 × 10−4). NRAS‐mutation(+) CRC preferentially occurred in elder patients (P = 0.02) and at the distal colon (P = 0.006), showed significantly less lymph vessel invasion (P = 0.002), and correlated with LME (P = 8 × 10−5). DNA methylation significantly accumulated at the proximal colon. NRAS‐mutation(+) CRC may constitute a different subgroup from KRAS‐mutation(+) CRC, showing significant correlation with LME, older age, distal colon, and relatively better prognosis.

Restriction-modification system with methyl-inhibited base excision and abasic-site cleavage activities

The restriction-modification systems use epigenetic modification to distinguish between self and nonself DNA. A modification enzyme transfers a methyl group to a base in a specific DNA sequence while its cognate restriction enzyme introduces breaks in DNA lacking this methyl group. So far, all the restriction enzymes hydrolyze phosphodiester bonds linking the monomer units of DNA. We recently reported that a restriction enzyme (R.PabI) of the PabI superfamily with half-pipe fold has DNA glycosylase activity that excises an adenine base in the recognition sequence (5′-GTAC). We now found a second activity in this enzyme: at the resulting apurinic/apyrimidinic (AP) (abasic) site (5′-GT#C, # = AP), its AP lyase activity generates an atypical strand break. Although the lyase activity is weak and lacks sequence specificity, its covalent DNA–R.PabI reaction intermediates can be trapped by NaBH4 reduction. The base excision is not coupled with the strand breakage and yet causes restriction because the restriction enzyme action can impair transformation ability of unmethylated DNA even in the absence of strand breaks in vitro. The base excision of R.PabI is inhibited by methylation of the target adenine base. These findings expand our understanding of genetic and epigenetic processes linking those in prokaryotes and eukaryotes.