Progress report of the Tohoku Medical Megabank Community-Based Cohort Study: Study profile of the repeated center-based survey during second period in Miyagi Prefecture

BACKGROUND

The purpose of this study was to report the basic profile of the Miyagi Prefecture part of a repeated center-based survey during the second period (2nd period survey) of the Tohoku Medical Megabank Community-Based Cohort Study (TMM CommCohort Study), as well as the participants' characteristics based on their participation type in the baseline survey.

METHODS

The 2nd period survey, conducted from June 2017 to March 2021, included participants of the TMM CommCohort Study (May 2013 to March 2016). In addition to the questionnaire, blood, urine, and physiological function tests were performed during the 2nd period survey. There were three main ways of participation in the baseline survey: Type 1, Type 1 additional, or Type 2 survey. The 2nd period survey was conducted in the same manner as the Type 2 survey, which was based on the community support center (CSC).

RESULTS

In Miyagi Prefecture, 29,383 (57.7%) of 50,967 participants participated in the 2nd period survey. The participation rate among individuals who had visited the CSC was approximately 80%. Although some factors differed depending on the participation type in the baseline survey, the 2nd period survey respondents in the Type 1 and Type 2 survey groups at baseline had similar traits.

CONCLUSIONS

The 2nd period survey of the TMM CommCohort Study provided detailed follow-up information. Following up on the health conditions of the participants will clarify the long-term effects of disasters and contribute to personalized prevention.

A Pilot Study for Return of Individual Pharmacogenomic Results to Population-Based Cohort Study Participants

Introduction:

Pharmacogenomic (PGx) testing results provide valuable information on drug selection and appropriate dosing, maximization of efficacy, and minimization of adverse effects. Although the number of large-scale, next-generation-sequencing-based PGx studies has recently increased, little is known about the risks and benefits of returning PGx results to ostensibly healthy individuals in research settings.

Methods:

Single-nucleotide variants of three actionable PGx genes, namely, MT-RNR1, CYP2C19, and NUDT15, were returned to 161 participants in a population-based Tohoku Medical Megabank project. Informed consent was obtained from the participants after a seminar on the outline of this study. The results were sent by mail alongside sealed information letter intended for clinicians. As an exception, genetic counseling was performed for the MT-RNR1 m.1555A > G variant carriers by a medical geneticist, and consultation with an otolaryngologist was encouraged. Questionnaire surveys (QSs) were conducted five times to evaluate the participants’ understanding of the topic, psychological impact, and attitude toward the study.

Results:

Whereas the majority of participants were unfamiliar with the term PGx, and none had undergone PGx testing before the study, more than 80% of the participants felt that they could acquire basic PGx knowledge sufficient to understand their genomic results and were satisfied with their potential benefit and use in future prescriptions. On the other hand, some felt that the PGx concepts or terminology was difficult to fully understand and suggested that in-person return of the results was desirable.

Conclusions:

These results collectively suggest possible benefits of returning preemptive PGx information to ostensibly healthy cohort participants in a research setting.

Japonica Array NEO with increased genome-wide coverage and abundant disease risk SNPs

Ethnic-specific SNP arrays are becoming more important to increase the power of genome-wide association studies in diverse population. In the Tohoku Medical Megabank Project, we have been developing a series of Japonica Arrays (JPA) for genotyping participants based on reference panels constructed from whole-genome sequence data of the Japanese population. Here, we designed a novel version of the SNP array for the Japanese population, called Japonica Array NEO (JPA NEO), comprising a total of 666,883 markers. Among them, 654,246 tag SNPs of autosomes and X chromosome were selected from an expanded reference panel of 3,552 Japanese, 3.5KJPNv2, using pairwise r² of linkage disequilibrium measures. Additionally, 28,298 markers were included for the evaluation of previously identified disease risk markers from the literature and databases, and those present in the Japanese population were extracted using the reference panel. Through genotyping 286 Japanese samples, we found that the imputation quality r² and INFO score in the minor allele frequency bin >2.5–5% were >0.9 and >0.8, respectively, and >12 million markers were imputed with an INFO score >0.8. From these results, JPA NEO is a promising tool for genotyping the Japanese population with genome-wide coverage, contributing to the development of genetic risk scores.

Oral Microbiome Analysis in Prospective Genome Cohort Studies of the Tohoku Medical Megabank Project

A baseline oral microbiome study of the Tohoku Medical Megabank Organization (TMM) was planned to characterize the profile of the oral microbiome in the Japanese population. The study also aimed to clarify risk factors for multifactorial diseases by integrated analysis of the oral microbiome and host genome/omics information. From 2013 to 2016, we collected three types of oral biospecimens, saliva, supragingival plaque, and tongue swab, from a total of 25,101 participants who had a dental examination in TMM. In this study, we used two independent cohorts; the Community-Based Cohort and Birth and Three-Generation Cohort as discovery and validation cohorts, respectively, and we selected participants examined by a single dentist. We found through the 16S ribosomal RNA gene sequencing analysis of 834 participants of the Community-Based Cohort Study that there are differences in the microbial composition and community structure between saliva and plaque. The species diversities in both saliva and plaque were increased in correlation with the severity of periodontal disease. These results were nicely reproduced in the analysis of 455 participants of the Birth and Three-Generation Cohort Study. In addition, strong positive and negative associations of microbial taxa in both plaque and saliva with periodontitis-associated biofilm formation were detected by co-occurrence network analysis. The classes Actinobacteria and Bacilli, including oral health-associated bacterial species, showed a positive correlation in saliva. These results revealed differences in microbial composition and community structure between saliva and plaque and a correlation between microbial species and the severity of periodontal disease. We expect that the large database of the oral microbiome in the TMM biobank will help in the discovery of novel targets for the treatment and prevention of oral diseases, as well as for the discovery of therapeutic and/or preventive targets of systemic diseases.

Subtype-specific gout susceptibility loci and enrichment of selection pressure on ABCG2 and ALDH2 identified by subtype genome-wide meta-analyses of clinically defined gout patients

Objectives

Genome-wide meta-analyses of clinically defined gout were performed to identify subtype-specific susceptibility loci. Evaluation using selection pressure analysis with these loci was also conducted to investigate genetic risks characteristic of the Japanese population over the last 2000–3000 years.

Methods

Two genome-wide association studies (GWASs) of 3053 clinically defined gout cases and 4554 controls from Japanese males were performed using the Japonica Array and Illumina Array platforms. About 7.2 million single-nucleotide polymorphisms were meta-analysed after imputation. Patients were then divided into four clinical subtypes (the renal underexcretion type, renal overload type, combined type and normal type), and meta-analyses were conducted in the same manner. Selection pressure analyses using singleton density score were also performed on each subtype.

Results

In addition to the eight loci we reported previously, two novel loci, PIBF1 and ACSM2B, were identified at a genome-wide significance level (p<5.0×10^–8) from a GWAS meta-analysis of all gout patients, and other two novel intergenic loci, CD2-PTGFRN and SLC28A3-NTRK2, from normal type gout patients. Subtype-dependent patterns of Manhattan plots were observed with subtype GWASs of gout patients, indicating that these subtype-specific loci suggest differences in pathophysiology along patients’ gout subtypes. Selection pressure analysis revealed significant enrichment of selection pressure on ABCG2 in addition to ALDH2 loci for all subtypes except for normal type gout.

Conclusions

Our findings on subtype GWAS meta-analyses and selection pressure analysis of gout will assist elucidation of the subtype-dependent molecular targets and evolutionary involvement among genotype, phenotype and subtype-specific tailor-made medicine/prevention of gout and hyperuricaemia.

Establishment of Integrated Biobank for Precision Medicine and Personalized Healthcare: The Tohoku Medical Megabank Project

The Tohoku Medical Megabank (TMM) project was established to provide creative reconstruction of the Tohoku area that suffered from a huge earthquake and ensuing tsunami (the Great East Japan Earthquake, GEJE). TMM aims to establish two large-scale genome cohorts and an integrated biobank managing biospecimen and related information. It supports community medicine by establishing next-generation medical systems through a combination of the prospective genome cohort studies with a total of 150,000 participants and genomic medicine. The strategies for genome analyses in TMM are to develop an elaborate genome reference panel by means of high-fidelity Japanese whole-genome sequence, to design custom single nucleotide polymorphism (SNP) arrays based on the reference panel, and to obtain genotype data for all the TMM cohort participants subsequently. Disease-associated genomic information and omics data, including metabolomics and microbiome analysis, provide an essential platform for precision medicine and personalized healthcare (PHC). Ethical, legal, and social issues (ELSI) and education are important for implementing genomic medicine. The major considerations of ELSI regarding each participant of the cohort studies are the respect for the autonomy and the protection of privacies. Moreover, developing and provide human resources not only for the TMM project but also for the social implementation of precision medicine and PHC is required. We started a pilot study of the return of genomic results for familial hypercholesterolemia (FH) as a target disease. TMM aims to establish solid platforms that support precision medicine and PHC based on the genomic and omics information and environmental and lifestyle factors of the individuals, which is one of the most advanced medical care beyond the evidenced-based medicine in the near future.

A training and education program for genome medical research coordinators in the genome cohort study of the Tohoku Medical Megabank Organization

BACKGROUND

Genome cohort studies are used to analyze interactions between genetic and environmental factors, providing valuable information for personalized healthcare. Large-scale and long-term cohort studies require a number of specially trained personnel, of whom those involved in obtaining informed consent play a vital role, especially during the initial phase of such studies. The Japanese Society of Human Genetics (JSHG) previously established a certification system for genome medical research coordinators (GMRCs) responsible for obtaining written consent via face-to-face explanation. Meanwhile, in the Tohoku Medical Megabank Organization (ToMMo), GMRCs are expected to play important roles not only in obtaining informed consent and conducting various assessments, but also in communicating with participants throughout the long-term follow-up. Based on the JSHG program, we therefore developed a specific education and training program for ToMMo GMRCs consisting of 17 lectures, one practical training session on the informed consent procedure, and written and interview examinations. Re-education workshops aimed at self-improvement are also carried out following certification. In this study, we evaluated the education and training program in terms of overall understanding, usefulness, and satisfaction using an anonymous questionnaire.

METHODS

An anonymous questionnaire addressing each aspect of the education and training program (understanding, usefulness, and satisfaction) was distributed among 152 qualified ToMMo GMRCs. Responses were received from 94 participants (61.8%).

RESULTS

There was a significant association between the level of overall understanding of lectures and medical qualification (nurse or clinical laboratory technologist), but not with age or educational background. The level of understanding and overall usefulness were lower in sessions related to genetics and epidemiology than those dealing with ToMMo practices. In the re-education workshops, GMRCs showed a preference for and hoped to learn more about both background knowledge and research progress in the ToMMo.

CONCLUSIONS

The results of our questionnaire suggest that not all ToMMo GMRCs are able to understand everything during the initial education and training program, especially in terms of genomic medicine. Continuous re-education is therefore vital in improving knowledge, skills and motivation, and preparing GMRCs for a specialist role in community-based personalized healthcare.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Regional genetic differences among Japanese populations and performance of genotype imputation using whole-genome reference panel of the Tohoku Medical Megabank Project

Background

Genotype imputation from single-nucleotide polymorphism (SNP) genotype data using a haplotype reference panel consisting of thousands of unrelated individuals from populations of interest can help to identify strongly associated variants in genome-wide association studies. The Tohoku Medical Megabank (TMM) project was established to support the development of precision medicine, together with the whole-genome sequencing of 1070 human genomes from individuals in the Miyagi region (Northeast Japan) and the construction of the 1070 Japanese genome reference panel (1KJPN). Here, we investigated the performance of 1KJPN for genotype imputation of Japanese samples not included in the TMM project and compared it with other population reference panels.

Results

We found that the 1KJPN population was more similar to other Japanese populations, Nagahama (south-central Japan) and Aki (Shikoku Island), than to East Asian populations in the 1000 Genomes Project other than JPT, suggesting that the large-scale collection (more than 1000) of Japanese genomes from the Miyagi region covered many of the genetic variations of Japanese in mainland Japan. Moreover, 1KJPN outperformed the phase 3 reference panel of the 1000 Genomes Project (1KGPp3) for Japanese samples, and IKJPN showed similar imputation rates for the TMM and other Japanese samples for SNPs with minor allele frequencies (MAFs) higher than 1%.

Conclusions

1KJPN covered most of the variants found in the samples from areas of the Japanese mainland outside the Miyagi region, implying 1KJPN is representative of the Japanese population’s genomes. 1KJPN and successive reference panels are useful genome reference panels for the mainland Japanese population. Importantly, the addition of whole genome sequences not included in the 1KJPN panel improved imputation efficiencies for SNPs with MAFs under 1% for samples from most regions of the Japanese archipelago.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4942-0) contains supplementary material, which is available to authorized users.

Omics research project on prospective cohort studies from the Tohoku Medical Megabank Project

Population-based prospective cohort studies are indispensable for modern medical research as they provide important knowledge on the influences of many kinds of genetic and environmental factors on the cause of disease. Although traditional cohort studies are mainly conducted using questionnaires and physical examinations, modern cohort studies incorporate omics and genomic approaches to obtain comprehensive physical information, including genetic information. Here, we report the design and midterm results of multi-omics analysis on population-based prospective cohort studies from the Tohoku Medical Megabank (TMM) Project. We have incorporated genomic and metabolomic studies in the TMM cohort study as both metabolome and genome analyses are suitable for high-throughput analysis of large-scale cohort samples. Moreover, an association study between the metabolome and genome show that metabolites are an important intermediate phenotype connecting genetic and lifestyle factors to physical and pathologic phenotypes. We apply our metabolome and genome analyses to large-scale cohort samples in the following studies.

GWAS of clinically defined gout and subtypes identifies multiple susceptibility loci that include urate transporter genes

Objective

A genome-wide association study (GWAS) of gout and its subtypes was performed to identify novel gout loci, including those that are subtype-specific.

Methods

Putative causal association signals from a GWAS of 945 clinically defined gout cases and 1213 controls from Japanese males were replicated with 1396 cases and 1268 controls using a custom chip of 1961 single nucleotide polymorphisms (SNPs). We also first conducted GWASs of gout subtypes. Replication with Caucasian and New Zealand Polynesian samples was done to further validate the loci identified in this study.

Results

In addition to the five loci we reported previously, further susceptibility loci were identified at a genome-wide significance level (p<5.0×10⁻⁸): urate transporter genes (SLC22A12 and SLC17A1) and HIST1H2BF-HIST1H4E for all gout cases, and NIPAL1 and FAM35A for the renal underexcretion gout subtype. While NIPAL1 encodes a magnesium transporter, functional analysis did not detect urate transport via NIPAL1, suggesting an indirect association with urate handling. Localisation analysis in the human kidney revealed expression of NIPAL1 and FAM35A mainly in the distal tubules, which suggests the involvement of the distal nephron in urate handling in humans. Clinically ascertained male patients with gout and controls of Caucasian and Polynesian ancestries were also genotyped, and FAM35A was associated with gout in all cases. A meta-analysis of the three populations revealed FAM35A to be associated with gout at a genome-wide level of significance (p_meta=3.58×10⁻⁸).

Conclusions

Our findings including novel gout risk loci provide further understanding of the molecular pathogenesis of gout and lead to a novel concept for the therapeutic target of gout/hyperuricaemia.

The structural origin of metabolic quantitative diversity

Relationship between structural variants of enzymes and metabolic phenotypes in human population was investigated based on the association study of metabolite quantitative traits with whole genome sequence data for 512 individuals from a population cohort. We identified five significant associations between metabolites and non-synonymous variants. Four of these non-synonymous variants are located in enzymes involved in metabolic disorders, and structural analyses of these moderate non-synonymous variants demonstrate that they are located in peripheral regions of the catalytic sites or related regulatory domains. In contrast, two individuals with larger changes of metabolite levels were also identified, and these individuals retained rare variants, which caused non-synonymous variants located near the catalytic site. These results are the first demonstrations that variant frequency, structural location, and effect for phenotype correlate with each other in human population, and imply that metabolic individuality and susceptibility for diseases may be elicited from the moderate variants and much more deleterious but rare variants.

Rare variant discovery by deep whole-genome sequencing of 1,070 Japanese individuals

The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of <1.0%. This detailed analysis detected signatures for purifying selection on regulatory elements as well as coding regions. We also catalogue structural variants, including 3.4 million insertions and deletions, and 25,923 genic copy-number variants. The 1KJPN was effective for imputing genotypes of the Japanese population genome wide. These data demonstrate the value of high-coverage sequencing for constructing population-specific variant panels, which covers 99.0% SNVs of minor allele frequency ≥0.1%, and its value for identifying causal rare variants of complex human disease phenotypes in genetic association studies.

The Tohoku Medical Megabank Organization establishes a biobank with detailed patient health care and genome information. Here the authors analyse whole-genome sequences of 1,070 Japanese individuals, allowing them to catalogue 21 million single-nucleotide variants including 12 million novel ones.

Abstract 1397: Identification of a hepatocellular carcinoma cell line capable of monitoring the differentiation of CD13(+)CD166(-) cancer stem cells and effects of sorafenib

Background: Recent cancer stem cell (CSC) research has highlighted the necessity to develop drugs targeting CSCs. However, only few in vitro models are available which show CSC hierarchy and allow monitoring of CSC differentiation and screening of drugs targeting CSC.

Aim: In the present study, we aimed to identify such hepatocellular carcinoma (HCC) cell line(s).

Methods: Based on a hypothesis that the proportion of CSCs would decrease during culture because of their differentiation into progenitors (population switch), we tested the positivity of several markers (CD13, EpCAM, CD133, CD44, CD90, CD24, CD166) in HCC cell lines (HuH-7, Li-7, PLC/PRF/5, HLF, HLE) before and after 2 months culture.

Results: Only Li-7 cells showed relevant population switch: the proportion of CD13 decreased, while that of CD166 increased during the culture. High tumorigenicity of Li-7 was lost after the population switch. CD13(+)/CD166(-) cells showed slow growth and reconstructed the bulk Li-7 populations composed of CD13(+)/CD166(-), CD13(-)/CD166(-) and CD13(-)/CD166(+) fractions, while CD13(-)/CD166(+) cells showed rapid growth but could not reproduce any other cell populations. CD13(+)/CD166(-) cells showed high ALDH activity, spheroid forming ability and resistance to 5FU. Microarray analysis demonstrated higher expression of stemness-related genes in CD166(-) than CD166(+) fraction. These results indicated the hierarchy of Li-7 cells, in which CD13(+)/CD166(-) and CD13(-)/CD166(+) cells serve as slow growing CSCs and rapid growing progenitors, respectively. Interestingly, sorafenib selectively targeted CD166(-) fraction including CD13(+) CSCs, which exhibited higher mRNA expression of FGF3 and FGF4, candidate biomarkers for sorafenib. 5-FU followed by sorafenib inhibited the growth of bulk Li-7 cells more effectively than reverse sequence or either alone.

Conclusions: We identified an useful cell line Li-7 capable of monitoring the CD13 (+) CSC differentiation and screening the drugs targeting the CSCs of HCC and found that sorafenib preferentially targeted CD13(+) CSCs in vitro. Screening cell lines for “population switch” may be useful to identify novel candidate model systems for screening drugs targeting CSCs.

Citation Format: Takeshi Yamada, Masato Abei, Inaho Danjoh, Ryoko Shirota, Taro Yamashita, Shinji Endo, Ichinosuke Hyodo, Yukio Nakamura. Identification of a hepatocellular carcinoma cell line capable of monitoring the differentiation of CD13(+)CD166(-) cancer stem cells and effects of sorafenib. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1397. doi:10.1158/1538-7445.AM2015-1397

Identification of a unique hepatocellular carcinoma line, Li-7, with CD13(+) cancer stem cells hierarchy and population change upon its differentiation during culture and effects of sorafenib

Backgrounds

Cancer stem cell (CSC) research has highlighted the necessity of developing drugs targeting CSCs. We investigated a hepatocellular carcinoma (HCC) cell line that not only has CSC hierarchy but also shows phenotypic changes (population changes) upon differentiation of CSC during culture and can be used for screening drugs targeting CSC.

Methods

Based on a hypothesis that the CSC proportion should decrease upon its differentiation into progenitors (population change), we tested HCC cell lines (HuH-7, Li-7, PLC/PRF/5, HLF, HLE) before and after 2 months culture for several markers (CD13, EpCAM, CD133, CD44, CD90, CD24, CD166). Tumorigenicity was tested using nude mice. To evaluate the CSC hierarchy, we investigated reconstructivity, proliferation, ALDH activity, spheroid formation, chemosensitivity and microarray analysis of the cell populations sorted by FACS.

Results

Only Li-7 cells showed a population change during culture: the proportion of CD13 positive cells decreased, while that of CD166 positive cells increased. The high tumorigenicity of the Li-7 was lost after the population change. CD13(+)/CD166(−) cells showed slow growth and reconstructed the bulk Li-7 populations composed of CD13(+)/CD166(−), CD13(−)/CD166(−) and CD13(−)/CD166(+) fractions, whereas CD13(−)/CD166(+) cells showed rapid growth but could not reproduce any other population. CD13(+)/CD166(−) cells showed high ALDH activity, spheroid forming ability and resistance to 5-fluorouracil. Microarray analysis demonstrated higher expression of stemness-related genes in CD166(−) than CD166(+) fraction. These results indicated a hierarchy in Li-7 cells, in which CD13(+)/CD166(−) and CD13(−)/CD166(+) cells serve as slow growing CSCs and rapid growing progenitors, respectively. Sorafenib selectively targeted the CD166(−) fraction, including CD13(+) CSCs, which exhibited higher mRNA expression for FGF3 and FGF4, candidate biomarkers for sorafenib. 5-fluorouracil followed by sorafenib inhibited the growth of bulk Li-7 cells more effectively than the reverse sequence or either alone.

Conclusions

We identified a unique HCC line, Li-7, which not only shows heterogeneity for a CD13(+) CSC hierarchy, but also undergoes a “population change” upon CSC differentiation. Sorafenib targeted the CSC in vitro, supporting the use of this model for screening drugs targeting the CSC. This type of “heterogeneous, unstable” cell line may prove more useful in the CSC era than conventional “homogeneous, stable” cell lines.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1297-7) contains supplementary material, which is available to authorized users.

Genome-wide association study of clinically defined gout identifies multiple risk loci and its association with clinical subtypes

OBJECTIVE

Gout, caused by hyperuricaemia, is a multifactorial disease. Although genome-wide association studies (GWASs) of gout have been reported, they included self-reported gout cases in which clinical information was insufficient. Therefore, the relationship between genetic variation and clinical subtypes of gout remains unclear. Here, we first performed a GWAS of clinically defined gout cases only.

METHODS

A GWAS was conducted with 945 patients with clinically defined gout and 1213 controls in a Japanese male population, followed by replication study of 1048 clinically defined cases and 1334 controls.

RESULTS

Five gout susceptibility loci were identified at the genome-wide significance level (p<5.0×10(-8)), which contained well-known urate transporter genes (ABCG2 and SLC2A9) and additional genes: rs1260326 (p=1.9×10(-12); OR=1.36) of GCKR (a gene for glucose and lipid metabolism), rs2188380 (p=1.6×10(-23); OR=1.75) of MYL2-CUX2 (genes associated with cholesterol and diabetes mellitus) and rs4073582 (p=6.4×10(-9); OR=1.66) of CNIH-2 (a gene for regulation of glutamate signalling). The latter two are identified as novel gout loci. Furthermore, among the identified single-nucleotide polymorphisms (SNPs), we demonstrated that the SNPs of ABCG2 and SLC2A9 were differentially associated with types of gout and clinical parameters underlying specific subtypes (renal underexcretion type and renal overload type). The effect of the risk allele of each SNP on clinical parameters showed significant linear relationships with the ratio of the case-control ORs for two distinct types of gout (r=0.96 [p=4.8×10(-4)] for urate clearance and r=0.96 [p=5.0×10(-4)] for urinary urate excretion).

CONCLUSIONS

Our findings provide clues to better understand the pathogenesis of gout and will be useful for development of companion diagnostics.

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population

BACKGROUND

Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously.

RESULTS

Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%.

CONCLUSIONS

Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Dominant expansion of a cryptic subclone with an abnormal karyotype in B lymphoblastoid cell lines during culture

Although B lymphoblastoid cell lines (B-LCLs) are thought to maintain their original genomic structures during long-term culture, there has been considerable disagreement on the actual genomic stability of these cells. This study was initiated to determine whether B-LCLs develop cell populations with abnormal genomes during culture and to search for factors important to the maintenance of the original genome. We established continuous cultures of B-LCLs for more than 6 months and analyzed the cells using array-based comparative genome hybridization (CGH) analysis, conventional karyotyping and analysis of V(D)J recombination in the immunoglobulin (Ig) gene. We found that one B-LCL acquired an extra chromosome 4 without any other genomic rearrangements at passage 16 of continuous culture. At the Ig light- and heavy-chain loci, analysis of the major cell population showed a difference between cultures at early and later passages. Another aneuploid line was detected among B-LCLs established elsewhere and deposited previously into the RIKEN Cell Bank. Our findings indicate that some of the genomic rearrangements in B-LCLs are not caused by gradual accumulation of mutations and rearrangements during the B-LCL establishment processes, but rather as a result of a change in the cell population from clones with a normal genome to clones with de novo rearrangements. It is therefore feasible to maintain B-LCLs with a normal genomic structure by cell cloning or similar treatment.

The Sonoda-Tajima Cell Collection: a human genetics research resource with emphasis on South American indigenous populations

The Sonoda–Tajima Cell Collection includes cell samples obtained from a range of ethnic minority groups across the world but in particular from South America. The collection is made all the more valuable by the fact that some of these ethnic populations have since died out, and thus it will be impossible to prepare a similar cell collection again. The collection was donated to our institute, a public cell bank in Japan, by Drs Sonoda and Tajima to make it available to researchers throughout the world. The original cell collection was composed of cryopreserved peripheral blood samples that would obviously have been rapidly exhausted if used directly. We, therefore, immortalized some samples with the Epstein–Barr virus and established B-lymphoblastoid cell lines (B-LCLs). As there is continuing controversy over whether the B-LCL genome is stably maintained, we performed an array comparative genomic hybridization (CGH) analysis to confirm the genomic stability of the cell lines. The array CGH analysis of the B-LCL lines and their parental B cells demonstrated that genomic stability was maintained in the long-term cell cultures. The B-LCLs of the Sonoda–Tajima Collection will therefore be made available to interested scientists around the world. At present, 512 B-LCLs have been developed, and we are willing to increase the number if there is sufficient demand.

Is parainfluenza virus a threatening virus for human cancer cell lines?

Immortalized cell lines, such as human cancer cell lines, are an indispensable experimental resource for many types of biological and medical research. However, unless the cell line has been authenticated prior to use, interpretation of experimental results may be problematic. The potential problems this may cause are illustrated by studies in which authentication of cell lines has not been carried out. For example, immortalized cell lines may unknowingly be infected with viruses that alter their characteristics. In fact, parainfluenza virus type 5 (PIV5) poses a threat to the use of immortalized cell lines in biological and medical research; PIV5 infection significantly alters cellular physiology associated with the response to interferon. If PIV5 infection is widespread in immortalized cell lines, then a very large number of published studies might have to be re-evaluated. Fortunately, analyses of a large number of immortalized cell lines indicate that PIV5 infection is not widespread.

IL-6-induced Bcl6 variant 2 supports IL-6-dependent myeloma cell proliferation and survival through STAT3

IL-6 is a growth and survival factor for myeloma cells, although the mechanism by which it induces myeloma cell proliferation through gene expression is largely unknown. Microarray analysis showed that some B-cell lymphoma-associated oncogenes such as Bcl6, which is absent in normal plasma cells, were upregulated by IL-6 in IL-6-dependent myeloma cell lines. We found that Bcl6 variant 2 was upregulated by STAT3. ChIP assay and EMSA showed that STAT3 bound to the upstream region of variant 2 DNA. Expression of p53, a direct target gene of Bcl6, was downregulated in the IL-6-stimulated cells, and this process was impaired by an HDAC inhibitor. Bcl6 was knocked down by introducing small hairpin RNA, resulting in decreased proliferation and increased sensitivity to a DNA damaging agent. Thus, STAT3-inducible Bcl6 variant 2 appears to generate an important IL-6 signal that supports proliferation and survival of IL-6-dependent myeloma cells.

Discovery of aberrant expression of R-RAS by cancer-linked DNA hypomethylation in gastric cancer using microarrays

Although hypomethylation was the originally identified epigenetic change in cancer, it was overlooked for many years in preference to hypermethylation. Recently, gene activation by cancer-linked hypomethylation has been rediscovered. However, in gastric cancer, genome-wide screening of the activated genes has not been found. By using microarrays, we identified 1,383 gene candidates reactivated in at least one cell line of eight gastric cancer cell lines after treatment with 5-aza-2'deoxycytidine and trichostatin A. Of the 1,383 genes, 159 genes, including oncogenes ELK1, FRAT2, R-RAS, RHOB, and RHO6, were further selected as gene candidates that are silenced by DNA methylation in normal stomach mucosa but are activated by DNA demethylation in a subset of gastric cancers. Next, we showed that demethylation of specific CpG sites within the first intron of R-RAS causes activation in more than half of gastric cancers. Introduction of siRNA into R-RAS-expressing cells resulted in the disappearance of the adhered cells, suggesting that functional blocking of the R-RAS-signaling pathway has great potential for gastric cancer therapy. Our extensive gene list provides other candidates for this class of oncogene.

Highly specific marker genes for detecting minimal gastric cancer cells in cytology negative peritoneal washings

Peritoneal wash cytology plays a pivotal role in the decision for gastric cancer treatment because advanced gastric cancer often turns out incurable with peritoneal metastasis. Molecular detection of minimal cancer cells from peritoneal washings may overcome the sensitivity boundary of conventional cytology and contribute to the prediction of the disease outcome. To select marker candidates out of ten thousands of genes, we performed microarray analyses in 12 gastric cell lines and 8 peritoneal washings of early stage cases. With 40 candidates selected by the above expression profiling, RT-PCR in 16 representative peritoneal wash samples was performed to identify genes specific to cytology positive samples. The finally selected five genes, CK20, FABP1, MUC2, TFF1, and TFF2, were then evaluated for their utility as a marker for minimal residual disease in 99 peritoneal wash samples. Nested RT-PCR using the five genes showed positive results highly specific to incurable cases (91-100%). With a high specificity, the combination of these five genes succeeded in identifying 6 out of 20 (30%) additional patients with all types of early recurrence that could not be predicted by the conventional method. The six newly identified recurrences included four non-peritoneal ones, showing that RT-PCR using the five genes without a real-time quantitative PCR technique contributes to the detection of minimal residual disease.

Ras-mediated signaling pathway regulates the expression of a low-molecular-weight heat-shock protein in fission yeast

In fission yeast, Schizosaccharomyces pombe, deficiency of ras1 gene causes an abnormal cell shape and abolishes mating ability. However, target genes of this signaling pathway are largely unknown because of the lack of an appropriate analysis system. To overcome this problem, we have started a novel project to categorize entire genes based on their expression levels under different growth conditions. Using this strategy, we screened genes whose expression levels were affected in the presence or absence of the ras1 gene product. For this purpose, we utilized high-density arrays of clones covering the entire genome of the fission yeast, and probed with labelled cDNA derived from various strains and growth conditions. Here, we demonstrate the detection of a low-molecular-weight heat-shock protein gene, hsp16, whose expression is very likely to be regulated by a ras-mediated signaling pathway, but not by the heat-shock response.

Identification and characterization of a novel trans-membrane protein gene, pdh1, from Schizosaccharomyces pombe

We have cloned a new gene, pdh1, from genomic DNA of fission yeast Schizosaccharomyces pombe. pdh1 is actively transcribed as 1400-nucleotide mRNA in vegetatively growing cells and can code for a 226 amino acid polypeptide (pdh1p). Computational structural prediction has revealed that the pdh1p is a highly hydrophobic protein with seven transmembrane domains. The prediction has also detected a possible C-kinase phosphorylation site within the longest hydrophilic loop.

Enzymic characterization of fission yeast farnesyl transferase: recognition of the -CAAL motif at the C-terminus

The enzyme farnesyl transferase (FTase) catalyzes the posttranslational modification of Ras and other Ras family proteins with a C15 farnesyl group. The target proteins have a consensus -CAAX motif (X, any amino acid except leucine) at the C-terminus. Since proteins that have leucine as the C-terminal amino acid X are modified with a C20 geranylgeranyl group, it is thought that the C-terminal leucine is the signal (-CAAL motif) for selection of isoprenoid molecules. Here, we report the presence of multiple FTase activities in the fission yeast Schizosaccharomyces pombe, each seeming to correspond to a particular protein known to be modified by the farnesyl group in vivo. Using enzymic activities specific to S. pombe Ras1, we found similar affinities for FTases in the wild-type (EVSTKCCVIC) and mutant Ras1 peptide, in which the C-terminal amino acid is replaced by leucine (EVSTKCCVIL). These results suggest that recognition and selection of the correct isoprenoid group by the FTases require other amino acid sequences of the target protein in addition to the C-terminal -CAAX motif.

Biochemical similarity of Schizosaccharomyces pombe ras1 protein with RAS2 protein of Saccharomyces cervisiae

Schizosaccharomyces pombe contains single ras oncogene homologue, ras1, that functions in the signal transduction pathway conducting the cell's mating processes. To understand the biochemical basis of yeast ras proteins, we have purified the ras1 protein and compared the major biochemical constants with those of RAS2 protein from Saccharomyces cerevisiae and mammalian ras proteins. The purified ras1 protein showed a remarkably high Kd value for GDP binding (178 nM) and for binding with ATP. In contrast, the Kd value for GTP binding and the rate of GTPase activity were 64 nM and 77 x 10(-6) s-1 at 37 degrees C, respectively; both were higher than normal p21ras protein, but at the same level as the RAS2 protein. We directly measured rate of GTP binding and GDP binding which were 3.9 x 10(-3) s-1 and 1.8 x 10(-3) s-1 at 30 degrees C, respectively. On the other hand, exchange rates between bound and free nucleotides remained almost constant throughout the tested combination of GTP and GDP, and were several-fold lower than the binding rate. These results suggest that the release of the guanine nucleotide is the rate-limiting step in the ras-GTP/GDP cycle. As a whole, the biochemical properties of the ras1 protein are close to those of the RAS2 protein, although these two proteins function differently in the signal transduction pathway in the cells.