Circular RNAs are temporospatially regulated throughout development and ageing in the rat.

Circular RNAs (circRNAs) are covalently closed structural isoforms of linear mRNA which have been observed across a broad range of species and tissues. Here, we provide a comprehensive circRNAs expression catalogue for the rat including 8 organs of both sexes during 4 developmental stages using a public RNAseq dataset. These analyses revealed thousands of circular RNA species, many expressed in an organ-specific manner along with their host genes which were enriched with tissue-specific biological functions. A large number of circRNAs also displayed a developmental-dependent expression pattern and are accumulated during ageing. CircRNAs also displayed some sexually dimorphic expression, with gender associated differences observed in various tissues and developmental stages. These observations suggest that circRNAs are dynamically expressed in a spatial-, temporal- and gender-specific manner in mammals, and may have important biological function in differentiation, development and aging.

MiR-137: an important player in neural development and neoplastic transformation.

MicroRNAs (miRNAs) represent an important class of small regulatory RNAs that control gene expression posttranscriptionally by targeting mRNAs for degradation or translation inhibition. Early studies have revealed a complex role for miRNAs in major biological processes such as development, differentiation, growth and metabolism. MiR-137 in particular, has been of great interest due to its critical role in brain function and putative involvement in the etiology of both neuropsychiatric disorders and cancer. Several lines of evidence suggest that development, differentiation and maturation of the nervous system is strongly linked to the expression of miR-137 and its regulation of a large number of downstream target genes in various pathways. Dysregulation of this molecule has also been implicated in major mental illnesses through its position in a variant allele highly associated with schizophrenia in the largest mega genome-wide association studies. Interestingly, miR-137 has also been shown to act as a tumor suppressor, with numerous studies finding reduced expression in neoplasia including brain tumor. Restoration of miR-137 expression has also been shown to inhibit cell proliferation, migration and metastasis, and induce cell cycle arrest, differentiation and apoptosis. These properties of miR-137 propose its potential for prognosis, diagnosis and as a therapeutic target for treatment of several human neurological and neoplastic disorders. In this review, we provide details on the discovery, targets, function, regulation and disease involvement of miR-137 with a broad look at recent discovery in this area.

Developmental suppression of schizophrenia-associated miR-137 alters sensorimotor function in zebrafish.

The neurodevelopmentally regulated microRNA miR-137 was strongly implicated as risk locus for schizophrenia in the most recent genome wide association study coordinated by the Psychiatric Genome Consortium (PGC). This molecule is highly conserved in vertebrates enabling the investigation of its function in the developing zebrafish. We utilized this model system to achieve overexpression and suppression of miR-137, both transiently and stably through transgenesis. While miR-137 overexpression was not associated with an observable specific phenotype, downregulation by antisense morpholino and/or transgenic expression of miR-sponge RNA induced significant impairment of both embryonic and larval touch-sensitivity without compromising overall anatomical development. We observed miR-137 expression and activity in sensory neurons including Rohon-Beard neurons and dorsal root ganglia, two neuronal cell types that confer touch-sensitivity in normal zebrafish, suggesting a role of these cell types in the observed phenotype. The lack of obvious anatomical or histological pathology in these cells, however, suggested that subtle axonal network defects or a change in synaptic function and neural connectivity might be responsible for the behavioral phenotype rather than a change in the cellular morphology or neuroanatomy.

Distinct miRNA expression in dorsal striatal subregions is associated with risk for addiction in rats.

Recently, we published data using an animal model that allowed us to characterize animals into two groups, addiction vulnerable and addiction resilient, where we identified that addiction/relapse vulnerability was associated with deficits in synaptic plasticity-associated gene expression in the dorsal striatum (DS). Notable was the strong reduction in expression for activity-regulated cytoskeleton-associated protein (Arc) considered a master regulator of synaptic plasticity. In the present study, we confirmed that Arc messenger RNA was significantly decreased in the DS, but importantly, we identified that this reduction was restricted to the dorsomedial (DMS) and not dorsolateral striatum (DLS). There is recent evidence of microRNA (miRNA)-associated posttranscriptional suppression of Arc and animal models of addiction have identified a key role for miRNA in the regulation of addiction-relevant genes. In further support of this link, we identified several differentially expressed miRNA with the potential to influence addiction-relevant plasticity genes, including Arc. A key study recently reported that miR-212 expression is protective against compulsive cocaine-seeking. Supporting this hypothesis, we found that miR-212 expression was significantly reduced in the DMS but not DLS of addiction-vulnerable animals. Together, our data provide strong evidence that miRNA promote ongoing plasticity deficits in the DS of addiction-vulnerable animals.

Alteration of imprinted Dlk1-Dio3 miRNA cluster expression in the entorhinal cortex induced by maternal immune activation and adolescent cannabinoid exposure.

A significant feature of the cortical neuropathology of schizophrenia is a disturbance in the biogenesis of short non-coding microRNA (miRNA) that regulate translation and stability of mRNA. While the biological origin of this phenomenon has not been defined, it is plausible that it relates to major environmental risk factors associated with the disorder such as exposure to maternal immune activation (MIA) and adolescent cannabis use. To explore this hypothesis, we administered the viral mimic poly I:C to pregnant rats and further exposed some of their maturing offsprings to daily injections of the synthetic cannabinoid HU210 for 14 days starting on postnatal day 35. Whole-genome miRNA expression analysis was then performed on the left and right hemispheres of the entorhinal cortex (EC), a region strongly associated with schizophrenia. Animals exposed to either treatment alone or in combination exhibited significant differences in the expression of miRNA in the left hemisphere, whereas the right hemisphere was less responsive. Hemisphere-associated differences in miRNA expression were greatest in the combined treatment and highly over-represented in a single imprinted locus on chromosome 6q32. This observation was significant as the syntenic 14q32 locus in humans encodes a large proportion of miRNAs differentially expressed in peripheral blood lymphocytes from patients with schizophrenia, suggesting that interaction of early and late environmental insults may affect miRNA expression, in a manner that is relevant to schizophrenia.

Genome-wide mRNA and miRNA analysis of peripheral blood mononuclear cells (PBMC) reveals different miRNAs regulating HIV/HCV co-infection.

Co-infection with human immunodeficiency virus (HIV) and hepatitis C virus (HCV) is common due to shared transmission routes. The genomic basis of HIV/HCV co-infection and its regulation by microRNA (miRNA) is unknown. Therefore, our objective was to investigate genome-wide mRNA expression and its regulation by miRNA in primary PBMCs derived from 27 patients (5 HCV - mono-infected, 5 HIV-mono-infected, 12 HCV/HIV co-infected, and 5 healthy controls). This revealed 27 miRNAs and 476 mRNAs as differentially expressed (DE) in HCV/HIV co-infection when compared to controls (adj p<0.05). Our study shows the first evidence of miRNAs specific for co-infection, several of which are correlated with key gene targets demonstrating functional relationships to pathways in cancer, immune-function, and metabolism. Notable was the up regulation of HCV-specific miR-122 in co-infection (FC>50, p=4.02E-06), which may have clinical/biological implications.

The long non-coding RNA Gomafu is acutely regulated in response to neuronal activation and involved in schizophrenia-associated alternative splicing.

Schizophrenia (SZ) is a complex disease characterized by impaired neuronal functioning. Although defective alternative splicing has been linked to SZ, the molecular mechanisms responsible are unknown. Additionally, there is limited understanding of the early transcriptomic responses to neuronal activation. Here, we profile these transcriptomic responses and show that long non-coding RNAs (lncRNAs) are dynamically regulated by neuronal activation, including acute downregulation of the lncRNA Gomafu, previously implicated in brain and retinal development. Moreover, we demonstrate that Gomafu binds directly to the splicing factors QKI and SRSF1 (serine/arginine-rich splicing factor 1) and dysregulation of Gomafu leads to alternative splicing patterns that resemble those observed in SZ for the archetypal SZ-associated genes DISC1 and ERBB4. Finally, we show that Gomafu is downregulated in post-mortem cortical gray matter from the superior temporal gyrus in SZ. These results functionally link activity-regulated lncRNAs and alternative splicing in neuronal function and suggest that their dysregulation may contribute to neurological disorders.

Decreased cortical muscarinic M1 receptors in schizophrenia are associated with changes in gene promoter methylation, mRNA and gene targeting microRNA.

Many studies have shown decreased cortical muscarinic M1 receptors (CHRM1) in schizophrenia (Sz), with one study showing Sz can be separated into two populations based on a marked loss of CHRM1 (-75%) in -25% of people (Def-Sz) with the disorder. To better understand the mechanism contributing to the loss of CHRM1 in Def-Sz, we measured specific markers of gene expression in the cortex of people with Sz as a whole, people differentiated into Def-Sz and people with Sz that do not have a deficit in cortical CHRM1 (Non-Def-Sz) and health controls. We now report that cortical CHRM1 gene promoter methylation and CHRM1 mRNA are decrease in Sz, Def-Sz and Non-Def-Sz but levels of the micro RNA (miR)-107, a CHRM1 targeting miR, are increased only in Def-Sz. We also report in vitro data strongly supporting the notion that miR-107 levels regulate CHRM1 expression. These data suggest there is a reversal of the expected inverse relationship between gene promoter methylation and CHRM1 mRNA in people with Sz and that a breakdown in gene promoter methylation control of CHRM1 expression is contributing to the global pathophysiology of the syndrome. In addition, our data argues that increased levels of at least one miR, miR-107, is contributing to the marked loss of cortical CHRM1 in Def-Sz and this may be a differentiating pathophysiology. These latter data continue to support the hypothesis that microRNAs (miRNA) have a role in the underlying neurobiology of Sz but argue they are differentially affected in subsets of people within that syndrome.

Genome-wide supported variant MIR137 and severe negative symptoms predict membership of an impaired cognitive subtype of schizophrenia.

Progress in determining the aetiology of schizophrenia (Sz) has arguably been limited by a poorly defined phenotype. We sought to delineate empirically derived cognitive subtypes of Sz to investigate the association of a genetic variant identified in a recent genome-wide association study with specific phenotypic characteristics of Sz. We applied Grade of Membership (GoM) analyses to 617 patients meeting ICD-10 criteria for Sz (n=526) or schizoaffective disorder (n=91), using cognitive performance indicators collected within the Australian Schizophrenia Research Bank. Cognitive variables included subscales from the Repeatable Battery for the Assessment of Neuropsychological Status, the Controlled Oral Word Association Test and the Letter Number Sequencing Test, and standardised estimates of premorbid and current intelligence quotient. The most parsimonious GoM solution yielded two subtypes of clinical cases reflecting those with cognitive deficits (CDs; N=294), comprising 47.6% of the sample who were impaired across all cognitive measures, and a cognitively spared group (CS; N=323) made up of the remaining 52.4% who performed relatively well on all cognitive tests. The CD subgroup were more likely to be unemployed, had an earlier illness onset, and greater severity of functional disability and negative symptoms than the CS group. Risk alleles on the MIR137 single-nucleotide polymorphism (SNP) predicted membership of CD subtype only in combination with higher severity of negative symptoms. These findings provide the first evidence for association of the MIR137 SNP with a specific Sz phenotype characterised by severe CDs and negative symptoms, consistent with the emerging role of microRNAs in the regulation of proteins responsible for neural development and function.

Imprinted DLK1-DIO3 region of 14q32 defines a schizophrenia-associated miRNA signature in peripheral blood mononuclear cells.

MicroRNAs (miRNAs) regulate gene expression at the post-transcriptional level and are important for coordinating nervous system development and neuronal function in the mature brain. We have recently identified schizophrenia-associated alteration of cortical miRNA biogenesis and expression in post-mortem brain tissue with implications for the dysregulation of schizophrenia candidate genes. Although these changes were observed in the central nervous system, it is plausible that schizophrenia-associated miRNA expression signatures may also be detected in non-neural tissue. To explore this possibility, we investigated the miRNA expression profile of peripheral blood mononuclear cells (PBMCs) from 112 patients with schizophrenia and 76 non-psychiatric controls. miRNA expression analysis of total RNA conducted using commercial miRNA arrays revealed that 33 miRNAs were significantly downregulated after correction for multiple testing with a false discovery rate (FDR) of 0%, which increased to 83 when we considered miRNA with an FDR<5%. Seven miRNAs altered in microarray analysis of schizophrenia were also confirmed to be downregulated by quantitative real-time reverse transcription-polymerase chain reaction. A large subgroup consisting of 17 downregulated miRNAs is transcribed from a single imprinted locus at the maternally expressed DLK1-DIO3 region on chromosome 14q32. This pattern of differentially expressed miRNA in PBMCs may be indicative of significant underlying genetic or epigenetic alteration associated with schizophrenia.

Schizophrenia is associated with an increase in cortical microRNA biogenesis.

MicroRNA expression profiling and quantitative reverse transcription-PCR analysis of the superior temporal gyrus and the dorsolateral prefrontal cortex revealed a significant schizophrenia-associated increase in global microRNA expression. This change was associated with an elevation of primary microRNA processing and corresponded with an increase in the microprocessor component DGCR8. The biological implications for this extensive increase in gene silencing are profound, and were exemplified by members of the miR-15 family and other related microRNA, which were significantly upregulated in both brain regions. This functionally convergent influence is overrepresented in pathways involved in synaptic plasticity and includes many genes and pathways associated with schizophrenia, some of which were substantiated in vitro by reporter gene assay. Given the magnitude of microRNA changes and their wide sphere of influence, this phenomenon could represent an important dimension in the pathogenesis of schizophrenia.

Physiological and pathophysiological influences on thirst.

Thirst motivates animals to seek fluid and drink it. It is regulated by the central nervous system and arises from neural and chemical signals from the periphery interacting in the brain to stimulate a drive to drink. Our research has focussed on the lamina terminalis and the manner in which osmotic and hormonal stimuli from the circulation are detected by neurons in this region and how that information is integrated with other neural signals to generate thirst. Our studies of osmoregulatory drinking in the sheep and rat have produced evidence that osmoreceptors for thirst exist in the dorsal cap of the organum vasculosum of the lamina terminalis (OVLT) and in the periphery of the subfornical organ, and possibly also in the median preoptic nucleus. In the rat, the hormones angiotensin II and relaxin act on neurons in the periphery of the subfornical organ to stimulate drinking. Studies of human thirst using functional magnetic resonance imaging (fMRI) techniques show that systemic hypertonicity activates the lamina terminalis and the anterior cingulate cortex, but the neural circuitry that connects sensors in the lamina terminalis to cortical regions subserving thirst remains to be determined. Regarding pathophysiological influences on thirst mechanisms, both excessive (polydipsia) and inadequate (hypodisia) water intake may have dire consequences. One of the most common primary polydipsias is that observed in some cases of schizophrenia. The neural mechanisms causing the excessive water intake in this disorder are unknown, so too are the factors that result in impaired thirst and inadequate fluid intake in some elderly humans.

Catalytic DNA: a novel tool for gene suppression.

RNA, as an intermediate in the production of every gene encoded protein and the genetic material of many pathogenic viruses, presents an attractive target for both biological and therapeutic manipulation. Despite its extensive involvement in living systems, its chemical diversity based on four units is relatively low compared with protein. This provides the opportunity for a generic approach to targeting with specificity based on primary structure rather than complex higher order structures. This form of recognition occurs naturally in complementary nucleic acids, due to an ability to bind their single stranded target through Watson-Crick interactions. The most established nucleic acid based approach to gene suppression at the RNA level is through antisense oligodeoxynucleotides (ODNs). These compounds form heteroduplex with target RNA which are thought to either block its function or mediate its destruction by activation of RNase H. Alternatively, RNA can be targeted by catalytic RNA such as the hammerhead ribozyme. Ribozymes have the advantage of being equipped with their own RNA cleavage apparatus and are therefore independent of host nuclear protein activity. At present, the utility of ribozyme oligonucleotides is restricted by the relative difficulty synthesising active molecules with sufficient resistance to nuclease degradation. Recently the power of in vitro selection has been used to evolve catalytic DNA sequences with RNA cleavage specificity and activity rivalling the very best ribozymes, while maintaining the more robust chemistry of an ODN. These deoxyribozymes or DNAzymes have tremendous potential as gene suppression agents for both target validation and therapeutic applications. A number of studies evaluating the biological activity of these compounds have shown promising results. However, as with other oligonucleotide based strategies, future exploitation of this approach may depend on accessory technology to assist with the accessibility of a target which is folded by its own secondary structure and hidden within the intracellular compartment.

The influence of arm length asymmetry and base substitution on the activity of the 10-23 DNA enzyme.

A small oligodeoxyribonucleotide derived from in vitro selection has been shown to be capable of efficient sequence-specific cleavage of RNA at purine-pyrimidine junctions. As the reaction readily takes place under simulated physiologic conditions, this molecule described as the 10-23 general purpose RNA-cleaving DNA enzyme, has potential as a therapeutic agent. To further explore the character of this prototype, we examined the influence of base substitution and binding arm length asymmetry on its RNA cleaving activity. Surprisingly, substitution of the proximal nucleotide on the 3'-arm, to allow nonstandard Watson-Crick interactions, was found in some instances to improve the cleavage reaction rate. Although the identity of the unpaired purine in the RNA substrate cleavage site was found to have only a subtle influence on the rate of catalysis, with a slight decrease observed when a G at this position was changed to an A, nucleotide substitution (G to C) in the core motif at position 14 was found to completely abolish catalysis. The effect of arm length reduction varied with RNA substrate sequence and extent of helix asymmetry. Where the cleavage rate of one substrate was impaired by truncation of the deoxyribozymes 5'-arm (6 bp), the same modification in reactions with a different sequence produced a rate enhancement. Truncation of the 3'-arm, however, had no effect on the reaction rate of the one substrate tested yet nearly halved the cleavage rate in another substrate.

Catalytic nucleic acids: from lab to applications.

Since the discovery of self-cleavage and ligation activity of the group I intron, the expansion of research interest in catalytic nucleic acids has provided a valuable nonprotein resource for manipulating biomolecules. Although a multitude of reactions can be enhanced by this class of catalyst, including trans-splicing activity of the group I intron (which could be applied to gene correction), RNA-cleaving RNA enzymes or "ribozymes" hold center stage because of their tremendous potential for mediating gene inactivation. This application has been driven predominantly by the "hammerhead" and "hairpin" ribozymes as they induce specific RNA cleavage from a very small catalytic domain, allowing delivery either as a transgene expression product or directly as a synthetic oligonucleotide. Although advances in the development of RNA modifications have improved the biological half-life of synthetic ribozymes, their use is restricted by the mechanistic dependence on conserved 2'OH-moieties. Recently a new class of catalytic nucleic acid made entirely of DNA has emerged through in vitro selection. DNA enzymes or deoxyribozyme with extraordinary RNA cleavage activity has already demonstrated their capacity for gene suppression both in vitro and in vivo. These new molecules, although rivaling the activity and stability of synthetic ribozymes, are limited equally by inefficient delivery to the intracellular target RNA. The challenge of in vivo delivery is being addressed with the assessment of a variety of approaches in animal models with the aim of bringing these compounds closer to the clinic.

Nucleic acid mutation analysis using catalytic DNA.

The sequence specificity of the '10-23' RNA-cleaving DNA enzyme (deoxyribozyme) was utilised to discriminate between subtle differences in nucleic acid sequence in a relatively conserved segment of the L1 gene from a number of different human papilloma virus (HPV) genotypes. DNA enzymes specific for the different HPV types were found to cleave their respective target oligoribonucleotide substrates with high efficiency compared with their unmatched counterparts, which were usually not cleaved or cleaved with very low efficiency. This specificity was achieved despite the existence of only very small differences in the sequence of one binding arm. As an example of how this methodology may be applied to mutation analysis of tissue samples, type-specific deoxyribozyme cleavable substrates were generated by genomic PCR using a chimeric primer containing three bases of RNA. The RNA component enabled each amplicon to be cleavable in the presence of its matching deoxyribozyme. In this format, the specificity of deoxyribozyme cleavage is defined by Watson-Crick interactions between one substrate-binding domain (arm I) and the polymorphic sequence which is amplified during PCR. Deoxy-ribozyme-mediated cleavage of amplicons generated by this method was used to examine the HPV status of genomic DNA derived from Caski cells, which are known to be positive for HPV16. This method is applicable to many types of nucleic acid sequence variation, including single nucleotide polymorphisms.

Suppression of smooth muscle cell proliferation by a c-myc RNA-cleaving deoxyribozyme.

A small catalytic DNA molecule targeting c-myc RNA was found to be a potent inhibitor of smooth muscle cell (SMC) proliferation. The catalytic domain of this molecule was based on that previously derived by in vitro selection (Santoro, S. W., and Joyce, G. F. (1997) Proc. Natl. Acad. Sci. U. S. A. 94, 4262-4266) and is known as the "10-23" general purpose RNA-cleaving deoxyribozyme. In addition to inhibiting SMC proliferation at low concentration, this molecule (targeting the translation initiation region of c-myc RNA) was found to efficiently cleave its full-length substrate in vitro and down-regulate c-myc gene expression in smooth muscle cells. The serum nuclease stability of this molecule was enhanced without substantial loss of kinetic efficiency by inclusion of a 3'-3'-internucleotide inversion at the 3'-terminal. The extent of SMC suppression was found to be influenced by the length of the substrate binding arms. This correlated to some extent with catalytic activity in both the short substrate under multiple turnover conditions and the full-length substrate under single turnover conditions, with the 9 + 9 base arm molecule producing the greatest activity.

Protein-DNA footprinting of the human epsilon-globin promoter in human intact cells using nitrogen mustard analogues and other DNA-damaging agents.

Nitrogen mustard analogues, bleomycin and dimethyl sulphate (DMS) have been used as probes of protein-DNA interactions in intact human cells. The sites of damage have been determined at base pair resolution in the single copy epsilon-globin gene promoter in erythroid K562 cells, non-erythroid HeLa cells and purified DNA. Exponential amplification of gene-specific damage fragments was achieved using the ligation-mediated polymerase chain reaction (LMPCR) technique and analysed on DNA sequencing gels. A comparison of the relative damage band intensities between purified DNA and intact cells revealed several significant differences - both protection (footprint) and enhancement. These differences occurred at putative transcription factor binding sites and hence are thought to be due to protein-DNA interactions. A major feature of the band intensity ratio plots was the footprint observed at the CCAAT box binding motif as revealed by nitrogen mustard analogues. Enhanced band intensity (hypersensitivity) was displayed at the 5'- and 3'-ends of the CCAAT box in K562 cells - this feature was absent in HeLa cells and in vitro reconstitutions. A footprint was found at the GATA-1 motif in K562 cells that was also absent in non-expressing HeLa cells. Footprints were also evident at the TATA box, CACC box and the epsilonF1 DNA binding motif in K562 cells.

Target site selection for an RNA-cleaving catalytic DNA.

A small catalytic DNA, known as the 10-23 DNA enzyme or deoxyribozyme, has been shown to efficiently hydrolyze RNA at purine-pyrimidine (R-Y) junctions in vitro. Although these potentially cleavable junctions are ubiquitous, they are often protected from deoxyribozyme activity by RNA secondary structure. We have developed a multiplex cleavage assay for screening the entire length of a target RNA molecule for deoxyribozyme cleavage sites that are efficient, both in terms of kinetics and accessibility. This strategy allowed us to simultaneously compare the RNA cleaving activity of 80 deoxyribozymes for a model target gene (HPV16 E6), and an additional 60 deoxyribozymes against the rat c-myc target. The human papilloma virus (HPV) target was used primarily to characterize the multiplex system and determine its validity. The c-myc target, coupled with a smooth muscle cell proliferation assay, allowed us to assess the relationship between in vitro cleavage efficiency and c-myc gene suppression in cell culture. The multiplex reaction approach streamlines the process of revealing effective deoxyribozymes in a functional assay and provides accessibility data that may also be applicable to site selection for other hybridization-based agents.

The DNA sequence specificity of hedamycin damage determined by ligation-mediated PCR and linear amplification.

The DNA sequence specificity of hedamycin damage was compared using two techniques: Ligation-mediated PCR (LMPCR) and Linear Amplification (LA). The electrophoretic mobility of LA products were made comparable with LMPCR products by using a primer with a 27 bp non-binding 5'-extension. By direct comparison of the damage intensity (as represented by each of the methods), considerable bias was found in the LMPCR system with respect to LA--resulting in the under representation of lesions with T as the base 3' to the damage site. In view of this observation some caution should be exercised in the interpretation of data for DNA damage/repair specificity derived by LMPCR. In addition the extended primer LA method used in these experiments, could be applied to generate a dideoxy sequencing ladder for analysis of LMPCR products. This would negate the need to prepare Maxam-Gilbert chemical sequencing fragments for amplification through LMPCR.