Zhijian Song

Recurrent gain-of-function USP8 mutations in Cushing's disease

Cushings disease, also known as adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (PAs) that cause excess cortisol production, accounts for up to 85% of corticotrophin-dependent Cushings syndrome cases. However, the genetic alterations in this disease are unclear. Here, we performed whole-exome sequencing of DNA derived from 12 ACTH-secreting PAs and matched blood samples, which revealed three types of somatic mutations in a candidate gene, USP8 (encoding ubiquitin-specific protease 8), exclusively in exon 14 in 8 of 12 ACTH-secreting PAs. We further evaluated somatic USP8 mutations in additional 258 PAs by Sanger sequencing. Targeted sequencing further identified a total of 17 types of USP8 variants in 67 of 108 ACTH-secreting PAs (62.04%). However, none of these mutations was detected in other types of PAs (n = 150). These mutations aggregate within the 14-3-3 binding motif of USP8 and disrupt the interaction between USP8 and 14-3-3 protein, resulting in an elevated capacity to protect EGFR from lysosomal degradation. Accordingly, PAs with mutated USP8 display a higher incidence of EGFR expression, elevated EGFR protein abundance and mRNA expression levels of POMC, which encodes the precursor of ACTH. PAs with mutated USP8 are significantly smaller in size and have higher ACTH production than wild-type PAs. In surgically resected primary USP8-mutated tumor cells, USP8 knockdown or blocking EGFR effectively attenuates ACTH secretion. Taken together, somatic gain-of-function USP8 mutations are common and contribute to ACTH overproduction in Cushings disease. Inhibition of USP8 or EGFR is promising for treating USP8-mutated corticotrophin adenoma. Our study highlights the potentially functional mutated gene in Cushings disease and provides insights into the therapeutics of this disease.

Genome-wide association analysis identifies 30 new susceptibility loci for schizophrenia

We conducted a genome-wide association study (GWAS) with replication in 36,180 Chinese individuals and performed further transancestry meta-analyses with data from the Psychiatry Genomics Consortium (PGC2). Approximately 95% of the genome-wide significant (GWS) index alleles (or their proxies) from the PGC2 study were overrepresented in Chinese schizophrenia cases, including ∼50% that achieved nominal significance and ∼75% that continued to be GWS in the transancestry analysis. The Chinese-only analysis identified seven GWS loci; three of these also were GWS in the transancestry analyses, which identified 109 GWS loci, thus yielding a total of 113 GWS loci (30 novel) in at least one of these analyses. We observed improvements in the fine-mapping resolution at many susceptibility loci. Our results provide several lines of evidence supporting candidate genes at many loci and highlight some pathways for further research. Together, our findings provide novel insight into the genetic architecture and biological etiology of schizophrenia.

The GSK3B gene confers risk for both major depressive disorder and schizophrenia in the Han Chinese population

BACKGROUND Glycogen synthease kinase-3B is a key gene encoding a protein kinase which is abundant in brain, and is involved in signal transduction cascades of neuronal cell development and energy metabolism. Previous researches proposed GSK3B as a potential region for schizophrenia. METHOD To validate the susceptibility of GSK3B to major depressive disorder, and to investigate the overlapping risk conferred by GSK3B for mental disorders, we performed a large-scale case-control study, analyzed 6 tag single nucleotide polymorphisms using TaqMan® technology in 1,045 major depressive disorder patients, 1,235 schizophrenia patients and 1,235 normal controls of Han Chinese origin. RESULTS We found rs334535 (Pallele=2.79E-03, Pgenotype=5.00E-03, OR=1.429) and rs2199503 (Pallele=0.020, Pgenotype= 0.040, OR=1.157) showed association with major depressive disorder before Bonferroni correction. rs6771023 (adjusted Pallele=1.64E-03, adjusted Pgenotype=6.00E-03, OR=0.701) and rs2199503 (adjusted Pallele=0.001, adjusted Pgenotype=0.002, OR=1.251) showed significant association with schizophrenia after Bonferroni correction. rs2199503 (adjusted Pallele=1.70E-03, adjusted Pgenotype=0.006, OR=1.208) remained to be significant in the combined cases of major depressive disorder and schizophrenia after Bonferroni correction. LIMITATIONS Further validations of our findings in samples with larger scale are suggested, and functional genomic study is needed to elucidate the role of GSK3B in signal pathway and psychiatric disorders. CONCLUSIONS Our results provide evidence that the GSK3B gene could be a promising region which contains genetic risk for both major depressive disorder and schizophrenia in the Han Chinese population. The study on variants conferring overlapping risk for multiple psychiatric disorders could be tangible pathogenesis support and clinical or diagnostic references.

SHEsisPlus, a toolset for genetic studies on polyploid species

Currently, algorithms and softwares for genetic analysis of diploid organisms with bi-allelic markers are well-established, while those for polyploids are limited. Here, we present SHEsisPlus, the online algorithm toolset for both dichotomous and quantitative trait genetic analysis on polyploid species (compatible with haploids and diploids, too). SHEsisPlus is also optimized for handling multiple-allele datasets. It’s free, open source and also designed to perform a range of analyses, including haplotype inference, linkage disequilibrium analysis, epistasis detection, Hardy-Weinberg equilibrium and single locus association tests. Meanwhile, we developed an accurate and efficient haplotype inference algorithm for polyploids and proposed an entropy-based algorithm to detect epistasis in the context of quantitative traits. A study of both simulated and real datasets showed that our haplotype inference algorithm was much faster and more accurate than existing ones. Our epistasis detection algorithm was the first try to apply information theory to characterizing the gene interactions in quantitative trait datasets. Results showed that its statistical power was significantly higher than conventional approaches. SHEsisPlus is freely available on the web at http://shesisplus.bio-x.cn/. Source code is freely available for download at https://github.com/celaoforever/SHEsisPlus.

Genome-wide Analysis of the Role of Copy Number Variation in Schizophrenia Risk in Chinese

BACKGROUND Compelling evidence suggested the role of copy number variations (CNVs) in schizophrenia susceptibility. Most of the evidence was from studies in populations with European ancestry. We tried to validate the associated CNV loci in a Han Chinese population and identify novel loci conferring risk of schizophrenia. METHODS We performed a genome-wide CNV analysis on 6588 patients with schizophrenia and 11,904 control subjects of Han Chinese ancestry. RESULTS Our data confirmed increased genome-wide CNV (>500 kb and <1%) burden in schizophrenia, and the increasing trend was more significant when only >1 Mb CNVs were considered. We also replicated several associated loci that were previously identified in European populations, including duplications at 16p11.2, 15q11.2-13.1, 7q11.23, and VIPR2 and deletions at 22q11.2, 1q21.1-q21.2, and NRXN1. In addition, we discovered three additional new potential loci (odds ratio >6, p < .05): duplications at 1p36.32, 10p12.1, and 13q13.3, involving many neurodevelopmental and synaptic related genes. CONCLUSIONS Our findings provide further support for the role of CNVs in the etiology of schizophrenia.

Genome-Wide Association Study of Bladder Cancer in a Chinese Cohort Reveals a New Susceptibility Locus at 5q12.3

Genome-wide association studies (GWAS) of bladder cancer have identified a number of susceptibility loci in European populations but have yet to uncover the genetic determinants underlying bladder cancer incidence among other ethnicities. Therefore, we performed the first GWAS in a Chinese cohort comprising 3,406 cases of bladder cancer and 4,645 controls. We identified a new susceptibility locus for bladder cancer at 5q12.3, located in the intron of CWC27 (rs2042329), that was significantly associated with disease risk (OR = 1.40; P = 4.61 × 10(-11)). However, rs2042329 was not associated with bladder cancer risk in patients of European descent. The rs2042329 risk allele was also related to significantly increased expression levels of CWC27 mRNA and protein in bladder cancer tissues from Chinese patients. Additional functional analyses suggested that CWC27 played an oncogenic role in bladder cancer by inducing cell proliferation and suppressing apoptosis. In conclusion, the identification of a risk-associated locus at 5q12.3 provides new insights into the inherited susceptibility to bladder cancer in Chinese populations and may help to identify high-risk individuals. Cancer Res; 76(11); 3277-84. ©2016 AACR.

Common variants at 10p12.31, 10q21.1 and 13q12.13 are associated with sporadic pituitary adenoma

Pituitary adenoma is one of the most common intracranial neoplasms, and its genetic basis remains largely unknown. To identify genetic susceptibility loci for sporadic pituitary adenoma, we performed a three-stage genome-wide association study (GWAS) in the Han Chinese population. We first analyzed genome-wide SNP data in 771 pituitary adenoma cases and 2,788 controls and then carried forward the promising variants for replication in another 2 independent sets (2,542 cases and 3,620 controls in total). We identified three new susceptibility loci below the genome-wide significance threshold (P < 5 × 10−8) in the combined analyses: 10p12.31 (rs2359536, Pmeta = 2.25 × 10−10 and rs10828088, Pmeta = 6.27 × 10−10), 10q21.1 (rs10763170, Pmeta = 6.88 × 10−10) and 13q12.13 (rs17083838, Pmeta = 1.89 × 10−8). This study is the first GWAS to our knowledge on sporadic pituitary adenoma, and our results provide insight into the genetic basis of this disease.

Endothelial Nitric Oxide Synthase (eNOS) T-786C, 4a4b, and G894T Polymorphisms and Male Infertility: Study for Idiopathic Asthenozoospermia and Meta-Analysis

ABSTRACT Recent studies on the eNOS gene and male infertility show that expression of eNOS regulates normal spermatogenesis in the testis, and the eNOS gene variants (T-786C, 4a4b, and G894T) are potentially involved in impairment of spermatogenesis and sperm function. Thus, we conducted this association and meta-analysis study to further validate whether variants of those three loci affected the risk of idiopathic asthenozoospermia (AZS) and male infertility. Approximately 340 Chinese idiopathic AZS patients and 342 healthy men were included for this case-control study, genotyped by gel electrophoresis analysis or direct sequencing of PCR products. The eNOS mRNA isolated from the semen of patients was further examined by quantitative real-time PCR. Also, a meta-analysis of association between eNOS gene polymorphisms and male infertility was performed. A significant association was identified on allelic level between 4a4b variant and AZS in our study (chi-squared = 7.53, corrected P = 0.018, odds ratio (OR) = 1.808), while there were no significant difference of T-786C and G894T for asthenozoospermia in both genotype and allele distributions. In addition, expression of eNOS was up-regulated in patients compared with controls (about 2.4-fold, P < 0.001). Furthermore, the results of the meta-analysis support the conclusion that the T-786C and 4a4b loci were associated with male infertility in both Asian and Caucasian populations. Our study provides genetic evidence for the eNOS gene being a risk factor for idiopathic AZS and male infertility. Considering genetic differences among populations and complex pathogenesis of male infertility, more validating studies using independent samples are suggested in the future.

ITIH family genes confer risk to schizophrenia and major depressive disorder in the Han Chinese population.

As a major extracellular matrix component, ITIHs played an important role in inflammation and carcinogenesis. Several genome-wide association studies have reported that some positive signals which were derived from the tight linkage disequilibrium region on chromosome 3p21 were associated with both schizophrenia and bipolar disorders in the Caucasian population. To further investigate whether this genomic region is also a susceptibility locus of schizophrenia and major depressive disorder in the Han Chinese population, we conducted this study by recruiting 1235 schizophrenia patients, 1045 major depressive disorder patients and 1235 healthy control subjects in the Han Chinese samples for a case-control study. We genotyped seven SNPs within this region using TaqMan® technology. We found that rs2710322 was significantly associated with schizophrenia (adjusted P(allele) = 0.0018, adjusted P(genotype) = 0.006, OR [95% CI] = 1.278 [1.117-1.462]) while rs1042779 was weakly associated with schizophrenia (adjusted P(allele) = 0.048, OR [95% CI] = 1.164 [1.040-1.303]) and major depressive disorder (adjusted P(allele) = 0.042, OR [95% CI] = 1.178 [1.047-1.326]); it was also our finding that rs3821831 was positively associated with major depressive disorder (adjusted P(allele) = 0.003, adjusted P(genotype) = 0.006, OR [95% CI] = 1.426 [1.156-1.760]). Furthermore, no haplotype was found to be associated with schizophrenia and major depressive disorder. Via the association analysis which combines the schizophrenia and major depressive disorder cases, we also notice that rs1042779 and rs3821831 were significantly associated with combined cases (rs1042779: adjusted P(allele) = 0.012, adjusted P(genotype) = 0.018, OR [95% CI] = 1.171 [1.060-1.292]; rs3821831:adjusted P(genotype) = 0.012, OR [95% CI] = 1.193 [1.010-1.410]). Our results revealed that the shared genetic risk factors of both schizophrenia and major depressive disorder exist in ITIH family genes in the Han Chinese population.

The NVL gene confers risk for both major depressive disorder and schizophrenia in the Han Chinese population.

NVL (nuclear VCP (valosin containing protein)/p97-Like), a member of the AAA-ATPase (ATPases associated with various cellular activities) family, encodes a novel hTERT (human telomerase reverse transcriptase)-interacting protein NVL2 which is a telomerase component essential for holoenzyme assembly. Previous researches have reported the impacts of telomerase activity on mental illness and the potential association between NVL and major depressive disorder. To validate the susceptibility of NVL to major depressive disorder, and to investigate the overlapping risk conferred by NVL for both major depressive disorder and schizophrenia, we analyzed 9 tag single nucleotide polymorphisms (tag SNPs) using TaqMan® technology, in 1045 major depressive disorder patients, 1235 schizophrenia patients and 1235 normal controls of Han Chinese origin. We found that rs10916583 (P(allele) = 0.020, P(genotype) = 0.028, OR = 1.156) and rs16846649 (adjusted P(allele) = 0.014, P(genotype) = 0.007, OR = 0.718) were associated with major depressive disorder, while rs10916583 (adjusted P(allele) = 1.08E-02, OR = 1.213), rs16846649 (adjusted P(allele) = 7.40E-06, adjusted P(genotype) = 8.07E-05, OR = 0.598) and rs10799541 (adjusted P(allele) = 8.10E-03, adjusted P(genotype) = 0.049, OR= 0.826) showed statistically significant association with schizophrenia after Bonferroni correction. Furthermore, rs10916583 (adjusted P(allele) = 9.00E-03, adjusted P(genotype) = 3.15E-02, OR = 1.187) and rs16846649 (adjusted P(allele) = 8.92E-06, adjusted P(genotype) = 8.84E-05, OR = 0.653) remained strongly associated with the analysis of combined cases of major depressive disorder and schizophrenia after Bonferroni correction. Our results indicated that the NVL gene may contain overlapping common genetic risk factors for major depressive disorder and schizophrenia in the Han Chinese population. The roles of NVL in telomerase biogenesis were also highlighted in psychiatric pathogenesis. The study on variants conferring overlapping risk for multiple psychiatric disorders could be tangible pathogenesis support and clinical or diagnostic references.