A rare functional noncoding variant at the GWAS-implicated MIR137/MIR2682 locus might confer risk to schizophrenia and bipolar disorder.

Schizophrenia (SZ) genome-wide association studies (GWASs) have identified common risk variants in >100 susceptibility loci; however, the contribution of rare variants at these loci remains largely unexplored. One of the strongly associated loci spans MIR137 (miR137) and MIR2682 (miR2682), two microRNA genes important for neuronal function. We sequenced ∼6.9 kb MIR137/MIR2682 and upstream regulatory sequences in 2,610 SZ cases and 2,611 controls of European ancestry. We identified 133 rare variants with minor allele frequency (MAF) <0.5%. The rare variant burden in promoters and enhancers, but not insulators, was associated with SZ (p = 0.021 for MAF < 0.5%, p = 0.003 for MAF < 0.1%). A rare enhancer SNP, 1:g.98515539A>T, presented exclusively in 11 SZ cases (nominal p = 4.8 × 10(-4)). We further identified its risk allele T in 2 of 2,434 additional SZ cases, 11 of 4,339 bipolar (BP) cases, and 3 of 3,572 SZ/BP study controls and 1,688 population controls; yielding combined p values of 0.0007, 0.0013, and 0.0001 for SZ, BP, and SZ/BP, respectively. The risk allele T of 1:g.98515539A>T reduced enhancer activity of its flanking sequence by >50% in human neuroblastoma cells, predicting lower expression of MIR137/MIR2682. Both empirical and computational analyses showed weaker transcription factor (YY1) binding by the risk allele. Chromatin conformation capture (3C) assay further indicated that 1:g.98515539A>T influenced MIR137/MIR2682, but not the nearby DPYD or LOC729987. Our results suggest that rare noncoding risk variants are associated with SZ and BP at MIR137/MIR2682 locus, with risk alleles decreasing MIR137/MIR2682 expression.

MicroRNA-9 and microRNA-326 regulate human dopamine D2 receptor expression, and the microRNA-mediated expression regulation is altered by a genetic variant.

The human dopamine receptor D2 (DRD2) has been implicated in the pathophysiology of schizophrenia and other neuropsychiatric disorders. Most antipsychotic drugs influence dopaminergic transmission through blocking dopamine receptors, primarily DRD2. We report here the post-transcriptional regulation of DRD2 expression by two brain-expressed microRNAs (miRs), miR-326 and miR-9, in an ex vivo mode, and show the relevance of miR-mediated DRD2 expression regulation in human dopaminergic neurons and in developing human brains. Both miRs targeted the 3'-UTR (untranslated region) of DRD2 in NT2 (neuron-committed teratocarcinoma, which endogenously expresses DRD2) and CHO (Chinese hamster ovary) cell lines, decreasing luciferase activity measured by a luciferase reporter gene assay. miR-326 overexpression reduced DRD2 mRNA and DRD2 receptor synthesis. Both antisense miR-326 and antisense miR-9 increased DRD2 protein abundance, suggesting an endogenous repression of DRD2 expression by both miRs. Furthermore, a genetic variant (rs1130354) within the DRD2 3'-UTR miR-targeting site interferes with miR-326-mediated repression of DRD2 expression. Finally, co-expression analysis identified an inverse correlation of DRD2 expression with both miR-326 and miR-9 in differentiating dopaminergic neurons derived from human induced pluripotent stem cells (iPSCs) and in developing human brain regions implicated in schizophrenia. Our study provides empirical evidence suggesting that miR-326 and miR-9 may regulate dopaminergic signaling, and miR-326 and miR-9 may be considered as potential drug targets for the treatment of disorders involving abnormal DRD2 function, such as schizophrenia.

Transcriptome study of differential expression in schizophrenia.

Schizophrenia genome-wide association studies (GWAS) have identified common SNPs, rare copy number variants (CNVs) and a large polygenic contribution to illness risk, but biological mechanisms remain unclear. Bioinformatic analyses of significantly associated genetic variants point to a large role for regulatory variants. To identify gene expression abnormalities in schizophrenia, we generated whole-genome gene expression profiles using microarrays on lymphoblastoid cell lines (LCLs) from 413 cases and 446 controls. Regression analysis identified 95 transcripts differentially expressed by affection status at a genome-wide false discovery rate (FDR) of 0.05, while simultaneously controlling for confounding effects. These transcripts represented 89 genes with functions such as neurotransmission, gene regulation, cell cycle progression, differentiation, apoptosis, microRNA (miRNA) processing and immunity. This functional diversity is consistent with schizophrenia's likely significant pathophysiological heterogeneity. The overall enrichment of immune-related genes among those differentially expressed by affection status is consistent with hypothesized immune contributions to schizophrenia risk. The observed differential expression of extended major histocompatibility complex (xMHC) region histones (HIST1H2BD, HIST1H2BC, HIST1H2BH, HIST1H2BG and HIST1H4K) converges with the genetic evidence from GWAS, which find the xMHC to be the most significant susceptibility locus. Among the differentially expressed immune-related genes, B3GNT2 is implicated in autoimmune disorders previously tied to schizophrenia risk (rheumatoid arthritis and Graves' disease), and DICER1 is pivotal in miRNA processing potentially linking to miRNA alterations in schizophrenia (e.g. MIR137, the second strongest GWAS finding). Our analysis provides novel candidate genes for further study to assess their potential contribution to schizophrenia.

No significant association of 14 candidate genes with schizophrenia in a large European ancestry sample: implications for psychiatric genetics.

OBJECTIVE: The authors carried out a genetic association study of 14 schizophrenia candidate genes (RGS4, DISC1, DTNBP1, STX7, TAAR6, PPP3CC, NRG1, DRD2, HTR2A, DAOA, AKT1, CHRNA7, COMT, and ARVCF). This study tested the hypothesis of association of schizophrenia with common single nucleotide polymorphisms (SNPs) in these genes using the largest sample to date that has been collected with uniform clinical methods and the most comprehensive set of SNPs in each gene. METHOD: The sample included 1,870 cases (schizophrenia and schizoaffective disorder) and 2,002 screened comparison subjects (i.e. controls), all of European ancestry, with ancestral outliers excluded based on analysis of ancestry-informative markers. The authors genotyped 789 SNPs, including tags for most common SNPs in each gene, SNPs previously reported as associated, and SNPs located in functional domains of genes such as promoters, coding exons (including nonsynonymous SNPs), 3' untranslated regions, and conserved noncoding sequences. After extensive data cleaning, 648 SNPs were analyzed for association of single SNPs and of haplotypes. RESULTS: Neither experiment-wide nor gene-wide statistical significance was observed in the primary single-SNP analyses or in secondary analyses of haplotypes or of imputed genotypes for additional common HapMap SNPs. Results in SNPs previously reported as associated with schizophrenia were consistent with chance expectation, and four functional polymorphisms in COMT, DRD2, and HTR2A did not produce nominally significant evidence to support previous evidence for association. CONCLUSIONS: It is unlikely that common SNPs in these genes account for a substantial proportion of the genetic risk for schizophrenia, although small effects cannot be ruled out.

Transcriptomic signatures of schizophrenia revealed by dopamine perturbation in an ex vivo model.

The dopaminergic hypothesis of schizophrenia (SZ) postulates that dopaminergic over activity causes psychosis, a central feature of SZ, based on the observation that blocking dopamine (DA) improves psychotic symptoms. DA is known to have both receptor- and non-receptor-mediated effects, including oxidative mechanisms that lead to apoptosis. The role of DA-mediated oxidative processes in SZ has been little studied. Here, we have used a cell perturbation approach and measured transcriptomic profiles by RNAseq to study the effect of DA exposure on transcription in B-cell transformed lymphoblastoid cell lines (LCLs) from 514 SZ cases and 690 controls. We found that DA had widespread effects on both cell growth and gene expression in LCLs. Overall, 1455 genes showed statistically significant differential DA response in SZ cases and controls. This set of differentially expressed genes is enriched for brain expression and for functions related to immune processes and apoptosis, suggesting that DA may play a role in SZ pathogenesis through modulating those systems. Moreover, we observed a non-significant enrichment of genes near genome-wide significant SZ loci and with genes spanned by SZ-associated copy number variants (CNVs), which suggests convergent pathogenic mechanisms detected by both genetic association and gene expression. The study suggests a novel role of DA in the biological processes of immune and apoptosis that may be relevant to SZ pathogenesis. Furthermore, our results show the utility of pathophysiologically relevant perturbation experiments to investigate the biology of complex mental disorders.

A novel carbohydrate-based therapeutic GCS-100 overcomes bortezomib resistance and enhances dexamethasone-induced apoptosis in multiple myeloma cells.

Human multiple myeloma is a presently incurable hematologic malignancy, and novel biologically based therapies are urgently needed. GCS-100 is a polysaccharide derived from citrus pectin in clinical development for the treatment of cancer. Here we show that GCS-100 induces apoptosis in various multiple myeloma cell lines, including those resistant to dexamethasone, melphalan, or doxorubicin. Examination of purified patient multiple myeloma cells showed similar results. Specifically, GCS-100 decreases viability of bortezomib/PS-341-resistant multiple myeloma patient cells. Importantly, GCS-100 inhibits multiple myeloma cell growth induced by adhesion to bone marrow stromal cells; overcome the growth advantage conferred by antiapoptotic protein Bcl-2, heat shock protein-27, and nuclear factor-kappaB; and blocks vascular endothelial growth factor-induced migration of multiple myeloma cells. GCS-100-induced apoptosis is associated with activation of caspase-8 and caspase-3 followed by proteolytic cleavage of poly(ADP-ribose) polymerase enzyme. Combined with dexamethasone, GCS-100 induces additive anti-multiple myeloma cytotoxicity associated with mitochondrial apoptotic signaling via release of cytochrome c and Smac followed by activation of caspase-3. Moreover, GCS-100 + dexamethasone-induced apoptosis in multiple myeloma cells is accompanied by a marked inhibition of an antiapoptotic protein Galectin-3, without significant alteration in Bcl-2 expression. Collectively, these findings provide the framework for clinical evaluation of GCS-100, either alone or in combination with dexamethasone, to inhibit tumor growth, overcome drug resistance, and improve outcome for patients with this universally fatal hematologic malignancy.

Maximum Entropy for the International Division of Labor.

As a result of the international division of labor, the trade value distribution on different products substantiated by international trade flows can be regarded as one country's strategy for competition. According to the empirical data of trade flows, countries may spend a large fraction of export values on ubiquitous and competitive products. Meanwhile, countries may also diversify their exports share on different types of products to reduce the risk. In this paper, we report that the export share distribution curves can be derived by maximizing the entropy of shares on different products under the product's complexity constraint once the international market structure (the country-product bipartite network) is given. Therefore, a maximum entropy model provides a good fit to empirical data. The empirical data is consistent with maximum entropy subject to a constraint on the expected value of the product complexity for each country. One country's strategy is mainly determined by the types of products this country can export. In addition, our model is able to fit the empirical export share distribution curves of nearly every country very well by tuning only one parameter.

Genome-wide survey of interindividual differences of RNA stability in human lymphoblastoid cell lines.

The extent to which RNA stability differs between individuals and its contribution to the interindividual expression variation remain unknown. We conducted a genome-wide analysis of RNA stability in seven human HapMap lymphoblastoid cell lines (LCLs) and analyzed the effect of DNA sequence variation on RNA half-life differences. Twenty-six percent of the expressed genes exhibited RNA half-life differences between LCLs at a false discovery rate (FDR) < 0.05, which accounted for ~ 37% of the gene expression differences between individuals. Nonsense polymorphisms were associated with reduced RNA half-lives. In genes presenting interindividual RNA half-life differences, higher coding GC3 contents (G and C percentages at the third-codon positions) were correlated with increased RNA half-life. Consistently, G and C alleles of single nucleotide polymorphisms (SNPs) in protein coding sequences were associated with enhanced RNA stability. These results suggest widespread interindividual differences in RNA stability related to DNA sequence and composition variation.