John Toubia

The RNA Binding Protein Quaking Regulates Formation of circRNAs

Circular RNAs (circRNAs), formed by non-sequential back-splicing of pre-mRNA transcripts, are a widespread form of non-coding RNA in animal cells. However, it is unclear whether the majority of circRNAs represent splicing by-products without function or are produced in a regulated manner to carry out specific cellular functions. We show that hundreds of circRNAs are regulated during human epithelial-mesenchymal transition (EMT) and find that the production of over one-third of abundant circRNAs is dynamically regulated by the alternative splicing factor, Quaking (QKI), which itself is regulated during EMT. Furthermore, by modulating QKI levels, we show the effect on circRNA abundance is dependent on intronic QKI binding motifs. Critically, the addition of QKI motifs is sufficient to induce de novo circRNA formation from transcripts that are normally linearly spliced. These findings demonstrate circRNAs are both purposefully synthesized and regulated by cell-type specific mechanisms, suggesting they play specific biological roles in EMT.

Dichotomy in the NRT Gene Families of Dicots and Grass Species

A large proportion of the nitrate (NO3 −) acquired by plants from soil is actively transported via members of the NRT families of NO3 − transporters. In Arabidopsis, the NRT1 family has eight functionally characterised members and predominantly comprises low-affinity transporters; the NRT2 family contains seven members which appear to be high-affinity transporters; and there are two NRT3 (NAR2) family members which are known to participate in high-affinity transport. A modified reciprocal best hit (RBH) approach was used to identify putative orthologues of the Arabidopsis NRT genes in the four fully sequenced grass genomes (maize, rice, sorghum, Brachypodium). We also included the poplar genome in our analysis to establish whether differences between Arabidopsis and the grasses may be generally applicable to monocots and dicots. Our analysis reveals fundamental differences between Arabidopsis and the grass species in the gene number and family structure of all three families of NRT transporters. All grass species possessed additional NRT1.1 orthologues and appear to lack NRT1.6/NRT1.7 orthologues. There is significant separation in the NRT2 phylogenetic tree between NRT2 genes from dicots and grass species. This indicates that determination of function of NRT2 genes in grass species will not be possible in cereals based simply on sequence homology to functionally characterised Arabidopsis NRT2 genes and that proper functional analysis will be required. Arabidopsis has a unique NRT3.2 gene which may be a fusion of the NRT3.1 and NRT3.2 genes present in all other species examined here. This work provides a framework for future analysis of NO3 − transporters and NO3 − transport in grass crop species.

HvZIP7 mediates zinc accumulation in barley (Hordeum vulgare) at moderately high zinc supply

High expression of zinc (Zn)-regulated, iron-regulated transporter-like protein (ZIP) genes increases root Zn uptake in dicots, leading to high accumulation of Zn in shoots. However, none of the ZIP genes tested previously in monocots could enhance shoot Zn accumulation. In this report, barley (Hordeum vulgare) HvZIP7 was investigated for its functions in Zn transport. The functions of HvZIP7 in planta were studied using in situ hybridization and transient analysis of subcellular localization with a green fluorescent protein (GFP) reporter. Transgenic barley lines overexpressing HvZIP7 were also generated to further understand the functions of HvZIP7 in metal transport. HvZIP7 is strongly induced by Zn deficiency, primarily in vascular tissues of roots and leaves, and its protein was localized in the plasma membrane. These properties are similar to its closely related homologs in dicots. Overexpression of HvZIP7 in barley plants increased Zn uptake when moderately high concentrations of Zn were supplied. Significantly, there was a specific enhancement of shoot Zn accumulation, with no measurable increase in iron (Fe), manganese (Mn), copper (Cu) or cadmium (Cd). HvZIP7 displays characteristics of low-affinity Zn transport. The unique function of HvZIP7 provides new insights into the role of ZIP genes in Zn homeostasis in monocots, and offers opportunities to develop Zn biofortification strategies in cereals.

Assessing the gene regulatory properties of Argonaute-bound small RNAs of diverse genomic origin.

High-throughput sequencing reveals an abundance of microRNA-sized fragments derived from larger non-coding RNAs. Roles for these small RNAs in gene silencing are suggested by their co-precipitation with Argonaute, the microRNA effector protein, though the extent to which they suppress gene expression endogenously remains unclear. To address this, we used luciferase reporters to determine the endogenous functionality of small RNAs from a diverse range of sources. We demonstrate small RNAs derived from snoRNAs have the capacity to act in a microRNA-like manner, though we note the vast majority of these are bound to Argonaute at levels below that required for detectable silencing activity. We show Argonaute exhibits a high degree of selectivity for the small RNAs with which it interacts and note that measuring Argonaute-associated levels is a better indicator of function than measuring total expression. Although binding to Argonaute at sufficient levels is necessary for demonstrating microRNA functionality in our reporter assay, this alone is not enough as some small RNAs derived from other non-coding RNAs (tRNAs, rRNAs, Y-RNAs) are associated with Argonaute at very high levels yet do not serve microRNA-like roles.

Targeting sphingosine kinase 1 induces MCL1 dependent cell death in acute myeloid leukemia

Acute myeloid leukemia (AML) is an aggressive malignancy where despite improvements in conventional chemotherapy and bone marrow transplantation, overall survival remains poor. Sphingosine kinase 1 (SPHK1) generates the bioactive lipid sphingosine 1-phosphate (S1P) and has established roles in tumor initiation, progression, and chemotherapy resistance in a wide range of cancers. The role and targeting of SPHK1 in primary AML, however, has not been previously investigated. Here we show that SPHK1 is overexpressed and constitutively activated in primary AML patient blasts but not in normal mononuclear cells. Subsequent targeting of SPHK1 induced caspase-dependent cell death in AML cell lines, primary AML patient blasts, and isolated AML patient leukemic progenitor/stem cells, with negligible effects on normal bone marrow CD34+ progenitors from healthy donors. Furthermore, administration of SPHK1 inhibitors to orthotopic AML patient-derived xenografts reduced tumor burden and prolonged overall survival without affecting murine hematopoiesis. SPHK1 inhibition was associated with reduced survival signaling from S1P receptor 2, resulting in selective downregulation of the prosurvival protein MCL1. Subsequent analysis showed that the combination of BH3 mimetics with either SPHK1 inhibition or S1P receptor 2 antagonism triggered synergistic AML cell death. These results support the notion that SPHK1 is a bona fide therapeutic target for the treatment of AML.

Clusters of genes encoding fructan biosynthesizing enzymes in wheat and barley

Fructans are soluble carbohydrates with health benefits and possible roles in plant adaptation. Fructan biosynthetic genes were isolated using comparative genomics and physical mapping followed by BAC sequencing in barley. Genes encoding sucrose:sucrose 1-fructosyltransferase (1-SST), fructan:fructan 1-fructosyltransferase (1-FFT) and sucrose:fructan 6-fructosyltransferase (6-SFT) were clustered together with multiple copies of vacuolar invertase genes and a transposable element on two barley BAC. Intron–exon structures of the genes were similar. Phylogenetic analysis of the fructosyltransferases and invertases in the Poaceae showed that the fructan biosynthetic genes may have evolved from vacuolar invertases. Quantitative real-time PCR was performed using leaf RNA extracted from three wheat cultivars grown under different conditions. The 1-SST, 1-FFT and 6-SFT genes had correlated expression patterns in our wheat experiment and in existing barley transcriptome database. Single nucleotide polymorphism (SNP) markers were developed and successfully mapped to a major QTL region affecting wheat grain fructan accumulation in two independent wheat populations. The alleles controlling high- and low- fructan in parental lines were also found to be associated in fructan production in a diverse set of 128 wheat lines. To the authors’ knowledge, this is the first report on the mapping and sequencing of a fructan biosynthetic gene cluster and in particular, the isolation of a novel 1-FFT gene from barley.

Naturally existing isoforms of miR-222 have distinct functions

Abstract Deep-sequencing reveals extensive variation in the sequence of endogenously expressed microRNAs (termed ‘isomiRs’) in human cell lines and tissues, especially in relation to the 3′ end. From the immunoprecipitation of the microRNA-binding protein Argonaute and the sequencing of associated small RNAs, we observe extensive 3′-isomiR variation, including for miR-222 where the majority of endogenously expressed miR-222 is extended by 1–5 nt compared to the canonical sequence. We demonstrate this 3′ heterogeneity has dramatic implications for the phenotype of miR-222 transfected cells, with longer isoforms promoting apoptosis in a size (but not 3′ sequence)-dependent manner. The transfection of longer miR-222 isomiRs did not induce an interferon response, but did downregulate the expression of many components of the pro-survival PI3K-AKT pathway including PIK3R3, a regulatory subunit whose knockdown phenocopied the expression of longer 222 isoforms in terms of apoptosis and the inhibition of other PI3K-AKT genes. As this work demonstrates the capacity for 3′ isomiRs to mediate differential functions, we contend more attention needs to be given to 3′ variance given the prevalence of this class of isomiR.

Maize maintains growth in response to decreased nitrate supply through a highly dynamic and developmental stage‐specific transcriptional response

Elucidation of the gene networks underlying the response to N supply and demand will facilitate the improvement of the N uptake efficiency of plants. We undertook a transcriptomic analysis of maize to identify genes responding to both a non-growth-limiting decrease in NO3- provision and to development-based N demand changes at seven representative points across the life cycle. Gene co-expression networks were derived by cluster analysis of the transcript profiles. The majority of NO3--responsive transcription occurred at 11 (D11), 18 (D18) and 29 (D29) days after emergence, with differential expression predominating in the root at D11 and D29 and in the leaf at D18. A cluster of 98 probe sets was identified, the expression pattern of which is similar to that of the high-affinity NO3- transporter (NRT2) genes across the life cycle. The cluster is enriched with genes encoding enzymes and proteins of lipid metabolism and transport, respectively. These are candidate genes for the response of maize to N supply and demand. Only a few patterns of differential gene expression were observed over the entire life cycle; however, the composition of the classes of the genes differentially regulated at individual time points was unique, suggesting tightly controlled regulation of NO3--responsive gene expression.

Don't go in circles: confounding factors in gene expression profiling

Quantification of gene expression is a crucial research tool in the life sciences, which makes it important to identify any factors that could compromise its accuracy. One of these factors are non‐polyadenylated (poly(A)−) transcripts, including circular RNAs (circRNAs) that can skew quantification of gene expression as they resemble messenger RNAs (mRNAs). Here, we highlight the impact circRNAs and other poly(A)− transcripts have on gene expression profiling and the biological conclusions drawn from such experiments. We also propose easily adoptable strategies to increase the accuracy of gene expression quantification.

miR‐200/375 control epithelial plasticity‐associated alternative splicing by repressing the RNA‐binding protein Quaking

Members of the miR‐200 family are critical gatekeepers of the epithelial state, restraining expression of pro‐mesenchymal genes that drive epithelial–mesenchymal transition (EMT) and contribute to metastatic cancer progression. Here, we show that miR‐200c and another epithelial‐enriched miRNA, miR‐375, exert widespread control of alternative splicing in cancer cells by suppressing the RNA‐binding protein Quaking (QKI). During EMT, QKI‐5 directly binds to and regulates hundreds of alternative splicing targets and exerts pleiotropic effects, such as increasing cell migration and invasion and restraining tumour growth, without appreciably affecting mRNA levels. QKI‐5 is both necessary and sufficient to direct EMT‐associated alternative splicing changes, and this splicing signature is broadly conserved across many epithelial‐derived cancer types. Importantly, several actin cytoskeleton‐associated genes are directly targeted by both QKI and miR‐200c, revealing coordinated control of alternative splicing and mRNA abundance during EMT. These findings demonstrate the existence of a miR‐200/miR‐375/QKI axis that impacts cancer‐associated epithelial cell plasticity through widespread control of alternative splicing.