Mitochondrial Bioenergetics of Functional Wound Closure is Dependent on Macrophage-Keratinocyte Exosomal Crosstalk.

Tissue nanotransfection (TNT)-based fluorescent labeling of cell-specific exosomes has shown that exosomes play a central role in physiological keratinocyte-macrophage (mϕ) crosstalk at the wound-site. Here, we report that during the early phase of wound reepithelialization, macrophage-derived exosomes (Exomϕ), enriched with the outer mitochondrial membrane protein TOMM70, are localized in leading-edge keratinocytes. TOMM70 is a 70 kDa adaptor protein anchored in the mitochondrial outer membrane and plays a critical role in maintaining mitochondrial function and quality. TOMM70 selectively recognizes cytosolic chaperones by its tetratricopeptide repeat (TPR) domain and facilitates the import of preproteins lacking a positively charged mitochondrial targeted sequence. Exosomal packaging of TOMM70 in mϕ was independent of mitochondrial fission. TOMM70-enriched Exomϕ compensated for the hypoxia-induced depletion of epidermal TOMM70, thereby rescuing mitochondrial metabolism in leading-edge keratinocytes. Thus, macrophage-derived TOMM70 is responsible for the glycolytic ATP supply to power keratinocyte migration. Blockade of exosomal uptake from keratinocytes impaired wound closure with the persistence of proinflammatory mϕ in the wound microenvironment, pointing toward a bidirectional crosstalk between these two cell types. The significance of such bidirectional crosstalk was established by the observation that in patients with nonhealing diabetic foot ulcers, TOMM70 is deficient in keratinocytes of wound-edge tissues.

Monoallelically expressed noncoding RNAs form nucleolar territories on NOR-containing chromosomes and regulate rRNA expression.

Out of the several hundred copies of rRNA genes arranged in the nucleolar organizing regions (NOR) of the five human acrocentric chromosomes, ~50% remain transcriptionally inactive. NOR-associated sequences and epigenetic modifications contribute to the differential expression of rRNAs. However, the mechanism(s) controlling the dosage of active versus inactive rRNA genes within each NOR in mammals is yet to be determined. We have discovered a family of ncRNAs, SNULs (Single NUcleolus Localized RNA), which form constrained sub-nucleolar territories on individual NORs and influence rRNA expression. Individual members of the SNULs monoallelically associate with specific NOR-containing chromosomes. SNULs share sequence similarity to pre-rRNA and localize in the sub-nucleolar compartment with pre-rRNA. Finally, SNULs control rRNA expression by influencing pre-rRNA sorting to the DFC compartment and pre-rRNA processing. Our study discovered a novel class of ncRNAs influencing rRNA expression by forming constrained nucleolar territories on individual NORs.

Human fetal dermal fibroblast-myeloid cell diversity is characterized by dominance of pro-healing Annexin1-FPR1 signaling.

Fetal skin achieves scarless wound repair. Dermal fibroblasts play a central role in extracellular matrix deposition and scarring outcomes. Both fetal and gingival wound repair share minimal scarring outcomes. We tested the hypothesis that compared to adult skin fibroblasts, human fetal skin fibroblast diversity is unique and partly overlaps with gingival skin fibroblasts. Human fetal skin (FS, n = 3), gingiva (HGG, n = 13), and mature skin (MS, n = 13) were compared at single-cell resolution. Dermal fibroblasts, the most abundant cluster, were examined to establish a connectome with other skin cells. Annexin1-FPR1 signaling pathway was dominant in both FS as well as HGG fibroblasts and related myeloid cells while scanty in MS fibroblasts. Myeloid-specific FPR1-ORF delivered in murine wound edge using tissue nanotransfection (TNT) technology significantly enhanced the quality of healing. Pseudotime analyses identified the co-existence of an HGG fibroblast subset with FPR1high myeloid cells of fetal origin indicating common underlying biological processes.

The Prolonged Terminal Phase of Human Life Induces Survival Response in the Skin Transcriptome.

Human death marks the end of organismal life under conditions such that the components of the human body continue to be alive. Such postmortem cellular survival depends on the nature (Hardy scale of slow-fast death) of human death. Slow and expected death typically results from terminal illnesses and includes a prolonged terminal phase of life. As such organismal death process unfolds, do cells of the human body adapt for postmortem cellular survival? Organs with low energy cost-of-living, such as the skin, are better suited for postmortem cellular survival. In this work, the effect of different durations of terminal phase of human life on postmortem changes in cellular gene expression was investigated using RNA sequencing data of 701 human skin samples from the Genotype-Tissue Expression (GTEx) database. Longer terminal phase (slow-death) was associated with a more robust induction of survival pathways (PI3K-Akt signaling) in postmortem skin. Such cellular survival response was associated with the upregulation of embryonic developmental transcription factors such as FOXO1 , FOXO3 , ATF4 and CEBPD . Upregulation of PI3K-Akt signaling was independent of sex or duration of death-related tissue ischemia. Analysis of single nucleus RNA-seq of post-mortem skin tissue specifically identified the dermal fibroblast compartment to be most resilient as marked by adaptive induction of PI3K-Akt signaling. In addition, slow death also induced angiogenic pathways in the dermal endothelial cell compartment of postmortem human skin. In contrast, specific pathways supporting functional properties of the skin as an organ were downregulated following slow death. Such pathways included melanogenesis and those representing the skin extracellular matrix (collagen expression and metabolism). Efforts to understand the significance of death as a biological variable (DABV) in influencing the transcriptomic composition of surviving component tissues has far-reaching implications including rigorous interpretation of experimental data collected from the dead and mechanisms involved in transplant-tissue obtained from dead donors.

Identification of a physiologic vasculogenic fibroblast state to achieve tissue repair.

Tissue injury to skin diminishes miR-200b in dermal fibroblasts. Fibroblasts are widely reported to directly reprogram into endothelial-like cells and we hypothesized that miR-200b inhibition may cause such changes. We transfected human dermal fibroblasts with anti-miR-200b oligonucleotide, then using single cell RNA sequencing, identified emergence of a vasculogenic subset with a distinct fibroblast transcriptome and demonstrated blood vessel forming function in vivo. Anti-miR-200b delivery to murine injury sites likewise enhanced tissue perfusion, wound closure, and vasculogenic fibroblast contribution to perfused vessels in a FLI1 dependent manner. Vasculogenic fibroblast subset emergence was blunted in delayed healing wounds of diabetic animals but, topical tissue nanotransfection of a single anti-miR-200b oligonucleotide was sufficient to restore FLI1 expression, vasculogenic fibroblast emergence, tissue perfusion, and wound healing. Augmenting a physiologic tissue injury adaptive response mechanism that produces a vasculogenic fibroblast state change opens new avenues for therapeutic tissue vascularization of ischemic wounds.

Neuron-oligodendrocyte potassium shuttling at nodes of Ranvier protects against inflammatory demyelination.

Multiple sclerosis (MS) is a progressive inflammatory demyelinating disease of the CNS. Increasing evidence suggests that vulnerable neurons in MS exhibit fatal metabolic exhaustion over time, a phenomenon hypothesized to be caused by chronic hyperexcitability. Axonal Kv7 (outward-rectifying) and oligodendroglial Kir4.1 (inward-rectifying) potassium channels have important roles in regulating neuronal excitability at and around the nodes of Ranvier. Here, we studied the spatial and functional relationship between neuronal Kv7 and oligodendroglial Kir4.1 channels and assessed the transcriptional and functional signatures of cortical and retinal projection neurons under physiological and inflammatory demyelinating conditions. We found that both channels became dysregulated in MS and experimental autoimmune encephalomyelitis (EAE), with Kir4.1 channels being chronically downregulated and Kv7 channel subunits being transiently upregulated during inflammatory demyelination. Further, we observed that pharmacological Kv7 channel opening with retigabine reduced neuronal hyperexcitability in human and EAE neurons, improved clinical EAE signs, and rescued neuronal pathology in oligodendrocyte-Kir4.1-deficient (OL-Kir4.1-deficient) mice. In summary, our findings indicate that neuron-OL compensatory interactions promoted resilience through Kv7 and Kir4.1 channels and identify pharmacological activation of nodal Kv7 channels as a neuroprotective strategy against inflammatory demyelination.

Genome-wide DNA hypermethylation opposes healing in patients with chronic wounds by impairing epithelial-mesenchymal transition.

An extreme chronic wound tissue microenvironment causes epigenetic gene silencing. An unbiased whole-genome methylome was studied in the wound-edge tissue of patients with chronic wounds. A total of 4,689 differentially methylated regions (DMRs) were identified in chronic wound-edge skin compared with unwounded human skin. Hypermethylation was more frequently observed (3,661 DMRs) in the chronic wound-edge tissue compared with hypomethylation (1,028 DMRs). Twenty-six hypermethylated DMRs were involved in epithelial-mesenchymal transition (EMT). Bisulfite sequencing validated hypermethylation of a predicted specific upstream regulator TP53. RNA-Seq analysis was performed to qualify findings from methylome analysis. Analysis of the downregulated genes identified the TP53 signaling pathway as being significantly silenced. Direct comparison of hypermethylation and downregulated genes identified 4 genes, ADAM17, NOTCH, TWIST1, and SMURF1, that functionally represent the EMT pathway. Single-cell RNA-Seq studies revealed that these effects on gene expression were limited to the keratinocyte cell compartment. Experimental murine studies established that tissue ischemia potently induces wound-edge gene methylation and that 5'-azacytidine, inhibitor of methylation, improved wound closure. To specifically address the significance of TP53 methylation, keratinocyte-specific editing of TP53 methylation at the wound edge was achieved by a tissue nanotransfection-based CRISPR/dCas9 approach. This work identified that reversal of methylation-dependent keratinocyte gene silencing represents a productive therapeutic strategy to improve wound closure.

P2R Inhibitors Prevent Antibody-Mediated Complement Activation in an Animal Model of Neuromyelitis Optica : P2R Inhibitors Prevent Autoantibody Injury.

Purinergic 2 receptors (P2Rs) contribute to disease-related immune cell signaling and are upregulated in various pathological settings, including neuroinflammation. P2R inhibitors have been used to treat inflammatory diseases and can protect against complement-mediated cell injury. However, the mechanisms behind these anti-inflammatory properties of P2R inhibitors are not well understood, and their potential in CNS autoimmunity is underexplored. Here, we tested the effects of P2R inhibitors on glial toxicity in a mouse model of neuromyelitis optica spectrum disorder (NMOSD). NMOSD is a destructive CNS autoimmune disorder, in which autoantibodies against astrocytic surface antigen Aquaporin 4 (AQP4) mediate complement-dependent loss of astrocytes. Using two-photon microscopy in vivo, we found that various classes of P2R inhibitors prevented AQP4-IgG/complement-dependent astrocyte death. In vitro, these drugs inhibited the binding of AQP4-IgG or MOG-IgG to their antigen in a dose-dependent manner. Size-exclusion chromatography and circular dichroism spectroscopy revealed a partial unfolding of antibodies in the presence of various P2R inhibitors, suggesting a shared interference with IgG antibodies leading to their conformational change. Our study demonstrates that P2R inhibitors can disrupt complement activation by direct interaction with IgG. This mechanism is likely to influence the role of P2R inhibitors in autoimmune disease models and their therapeutic impact in human disease.

CASowary: CRISPR-Cas13 guide RNA predictor for transcript depletion.

BACKGROUND: Recent discovery of the gene editing system - CRISPR (Clustered Regularly Interspersed Short Palindromic Repeats) associated proteins (Cas), has resulted in its widespread use for improved understanding of a variety of biological systems. Cas13, a lesser studied Cas protein, has been repurposed to allow for efficient and precise editing of RNA molecules. The Cas13 system utilizes base complementarity between a crRNA/sgRNA (crispr RNA or single guide RNA) and a target RNA transcript, to preferentially bind to only the target transcript. Unlike targeting the upstream regulatory regions of protein coding genes on the genome, the transcriptome is significantly more redundant, leading to many transcripts having wide stretches of identical nucleotide sequences. Transcripts also exhibit complex three-dimensional structures and interact with an array of RBPs (RNA Binding Proteins), both of which may impact the effectiveness of transcript depletion of target sequences. However, our understanding of the features and corresponding methods which can predict whether a specific sgRNA will effectively knockdown a transcript is very limited. RESULTS: Here we present a novel machine learning and computational tool, CASowary, to predict the efficacy of a sgRNA. We used publicly available RNA knockdown data from Cas13 characterization experiments for 555 sgRNAs targeting the transcriptome in HEK293 cells, in conjunction with transcriptome-wide protein occupancy information. Our model utilizes a Decision Tree architecture with a set of 112 sequence and target availability features, to classify sgRNA efficacy into one of four classes, based upon expected level of target transcript knockdown. After accounting for noise in the training data set, the noise-normalized accuracy exceeds 70%. Additionally, highly effective sgRNA predictions have been experimentally validated using an independent RNA targeting Cas system - CIRTS, confirming the robustness and reproducibility of our model's sgRNA predictions. Utilizing transcriptome wide protein occupancy map generated using POP-seq in HeLa cells against publicly available protein-RNA interaction map in Hek293 cells, we show that CASowary can predict high quality guides for numerous transcripts in a cell line specific manner. CONCLUSIONS: Application of CASowary to whole transcriptomes should enable rapid deployment of CRISPR/Cas13 systems, facilitating the development of therapeutic interventions linked with aberrations in RNA regulatory processes.

Endothelial Phospholipase Cγ2 Improves Outcomes of Diabetic Ischemic Limb Rescue Following VEGF Therapy.

Therapeutic vascular endothelial growth factor (VEGF) replenishment has met with limited success for the management of critical limb-threatening ischemia. To improve outcomes of VEGF therapy, we applied single-cell RNA sequencing (scRNA-seq) technology to study the endothelial cells of the human diabetic skin. Single-cell suspensions were generated from the human skin followed by cDNA preparation using the Chromium Next GEM Single-cell 3' Kit v3.1. Using appropriate quality control measures, 36,487 cells were chosen for downstream analysis. scRNA-seq studies identified that although VEGF signaling was not significantly altered in diabetic versus nondiabetic skin, phospholipase Cγ2 (PLCγ2) was downregulated. The significance of PLCγ2 in VEGF-mediated increase in endothelial cell metabolism and function was assessed in cultured human microvascular endothelial cells. In these cells, VEGF enhanced mitochondrial function, as indicated by elevation in oxygen consumption rate and extracellular acidification rate. The VEGF-dependent increase in cell metabolism was blunted in response to PLCγ2 inhibition. Follow-up rescue studies therefore focused on understanding the significance of VEGF therapy in presence or absence of endothelial PLCγ2 in type 1 (streptozotocin-injected) and type 2 (db/db) diabetic ischemic tissue. Nonviral topical tissue nanotransfection technology (TNT) delivery of CDH5 promoter-driven PLCγ2 open reading frame promoted the rescue of hindlimb ischemia in diabetic mice. Improvement of blood flow was also associated with higher abundance of VWF+/CD31+ and VWF+/SMA+ immunohistochemical staining. TNT-based gene delivery was not associated with tissue edema, a commonly noted complication associated with proangiogenic gene therapies. Taken together, our study demonstrates that TNT-mediated delivery of endothelial PLCγ2, as part of combination gene therapy, is effective in diabetic ischemic limb rescue.

Lantern: an integrative repository of functional annotations for lncRNAs in the human genome.

BACKGROUND: With advancements in omics technologies, the range of biological processes where long non-coding RNAs (lncRNAs) are involved, is expanding extensively, thereby generating the need to develop lncRNA annotation resources. Although, there are a plethora of resources for annotating genes, despite the extensive corpus of lncRNA literature, the available resources with lncRNA ontology annotations are rare. RESULTS: We present a lncRNA annotation extractor and repository (Lantern), developed using PubMed's abstract retrieval engine and NCBO's recommender annotation system. Lantern's annotations were benchmarked against lncRNAdb's manually curated free text. Benchmarking analysis suggested that Lantern has a recall of 0.62 against lncRNAdb for 182 lncRNAs and precision of 0.8. Additionally, we also annotated lncRNAs with multiple omics annotations, including predicted cis-regulatory TFs, interactions with RBPs, tissue-specific expression profiles, protein co-expression networks, coding potential, sub-cellular localization, and SNPs for ~ 11,000 lncRNAs in the human genome, providing a one-stop dynamic visualization platform. CONCLUSIONS: Lantern integrates a novel, accurate semi-automatic ontology annotation engine derived annotations combined with a variety of multi-omics annotations for lncRNAs, to provide a central web resource for dissecting the functional dynamics of long non-coding RNAs and to facilitate future hypothesis-driven experiments. The annotation pipeline and a web resource with current annotations for human lncRNAs are freely available on sysbio.lab.iupui.edu/lantern.

Transcriptome-wide high-throughput mapping of protein-RNA occupancy profiles using POP-seq.

Interaction between proteins and RNA is critical for post-transcriptional regulatory processes. Existing high throughput methods based on crosslinking of the protein-RNA complexes and poly-A pull down are reported to contribute to biases and are not readily amenable for identifying interaction sites on non poly-A RNAs. We present Protein Occupancy Profile-Sequencing (POP-seq), a phase separation based method in three versions, one of which does not require crosslinking, thus providing unbiased protein occupancy profiles on whole cell transcriptome without the requirement of poly-A pulldown. Our study demonstrates that ~ 68% of the total POP-seq peaks exhibited an overlap with publicly available protein-RNA interaction profiles of 97 RNA binding proteins (RBPs) in K562 cells. We show that POP-seq variants consistently capture protein-RNA interaction sites across a broad range of genes including on transcripts encoding for transcription factors (TFs), RNA-Binding Proteins (RBPs) and long non-coding RNAs (lncRNAs). POP-seq identified peaks exhibited a significant enrichment (p value < 2.2e-16) for GWAS SNPs, phenotypic, clinically relevant germline as well as somatic variants reported in cancer genomes, suggesting the prevalence of uncharacterized genomic variation in protein occupied sites on RNA. We demonstrate that the abundance of POP-seq peaks increases with an increase in expression of lncRNAs, suggesting that highly expressed lncRNA are likely to act as sponges for RBPs, contributing to the rewiring of protein-RNA interaction network in cancer cells. Overall, our data supports POP-seq as a robust and cost-effective method that could be applied to primary tissues for mapping global protein occupancies.

Manganese-induced oxidative stress, ultrastructural changes, and proteomics studies in rice plants.

Manganese (Mn) is an essential element for plant growth but it becomes phytotoxic at higher concentrations. The effect of Mn-excess in hydroponics medium was examined on growth, oxidative stress, and ultrastructural changes in chloroplasts and mitochondria as well proteomic alterations in rice (Oryza sativa L.) seedlings. Seedlings grown with 1 mM and 2 mM Mn in nutrient medium for 8 days showed decline in length and fresh biomass, and decline in net photosynthetic rate, transpiration rate, and stomatal conductance. Shoots of the seedlings had higher Mn content than roots. Mn-treated seedlings showed increased production of O2·-, H2O2, and .OH, increased lipid peroxidation, increased carbonylation of proteins, and increased proteolytic activity compared to untreated seedlings. Mn-treated seedlings showed disorganization and swelling of chloroplasts with appearance of plastoglobuli in TEM images and deformity in shape of mitochondria. Using confocal microscopy depolarization of mitochondrial membrane was observed marked by green fluorescence of JC-1 dye monomers in Mn-treated roots. Proteomics studies from leaves of Mn-treated seedlings involving 2DE and PDQuest analysis showed differential expression of 23 proteins, among which MALDI-TOF/TOF mass spectrometry analysis revealed Mn-led downregulation of photosynthesis-related proteins, namely oxygen-evolving complex protein associated with PSII, PAP-3, enzyme involved in protein folding peptidyl-prolyl cis-trans isomerase (PPIase) and carbohydrate metabolizing enzymes hydrolase, fructose-bisphosphate aldolase, transketolase, and isocitrate dehydrogenase, whereas ATP-dependent Clp protease, peroxidase, and nucleic acid-binding proteins were downregulated due to Mn treatment. Results indicate that Mn-excess inhibits growth of rice plants with induction of oxidative stress, causing structural alterations in chloroplasts, mitochondria, inhibiting photosynthesis, and downregulating many photosynthesis and carbohydrate metabolism-related proteins.

Role of SARS-CoV-2 in Altering the RNA-Binding Protein and miRNA-Directed Post-Transcriptional Regulatory Networks in Humans.

The outbreak of a novel coronavirus SARS-CoV-2 responsible for the COVID-19 pandemic has caused a worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2-mediated alterations on post-transcriptional gene regulations across human tissues remain elusive. In this study, we analyzed publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of the human post-transcriptional regulatory networks governed by RNA-binding proteins (RBPs) and micro-RNAs (miRs) due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2-encoded proteins directly interacted with 51 human RBPs, of which the majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2-infected lung cells that revealed enrichment for the immune response, cytokine-mediated signaling, and metabolism-associated genes. This study also characterized the alternative splicing events in SARS-CoV-2-infected cells compared to the control, demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to the viral infection. A motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs, suggesting the sponging of RBPs by the SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate the host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provided a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs across tissue types during a SARS-CoV-2 infection.

Granzyme A-producing T helper cells are critical for acute graft-versus-host disease.

Acute graft-versus-host disease (aGVHD) can occur after hematopoietic cell transplant in patients undergoing treatment for hematological malignancies or inborn errors. Although CD4+ T helper (Th) cells play a major role in aGVHD, the mechanisms by which they contribute, particularly within the intestines, have remained elusive. We have identified a potentially novel subset of Th cells that accumulated in the intestines and produced the serine protease granzyme A (GrA). GrA+ Th cells were distinct from other Th lineages and exhibited a noncytolytic phenotype. In vitro, GrA+ Th cells differentiated in the presence of IL-4, IL-6, and IL-21 and were transcriptionally unique from cells cultured with either IL-4 or the IL-6/IL-21 combination alone. In vivo, both STAT3 and STAT6 were required for GrA+ Th cell differentiation and played roles in maintenance of the lineage identity. Importantly, GrA+ Th cells promoted aGVHD-associated morbidity and mortality and contributed to crypt destruction within intestines but were not required for the beneficial graft-versus-leukemia effect. Our data indicate that GrA+ Th cells represent a distinct Th subset and are critical mediators of aGVHD.

Role of SARS-CoV-2 in altering the RNA binding protein and miRNA directed post-transcriptional regulatory networks in humans.

The outbreak of a novel coronavirus SARS-CoV-2 responsible for COVID-19 pandemic has caused worldwide public health emergency. Due to the constantly evolving nature of the coronaviruses, SARS-CoV-2 mediated alteration on post-transcriptional gene regulation across human tissues remains elusive. In this study, we analyze publicly available genomic datasets to systematically dissect the crosstalk and dysregulation of human post-transcriptional regulatory networks governed by RNA binding proteins (RBPs) and micro-RNAs (miRs), due to SARS-CoV-2 infection. We uncovered that 13 out of 29 SARS-CoV-2 encoded proteins directly interact with 51 human RBPs of which majority of them were abundantly expressed in gonadal tissues and immune cells. We further performed a functional analysis of differentially expressed genes in mock-treated versus SARS-CoV-2 infected lung cells that revealed enrichment for immune response, cytokine-mediated signaling, and metabolism associated genes. This study also characterized the alternative splicing events in SARS-CoV-2 infected cells compared to control demonstrating that skipped exons and mutually exclusive exons were the most abundant events that potentially contributed to differential outcomes in response to viral infection. Motif enrichment analysis on the RNA genomic sequence of SARS-CoV-2 clearly revealed the enrichment for RBPs such as SRSFs, PCBPs, ELAVs, and HNRNPs suggesting the sponging of RBPs by SARS-CoV-2 genome. A similar analysis to study the interactions of miRs with SARS-CoV-2 revealed functionally important miRs that were highly expressed in immune cells, suggesting that these interactions may contribute to the progression of the viral infection and modulate host immune response across other human tissues. Given the need to understand the interactions of SARS-CoV-2 with key post-transcriptional regulators in the human genome, this study provides a systematic computational analysis to dissect the role of dysregulated post-transcriptional regulatory networks controlled by RBPs and miRs, across tissues types during SARS-CoV-2 infection.

Targeting Bim via a lncRNA Morrbid Regulates the Survival of Preleukemic and Leukemic Cells.

Inhibition of anti-apoptotic proteins BCL-2 and MCL-1 to release pro-apoptotic protein BIM and reactivate cell death could potentially be an efficient strategy for the treatment of leukemia. Here, we show that a lncRNA, MORRBID, a selective transcriptional repressor of BIM, is overexpressed in human acute myeloid leukemia (AML), which is associated with poor overall survival. In both human and animal models, MORRBID hyperactivation correlates with two recurrent AML drivers, TET2 and FLT3ITD. Mice with individual mutations of Tet2 or Flt3ITD develop features of chronic myelomonocytic leukemia (CMML) and myeloproliferative neoplasm (MPN), respectively, and combined presence results in AML. We observe increased levels of Morrbid in murine models of CMML, MPN, and AML. Functionally, loss of Morrbid in these models induces increased expression of Bim and cell death in immature and mature myeloid cells, which results in reduced infiltration of leukemic cells in tissues and prolongs the survival of AML mice.

Embryonic ethanol exposure alters expression of sox2 and other early transcripts in zebrafish, producing gastrulation defects.

Ethanol exposure during prenatal development causes fetal alcohol spectrum disorder (FASD), the most frequent preventable birth defect and neurodevelopmental disability syndrome. The molecular targets of ethanol toxicity during development are poorly understood. Developmental stages surrounding gastrulation are very sensitive to ethanol exposure. To understand the effects of ethanol on early transcripts during embryogenesis, we treated zebrafish embryos with ethanol during pre-gastrulation period and examined the transcripts by Affymetrix GeneChip microarray before gastrulation. We identified 521 significantly dysregulated genes, including 61 transcription factors in ethanol-exposed embryos. Sox2, the key regulator of pluripotency and early development was significantly reduced. Functional annotation analysis showed enrichment in transcription regulation, embryonic axes patterning, and signaling pathways, including Wnt, Notch and retinoic acid. We identified all potential genomic targets of 25 dysregulated transcription factors and compared their interactions with the ethanol-dysregulated genes. This analysis predicted that Sox2 targeted a large number of ethanol-dysregulated genes. A gene regulatory network analysis showed that many of the dysregulated genes are targeted by multiple transcription factors. Injection of sox2 mRNA partially rescued ethanol-induced gene expression, epiboly and gastrulation defects. Additional studies of this ethanol dysregulated network may identify therapeutic targets that coordinately regulate early development.

SliceIt: A genome-wide resource and visualization tool to design CRISPR/Cas9 screens for editing protein-RNA interaction sites in the human genome.

Several protein-RNA cross linking protocols have been established in recent years to delineate the molecular interaction of an RNA Binding Protein (RBP) and its target RNAs. However, functional dissection of the role of the RBP binding sites in modulating the post-transcriptional fate of the target RNA remains challenging. CRISPR/Cas9 genome editing system is being commonly employed to perturb both coding and noncoding regions in the genome. With the advancements in genome-scale CRISPR/Cas9 screens, it is now possible to not only perturb specific binding sites but also probe the global impact of protein-RNA interaction sites across cell types. Here, we present SliceIt (http://sliceit.soic.iupui.edu/), a database of in silico sgRNA (single guide RNA) library to facilitate conducting such high throughput screens. SliceIt comprises of ~4.8 million unique sgRNAs with an estimated range of 2-8 sgRNAs designed per RBP binding site, for eCLIP experiments of >100 RBPs in HepG2 and K562 cell lines from the ENCODE project. SliceIt provides a user friendly environment, developed using advanced search engine framework, Elasticsearch. It is available in both table and genome browser views facilitating the easy navigation of RBP binding sites, designed sgRNAs, exon expression levels across 53 human tissues along with prevalence of SNPs and GWAS hits on binding sites. Exon expression profiles enable examination of locus specific changes proximal to the binding sites. Users can also upload custom tracks of various file formats directly onto genome browser, to navigate additional genomic features in the genome and compare with other types of omics profiles. All the binding site-centric information is dynamically accessible via "search by gene", "search by coordinates" and "search by RBP" options and readily available to download. Validation of the sgRNA library in SliceIt was performed by selecting RBP binding sites in Lipt1 gene and designing sgRNAs. Effect of CRISPR/Cas9 perturbations on the selected binding sites in HepG2 cell line, was confirmed based on altered proximal exon expression levels using qPCR, further supporting the utility of the resource to design experiments for perturbing protein-RNA interaction networks. Thus, SliceIt provides a one-stop repertoire of guide RNA library to perturb RBP binding sites, along with several layers of functional information to design both low and high throughput CRISPR/Cas9 screens, for studying the phenotypes and diseases associated with RBP binding sites.