DNA Hypermethylation of Fgf16 and Tbx22 Associated with Cleft Palate During Palatal Fusion
Total Page:16
File Type:pdf, Size:1020Kb
Original Article http://dx.doi.org/10.1590/1678-7757-2018-0649 DNA hypermethylation of Fgf16 and Tbx22 associated with cleft palate during palatal fusion Abstract Xuan SHU1 Objective: Cleft palate (CP) is a congenital birth defect caused by the failure of palatal fusion. Little is known about the potential role of DNA methylation Zejun DONG1 in the pathogenesis of CP. This study aimed to explore the potential role of Liuhanghang CHENG1 DNA methylation in the mechanism of CP. Methodology: We established an Shenyou SHU1 all-trans retinoic acid (ATRA)-induced CP model in C57BL/6J mice and used methylation-dependent restriction enzymes (MethylRAD, FspEI) combined with high-throughput sequencing (HiSeq X Ten) to compare genome-wide DNA methylation profiles of embryonic mouse palatal tissues, between embryos from ATRA-treated vs. untreated mice, at embryonic gestation day 14.5 (E14.5) (n=3 per group). To confirm differentially methylated levels of susceptible genes, real-time quantitative PCR (qPCR) was used to correlate expression of differentially methylated genes related to CP. Results: We identified 196 differentially methylated genes, including 17,298 differentially methylated CCGG sites between ATRA-treated vs. untreated embryonic mouse palatal tissues (P<0.05, log2FC>1). The CP-related genes Fgf16 (P=0.008, log2FC=1.13) and Tbx22 (P=0.011, log2FC=1.64,) were hypermethylated. Analysis of Fgf16 and Tbx22, using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), identified 3 GO terms and 1 KEGG pathway functionally related to palatal fusion. The qPCR showed that changes in expression level negatively correlated with methylation levels. Conclusions: Taken together, these results suggest that hypermethylation of Fgf16 and Tbx22 is associated with decreased gene expression, which might be responsible for developmental failure of palatal fusion, eventually resulting in the formation of CP. Keywords: DNA methylation. Cleft palate. Gene expression. Submitted: November 1, 2018 Modification: February 17, 2019 Accepted: March 12, 2019 Corresponding address: 1Second Affiliated Hospital of Shantou University Medical College, Cleft Lip and Palate Treatment Shenyou Shu Center, Shantou, Guangdong, China. Cleft Lip and Palate Treatment Center, Second Affiliated Hospital of Shantou University Medical College 69 Dongxia North Road, Jinping District, Shantou 515041 - China. Phone: +86-18023235288 - Fax: +86-0754-83141156 e-mail: [email protected] J Appl Oral Sci. 1/8 2019;27:e20180649 DNA hypermethylation of Fgf16 and Tbx22 associated with cleft palate during palatal fusion Introduction and humans. We integrated DNA methylation analysis of Fgf16 and Tbx22, by MethylRAD-seq, with Cleft palate (CP) is a congenital birth defect their biological characteristics identified by GO and caused by both environmental and genetic factors.1 KEGG pathway enrichment analysis. The restriction It is universally acknowledged that palatal fusion is enzyme FspEI recognizes 5-methylcytosine and the most crucial process in the palate formation. In 5-hydroxymethylcytosine in CCGG and CCWGG (where 10 a mouse, palatal shelves appose at the midline and W=A or T), and cleaves DNA bilaterally to generate palatal fusion occurs at E 14.5, and an imbalance of 32 bp fragments at a methylated CCGG site and 31 embryonic palatal mesenchyme cell proliferation and bp fragments at a methylated CCWGG site, with either apoptosis might result in CP.2 Previous studies have FspE1 recognition site in the middle. The qPCR was demonstrated that Fgf16 and Tbx22 participate in used to correlate Fgf16 and Tbx22 expression levels murine palate development.3 However, the role of with methylation levels to elucidate the molecular Fgf16 and Tbx22 during palatal fusion is still unknown. regulatory mechanisms underlying the development DNA methylation is an important epigenetic of CP. modification and plays a crucial role in many biological processes, such as embryogenesis, cellular differentiation, X-chromosome inactivation, genomic Methodology imprinting and transcriptional silencing.4 Methylation patterns of specific genes have been found to undergo Ethics, animals and treatment dynamic changes in embryonic development and This study was approved by the Laboratory Animal contribute to tissue-specific gene expression.5 The DNA Ethical Committee of the Medical College of Shantou methylation pattern of a mouse undergoes dynamic University (SUMC2015-106; Guangdong, China). and widespread alterations during palatogenesis, C57BL/6J mice of 20–28 g in body weight and 8 to and failing to establish correct methylation patterns 10 weeks of age were purchased from the Beijing can result in CP,6 which suggests that CP-susceptible Vital River Laboratory Animal Technology Co. Ltd. genes (e.g. Fgf16 and Tbx22) in embryonic mice (Beijing, China). In this study, female mice were mated provide new clues to epigenetic markers involved in with male mice of similar weight and age overnight. CP. Nevertheless, details on the methylation patterns Embryonic gestation day 0.5 (E0.5) was designated of CP-susceptible genes during palatal fusion are to be at 8 AM the next day when a vaginal plug was very limited, and the methylation pattern of Fgf16 observed. Pregnant mice at E10.5 were randomly and Tbx22 underlying palate development and its divided into an all-trans retinoic acid (ATRA, Sigma- contribution to CP is still unclear. Aldrich, St. Louis, MO, USA)-treated and control To explore the potential involvement of DNA (mock-treated) group. Mice in the ATRA-treated group methylation in regulating palatal fusion, we previously were treated, via oral gavage, with ATRA at 70 mg/ established a CP model in which pregnant C57BL/6J kg dissolved in corn oil as described previously.7 The mice are treated with all-trans retinoic acid (ATRA) untreated group was given an equivalent volume of to introduce CP in the embryos.7 ATRA is a vitamin A corn oil. At E14.5, mice were euthanized, and the metabolite and functions to support normal pattern embryonic palatal shelves (3 ATRA-treated samples formation during embryogenesis. Abnormally high vs. 3 untreated samples) were resected and stored concentrations of ATRA may affect palatogenesis at -80°C until use. by interfering in medial edge epithelia (MEE) cell differentiation and apoptosis, including inhibition of DNA preparation, library construction and mesenchymal proliferation and signaling growth factors MethylRAD-seq (transforming growth factor-β and platelet-derived Genomic DNA was extracted from palatal shelf growth factor).8 Cuervo, et al.9 (2002) reported that tissues using the conventional cetyltrimethyl high ATRA concentration blocked palatal shelf fusion ammonium bromide (CTAB) method following the and increased apoptosis within the MEE cell adhesion manufacturer’s instructions (AMRESCO Inc, Solon, process, and which may result in fetal malformations, OH, USA), and MethylRAD library construction and including cleft palate, in both experimental animals sequencing were the same as previously reported.10 J Appl Oral Sci. 2/8 2019;27:e20180649 SHU X, DONG Z, CHENG L, SHU S Paired-end sequencing was performed on a HiSeq X Protein-protein interaction (PPI) network Ten platform (100-150 bp) (Illumina Inc, San Diego, construction CA, USA), according to the manufacturer’s protocol, In this study, we used the online database, The by Shanghai Oebiotech Co. Ltd. (Shanghai, China). Search Tool for the Retrieval of Interacting Genes (STRING) (https://string-db.org/cgi/input.pl), to MethylRAD-seq analysis construct the PPI network of differentially methylated After quality-control and filtering of the original genes. Then the PPI network was constructed and reads, high-quality reads were mapped against these visualized using Cytoscape 3.5.1. reference sites (ftp://ftp.ensembl.org/pub/release-84/ fasta/mus_musculus/dna/Mus_musculus.GRCm38. Validation of susceptibility gene expression dna.toplevel.fa.gz) using the SOAP program.11 In by qPCR order to improve accuracy in the follow-up analysis, The same samples used for MethylRAD-seq were Pear software (v0.9.6)12was used to re-filter paired- used for reverse transcription. Briefly, total RNA end sequencing by removing: (i) low quality reads was isolated from mouse palatal shelf tissues and (Phred quality score lower than 30) and (ii) sequences reverse-transcribed into cDNA using TRIzol reagent containing too many N bases (sites containing more (Invitrogen Inc, Carlsbad, CA USA) and the Thermo than 8% of N bases). Signatures containing FspEI First cDNA Synthesis Kit (Thermo Scientific™, sites were extracted from the genome as the reference Waltham, MA, USA), respectively, according to the sequence. Sites covered by at least three reads were manufacturers’ protocols. In each qPCR tube, a 20 μl regarded as authentic methylated sites. The number of reaction mix was prepared using 2×SG Green qPCR Mix methylated sites and the depth of signature coverage (with ROX) (Sino Gene, Beijing, China). The thermal of each sample were then calculated. The untranslated profile was 40 cycles of 95°C for 10 seconds and region (UTR) was calculated using snpEff software 60°C for 30 seconds, followed by a dissociation curve (version: 4.3p),13 and was counted using bed tools check. The qPCR primer sequences for Fgf16 were software (v2.25.0)14 according to the annotation as follows: forward, 5′ACGTGAATGTGTTTTCCGGG3′, document and the distribution of methylation sites in reverse, 5′CCGTCTTTATTCAGGGCCAC3′. The qPCR the different gene elements