Genome-Wide Association Study of Susceptibility Loci for Radiation-Induced Brain Injury.

BACKGROUND: Radiation-induced brain injury is a nonnegligible issue in the management of cancer patients treated by partial or whole brain irradiation. In particular, temporal lobe injury (TLI), a deleterious late complication in nasopharyngeal carcinoma, greatly affects the long-term life quality of these patients. Although genome-wide association studies (GWASs) have successfully identified single nucleotide polymorphisms (SNPs) associated with radiation toxicity, genetic variants contributing to the radiation-induced brain injury have not yet been assessed. METHODS: We recruited and performed follow-up for a prospective observational cohort, Genetic Architecture of Radiotherapy Toxicity and Prognosis, using magnetic resonance imaging for TLI diagnosis. We conducted genome-wide association analysis in 1082 patients and validated the top associations in two independent cohorts of 1119 and 741 patients, respectively. All statistical tests were two-sided. RESULTS: We identified a promoter variant rs17111237 (A > G, minor allele frequency [MAF] = 0.14) in CEP128 associated with TLI risk (hazard ratio = 1.45, 95% confidence interval = 1.26 to 1.66, Pcombined=3.18 × 10-7) which is in moderate linkage disequilibrium (LD) with rs162171 (MAF = 0.18, R2 = 0.69), the top signal in CEP128 (hazard ratio = 1.46, 95% confidence interval = 1.29-1.66, Pcombined= 6.17 × 10-9). Combining the clinical variables with the top SNP, we divided the patients into different subgroups with varying risk with 5-year TLI-free rates ranging from 33.7% to 95.5%. CEP128, a key component of mother centriole, tightly interacts with multiple radiation-resistant genes and plays an important role in maintaining the functional cilia, which otherwise will lead to a malfunction of the neural network. We found that A > G alteration at rs17111237 impaired the promoter activity of CEP128 and knockdown of CEP128 decreased the clonogenic cell survival of U87 cells under radiation. Noteworthy, 12.7% (27/212) of the GWAS-based associated genes (P < .001) were enriched in the neurogenesis pathway. CONCLUSIONS: This three-stage study is the first GWAS of radiation-induced brain injury that implicates the genetic susceptibility gene CEP128 involved in TLI development and provides the novel insight into the underlying mechanisms of radiation-induced brain injury.

Bagging Nearest-Neighbor Prediction independence Test: an efficient method for nonlinear dependence of two continuous variables.

Testing dependence/correlation of two variables is one of the fundamental tasks in statistics. In this work, we proposed an efficient method for nonlinear dependence of two continuous variables (X and Y). We addressed this research question by using BNNPT (Bagging Nearest-Neighbor Prediction independence Test, software available at https://sourceforge.net/projects/bnnpt/). In the BNNPT framework, we first used the value of X to construct a bagging neighborhood structure. We then obtained the out of bag estimator of Y based on the bagging neighborhood structure. The square error was calculated to measure how well Y is predicted by X. Finally, a permutation test was applied to determine the significance of the observed square error. To evaluate the strength of BNNPT compared to seven other methods, we performed extensive simulations to explore the relationship between various methods and compared the false positive rates and statistical power using both simulated and real datasets (Rugao longevity cohort mitochondrial DNA haplogroups and kidney cancer RNA-seq datasets). We concluded that BNNPT is an efficient computational approach to test nonlinear correlation in real world applications.

An Integrative Computational Approach to Evaluate Genetic Markers for Bipolar Disorder.

Studies to date have reported hundreds of genes connected to bipolar disorder (BP). However, many studies identifying candidate genes have lacked replication, and their results have, at times, been inconsistent with one another. This paper, therefore, offers a computational workflow that can curate and evaluate BP-related genetic data. Our method integrated large-scale literature data and gene expression data that were acquired from both postmortem human brain regions (BP case/control: 45/50) and peripheral blood mononuclear cells (BP case/control: 193/593). To assess the pathogenic profiles of candidate genes, we conducted Pathway Enrichment, Sub-Network Enrichment, and Gene-Gene Interaction analyses, with 4 metrics proposed and validated for each gene. Our approach developed a scalable BP genetic database (BP_GD), including BP related genes, drugs, pathways, diseases and supporting references. The 4 metrics successfully identified frequently-studied BP genes (e.g. GRIN2A, DRD1, DRD2, HTR2A, CACNA1C, TH, BDNF, SLC6A3, P2RX7, DRD3, and DRD4) and also highlighted several recently reported BP genes (e.g. GRIK5, GRM1 and CACNA1A). The computational biology approach and the BP database developed in this study could contribute to a better understanding of the current stage of BP genetic research and assist further studies in the field.

Altered expression of mRNA profiles in blood of early-onset schizophrenia.

To identify gene expression abnormalities in schizophrenia (SZ), we generated whole-genome gene expression profiles using microarrays on peripheral blood mononuclear cells (PBMCs) from 18 early-onset SZ cases and 12 controls. We detected 84 transcripts differentially expressed by diagnostic status, with 82 genes being upregulated and 2 downregulated. We identified two SZ associated gene coexpression modules (green and red), including 446 genes . The green module is positively correlated with SZ, encompassing predominantly up-regulated genes in SZ; while the red module was negatively correlated with disease status, involving mostly nominally down-regulated genes in SZ. The olfactory transduction pathway was the most enriched pathways for the genes within the two modules. The expression levels of several hub genes, including AKT1, BRCA1, CCDC134, UBD, and ZIC2 were validated using real-time quantitative PCR. Our findings indicate that mRNA coexpression abnormalities may serve as a promising mechanism underlying the development of SZ.

Exploring Transcription Factors-microRNAs Co-regulation Networks in Schizophrenia.

BACKGROUND: Transcriptional factors (TFs) and microRNAs (miRNAs) have been recognized as 2 classes of principal gene regulators that may be responsible for genome coexpression changes observed in schizophrenia (SZ). METHODS: This study aims to (1) identify differentially coexpressed genes (DCGs) in 3 mRNA expression microarray datasets; (2) explore potential interactions among the DCGs, and differentially expressed miRNAs identified in our dataset composed of early-onset SZ patients and healthy controls; (3) validate expression levels of some key transcripts; and (4) explore the druggability of DCGs using the curated database. RESULTS: We detected a differential coexpression network associated with SZ and found that 9 out of the 12 regulators were replicated in either of the 2 other datasets. Leveraging the differentially expressed miRNAs identified in our previous dataset, we constructed a miRNA-TF-gene network relevant to SZ, including an EGR1-miR-124-3p-SKIL feed-forward loop. Our real-time quantitative PCR analysis indicated the overexpression of miR-124-3p, the under expression of SKIL and EGR1 in the blood of SZ patients compared with controls, and the direction of change of miR-124-3p and SKIL mRNA levels in SZ cases were reversed after a 12-week treatment cycle. Our druggability analysis revealed that many of these genes have the potential to be drug targets. CONCLUSIONS: Together, our results suggest that coexpression network abnormalities driven by combinatorial and interactive action from TFs and miRNAs may contribute to the development of SZ and be relevant to the clinical treatment of the disease.

Increased Variability of Genomic Transcription in Schizophrenia.

Schizophrenia (SZ) is a severe chronic mental disorder with a high heritability. Current microarray analyses typically focus on identifying differentially expressed genes or enriched pathways relevant to phenotypes. Whether there is a variability change of the genomic transcription in diseases has rarely been explored. In this study, we compared coefficient of variation (CV, the ratio of the standard deviation to the mean) of genome transcription of early-onset SZ (EOS) patients with controls in a blood mRNA microarray dataset and a blood microRNA (miRNA) microarray dataset. Furthermore, we compared CV of the expression levels of 17 genes in blood of the 30 patients before and after the 12-week treatment using real-time quantitative PCR (RT-qPCR) analysis. Our results indicated a significant increase of CV of genome transcription in patients compared with controls in both the mRNA and the miRNA datasets. The 30 after-treatment patients showed a significant decrease of CV of gene expression compared with the before-treatment patients. Our study may implicate the blood gene expression variability in SZ, providing further evidence supporting the abnormality of peripheral blood transcriptome in SZ. Given that peripheral blood can be easily collected from patients and followed longitudinally, our results may indicate a new way to facilitate the identification of the signatures of clinical subtypes, their prognosis and treatment response.

Heterogeneity revealed through meta-analysis might link geographical differences with nasopharyngeal carcinoma incidence in Han Chinese populations.

BACKGROUND: Nasopharyngeal carcinoma (NPC) is an epithelial malignancy highly prevalent in southern China, and incidence rates among Han Chinese people vary according to geographic region. Recently, three independent genome-wide association studies (GWASs) confirmed that HLA-A is the main risk gene for NPC. However, the results of studies conducted in regions with dissimilar incidence rates contradicted the claims that HLA-A is the sole risk gene and that the association of rs29232 is independent of the HLA-A effect in the chromosome 6p21.3 region. METHODS: We performed a meta-analysis, selecting five single-nucleotide polymorphisms (SNPs) in chromosome 6p21.3 mapped in three published GWASs and four case-control studies. The studies involved 8994 patients with NPC and 11,157 healthy controls, all of whom were Han Chinese. RESULTS: The rs2517713 SNP located downstream of HLA-A was significantly associated with NPC (P = 1.08 × 10(-91), odds ratio [OR] = 0.58, 95 % confidence interval [CI] = 0.55-0.61). The rs29232 SNP exhibited a moderate level of heterogeneity (I(2) = 47 %) that disappeared (I(2) = 0 %) after stratification by moderate- and high-incidence NPC regions. CONCLUSIONS: Our results suggested that the HLA-A gene is strongly associated with NPC risk. In addition, the heterogeneity revealed by the meta-analysis of rs29232 might be associated with regional differences in NPC incidence among Han Chinese people. The higher OR of rs29232 and the fact that rs29232 was independent of the HLA-A effect in the moderate-incidence population suggested that rs29232 might have greater relevance to NPC incidence in a moderate-incidence population than in a high-incidence population.

How genome-wide SNP-SNP interactions relate to nasopharyngeal carcinoma susceptibility.

This study is the first to use genome-wide association study (GWAS) data to evaluate the multidimensional genetic architecture underlying nasopharyngeal cancer. Since analysis of data from GWAS confirms a close and consistent association between elevated risk for nasopharyngeal carcinoma (NPC) and major histocompatibility complex class 1 genes, our goal here was to explore lesser effects of gene-gene interactions. We conducted an exhaustive genome-wide analysis of GWAS data of NPC, revealing two-locus interactions occurring between single nucleotide polymorphisms (SNPs), and identified a number of suggestive interaction loci which were missed by traditional GWAS analyses. Although none of the interaction pairs we identified passed the genome-wide Bonferroni-adjusted threshold for significance, using independent GWAS data from the same population (Stage 2), we selected 66 SNP pairs in 39 clusters with P<0.01. We identified that in several chromosome regions, multiple suggestive interactions group to form a block-like signal, effectively reducing the rate of false discovery. The strongest cluster of interactions involved the CREB5 gene and a SNP rs1607979 on chromosome 17q22 (P = 9.86×10(-11)) which also show trans-expression quantitative loci (eQTL) association in Chinese population. We then detected a complicated cis-interaction pattern around the NPC-associated HLA-B locus, which is immediately adjacent to copy-number variations implicated in male susceptibility for NPC. While it remains to be seen exactly how and to what degree SNP-SNP interactions such as these affect susceptibility for nasopharyngeal cancer, future research on these questions holds great promise for increasing our understanding of this disease's genetic etiology, and possibly also that of other gene-related cancers.

A population-based study of hedgehog pathway gene variants in relation to the dual risk of basal cell carcinoma plus another cancer.

INTRODUCTION: A personal history of basal cell carcinoma (BCC) is associated with increased risk of other malignancies, but the reason is unknown. The hedgehog pathway is critical to the etiology of BCC, and is also believed to contribute to susceptibility to other cancers. This study tested the hypothesis that hedgehog pathway and pathway-related gene variants contribute to the increased risk of subsequent cancers among those with a history of BCC. METHODS: The study was nested within the ongoing CLUE II cohort study, established in 1989 in Washington County, Maryland, USA. The study consisted of a cancer-free control group (n=2296) compared to three different groups of cancer cases ascertained through 2007, those diagnosed with: (1) Other (non-BCC) cancer only (n=2349); (2) BCC only (n=534); and (3) BCC plus other cancer (n=446). The frequencies of variant alleles were compared among these four groups for 20 single nucleotide polymorphisms (SNPs) in 6 hedgehog pathway genes (SHH, IHH, PTCH2, SMO, GLI1, SUFU), and also 22 SNPs in VDR and 8 SNPs in FAS, which have cross-talk with the hedgehog pathway. RESULTS: Comparing those with both BCC and other cancer versus those with no cancer, no significant associations were observed for any of the hedgehog pathway SNPs, or for the FAS SNPs. One VDR SNP was nominally significantly associated with the BCC cancer-prone phenotype, rs11574085 [per minor allele odds ratio (OR) 1.38, 95% confidence interval (CI) 1.05-1.82; p-value=0.02]. CONCLUSION: The hedgehog pathway gene SNPs studied, along with the VDR and FAS SNPs studied, are not strongly associated with the BCC cancer-prone phenotype.

A gender-specific association of CNV at 6p21.3 with NPC susceptibility.

Copy number variations (CNVs), a major source of human genetic polymorphism, have been suggested to have an important role in genetic susceptibility to common diseases such as cancer, immune diseases and neurological disorders. Nasopharyngeal carcinoma (NPC) is a multifactorial tumor closely associated with genetic background and with a male preponderance over female (3:1). Previous genome-wide association studies have identified single-nucleotide polymorphisms (SNPs) that are associated with NPC susceptibility. Here, we sought to explore the possible association of CNVs with NPC predisposition. Utilizing genome-wide SNP-based arrays and five CNV-prediction algorithms, we identified eight regions with CNV that were significantly overrepresented in NPC patients compared with healthy controls. These CNVs included six deletions (on chromosomes 3, 6, 7, 8 and 19), and two duplications (on chromosomes 7 and 12). Among them, the CNV located at chromosome 6p21.3, with single-copy deletion of the MICA and HCP5 genes, showed the highest association with NPC. Interestingly, it was more specifically associated with an increased NPC risk among males. This gender-specific association was replicated in an independent case-control sample using a self-established deletion-specific polymerase chain reaction strategy. To the best of our knowledge, this is the first study to explore the role of constitutional CNVs in NPC, using a genome-wide platform. Moreover, we identified eight novel candidate regions with CNV that merit future investigation, and our results suggest that similar to neuroblastoma and prostate cancer, genetic structural variations might contribute to NPC predisposition.

Traditional Cantonese diet and nasopharyngeal carcinoma risk: a large-scale case-control study in Guangdong, China.

BACKGROUND: Nasopharyngeal carcinoma (NPC) is rare in most parts of the world but is a common malignancy in southern China, especially in Guangdong. Dietary habit is regarded as an important modifier of NPC risk in several endemic areas and may partially explain the geographic distribution of NPC incidence. In China, rapid economic development during the past few decades has changed the predominant lifestyle and dietary habits of the Chinese considerably, requiring a reassessment of diet and its potential influence on NPC risk in this NPC-endemic area. METHODS: To evaluate the association between dietary factors and NPC risk in Guangdong, China, a large-scale, hospital-based case-control study was conducted. 1387 eligible cases and 1459 frequency matched controls were recruited. Odds ratios (ORs) and the corresponding 95% confidence intervals (CIs) were estimated using a logistic regression model, adjusting for age, sex, education, dialect, and habitation household type. RESULTS: Observations made include the following: 1) consumption of canton-style salted fish, preserved vegetables and preserved/cured meat were significantly associated with increased risk of NPC, with enhanced odds ratios (OR) of 2.45 (95% CI: 2.03-2.94), 3.17(95% CI: 2.68-3.77) and 2.09 (95% CI: 1.22-3.60) respectively in the highest intake frequency stratum during childhood; 2) consumption of fresh fruit was associated with reduced risk with a dose-dependent relationship (p = 0.001); and 3) consumption of Canton-style herbal tea and herbal slow-cooked soup was associated with decreased risk, with ORs of 0.84 (95% CI: 0.68-1.03) and 0.58 (95% CI: 0.47-0.72) respectively in the highest intake frequency stratum. In multivariate analyses, these associations remained significant. CONCLUSIONS: It can be inferred that previously established dietary risk factors in the Cantonese population are still stable and have contributed to the incidence of NPC.

Replication study supports evidence for linkage to 9p24 in obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a severe psychiatric illness that is characterized by intrusive and senseless thoughts and impulses (obsessions) and by repetitive behaviors (compulsions). Family, twin, and segregation studies support the presence of both genetic and environmental susceptibility factors, and the only published genome scan for OCD identified a candidate region on 9p24 at marker D9S288 that met criteria for suggestive significance (Hanna et al. 2002). In an attempt to replicate this finding, we genotyped 50 pedigrees with OCD, using microsatellite markers spanning the 9p24 candidate region, and analyzed the data, using parametric and nonparametric linkage analyses under both a narrow phenotype model (DSM-IV OCD definite; 41 affected sib pairs) and a broad phenotype model (DSM-IV OCD definite and probable; 50 affected sib pairs). Similar to what was described by Hanna et al. (2002), our strongest findings came with the dominant parameters and the narrow phenotype model: the parametric signal peaked at marker D9S1792 with an HLOD of 2.26 ( alpha =0.59), and the nonparametric linkage signal (NPL) peaked at marker D9S1813 with an NPL of 2.52 (P=.006). These findings are striking in that D9S1813 and D9S1792 lie within 0.5 cM (<350 kb) of the original 9p24 linkage signal at D9S288; furthermore, pedigree-based association analyses also implicated the 9p24 candidate region by identifying two markers (D9S288 and GATA62F03) with modest evidence (P=.046 and .02, respectively) for association.