Cancer-specific functional profiling in microsatellite-unstable (MSI) colon and endometrial cancers using combined differentially expressed genes and biclustering analysis.

Microsatellite-unstable (MSI) cancers have distinct genetic and clinical features from microsatellite-stable cancers, but the molecular functional differences between MSI cancers originating from different tissues or organs have not been well studied because the application of usual differentially expressed gene (DEG) analysis is error-prone, producing too many noncancer-specific normally functioning genes. To maximize therapeutic efficacy, biomarkers reflecting cancer-specific differences between MSI cancers of different tissue origins should be identified. To identify functional differences between MSI colon and endometrial cancers, we combined DEG analysis and biclustering instead of DEG analysis alone and refined functionally relevant biclusters reflecting genuine functional differences between the 2 tumors. Specifically, using The Cancer Genome Atlas and genome-tissue expression as data sources, gene ontology (GO) enrichment tests were performed after routinely identifying DEGs between the 2 tumors with the exclusion of DEGs identified in their normal counterparts. Cancer-specific biclusters and associated enriched GO terms were obtained by biclustering with enrichment tests for the preferences for cancer type (either colon or endometrium) and GO enrichment tests for each cancer-specific bicluster, respectively. A novel childness score was developed to select functionally relevant biclusters among cancer-specific biclusters based on the extent to which the enriched GO terms of the biclusters tended to be child terms of the enriched GO terms in DEGs. The selected biclusters were tested using survival analysis to validate their clinical significance. We performed multiple sequential analyses to produce functionally relevant biclusters from the RNA sequencing data of MSI colon and endometrial cancer samples and their normal counterparts. We identified 3066 cancer-specific DEGs. Biclustering analysis revealed 153 biclusters and 41 cancer-specific biclusters were selected using Fisher exact test. A mean childness score over 0.6 was applied as the threshold and yielded 8 functionally relevant biclusters from cancer-specific biclusters. Functional differences appear to include gland cavitation and the TGF-ß receptor, G protein, and cytokine pathways. In the survival analysis, 6 of the 8 functionally relevant biclusters were statistically significant. By attenuating noise and applying a synergistic contribution of DEG results, we refined candidate biomarkers to complement tissue-specific features of MSI tumors.

Exome sequencing identified five novel USH2A variants in Korean patients with retinitis pigmentosa.

BACKGROUND: Retinitis pigmentosa (RP) is an inherited disorder that causes progressive loss of vision. This study aimed to describe the possible causative variants of the USH2A gene in Korean RP families and their associated phenotypes. MATERIALS AND METHODS: We recruited 94 RP families (220 subjects, including 94 probands and 126 family members) in a Korean cohort, and analyzed USH2A gene variants through whole-exome sequencing. The pathogenicity of the variants was classified according to American College of Medical Genetics and Genomics and Association for Molecular Pathology guidelines. RESULTS: We found 14 USH2A disease-causing variants, including 5 novel variants. Disease causing variants were identified in 10 probands with RP, accounting for 10.6% (10/94) of the Korean RPs in the cohort. To visually represent the structural changes induced by novel variants, we modeled the three-dimensional structures of the wild-type and mutant proteins. CONCLUSIONS: This study expands the spectrum of USH2A variants and provides information for future therapeutic strategies for RP.

Brain Region-Dependent Alternative Splicing of Alzheimer Disease (AD)-Risk Genes Is Associated With Neuropathological Features in AD.

PURPOSE: Alzheimer disease (AD) is one of the most complex diseases and is characterized by AD-related neuropathological features, including accumulation of amyloid-ß plaques and tau neurofibrillary tangles. Dysregulation of alternative splicing (AS) contributes to these features, and there is heterogeneity in features across brain regions between AD patients, leading to different severity and progression rates; however, brain region-specific AS mechanisms still remain unclear. Therefore, we aimed to systemically investigate AS in multiple brain regions of AD patients and how they affect clinical features. METHODS: We analyzed RNA sequencing (RNA-Seq) data obtained from brain regions (frontal and temporal) of AD patients. Reads were mapped to the hg19 reference genome using the STAR aligner, and exon skipping (ES) rates were estimated as percent spliced in (PSI) by rMATs. We focused on AD-risk genes discovered by genome-wide association studies, and accordingly evaluated associations between PSI of skipped exons in AD-risk genes and Braak stage and plaque density mean (PM) for each brain region. We also integrated whole-genome sequencing data of the ascertained samples with RNA-Seq data to identify genetic regulators of feature-associated ES. RESULTS: We identified 26 and 41 ES associated with Braak stage in frontal and temporal regions, respectively, and 10 and 50 ES associated with PM. Among those, 10 were frontal-specific (CLU and NTRK2), 65 temporal-specific (HIF1A and TRPC4AP), and 26 shared ES (APP) that accompanied functional Gene Ontology terms, including axonogenesis in shared-ES genes. We further identified genetic regulators that account for 44 ES (44% of the total). Finally, we present as a case study the systematic regulation of an ES in APP, which is important in AD pathogenesis. CONCLUSION: This study provides new insights into brain region-dependent AS regulation of the architecture of AD-risk genes that contributes to AD pathologies, ultimately allowing identification of a treatment target and region-specific biomarkers for AD.

Genome-wide gene and serum ferritin interaction in the development of type 2 diabetes in adults aged 40 years or older.

BACKGROUND AND AIMS: Elevated serum ferritin is associated with incident Type 2 diabetes (T2D), but the interactions between serum ferritin and genetic factors which may improve understanding underlying mechanism in the development of T2D are still unclear. We determined the gene-ferritin interactions on the development of T2D by genome-wide gene-ferritin interaction analyses. METHODS AND RESULTS: A total of 3405 participants from two prospective cohorts of community living residents were included, and the median follow-time was 3.99 years. Genome-wide gene-ferritin interactions were analyzed using the joint test with two degrees of freedom and the interaction test with one degree of freedom. There were 18 SNPs selected in the joint test. Finally, four independent variants [rs355140 (LINC00312), rs4075576 (nearby PDGFA), rs1332202 (PTPRD), and rs713157 (nearby LINC00900)] with low pairwise linkage disequilibrium (r2<0.2) and located at least 1000 kb from the index SNP showed interactions with serum ferritin level. In the association analyses between serum ferritin levels (tertiles of ferritin and ferritin status) and the incidence of T2D according to genotype, the Incidence Rate Ratios (IRRs) in the highest tertile of ferritin level (vs. the lowest tertile) were greater for participants with heterozygotes of risk alleles of each of the four SNP than IRRs for those with wild type. Compared with the normal group, the elevated ferritin group also had a higher risk of T2D for all genetic variants of risk alleles, particularly its homozygotes. CONCLUSION: Serum ferritin level interacts with genetic variants (rs355140, rs4075576, rs1332202, and rs713157) in the development of T2D.

Alternative Splicing Regulation of Low-Frequency Genetic Variants in Exon 2 of TREM2 in Alzheimer's Disease by Splicing-Based Aggregation.

TREM2 is among the most well-known Alzheimer's disease (AD) risk genes; however, the functional roles of its AD-associated variants remain to be elucidated, and most known risk alleles are low-frequency variants whose investigation is challenging. Here, we utilized a splicing-guided aggregation method in which multiple low-frequency TREM2 variants were bundled together to investigate the functional impact of those variants on alternative splicing in AD. We analyzed whole genome sequencing (WGS) and RNA-seq data generated from cognitively normal elderly controls (CN) and AD patients in two independent cohorts, representing three regions in the frontal lobe of the human brain: the dorsolateral prefrontal cortex (CN = 213 and AD = 376), frontal pole (CN = 72 and AD = 175), and inferior frontal (CN = 63 and AD = 157). We observed an exon skipping event in the second exon of TREM2, with that exon tending to be more frequently skipped (p = 0.0012) in individuals having at least one low-frequency variant that caused loss-of-function for a splicing regulatory element. In addition, genes differentially expressed between AD patients with high vs. low skipping of the second exon (i.e., loss of a TREM2 functional domain) were significantly enriched in immune-related pathways. Our splicing-guided aggregation method thus provides new insight into the regulation of alternative splicing of the second exon of TREM2 by low-frequency variants and could be a useful tool for further exploring the potential molecular mechanisms of multiple, disease-associated, low-frequency variants.

ADAS-viewer: web-based application for integrative analysis of multi-omics data in Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative disorder and is represented by complicated biological mechanisms and complexity of brain tissue. Our understanding of the complicated molecular architecture that contributes to AD progression benefits from performing comprehensive and systemic investigations with multi-layered molecular and biological data from different brain regions. Since recently different independent studies generated various omics data in different brain regions of AD patients, multi-omics data integration can be a useful resource for better comprehensive understanding of AD. Here we present a web platform, ADAS-viewer, that provides researchers with the ability to comprehensively investigate and visualize multi-omics data from multiple brain regions of AD patients. ADAS-viewer offers means to identify functional changes in transcript and exon expression (i.e., alternative splicing) along with associated genetic or epigenetic regulatory effects. Specifically, it integrates genomic, transcriptomic, methylation, and miRNA data collected from seven different brain regions (cerebellum, temporal cortex, dorsolateral prefrontal cortex, frontal pole, inferior frontal gyrus, parahippocampal gyrus, and superior temporal gyrus) across three independent cohort datasets. ADAS-viewer is particularly useful as a web-based application for analyzing and visualizing multi-omics data across multiple brain regions at both transcript and exon level, allowing the identification of candidate biomarkers of Alzheimer's disease.

Independent and interactive associations of season, dietary vitamin D, and vitamin D-related genetic variants with serum 25(OH)D in Korean adults aged 40 years or older.

Only limited information is available on the inter-relationships between genetic and non-genetic factors such as diet and sunlight exposure with serum 25-hydroxyvitamin D [25(OH)D] concentration. This cross-sectional study aimed to examine the independent and interactive associations of season, dietary vitamin D intake, and SNPs of 11 vitamin D-related candidate genes with serum 25(OH)D concentration among 2,721 adults aged ≥40 years at baseline from the Yangpyeong cohort, a part of the Korean Genome Epidemiology Study (KoGES). The interactions between season or dietary vitamin D and 556 SNPs were evaluated using 2-degree of freedom joint tests. Season was strongly (pdifference = 1.00 × 10-12) and dietary vitamin D intake was slightly but significantly associated with serum 25(OH)D concentration (pdifference = 0.0119). Among five SNPs (rs11723621-GC, rs7041-GC, rs10500804-CYP2R1, rs7129781-CYP2R1, and rs2852853-DHCR7) identified in the screening steps, only one, rs10500804-CYP2R1, significantly interacted with season (pinteraction = 8.01 × 10-5). The inverse association between number of minor alleles of rs10500804-CYP2R1 and concentration of 25(OH)D was significant only in summer/fall. Conversely, dietary vitamin D intake was positively associated only in winter/spring. In conclusion, season, dietary vitamin D intake, and four SNPs in GC, CYP2R1, and DHCR7 are independently and rs10500804-CYP2R1 is interactively associated with serum 25(OH)D concentration. Serum 25(OH)D is influenced by genotype of rs10500804-CYP2R1 in summer/fall when sunlight exposure is high, while dietary vitamin D intake is an important determinant of serum 25(OH)D during the seasons with low cutaneous vitamin D synthesis.

Association of an APBA3 Missense Variant with Risk of Premature Ovarian Failure in the Korean Female Population.

Premature ovarian failure (POF) is a complex disease of which the etiology is influenced by numerous genetic variations. Several POF candidate genes have been reported. However, no causal genes with high odds ratio (OR) have yet been discovered. This study included 564 females of Korean ethnicity, comprising 60 patients with POF and 182 controls in the discovery set and 105 patients with POF and 217 controls in the replication set. We conducted genome-wide association analysis to search for novel candidate genes predicted to influence POF development using Axiom Precision Medicine Research Arrays and additive model logistic regression analysis. One statistically significant single nucleotide polymorphism (SNP), rs55941146, which encodes a missense alteration (Val > Gly) in the APBA3 gene, was identified with OR values for association with POF of 13.33 and 4.628 in the discovery and replication sets, respectively. No rs55941146 minor allele homozygotes were present in either cases or controls. The APBA3 protein binds FIH-1 that inhibits hypoxia inducible factor-1α (HIF-1α). HIF-1α contributes to granulosa cell proliferation, which is crucial for ovarian follicle growth, by regulating cell proliferation factors and follicle stimulating hormone-mediated autophagy. Our data demonstrate that APBA3 is a candidate novel causal gene for POF.

Analysis of unmapped regions associated with long deletions in Korean whole genome sequences based on short read data.

While studies aimed at detecting and analyzing indels or single nucleotide polymorphisms within human genomic sequences have been actively conducted, studies on detecting long insertions/deletions are not easy to orchestrate. For the last 10 years, the availability of long read data of human genomes from PacBio or Nanopore platforms has increased, which makes it easier to detect long insertions/deletions. However, because long read data have a critical disadvantage due to their relatively high cost, many next generation sequencing data are produced mainly by short read sequencing machines. Here, we constructed programs to detect so-called unmapped regions (UMRs, where no reads are mapped on the reference genome), scanned 40 Korean genomes to select UMR long deletion candidates, and compared the candidates with the long deletion break points within the genomes available from the 1000 Genomes Project (1KGP). An average of about 36,000 UMRs were found in the 40 Korean genomes tested, 284 UMRs were common across the 40 genomes, and a total of 37,943 UMRs were found. Compared with the 74,045 break points provided by the 1KGP, 30,698 UMRs overlapped. As the number of compared samples increased from 1 to 40, the number of UMRs that overlapped with the break points also increased. This eventually reached a peak of 80.9% of the total UMRs found in this study. As the total number of overlapped UMRs could probably grow to encompass 74,045 break points with the inclusion of more Korean genomes, this approach could be practically useful for studies on long deletions utilizing short read data.

Differential Expression Profiling of Long Noncoding RNA and mRNA during Osteoblast Differentiation in Mouse.

Long noncoding RNAs (lncRNAs) are emerging as an important controller affecting metabolic tissue development, signaling, and function. However, little is known about the function and profile of lncRNAs in osteoblastic differentiation in mice. Here, we analyzed the RNA-sequencing (RNA-Seq) datasets obtained for 18 days in two-day intervals from neonatal mouse calvarial pre-osteoblast-like cells. Over the course of osteoblast differentiation, 4058 mRNAs and 3948 lncRNAs were differentially expressed, and they were grouped into 12 clusters according to the expression pattern by fuzzy c-means clustering. Using weighted gene coexpression network analysis, we identified 9 modules related to the early differentiation stage (days 2-8) and 7 modules related to the late differentiation stage (days 10-18). Gene ontology and KEGG pathway enrichment analysis revealed that the mRNA and lncRNA upregulated in the late differentiation stage are highly associated with osteogenesis. We also identified 72 mRNA and 89 lncRNAs as potential markers including several novel markers for osteoblast differentiation and activation. Our findings provide a valuable resource for mouse lncRNA study and improves our understanding of the biology of osteoblastic differentiation in mice.

KoVariome: Korean National Standard Reference Variome database of whole genomes with comprehensive SNV, indel, CNV, and SV analyses.

High-coverage whole-genome sequencing data of a single ethnicity can provide a useful catalogue of population-specific genetic variations, and provides a critical resource that can be used to more accurately identify pathogenic genetic variants. We report a comprehensive analysis of the Korean population, and present the Korean National Standard Reference Variome (KoVariome). As a part of the Korean Personal Genome Project (KPGP), we constructed the KoVariome database using 5.5 terabases of whole genome sequence data from 50 healthy Korean individuals in order to characterize the benign ethnicity-relevant genetic variation present in the Korean population. In total, KoVariome includes 12.7M single-nucleotide variants (SNVs), 1.7M short insertions and deletions (indels), 4K structural variations (SVs), and 3.6K copy number variations (CNVs). Among them, 2.4M (19%) SNVs and 0.4M (24%) indels were identified as novel. We also discovered selective enrichment of 3.8M SNVs and 0.5M indels in Korean individuals, which were used to filter out 1,271 coding-SNVs not originally removed from the 1,000 Genomes Project when prioritizing disease-causing variants. KoVariome health records were used to identify novel disease-causing variants in the Korean population, demonstrating the value of high-quality ethnic variation databases for the accurate interpretation of individual genomes and the precise characterization of genetic variations.

Photosynthetic functions of Synechococcus in the ocean microbiomes of diverse salinity and seasons.

Synechococcus is an important photosynthetic picoplankton in the temperate to tropical oceans. As a photosynthetic bacterium, Synechococcus has an efficient mechanism to adapt to the changes in salinity and light intensity. The analysis of the distributions and functions of such microorganisms in the ever changing river mouth environment, where freshwater and seawater mix, should help better understand their roles in the ecosystem. Toward this objective, we have collected and sequenced the ocean microbiome in the river mouth of Kwangyang Bay, Korea, as a function of salinity and temperature. In conjunction with comparative genomics approaches using the sequenced genomes of a wide phylogeny of Synechococcus, the ocean microbiome was analyzed in terms of their composition and clade-specific functions. The results showed significant differences in the compositions of Synechococcus sampled in different seasons. The photosynthetic functions in such enhanced Synechococcus strains were also observed in the microbiomes in summer, which is significantly different from those in other seasons.

Idiopathic hypereosinophilia is clonal disorder? Clonality identified by targeted sequencing.

Idiopathic hypereosinophilia (IHE)/idiopathic hypereosinophilic syndrome (IHES) has been defined by a persistent elevation of the blood eosinophil count exceeding 1.5×103/µL, without evidence of reactive or clonal causes. While T-cell clonality assessment has been recommended for unexplained hypereosinophilia, this approach is not often applied to routine practice in the clinic. We hypothesized that the clonality would exist in a subset of IHE/IHES patients. We aimed to investigate the candidate mutations and T-cell clonality in IHE/IHES and to explore the role of mutations in eosinophil proliferation. We performed targeted capture sequencing for 88 genes using next-generation sequencing, T-cell receptor (TCR) gene rearrangement assays, and pathway network analysis in relation to eosinophil proliferation. By targeted sequencing, 140 variants in 59 genes were identified. Sixteen out of 30 patients (53.3%) harbored at least one candidate mutation. The most frequently affected genes were NOTCH1 (26.7%), SCRIB and STAG2 (16.7%), and SH2B3 (13.3%). Network analysis revealed that our 21 candidate genes (BIRC3, BRD4, CSF3R, DNMT3A, EGR2, EZH2, FAT4, FLT3, GATA2, IKZF, JAK2, MAPK1, MPL, NF1, NOTCH1, PTEN, RB1, RUNX1, TET2, TP53 and WT1) are functionally linked to the eosinophilopoietic pathway. Among the 21 candidate genes, five genes (MAPK1, RUNX1, GATA2, NOTCH1 and TP53) with the highest number of linkages were considered major genes. A TCR assay revealed that four patients (13.3%) had a clonal TCR rearrangement. NOTCH1 was the most frequently mutated gene and was shown to be a common node for eosinophilopoiesis in our network analysis, while the possibility of hidden T cell malignancy was indwelling in the presence of NOTCH1 mutation, though not revealed by aberrant T cell study. Collectively, these results provide new evidence that mutations affecting eosinophilopoiesis underlie a subset of IHE/IHES, and the candidate genes are inferred to act their potential roles in the eosinophilopoietic pathway.

Pan-genome analysis of Bacillus for microbiome profiling.

Recent advances in high-throughput sequencing technology allow for in-depth studies on microbial genomes and their communities. While multiple strains of the same species could display genomic variations with different gene contents in diverse habitats and hosts, the essential functions for a specific species are conserved as core genes that are shared among strains. We have comprehensively analyzed 238 strains of five different Bacillus species to identify the properties of core and strain-specific genes. Core and strain-specific genes in each Bacillus species show significant differences in their functions and genomic signatures. Using the core genes defined in this study, we have precisely identified the Bacillus species that exist in food microbiomes. Without resorting to culture-based whole genome sequencing, an unexpectedly large portion of the core genes, 98.22% of core genes in B. amyloliquefaciens and 97.77% of B. subtilis, were reconstructed from the microbiome. We have performed a pan-genome analysis on the core gene data of multiple Bacillus species to investigate the Bacillus species in food microbiome. Our findings provide a comprehensive genetic landscape of the Bacillus species, which is also consistent with previous studies on a limited number of strains and species. Analysis based on comprehensive core genes should thus serve as a powerful profiling tool to better understand major constituents in fermented food microbiomes.

RNA-Seq for Gene Expression Profiling of Human Necrotizing Enterocolitis: a Pilot Study.

Necrotizing enterocolitis (NEC) characterized by inflammatory intestinal necrosis is a major cause of mortality and morbidity in newborns. Deep RNA sequencing (RNA-Seq) has recently emerged as a powerful technology enabling better quantification of gene expression than microarrays with a lower background signal. A total of 10 transcriptomes from 5 pairs of NEC lesions and adjacent normal tissues obtained from preterm infants with NEC were analyzed. As a result, a total of 65 genes (57 down-regulated and 8 up-regulated) revealed significantly different expression levels in the NEC lesion compared to the adjacent normal region, based on a significance at fold change ≥ 1.5 and P ≤ 0.05. The most significant gene, DPF3 (P < 0.001), has recently been reported to have differential expressions in colon segments. Our gene ontology analysis between NEC lesion and adjacent normal tissues showed that down-regulated genes were included in nervous system development with the most significance (P = 9.3 × 10⁻7; P(corr) = 0.0003). In further pathway analysis using Pathway Express based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, genes involved in thyroid cancer and axon guidance were predicted to be associated with different expression (P(corr) = 0.008 and 0.020, respectively). Although further replications using a larger sample size and functional evaluations are needed, our results suggest that altered gene expression and the genes' involved functional pathways and categories may provide insight into NEC development and aid in future research.

Mitochondrial DNA copy number augments performance of A1C and oral glucose tolerance testing in the prediction of type 2 diabetes.

Here, we tested the performance of the mitochondrial DNA copy number (mtDNA-CN) in predicting future type 2 diabetes (n = 1108). We used the baseline clinical data (age, sex, body mass index, waist-to-hip ratio, systolic and diastolic blood pressure) and the mtDNA-CN, hemoglobin A1c (A1C) levels and results of oral glucose tolerance test (OGTT) including fasting plasma glucose, 1-hour glucose, and 2-hour glucose levels, to predict future diabetes. We built a prediction model using the baseline data and the diabetes status at biannual follow-up of 8 years. The mean area under curve (AUC) for all follow-ups of the full model including all variables was 0.92 ± 0.04 (mean ± standard deviation), while that of the model excluding the mtDNA-CN was 0.90 ± 0.03. The sensitivity of the f4ull model was much greater than that of the model not including mtDNA-CN: the mean sensitivities of the model with and without mtDNA-CN were 0.60 ± 0.06 and 0.53 ± 0.04, respectively. We found that the mtDNA-CN of peripheral leukocytes is a biomarker that augments the predictive power for future diabetes of A1C and OGTT. We believe that these results could provide invaluable information for developing strategies for the management of diabetes.

[Corrigendum] Promoter methylation of RASSF1A modulates the effect of the microtubule-targeting agent docetaxel in breast cancer.

Following the publication of this article, an interested reader drew to our attention that there were possible anomalies in the presentation of Fig. 5B in the above article. After having examined the figure, we recognized that several errors had indeed occurred during the process of compiling the figure. A corrected version of Fig. 5 is shown below, containing new data for Fig. 5B, after our having re-performed the western blot experiment according to the identical procedure detailed in the paper. We obtained broadly similar results to those featured originally in the article; therefore, the revision of this figure does not affect the conclusions reported in the study. We thank the reader of our article who drew this matter to our attention. [the original article was published in the International Journal of Oncology 41: 611-620, 2012; DOI: 10.3892/ijo.2012.1470].

A genome-wide association analysis of chromosomal aberrations and Hirschsprung disease.

Hirschsprung disease (HSCR) is a neurocristopathy characterized by the absence of intramural ganglion cells along variable lengths of the gastrointestinal tract. Although the RET proto-oncogene is considered to be the main risk factor for HSCR, only about 30% of the HSCR cases can be explained by variations in previously known genes including RET. Recently, copy number variation (CNV) and loss of heterozygosity (LOH) have emerged as new ways to understand human genomic variation. The goal of this present study is to identify new HSCR genetic factors related to CNV in Korean patients. In the genome-wide genotyping, using Illumina's HumanOmni1-Quad BeadChip (1,140,419 markers), of 123 HSCR patients and 432 unaffected subjects (total n = 555), a total of 8,188 CNVs (1 kb âˆ¼ 1 mb) were identified by CNVpartition. As a result, 16 CNV regions and 13 LOH regions were identified as associated with HSCR (minimum P = 0.0005). Two top CNV regions (deletions at chr6:32675155-32680480 and chr22:20733495-21607293) were successfully validated by additional real-time quantitative polymerase chain reaction analysis. In addition, 2 CNV regions (6p21.32 and 22q11.21) and 2 LOH regions (3p22.2 and 14q23.3) were discovered to be unique to the HSCR patients group. Regarding the large-scale chromosomal aberrations (>1 mb), 11 large aberrations in the HSCR patients group were identified, which suggests that they may be a risk factor for HSCR. Although further replication in a larger cohort is needed, our findings may contribute to the understanding of the etiology of HSCR.

Association Analysis of SLC6A20 Polymorphisms With Hirschsprung Disease.

PURPOSE: Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder, which is caused by no neuronal ganglion cells in part or all of distal gastrointestinal tract. Recently, our genome-wide association study has identified solute carrier family 6, proline IMINO transporter, member 20 (SLC6A20) as one of the potential risk factors for HSCR development. This study performed a replication study for the association of SLC6A20 polymorphisms with HSCR and an extended analysis to investigate further associations for subgroups and haplotypes. METHODS: For the replication study, a total of 40 single nucleotide polymorphisms (SNPs) of SLC6A20 were genotyped in 187 HSCR subjects composed of 121 short-segment HSCR, 45 long-segment HSCR (L-HSCR), 21 total colonic aganglionosis, and 283 unaffected controls. Imputation was performed using genotype data from our genome-wide association study and this replication study. RESULTS: Imputed meta-analysis revealed that 13 SLC6A20 SNPs (minimum P = 0.0002 at rs6770261) were significantly associated with HSCR even after correction for multiple comparisons using false discovery rate (FDR) (minimum PFDR =  .005). In further subgroup analysis, SLC6A20 polymorphisms appeared to have increased associations with L-HSCR. Moreover, haplotype analysis also showed significant associations between 2 haplotypes (BL3_ht2 and BL4_ht2) and HSCR susceptibility (PFDR <  .05). CONCLUSIONS: Although further replications and functional evaluations are required, our results suggest that SLC6A20 may have roles in HSCR development and in the extent of aganglionic segment during enteric nervous system development.

Replication of genome wide association studies on hepatocellular carcinoma susceptibility loci of STAT4 and HLA-DQ in a Korean population.

A recent genome-wide association study (GWAS) for hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) identified two loci (rs7574865 in STAT4 and rs9275319 in HLA-DQ) in a Chinese population. We attempted to replicate the associations between the two SNP loci and the risk of HCC in a Korean population. The rs7574865 in STAT4 and rs9275319 in HLA-DQ were genotyped in a total of 3838 Korean subjects composed of 287 HBV-related hepatocellular carcinoma patients, 671 chronic hepatitis B virus (CHB) patients, and 2880 population controls using TaqMan genotyping assay. Gene expression was measured by microarray. A logistic regression analysis revealed that rs7574865 in STAT4 and rs9275319 in HLA-DQ were associated with the risk of CHB (OR = 1.25, P = 0.0002 and OR = 1.57, P= 1.44 × 10(-10), respectively). However, these loci were no association with the risk of HBV-related HCC among CHB patients. In the gene expression analyses, although no significant differences in mRNA expression of nearby genes according to genotypes were detected, a significantly decreased mRNA expression in HCC subjects was observed in STAT4, HLA-DQA1, and HLA-DQB1. Although the genetic effects of two HCC susceptibility loci were not replicated, the two loci were found to exert susceptibility effects on the risk of CHB in a Korean population. In addition, the decreased mRNA expression of STAT4, HLA-DQA1, and HLA-DQB1 in HCC tissue might provide a clue to understanding their role in the progression to HCC.