mRNA transcript diversity creates new opportunities for pharmacological intervention

SARS-CoV-2 infection severity and mortality is modulated by repeat-mediated regulation of alternative splicing

Like single-stranded RNA viruses, SARS-CoV-2 hijacks the host transcriptional machinery for its own replication. Numerous traditional differential gene expression-based investigations have examined the diverse clinical symptoms caused by SARS-CoV-2 infection. The virus, on the other hand, also affects the host splicing machinery, causing host transcriptional dysregulation, which can lead to diverse clinical outcomes. Hence, in this study, we performed host transcriptome sequencing of 125 hospital-admitted COVID-19 patients to understand the transcriptomic differences between the severity sub-phenotypes of mild, moderate, severe, and mortality. We performed transcript-level differential expression analysis, investigated differential isoform usage, looked at the splicing patterns within the differentially expressed transcripts (DET), and elucidated the possible genome regulatory features. Our DTE analysis showed evidence of diminished transcript length and diversity as well as altered promoter site usage in the differentially expressed protein-coding transcripts in the COVID-19 mortality patients. We also investigated the potential mechanisms driving the alternate splicing and discovered a compelling differential enrichment of repeats in the promoter region and a specific enrichment of SINE (Alu) near the splicing sites of differentially expressed transcripts. These findings suggested a repeat-mediated plausible regulation of alternative splicing as a potential modulator of COVID-19 disease severity. In this work, we emphasize the role of scarcely elucidated functional role of alternative splicing in influencing COVID-19 disease severity sub-phenotypes, clinical outcomes, and its putative mechanism. IMPORTANCE The wide range of clinical symptoms reported during the COVID-19 pandemic inherently highlights the numerous factors that influence the progression and prognosis of SARS-CoV-2 infection. While several studies have investigated the host response and discovered immunological dysregulation during severe infection, most of them have the common theme of focusing only up to the gene level. Viruses, especially RNA viruses, are renowned for hijacking the host splicing machinery for their own proliferation, which inadvertently puts pressure on the host transcriptome, exposing another side of the host response to the pathogen challenge. Therefore, in this study, we examine host response at the transcript-level to discover a transcriptional difference that culminates in differential gene-level expression. Importantly, this study highlights diminished transcript diversity and possible regulation of transcription by differentially abundant repeat elements near the promoter region and splicing sites in COVID-19 mortality patients, which together with differentially expressed isoforms hold the potential to elaborate disease severity and outcome.

Network-based method for drug target discovery at the isoform level

Identification of primary targets associated with phenotypes can facilitate exploration of the underlying molecular mechanisms of compounds and optimization of the structures of promising drugs. However, the literature reports limited effort to identify the target major isoform of a single known target gene. The majority of genes generate multiple transcripts that are translated into proteins that may carry out distinct and even opposing biological functions through alternative splicing. In addition, isoform expression is dynamic and varies depending on the developmental stage and cell type. To identify target major isoforms, we integrated a breast cancer type-specific isoform coexpression network with gene perturbation signatures in the MCF7 cell line in the Connectivity Map database using the ‘shortest path’ drug target prioritization method. We used a leukemia cancer network and differential expression data for drugs in the HL-60 cell line to test the robustness of the detection algorithm for target major isoforms. We further analyzed the properties of target major isoforms for each multi-isoform gene using pharmacogenomic datasets, proteomic data and the principal isoforms defined by the APPRIS and STRING datasets. Then, we tested our predictions for the most promising target major protein isoforms of DNMT1, MGEA5 and P4HB4 based on expression data and topological features in the coexpression network. Interestingly, these isoforms are not annotated as principal isoforms in APPRIS. Lastly, we tested the affinity of the target major isoform of MGEA5 for streptozocin through in silico docking. Our findings will pave the way for more effective and targeted therapies via studies of drug targets at the isoform level.

A novel role of U1 snRNP: Splice site selection from a distance

Removal of introns by pre-mRNA splicing is fundamental to gene function in eukaryotes. However, understanding the mechanism by which exon-intron boundaries are defined remains a challenging endeavor. Published reports support that the recruitment of U1 snRNP at the 5'ss marked by GU dinucleotides defines the 5'ss as well as facilitates 3'ss recognition through cross-exon interactions. However, exceptions to this rule exist as U1 snRNP recruited away from the 5'ss retains the capability to define the splice site, where the cleavage takes place. Independent reports employing exon 7 of Survival Motor Neuron (SMN) genes suggest a long-distance effect of U1 snRNP on splice site selection upon U1 snRNP recruitment at target sequences with or without GU dinucleotides. These findings underscore that sequences distinct from the 5'ss may also impact exon definition if U1 snRNP is recruited to them through partial complementarity with the U1 snRNA. In this review we discuss the expanded role of U1 snRNP in splice-site selection due to U1 ability to be recruited at more sites than predicted solely based on GU dinucleotides.

Voltage-Dependent Sarcolemmal Ion Channel Abnormalities in the Dystrophin-Deficient Heart

Mutations in the gene encoding for the intracellular protein dystrophin cause severe forms of muscular dystrophy. These so-called dystrophinopathies are characterized by skeletal muscle weakness and degeneration. Dystrophin deficiency also gives rise to considerable complications in the heart, including cardiomyopathy development and arrhythmias. The current understanding of the pathomechanisms in the dystrophic heart is limited, but there is growing evidence that dysfunctional voltage-dependent ion channels in dystrophin-deficient cardiomyocytes play a significant role. Herein, we summarize the current knowledge about abnormalities in voltage-dependent sarcolemmal ion channel properties in the dystrophic heart, and discuss the potentially underlying mechanisms, as well as their pathophysiological relevance.

Gene isoforms as expression-based biomarkers predictive of drug response in vitro

Next-generation sequencing technologies have recently been used in pharmacogenomic studies to characterize large panels of cancer cell lines at the genomic and transcriptomic levels. Among these technologies, RNA-sequencing enable profiling of alternatively spliced transcripts. Given the high frequency of mRNA splicing in cancers, linking this feature to drug response will open new avenues of research in biomarker discovery. To identify robust transcriptomic biomarkers for drug response across studies, we develop a meta-analytical framework combining the pharmacological data from two large-scale drug screening datasets. We use an independent pan-cancer pharmacogenomic dataset to test the robustness of our candidate biomarkers across multiple cancer types. We further analyze two independent breast cancer datasets and find that specific isoforms of IGF2BP2, NECTIN4, ITGB6, and KLHDC9 are significantly associated with AZD6244, lapatinib, erlotinib, and paclitaxel, respectively. Our results support isoform expressions as a rich resource for biomarkers predictive of drug response.

Advancing psychiatric pharmacogenomics using drug development paradigms

Drugs used to treat psychiatric disorders, even when taken as directed, fail to provide adequate relief for a sizeable proportion of patients. Despite our advancements in understanding human genetics and development of high-throughput tools to probe variation, pharmacogenomics has yielded marginal ability to predict drug response for psychiatric disorders. Here, I review the current pharmacogenomics paradigm, identifying opportunities to incorporate drug development strategies designed to increase the probability of delivering a successful molecule to the clinic. This includes using in-depth pharmacokinetic profiles, clear measures of target engagement and target-specific pharmacodynamic responses orthogonal to clinical response. The complex pharmacological profiles psychiatric drugs require re-examination of simplified clinical response-oriented pharmacogenetic hypotheses, in favor of a more complete patient profile.

Genome annotation for clinical genomic diagnostics: strengths and weaknesses

The Human Genome Project and advances in DNA sequencing technologies have revolutionized the identification of genetic disorders through the use of clinical exome sequencing. However, in a considerable number of patients, the genetic basis remains unclear. As clinicians begin to consider whole-genome sequencing, an understanding of the processes and tools involved and the factors to consider in the annotation of the structure and function of genomic elements that might influence variant identification is crucial. Here, we discuss and illustrate the strengths and weaknesses of approaches for the annotation and classification of important elements of protein-coding genes, other genomic elements such as pseudogenes and the non-coding genome, comparative-genomic approaches for inferring gene function, and new technologies for aiding genome annotation, as a practical guide for clinicians when considering pathogenic sequence variation. Complete and accurate annotation of structure and function of genome features has the potential to reduce both false-negative (from missing annotation) and false-positive (from incorrect annotation) errors in causal variant identification in exome and genome sequences. Re-analysis of unsolved cases will be necessary as newer technology improves genome annotation, potentially improving the rate of diagnosis.

Hiding in the Shadows: CPOX Expression and 5-ALA Induced Fluorescence in Human Glioma Cells

Protoporphyrin IX (PpIX) is widely used in photodynamic diagnosis. To date, the details of molecular mechanisms underlying PpIX accumulation in malignant cells after 5-ALA administration remain unclear. The fluorescence of PpIX was studied in human glioma cells. Several cell cultures were established from glioma tumor tissue to study the differences between fluorescence-positive and fluorescence-negative human glioma tumors. The cell cultures demonstrated fluorescence profiles similar to those of source tumor tissues, which allows us to use these cultures in experimental research. Dynamics of the rates of synthesis and degradation of fluorescent protoporphyrin IX was studied in the cultures obtained. In addition, the expression of CPOX, an enzyme involved in PpIX synthesis, was evaluated. mRNA levels of heme biosynthesis enzymes were analyzed, and PpIX fluorescence proved to correlate with the CPOX protein level, whereas no such correlation was observed at the mRNA level. Fluorescence intensity decreased at low levels of the enzyme, which indicates its critical role in PpIX fluorescence. Finally, the fluorescence intensity proved to correlate with the proliferative activity.

Regulation and Functional Implications of Opioid Receptor Splicing in Opioid Pharmacology and HIV Pathogenesis

Despite the identification and characterization of four opioid receptor subtypes and the genes from which they are encoded, pharmacological data does not conform to the predications of a four opioid receptor model. Instead, current studies of opioid pharmacology suggest the existence of additional receptor subtypes; however, no additional opioid receptor subtype has been identified to date. It is now understood that this discrepancy is due to the generation of multiple isoforms of opioid receptor subtypes. While several mechanisms are utilized to generate these isoforms, the primary mechanism involves alternative splicing of the pre-mRNA transcript. Extensive alternative splicing patterns for opioid receptors have since been identified and discrepancies in opioid pharmacology are now partially attributed to variable expression of these isoforms. Recent studies have been successful in characterizing the localization of these isoforms as well as their specificity in ligand binding; however, the regulation of opioid receptor splicing specificity is poorly characterized. Furthermore, the functional significance of individual receptor isoforms and the extent to which opioid- and/or HIV-mediated changes in the opioid receptor isoform profile contributes to altered opioid pharmacology or the well-known physiological role of opioids in the exacerbation of HIV neurocognitive dysfunction is unknown. As such, the current review details constitutive splicing mechanisms as well as the specific architecture of opioid receptor genes, transcripts, and receptors in order to highlight the current understanding of opioid receptor isoforms, potential mechanisms of their regulation and signaling, and their functional significance in both opioid pharmacology and HIV-associated neuropathology.

A reliable method for quantification of splice variants using RT-qPCR

Background

The majority of protein isoforms arise from alternative splicing of the encoding primary RNA transcripts. To understand the significance of single splicing events, reliable techniques are needed to determine their incidence. However, existing methods are labour-intensive, error-prone or of limited use.

Results

Here, we present an improved method to determine the relative incidence of transcripts that arise from alternative splicing at a single site. Splice variants were quantified within a single sample using one-step reverse transcription quantitative PCR. Amplification products obtained with variant specific primer pairs were compared to those obtained with primer pairs common to both variants. The identities of variant specific amplicons were simultaneously verified by melt curve analysis. Independent calculations of the relative incidence of each variant were performed. Since the relative incidences of variants have to add upto 100 %, the method provides an internal control to monitor experimental errors and uniform reverse transcription. The reliability of the method was tested using mixtures of cDNA templates as well as RNA samples from different sources.

Conclusion

The method described here, is easy to set up and does not need unrelated reference genes and time consuming, error-prone standard curves. It provides a reliable and precise technique to distinguish small differences of the relative incidence of two splice variants.

Experimental approaches to studying the nature and impact of splicing variation in zebrafish

From a fixed number of genes carried in all cells, organisms create considerable diversity in cellular phenotype through differential regulation of gene expression. One prevalent source of transcriptome diversity is alternative pre-mRNA splicing, which is manifested in many different forms. Zebrafish models of splicing dysfunction due to mutated spliceosome components provide opportunity to link biochemical analyses of spliceosome structure and function with whole organism phenotypic outcomes. Drawing from experience with two zebrafish mutants: cephalophŏnus (a prpf8 mutant, isolated for defects in granulopoiesis) and caliban (a rnpc3 mutant, isolated for defects in digestive organ development), we describe the use of glycerol gradient sedimentation and native gel electrophoresis to resolve components of aberrant splicing complexes. We also describe how RNAseq can be employed to examine relatively rare alternative splicing events including intron retention. Such experimental approaches in zebrafish can promote understanding of how splicing variation and dysfunction contribute to phenotypic diversity and disease pathogenesis.

Recent advances in therapeutic strategies that focus on the regulation of ion channel expression

A number of different ion channel types are involved in cell signaling networks, and homeostatic regulatory mechanisms contribute to the control of ion channel expression. Profiling of global gene expression using microarray technology has recently provided novel insights into the molecular mechanisms underlying the homeostatic and pathological control of ion channel expression. It has demonstrated that the dysregulation of ion channel expression is associated with the pathogenesis of neural, cardiovascular, and immune diseases as well as cancers. In addition to the transcriptional, translational, and post-translational regulation of ion channels, potentially important evidence on the mechanisms controlling ion channel expression has recently been accumulated. The regulation of alternative pre-mRNA splicing is therefore a novel therapeutic strategy for the treatment of dominant-negative splicing disorders. Epigenetic modification plays a key role in various pathological conditions through the regulation of pluripotency genes. Inhibitors of pre-mRNA splicing and histone deacetyalase/methyltransferase have potential as potent therapeutic drugs for cancers and autoimmune and inflammatory diseases. Moreover, membrane-anchoring proteins, lysosomal and proteasomal degradation-related molecules, auxiliary subunits, and pharmacological agents alter the protein folding, membrane trafficking, and post-translational modifications of ion channels, and are linked to expression-defect channelopathies. In this review, we focused on recent insights into the transcriptional, spliceosomal, epigenetic, and proteasomal regulation of ion channel expression: Ca(2+) channels (TRPC/TRPV/TRPM/TRPA/Orai), K(+) channels (voltage-gated, KV/Ca(2+)-activated, KCa/two-pore domain, K2P/inward-rectifier, Kir), and Ca(2+)-activated Cl(-) channels (TMEM16A/TMEM16B). Furthermore, this review highlights expression of these ion channels in expression-defect channelopathies.

Transcriptomic variation of pharmacogenes in multiple human tissues and lymphoblastoid cell lines

Variation in the expression level and activity of genes involved in drug disposition and action (‘pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20–45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Hypertension pharmacogenomics: in search of personalized treatment approaches

Cardiovascular and renal diseases are associated with many risk factors, of which hypertension is one of the most prevalent. Worldwide, blood pressure control is only achieved in ∼50% of those treated for hypertension, despite the availability of a considerable number of antihypertensive drugs from different pharmacological classes. Although many reasons exist for poor blood pressure control, a likely contributor is the inability to predict to which antihypertensive drug an individual is most likely to respond. Hypertension pharmacogenomics and other 'omics' technologies have the potential to identify genetic signals that are predictive of response or adverse outcome to particular drugs, and guide selection of hypertension treatment for a given individual. Continued research in this field will enhance our understanding of how to maximally deploy the various antihypertensive drug classes to optimize blood pressure response at the individual level. This Review summarizes the available literature on the most convincing genetic signals associated with antihypertensive drug responses and adverse cardiovascular outcomes. Future research in this area will be facilitated by enhancing collaboration between research groups through consortia such as the International Consortium for Antihypertensives Pharmacogenomics Studies, with the goal of translating replicated findings into clinical implementation.

Genomic structure and expression of the human serotonin 2A receptor gene (HTR2A) locus: identification of novel HTR2A and antisense (HTR2A-AS1) exons

Background

The serotonin 2A receptor is widely implicated in genetic association studies and remains an important drug target for psychiatric, neurological, and cardiovascular conditions. RNA sequencing redefined the architecture of the serotonin 2A receptor gene (HTR2A), revealing novel mRNA transcript isoforms utilizing unannotated untranslated regions of the gene. Expression of these untranslated regions is modulated by common single nucleotide polymorphisms (SNPs), namely rs6311. Previous studies did not fully capture the complexity of the sense- and antisense-encoded transcripts with respect to novel exons in the HTR2A gene locus. Here, we comprehensively catalogued exons and RNA isoforms for both HTR2A and HTR2A-AS1 using RNA-Seq from human prefrontal cortex and multiple mouse tissues. We subsequently tested associations between expression of newfound gene features and common SNPs in humans.

Results

We find that the human HTR2A gene spans ~66 kilobases and consists of 7, rather than 4 exons. Furthermore, the revised human HTR2A-AS1 gene spans ~474 kilobases and consists of 18, rather than 3 exons. Three HTR2A exons directly overlap with HTR2A-AS1 exons, suggesting potential for complementary nucleotide interactions. The repertoire of possible mouse Htr2a splice isoforms is remarkably similar to humans and we also find evidence for overlapping sense-antisense transcripts in the same relative positions as the human transcripts. rs6311 and SNPs in high linkage disequilibrium are associated with HTR2A-AS1 expression, in addition to previously described associations with expression of the extended 5’ untranslated region of HTR2A.

Conclusions

Our proposed HTR2A and HTR2A-AS1 gene structures dramatically differ from current annotations, now including overlapping exons on the sense and anti-sense strands. We also find orthologous transcript isoforms expressed in mice, providing opportunities to elucidate the biological roles of the human isoforms using a model system. Associations between rs6311 and expression of HTR2A and HTR2A-AS1 suggest this polymorphism is capable of modulating the expression of the sense or antisense transcripts. Still unclear is whether these SNPs act directly on the expression of the sense or antisense transcripts and whether overlapping exons are capable of interacting through complimentary base-pairing. Additional studies are necessary to determine the extent and nature of interactions between the SNPs and the transcripts prior to interpreting these findings in the context of phenotypes associated with HTR2A.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0325-6) contains supplementary material, which is available to authorized users.

Molecular identification of the dominant-negative, splicing isoform of the two-pore domain K(+) channel K(2P)5.1 in lymphoid cells and enhancement of its expression by splicing inhibition

The two-pore domain background K(+) channel K2P5.1 is expected as a possible therapeutic target for autoimmune and inflammatory disorders and cancers because it plays an important role in maintaining the resting membrane potential and regulation of Ca(2+) signaling in T lymphocytes and cancer cells. However, the lack of selective K2P5.1 blockers has led to difficulties conducting experimental studies on this K(+) channel. We identified a novel splicing isoform of K2P5.1, K2P5.1B from the mammalian spleen, which lacked the N-terminus of full-length K2P5.1A. A co-immunoprecipitation assay using mice spleen lysates revealed an interaction between K2P5.1A and K2P5.1B in the cytoplasmic C-terminal domain. In a heterologous HEK293 expression system, K2P5.1B inhibited the trafficking of K2P5.1A to the plasma membrane. The alkaline pHe-induced hyperpolarizing response was significantly suppressed in K2P5.1B-transfected human leukemia K562 cells. Enhancement in cell proliferation by the overexpression of K2P5.1A in K562 was significantly prevented by the transfection of K2P5.1B. The spliceosome inhibitor pladienolide B significantly enhanced the relative expression of K2P5.1B in K562, resulting in decreases in the activity of K2P5.1A. K2P5.1B suppresses the function of the K2P5.1 K(+) channel in a dominant-negative manner, suggesting that the mRNA splicing mechanisms underlying the transcriptional regulation of K2P5.1B may be a new therapeutic strategy for autoimmune and inflammatory disorders and cancers.

Cosplicing network analysis of mammalian brain RNA-Seq data utilizing WGCNA and Mantel correlations

Across species and tissues and especially in the mammalian brain, production of gene isoforms is widespread. While gene expression coordination has been previously described as a scale-free coexpression network, the properties of transcriptome-wide isoform production coordination have been less studied. Here we evaluate the system-level properties of cosplicing in mouse, macaque, and human brain gene expression data using a novel network inference procedure. Genes are represented as vectors/lists of exon counts and distance measures sensitive to exon inclusion rates quantifies differences across samples. For all gene pairs, distance matrices are correlated across samples, resulting in cosplicing or cotranscriptional network matrices. We show that networks including cosplicing information are scale-free and distinct from coexpression. In the networks capturing cosplicing we find a set of novel hubs with unique characteristics distinguishing them from coexpression hubs: heavy representation in neurobiological functional pathways, strong overlap with markers of neurons and neuroglia, long coding lengths, and high number of both exons and annotated transcripts. Further, the cosplicing hubs are enriched in genes associated with autism spectrum disorders. Cosplicing hub homologs across eukaryotes show dramatically increasing intronic lengths but stable coding region lengths. Shared transcription factor binding sites increase coexpression but not cosplicing; the reverse is true for splicing-factor binding sites. Genes with protein-protein interactions have strong coexpression and cosplicing. Additional factors affecting the networks include shared microRNA binding sites, spatial colocalization within the striatum, and sharing a chromosomal folding domain. Cosplicing network patterns remain relatively stable across species.

The current and potential impact of genetics and genomics on neuropsychopharmacology

One justification for the major scientific and financial investments in genetic and genomic studies in medicine is their therapeutic potential, both for revealing novel targets for drugs which treat the disease process, as well as allowing for more effective and safe use of existing medications. This review considers the extent to which this promise has yet been realised within psychopharmacology, how things are likely to develop in the foreseeable future, and the key issues involved. It draws primarily on examples from schizophrenia and its treatments. One observation is that there is evidence for a range of genetic influences on different aspects of psychopharmacology in terms of discovery science, but far less evidence that meets the standards required before such discoveries impact upon clinical practice. One reason is that results reveal complex genetic influences that are hard to replicate and usually of very small effect. Similarly, the slow progress being made in revealing the genes that underlie the major psychiatric syndromes hampers attempts to apply the findings to identify novel drug targets. Nevertheless, there are some intriguing positive findings of various kinds, and clear potential for genetics and genomics to play an increasing and major role in psychiatric drug discovery.

Pharmacologic management of Duchenne muscular dystrophy: target identification and preclinical trials

Duchenne muscular dystrophy (DMD) is an X-linked human disorder in which absence of the protein dystrophin causes degeneration of skeletal and cardiac muscle. For the sake of treatment development, over and above definitive genetic and cell-based therapies, there is considerable interest in drugs that target downstream disease mechanisms. Drug candidates have typically been chosen based on the nature of pathologic lesions and presumed underlying mechanisms and then tested in animal models. Mammalian dystrophinopathies have been characterized in mice (mdx mouse) and dogs (golden retriever muscular dystrophy [GRMD]). Despite promising results in the mdx mouse, some therapies have not shown efficacy in DMD. Although the GRMD model offers a higher hurdle for translation, dogs have primarily been used to test genetic and cellular therapies where there is greater risk. Failed translation of animal studies to DMD raises questions about the propriety of methods and models used to identify drug targets and test efficacy of pharmacologic intervention. The mdx mouse and GRMD dog are genetically homologous to DMD but not necessarily analogous. Subcellular species differences are undoubtedly magnified at the whole-body level in clinical trials. This problem is compounded by disparate cultures in clinical trials and preclinical studies, pointing to a need for greater rigor and transparency in animal experiments. Molecular assays such as mRNA arrays and genome-wide association studies allow identification of genetic drug targets more closely tied to disease pathogenesis. Genes in which polymorphisms have been directly linked to DMD disease progression, as with osteopontin, are particularly attractive targets.

Recent genetic findings in schizophrenia and their therapeutic relevance

Over 100 loci are now associated with schizophrenia risk as identified by single nucleotide polymorphisms (SNPs) in genome-wide association studies. These findings mean that 'genes for schizophrenia' have unquestionably been found. However, many questions remain unanswered, including several which affect their therapeutic significance. The SNPs individually have minor effects, and even cumulatively explain only a modest fraction of the genetic predisposition. The remainder likely results from many more loci, from rare variants, and from gene-gene and gene-environment interactions. The risk SNPs are almost all non-coding, meaning that their biological significance is unclear; probably their effects are mediated via an influence on gene regulation, and emerging evidence suggests that some key molecular events occur during early brain development. The loci include novel genes of unknown function as well as genes and pathways previously implicated in the pathophysiology of schizophrenia, e.g. NMDA receptor signalling. Genes in the latter category have the clearer therapeutic potential, although even this will be a challenging process because of the many complexities concerning the genetic architecture and mediating mechanisms. This review summarises recent schizophrenia genetic findings and some key issues they raise, particularly with regard to their implications for identifying and validating novel drug targets.

Expression of ZNF804A in human brain and alterations in schizophrenia, bipolar disorder, and major depressive disorder: a novel transcript fetally regulated by the psychosis risk variant rs1344706

IMPORTANCE

The single-nucleotide polymorphism rs1344706 in the zinc finger protein 804A gene (ZNF804A) shows genome-wide association with schizophrenia and bipolar disorder. Little is known regarding the expression of ZNF804A and the functionality of rs1344706.

OBJECTIVES

To characterize ZNF804A expression in human brain and to investigate how it changes across the life span and how it is affected by rs1344706, schizophrenia, bipolar disorder, and major depressive disorder.

DESIGN, SETTING, AND PARTICIPANTS

Molecular and immunochemical methods were used to study ZNF804A messenger RNA (mRNA) and ZNF804A protein, respectively. ZNF804A transcripts were investigated using next-generation sequencing and polymerase chain reaction-based methods, and ZNF804A protein was investigated using Western blots and immunohistochemistry. Samples of dorsolateral prefrontal cortex and inferior parietal lobe tissue were interrogated from 697 participants between 14 weeks' gestational age and age 85 years, including patients with schizophrenia, bipolar disorder, or major depressive disorder.

MAIN OUTCOMES AND MEASURES

Quantitative measurements of ZNF804A mRNA and immunoreactivity, and the effect of diagnosis and rs1344706 genotype.

RESULTS

ZNF804A was expressed across the life span, with highest expression prenatally. An abundant and developmentally regulated truncated ZNF804A transcript was identified, missing exons 1 and 2 (ZNF804AE3E4) and predicted to encode a protein lacking the zinc finger domain. rs1344706 influenced expression of ZNF804AE3E4 mRNA in fetal brain (P = .02). In contrast, full-length ZNF804A showed no association with genotype (P > .05). ZNF804AE3E4 mRNA expression was decreased in patients with schizophrenia (P = .006) and increased in those with major depressive disorder (P < .001), and there was a genotype-by-diagnosis interaction in bipolar disorder (P = .002). ZNF804A immunoreactivity was detected in fetal and adult human cerebral cortex. It was localized primarily to pyramidal neurons, with cytoplasmic as well as dendritic and nuclear staining. No differences in ZNF804A-immunoreactive neurons were seen in schizophrenia or related to rs1344706 (P > .05).

CONCLUSIONS AND RELEVANCE

rs1344706 influences the expression of ZNF804AE3E4, a novel splice variant. The effect is limited to fetal brain and to this isoform. It may be part of the mechanism by which allelic variation in ZNF804A affects risk of psychosis. ZNF804A is translated in human brain, where its functions may extend beyond its predicted role as a transcription factor.

Altered gamma-aminobutyric acid type B receptor subunit 1 splicing in alcoholics

BACKGROUND

Chronic alcohol exposure can change splice variant expression. The gamma-aminobutyric acid type B (GABAB) receptor undergoes splicing and is an alcoholism treatment target, but there is little information about splicing changes in this receptor in alcoholics. We studied GABAB receptor subunit 1 (GABAB1) splicing in alcoholic postmortem brains.

METHODS

To maximize GABAB1 splice junction identification, we combined gene specific libraries with RNA-seq. Splice junctions and mapped reads were also found from intronic and intergenic regions. We compared GABAB1 splice junctions in prefrontal cortices from 14 alcoholic and 15 control subjects and introduced new strategies, reads per kilobase of splice junction model per million mapped reads and reads per kilobase of gene model per million mapped reads, for quantitating splice junction and gene expression.

RESULTS

Novel splice junction detection indicated that the GABAB1 gene is at least two times longer than the previously reported gene length. GABAB1 exon and intron expression data showed low expression at the 5' end exons and exon grouping. This indicated that there are short splicing variants in addition to GABAB receptor subunit GABAB1a, the longest known major transcript. We found that chronic alcohol altered exon/intron expression and splice junction levels. Decreased expression of the gamma-aminobutyric acid binding site, a transmembrane domain and a microRNA binding site may decrease normal GABAB1 transcript population and thereby decrease normal signal transduction in alcoholics.

CONCLUSIONS

We discovered novel, complex splicing of GABAB1 in human brain and showed that chronic alcohol produces additional splicing complexity.

Coexpression and cosplicing network approaches for the study of mammalian brain transcriptomes

Next-generation sequencing experiments have demonstrated great potential for transcriptome profiling. While transcriptome sequencing greatly increases the level of biological detail, system-level analysis of these high-dimensional datasets is becoming essential. We illustrate gene network approaches to the analysis of transcriptional data, with particular focus on the advantage of RNA-Seq technology compared to microarray platforms. We introduce a novel methodology for constructing cosplicing networks, based on distance measures combined with matrix correlations. We find that the cosplicing network is distinct and complementary to the coexpression network, although it shares the scale-free properties. In the cosplicing network, we find a set of novel hubs that have unique characteristics distinguishing them from coexpression hubs: they are heavily represented in neurobiological functional pathways and have strong overlap with markers of neurons and neuroglia, long-coding lengths, and high number of both exons and annotated transcripts. We also find that gene networks are plastic in the face of genetic and environmental pressures.

Mechano-regulation of alternative splicing

Alternative splicing contributes to the complexity of proteome by producing multiple mRNAs from a single gene. Affymetrix exon arrays and experiments in vivo or in vitro demonstrated that alternative splicing was regulated by mechanical stress. Expression of mechano-growth factor (MGF) which is the splicing isoform of insulin-like growth factor 1(IGF-1) and vascular endothelial growth factor (VEGF) splicing variants such as VEGF₁₂₁, VEGF_165, VEGF₂₀₆, VEGF₁₈₉, VEGF₁₆₅ and VEGF₁₄₅ are regulated by mechanical stress. However, the mechanism of this process is not yet clear. Increasing evidences showed that the possible mechanism is related to Ca²⁺ signal pathway and phosphorylation signal pathway. This review proposes possible mechanisms of mechanical splicing regulation. This will contribute to the biomechanical study of alternative splicing.

RNA-seq characterization of spinal cord injury transcriptome in acute/subacute phases: a resource for understanding the pathology at the systems level

Spinal cord injury (SCI) is a devastating neurological disease without effective treatment. To generate a comprehensive view of the mechanisms involved in SCI pathology, we applied RNA-Sequencing (RNA-Seq) technology to characterize the temporal changes in global gene expression after contusive SCI in mice. We sequenced tissue samples from acute and subacute phases (2 days and 7 days after injury) and systematically characterized the transcriptomes with the goal of identifying pathways and genes critical in SCI pathology. The top enriched functional categories include “inflammation response,” “neurological disease,” “cell death and survival” and “nervous system development.” The top enriched pathways include LXR/RXR Activation and Atherosclerosis Signaling, etc. Furthermore, we developed a systems-based analysis framework in order to identify key determinants in the global gene networks of the acute and sub-acute phases. Some candidate genes that we identified have been shown to play important roles in SCI, which demonstrates the validity of our approach. There are also many genes whose functions in SCI have not been well studied and can be further investigated by future experiments. We have also incorporated pharmacogenomic information into our analyses. Among the genes identified, the ones with existing drug information can be readily tested in SCI animal models. Therefore, in this study we have described an example of how global gene profiling can be translated to identifying genes of interest for functional tests in the future and generating new hypotheses. Additionally, the RNA-Seq enables splicing isoform identification and the estimation of expression levels, thus providing useful information for increasing the specificity of drug design and reducing potential side effect. In summary, these results provide a valuable reference data resource for a better understanding of the SCI process in the acute and sub-acute phases.

Alternative splicing for diseases, cancers, drugs, and databases

Alternative splicing is a major diversification mechanism in the human transcriptome and proteome. Several diseases, including cancers, have been associated with dysregulation of alternative splicing. Thus, correcting alternative splicing may restore normal cell physiology in patients with these diseases. This paper summarizes several alternative splicing-related diseases, including cancers and their target genes. Since new cancer drugs often target spliceosomes, several clinical drugs and natural products or their synthesized derivatives were analyzed to determine their effects on alternative splicing. Other agents known to have modulating effects on alternative splicing during therapeutic treatment of cancer are also discussed. Several commonly used bioinformatics resources are also summarized.

RNA editing and drug discovery for cancer therapy

RNA editing is vital to provide the RNA and protein complexity to regulate the gene expression. Correct RNA editing maintains the cell function and organism development. Imbalance of the RNA editing machinery may lead to diseases and cancers. Recently, RNA editing has been recognized as a target for drug discovery although few studies targeting RNA editing for disease and cancer therapy were reported in the field of natural products. Therefore, RNA editing may be a potential target for therapeutic natural products. In this review, we provide a literature overview of the biological functions of RNA editing on gene expression, diseases, cancers, and drugs. The bioinformatics resources of RNA editing were also summarized.

Human osteopontin splicing isoforms: known roles, potential clinical applications and activated signaling pathways

Human osteopontin is subject to alternative splicing, which generates three isoforms, termed OPNa, OPNb and OPNc. These variants show specific expression and roles in different cell contexts. We present an overview of current knowledge of the expression profile of human OPN splicing isoforms (OPN-SIs), their tissue-specific roles, and the pathways mediating their functional properties in different pathophysiological conditions. We also describe their putative application as biomarkers, and their potential use as therapeutic targets by using antibodies, oligonucleotides or siRNA molecules. This synthesis provides new clues for a better understanding of human OPN splice variants, their roles in normal and pathological conditions, and their possible clinical applications.

Dynorphin--still an extraordinarily potent opioid peptide

This issue of Molecular Pharmacology is dedicated to Dr. Avram Goldstein, the journal's founding editor and one of the leaders in the development of modern pharmacology. This article focuses on his contributions to the discovery of the dynorphins and evidence that members of this family of opioid peptides are endogenous agonists for the kappa opioid receptor. In his original publication describing the purification and sequencing of dynorphin A, Avram described this peptide as "extraordinarily potent" ("dyn" from the Greek, dynamis = power and "orphin" for endogenous morphine peptide). The name originally referred to its high affinity and great potency in the bioassay that was used to follow its activity during purification, but the name has come to have a second meaning: studies of its physiologic function in brain continue to provide powerful insights to the molecular mechanisms controlling mood disorders and drug addiction. During the 30 years since its discovery, we have learned that the dynorphin peptides are released in brain during stress exposure. After they are released, they activate kappa opioid receptors distributed throughout the brain and spinal cord, where they trigger cellular responses resulting in different stress responses: analgesia, dysphoria-like behaviors, anxiety-like responses, and increased addiction behaviors in experimental animals. Avram predicted that a detailed molecular analysis of opiate drug actions would someday lead to better treatments for drug addiction, and he would be gratified to know that subsequent studies enabled by his discovery of the dynorphins resulted in insights that hold great promise for new treatments for addiction and depressive disorders.

Omics and drug response

A new generation of technologies commonly named omics permits assessment of the entirety of the components of biological systems and produces an explosion of data and a major shift in our concepts of disease. These technologies will likely shape the future of health care. One aspect of these advances is that the data generated document the uniqueness of each human being in regard to disease risk and treatment response. These developments have reemphasized the concept of personalized medicine. Here we review the impact of omics technologies on one key aspect of personalized medicine: the individual drug response. We describe how knowledge of different omics may affect treatment decisions, namely drug choice and drug dose, and how it can be used to improve clinical outcomes.