==== Front Neuropsychopharmacol Rep Neuropsychopharmacol Rep 10.1002/(ISSN)2574-173X NPR2 Neuropsychopharmacology Reports 2574-173X John Wiley and Sons Inc. Hoboken 33037870 10.1002/npr2.12145 NPR212145 Micro Report Micro Reports Use of the Illumina EPIC methylation array for epigenomic research in the crab‐eating macaque (Macaca fascicularis) NAKACHI et al.Nakachi Yutaka https://orcid.org/0000-0003-1964-2495 1 Ishii Kazuhiro 2 Bundo Miki 1 Masuda Tomoyuki https://orcid.org/0000-0001-6429-9468 2 tmasu@md.tsukuba.ac.jp Iwamoto Kazuya https://orcid.org/0000-0002-1780-692X 1 tmasu@md.tsukuba.ac.jpiwamotok@kumamoto-u.ac.jp 1 Department of Molecular Brain Science Graduate School of Medical Sciences Kumamoto University Kumamoto Japan 2 Department of Neurology Faculty of Medicine University of Tsukuba Ibaraki Japan * Correspondence Tomoyuki Masuda, Department of Neurology, Faculty of Medicine, University of Tsukuba, Ibaraki 305‐8575, Japan. Email: tmasu@md.tsukuba.ac.jp Kazuya Iwamoto, Department of Molecular Brain Science, Graduate School of Medical Sciences, Kumamoto University, Kumamoto 860‐8556, Japan. Email: iwamotok@kumamoto-u.ac.jp 10 10 2020 12 2020 40 4 10.1111/npr2.v40.4423 426 12 7 2020 26 8 2020 10 9 2020 © 2020 The Authors. Neuropsychopharmacology Reports published by John Wiley & Sons Australia, Ltd on behalf of The Japanese Society of Neuropsycho PharmacologyThis is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.Abstract Background Commercially available Illumina DNA methylation arrays (HumanMethylation 27K, HumanMethylation450, and MethylationEPIC BeadChip) can be used for comprehensive DNA methylation analyses of not only the human genome but also other mammalian genomes, ranging from those of nonhuman primates to those of rodents. However, practical application of the EPIC array to the crab‐eating macaque has not been reported. Methods Through bioinformatic analyses involving cross‐species comparison and consideration of probe performance, we selected array probes that can be reliably used for the crab‐eating macaque genome. A DNA methylation assay using an EPIC array was performed on genomic DNA extracted from the brains of five crab‐eating macaques. The obtained DNA methylation data were compared with a publicly available dataset. Results Among the 865 918 probes in the EPIC array, a total of 183 509 probes (21.2%) were selected as high‐confidence array probes in the crab‐eating macaque. Subsequent comparisons revealed that the data from these probes showed good concordance with other DNA methylation datasets of the crab‐eating macaque. Conclusion The selected high‐confidence array probes would be useful for high‐throughput DNA methylation assays of the crab‐eating macaque. Epigenetic research in the non‐human primates, such as crab‐eating macaque, will be important to understand the pathophysiology of psychiatric disorders. Among the methylation array probes for human genome, the probes that can reliably measure DNA methylation levels of the crab‐eating macaque are reported. bisulfite modificationcynomolgus monkeycytosineDNA methylationnonhuman primateJSPS KAKENHI16H0639917K1028318H0275318H0542818H0543018K0756720K07240Japan Agency for Medical Research and Development 10.13039/100009619JP20dm0107123JP20dm0207074 source-schema-version-number2.0cover-dateDecember 2020details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:5.9.5 mode:remove_FC converted:08.12.2020 Nakachi Y , Ishii K , Bundo M , Masuda T , Iwamoto K . Use of the Illumina EPIC methylation array for epigenomic research in the crab‐eating macaque (Macaca fascicularis) . Neuropsychopharmacol Rep 2020 ;40 :423 –426 . 10.1002/npr2.12145 33037870 ==== Body 1 INTRODUCTION DNA methylation is a central epigenetic mechanism that regulates long‐lasting gene expression and is considered to reflect the complex interaction between genes and environments. Altered DNA methylation in the brain has been reported in a wide range of neuropsychiatric disorders and their animal models. Thus, elucidating the comprehensive DNA methylation profile may be quite important for understanding pathophysiology. 1 , 2 The commercially available array‐based DNA methylation assay, which allows quantification of DNA methylation levels at hundreds of thousands of CpG sites, has been widely used in humans. Although the assay was specifically designed for the human genome, previous studies have reported its application to the genomes of other species, ranging from nonhuman primates to rodents. 3 , 4 , 5 , 6 , 7 , 8 , 9 Typically, this was done by assessing the conservation of genomic context around the probe region between humans and the species of interest by a bioinformatic approach. The crab‐eating macaque (cynomolgus monkey, Macaca fascicularis) is a cercopithecine primate in the group of Old World monkeys and is widely distributed in Southeast Asia. Due to its ease of rearing and propagation, as well as the close relatedness of its brain structure and function to those of humans, this macaque has become one of the most widely used nonhuman primates in biomedical research. Previously, applications of HumanMethylation27K and HumanMethylation450K arrays (Illumina), which contain approximately 27 000 and 450 000 array probes, respectively, to the crab‐eating macaque have been reported. 3 , 8 Here, we report the application of the updated version of the array, MethylationEPIC (Illumina), containing approximately 850 000 probes, to the crab‐eating macaque. 2 MATERIALS AND METHODS 2.1 Animals Five cynomolgus monkeys (Macaca fascicularis; 2 males and 3 females, 3‐4 years old) were purchased from HAMRI Co., Ltd. All procedures for animal care and experimentation were approved by the University of Tsukuba Animal Experiment Committee and were carried out in accordance with the Guidelines for Proper Conduct of Animal Experiments by the Science Council of Japan. 2.2 Probe filtering and selection strategy Probe selection was performed according to our previous report. 9 In brief, we used the manifest file provided by the manufacturer (Illumina). The information for SourceSeq, ProbeSeq, MapInfo, and AlleleA ProbeSeq was provided in the manifest file (Infinium MethylationEPIC v1.0 B4), which indicated the bisulfite‐unconverted design sequence, the probe sequence, the coordinates of the target site, and the sequence for probe A, respectively. We retrieved human bisulfite‐unconverted probe sequences (referred to as modified SourceSeq) based on the MapInfo and AlleleA ProbeSeq. Illumina arrays include two types of array probes. Type I probe contains two probes designed for methylated and unmethylated CpG at one CpG site, whereas type II probe contains single probe for calculation of methylation level at one CpG site. For type I and type II probes, we retrieved 50 and 51 bp of the human sequence, respectively. We mapped the modified SourceSeq to the reference genome of the crab‐eating macaque, Macaca_fascicularis_5.0 (macFas5, GCA_000364345.1 without the mitochondrial sequence). We used Blat 10 for mapping and selected probes with the following criteria: (a) no insertions/deletions; (b) only one best hit and no suboptimal hits; (c) no mismatches within 5 bp of the target site, with up to two mismatches allowed in the rest of the sequence; and (d) no ambiguous bases in the mapped reference genome. All location descriptions of mapped probes and the relationship of orthologous genes between humans and the crab‐eating macaque were confirmed by using the Ensembl release 96 dataset corresponding to the Macaca_fascicularis_5.0 genome. Selected probes were visualized with a custom viewer (CHARANGO viewer; CHromosomal Alignment Representation and ANnotation Graphical Omics viewer). Scripts are available upon request. 2.3 DNA methylation assay Genomic DNA was extracted from frozen occipital lobes using the AllPrep DNA/RNA Mini Kit (QIAGEN). DNA methylation profiles were obtained using an Illumina MethylationEPIC kit (Illumina). All experimental methods were performed strictly according to the manufacturer's protocols (Illumina). 2.4 Data analysis The IDAT file, a raw data file containing the intensities of the probes in the array, was analyzed using the minfi package 11 of R/Bioconductor and noob background correction. 12 The probes for which detection P values were <.01 across all samples (N = 5) were used. Pairwise Pearson's correlations of beta values, which are proxies for DNA methylation level and range from 0 to 1, were calculated. Average beta values were used for comparison with the published dataset. 13 3 RESULTS AND DISCUSSION Among the 865 918 probes in the HumanMethylationEPIC array, we selected the probes that could hybridize with the bisulfite‐converted crab‐eating macaque genome (macFas5). As a result, we obtained 183 509 high‐confidence probes in total (Figure 1, Table S1). These probes account for 21.2% of all the probes in the array. In our similar calculations in the common marmoset 9 and the chimpanzee (Nakachi et al, unpublished data), we obtained approximately 9% and 44% high‐confidence probes, respectively. Considering the evolutionary relationships, the number of obtained probes in this study seems reasonable. We added annotation based on macFas5 to each probe, by which one can consider orthologous relationships between humans and the crab‐eating macaque (Table S1). FIGURE 1 Chromosomal distribution of high‐confidence EPIC array probes. Chromosomal locations of high‐confidence probes (N = 183 509) are shown. Both type I (N = 37 214, left‐hand side of each chromosome) and type II (N = 146 295, right‐hand side of each chromosome) probes are illustrated The high‐confidence probes were mapped to the entire macaque genome for both type I (N = 37 214) and type II (N = 146 295) probes (Figure 1). To validate the performance of the high‐confidence probes, we profiled DNA methylation patterns of the crab‐eating macaque brains using the EPIC array (occipital lobes, N = 5). In this analysis, we used autosomal CpG probes (N = 174 611) that passed detection P values (<.01) across five samples. The finally selected probes (N = 172 688) showed high concordance not only for all probes (total average R = 0.980) but also for both type I (R = 0.989) and type II (R = 0.974) probes among five biological replicates. We then compared the DNA methylation values of high‐confidence probes and independent DNA methylation datasets for the crab‐eating macaque. Because the DNA methylation study using the occipital lobe of the crab‐eating macaque has not been released, we utilized a previous study that reported methylated regions of seven different brain regions of the crab‐eating macaque. 13 Within the 1310 methylated regions reported, 13 a total of 2844 high‐confidence CpG probes were included. Comparison of DNA methylation levels revealed high consistency with different brain regions (R = 0.709 for the hippocampus to 0.743 for the inferior temporal lobe, average R for the seven brain regions was 0.733) (Figure 2). Similarly, both type I and type II probes also showed consistency, with average R values = .788 and .691, respectively. Because we utilized the published dataset of different brain regions of different macaques, and we did not perform validation experiments using the same brain region of same macaques, correlations remained the modest levels. FIGURE 2 Comparison of DNA methylation levels between high‐confidence probes and other datasets. Average beta values (N = 5) of high‐confidence probes and methylation levels reported in other datasets 13 are plotted. Among the comparisons with seven different brain regions, data on the inferior temporal lobe are shown as an example We provided high‐confidence Illumina EPIC array probes that can be used in the crab‐eating macaque genome. Using the well‐established, low‐cost platform, epigenomic research related to neuropsychiatric disorders will be facilitated in the crab‐eating macaque. In particular, given the high conservation with the human genome in each array probe, the probes can be readily used for interspecies comparisons. In this study, we used the occipital lobes, however, other tissues such as blood can be applied similarly, as we demonstrated previously in other primate. 9 The limitations of DNA methylation assay using the methylation array are (a) it can analyze a limited number of the CpGs, compared to other approaches, such as whole‐genome bisulfite sequencing, (b) it cannot discriminate the different cytosine modifications solely, such as methylcytosine and hydroxymethyl cytosine, and (c) its quantitative data can be affected the presence of the SNPs on the array probes. It should be noted that a lack of polymorphism information in nonhuman primates might cause false methylation signals due to hybridization artifacts (Nakachi et al, unpublished data). Further refinement of the high‐confidence probes awaits the accumulation of more detailed genomic information on the crab‐eating macaque. CONFLICT OF INTEREST The authors declare no conflicts of interest. AUTHOR CONTRIBUTIONS TM, MB, and K. Iwamoto designed the research. K. Ishii and TM performed the experiment. YN analyzed the data. YN, MB, TM, and K. Iwamoto prepared the manuscript. Supporting information Tab S1 Click here for additional data file. ACKNOWLEDGMENTS We thank M. Iseki for her assistance. This research was supported by a Grant‐in‐Aid for research on study (Toxic study group) of the Health Effects on Diphenylarsinic Acid supported by the Ministry of the Environment, Japan, and in part by JSPS KAKENHI grant numbers 17K10283, 18H05428, 18H05430, 18H02753, 16H06399, 18K07567, and 20K07240. This research was also partly supported by AMED under grant numbers JP20dm0107123 and JP20dm0207074. DATA AVAILABILITY STATEMENT Raw data have been deposited in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) and are accessible through GSE154078 14 . ==== Refs REFERENCES 1 Iwamoto K . Understanding the epigenetic architecture of psychiatric disorders: modifications and beyond . Psychiatry Clin Neurosci . 2018 ;72 :194 .29633529 2 Kular L , Kular S . Epigenetics applied to psychiatry: clinical opportunities and future challenges . Psychiatry Clin Neurosci . 2018 ;72 :195 –211 .29292553 3 Howard TD , Ho SM , Zhang L , Chen J , Cui W , Slager R , et al. Epigenetic changes with dietary soy in cynomolgus monkeys . PLoS ONE . 2011 ;6 :e26791.22046358 4 Iwamoto K , Bundo M , Ueda J , Oldham MC , Ukai W , Hashimoto E , et al. Neurons show distinctive DNA methylation pro‐file and higher interindividual variations compared with non‐neurons . Genome Res . 2011 ;21 :688 –96 .21467265 5 Hernando‐Herraez I , Prado‐Martinez J , Garg P , Fernandez‐Callejo M , Heyn H , Hvilsom C , et al. Dynamics of DNA methylation in recent human and great ape evolution . PLoS Genet . 2013 ;9 :e1003763.24039605 6 Wong NC , Ng J , Hall NE , Lunke S , Salmanidis M , Brumatti G , et al. Exploring the utility of human DNA methylation arrays for profiling mouse genomic DNA . Genomics . 2013 ;102 :38 –46 .23639479 7 Chopra P , Papale LA , White AT , Hatch A , Brown RM , Garthwaite MA , et al. Array‐based assay detects genome‐wide 5‐mC and 5‐hmC in the brains of humans, non‐human primates, and mice . BMC Genom . 2014 ;15 :131 . 8 Ong ML , Tan PY , MacIsaac JL , Mah SM , Buschdorf JP , Cheong CY , et al. Infinium monkeys: Infinium 450K array for the Cynomolgus macaque (Macaca fascicularis) . G3 (Bethesda) . 2014 ;4 :1227 –34 .24815017 9 Ueda J , Murata Y , Bundo M , Oh‐Nishi A , Kasai H , Ikegame T , et al. Use of human methylation arrays for epigenome research in the common marmoset (Callithrix jacchus) . Neurosci Res . 2017 ;120 :60 –5 .28215819 10 Kent WJ . BLAT ‐ The BLAST‐Like alignment tool . Genome Res . 2002 ;12 :656 –64 .11932250 11 Aryee MJ , Jaffe AE , Corrada‐Bravo H , Ladd‐Acosta C , Feinberg AP , Hansen KD , et al. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays . Bioinformatics . 2014 ;30 :1363 –9 .24478339 12 Triche TJ , Weisenberger DJ , Van Den Berg D , Laird PW , Siegmund KD . Low‐level processing of Illumina Infinium DNA Methylation BeadArrays . Nucleic Acids Res . 2013 ;41 :e90.23476028 13 Lim WJ , Kim KH , Kim JY , Kim HJ , Kim M , Park JL , et al. Investigation of gene expression and DNA methylation from seven different brain regions of a crab‐eating monkey as determined by RNA‐Seq and whole‐genome bisulfite sequencing . Front Genet . 2019 ;10 :694 .31428131 14 Nakachi Y , Ishii K , Bundo M , Masuda T , Iwamoto K . Use of the Illumina EPIC methylation array for epigenomic research in the crab‐eating macaque (Macaca fascicularis); Gene Expression Omnibus (GEO); GSE154078 ; 2020 .