NRF2 Promotes Tumor Maintenance by Modulating mRNA Translation in Pancreatic Cancer.

Pancreatic cancer is a deadly malignancy that lacks effective therapeutics. We previously reported that oncogenic Kras induced the redox master regulator Nfe2l2/Nrf2 to stimulate pancreatic and lung cancer initiation. Here, we show that NRF2 is necessary to maintain pancreatic cancer proliferation by regulating mRNA translation. Specifically, loss of NRF2 led to defects in autocrine epidermal growth factor receptor (EGFR) signaling and oxidation of specific translational regulatory proteins, resulting in impaired cap-dependent and cap-independent mRNA translation in pancreatic cancer cells. Combined targeting of the EGFR effector AKT and the glutathione antioxidant pathway mimicked Nrf2 ablation to potently inhibit pancreatic cancer ex vivo and in vivo, representing a promising synthetic lethal strategy for treating the disease.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Organoid models of human and mouse ductal pancreatic cancer.

Pancreatic cancer is one of the most lethal malignancies due to its late diagnosis and limited response to treatment. Tractable methods to identify and interrogate pathways involved in pancreatic tumorigenesis are urgently needed. We established organoid models from normal and neoplastic murine and human pancreas tissues. Pancreatic organoids can be rapidly generated from resected tumors and biopsies, survive cryopreservation, and exhibit ductal- and disease-stage-specific characteristics. Orthotopically transplanted neoplastic organoids recapitulate the full spectrum of tumor development by forming early-grade neoplasms that progress to locally invasive and metastatic carcinomas. Due to their ability to be genetically manipulated, organoids are a platform to probe genetic cooperation. Comprehensive transcriptional and proteomic analyses of murine pancreatic organoids revealed genes and pathways altered during disease progression. The confirmation of many of these protein changes in human tissues demonstrates that organoids are a facile model system to discover characteristics of this deadly malignancy.

Suppression of protein tyrosine phosphatase N23 predisposes to breast tumorigenesis via activation of FYN kinase.

Disruption of the balanced modulation of reversible tyrosine phosphorylation has been implicated in the etiology of various human cancers, including breast cancer. Protein Tyrosine Phosphatase N23 (PTPN23) resides in chromosomal region 3p21.3, which is hemizygously or homozygously lost in some breast cancer patients. In a loss-of-function PTPome screen, our laboratory identified PTPN23 as a suppressor of cell motility and invasion in mammary epithelial and breast cancer cells. Now, our TCGA (The Cancer Genome Atlas) database analyses illustrate a correlation between low PTPN23 expression and poor survival in breast cancers of various subtypes. Therefore, we investigated the tumor-suppressive function of PTPN23 in an orthotopic transplantation mouse model. Suppression of PTPN23 in Comma 1Dß cells induced breast tumors within 56 wk. In PTPN23-depleted tumors, we detected hyperphosphorylation of the autophosphorylation site tyrosine in the SRC family kinase (SFK) FYN as well as Tyr142 in ß-catenin. We validated the underlying mechanism of PTPN23 function in breast tumorigenesis as that of a key phosphatase that normally suppresses the activity of FYN in two different models. We demonstrated that tumor outgrowth from PTPN23-deficient BT474 cells was suppressed in a xenograft model in vivo upon treatment with AZD0530, an SFK inhibitor. Furthermore, double knockout of FYN and PTPN23 via CRISPR/CAS9 also attenuated tumor outgrowth from PTPN23 knockout Cal51 cells. Overall, this mechanistic analysis of the tumor-suppressive function of PTPN23 in breast cancer supports the identification of FYN as a therapeutic target for breast tumors with heterozygous or homozygous loss of PTPN23.

piRNA-directed cleavage of meiotic transcripts regulates spermatogenesis.

MIWI catalytic activity is required for spermatogenesis, indicating that piRNA-guided cleavage is critical for germ cell development. To identify meiotic piRNA targets, we augmented the mouse piRNA repertoire by introducing a human meiotic piRNA cluster. This triggered a spermatogenesis defect by inappropriately targeting the piRNA machinery to mouse mRNAs essential for germ cell development. Analysis of such de novo targets revealed a signature for pachytene piRNA target recognition. This enabled identification of both transposable elements and meiotically expressed protein-coding genes as targets of native piRNAs. Cleavage of genic targets began at the pachytene stage and resulted in progressive repression through meiosis, driven at least in part via the ping-pong cycle. Our data support the idea that meiotic piRNA populations must be strongly selected to enable successful spermatogenesis, both driving the response away from essential genes and directing the pathway toward mRNA targets that are regulated by small RNAs in meiotic cells.

TGF-ß/Smad signaling through DOCK4 facilitates lung adenocarcinoma metastasis.

The mechanisms by which TGF-ß promotes lung adenocarcinoma (ADC) metastasis are largely unknown. Here, we report that in lung ADC cells, TGF-ß potently induces expression of DOCK4, but not other DOCK family members, via the Smad pathway and that DOCK4 induction mediates TGF-ß's prometastatic effects by enhancing tumor cell extravasation. TGF-ß-induced DOCK4 stimulates lung ADC cell protrusion, motility, and invasion without affecting epithelial-to-mesenchymal transition. These processes, which are fundamental to tumor cell extravasation, are driven by DOCK4-mediated Rac1 activation, unveiling a novel link between TGF-ß and Rac1. Thus, our findings uncover the atypical Rac1 activator DOCK4 as a key component of the TGF-ß/Smad pathway that promotes lung ADC cell extravasation and metastasis.

RNF17 blocks promiscuous activity of PIWI proteins in mouse testes.

PIWI proteins and their associated piRNAs protect germ cells from the activity of mobile genetic elements. Two classes of piRNAs­primary and secondary­are defined by their mechanisms of biogenesis. Primary piRNAs are processed directly from transcripts of piRNA cluster loci, whereas secondary piRNAs are generated in an adaptive amplification loop, termed the ping-pong cycle. In mammals, piRNA populations are dynamic, shifting as male germ cells develop. Embryonic piRNAs consist of both primary and secondary species and are mainly directed toward transposons. In meiotic cells, the piRNA population is transposon-poor and largely restricted to primary piRNAs derived from pachytene piRNA clusters. The transition from the embryonic to the adult piRNA pathway is not well understood. Here we show that RNF17 shapes adult meiotic piRNA content by suppressing the production of secondary piRNAs. In the absence of RNF17, ping-pong occurs inappropriately in meiotic cells. Ping-pong initiates piRNA responses against not only transposons but also protein-coding genes and long noncoding RNAs, including genes essential for germ cell development. Thus, the sterility of Rnf17 mutants may be a manifestation of a small RNA-based autoimmune reaction.

Single-cell RNA-seq analysis identifies markers of resistance to targeted BRAF inhibitors in melanoma cell populations.

Single-cell RNA-seq's (scRNA-seq) unprecedented cellular resolution at a genome-wide scale enables us to address questions about cellular heterogeneity that are inaccessible using methods that average over bulk tissue extracts. However, scRNA-seq data sets also present additional challenges such as high transcript dropout rates, stochastic transcription events, and complex population substructures. Here, we present a single-cell RNA-seq analysis and klustering evaluation (SAKE), a robust method for scRNA-seq analysis that provides quantitative statistical metrics at each step of the analysis pipeline. Comparing SAKE to multiple single-cell analysis methods shows that most methods perform similarly across a wide range of cellular contexts, with SAKE outperforming these methods in the case of large complex populations. We next applied the SAKE algorithms to identify drug-resistant cellular populations as human melanoma cells respond to targeted BRAF inhibitors (BRAFi). Single-cell RNA-seq data from both the Fluidigm C1 and 10x Genomics platforms were analyzed with SAKE to dissect this problem at multiple scales. Data from both platforms indicate that BRAF inhibitor-resistant cells can emerge from rare populations already present before drug application, with SAKE identifying both novel and known markers of resistance. These experimentally validated markers of BRAFi resistance share overlap with previous analyses in different melanoma cell lines, demonstrating the generality of these findings and highlighting the utility of single-cell analysis to elucidate mechanisms of BRAFi resistance.

Ten simple rules for running a successful women-in-STEM organization on an academic campus.

The current academic culture facing women in science, technology, engineering, and math (STEM) fields in the United States has sparked the formation of grassroots advocacy groups to empower female scientists in training. However, the impact of these initiatives often goes unmeasured and underappreciated. Our Women in Science and Engineering (WiSE) organization serves postdoctoral researchers, graduate students, and research technicians (trainees) at a private research institute for biological sciences. Here we propose the following guidelines for cultivating a successful women-in-STEM-focused group based upon survey results from our own scientific community as well as the experience of our WiSE group leaders. We hope these recommendations can provide guidance to advocacy groups at other research and academic organizations that wish to strengthen their efforts. Whereas our own group specifically focuses on the underrepresented state of women in science, we hope these guidelines may be adapted and applied to groups that advocate for any minority group within the greater scientific community (i.e., those of gender, race/ethnicity, socioeconomic background, sexual orientation, etc.).

Multiple roles for Piwi in silencing Drosophila transposons.

Silencing of transposons in the Drosophila ovary relies on three Piwi family proteins--Piwi, Aubergine (Aub), and Ago3--acting in concert with their small RNA guides, the Piwi-interacting RNAs (piRNAs). Aub and Ago3 are found in the germ cell cytoplasm, where they function in the ping-pong cycle to consume transposon mRNAs. The nuclear Piwi protein is required for transposon silencing in both germ and somatic follicle cells, yet the precise mechanisms by which Piwi acts remain largely unclear. We investigated the role of Piwi by combining cell type-specific knockdowns with measurements of steady-state transposon mRNA levels, nascent RNA synthesis, chromatin state, and small RNA abundance. In somatic cells, Piwi loss led to concerted effects on nascent transcripts and transposon mRNAs, indicating that Piwi acts through transcriptional gene silencing (TGS). In germ cells, Piwi loss showed disproportionate impacts on steady-state RNA levels, indicating that it also exerts an effect on post-transcriptional gene silencing (PTGS). Piwi knockdown affected levels of germ cell piRNAs presumably bound to Aub and Ago3, perhaps explaining its post-transcriptional impacts. Overall, our results indicate that Piwi plays multiple roles in the piRNA pathway, in part enforcing transposon repression through effects on local chromatin states and transcription but also participating in germ cell piRNA biogenesis.

Retrotransposon activation contributes to neurodegeneration in a Drosophila TDP-43 model of ALS.

Amyotrophic lateral sclerosis (ALS) and frontotemporal lobar degeneration (FTLD) are two incurable neurodegenerative disorders that exist on a symptomological spectrum and share both genetic underpinnings and pathophysiological hallmarks. Functional abnormality of TAR DNA-binding protein 43 (TDP-43), an aggregation-prone RNA and DNA binding protein, is observed in the vast majority of both familial and sporadic ALS cases and in ~40% of FTLD cases, but the cascade of events leading to cell death are not understood. We have expressed human TDP-43 (hTDP-43) in Drosophila neurons and glia, a model that recapitulates many of the characteristics of TDP-43-linked human disease including protein aggregation pathology, locomotor impairment, and premature death. We report that such expression of hTDP-43 impairs small interfering RNA (siRNA) silencing, which is the major post-transcriptional mechanism of retrotransposable element (RTE) control in somatic tissue. This is accompanied by de-repression of a panel of both LINE and LTR families of RTEs, with somewhat different elements being active in response to hTDP-43 expression in glia versus neurons. hTDP-43 expression in glia causes an early and severe loss of control of a specific RTE, the endogenous retrovirus (ERV) gypsy. We demonstrate that gypsy causes the degenerative phenotypes in these flies because we are able to rescue the toxicity of glial hTDP-43 either by genetically blocking expression of this RTE or by pharmacologically inhibiting RTE reverse transcriptase activity. Moreover, we provide evidence that activation of DNA damage-mediated programmed cell death underlies both neuronal and glial hTDP-43 toxicity, consistent with RTE-mediated effects in both cell types. Our findings suggest a novel mechanism in which RTE activity contributes to neurodegeneration in TDP-43-mediated diseases such as ALS and FTLD.

Novel DICER-LIKE1 siRNAs Bypass the Requirement for DICER-LIKE4 in Maize Development.

Dicer enzymes function at the core of RNA silencing to defend against exogenous RNA or to regulate endogenous genes. Plant DICER-LIKE4 (DCL4) performs dual functions, acting in antiviral defense and in development via the biogenesis of trans-acting short-interfering RNAs (siRNAs) termed tasiR-ARFs. These small RNAs play an essential role in the grasses, spatially defining the expression domain of AUXIN RESPONSE FACTOR3 (ARF3) transcription factors. However, contrary to tasiR-ARFs' essential function in development, DCL4 proteins exhibit strong evidence of recurrent adaptation typical of host factors involved in antiviral immunity. Here, we address how DCL4 balances its role in development with pressures to diversify in response to viral attack. We show that, in contrast to other tasiR-ARF biogenesis mutants, dcl4 null alleles have an uncharacteristically mild phenotype, correlated with normal expression of select arf3 targets. Loss of DCL4 activity yields a class of 22-nucleotide tasiR-ARF variants associated with the processing of arf3 transcripts into 22-nucleotide secondary siRNAs by DCL1. Our findings reveal a DCL1-dependent siRNA pathway that bypasses the otherwise adverse developmental effects of mutations in DCL4. This pathway is predicted to have important implications for DCL4's role in antiviral defense by reducing the selective constraints on DCL4 and allowing it to diversify in response to viral suppressors.

G-protein-coupled receptor GPR161 is overexpressed in breast cancer and is a promoter of cell proliferation and invasion.

Triple-negative breast cancer (TNBC) accounts for 20% of breast cancer in women and lacks an effective targeted therapy. Therefore, finding common vulnerabilities in these tumors represents an opportunity for more effective treatment. Despite the growing appreciation of G-protein-coupled receptor (GPCR)-mediated signaling in cancer pathogenesis, very little is known about the role GPCRs play in TNBC. Using genomic information of human breast cancer, we have discovered that the orphan GPCR, G-protein-coupled receptor 161 (GPR161) is overexpressed specifically in TNBC and correlates with poor prognosis. Knockdown of GPR161 impairs proliferation of human basal breast cancer cell lines. Overexpression of GPR161 in human mammary epithelial cells increases cell proliferation, migration, intracellular accumulation of E-cadherin, and formation of multiacinar structures in 3D culture. GPR161 forms a signaling complex with the scaffold proteins ß-arrestin 2 and Ile Gln motif containing GTPase Activating Protein 1, a regulator of mammalian target of rapamycin complex 1 and E-cadherin. Consistently, GPR161 amplified breast tumors and cells overexpressing GPR161 activate mammalian target of rapamycin signaling and decrease Ile Gln motif containing GTPase Activating Protein 1 phosphorylation. Thus, we identify the orphan GPCR, GPR161, as an important regulator and a potential drug target for TNBC.

Genome-wide analysis of leafbladeless1-regulated and phased small RNAs underscores the importance of the TAS3 ta-siRNA pathway to maize development.

Maize leafbladeless1 (lbl1) encodes a key component in the trans-acting short-interfering RNA (ta-siRNA) biogenesis pathway. Correlated with a great diversity in ta-siRNAs and the targets they regulate, the phenotypes conditioned by mutants perturbing this small RNA pathway vary extensively across species. Mutations in lbl1 result in severe developmental defects, giving rise to plants with radial, abaxialized leaves. To investigate the basis for this phenotype, we compared the small RNA content between wild-type and lbl1 seedling apices. We show that LBL1 affects the accumulation of small RNAs in all major classes, and reveal unexpected crosstalk between ta-siRNA biogenesis and other small RNA pathways regulating transposons. Interestingly, in contrast to data from other plant species, we found no evidence for the existence of phased siRNAs generated via the one-hit model. Our analysis identified nine TAS loci, all belonging to the conserved TAS3 family. Information from RNA deep sequencing and PARE analyses identified the tasiR-ARFs as the major functional ta-siRNAs in the maize vegetative apex where they regulate expression of AUXIN RESPONSE FACTOR3 (ARF3) homologs. Plants expressing a tasiR-ARF insensitive arf3a transgene recapitulate the phenotype of lbl1, providing direct evidence that deregulation of ARF3 transcription factors underlies the developmental defects of maize ta-siRNA biogenesis mutants. The phenotypes of Arabidopsis and Medicago ta-siRNA mutants, while strikingly different, likewise result from misexpression of the tasiR-ARF target ARF3. Our data indicate that diversity in TAS pathways and their targets cannot fully account for the phenotypic differences conditioned by ta-siRNA biogenesis mutants across plant species. Instead, we propose that divergence in the gene networks downstream of the ARF3 transcription factors or the spatiotemporal pattern during leaf development in which these proteins act constitute key factors underlying the distinct contributions of the ta-siRNA pathway to development in maize, Arabidopsis, and possibly other plant species as well.

An atlas of chromatoid body components.

The genome of male germ cells is actively transcribed during spermatogenesis to produce phase-specific protein-coding mRNAs and a considerable amount of different noncoding RNAs. Ribonucleoprotein (RNP) granule-mediated RNA regulation provides a powerful means to secure the quality and correct expression of the requisite transcripts. Haploid spermatids are characterized by a unique, unusually large cytoplasmic granule, the chromatoid body (CB), which emerges during the switch between the meiotic and post-meiotic phases of spermatogenesis. To better understand the role of the CB in male germ cell differentiation, we isolated CBs from mouse testes and revealed its full RNA and protein composition. We showed that the CB is mainly composed of RNA-binding proteins and other proteins involved RNA regulation. The CB was loaded with RNA, including pachytene piRNAs, a diverse set of mRNAs, and a number of uncharacterized long noncoding transcripts. The CB was demonstrated to accumulate nascent RNA during all the steps of round spermatid differentiation. Our results revealed the CB as a large germ cell-specific RNP platform that is involved in the control of the highly complex transcriptome of haploid male germ cells.

TEtranscripts: a package for including transposable elements in differential expression analysis of RNA-seq datasets.

MOTIVATION: Most RNA-seq data analysis software packages are not designed to handle the complexities involved in properly apportioning short sequencing reads to highly repetitive regions of the genome. These regions are often occupied by transposable elements (TEs), which make up between 20 and 80% of eukaryotic genomes. They can contribute a substantial portion of transcriptomic and genomic sequence reads, but are typically ignored in most analyses. RESULTS: Here, we present a method and software package for including both gene- and TE-associated ambiguously mapped reads in differential expression analysis. Our method shows improved recovery of TE transcripts over other published expression analysis methods, in both synthetic data and qPCR/NanoString-validated published datasets. AVAILABILITY AND IMPLEMENTATION: The source code, associated GTF files for TE annotation, and testing data are freely available at http://hammelllab.labsites.cshl.edu/software. CONTACT: mhammell@cshl.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

TDP-43 dysregulation of polyadenylation site selection is a defining feature of RNA misprocessing in ALS/FTD and related disorders.

Nuclear clearance and cytoplasmic aggregation of the RNA-binding protein TDP-43 are observed in many neurodegenerative disorders, including amyotrophic lateral sclerosis (ALS) and fronto- temporal dementia (FTD). Although TDP-43 dysregulation of splicing has emerged as a key event in these diseases, TDP-43 can also regulate polyadenylation; yet, this has not been adequately studied. Here, we applied the dynamic analysis of polyadenylation from RNA-seq (DaPars) tool to ALS/FTD transcriptome datasets, and report extensive alternative polyadenylation (APA) upon TDP-43 alteration in ALS/FTD cell models and postmortem ALS/FTD neuronal nuclei. Importantly, many identified APA genes highlight pathways implicated in ALS/FTD pathogenesis. To determine the functional significance of APA elicited by TDP-43 nuclear depletion, we examined microtubule affinity regulating kinase 3 (MARK3). Nuclear loss of TDP-43 yielded increased expression of MARK3 transcripts with longer 3'UTRs, resulting in greater transcript stability and elevated MARK3 protein levels, which promotes increased neuronal tau S262 phosphorylation. Our findings define changes in polyadenylation site selection as a previously unrecognized feature of TDP-43-driven disease pathology in ALS/FTD and highlight a potentially novel mechanistic link between TDP-43 dysfunction and tau regulation.

Effect of life history on microRNA expression during C. elegans development.

Animals have evolved mechanisms to ensure the robustness of developmental outcomes to changing environments. MicroRNA expression may contribute to developmental robustness because microRNAs are key post-transcriptional regulators of developmental gene expression and can affect the expression of multiple target genes. Caenorhabditis elegans provides an excellent model to study developmental responses to environmental conditions. In favorable environments, C. elegans larvae develop rapidly and continuously through four larval stages. In contrast, in unfavorable conditions, larval development may be interrupted at either of two diapause stages: The L1 diapause occurs when embryos hatch in the absence of food, and the dauer diapause occurs after the second larval stage in response to environmental stimuli encountered during the first two larval stages. Dauer larvae are stress resistant and long lived, permitting survival in harsh conditions. When environmental conditions improve, dauer larvae re-enter development, and progress through two post-dauer larval stages to adulthood. Strikingly, all of these life history options (whether continuous or interrupted) involve an identical pattern and sequence of cell division and cell fates. To identify microRNAs with potential functions in buffering development in the context of C. elegans life history options, we used multiplex real-time PCR to assess the expression of 107 microRNAs throughout development in both continuous and interrupted life histories. We identified 17 microRNAs whose developmental profile of expression is affected by dauer life history and/or L1 diapause, compared to continuous development. Hence these microRNAs could function to regulate gene expression programs appropriate for different life history options in the developing worm.

Single-cell transcriptomics of human traumatic brain injury reveals activation of endogenous retroviruses in oligodendroglia.

Traumatic brain injury (TBI) is a leading cause of chronic brain impairment and results in a robust, but poorly understood, neuroinflammatory response that contributes to the long-term pathology. We used single-nuclei RNA sequencing (snRNA-seq) to study transcriptomic changes in different cell populations in human brain tissue obtained acutely after severe, life-threatening TBI. This revealed a unique transcriptional response in oligodendrocyte precursors and mature oligodendrocytes, including the activation of a robust innate immune response, indicating an important role for oligodendroglia in the initiation of neuroinflammation. The activation of an innate immune response correlated with transcriptional upregulation of endogenous retroviruses in oligodendroglia. This observation was causally linked in vitro using human glial progenitors, implicating these ancient viral sequences in human neuroinflammation. In summary, this work provides insight into the initiating events of the neuroinflammatory response in TBI, which has therapeutic implications.

LINE-1 retrotransposons drive human neuronal transcriptome complexity and functional diversification.

The genetic mechanisms underlying the expansion in size and complexity of the human brain remain poorly understood. Long interspersed nuclear element-1 (L1) retrotransposons are a source of divergent genetic information in hominoid genomes, but their importance in physiological functions and their contribution to human brain evolution are largely unknown. Using multiomics profiling, we here demonstrate that L1 promoters are dynamically active in the developing and the adult human brain. L1s generate hundreds of developmentally regulated and cell type-specific transcripts, many that are co-opted as chimeric transcripts or regulatory RNAs. One L1-derived long noncoding RNA, LINC01876, is a human-specific transcript expressed exclusively during brain development. CRISPR interference silencing of LINC01876 results in reduced size of cerebral organoids and premature differentiation of neural progenitors, implicating L1s in human-specific developmental processes. In summary, our results demonstrate that L1-derived transcripts provide a previously undescribed layer of primate- and human-specific transcriptome complexity that contributes to the functional diversification of the human brain.