Scan for Motifs: a webserver for the analysis of post-transcriptional regulatory elements in the 3' untranslated regions (3' UTRs) of mRNAs

A second hotspot for pathogenic exon-skipping variants in CDC45

Biallelic pathogenic variants in CDC45 are associated with Meier-Gorlin syndrome with craniosynostosis (MGORS type 7), which also includes short stature and absent/hypoplastic patellae. Identified variants act through a hypomorphic loss of function mechanism, to reduce CDC45 activity and impact DNA replication initiation. In addition to missense and premature termination variants, several pathogenic synonymous variants have been identified, most of which cause increased exon skipping of exon 4, which encodes an essential part of the RecJ-orthologue's DHH domain. Here we have identified a second cohort of families segregating CDC45 variants, where patients have craniosynostosis and a reduction in height, alongside common facial dysmorphisms, including thin eyebrows, consistent with MGORS7. Skipping of exon 15 is a consequence of two different variants, including a shared synonymous variant that is enriched in individuals of East Asian ancestry, while other variants in trans are predicted to alter key intramolecular interactions in α/β domain II, or cause retention of an intron within the 3'UTR. Our cohort and functional data confirm exon skipping is a relatively common pathogenic mechanism in CDC45, and highlights the need for alternative splicing events, such as exon skipping, to be especially considered for variants initially predicted to be less likely to cause the phenotype, particularly synonymous variants.

RNA localization during early development of the axolotl

The asymmetric localization of biomolecules is critical for body plan development. One of the most popular model organisms for early embryogenesis studies is Xenopus laevis but there is a lack of information in other animal species. Here, we compared the early development of two amphibian species-the frog X. laevis and the axolotl Ambystoma mexicanum. This study aimed to identify asymmetrically localized RNAs along the animal-vegetal axis during the early development of A. mexicanum. For that purpose, we performed spatial transcriptome-wide analysis at low resolution, which revealed dynamic changes along the animal-vegetal axis classified into the following categories: profile alteration, de novo synthesis and degradation. Surprisingly, our results showed that many of the vegetally localized genes, which are important for germ cell development, are degraded during early development. Furthermore, we assessed the motif presence in UTRs of degraded mRNAs and revealed the enrichment of several motifs in RNAs of germ cell markers. Our results suggest novel reorganization of the transcriptome during embryogenesis of A. mexicanum to converge to the similar developmental pattern as the X. laevis.

Evolutionary conservation of maternal RNA localization in fishes and amphibians revealed by TOMO-Seq

Asymmetrical localization of biomolecules inside the egg, results in uneven cell division and establishment of many biological processes, cell types and the body plan. However, our knowledge about evolutionary conservation of localized transcripts is still limited to a few models. Our goal was to compare localization profiles along the animal-vegetal axis of mature eggs from four vertebrate models, two amphibians (Xenopus laevis, Ambystoma mexicanum) and two fishes (Acipenser ruthenus, Danio rerio) using the spatial expression method called TOMO-Seq. We revealed that RNAs of many known important transcripts such as germ layer determinants, germ plasm factors and members of key signalling pathways, are localized in completely different profiles among the models. It was also observed that there was a poor correlation between the vegetally localized transcripts but a relatively good correlation between the animally localized transcripts. These findings indicate that the regulation of embryonic development within the animal kingdom is highly diverse and cannot be deduced based on a single model.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the NAM pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

Nucleotide Variation in Drosophila cryptochrome Is Linked to Circadian Clock Function: An Association Analysis

Cryptochrome (CRY) is a conserved protein associated with the circadian clock in a broad range of organisms, including plants, insects, and mammals. In Drosophila, cry is a pleiotropic gene that encodes a blue light-dedicated circadian photoreceptor, as well as an electromagnetic field sensor and a geotaxis behavior regulator. We have generated a panel of nearly-isogenic strains that originated from various wild populations and which carry different natural alleles of cry. Sequencing of these alleles revealed substantial polymorphism, the functional role of which was elusive. To link this natural molecular diversity to gene function, we relied on association mapping. Such analysis revealed two major haplogroups consisting of six linked nucleotides associated with circadian phase (haplotypes All1/All2). We also generated a maximum-likelihood gene-tree that uncovered an additional pair of haplogroups (B1/B2). Behavioral analysis of the different haplotypes indicated significant effect on circadian phase and period, as well on the amount of activity and sleep. The data also suggested substantial epistasis between the All and B haplogroups. Intriguingly, circadian photosensitivity, assessed by light-pulse experiments, did not differ between the genotypes. Using CRISPR-mediated transgenic flies, we verified the effect of B1/B2 polymorphism on circadian phase. The transgenic flies also exhibited substantially different levels of cry transcription. We, moreover, analyzed the geographical distribution of the B1/B2 haplotypes, focusing on a 12 bp insertion/deletion polymorphism that differentiates the two haplotypes. Analysis of cry sequences in wild populations across Europe revealed a geographical cline of B1/B2 indel frequency, which correlated with seasonal bioclimatic variables. This spatial distribution of cry polymorphism reinforces the functional importance of these haplotypes in the circadian system and local adaptation.

The eukaryotic translation initiation factor eIF4E elevates steady-state m7G capping of coding and noncoding transcripts

Methyl-7-guanosine (m7G) "capping" of coding and some noncoding RNAs is critical for their maturation and subsequent activity. Here, we discovered that eukaryotic translation initiation factor 4E (eIF4E), itself a cap-binding protein, drives the expression of the capping machinery and increased capping efficiency of ∼100 coding and noncoding RNAs. To quantify this, we developed enzymatic (cap quantification; CapQ) and quantitative cap immunoprecipitation (CapIP) methods. The CapQ method has the further advantage that it captures information about capping status independent of the type of 5' cap, i.e., it is not restricted to informing on m7G caps. These methodological advances led to unanticipated revelations: 1) Many RNA populations are inefficiently capped at steady state (∼30 to 50%), and eIF4E overexpression increased this to ∼60 to 100%, depending on the RNA; 2) eIF4E physically associates with noncoding RNAs in the nucleus; and 3) approximately half of eIF4E-capping targets identified are noncoding RNAs. eIF4E's association with noncoding RNAs strongly positions it to act beyond translation. Coding and noncoding capping targets have activities that influence survival, cell morphology, and cell-to-cell interaction. Given that RNA export and translation machineries typically utilize capped RNA substrates, capping regulation provides means to titrate the protein-coding capacity of the transcriptome and, for noncoding RNAs, to regulate their activities. We also discovered a cap sensitivity element (CapSE) which conferred eIF4E-dependent capping sensitivity. Finally, we observed elevated capping for specific RNAs in high-eIF4E leukemia specimens, supporting a role for cap dysregulation in malignancy. In all, levels of capping RNAs can be regulated by eIF4E.

Tandem RNA isolation reveals functional rearrangement of RNA-binding proteins on CDKN1B/p27Kip1 3'UTRs in cisplatin treated cells

Post-transcriptional control of gene expression is mediated via RNA-binding proteins (RBPs) that interact with mRNAs in a combinatorial fashion. While recent global RNA interactome capture experiments expanded the repertoire of cellular RBPs quiet dramatically, little is known about the assembly of RBPs on particular mRNAs; and how these associations change and control the fate of the mRNA in drug-treatment conditions. Here we introduce a novel biochemical approach, termed tobramycin-based tandem RNA isolation procedure (tobTRIP), to quantify proteins associated with the 3ʹUTRs of cyclin-dependent kinase inhibitor 1B (CDKN1B/p27^Kip1) mRNAs in vivo. P27^Kip1 plays an important role in mediating a cell’s response to cisplatin (CP), a widely used chemotherapeutic cancer drug that induces DNA damage and cell cycle arrest. We found that p27^Kip1 mRNA is stabilized upon CP treatment of HEK293 cells through elements in its 3ʹUTR. Applying tobTRIP, we further compared the associated proteins in CP and non-treated cells, and identified more than 50 interacting RBPs, many functionally related and evoking a coordinated response. Knock-downs of several of the identified RBPs in HEK293 cells confirmed their involvement in CP-induced p27 mRNA regulation; while knock-down of the KH-type splicing regulatory protein (KHSRP) further enhanced the sensitivity of MCF7 adenocarcinoma cancer cells to CP treatment. Our results highlight the benefit of specific in vivo mRNA-protein interactome capture to reveal post-transcriptional regulatory networks implicated in cellular drug response and adaptation.

Evolutionary and functional implications of 3' untranslated region length of mRNAs by comprehensive investigation among four taxonomically diverse metazoan species

Background

In eukaryotic organisms, it has been well acknowledged that 3′ untranslated regions (3′ UTRs) of mRNA are actively involved in post-transcriptional regulations of gene expression. Although both shortening and lengthening of 3′ UTRs of specific candidate genes were explicitly documented to have functional consequences, landscape of 3′ UTR lengths in relation to evolutionary dynamics and biological meanings remains to be elucidated when large-scale data become available.

Objectives

The primary objective of this study was to revealed different inter- and intra-species patterns on length distribution of 3′ UTRs in comparison with 5′ UTRs and coding regions.

Methods

In the present study, we investigated 3′ UTR lengths in a highly curated set of 57,135 mRNA sequences among four well-studied and taxonomically diverse metazoan species (fruit fly, zebrafish, mouse and human).

Results

The average length ratio of 3′–5′ UTRs considerably increased from fruit fly (twofold) to human (fivefold). Moreover, genes can be characterized by the obviously different patterns of evolutionary change on 3′ UTR lengths. By utilizing the Gene Ontology annotations, genes with differential lengths of 3′ UTRs were suggested to have the divergent functional implications. In human, we further revealed that ubiquitously transcribed genes had higher median lengths of 3′ UTRs than the genes of tissue-restricted expressions.

Conclusion

We conducted a comprehensive analysis and provided an overview regarding 3′ UTRs biology of mRNAs in animals, whereas the mechanistic explanations of 3′ UTRs length variation in correlation to regulation of gene expression still remain to be further studied.

Electronic supplementary material

The online version of this article (10.1007/s13258-019-00808-8) contains supplementary material, which is available to authorized users.

The physiological role and toxicological significance of the non-metal-selective cadmium/copper-metallothionein isoform differ between embryonic and adult helicid snails

Metal regulation is essential for terrestrial gastropods to survive. In helicid snails, two metal-selective metallothionein (MT) isoforms with different functions are expressed. A cadmium-selective isoform (CdMT) plays a major role in Cd²⁺ detoxification and stress response, whereas a copper-selective MT (CuMT) is involved in Cu homeostasis and hemocyanin synthesis. A third, non-metal-selective isoform, called Cd/CuMT, was first characterized in Cantareus aspersus. The aim of this study was to quantify the transcriptional activity of all three MT genes in unexposed and metal-exposed (Cd, Cu) embryonic Roman snails. In addition, the complete Cd/CuMT mRNA of the Roman snail (Helix pomatia) was characterized, and its expression quantified in unexposed and Cd-treated adult individuals. In embryos of Helix pomatia, the Cd/CuMT gene was induced upon Cu exposure. Its transcription levels were many times higher than that of the other two MT genes, and also exceeded by far the Cd/CuMT mRNA concentrations of adult snails. In the hepatopancreas of adult Roman snails, no Cd/CuMT could be detected at the protein level, irrespective of whether the snails had been exposed to Cd or not. This contrasts with the situation in the near relative, Cantareus aspersus. It appeared that the 3'-UTR of the Cd/CuMT mRNA differed largely between Cantareus aspersus and Helix pomatia, being larger in the latter species, with a number of putative binding sites for proteins and miRNAs known to inhibit mRNA translation. We suggest this as a possible mechanism responsible for the lack of Cd/CuMT protein expression in adult Roman snails.

Characterization of a low expression haplotype in canine glutathione S-transferase (GSTT1) and its prevalence in golden retrievers

Glutathione S-transferase-theta (GSTT1) is a carcinogen detoxification enzyme, and low activity variants are associated with lymphoma in humans. We recently found a variant in the 3' untranslated region (UTR) of canine GSTT1, *101_102insT, which was predicted to change miRNA binding and was found in 5 of 17 golden retriever (GR) dogs with lymphoma but none of 14 healthy GRs. The aim of this study was to determine whether this variant led to decreased GSTT1 expression and was a discernible risk factor for lymphoma within the GR breed. On resequencing, *101_102insT appeared to be in complete linkage disequilibrium with 3 additional 3'UTR variants, leading to the inferred haplotype *3T>C; *101_102insT; *190C>A; *203T>C. In canine livers that were heterozygous for this variant haplotype, GSTT1 protein expression was significantly lower compared to the reference haplotype (densitometry .40 vs .64, P = .022), and GSTT1 transcript levels by qPCR were also significantly lower (fold difference .52, P = .012), without evidence of substantial allelic expression imbalance. The variant haplotype led to >50% decrease in expression in vitro (.31 ± .07 vs .64 ± .19; P = .019). We found no significant difference in minor allele frequencies between 71 GR dogs with lymphoma (MAF .162) and 33 healthy age-matched controls (MAF .136, P = .69). Our results indicate that the variant GSTT1 3'UTR haplotype containing *101_102insT reduces gene expression, which could lead to impaired carcinogen detoxification, but was not a detectable risk factor for lymphoma in GR dogs.

Identification of single nucleotide polymorphisms in carnosine-related genes and effects of genotypes on pork meat quality attributes

Carnosine has pH-buffering and antioxidant properties that may bring advantages in terms of meat quality attributes. This study aimed at identifying polymorphisms in carnosine-related genes (CARNS1, SLC6A6, SLC15A3, SLC15A4) that might associate with muscle carnosine content and meat quality traits in pigs (Duroc, Landrace, Yorkshire). Twenty seven SNPs were identified and association analyses performed for SLC15A3 c.*35C>T and c.*52C>T (3' UTR region), and SLC15A4 c.658A>G (Ile220Val) and c.818G>A (Ser273Asn) SNPs. Associations were observed for SNP c.658A>G with carnosine content, color b* and L*, drip and cooking losses, pH24h and glycolytic potential values (P≤0.05). The same associations were observed for SNP c.818G>A, but they were not significant after FDR correction. Results suggest that specific SLC15A4 gene variants might increase muscle carnosine content and improve meat quality. With a minor allele frequency of 0.17 for SNP c.658A>G in Yorkshire pigs, selection in favor of the c.658A allele may be considered as a mean to improve pork quality attributes.

Introduction to Bioinformatics Resources for Post-transcriptional Regulation of Gene Expression

Untranslated regions (UTRs) and, to a lesser extent, coding sequences of mRNAs are involved in defining the fate of the mature transcripts through the modulation of three primary control processes, mRNA localization, degradation and translation; the action of trans-factors such as RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) combined with the presence of defined sequence and structural cis-elements ultimately determines translation levels. Identifying functional regions in UTRs and uncovering post-transcriptional regulators acting upon these regions is thus of paramount importance to understand the spectrum of regulatory possibilities for any given mRNA. This tasks can now be approached computationally, to reduce the space of testable hypotheses and to drive experimental validation.This chapter focuses on presenting databases and tools allowing to study the various aspects of post-transcriptional regulation, including motif search (sequence and secondary structure), prediction of regulatory networks (e.g., RBP and ncRNA binding sites), profiling of the mRNAs translational state, and other aspects of this level of gene expression regulation. Two analysis pipelines are also presented as practical examples of how the described tools could be integrated and effectively employed.

Combinatorial Control of mRNA Fates by RNA-Binding Proteins and Non-Coding RNAs

Post-transcriptional control of gene expression is mediated by RNA-binding proteins (RBPs) and small non-coding RNAs (e.g., microRNAs) that bind to distinct elements in their mRNA targets. Here, we review recent examples describing the synergistic and/or antagonistic effects mediated by RBPs and miRNAs to determine the localisation, stability and translation of mRNAs in mammalian cells. From these studies, it is becoming increasingly apparent that dynamic rearrangements of RNA-protein complexes could have profound implications in human cancer, in synaptic plasticity, and in cellular differentiation.