Search

Human metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is a nuclear long noncoding RNA (lncRNA) that is highly overexpressed in many cancer tissues and plays important roles in tumor progression and metastasis. The MALAT1 primary transcript contains evolutionarily conserved structural elements in its 3′-terminal region: a triple helix forming element called element for nuclear expression (ENE) and a downstream tRNA-like structure called mascRNA. Instead of being polyadenylated, mature MALAT1 is generated by recognition and processing of the mascRNA by RNase P. A genomically encoded A-rich tract at the new 3′ end of MALAT1, which is generated upon RNase P cleavage, forms a triple helical structure with the upstream ENE. Triplex formation is vital for stabilization of the mature transcript and for subsequent accumulation and oncogenic activity of MALAT1. Here, we demonstrate that efficient 3′-end maturation of MALAT1 is dependent on an interaction between the A-rich tract and the mascRNA 3′ trailer. Using mutational analyses of cell-based reporter accumulation, we show that an extended mascRNA acceptor stem and formation of a single bulged A 5′ to the RNase P cleavage site are required for efficient maturation of the nascent MALAT1 3′ end. Our results should benefit the development of therapeutic approaches to cancer through targeting MALAT1.

Lpu-miR171 participates in lily somatic embryogenesis by regulating SCARECROW-LIKE 6 transcription factors. However, the regulatory mechanism of lpu-miR171 remains unclear. In this study, we identified the primary transcript of lpu-miR171a and analyzed the factors influencing its expression in lily somatic embryogenesis. The full-length lpu-miR171a promoter and five 5′-deletion fragments were cloned, fused to GUS, and transferred into Arabidopsis thaliana to evaluate the promoter function. The results showed that the 1248 bp sequence located upstream of the MIR171a TSS might drive GUS expression. GUS activity was detected in only the cotyledons and radicles of transgenic A. thaliana. The GUS content was highest during seed germination and then gradually declined as growth continued. The lpu-miR171a promoter responded to light, abiotic stress, and stress-related exogenous hormone signals, including ABA, SA, and MeJA. Factors affecting lpu-miR171 family member expression were also evaluated at the posttranscriptional level, and endogenous target mimics regulating lpu-miR171a and lpu-miR171b expression were predicted using the transcriptome database. The qRT-PCR results indicated that the expression levels of lpu-miR171a and lpu-miR171b were affected by competing endogenous RNAs. This study provides a molecular basis for exploring the regulatory mechanism of lpu-miR171 during somatic embryogenesis in Lilium. For the first time, we identified and analyzed the regulatory factors of lpu-miR171 at the transcriptional and posttranscriptional levels during lily somatic embryogenesis. The 1248 bp sequence located upstream of the MIR171a TSS may drive GUS expression in the cotyledons and radicles of Arabidopsis thaliana. The lpu-miR171a promoter responded to light, abiotic stress and stress-related exogenous hormone signals, including ABA, SA and MeJA. The expression of miR171a and miR171b is influenced by endogenous target mimics at the posttranscriptional level.

The microprocessor complex cleaves the primary transcript of microRNA (pri-miRNA) to initiate miRNA maturation. Microprocessor is known to consist of RNase III DROSHA and dsRNA-binding DGCR8. Here, we identify Enhancer of Rudimentary Homolog (ERH) as a new component of Microprocessor. Through a crystal structure and biochemical experiments, we reveal that ERH uses its hydrophobic groove to bind to a conserved region in the N-terminus of DGCR8, in a 2:2 stoichiometry. Knock-down of ERH or deletion of the DGCR8 N-terminus results in a reduced processing of suboptimal pri-miRNAs in polycistronic miRNA clusters. ERH increases the processing of suboptimal pri-miR-451 in a manner dependent on its neighboring pri-miR-144. Thus, the ERH dimer may mediate ‘cluster assistance’ in which Microprocessor is loaded onto a poor substrate with help from a high-affinity substrate in the same cluster. Our study reveals a role of ERH in the miRNA biogenesis pathway.

miR-18a is a member of primary transcript called miR-17-92a (C13orf25 or MIR17HG) which also contains five other miRNAs: miR-17, miR-19a, miR-20a, miR-19b and miR-92a. This cluster as a whole shows specific characteristics, where miR-18a seems to be unique. In contrast to the other members, the expression of miR-18a is additionally controlled and probably functions as its own internal controller of the cluster. miR-18a regulates many genes involved in proliferation, cell cycle, apoptosis, response to different kinds of stress, autophagy and differentiation. The disturbances of miR-18a expression are observed in cancer as well as in different diseases or pathological states. The miR-17-92a cluster is commonly described as oncogenic and it is known as ‘oncomiR-1’, but this statement is a simplification because miR-18a can act both as an oncogene and a suppressor. In this review we summarize the current knowledge about miR-18a focusing on its regulation, role in cancer biology and utility as a potential biomarker.

Human immunodeficiency virus 1 (HIV-1) is a retrovirus with a ten-kilobase single-stranded RNA genome. HIV-1 must express all of its gene products from a single primary transcript, which undergoes alternative splicing to produce diverse protein products that include structural proteins and regulatory factors1,2. Despite the critical role of alternative splicing, the mechanisms that drive the choice of splice site are poorly understood. Synonymous RNA mutations that lead to severe defects in splicing and viral replication indicate the presence of unknown cis-regulatory elements3. Here we use dimethyl sulfate mutational profiling with sequencing (DMS-MaPseq) to investigate the structure of HIV-1 RNA in cells, and develop an algorithm that we name 'detection of RNA folding ensembles using expectation-maximization' (DREEM), which reveals the alternative conformations that are assumed by the same RNA sequence. Contrary to previous models that have analysed population averages4, our results reveal heterogeneous regions of RNA structure across the entire HIV-1 genome. In addition to confirming that in vitro characterized5 alternative structures for the HIV-1 Rev responsive element also exist in cells, we discover alternative conformations at critical splice sites that influence the ratio of transcript isoforms. Our simultaneous measurement of splicing and intracellular RNA structure provides evidence for the long-standing hypothesis6-8 that heterogeneity in RNA conformation regulates splice-site use and viral gene expression.

Many RNAs are processed into biologically active transcripts, the aberrant expression of which can contribute to disease phenotypes. For example, the primary microRNA-17-92 (pri-miR-17-92) cluster contains six microRNAs (miRNAs) that collectively act in several disease settings. Herein, we used sequence-based design of structure-specific ligands to target a common structure in the Dicer processing sites of three miRNAs in the cluster, miR-17, miR-18a, and miR-20a, thereby inhibiting their biogenesis. The compound was optimized to afford a dimeric molecule that binds the Dicer processing site and an adjacent bulge, affording a 100-fold increase in potency. The dimer's mode of action was then extended from simple binding to direct cleavage by conjugation to bleomycin A5 in a manner that imparts RNA-selective cleavage or to indirect cleavage by recruiting an endogenous nuclease, or a ribonuclease targeting chimera (RIBOTAC). Interestingly, the dimer-bleomycin conjugate cleaves the entire pri-miR-17-92 cluster and hence functionally inhibits all six miRNAs emanating from it. The compound selectively reduced levels of the cluster in three disease models: polycystic kidney disease, prostate cancer, and breast cancer, rescuing disease-associated phenotypes in the latter two. Further, the bleomycin conjugate exerted selective effects on the miRNome and proteome in prostate cancer cells. In contrast, the RIBOTAC only depleted levels of pre- and mature miR-17, -18a, and 20a, with no effect on the primary transcript, in accordance with the cocellular localization of RNase L, the pre-miRNA targets, and the compound. These studies demonstrate a strategy to tune RNA structure-targeting compounds to the cellular localization of the target.

Glioblastoma (GBM) is the most malignant primary brain tumor in adults. Due to its invasive nature, it cannot be thoroughly eliminated. WD repeat domain 12 (WDR12) processes the 32S precursor rRNA but cannot affect the synthesis of the 45S/47S primary transcript. In this study, we found that WDR12 is highly expressed in GBM according to the analysis results of mRNA expression by The Cancer Genome Atlas database. The high expression level of WDR12 is dramatically related to shorter overall survival and reduced disease-free survival. Next, we knocked down WDR12 and found that knockdown of WDR12 promoted the apoptosis and inhibited the proliferation by cell biology experiments. Differential expression genes in gene-chip revealed that WDR12 knockdown mainly inhibited cell cycle. Finally, we also found that WDR12 is associated with PLK1 and EZH2 in cell proliferation of GBM. Resumptively, this report showed a possible evidence that WDR12 drove malignant behavior of GBM, whose expression may present a neoteric independent prognostic biomarker in GBM.

In Saccharomyces cerevisiae, the mitoribosomal RNA of the minor subunit,15S rRNA, is transcribed as a bicistronic transcript along with tRNA(W). 5’ and 3’ sequences flanking the mature transcript must be removed by cleavage at the respective junctions before incorporating it into the mitoribosome. An in vivo dose-response triphasic system was created to elucidate the role of Ccm1p in the processing of 15S rRNA: Ccm1p supply (“On”), deprivation (“Off”), and resupply (“Back on”). After 72 h under “Off” status, the cells started to exhibit a complete mutant phenotype as assessed by their lack of growth in glycerol medium, while keeping their mitochondrial DNA integrity (ρ(+)). Full functionality of mitochondria was reacquired upon “Back on.” 15S rRNA levels and phenotype followed the Ccm1p intramitochondrial concentrations throughout the “On-Off-Back on” course. Under “Off” status, cells gradually accumulated unprocessed 5’ and 3’ junctions, which reached significant levels at 72–96 h, probably due to a saturation of the mitochondrial degradosome (mtEXO). The Ccm1p/mtEXO mutant (Δccm1/Δdss1) showed a copious accumulation of 15S rRNA primary transcript forms, which were cleaved upon Ccm1p resupply. The gene that codes for the RNA component of RNase P was conserved in wild-type and mutant strains. Our results indicate that Ccm1p is crucial in processing the 15S rRNA primary transcript and does not stabilize the already mature 15S rRNA. Consequently, failure of this function in Δccm1 cells results, as it happens to any other unprocessed primary transcripts, in total degradation of 15S rRNA by mtEXO, whose mechanism of action is discussed.

Alternative splicing (AS) of mRNA precursors allows the synthesis of multiple mRNAs from a single primary transcript, significantly expanding the information content and regulatory possibilities of higher eukaryotic genomes. During mammalian development, AS drives certain decisive changes in different physiological processes. As development progresses, the maternal-to-zygotic transition (MZT) will trigger two processes: elimination of a subset of maternal mRNA and transcription of the zygote genome begins. Recent high-throughput technological advancements have facilitated genome-wide AS, whereas its analysis in mouse oocyte transition to the zygote stage has not been reported. We present a high-resolution global analysis of AS transitions and discovered extensive AS transitions between mouse oocyte and zygote. The difference of AS patterns was further confirmed using reverse transcription-polymerase chain reaction analysis. Many genes with specific AS events in mouse oocytes are differentially expressed between oocyte and zygote, but only a few genes with specific AS events in zygote are differentially expressed between oocyte and zygote. We provide a landscape of AS events in mouse oocyte and zygote. Our results advance the understanding of AS transitions during mouse fertilization and its potential functions for MZT and further development.

The SGI1-family elements that are specifically mobilized by the IncA- and IncC-family plasmids are important vehicles of antibiotic resistance among enteric bacteria. Although SGI1 exploits many plasmid-derived conjugation and regulatory functions, the basic mobilization module of the island is unrelated to that of IncC plasmids. This module contains the oriT and encodes the mobilization proteins MpsA and MpsB, which belong to the tyrosine recombinases and not to relaxases. Here we report an additional, essential transfer factor of SGI1. This is a small RNA deriving from the 3'-end of a primary RNA that can also serve as mRNA of ORF S022. The functional domain of this sRNA named sgm-sRNA is encoded between the mpsA gene and the oriT of SGI1. Terminator-like sequence near the promoter of the primary transcript possibly has a regulatory function in controlling the amount of full-length primary RNA, which is converted to the active sgm-sRNA through consecutive maturation steps influenced by the 5'-end of the primary RNA. The mobilization module of SGI1 seems unique due to its atypical relaxase and the newly identified sgm-sRNA, which is required for the horizontal transfer of the island but appears to act differently from classical regulatory sRNAs.

Knowledge of mitochondrial transcription start sites and promoter sequences is key to understanding mechanisms of transcription initiation in plant mitochondria. Transcription start sites can be straightforwardly determined by the mapping of primary transcript 5' ends. This chapter describes a next-generation sequencing-based protocol for the mitochondrial genome-wide mapping of transcription start sites in Arabidopsis thaliana. Like other strategies aiming at the determination of primary transcript 5' ends, this protocol exploits that only primary but not processed transcripts are 5'-triphosphorylated and, based on this property, can be enzymatically selected for. However, it uses nascent transcripts, in order to (1) enhance mitochondrial coverage compared with other compartments, (2) reduce rRNA and other background, and (3) also capture the primary 5' ends of rapidly degraded or processed transcripts.

Simple Summary Alternatively spliced tissue factor (asTF) is a naturally occurring isoform of tissue factor (TF) generated via the omission of exon 5 during the processing of TF’s primary transcript. In human and mouse, asTF protein features a unique C-terminus that lacks a transmembrane domain, rendering it soluble. asTF protein is able to associate with a subset on integrins on cell surfaces, which can trigger outside-in signaling programs in a variety of cell types. In this review, we discuss recent findings on asTF’s proto-oncogenic effects, regulatory mechanisms enabling asTF’s biosynthesis, and asTF’s potential as a biomarker and therapeutic target. Abstract In human and mouse, alternative splicing of tissue factor’s primary transcript yields two mRNA species: one features all six TF exons and encodes full-length tissue factor (flTF), and the other lacks exon 5 and encodes alternatively spliced tissue factor (asTF). flTF, which is oftentimes referred to as “TF”, is an integral membrane glycoprotein due to the presence of an alpha-helical domain in its C-terminus, while asTF is soluble due to the frameshift resulting from the joining of exon 4 directly to exon 6. In this review, we focus on asTF—the more recently discovered isoform of TF that appears to significantly contribute to the pathobiology of several solid malignancies. There is currently a consensus in the field that asTF, while dispensable to normal hemostasis, can activate a subset of integrins on benign and malignant cells and promote outside-in signaling eliciting angiogenesis; cancer cell proliferation, migration, and invasion; and monocyte recruitment. We provide a general overview of the pioneering, as well as more recent, asTF research; discuss the current concepts of how asTF contributes to cancer progression; and open a conversation about the emerging utility of asTF as a biomarker and a therapeutic target.

RNA processing converts primary transcript RNA into mature RNA. Altered RNA processing drives tumor initiation and maintenance, and may generate novel therapeutic opportunities. However, the role of RNA processing factors in gastric cancer (GC) has not been clearly elucidated. This study presents a comprehensive analysis exploring the clinical, molecular, immune, and drug response features underlying the RNA processing factors in GC. This study included 1079 GC cases from The Cancer Genome Atlas (TCGA, training set), our hospital cohort, and two other external validation sets (GSE15459, GSE62254). We developed an RNA processing-related prognostic signature using Cox regression with the least absolute shrinkage and selection operator (LASSO) penalty. The prognostic value of the signature was evaluated using a multiple-method approach. The genetic variants, pathway activation, immune heterogeneity, drug response, and splicing features associated with the risk signature were explored using bioinformatics methods. Among the tested 819 RNA processing genes, we identified five distinct RNA processing patterns with specific clinical outcomes and biological features. A 10-gene RNA processing-related prognostic signature, involving ZBTB7A, METTL2B, CACTIN, TRUB2, POLDIP3, TSEN54, SUGP1, RBMS1, TGFB1, and PWP2, was further identified. The signature was a powerful and robust prognosis factor in both the training and validation datasets. Notably, it could stratify the survival of patients with GC in specific tumor-node-metastasis (TNM) classification subgroups. We constructed a composite prognostic nomogram to facilitate clinical practice by integrating this signature with other clinical variables (TNM stage, age). Patients with low-risk scores were characterized with good clinical outcomes, proliferation, and metabolism hallmarks. Conversely, poor clinical outcome, invasion, and metastasis hallmarks were enriched in the high-risk group. The RNA processing signature was also involved in tumor microenvironment reprogramming and regulating alternative splicing, causing different drug response features between the two risk groups. The low-risk subgroup was characterized by high genomic instability, high alternative splicing and might benefit from the immunotherapy. Our findings highlight the prognostic value of RNA processing factors for patients with GC and provide insights into the specific clinical and molecular features underlying the RNA processing-related signature, which may be important for patient management and targeting treatment.

SUMMARYCis-regulatory elements (CREs) play a critical role in the development, maintenance, and disease-states of all human cell types. In the human retina, CREs have been implicated in a variety of inherited retinal disorders. To characterize cell-class-specific CREs in the human retina and elucidate their potential functions in development and disease, we performed single-nucleus (sn)ATAC-seq and snRNA-seq on the developing and adult human retina and on human retinal organoids. These analyses allowed us to identify cell-class-specific CREs, enriched transcription factor binding motifs, putative target genes, and to examine how these features change over development. By comparing DNA accessibility between the human retina and retinal organoids we found that CREs in organoids are highly correlated at the single-cell level, validating the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we studied the function of a disease-associated CRE at 5q14.3 in organoids, identifying its principal target gene as the miR-9-2 primary transcript and demonstrating a dual role for this CRE in regulating neurogenesis and gene regulatory programs in mature glia. This study provides a rich resource for characterizing cell-class-specific CREs in the human retina and showcases retinal organoids as a model in which to study the function of retinal CREs that influence retinal development and disease.HIGHLIGHTSSingle-cell map of cis-regulatory elements in developing and adult human retina.Correlation of single-cell DNA accessibility between human retina and retinal organoids.Association of disease risk loci with cell-class-specific accessibility.Modeling of enhancer function at the 5q14.3 retinal disease-risk locus.

Nascent transcript measurements derived from run-on sequencing experiments are critical for the investigation of transcriptional mechanisms and regulatory networks. However, conventional gene annotations specify the boundaries of mRNAs, which significantly differ from the boundaries of primary transcripts. Moreover, transcript isoforms with distinct transcription start and end coordinates can vary between cell types. Therefore, new primary transcript annotations are needed to accurately interpret run-on data. We developed the primaryTranscriptAnnotation R package to infer the transcriptional start and termination sites of annotated genes from genomic run-on data. We then used these inferred co-ordinates to annotate transcriptional units identified de novo. Hence, this package provides the novel utility to integrate data-driven primary transcript annotations with transcriptional unit coordinates identified in an unbiased manner. Our analyses demonstrated that this new methodology increases the sensitivity for detecting differentially expressed transcripts and provides more accurate quantification of RNA polymerase pause indices, consistent with the importance of using accurate primary transcript coordinates for interpreting genomic nascent transcription data.Availabilityhttps://github.com/WarrenDavidAnderson/genomicsRpackage/tree/master/primaryTranscriptAnnotation

The RNA guanine-N7 methyltransferase (RNMT) in complex with RNMT-activating miniprotein (RAM) catalyses the formation of a N7-methylated guanosine cap structure on the 5′ end of nascent RNA polymerase II transcripts. The mRNA cap protects the primary transcript from exonucleases and recruits cap-binding complexes that mediate RNA processing, export and translation. By using microsecond standard and accelerated molecular dynamics simulations, we provide for the first time a detailed molecular mechanism of allosteric regulation of RNMT by RAM. We show that RAM selects the RNMT active site conformations that are optimal for binding of substrates (AdoMet and the cap), thus enhancing their affinity. Furthermore, our results strongly suggest the likely scenario in which the cap binding promotes the subsequent AdoMet binding, consistent with the previously suggested cooperative binding model. By employing the network community analyses, we revealed the underlying long-range allosteric networks and paths that are crucial for allosteric regulation by RAM. Our findings complement and explain previous experimental data on RNMT activity. Moreover, this study provides the most complete description of the cap and AdoMet binding poses and interactions within the enzyme’s active site. This information is critical for the drug discovery efforts that consider RNMT as a promising anti-cancer target.

© The Author(s) 2019., Alternative splicing is facilitated by accessory proteins that guide spliceosome subunits to the primary transcript. Many of these splicing factors recognize the RNA polymerase II tail, but SFPQ is a notable exception even though essential for mammalian RNA processing. This study reveals a novel role for Dido3, one of three Dido gene products, in alternative splicing. Binding of the Dido3 amino terminus to histones and to the polymerase jaw domain was previously reported, and here we show interaction between its carboxy terminus and SFPQ. We generated a mutant that eliminates Dido3 but preserves other Dido gene products, mimicking reduced Dido3 levels in myeloid neoplasms. Dido mutation suppressed SFPQ binding to RNA and increased skipping for a large group of exons. Exons bearing recognition sequences for alternative splicing factors were nonetheless included more efficiently. Reduced SFPQ recruitment may thus account for increased skipping of SFPQ-dependent exons, but could also generate a splicing factor surplus that becomes available to competing splice sites. Taken together, our data indicate that Dido3 is an adaptor that controls SFPQ utilization in RNA splicing. Distributing splicing factor recruitment over parallel pathways provides mammals with a simple mechanism to regulate exon usage while maintaining RNA splicing efficiency., Spanish Ministry of Economics, Industry, and Competitiveness [SAF2016-75456-R, AEI/FEDER, EU]; Comunidad Autónoma de Madrid [B2017/BMD-3703 MITIC]; Foundation Alfonso Martin Escudero (to C.M.A.); CSIC predoctoral fellow FPI grant [BES-2014-068580 to C.M.G.]. Funding for open access charge: Project grants (SAF2016 and B2017) as listed under Funding.

Background Gene co-expression studies can provide important insights into molecular and cellular signaling pathways. The GeneNetwork database is a unique resource for co-expression analysis using data from a variety of tissues across genetically distinct inbred mice. However, extraction of biologically meaningful co-expressed gene sets is challenging due to variability in microarray platforms, probe quality, normalization methods, and confounding biological factors. In this study, we tested whether literature derived functional cohesion could be used as an objective metric in lieu of ‘ground truth’ to evaluate the quality of probes and microarray datasets. Results We examined Sirtuin-3 (Sirt3) co-expressed gene sets extracted from either liver or brain tissues of BXD recombinant inbred mice in the GeneNetwork database. Depending on the microarray platform, there were as many as 26 probes that targeted different regions of Sirt3 primary transcript. Co-expressed gene sets (ranging from 100–1000 genes) associated with each Sirt3 probe were evaluated using the previously developed literature-derived cohesion p-value (LPv) and benchmarked against ‘gold standards’ derived from proteomic studies or Gene Ontology classifications. We found that the maximal F-measure was obtained at an average window size of 535 genes. Using set size of 500 genes, the Pearson correlations between LPv and F-measure as well as between LPv and mitochondrial gene enrichment p-values were 0.90 and 0.93, respectively. Importantly, we found that the LPv approach can distinguish high quality Sirt3 probes. Analysis of the most functionally cohesive Sirt3 co-expressed gene set revealed core metabolic pathways that were shared between hippocampus and liver as well as distinct pathways which were unique to each tissue. These results are consistent with other studies that suggest Sirt3 is a key metabolic regulator and has distinct functions in energy-producing vs. energy-demanding tissues. Conclusions Our results provide proof-of-concept that literature cohesion analysis is useful for evaluating the quality of probes and microarray datasets, particularly when experimentally derived gold standards are unavailable. Our approach would enable researchers to rapidly identify biologically meaningful co-expressed gene sets and facilitate discovery from high throughput genomic data. Electronic supplementary material The online version of this article (10.1186/s12859-019-2621-z) contains supplementary material, which is available to authorized users.

Proper expression of the PPARG gene, which encodes a key transcription factor of adipogenesis, is indispensable in the formation of mature adipocytes. The positioning of a gene within the nuclear space has been implicated in gene regulation. We here report on the significance of the PPARG gene’s nuclear positioning for its activity during in vitro adipogenesis in the pig. We used an established system of differentiation of mesenchymal stem cells derived from bone marrow and adipose tissue into adipocytes. The differentiation process was carried out for 7 days, and the cells were examined using the 3D DNA/immuno-FISH and RNA/DNA-FISH approaches. PPARG transcript level was measured using real-time PCR, and PPARγ activity was detected with colorimetric assay. Changes in the nuclear location of the PPARG gene were observed when we compared undifferentiated mesenchymal stem cells with mature adipocytes. The gene moved from the nuclear periphery to the nuclear center as its transcriptional activity increased. The RNA/DNA-FISH approach shows that differences in primary transcript production correlated with the allele’s nuclear positioning. Transcriptionally active alleles preferentially occupy the central part of the nucleus, while inactive alleles are found on the nuclear periphery. We also show that transcription of PPARG begins with one allele, but that both alleles are active in later stages of differentiation. Our results provide evidence that functionally distinct alleles of the PPARG gene are positioned in different parts of the cell nucleus. This confirms the importance of nuclear architecture to the regulation of PPARG gene transcription, and thus to the fate of the adipose cell. Electronic supplementary material The online version of this article (10.1007/s10577-019-09604-2) contains supplementary material, which is available to authorized users.

Three out of the seven ribosomal RNA operons in Escherichia coli end in dual terminator structures. Between the two terminators of each operon is a short sequence that we report here to be an sRNA gene, transcribed as part of the ribosomal RNA primary transcript by read-through of the first terminator. The sRNA genes (rrA, rrB and rrF) from the three operons (rrnA, rrnB and rrnD) are more than 98% identical, and pull-down experiments show that their transcripts interact with Hfq and CsrA. Deletion of rrA, B, F, as well as overexpression of rrB, only modestly affect known CsrA-regulated phenotypes like biofilm formation, pgaA translation and glgC translation, and the role of the sRNAs in vivo may not yet be fully understood. Since RrA, B, F are short-lived and transcribed along with the ribosomal RNA components, their concentration reflect growth-rate regulation at the ribosomal RNA promoters and they could function to fine-tune other growth-phase-dependent processes in the cell. The primary and secondary structure of these small RNAs are conserved among species belonging to different genera of Enterobacteriales.

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC—seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

Summary Maturation of canonical microRNA (miRNA) is initiated by DROSHA that cleaves the primary transcript (pri-miRNA). More than 1,800 miRNA loci are annotated in humans, but it remains largely unknown whether and at which sites pri-miRNAs are cleaved by DROSHA. Here, we performed in vitro processing on a full set of human pri-miRNAs (miRBase version 21) followed by sequencing. This comprehensive profiling enabled us to classify miRNAs on the basis of DROSHA dependence and map their cleavage sites with respective processing efficiency measures. Only 758 pri-miRNAs are confidently processed by DROSHA, while the majority may be non-canonical or false entries. Analyses of the DROSHA-dependent pri-miRNAs show key cis-elements for processing. We observe widespread alternative processing and unproductive cleavage events such as “nick” or “inverse” processing. SRSF3 is a broad-acting auxiliary factor modulating alternative processing and suppressing unproductive processing. The profiling data and methods developed in this study will allow systematic analyses of miRNA regulation.

Cis-regulatory elements (CREs) play a critical role in the development, maintenance, and disease-states of all human cell types. In the human retina, CREs have been implicated in a variety of inherited retinal disorders. To characterize cell-class-specific CREs in the human retina and elucidate their potential functions in development and disease, we performed single-nucleus (sn)ATAC-seq and snRNA-seq on the developing and adult human retina and on human retinal organoids. These analyses allowed us to identify cell-class-specific CREs, enriched transcription factor binding motifs, putative target genes, and to examine how these features change over development. By comparing DNA accessibility between the human retina and retinal organoids we found that CREs in organoids are highly correlated at the single-cell level, validating the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we studied the function of a disease-associated CRE at 5q14.3 in organoids, identifying its principal target gene as the miR-9-2 primary transcript and demonstrating a dual role for this CRE in regulating neurogenesis and gene regulatory programs in mature glia. This study provides a rich resource for characterizing cell-class-specific CREs in the human retina and showcases retinal organoids as a model in which to study the function of retinal CREs that influence retinal development and disease.

Transcriptomics is a developing field with new methods of analysis being produced which may hold advantages in price, accuracy, or information output. QuantSeq is a form of 3' sequencing produced by Lexogen which aims to obtain similar gene-expression information to RNA-seq with significantly fewer reads, and therefore at a lower cost. QuantSeq is also able to provide information on differential polyadenylation. We applied both QuantSeq at low read depth and total RNA-seq to the same two sets of mouse spinal cord RNAs, each comprised by four controls and four mutants related to the neurodegenerative disease amyotrophic lateral sclerosis. We found substantial differences in which genes were found to be significantly differentially expressed by the two methods. Some of this difference likely due to the difference in number of reads between our QuantSeq and RNA-seq data. Other sources of difference can be explained by the differences in the way the two methods handle genes with different primary transcript lengths and how likely each method is to find a gene to be differentially expressed at different levels of overall gene expression. This work highlights how different methods aiming to assess expression difference can lead to different results.

Duchenne muscular dystrophy (DMD) is a life-shortening X-linked genetic disorder characterized by progressive wasting and weakening of muscles in boys. Loss-of-function mutations in the DMD gene, which codes for dystrophin, lead to this disease. The majority of mutations in this gene result in the exclusion of one or more exons from the transcript, eventually causing the remaining exons not to fit together correctly (i.e., out-of-frame mutations). Antisense oligonucleotides, e.g., phosphorodiamidate morpholino oligomers (PMOs), can induce therapeutic exon skipping during pre-mRNA processing to restore the reading frame of the primary transcript of DMD. As a result, truncated but partially functional dystrophin is produced, potentially slowing down the disease progression. Golodirsen is a provisionally approved PMO-based drug for approx. 8% of all DMD patients amenable to exon 53 skipping. This article summarizes golodirsen's pharmacology, efficacy and safety information. It also discusses some controversies that golodirsen met after the approval.

TheDrosophila shaggygene (sgg, GSK-3) encodes multiple protein isoforms with serine/threonine kinase activity and is a key player in diverse developmental signalling pathways. Currently it is unclear whether different Sgg proteoforms are similarly involved in signalling or if different proteoforms have distinct functions. We used CRISPR/Cas9 genome engineering to tag eight different Sgg proteoform classes and determined their localization during embryonic development. We performed proteomic analysis of the two major proteoform classes and generated mutant lines for both of these for transcriptomic and phenotypic analysis. We uncovered distinct tissue-specific localization patterns for all of the tagged proteoforms we examined, most of which have not previously been characterised directly at the protein level, including one proteoform initiating with a non-standard codon. Collectively this suggests complex developmentally regulated splicing of thesggprimary transcript. Further, affinity purification followed by mass spectrometric analyses indicate a different repertoire of interacting proteins for the two major proteoform classes we examined, one with ubiquitous expression (Sgg-PB) and one with nervous system specific expression (Sgg-PA). Specific mutation of these proteoforms shows that Sgg-PB performs the well characterised maternal and zygotic segmentations functions of thesgglocus, while Sgg-PA mutants show adult lifespan and locomotor defects consistent with its nervous system localisation. Our findings provide new insights into the role of GSK-3 proteoforms and intriguing links with the GSK-3α and GSK-3β encoded by independent vertebrate genes. Our analysis suggests that different protein isoforms generated by alternative splicing perform distinct functions.

Comprehensive research in various plants shows that the metabolic pathway of anthocyanin biosynthesis is affected by environmental factors and regulated by microRNAs through post-transcriptional regulation. In seedlings of Brassica rapa Tsuda, the accumulation of anthocyanin is induced by light. However, the roles of BrmiR828 in the light-induced synthesis of anthocyanin in Brassica rapa remain to be explored. Here, a primary transcript of BrmiR828 was identified to be located on the chromosomes of the A03 sub-genome. Five candidate MYB family genes were predicted as targets of BrmiR828 in the database of Brassica rapa (BRAD, V1.1) by using psRNATarget. The transcript abundance of mature BrmiR828 was reduced in seedlings of Brassica rapa Tsuda under blue light irradiation comparing with dark treatment. However, Real-time PCR showed the transcript level of the five candidate targets, Bra004162, Bra022602, Bra001917, Bra029113, and Bra039763 was up-regulated when the seedlings exposed to blue or UV-A light. Trans-acting siRNA gene 4 (BrTAS4) was also identified to have a higher transcript level under blue and UV-A light irradiation than that in dark treatment. RNA ligase mediated 5ˊamplification of cDNA ends (RLM-5&prime, RACE) showed that BrmiR828 can splice the mRNA of Bra039763, Bra022602, and BrTAS4 on binding sites. Phylogenetic analysis of candidate BrMYBs targets along with MYBs from Arabidopsis thaliana showed that Bra039763, Bra004162, Bra001917, Bra029113, and Bra022602 are classified to the same group with AtMYB75, AtMYB114, AtMYB90, AtMYB113, and AtMYB82 which are involved in the anthocyanin biosynthetic pathway. As a result, light-induced down-regulation of BrmiR828 can target BrTAS4, BrPAP1 (Bra039763), MYB82 (Bra022602) to negatively regulate their transcript levels leading to the accumulation of MYB transcription factors that positively regulate anthocyanin biosynthesis in light-exposed seedlings of Brassica rapa.

Neurofilaments (NFs), a major cytoskeletal component of motor neurons, play a key role in the differentiation, establishment and maintenance of their morphology and mechanical strength. The de novo assembly of these neuronal intermediate filaments requires the presence of the neurofilament light subunit (NEFL), whose expression is reduced in motor neurons in amyotrophic lateral sclerosis (ALS). This study used zebrafish as a model to characterize the NEFL homologue neflb, which encodes two different isoforms via a splicing of the primary transcript (neflbE4 and neflbE3). In vivo imaging showed that neflb is crucial for proper neuronal development, and that disrupting the balance between its two isoforms specifically affects the NF assembly and motor axon growth, with resultant motor deficits. This equilibrium is also disrupted upon the partial depletion of TDP-43 (TAR DNA-binding protein 43), an RNA-binding protein encoded by the gene TARDBP that is mislocalized into cytoplasmic inclusions in ALS. The study supports the interaction of the NEFL expression and splicing with TDP-43 in a common pathway, both biologically and pathogenetically.

The Microprocessor complex cleaves the primary transcript of microRNA (pri-miRNA) to initiate miRNA maturation. Microprocessor is known to consist of RNase III DROSHA and dsRNA-binding DGCR8. Here we identify Enhancer of Rudimentary Homolog (ERH) as a new component of the Microprocessor. ERH binds to a conserved region in the N-terminus of DGCR8. Knockdown of ERH or deletion of the DGCR8 N-terminus results in a decrease of processing of primary miRNAs with suboptimal hairpin structures that reside in polycistronic miRNA clusters. ERH increases the processing of suboptimal pri-miR-451 in a manner dependent on its neighboring pri-miR-144. Thus, the ERH dimer may mediate “cluster assistance” in which the Microprocessor is loaded onto a poor substrate with help from a high-affinity substrate in the same cluster. Our study reveals a role of ERH in the miRNA pathway.

Neurofilaments (NFs), a major cytoskeletal component of motor neurons, play a key role in their differentiation, establishment and maintenance of their morphology and mechanical strength. Thede novoassembly of these neuronal intermediate filaments requires the presence of the neurofilament light subunit, NEFL, which expression is reduced in motor neurons in Amyotrophic Lateral Sclerosis (ALS). This study used zebrafish as a model to characterize the NEFL homologueneflb, which encodes two different isoforms via splicing of the primary transcript (neflbE4andneflbE3).In vivoimaging showed thatneflbis crucial for proper neuronal development, and that disrupting the balance between its two isoforms specifically affects NF assembly and motor axon growth, with resulting motor deficits. This equilibrium is also disrupted upon partial depletion of TDP-43, a RNA binding protein that is mislocalized into cytoplasmic inclusions in ALS. The study supports interaction of NEFL expression and splicing with TDP-43 in a common pathway, both biologically and pathogenetically.

MicroRNAs (miRNAs) are processed products of primary miRNAs (pri-miRNAs) and regulate the target gene expression. Though the regulatory roles of the several mature plant miRNAs have been studied in detail, the functions of other regions of the pri-miRNAs are still unrecognized. Recent studies suggest that a few pri-miRNAs may encode small peptides, miRNA-encoded peptides (miPEPs); however, the functions of these peptides have not been studied in detail. We report that the pri-miR858a of Arabidopsis thaliana encodes a small peptide, miPEP858a, which regulates the expression of pri-miR858a and associated target genes. miPEP858a-edited and miPEP858a-overexpressing lines showed altered plant development and accumulated modulated levels of flavonoids due to changes in the expression of genes associated with the phenylpropanoid pathway and auxin signalling. The exogenous treatment of the miPEP858a-edited plants with synthetic miPEP858a complemented the phenotypes and the gene function. This study suggests the importance of miPEP858a in exerting control over plant development and the phenylpropanoid pathway.

Structural changes in RNAs are an important contributor to controlling gene expression not only at the post-transcriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5’ region of the primary transcript regulates the downstream functional unit – usually an ORF – through premature termination of transcription. Such elements not only occur naturally but they are also attractive devices in synthetic biology. The possibility to design such riboswitches or RNA thermometers is thus of considerable practical interest. Since these functional RNA elements act already during transcription, it is important to model and understand the dynamics of folding and, in particular, the formation of intermediate structures concurrently with transcription. Cotranscriptional folding simulations are therefore an important step to verify the functionality of design constructs before conducting expensive and labour-intensive wet lab experiments. For RNAs, full-fledged molecular dynamics simulations are far beyond practical reach both because of the size of the molecules and the time scales of interest. Even at the simplified level of secondary structures further approximations are necessary. The BarMap approach is based on representing the secondary structure landscape for each individual transcription step by a coarse-grained representation that only retains a small set of low-energy local minima and the energy barriers between them. The folding dynamics between two transcriptional elongation steps is modeled as a Markov process on this representation. Maps between pairs of consecutive coarse-grained landscapes make it possible to follow the folding process as it changes in response to transcription elongation.In its original implementation, the BarMap software provides a general framework to investigate RNA folding dynamics on temporally changing landscapes. It is, however, difficult to use in particular for specific scenarios such as cotranscriptional folding. To overcome this limitation, we developed the user-friendly BarMap-QA pipeline described in detail in this contribution. It is illustrated here by an elaborate example that emphasizes the careful monitoring of several quality measures. Using an iterative workflow, a reliable and complete kinetics simulation of a synthetic, transcription regulating riboswitch is obtained using minimal computational resources. All programs and scripts used in this contribution are free software and available for download as a source distribution for Linux®, or as a platform-independent Docker® image including support for Apple macOS® and Microsoft Windows®.

For decades, cardiovascular disease (CVD) has been the leading cause of death throughout most developed countries. Several studies relate RNA splicing, and more recently also circular RNAs (circRNAs), to CVD. CircRNAs originate from linear transcripts and have been shown to exhibit tissue-specific expression profiles. Here, we present an in-depth analysis of sequence, structure, modification, and cardiac circRNA interactions. We used human induced pluripotent stem cell-derived cardiac myocytes (hiPSC-CMs), human healthy and diseased (ischemic cardiomyopathy, dilated cardiomyopathy) cardiac tissue, and human umbilical vein endothelial cells (HUVECs) to profile circRNAs. We identified shared circRNAs across all samples, as well as model-specific circRNA signatures. Based on these circRNAs, we identified 63 positionally conserved and expressed circRNAs in human, pig, and mouse hearts. Furthermore, we found that the sequence of circRNAs can deviate from the sequence derived from the genome sequence, an important factor in assessing potential functions. Integration of additional data yielded evidence for m6A-methylation of circRNAs, potentially linked to translation, as well as, circRNAs overlapping with potential Argonaute 2 binding sites, indicating potential association with the RISC complex. Moreover, we describe, for the first time in cardiac model systems, a sub class of circRNAs containing the start codon of their primary transcript (AUG circRNAs) and observe an enrichment for m6A-methylation for AUG circRNAs.

The main purpose of this review is to recall for investigators - and in particular students -, some of the early data and concepts in molecular genetics and biology that are rarely cited in the current literature and are thus invariably overlooked. There is a growing tendency among editors and reviewers to consider that only data produced in the last 10–20 years or so are pertinent. However this is not the case. In exact science, sound data and lucid interpretation never become obsolete, and even if forgotten, will resurface sooner or later. In the field of gene expression, covered in the present review, recent post-genomic data have indeed confirmed many of the earlier results and concepts developed in the mid-seventies, well before the start of the recombinant DNA revolution. Human brains and even the most powerful computers, have difficulty in handling and making sense of the overwhelming flow of data generated by recent high-throughput technologies. This was easier when low throughput, more integrative methods based on biochemistry and microscopy dominated biological research. Nowadays, the need for organising concepts is ever more important, otherwise the mass of available data can generate only "building ruins" - the bricks without an architect. Concepts such as pervasive transcription of genomes, large genomic domains, full domain transcripts (FDTs) up to 100 kb long, the prevalence of post-transcriptional events in regulating eukaryotic gene expression, and the 3D-genome architecture, were all developed and discussed before 1990, and are only now coming back into vogue. Thus, to review the impact of earlier concepts on later developments in the field, I will confront former and current data and ideas, including a discussion of old and new methods. Whenever useful, I shall first briefly report post-genomic developments before addressing former results and interpretations. Equally important, some of the terms often used sloppily in scientific discussions will be clearly defined. As a basis for the ensuing discussion, some of the issues and facts related to eukaryotic gene expression will first be introduced. In chapter 2 the evolution in perception of biology over the last 60 years and the impact of the recombinant DNA revolution will be considered. Then, in chapter 3 data and theory concerning the genome, gene expression and genetics will be reviewed. The experimental and theoretical definition of the gene will be discussed before considering the 3 different types of genetic information - the "Triad" - and the importance of post-transcriptional regulation of gene expression in the light of the recent finding that 90% of genomic DNA seems to be transcribed. Some previous attempts to provide a conceptual framework for these observations will be recalled, in particular the "Cascade Regulation Hypothesis" (CRH) developed in 1967–85, and the "Gene and Genon" concept proposed in 2007. A knowledge of the size of primary transcripts is of prime importance, both for experimental and theoretical reasons, since these molecules represent the primary units of the "RNA genome" on which most of the post-transcriptional regulation of gene expression occurs. In chapter 4, I will first discuss some current post-genomic topics before summarising the discovery of the high Mr-RNA transcripts, and the investigation of their processing spanning the last 50 years. Since even today, a consensus concerning the real form of primary transcripts in eukaryotic cells has not yet been reached, I will refer to the viral and specialized cellular models which helped early on to understand the mechanisms of RNA processing and differential splicing which operate in cells and tissues. As a well-studied example of expression and regulation of a specific cellular gene in relation to differentiation and pathology, I will discuss the early and recent work on expression of the globin genes in nucleated avian erythroblasts. An important concept is that the primary transcript not only embodies protein-coding information and regulation of its expression, but also the 3D-structure of the genomic DNA from which it was derived. The wealth of recent post-genomic data published in this field emphasises the importance of a fundamental principle of genome organisation and expression that has been overlooked for years even though it was already discussed in the 1970–80ties. These issues are addressed in chapter 5 which focuses on the involvement of the nuclear matrix and nuclear architecture in DNA and RNA biology. This section will make reference to the Unified Matrix Hypothesis (UMH), which was the first molecular model of the 3D organisation of DNA and RNA. The chapter on the "RNA-genome and peripheral memories" discusses experimental data on the ribonucleoprotein complexes containing pre-mRNA (pre-mRNPs) and mRNA (mRNPs) which are organised in nuclear and cytoplasmic spaces respectively. Finally, "Outlook " will enumerate currently unresolved questions in the field, and will propose some ideas that may encourage further investigation, and comprehension of available experimental data still in need of interpretation. In chapter 8, some propositions and paradigms basic to the authors own analysis are discussed. "In conclusion" the raison d'etre of this review is recalled and positioned within the overall framework of scientific endeavour.

Ribosome biogenesis is fundamental to growth and development in eukaryotes and is linked to human diseases and cancer. Arabidopsis thaliana MORPHOLOGY OF ARGONAUTE1-52 SUPPRESSED 2 (MAS2) participates in splicing and 45S ribosomal DNA (rDNA) expression. In a screen for MAS2 interactors, we identified RIBOSOMAL RNA PROCESSING 7 (RRP7), an ortholog of yeast rRNA processing protein 7 (Rrp7), which is required for 18S ribosomal RNA (rRNA) maturation. Arabidopsis rrp7 mutants exhibit a pleiotropic phenotype including slow growth, altered shoot phyllotaxy, aberrant venation in lateral organs, partial infertility, and abscisic acid hypersensitivity in seedlings. In Arabidopsis, RRP7 localizes mainly to the nucleolus, the site of the 45S rDNA transcription that produces a 45S pre-rRNA primary transcript, precursor of the 25S, 18S and 5.8S rRNAs. Lack of RRP7 function perturbs 18S rRNA maturation, causes nucleolar hypertrophy, and results in an increased 25S/18S rRNA ratio. Arabidopsis contains hundreds of 45S rDNA genes whose expression is epigenetically regulated, and deregulated, in rrp7 mutants. Double mutant analysis revealed synergistic interactions between RRP7 alleles and alleles of MAS2, NUCLEOLIN1 (NUC1), and HISTONE DEACETYLASE 6 (HDA6), which encode epigenetic regulators of 45S rDNA transcription. Our results reveal the evolutionarily conserved but divergent roles of RRP7 as a ribosome biogenesis factor.

Prp4 kinase (Prp4k) is the first spliceosome-associated kinase shown to regulate splicing in fungi and metazoans, but nothing is yet known about its functions in plants. Here, Kanno and Venhuizen et al. report..., Splicing of precursor messenger RNAs (pre-mRNAs) is an essential step in the expression of most eukaryotic genes. Both constitutive splicing and alternative splicing, which produces multiple messenger RNA (mRNA) isoforms from a single primary transcript, are modulated by reversible protein phosphorylation. Although the plant splicing machinery is known to be a target for phosphorylation, the protein kinases involved remain to be fully defined. We report here the identification of pre-mRNA processing 4 (PRP4) KINASE A (PRP4KA) in a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana (Arabidopsis). Prp4 kinase is the first spliceosome-associated kinase shown to regulate splicing in fungi and mammals but it has not yet been studied in plants. In the same screen we identified mutants defective in SAC3A, a putative mRNA export factor that is highly coexpressed with PRP4KA in Arabidopsis. Whereas the sac3a mutants appear normal, the prp4ka mutants display a pleiotropic phenotype featuring atypical rosettes, late flowering, tall final stature, reduced branching, and lowered seed set. Analysis of RNA-sequencing data from prp4ka and sac3a mutants identified widespread and partially overlapping perturbations in alternative splicing in the two mutants. Quantitative phosphoproteomic profiling of a prp4ka mutant detected phosphorylation changes in several serine/arginine-rich proteins, which regulate constitutive and alternative splicing, and other splicing-related factors. Tests of PRP4KB, the paralog of PRP4KA, indicated that the two genes are not functionally redundant. The results demonstrate the importance of PRP4KA for alternative splicing and plant phenotype, and suggest that PRP4KA may influence alternative splicing patterns by phosphorylating a subset of splicing regulators.

Messenger RNAs (mRNAs) are molecules that represent the intermediate step in converting the genetic information carried in a cell's DNA to functional proteins. Structural features of mRNAs control both how effectively their contained data are translated into functional proteins and how rapidly they are destroyed. Keywords: transcription; translation; ribonucleases; polyadenylation; decay

The morphological and the primary transcript levels of the genes for aquaporin AaAQP4 and the membrane transporter AaV-H+-ATPase in the anal papillae of A. aegypti larvae exposed to squamocin were evaluated. Squamocin caused great vacuolization in the anal papillae wall, basal lamina disruption with loss of canalicular spaces, and disorganization of the cuticular layers. There was also decrease of AaAQP4 transcript levels and inhibition of AaV-H+-ATPase in the anal papillae of larvae. Morphological together with decrease in the AaAQP4 and AaV-H+-ATPase primary transcript levels suggest that squamocin may affect osmoregulation and ion-regulation of this insect followed by death.