Landscape of Actionable Genetic Alterations Profiled from 1,071 Tumor Samples in Korean Cancer Patients / Journal of the Korean Cancer Association, 대한암학회지

PURPOSE: With the emergence of next-generation sequencing (NGS) technology, profiling a wide range of genomic alterations has become a possibility resulting in improved implementation of targeted cancer therapy. In Asian populations, the prevalence and spectrum of clinically actionable genetic alterations has not yet been determined because of a lack of studies examining high-throughput cancer genomic data. MATERIALS AND METHODS: To address this issue, 1,071 tumor samples were collected from five major cancer institutes in Korea and analyzed using targeted NGS at a centralized laboratory. Samples were either fresh frozen or formalin-fixed, paraffin embedded (FFPE) and the quality and yield of extracted genomic DNA was assessed. In order to estimate the effect of sample condition on the quality of sequencing results, tissue preparation method, specimen type (resected or biopsied) and tissue storage time were compared. RESULTS: We detected 7,360 non-synonymous point mutations, 1,164 small insertions and deletions, 3,173 copy number alterations, and 462 structural variants. Fifty-four percent of tumors had one or more clinically relevant genetic mutation. The distribution of actionable variants was variable among different genes. Fresh frozen tissues, surgically resected specimens, and recently obtained specimens generated superior sequencing results over FFPE tissues, biopsied specimens, and tissues with long storage duration. CONCLUSION: In order to overcome, challenges involved in bringing NGS testing into routine clinical use, a centralized laboratory model was designed that could improve the NGS workflows, provide appropriate turnaround times and control costs with goal of enabling precision medicine.

The Effect of Probiotics, Antibiotics, and Antipyretic Analgesics on Gut Microbiota Modification

Human gut microbial community is playing a critical role in human health and associated with different human disease. In parallel, probiotics, antibiotics, and antipyretic analgesics (AAs) were developed to improve human health or cure human diseases. We therefore examined how probiotics, antibiotics, and AAs influence to the gut microbiota. Three independent case/control studies were designed from the cross-sectional cohort data of 1,463 healthy Koreans. The composition of the gut microbiota in each case and control group was determined via 16S ribosomal RNA Illumina next-generation sequencing. The correlation between microbial taxa and the consumption of each drug was tested using zero-inflated Gaussian mixture models, with covariate adjustment of age, sex, and body mass index (BMI). Probiotics, antibiotics, and AAs consumption yielded the significant differences in the gut microbiota, represented the lower abundance of Megasphaera in probiotics, the higher abundance of Fusobacteria in antibiotics, and the higher abundance of Butyrivibrio and Verrucomicrobia in AAs, compared to each control group. The reduction of Erysipelotrichaceae family was common in three drugs consumption.

Pathway Analysis of Metabolic Syndrome Using a Genome-Wide Association Study of Korea Associated Resource (KARE) Cohorts

Metabolic syndrome (MetS) is a complex disorder related to insulin resistance, obesity, and inflammation. Genetic and environmental factors also contribute to the development of MetS, and through genome-wide association studies (GWASs), important susceptibility loci have been identified. However, GWASs focus more on individual single-nucleotide polymorphisms (SNPs), explaining only a small portion of genetic heritability. To overcome this limitation, pathway analyses are being applied to GWAS datasets. The aim of this study is to elucidate the biological pathways involved in the pathogenesis of MetS through pathway analysis. Cohort data from the Korea Associated Resource (KARE) was used for analysis, which include 8,842 individuals (age, 52.2 +/- 8.9 years; body mass index, 24.6 +/- 3.2 kg/m2). A total of 312,121 autosomal SNPs were obtained after quality control. Pathway analysis was conducted using Meta-analysis Gene-Set Enrichment of Variant Associations (MAGENTA) to discover the biological pathways associated with MetS. In the discovery phase, SNPs from chromosome 12, including rs11066280, rs2074356, and rs12229654, were associated with MetS (p < 5 x 10(-6)), and rs11066280 satisfied the Bonferroni-corrected cutoff (unadjusted p < 1.38 x 10(-7), Bonferroni-adjusted p < 0.05). Through pathway analysis, biological pathways, including electron carrier activity, signaling by platelet-derived growth factor (PDGF), the mitogen-activated protein kinase kinase kinase cascade, PDGF binding, peroxisome proliferator-activated receptor (PPAR) signaling, and DNA repair, were associated with MetS. Through pathway analysis of MetS, pathways related with PDGF, mitogen-activated protein kinase, and PPAR signaling, as well as nucleic acid binding, protein secretion, and DNA repair, were identified. Further studies will be needed to clarify the genetic pathogenesis leading to MetS.

Trends and Directions in Personality Genetic Studies

How personality forms and whether personality genes exist are long-studied questions. Various concepts and theories have been presented for centuries. Personality is a complex trait and is developed through the interaction of genes and the environment. Twin and family studies have found that there are critical genetic and environmental components in the inheritance of personality traits, and modern advances in genetics are making it possible to identify specific variants for personality traits. Although genes that were found in studies on personality have not provided replicable association between genetic and personality variability, more and more genetic variants associated with personality traits are being discovered. Here, we present the current state of the art on genetic research in the personality field and finally list several of the recently published research highlights. First, we briefly describe the commonly used self-reported measures that define personality traits. Then, we summarize the characteristics of the candidate genes for personality traits and investigate gene variants that have been suggested to be associated with personality traits.

Genome-wide scan of granular corneal dystrophy, type II: confirmation of chromosome 5q31 and identification of new co-segregated loci on chromosome 3q26.3

Granular corneal dystrophy, type II (CGD2; Avellino corneal dystrophy) is the most common corneal dystrophy among Koreans, but its pathophysiology is still poorly understood. Many reports showed that even though the causative mutation is the same TGFBI R124H mutation, there are severe and mild phenotypes of the corneal dystrophy. We also observed the phenotype differences in our samples. For this reason, we focused our effort on the identification of unknown genetic factor related to phenotype variation. A total 551 individuals from 59 families were genotyped with SNP chip and used in genome-wide linkage analysis. From single-point linkage analyses, we confirmed the known 5q31 region for TGFBI gene, and selected novel nine candidate loci for CGD2. In simulation analysis, the only 3q26.3 region including neuroligin 1 gene (NLGN1) was supported by empirical statistic significance. To investigate the effect of genetic heterogeneity in linkage analysis, we classified CGD2 families into two subgroups. Although we could not find a significant evidence for correlation between the 3q26.3 region and CGD2 phenotypes, this first genome-wide analysis with CGD2 families in Korea has a very important value for offering insights in genetics of CGD2. In addition, the co-segregating loci with CGD2 including 3q26.3 would be a good target for further study to understand the pathophysiology of CGD2.

Genome-Wide Association Studies of the Korea Association REsource (KARE) Consortium

During the last decade, large community cohorts have been established by the Korea National Institutes of Health (KNIH), and enormous epidemiological and clinical data have been accumulated. Using these information and samples in the cohorts, KNIH set out to do a large-scale genome-wide association study (GWAS) in 2007, and the Korea Association REsource (KARE) consortium was launched to analyze the data to identify the underlying genetic risk factors of diseases and diverse health indexes, such as blood pressure, obesity, bone density, and blood biochemical traits. The consortium consisted of 6 research divisions, formed by 25 principal investigators in 19 organizations, including 18 universities, 2 institutes, and 1 company. Each division focused on one of the following subjects: the identification of genetic factors, the statistical analysis of gene-gene interactions, the genetic epidemiology of gene-environment interactions, copy number variation, the bioinformatics related to a GWAS, and a GWAS of nutrigenomics. In this special issue, the study results of the KARE consortium are provided as 9 articles. We hope that this special issue might encourage the genomics community to share data and scientists, including clinicians, to analyze the valuable Korean data of KARE.

Analysis of genetic and non-genetic factors that affect the QTc interval in a Mongolian population: the GENDISCAN study

The QTc interval is a complex quantitative trait and a strong prognostic indicator of cardiovascular mortality in general, healthy people. The aim of this study was to identify non-genetic factors and quantitative trait loci that govern the QTc interval in an isolated Mongolian population. We used multiple regression analysis to determine the relationship between the QTc interval and non-genetic factors including height, blood pressure, and the plasma lipid level. Whole genome linkage analyses were performed to reveal quantitative trait loci for the QTc interval with 349 microsatellite markers from 1,080 Mongolian subjects. Among many factors previously known for association with the QTc interval, age, sex, heart rate, QRS duration of electrocardiogram and systolic blood pressure were also found to have influence on the QTc interval. A genetic effect for the QTc interval was identified based on familial correlation with a heritability value of 0.31. In a whole genome linkage analysis, we identified the four potential linkage regions 7q31-34, 5q21, 4q28, and 2q36.

Erratum

The authors regret an error in discussion, the authors wrote that "We also performed an additional linkage analysis using only the adult population (> or = 20 years old)(Table 7, Figure 2 and 3)." In this sentence, Table 7 should be changed to Table 6.

Socioeconomic burden of the polycystic ovary syndrome / 대한산부인과학회지

OBJECTIVE: The objective of the study was to estimate socioeconomic burden of polycystic ovary syndrome (PCOS) during the reproductive life span using current definitions and prevalence or incidence data. METHODS: Questionnaires were given to 8,588 reproductive women reviewed at Ewha Womans University Mokdong hospital. The PCOS affected approximately 10.4% of reproductive-aged women (11 million women in Korea, prevalence rate according to 1990 National Institutes of Health PCOS diagnosis criteria). We tied general societal cost data for the different health consequences to reproductive-age PCOS costs, using prevalence data. RESULTS: We estimated the mean annual cost of the initial evaluation to be 76 hundred million won, that of hormonally treating menstrual dysfunction, providing infertility care, diagnosis/treatment of endometrial hyperplasia, GDM, type 2 DM, and hypertension to be 280 billion won. The total annual socioeconomic cost of evaluating and providing care to reproductive-aged PCOS women in Korea is 350 billion won. CONCLUSION: Because the cost of the diagnostic evaluation accounted for a relatively minor part of the total socioeconomic costs, more widespread screening for PCOS appears be a cost-effective strategy, leading to earlier diagnosis and intervention and possibly the amelioration and prevention of serious sequelae.

A comprehensive profile of DNA copy number variations in a Korean population: identification of copy number invariant regions among Koreans

To examine copy number variations among the Korean population, we compared individual genomes with the Korean reference genome assembly using the publicly available Korean HapMap SNP 50 k chip data from 90 individuals. Korean individuals exhibited 123 copy number variation regions (CNVRs) covering 27.2 mb, equivalent to 1.0% of the genome in the copy number variation (CNV) analysis using the combined criteria of P value (P or = 0.25) among study subjects. In contrast, when compared to the Affymetrix reference genome assembly from multiple ethnic groups, considerably more CNVRs (n = 643) were detected in larger proportions (5.0%) of the genome covering 135.1 mb even by more stringent criteria (P or = 0.25), reflecting ethnic diversity of structural variations between Korean and other populations. Some CNVRs were validated by the quantitative multiplex PCR of short fluorescent fragment (QMPSF) method, and then copy number invariant regions were detected among the study subjects. These copy number invariant regions would be used as good internal controls for further CNV studies. Lastly, we demonstrated that the CNV information could stratify even a single ethnic population with a proper reference genome assembly from multiple heterogeneous populations.

Association of polymorphisms in the Interleukin 23 receptor gene with osteonecrosis of femoral head in Korean population

Osteonecrosis of the femoral head (ONFH) is known as death of the cellular portion of the femoral head due to an interruption in the vascular supply. The underlying pathophysiology regarding bone cell death remains uncertain. Recently, several studies have shown that autoimmune disorders were related to the development of osteonecrosis. This study investigated the genetic effects of Interleukin 23 receptor (IL23R) polymorphisms regarding the risk of ONFH. Ten SNPs were selected and genotyped in 443 ONFH patients and 273 control subjects in order to perform the genetic association analysis. It was found that polymorphisms of the IL23R gene (rs4655686, rs1569922 and rs7539625) were significantly associated with an increased risk of ONFH (P values; 0.0198-0.0447, OR; 1.30-1.49). Particularly, a stratified analysis based on etiology (alcohol, steroid or idiopathic) showed that the associations between these polymorphisms and ONFH were most significant in idiopathic ONFH patients (P values; 0.0001-0.0150, OR; 1.45-2.17). These results suggest that IL23R polymorphisms may play an important role in the development of ONFH.

Prevalence of Atrial Fibrillation in Middle-Aged People in Korea: The Korean Genome and Epidemiology Study

BACKGROUND AND OBJECTIVES: The prevalence of atrial fibrillation (AF), the most common sustained arrhythmia, is expected to rise with the aging population, but very few studies have reported on the prevalence and risk factors of AF in Korea. SUBJECTS AND METHODS: We analyzed 10,012 Korean adults (4,750 men and 5,262 women), 40-69 years old, who were enrolled in the Korean Genome and Epidemiology Study. AF was diagnosed by single electrocardiogram recording in a baseline survey (2001-2003). RESULTS: The estimated prevalence of AF was 0.4% {95% confidence interval (CI), 0.28-0.52} in adults 40-69 years old, and increased to 1.0% in individuals 60-69 years old. The prevalence rate for men (0.6%) was higher than for women (0.2%) across all age groups. In multiple logistic regression analysis, AF was significantly associated with old age {odds ratio (OR), 8.15; 95% CI, 3.06-21.71}, male gender (OR, 4.04; 95% CI, 1.90-8.61), diabetes mellitus (OR, 2.15; 95% CI, 1.05-4.44), and congestive heart failure (OR, 14.11; 95% CI, 2.56-77.70). Obesity, however, did not show an association with AF. CONCLUSION: The prevalence of AF in Korean adults aged 40-69 years is approximately 0.4%, lower than that in Western populations. Age, male gender, diabetes, and heart failure are associated with AF. Further research in a larger population is necessary to verify for our results.

Genome-wide Linkage Study for Plasma HDL Cholesterol Level in an Isolated Population of Mongolia

High-density lipoprotein (HDL) whose primary role is to transport cholesterol from peripheral tissues to the liver, is associated with the incidence of coronary heart disease. We analyzed HDL cholesterol levels in a genetically isolated population of extended Mongolian families. A total of 1002 individuals (54.5% women) from 95 families were enrolled. After genotyping by use of 1000 microsatellite markers, we performed a genome-wide linkage search with variance component analysis. The estimated heritability of HDL cholesterol was 0.45, revealing that HDL cholesterol was under significant genetic influence. We found peak evidence of linkage (LOD score=1.88) for HDL cholesterol level on chromosome 6(nearest marker D6S1660) and potential evidences for linkage on chromosomes 1, 12 and 19 with the LOD scores of 1.32, 1.44 and 1.14, respectively. These results should pave the way for the discovery of the relevant genes by fine mapping and association analysis.

Genome-wide Linkage Study for Plasma HDL Cholesterol Level in an Isolated Population of Mongolia

High-density lipoprotein (HDL) whose primary role is to transport cholesterol from peripheral tissues to the liver, is associated with the incidence of coronary heart disease. We analyzed HDL cholesterol levels in a genetically isolated population of extended Mongolian families. A total of 1002 individuals (54.5% women) from 95 families were enrolled. After genotyping by use of 1000 microsatellite markers, we performed a genome-wide linkage search with variance component analysis. The estimated heritability of HDL cholesterol was 0.45, revealing that HDL cholesterol was under significant genetic influence. We found peak evidence of linkage (LOD score=1.88) for HDL cholesterol level on chromosome 6(nearest marker D6S1660) and potential evidences for linkage on chromosomes 1, 12 and 19 with the LOD scores of 1.32, 1.44 and 1.14, respectively. These results should pave the way for the discovery of the relevant genes by fine mapping and association analysis.

Heterogeneous Composition of Voltage-Dependent K+ Currents in Hepatic Stellate Cells

PURPOSE: Hepatic stellate cells (HSC) are a type of pericyte with varying characteristics according to their location. However, the electrophysiological properties of HSC are not completely understood. Therefore, this study investigated the difference in the voltage-dependent K(+) currents in HSC. MATERIALS AND METHODS: The voltage-dependent K(+) currents in rat HSC were evaluated using the whole cell configuration of the patch-clamp technique. RESULTS: Four different types of voltage-dependent K(+) currents in HSC were identified based on the outward and inward K(+) currents. Type D had the dominant delayed rectifier K(+) current, and type A had the dominant transient outward K(+) current. Type I had an inwardly rectifying K(+) current, whereas the non-type I did not. TEA (5mM) and 4-AP (2mM) suppressed the outward K(+) currents differentially in type D and A. Changing the holding potential from -80 to -40mV reduced the amplitude of the transient outward K(+) currents in type A. The inwardly rectifying K(+) currents either declined markedly or were sustained in type I during the hyperpolarizing step pulses from -120 to -150mV. CONCLUSION: There are four different configurations of voltage-dependent K(+) currents expressed in cultured HSC. These results are expected to provide information that will help determine the properties of the K(+) currents in HSC as well as the different type HSC populations.

FCAnalyzer: A Functional Clustering Analysis Tool for Predicted Transcription Regulatory Elements and Gene Ontology Terms

Numerous studies have reported that genes with similar expression patterns are co-regulated. From gene expression data, we have assumed that genes having similar expression pattern would share similar transcription factor binding sites (TFBSs). These function as the binding regions for transcription factors (TFs) and thereby regulate gene expression. In this context, various analysis tools have been developed. However, they have shortcomings in the combined analysis of expression patterns and significant TFBSs and in the functional analysis of target genes of significantly overrepresented putative regulators. In this study, we present a web-based A Functional Clustering Analysis Tool for Predicted Transcription Regulatory Elements and Gene Ontology Terms (FCAnalyzer). This system integrates microarray clustering data with similar expression patterns, and TFBS data in each cluster. FCAnalyzer is designed to perform two independent clustering procedures. The first process clusters gene expression profiles using the K-means clustering method, and the second process clusters predicted TFBSs in the upstream region of previously clustered genes using the hierarchical biclustering method for simultaneous grouping of genes and samples. This system offers retrieved information for predicted TFBSs in each cluster using Match(TM) in the TRANSFAC database. We used gene ontology term analysis for functional annotation of genes in the same cluster. We also provide the user with a combinatorial TFBS analysis of TFBS pairs. The enrichment of TFBS analysis and GO term analysis is statistically by the calculation of P values based on Fisher's exact test, hypergeometric distribution and Bonferroni correction. FCAnalyzer is a web-based, user-friendly functional clustering analysis system that facilitates the transcriptional regulatory analysis of co-expressed genes. This system presents the analyses of clustered genes, significant TFBSs, significantly enriched TFBS combinations, their target genes and TFBS-TF pairs.

Gene Expression Profiling of 6-MP (6-mercaptopurine) in Liver

The KFDA (Korea Food & Drug Administration) has performed a collaborative toxicogenomics project since 2003. Its aim is to construct a toxicology database of 12 compounds administered to mice at initial phase. We chose 6-MP (6-mercaptopurine) which has been used in the treatment of childhood leukemia. It was administered at low (0.224 mg/kg) and at high (2.24 mg/kg) dose (5 mice per group) intraperitonealy to the postnatal 6 weeks mice, then the serum and liver were collected at the indicated time (6, 24 and 72 h) after scarification. Serum biochemical markers for liver toxicity were measured and histopathologic studies also were carried out. The gene expression profiling was carried out by using Applied Biosystems 1700 Full Genome Expression Mouse. By self-organization maps (SOM), we identified groups with unique gene expression patterns, some of them are supposed to be related to 6-MP induced toxicity, including lipid metabolism abnormality, inflammatory response, oxidative stress, ATP depletion and cell death. The potential toxic effects appearing as gene expression changes are dependent of the time of 6-MP but independent of the dosage of it. This study would contribute to establishment of international database as well as national one about hepatotoxicity.

Gene Expression Analysis of Megakaryocytes Derived from Human Umbilical Cord CD34(+) Cells by Thrombopoietin

Although much is known about the molecular biology of platelets, the megakaryocytes' (MKs) molecular biology was not understood so well because of their rareness. By the cloning and characterization of thrombopoietin (TPO), which is the principal regulator of the growth and development of the MKs, researches on the MKs have been growing rapidly. To understand megakaryocytopoiesis, we investigated the gene expression profile of the MKs using oligonucleotide microarray where 10, 108 unique genes were spotted. Comparing the fluorescence intensities of which ratio is > or = |2|, 372 genes were up-regulated and 541 genes were down-regulated in MKs. For confirmatory expression, RNase protection assay (RPA) establishing abundant apoptotic gene expression was carried out. In MKs, many of the known genes, including several platelet related genes, GATA binding protein were highly expressed. Particularly, TGF beta, clusterin (complement lysis inhibitor), and thymosin beta 4 (actin-sequestering molecules) were expressed highly in MKs. As MKs specific expressed genes may regulate normal and pathologic platelet (and/or MK) functions, the transcript profiling using microarray was useful on molecular understanding of MKs.

Identification of the Marker-Genes for Dioxin(2,3,7,8-tetradibenzo-p-dioxin)-Induced Immune Dysfunction by Using the High-Density Oligonucleotide Microarray

In a variety of animal species, the perinatal exposure of experimental animals to the 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) leads to the immune dysfunction, which is more severe and persistent than that caused by adult exposure. We report here the changes of gene expression and the identification of the marker-genes representing the dioxin exposure. The expressions of the transcripts were analyzed using the 11K oligonucleotide- microarray from the bone marrow cells of male C57BL/6J mice after an intraperitoneal injection of 1 microgram TCDD/kg body weight at various time intervals: gestational 6.5 day(G6.5), 13.5 day(G13.5), 18.5 day(G18.5), and postnatal 3 (P3W)and 6 week (P6W). The type of self-organizing maps(SOM) representing the specific exposure dioxin could be identified as follows; G6.5D(C14), G13.5D(C0, C5, C10, C18), G18.5D(7), P3W(C2, C21), and P6W(C4, C15, C20). The candidate marker-genes were restricted to the transcripts, which could be consistently expressed greater than +/-2-fold in three experiments. The resulting candidates were 85 genes, the characteristics of that were involved in cell physiology and cell functions such as cell proliferation and immune function. We identified the biomarker-genes for dioxin exposure: smc -like 2 from SOM C14 for the dioxin exposure at G6.5D, focal adhesion kinase and 6 other genes from C0, and protein tyrosine phosphatase 4a2 and 3 other genes from C5 for G13.5D, platelet factor 4 from C7 for G18.5D, fos from C2 for P3W.

Comparison of oligonucleotide-microarray and serial analysis of gene expression (SAGE) in transcript profiling analysis of megakaryocytes derived from CD34+ cells

For the comprehensive analysis of transcript expression, the array-based hybridization analysis and the serial analysis of gene expression (SAGE) are commonly used platforms. The SAGE is based on a high-throughput sequencing of ditags derived from the transcript. DNA microarrays are a powerful tool for monitoring thousands of transcripts simultaneously, whereas the Genechip (Affimatrix microarray) technology is based on the hybridization of a single probe or other manufacturer's microarrays (cDNA- or oligonucleotide-microarray) procedures include the competitive hybridization of two probes. In this study, the quantitative accuracy of expression using oligonucleotide-microarray was determined by comparing data set from the SAGE. In previous study the microSAGE was performed for the megakaryocytes and non- megakaryocytes derived from human cord blood CD34(+)cells by ex vivo expansion using thrombopoietin, and a total of 38,909 tags representing 8,976 unique genes were obtained. On the identical RNA, expression profiling was also carried out using oligonucleotide-microarray (MAGIC II 10K chip, Macrogen). The most frequently expressed genes in human megakaryocytes were identified as platelet factor 1 followed by annexin A1, ribosomal protein S23. The majority of the 50 most highly expressed genes in the CD34(+)-derived megakaryocytes were those involved in protein synthesis, e.g., ribosomal proteins. The expression level through the single channel of oligonucleotide-microarray and SAGE have a fairly good correlation in terms of absolute analyses and that the correlation is higher for the genes with higher expression levels.