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.

Nuclear export of circular RNA.

Circular RNAs (circRNAs), which are increasingly being implicated in a variety of functions in normal and cancerous cells1-5, are formed by back-splicing of precursor mRNAs in the nucleus6-10. circRNAs are predominantly localized in the cytoplasm, indicating that they must be exported from the nucleus. Here we identify a pathway that is specific for the nuclear export of circular RNA. This pathway requires Ran-GTP, exportin-2 and IGF2BP1. Enhancing the nuclear Ran-GTP gradient by depletion or chemical inhibition of the major protein exporter CRM1 selectively increases the nuclear export of circRNAs, while reducing the nuclear Ran-GTP gradient selectively blocks circRNA export. Depletion or knockout of exportin-2 specifically inhibits nuclear export of circRNA. Analysis of nuclear circRNA-binding proteins reveals that interaction between IGF2BP1 and circRNA is enhanced by Ran-GTP. The formation of circRNA export complexes in the nucleus is promoted by Ran-GTP through its interactions with exportin-2, circRNA and IGF2BP1. Our findings demonstrate that adaptors such as IGF2BP1 that bind directly to circular RNAs recruit Ran-GTP and exportin-2 to export circRNAs in a mechanism that is analogous to protein export, rather than mRNA export.

A Prox1 enhancer represses haematopoiesis in the lymphatic vasculature.

Transcriptional enhancer elements are responsible for orchestrating the temporal and spatial control over gene expression that is crucial for programming cell identity during development1-3. Here we describe a novel enhancer element that is important for regulating the expression of Prox1 in lymphatic endothelial cells. This evolutionarily conserved enhancer is bound by key lymphatic transcriptional regulators including GATA2, FOXC2, NFATC1 and PROX1. Genome editing of the enhancer to remove five nucleotides encompassing the GATA2-binding site resulted in perinatal death of homozygous mutant mice due to profound lymphatic vascular defects. Lymphatic endothelial cells in enhancer mutant mice exhibited reduced expression of genes characteristic of lymphatic endothelial cell identity and increased expression of genes characteristic of haemogenic endothelium, and acquired the capacity to generate haematopoietic cells. These data not only reveal a transcriptional enhancer element important for regulating Prox1 expression and lymphatic endothelial cell identity but also demonstrate that the lymphatic endothelium has haemogenic capacity, ordinarily repressed by Prox1.

Genome-wide Analyses of Chromatin State in Human Mast Cells Reveal Molecular Drivers and Mediators of Allergic and Inflammatory Diseases.

Mast cells (MCs) are versatile immune cells capable of rapidly responding to a diverse range of extracellular cues. Here, we mapped the genomic and transcriptomic changes in human MCs upon diverse stimuli. Our analyses revealed broad H3K4me3 domains and enhancers associated with activation. Notably, the rise of intracellular calcium concentration upon immunoglobulin E (IgE)-mediated crosslinking of the high-affinity IgE receptor (FcεRI) resulted in genome-wide reorganization of the chromatin landscape and was associated with a specific chromatin signature, which we term Ca2+-dependent open chromatin (COC) domains. Examination of differentially expressed genes revealed potential effectors of MC function, and we provide evidence for fibrinogen-like protein 2 (FGL2) as an MC mediator with potential relevance in chronic spontaneous urticaria. Disease-associated single-nucleotide polymorphisms mapped onto cis-regulatory regions of human MCs suggest that MC function may impact a broad range of pathologies. The datasets presented here constitute a resource for the further study of MC function.

ESRP1 controls biogenesis and function of a large abundant multiexon circRNA.

While the majority of circRNAs are formed from infrequent back-splicing of exons from protein coding genes, some can be produced at quite high level and in a regulated manner. We describe the regulation, biogenesis and function of circDOCK1(2-27), a large, abundant circular RNA that is highly regulated during epithelial-mesenchymal transition (EMT) and whose formation depends on the epithelial splicing regulator ESRP1. CircDOCK1(2-27) synthesis in epithelial cells represses cell motility both by diverting transcripts from DOCK1 mRNA production to circRNA formation and by direct inhibition of migration by the circRNA. HITS-CLIP analysis and CRISPR-mediated deletions indicate ESRP1 controls circDOCK1(2-27) biosynthesis by binding a GGU-containing repeat region in intron 1 and detaining its splicing until Pol II completes its 157 kb journey to exon 27. Proximity-dependent biotinylation (BioID) assay suggests ESRP1 may modify the RNP landscape of intron 1 in a way that disfavours communication of exon 1 with exon 2, rather than physically bridging exon 2 to exon 27. The X-ray crystal structure of RNA-bound ESRP1 qRRM2 domain reveals it binds to GGU motifs, with the guanines embedded in clamp-like aromatic pockets in the protein.

Dual inhibition of airway inflammation and fibrosis by common ß cytokine receptor blockade.

BACKGROUND: Patients with severe asthma can present with eosinophilic type 2 (T2), neutrophilic, or mixed inflammation that drives airway remodeling and exacerbations and represents a major treatment challenge. The common ß (ßc) receptor signals for 3 cytokines, GM-CSF, IL-5, and IL-3, which collectively mediate T2 and neutrophilic inflammation. OBJECTIVE: To determine the pathogenesis of ßc receptor-mediated inflammation and remodeling in severe asthma and to investigate ßc antagonism as a therapeutic strategy for mixed granulocytic airway disease. METHODS: ßc gene expression was analyzed in bronchial biopsy specimens from patients with mild-to-moderate and severe asthma. House dust mite extract and Aspergillus fumigatus extract (ASP) models were used to establish asthma-like pathology and airway remodeling in human ßc transgenic mice. Lung tissue gene expression was analyzed by RNA sequencing. The mAb CSL311 targeting the shared cytokine binding site of ßc was used to block ßc signaling. RESULTS: ßc gene expression was increased in patients with severe asthma. CSL311 potently reduced lung neutrophils, eosinophils, and interstitial macrophages and improved airway pathology and lung function in the acute steroid-resistant house dust mite extract model. Chronic intranasal ASP exposure induced airway inflammation and fibrosis and impaired lung function that was inhibited by CSL311. CSL311 normalized the ASP-induced fibrosis-associated extracellular matrix gene expression network and strongly reduced signatures of cellular inflammation in the lung. CONCLUSIONS: ßc cytokines drive steroid-resistant mixed myeloid cell airway inflammation and fibrosis. The anti-ßc antibody CSL311 effectively inhibits mixed T2/neutrophilic inflammation and severe asthma-like pathology and reverses fibrosis gene signatures induced by exposure to commonly encountered environmental allergens.

An orthotopic syngeneic mouse model of bortezomib-resistant multiple myeloma.

While bortezomib has significant benefits in multiple myeloma (MM) therapy, the disease remains incurable due to the invariable development of bortezomib resistance. This emphasises the need for advanced models for preclinical evaluation of new therapeutic approaches for bortezomib-resistant MM. Here, we describe the development of an orthotopic syngeneic bortezomib-resistant MM mouse model based on the most well-characterised syngeneic MM mouse model derived from spontaneous MM-forming C57BL/KaLwRij mice. Using bortezomib-resistant 5TGM1 cells, we report and characterise a robust syngeneic mouse model of bortezomib-resistant MM that is well suited to the evaluation of new therapeutic approaches for proteasome inhibitor-resistant MM.

capCLIP: a new tool to probe translational control in human cells through capture and identification of the eIF4E-mRNA interactome.

Translation of eukaryotic mRNAs begins with binding of their m7G cap to eIF4E, followed by recruitment of other translation initiation factor proteins. We describe capCLIP, a novel method to comprehensively capture and quantify the eIF4E (eukaryotic initiation factor 4E) 'cap-ome' and apply it to examine the biological consequences of eIF4E-cap binding in distinct cellular contexts. First, we use capCLIP to identify the eIF4E cap-omes in human cells with/without the mTORC1 (mechanistic target of rapamycin, complex 1) inhibitor rapamycin, there being an emerging consensus that rapamycin inhibits translation of TOP (terminal oligopyrimidine) mRNAs by displacing eIF4E from their caps. capCLIP reveals that the representation of TOP mRNAs in the cap-ome is indeed systematically reduced by rapamycin, thus validating our new methodology. capCLIP also refines the requirements for a functional TOP sequence. Second, we apply capCLIP to probe the consequences of phosphorylation of eIF4E. We show eIF4E phosphorylation reduces overall eIF4E-mRNA association and, strikingly, causes preferential dissociation of mRNAs with short 5'-UTRs. capCLIP is a valuable new tool to probe the function of eIF4E and of other cap-binding proteins such as eIF4E2/eIF4E3.

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.

Extensive transcriptional responses are co-ordinated by microRNAs as revealed by Exon-Intron Split Analysis (EISA).

Epithelial-mesenchymal transition (EMT) has been a subject of intense scrutiny as it facilitates metastasis and alters drug sensitivity. Although EMT-regulatory roles for numerous miRNAs and transcription factors are known, their functions can be difficult to disentangle, in part due to the difficulty in identifying direct miRNA targets from complex datasets and in deciding how to incorporate 'indirect' miRNA effects that may, or may not, represent biologically relevant information. To better understand how miRNAs exert effects throughout the transcriptome during EMT, we employed Exon-Intron Split Analysis (EISA), a bioinformatic technique that separates transcriptional and post-transcriptional effects through the separate analysis of RNA-Seq reads mapping to exons and introns. We find that in response to the manipulation of miRNAs, a major effect on gene expression is transcriptional. We also find extensive co-ordination of transcriptional and post-transcriptional regulatory mechanisms during both EMT and mesenchymal to epithelial transition (MET) in response to TGF-ß or miR-200c respectively. The prominent transcriptional influence of miRNAs was also observed in other datasets where miRNA levels were perturbed. This work cautions against a narrow approach that is limited to the analysis of direct targets, and demonstrates the utility of EISA to examine complex regulatory networks involving both transcriptional and post-transcriptional mechanisms.

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.

Naturally existing isoforms of miR-222 have distinct functions.

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.

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.

Nitrogen assimilation system in maize is regulated by developmental and tissue-specific mechanisms.

KEY MESSAGE: We found metabolites, enzyme activities and enzyme transcript abundances vary significantly across the maize lifecycle, but weak correlation exists between the three groups. We identified putative genes regulating nitrate assimilation. Progress in improving nitrogen (N) use efficiency (NUE) of crop plants has been hampered by the complexity of the N uptake and utilisation systems. To understand this complexity we measured the activities of seven enzymes and ten metabolites related to N metabolism in the leaf and root tissues of Gaspe Flint maize plants grown in 0.5 or 2.5 mM NO3 (-) throughout the lifecycle. The amino acids had remarkably similar profiles across the lifecycle except for transient responses, which only appeared in the leaves for aspartate or in the roots for asparagine, serine and glycine. The activities of the enzymes for N assimilation were also coordinated to a certain degree, most noticeably with a peak in root activity late in the lifecycle, but with wide variation in the activity levels over the course of development. We analysed the transcriptional data for gene sets encoding the measured enzymes and found that, unlike the enzyme activities, transcript levels of the corresponding genes did not exhibit the same coordination across the lifecycle and were only weakly correlated with the levels of various amino acids or individual enzyme activities. We identified gene sets which were correlated with the enzyme activity profiles, including seven genes located within previously known quantitative trait loci for enzyme activities and hypothesise that these genes are important for the regulation of enzyme activities. This work provides insights into the complexity of the N assimilation system throughout development and identifies candidate regulatory genes, which warrant further investigation in efforts to improve NUE in crop plants.

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.

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.

Increased prevalence of germline pathogenic CHEK2 variants in individuals with pituitary adenomas.

CONTEXT: CHEK2 is a cell cycle checkpoint regulator gene with a long-established role as a clinically relevant, moderate risk breast cancer predisposition gene, with greater risk ascribed to truncating variants than missense variants. METHODS: We assessed 165 individuals with pituitary adenomas for CHEK2 variants. The study consisted of a primary cohort of 29 individuals who underwent germline and tumour whole exome sequencing, and a second, independent cohort of 136 individuals who had a targeted next-generation sequencing panel performed on both germline and tumour DNA (n=52) or germline DNA alone (n=84). RESULTS: We identified rare, coding, non-synonymous germline CHEK2 variants amongst 3/29 (10.3%) patients in our primary cohort and 5/165 (3.0%) patients overall, with affected patients having a range of hormone secretion types (prolactinoma, thyrotrophinoma, somatotrophinoma and non-functioning pituitary adenoma). No somatic variants were identified. Two variants were definitive null variants (c.1100delC, c.444+1G>A), classified as pathogenic. Two variants were missense variants (p.Asn186His, p.Thr476Met), classified as likely pathogenic. Even when considering the null variants only, the rate of CHEK2 variants was higher in our cohort compared to national control data (1.8% vs. 0.5%, P=0.049). CONCLUSIONS: This is the first study to suggest a role for the breast cancer predisposition gene, CHEK2, in pituitary tumorigenesis, with pathogenic/likely pathogenic variants found in 3% of patients with pituitary adenomas. As pituitary adenomas are relatively common and typically lack classical autosomal dominant family histories, risk alleles - such as these variants found in CHEK2 - might be a significant contributor to pituitary adenoma risk in the general population.

PAM variants in patients with thyrotrophinomas, cyclical Cushing's disease and prolactinomas.

Introduction: Germline loss-of-function variants in PAM, encoding peptidylglycine α-amidating monooxygenase (PAM), were recently discovered to be enriched in conditions of pathological pituitary hypersecretion, specifically: somatotrophinoma, corticotrophinoma, and prolactinoma. PAM is the sole enzyme responsible for C-terminal amidation of peptides, and plays a role in the biosynthesis and regulation of multiple hormones, including proopiomelanocortin (POMC). Methods: We performed exome sequencing of germline and tumour DNA from 29 individuals with functioning pituitary adenomas (12 prolactinomas, 10 thyrotrophinomas, 7 cyclical Cushing's disease). An unfiltered analysis was undertaken of all PAM variants with population prevalence <5%. Results: We identified five coding, non-synonymous PAM variants of interest amongst seven individuals (six germline, one somatic). The five variants comprised four missense variants and one truncating variant, all heterozygous. Each variant had some evidence of pathogenicity based on population prevalence, conservation scores, in silico predictions and/or prior functional studies. The yield of predicted deleterious PAM variants was thus 7/29 (24%). The variants predominated in individuals with thyrotrophinomas (4/10, 40%) and cyclical Cushing's disease (2/7, 29%), compared to prolactinomas (1/12, 8%). Conclusion: This is the second study to demonstrate a high yield of suspected loss-of-function, predominantly germline, PAM variants in individuals with pathological pituitary hypersecretion. We have extended the association with corticotrophinoma to include the specific clinical entity of cyclical Cushing's disease and demonstrated a novel association between PAM variants and thyrotrophinoma. PAM variants might act as risk alleles for pituitary adenoma formation, with a possible genotype-phenotype relationship between truncating variants and altered temporal secretion of cortisol.