1
|
Rodney AR, Buckley RM, Fulton RS, Fronick C, Richmond T, Helps CR, Pantke P, Trent DJ, Vernau KM, Munday JS, Lewin AC, Middleton R, Lyons LA, Warren WC. A domestic cat whole exome sequencing resource for trait discovery. Sci Rep 2021; 11:7159. [PMID: 33785770 DOI: 10.1038/s41598-021-86200-7] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2020] [Accepted: 02/17/2021] [Indexed: 12/13/2022] Open
Abstract
Over 94 million domestic cats are susceptible to cancers and other common and rare diseases. Whole exome sequencing (WES) is a proven strategy to study these disease-causing variants. Presented is a 35.7 Mb exome capture design based on the annotated Felis_catus_9.0 genome assembly, covering 201,683 regions of the cat genome. Whole exome sequencing was conducted on 41 cats with known and unknown genetic diseases and traits, of which ten cats had matching whole genome sequence (WGS) data available, used to validate WES performance. At 80 × mean exome depth of coverage, 96.4% of on-target base coverage had a sequencing depth > 20-fold, while over 98% of single nucleotide variants (SNVs) identified by WGS were also identified by WES. Platform-specific SNVs were restricted to sex chromosomes and a small number of olfactory receptor genes. Within the 41 cats, we identified 31 previously known causal variants and discovered new gene candidate variants, including novel missense variance for polycystic kidney disease and atrichia in the Peterbald cat. These results show the utility of WES to identify novel gene candidate alleles for diseases and traits for the first time in a feline model.
Collapse
|
2
|
Kim H, Yoshihara M, Suyama M. Comparative genomic analysis of inbred rat strains reveals the existence of ancestral polymorphisms. Mamm Genome 2020; 31:86-94. [PMID: 32166433 PMCID: PMC7200647 DOI: 10.1007/s00335-020-09831-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Accepted: 03/02/2020] [Indexed: 11/25/2022]
Abstract
In an alignment of closely related genomic sequences, the existence of discordant mutation sites, which do not reflect the phylogenetic relationship of the genomes, is often observed. Although these discordant mutation sites are thought to have emerged by ancestral polymorphism or gene flow, their frequency and distribution in the genome have not yet been analyzed in detail. Using the genome sequences of all protein coding genes of 25 inbred rat strains, we analyzed the frequency and genome-wide distribution of the discordant mutation sites. From the comparison of different substrains, it was found that these loci are not substrain specific, but are common among different groups of substrains, suggesting that the discordant sites might have mainly emerged through ancestral polymorphism. It was also revealed that the discordant sites are not uniformly distributed along chromosomes, but are concentrated at certain genomic loci, such as RT1, major histocompatibility complex of rats, and olfactory receptors, indicating that genes known to be highly polymorphic tend to have more discordant sites. Our results also showed that loci with a high density of discordant sites are also rich in heterozygous variants, even though these are inbred strains.
Collapse
Affiliation(s)
- Hyeonjeong Kim
- Division of Bioinformatics, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka, 812-8582, Japan
| | - Minako Yoshihara
- Division of Bioinformatics, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka, 812-8582, Japan
| | - Mikita Suyama
- Division of Bioinformatics, Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka, 812-8582, Japan.
| |
Collapse
|
3
|
Hoff SNK, Baalsrud HT, Tooming-Klunderud A, Skage M, Richmond T, Obernosterer G, Shirzadi R, Tørresen OK, Jakobsen KS, Jentoft S. Long-read sequence capture of the haemoglobin gene clusters across codfish species. Mol Ecol Resour 2018; 19:245-259. [PMID: 30329222 PMCID: PMC7379720 DOI: 10.1111/1755-0998.12955] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2018] [Revised: 10/05/2018] [Accepted: 10/09/2018] [Indexed: 11/30/2022]
Abstract
Combining high-throughput sequencing with targeted sequence capture has become an attractive tool to study specific genomic regions of interest. Most studies have so far focused on the exome using short-read technology. These approaches are not designed to capture intergenic regions needed to reconstruct genomic organization, including regulatory regions and gene synteny. Here, we demonstrate the power of combining targeted sequence capture with long-read sequencing technology for comparative genomic analyses of the haemoglobin (Hb) gene clusters across eight species separated by up to 70 million years. Guided by the reference genome assembly of the Atlantic cod (Gadus morhua) together with genome information from draft assemblies of selected codfishes, we designed probes covering the two Hb gene clusters. Use of custom-made barcodes combined with PacBio RSII sequencing led to highly continuous assemblies of the LA (~100 kb) and MN (~200 kb) clusters, which include syntenic regions of coding and intergenic sequences. Our results revealed an overall conserved genomic organization of the Hb genes within this lineage, yet with several, lineage-specific gene duplications. Moreover, for some of the species examined, we identified amino acid substitutions at two sites in the Hbb1 gene as well as length polymorphisms in its regulatory region, which has previously been linked to temperature adaptation in Atlantic cod populations. This study highlights the use of targeted long-read capture as a versatile approach for comparative genomic studies by generation of a cross-species genomic resource elucidating the evolutionary history of the Hb gene family across the highly divergent group of codfishes.
Collapse
Affiliation(s)
- Siv Nam Khang Hoff
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Helle T Baalsrud
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Ave Tooming-Klunderud
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Morten Skage
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | | | | | | | - Ole Kristian Tørresen
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Kjetill S Jakobsen
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Sissel Jentoft
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| |
Collapse
|
4
|
Förster DW, Bull JK, Lenz D, Autenrieth M, Paijmans JLA, Kraus RHS, Nowak C, Bayerl H, Kuehn R, Saveljev AP, Sindičić M, Hofreiter M, Schmidt K, Fickel J. Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species. Mol Ecol Resour 2018; 18:1356-1373. [PMID: 29978939 DOI: 10.1111/1755-0998.12924] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2017] [Revised: 05/28/2018] [Accepted: 06/05/2018] [Indexed: 11/29/2022]
Abstract
Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.
Collapse
Affiliation(s)
- Daniel W Förster
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
| | - James K Bull
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
| | - Dorina Lenz
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
| | - Marijke Autenrieth
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany.,Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany
| | | | - Robert H S Kraus
- Department of Biology, University of Konstanz, Konstanz, Germany.,Department of Migration and Immuno-Ecology, Max Planck Institute for Ornithology, Radolfzell, Germany
| | - Carsten Nowak
- Conservation Genetics Group, Senckenberg Research Institute and Natural History Museum Frankfurt, Gelnhausen, Germany
| | - Helmut Bayerl
- Unit of Molecular Zoology, Chair of Zoology, Department of Animal Science, Technical University of Munich, Freising, Germany
| | - Ralph Kuehn
- Unit of Molecular Zoology, Chair of Zoology, Department of Animal Science, Technical University of Munich, Freising, Germany.,Department of Fish, Wildlife and Conservation Ecology, New Mexico State University, Las Cruces, New Mexico
| | - Alexander P Saveljev
- Department of Animal Ecology, Russian Research Institute of Game Management and Fur Farming, Kirov, Russia
| | - Magda Sindičić
- Faculty of Veterinary Medicine, University of Zagreb, Zagreb, Croatia
| | - Michael Hofreiter
- Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany
| | - Krzysztof Schmidt
- Mammal Research Institute, Polish Academy of Sciences, Białowieza, Poland
| | - Jörns Fickel
- Department of Evolutionary Genetics, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany.,Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany
| |
Collapse
|
5
|
Foley JF, Phadke DP, Hardy O, Hardy S, Miller V, Madan A, Howard K, Kruse K, Lord C, Ramaiahgari S, Solomon GG, Shah RR, Pandiri AR, Herbert RA, Sills RC, Merrick BA. Whole exome sequencing in the rat. BMC Genomics 2018; 19:487. [PMID: 29925311 DOI: 10.1186/s12864-018-4858-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2018] [Accepted: 06/06/2018] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND The rat genome was sequenced in 2004 with the aim to improve human health altered by disease and environmental influences through gene discovery and animal model validation. Here, we report development and testing of a probe set for whole exome sequencing (WES) to detect sequence variants in exons and UTRs of the rat genome. Using an in-silico approach, we designed probes targeting the rat exome and compared captured mutations in cancer-related genes from four chemically induced rat tumor cell lines (C6, FAT7, DSL-6A/C1, NBTII) to validated cancer genes in the human database, Catalogue of Somatic Mutations in Cancer (COSMIC) as well as normal rat DNA. Paired, fresh frozen (FF) and formalin-fixed, paraffin-embedded (FFPE) liver tissue from naive rats were sequenced to confirm known dbSNP variants and identify any additional variants. RESULTS Informatics analysis of available gene annotation from rat RGSC6.0/rn6 RefSeq and Ensembl transcripts provided 223,636 unique exons representing a total of 26,365 unique genes and untranslated regions. Using this annotation and the Rn6 reference genome, an in-silico probe design generated 826,878 probe sequences of which 94.2% were uniquely aligned to the rat genome without mismatches. Further informatics analysis revealed 25,249 genes (95.8%) covered by at least one probe and 23,603 genes (93.5%) had every exon covered by one or more probes. We report high performance metrics from exome sequencing of our probe set and Sanger validation of annotated, highly relevant, cancer gene mutations as cataloged in the human COSMIC database, in addition to several exonic variants in cancer-related genes. CONCLUSIONS An in-silico probe set was designed to enrich the rat exome from isolated DNA. The platform was tested on rat tumor cell lines and normal FF and FFPE liver tissue. The method effectively captured target exome regions in the test DNA samples with exceptional sensitivity and specificity to obtain reliable sequencing data representing variants that are likely chemically induced somatic mutations. Genomic discovery conducted by means of high throughput WES queries should benefit investigators in discovering rat genomic variants in disease etiology and in furthering human translational research.
Collapse
|
6
|
Meckes JK, Lim PH, Wert SL, Luo W, Gacek SA, Platt D, Jankord R, Saar K, Redei EE. Brain region-specific expression of genes mapped within quantitative trait loci for behavioral responsiveness to acute stress in Fisher 344 and Wistar Kyoto male rats. PLoS One 2018; 13:e0194293. [PMID: 29529077 PMCID: PMC5847310 DOI: 10.1371/journal.pone.0194293] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2017] [Accepted: 02/28/2018] [Indexed: 12/27/2022] Open
Abstract
Acute stress responsiveness is a quantitative trait that varies in severity from one individual to another; however, the genetic component underlying the individual variation is largely unknown. Fischer 344 (F344) and Wistar Kyoto (WKY) rat strains show large differences in behavioral responsiveness to acute stress, such as freezing behavior in response to footshock during the conditioning phase of contextual fear conditioning (CFC). Quantitative trait loci (QTL) have been identified for behavioral responsiveness to acute stress in the defensive burying (DB) and open field test (OFT) from a reciprocal F2 cross of F344 and WKY rat strains. These included a significant QTL on chromosome 6 (Stresp10). Here, we hypothesized that the Stresp10 region harbors genes with sequence variation(s) that contribute to differences in multiple behavioral response phenotypes between the F344 and WKY rat strains. To test this hypothesis, first we identified differentially expressed genes within the Stresp10 QTL in the hippocampus, amygdala, and frontal cortex of F344 and WKY male rats using genome-wide microarray analyses. Genes with both expression differences and non-synonymous sequence variations in their coding regions were considered candidate quantitative trait genes (QTGs). As a proof-of-concept, the F344.WKY-Stresp10 congenic strain was generated with the Stresp10 WKY donor region into the F344 recipient strain. This congenic strain showed behavioral phenotypes similar to those of WKYs. Expression patterns of Gpatch11 (G-patch domain containing 11), Cdkl4 (Cyclin dependent kinase like 4), and Drc1 (Dynein regulatory complex subunit 1) paralleled that of WKY in the F344.WKY-Stresp10 strain matching the behavioral profiles of WKY as opposed to F344 parental strains. We propose that these genes are candidate QTGs for behavioral responsiveness to acute stress.
Collapse
Affiliation(s)
- Jeanie K. Meckes
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Patrick H. Lim
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Stephanie L. Wert
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Wendy Luo
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Stephanie A. Gacek
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Dana Platt
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| | - Ryan Jankord
- Applied Neuroscience, Warfighter Interface Division, 711th Human Performance Wing, Air Force Research Laboratory, Wright-Patterson Air Force Base, OH, United States of America
| | - Kathrin Saar
- Max Delbrück Center for Molecular Medicine in the Helmholtz Association, Berlin-Buch, Germany
| | - Eva E. Redei
- Department of Psychiatry and Behavioral Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States of America
| |
Collapse
|
7
|
Kuramoto T, Voigt B, Nakanishi S, Kitada K, Nakamura T, Wakamatsu K, Yoshihara M, Suyama M, Uemura R, Tanaka M, Kuwamura M, Shimizu S, Ohno Y, Sasa M, Serikawa T. Identification of Candidate Genes for Generalized Tonic-Clonic Seizures in Noda Epileptic Rat. Behav Genet 2017; 47:609-19. [PMID: 28936718 DOI: 10.1007/s10519-017-9870-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2017] [Accepted: 09/08/2017] [Indexed: 12/31/2022]
Abstract
The Noda epileptic rat (NER) exhibits generalized tonic-clonic seizures (GTCS). A genetic linkage analysis identified two GTCS-associated loci, Ner1 on Chr 1 and Ner3 on Chr 5. The wild-type Ner1 and Ner3 alleles suppressed GTCS when combined in double-locus congenic lines, but not when present in single-locus congenic lines. Global expression analysis revealed that cholecystokinin B receptor (Cckbr) and suppressor of tumorigenicity 5 (St5), which map within Ner1, and PHD finger protein 24 (Phf24), which maps within Ner3, were significantly downregulated in NER. De novo BAC sequencing detected an insertion of an endogenous retrovirus sequence in intron 2 of the Phf24 gene in the NER genome, and PHF24 protein was almost absent in the NER brain. Phf24 encodes a Gαi-interacting protein involved in GABAB receptor signaling pathway. Based on these findings, we conclude that Cckbr, St5, and Phf24 are strong candidate genes for GTCS in NER.
Collapse
|
8
|
Yoshihara M, Sato T, Saito D, Ohara O, Kuramoto T, Suyama M. A deletion in the intergenic region upstream of Ednrb causes head spot in the rat strain KFRS4/Kyo. BMC Genet 2017; 18:29. [PMID: 28356074 PMCID: PMC5372274 DOI: 10.1186/s12863-017-0497-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2017] [Accepted: 03/25/2017] [Indexed: 11/17/2022] Open
Abstract
Background Head spot is one of the phenotypes identified in the KFRS4/Kyo rat strain. Although previous linkage analysis suggested that Ednrb, which is frequently involved in coat color variations in various animals, could be the gene responsible for this phenotype, no mutations have been identified in its coding region. Results To identify mutations causative of this phenotype in KFRS4/Kyo, we analyzed target capture sequencing data that we recently generated. Our target capture method has a unique feature, i.e., it covers not only exonic regions but also conserved non-coding sequences (CNSs) among vertebrates; therefore, it has the potential to detect regulatory mutations. We identified a deletion of approximately 50 kb in length approximately 50 kb upstream of Ednrb. A comparative analysis with the epigenomic data in the corresponding region in humans and mice showed that one of the CNSs might be an enhancer. Further comparison with Hi-C data, which provide information about chromosome conformation, indicated that the putative enhancer is spatially close to the promoter of Ednrb, suggesting that it acts as an enhancer of Ednrb. Conclusions These in silico data analyses strongly suggest that the identified deletion in the intergenic region upstream of Ednrb, which might contain a melanocyte-specific enhancer, is the mutation causative of the head spot phenotype in the KFRS4/Kyo rat strain. Electronic supplementary material The online version of this article (doi:10.1186/s12863-017-0497-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Minako Yoshihara
- Medical Institute of Bioregulation, Kyushu University, Fukuoka, 812-8582, Japan.,AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka, 812-8582, Japan
| | - Tetsuya Sato
- Medical Institute of Bioregulation, Kyushu University, Fukuoka, 812-8582, Japan.,AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka, 812-8582, Japan
| | - Daisuke Saito
- Medical Institute of Bioregulation, Kyushu University, Fukuoka, 812-8582, Japan.,AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka, 812-8582, Japan
| | - Osamu Ohara
- Department of Technology Development, Kazusa DNA Research Institute, Kisarazu, Chiba, 292-0818, Japan
| | - Takashi Kuramoto
- Institute of Laboratory Animals, Graduate School of Medicine, Kyoto University, Kyoto, 606-8501, Japan
| | - Mikita Suyama
- Medical Institute of Bioregulation, Kyushu University, Fukuoka, 812-8582, Japan. .,AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka, 812-8582, Japan.
| |
Collapse
|
9
|
Yoshihara M, Sato T, Saito D, Ohara O, Kuramoto T, Suyama M. Application of target capture sequencing of exons and conserved non-coding sequences to 20 inbred rat strains. Genom Data 2016; 10:155-157. [PMID: 27882299 PMCID: PMC5114524 DOI: 10.1016/j.gdata.2016.11.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 10/28/2016] [Accepted: 11/09/2016] [Indexed: 11/23/2022]
Abstract
We report sequence data obtained by our recently devised target capture method TargetEC applied to 20 inbred rat strains. This method encompasses not only all annotated exons but also highly conserved non-coding sequences shared among vertebrates. The total length of the target regions covers 146.8 Mb. On an average, we obtained 31.7 × depth of target coverage and identified 154,330 SNVs and 24,368 INDELs for each strain. This corresponds to 470,037 unique SNVs and 68,652 unique INDELs among the 20 strains. The sequence data can be accessed at DDBJ/EMBL/GenBank under accession number PRJDB4648, and the identified variants have been deposited at http://bioinfo.sls.kyushu-u.ac.jp/rat_target_capture/20_strains.vcf.gz.
Collapse
Affiliation(s)
- Minako Yoshihara
- Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan; AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka 812-8582, Japan
| | - Tetsuya Sato
- Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan; AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka 812-8582, Japan
| | - Daisuke Saito
- Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan; AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka 812-8582, Japan
| | - Osamu Ohara
- Department of Technology Development, Kazusa DNA Research Institute, Kisarazu, Chiba 292-0818, Japan
| | - Takashi Kuramoto
- Institute of Laboratory Animals, Graduate School of Medicine, Kyoto University, Kyoto 606-8501, Japan
| | - Mikita Suyama
- Medical Institute of Bioregulation, Kyushu University, Maidashi 3-1-1, Higashi-ku, Fukuoka 812-8582, Japan; AMED-CREST, Japan Agency for Medical Research and Development, Fukuoka 812-8582, Japan
| |
Collapse
|