Search

Over the last few years, the number of microRNAs in the human genome has become a controversially debated issue. Several publications reported thousands of putative novel microRNAs not included in the curated microRNA gene database MirGeneDB and the repository miRBase. Recently, by using sequencing of ∼300 human tissues and cell lines, the human RNA atlas, an expanded inventory of human RNA annotations, was published, reporting thousands of putative microRNAs. We, the developers of established microRNA prediction tools and hosts of MirGeneDB, raise concerns about the frequently applied prediction and functional validation strategies, briefly discussing the drawbacks of false positive detections. By means of quantifying well-established biogenesis-derived features, we show that the reported novel microRNAs essentially represent false-positives and argue that the human microRNA complement, at about 550 microRNA genes, is already near complete. Output of available tools must be curated as false predictions will misguide scientists looking for biomarkers or therapeutic targets.

Malaria inflicts the highest rate of morbidity and mortality among the vector-borne diseases. The dramatic bottleneck of parasite numbers that occurs in the gut of the obligatory mosquito vector provides a promising target for novel control strategies. Using single-cell transcriptomics, we analyzed Plasmodium falciparum development in the mosquito gut, from unfertilized female gametes through the first 20 hours post blood feeding, including the zygote and ookinete stages. This study revealed the temporal gene expression of the ApiAP2 family of transcription factors, and of parasite stress genes in response to the harsh environment of the mosquito midgut. Further, employing structural protein prediction analyses we found several upregulated genes predicted to encode intrinsically disordered proteins (IDPs), a category of proteins known for their importance in regulation of transcription, translation and protein-protein interactions. IDPs are known for their antigenic properties and may serve as suitable targets for antibody or peptide-based transmission suppression strategies. In total, this study uncovers the P. falciparum transcriptome from early-to-late parasite development in the mosquito midgut, inside its natural vector, which provides an important resource for future malaria transmission-blocking initiatives. This repository contains original data such as raw sequencing data (fastq files), count matrices, metadata files, images parameters, GO term files, etc.isolated time points of P.falciparum single parasite from mosquito midgut experiments usingmicro-manipulation performed in this study.These folders contain input data for the R markdown script, which can be found in this GitHub repository:Github link. The filescontains all data necessary to reproduce the results of the study and raw_seq_filesis a directory for all of the fastq files and their metadata. processed_datafiles.zip have the expression matrix reads counts and metadata toreproduce the data of this study can also be found in the same GitHub repository (https://github.com/ANKARKLEVLAB/Single-cell-P.falciparum-midgut). for data interactive Visualizationplease go toScRNA-seq data Visualization Experiment design:Plasmodium falciparum NF54 gametocyte producing cell line were committed to sexual stages (gametocytes) and pooled to 1% gametocytaemia and 50% hematocrit. Gametocyte were membrane fed to Anopheles gambiae adult mosquito and dissected at time points; 2, 4, 8, 12 and 20 hours post infection, parasites were then isolated using Micro-maniupulation based on Anti-Pfs25 antibody conjugated to Alexa-488 flour and Hoechststain. cDNA library were prepared using SMART-seq2 and illumina sequenced. ScRNA-seq data were then analyzed using various computational tools, including Seurat toolkit and Slingshot for Pseudotime analyses. directories structure: Processed_datafiles.zip: metadata.txt (cell IDs + clusters + Time points + Pseudotime modules) zygote_raw_expressionmatrix.txt (columns = cell Ids , rows = gene accession number) img_data.txt ( imaging data ,The fluorescence intensity + Volume) metadata_NFFG.txt ( cell IDs + clusters + Time points including unfertilized female gametes) metadata_without_NFFG.txt (cell IDs + clusters + Pseudotime modules without unfertilized female gametes) raw seq data files and metadata batch_1.zip(RNA-seq fastq files linked to metadata_batch_1.csv) batch_2.zip(RNA-seq fastq files linked to metadata_batch_2.csv) raw_seq_metadata.zip folder RNA-seq metadata files information sheet: Sc_Pfalciparum_midgut_stages_Mohammed et al.xlsx (information about the study and experimental design)&nbsp

MicroRNAs are gene regulatory molecules that play important roles in numerous biological processes including human health. The function of a given microRNA is defined by its selection of target transcripts, yet current state-of-the-art experimental methods to identify microRNA targets are laborious and require millions of cells. We have overcome these limitations by fusing the microRNA effector protein Argonaute2 to the RNA editing domain of ADAR2, allowing for the first time the detection of microRNA targets transcriptome-wide in single cells. Our agoTRIBE method reports functional microRNA targets which are additionally supported by evolutionary sequence conservation. As a proof-of-principle, we study microRNA interactions in single cells, and find substantial differential targeting across the cell cycle. Lastly, agoTRIBE additionally provides transcriptome-wide measurements of RNA abundance and will allow the deconvolution of microRNA targeting in complex samples such as tissues at the single-cell level.

Respiratory syncytial virus (RSV) causes acute lower respiratory tract infection in infants, immunocompromised individuals and the elderly. As the only current specific treatment options for RSV are monoclonal antibodies, there is a need for efficacious antiviral treatments against RSV to be developed. We have previously shown that a group of synthetic non-coding single-stranded DNA oligonucleotides with lengths of 25–40 nucleotides can inhibit RSV infection in vitro and in vivo. Based on this, herein, we investigate whether naturally occurring single-stranded small non-coding RNA (sncRNA) fragments present in the airways have antiviral effects against RSV infection. From publicly available sequencing data, we selected sncRNA fragments such as YRNAs, tRNAs and rRNAs present in human bronchoalveolar lavage fluid (BALF) from healthy individuals. We utilized a GFP-expressing RSV to show that pre-treatment with the selected sncRNA fragments inhibited RSV infection in A549 cells in vitro. Furthermore, by using a flow cytometry-based binding assay, we demonstrate that these naturally occurring sncRNAs fragments inhibit viral infection most likely by binding to the RSV entry receptor nucleolin and thereby preventing the virus from binding to host cells, either directly or via steric hindrance. This finding highlights a new function of sncRNAs and displays the possibility of using naturally occurring sncRNAs as treatments against RSV.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition, or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 yr, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by (i) the presence of canid/Caniformia-specific sequences that never evolved outside of this clade, (ii) tissue-specific expression patterns (cartilage, liver, and muscle) that resemble those of modern dogs, and (iii) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition, and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing.

Accumulating behavioural data indicate that aggregation pheromones may mediate the formation and maintenance of mosquito swarms. However, chemical cues possibly luring mosquitoes to swarms have not been adequately investigated, and the likely molecular incitants of these complex reproductive behaviours remain unknown. Here we show that males of the important malaria vector species Anopheles arabiensis and An. gambiae produce and release aggregation pheromones that attract individuals to the swarm and enhance mating success. We found that males of both species released significantly higher amounts of 3-hydroxy-2-butanone (acetoin), 6-methyl-5-hepten-2-one (sulcatone), octanal, nonanal and decanal during swarming in the laboratory. Feeding males with stable-isotope-labelled glucose revealed that the males produced these five compounds. A blend composed of synthetic analogues to these swarming odours proved highly attractive to virgin males and females of both species under laboratory conditions and substantially increased mating in five African malaria vectors (An. gambiae, An. coluzzii, An. arabiensis, An. merus and An. funestus) in semi-field experiments. Our results not only narrow a conspicuous gap in understanding a vital aspect of the chemical ecology of male mosquitoes but also demonstrate fundamental roles of rhythmic and metabolic genes in the physiology and behavioural regulation of these vectors. These identified aggregation pheromones have great potential for exploitation against these highly dangerous insects. Manipulating such pheromones could increase the efficacy of malaria-vector control programmes.

Respiratory syncytial virus (RSV) causes acute lower respiratory tract infection in infants, immunocompromised individuals and the elderly. As the only current specific treatment options for RSV are monoclonal antibodies, there is a need for efficacious antiviral treatments against RSV to be developed. We have previously shown that a group of synthetic non-coding single-stranded DNA oligonucleotides with lengths of 25-40 nucleotides can inhibit RSV infection in vitro and in vivo. Based on this, herein, we investigate whether naturally occurring single-stranded small non-coding RNA (sncRNA) fragments present in the airways have antiviral effects against RSV infection. From publicly available sequencing data, we selected sncRNA fragments such as YRNAs, tRNAs and rRNAs present in human bronchoalveolar lavage fluid (BALF) from healthy individuals. We utilized a GFP-expressing RSV to show that pre-treatment with the selected sncRNA fragments inhibited RSV infection in A549 cells in vitro. Furthermore, by using a flow cytometry-based binding assay, we demonstrate that these naturally occurring sncRNAs fragments inhibit viral infection most likely by binding to the RSV entry receptor nucleolin and thereby preventing the virus from binding to host cells, either directly or via steric hindrance. This finding highlights a new function of sncRNAs and displays the possibility of using naturally occurring sncRNAs as treatments against RSV.

miRNAs are short noncoding RNA molecules that regulate gene expression by inhibiting translation or inducing degradation of target mRNAs. miRNAs are often expressed as polycistronic transcripts, so-called miRNA clusters, containing several miRNA precursors. The largest mammalian miRNA cluster, the miR-379-410 cluster, is expressed primarily during embryonic development and in the adult brain; however, downstream regulation of this cluster is not well understood. Here, we investigated adenosine deamination to inosine (RNA editing) in the miR-379-410 cluster by adenosine deaminase acting on RNA (ADAR) enzymes as a possible mechanism modulating the expression and activity of these miRNAs in a brain-specific manner. We show that the levels of editing in the majority of mature miRNAs are lower than the editing levels of the corresponding site in primary miRNA precursors. However, for one miRNA, miR-376b-3p, editing was significantly higher in the mature form than in the primary precursor. We found miR-376b-3p maturation is negatively regulated by ADAR2 in an editing activity-independent manner, whereas ADAR1-mediated and ADAR2-mediated editing were observed to be competitive. In addition, the edited miR-376b-3p targets a different set of mRNAs than unedited miR-376b-3p, including 4-aminobutyrate aminotransferase, encoding the enzyme responsible for the catabolism of the neurotransmitter gamma aminobutyric acid (GABA). Expression of edited miR-376b-3p led to increased intracellular GABA levels as well as increased cell surface presentation of GABA type A receptors. Our results indicate that both editing and editing-independent effects modulate the expression of miR-376b-3p, with the potential to regulate GABAergic signaling in the brain.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition, or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 yr, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by (i) the presence of canid/Caniformia-specific sequences that never evolved outside of this clade, (ii) tissue-specific expression patterns (cartilage, liver, and muscle) that resemble those of modern dogs, and (iii) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition, and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing.

SummaryIn early development, the environment triggers mnemonic epigenomic programs resulting in memory and learning experiences to confer cognitive phenotypes into adulthood. To uncover how environmental stimulation impacts the epigenome and genome organization, we used the paradigm of environmental enrichment (EE) in young mice constantly receiving novel stimulation. We profiled epigenome and chromatin architecture in whole cortex and sorted neurons by deep-sequencing techniques. Specifically, we studied chromatin accessibility, gene and protein regulation, and 3D genome conformation, combined with predicted enhancer and chromatin interactions. We identified increased chromatin accessibility, transcription factor binding including CTCF-mediated insulation, differential occupancy of H3K36me3 and H3K79me2, and changes in transcriptional programs required for neuronal development. EE stimuli led to local genome re-organization by inducing increased contacts between chromosomes 7 and 17 (inter-chromosomal). Our findings support the notion that EE-induced learning and memory processes are directly associated with the epigenome and genome organization.Highlights-Environmental enrichment (EE) alters chromatin conformation, CTCF binding, and spatially 3D genome changes, thereby regulating cognitive function during the first steps of life after birth.-Transcription-associated gene body marks H3K79me2 and H3K36me3 are differently influenced by EE in cortical brain cells and binding is exacerbated upon stimulation in an age-dependent manner.-EE-induced changes of 3D genome organization increase inter-chromosomal interactions of genes associated with synaptic transmission and AMPA receptor genes on chromosomes 7 and 17.

Objective The liver performs a central role in regulating energy homeostasis by increasing glucose output during fasting. Recent studies on Argonaute2 (Ago2), a key RNA-binding protein mediating the microRNA pathway, have illustrated its role in adaptive mechanisms according to changes in metabolic demand. Here we sought to characterize the functional role of Ago2 in the liver in the maintenance of systemic glucose homeostasis. Methods We first analyzed Ago2 expression in mouse primary hepatocyte cultures after modulating extracellular glucose concentrations and in the presence of activators or inhibitors of glucokinase activity. We then characterized a conditional loss-of-function mouse model of Ago2 in liver for alterations in systemic energy metabolism. Results Here we show that Ago2 expression in liver is directly correlated to extracellular glucose concentrations and that modulating glucokinase activity is adequate to affect hepatic Ago2 levels. Conditional deletion of Ago2 in liver resulted in decreased fasting glucose levels in addition to reducing hepatic glucose production. Moreover, loss of Ago2 promoted hepatic expression of AMP-activated protein kinase α1 (AMPKα1) by de-repressing its targeting by miR-148a, an abundant microRNA in the liver. Deletion of Ago2 from hyperglycemic, obese, and insulin-resistant Lepob/ob mice reduced both random and fasted blood glucose levels and body weight and improved insulin sensitivity. Conclusions These data illustrate a central role for Ago2 in the adaptive response of the liver to fasting. Ago2 mediates the suppression of AMPKα1 by miR-148a, thereby identifying a regulatory link between non-coding RNAs and a key stress regulator in the hepatocyte., Highlights • Hepatic Ago2 levels correlate with changes in extracellular glucose concentrations. • Conditional deletion of Ago2 in liver decreased fasting glucose levels. • Loss of Ago2 promoted AMPKα1 by de-repressing its targeting by miR-148a. • Deletion of Ago2 from Lepob/ob mice improved glycemia and insulin sensitivity.

The molecular basis of many health outcomes attributed to solar ultraviolet radiation (UVR) is unknown. We tested the hypothesis that they may originate from transcriptional changes in blood cells. This was determined by assessing the effect of fluorescent solar simulated radiation (FSSR) on the transcriptional profile of peripheral blood pre- and 6 h, 24 h and 48 h post-exposure in nine healthy volunteers. Expression of 20 genes was down-regulated and one was up-regulated at 6 h after FSSR. All recovered to baseline expression at 24 h or 48 h. These genes have been associated with immune regulation, cancer and blood pressure; health effects attributed to vitamin D via solar UVR exposure. Plasma 25-hydroxyvitamin D3 [25OHD3] levels increased over time after FSSR and were maximal at 48 h. The increase was more pronounced in participants with low basal 25OHD3 levels. Mediation analyses suggested that changes in gene expression due to FSSR were independent of 25OHD3 and blood cell subpopulations.

We have identified and validated a spaceflight-associated microRNA (miRNA) signature that is shared by rodents and humans in response to simulated and real spaceflight and regulates vascular damage by simulated deep space radiation. In previous studies, we identified miRNAs that are predicted to regulate rodent responses to spaceflight in low-Earth orbit. Here we have validated the expression of key spaceflight-associated miRNAs in rodents reacting to simulated spaceflight conditions (exposure to ionizing radiation combined with simulated microgravity) and in astronaut samples from the NASA Twin Study. To demonstrate their physiological relevance, we have utilized antagomirs to target a subset of miRNAs from the spaceflight associated signature to successfully rescue simulated deep space radiation-mediated damage in human 3D vascular constructs.

Single-cell RNA sequencing studies on gene co-expression patterns could yield important regulatory and functional insights, but have so far been limited by the confounding effects of differentiation and cell cycle. We apply a tailored experimental design that eliminates these confounders, and report thousands of intrinsically covarying gene pairs in mouse embryonic stem cells. These covariations form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariations and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariations. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex formation. Our results lend support to the concept of post-transcriptional RNA operons, but we further present evidence that nuclear proximity of genes may provide substantial functional regulation in mammalian single cells., Gene expression covariation can be studied by single-cell RNA sequencing. Here the authors analyze intrinsically covarying gene pairs by eliminating the confounding effects in single-cell experiments and observe covariation of proximal genes and miRNA-induced covariation of target mRNAs.

Ancient DNA sequencing is the key technology for paleogenomic studies and today a routine method in many laboratories. Recent analyses have shown that, under favoring conditions, also RNA can be sequenced from historical and even ancient samples. We have re-analyzed ancient RNA data from a Pleistocene canid and find - in addition to the previously described messenger RNA fragments - intact microRNAs, which are short transcripts with important gene regulatory functions. With an extraordinary age of 14,300 years, the canid microRNA profiles are the oldest ever reported. Despite their age, we show that the microRNA profiles are conclusive of taxonomic origin, tissue identity with organ- and cell-type specific signatures, and that they yield glimpses into gene regulatory activity and biological processes from the Pleistocene. In summary, we here show that straightforward microRNA analyses hold great promise for deeper insights into gene regulation in extinct animals.

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes. The function of the SWI/SNF complexes in transcriptional initiation has been well studied, while a function in alternative splicing has only been studied for a few cases for BRM-containing SWI/SNF complexes. Here, we have expressed BRG1 in C33A cells, a BRG1 and BRM-deficient cell line, and we have analysed the effects on the transcriptome by RNA sequencing. We have shown that BRG1 expression affects the splicing of a subset of genes. For some, BRG1 expression favours exon inclusion and for others, exon skipping. Some of the changes in alternative splicing induced by BRG1 expression do not require the ATPase activity of BRG1. Among the exons regulated through an ATPase-independent mechanism, the included exons had signatures of high GC-content and lacked a positioned nucleosome at the exon. By investigating three genes in which the expression of either wild-type BRG1 or a BRG1-ATPase-deficient variant favoured exon inclusion, we showed that expression of the ATPases promotes the local recruitment of RNA binding factors to chromatin and RNA in a differential manner. The hnRNPL, hnRNPU and SAM68 proteins associated to chromatin in C33A cells expressing BRG1 or BRM, but their association with RNA varied. We propose that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and altering their binding to the nascent pre-mRNA, which changes RNP structure.Author summarySplicing, in particular alternative splicing, is a combinatorial process which involves splicing factor complexes and many RNA binding splicing regulatory proteins in different constellations. Most splicing events occur during transcription, which also makes the DNA sequence, the chromatin state and the transcription rate at the exons important components that influence the splicing outcome. We show here that the ATP-dependent chromatin remodelling complex SWI/SNF influences the interactions of splicing regulatory factors with RNA during transcription on certain exons that have a high GC-content. The splicing on this type of exon rely on the ATPase BRG1 and favour inclusion of alternative exons in an ATP-independent manner. SWI/SNF complexes are known to alter the chromatin structure at promoters in transcription initiation, and have been previously shown to alter the transcription rate or nucleosome position in splicing. Our results suggests a further mechanism for chromatin remodelling proteins in splicing: to change the interaction patterns of RNA binding splicing regulatory factors at alternative exons to alter the splicing outcome.

Single-cell RNA sequencing studies into gene co-expression patterns could yield important new regulatory and functional insights, but have so far been limited by the confounding effects of cell differentiation and the cell cycle. We apply a tailored experimental design that eliminates these confounders, and report >80,000 intrinsically covarying gene pairs in mouse embryonic stem cells. These covariances form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariances, and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariances. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex stoichiometry. Our results lend support to the concept of post-transcriptional ‘RNA operons’, but we further present evidence that nuclear proximity of genes on the same or even distinct chromosomes also provides substantial functional regulation in mammalian single cells.

Summary Background Terrestrial ultraviolet (UV) radiation causes erythema, oxidative stress, DNA mutations and skin cancer. Skin can adapt to these adverse effects by DNA repair, apoptosis, keratinization and tanning. Objectives To investigate the transcriptional response to fluorescent solar‐simulated radiation (FSSR) in sun‐sensitive human skin in vivo. Methods Seven healthy male volunteers were exposed to 0, 3 and 6 standard erythemal doses (SED). Skin biopsies were taken at 6 h and 24 h after exposure. Gene and microRNA expression were quantified with next generation sequencing. A set of candidate genes was validated by quantitative polymerase chain reaction (qPCR); and wavelength dependence was examined in other volunteers through microarrays. Results The number of differentially expressed genes increased with FSSR dose and decreased between 6 and 24 h. Six hours after 6 SED, 4071 genes were differentially expressed, but only 16 genes were affected at 24 h after 3 SED. Genes for apoptosis and keratinization were prominent at 6 h, whereas inflammation and immunoregulation genes were predominant at 24 h. Validation by qPCR confirmed the altered expression of nine genes detected under all conditions; genes related to DNA repair and apoptosis; immunity and inflammation; pigmentation; and vitamin D synthesis. In general, candidate genes also responded to UVA1 (340–400 nm) and/or UVB (300 nm), but with variations in wavelength dependence and peak expression time. Only four microRNAs were differentially expressed by FSSR. Conclusions The UV radiation doses of this acute study are readily achieved daily during holidays in the sun, suggesting that the skin transcriptional profile of ‘typical’ holiday makers is markedly deregulated. What's already known about this topic? The skin's transcriptional profile underpins its adverse (i.e. inflammation) and adaptive molecular, cellular and clinical responses (i.e. tanning, hyperkeratosis) to solar ultraviolet radiation.Few studies have assessed microRNA and gene expression in vivo in humans, and there is a lack of information on dose, time and waveband effects. What does this study add? Acute doses of fluorescent solar‐simulated radiation (FSSR), of similar magnitude to those received daily in holiday situations, markedly altered the skin's transcriptional profiles.The number of differentially expressed genes was FSSR‐dose‐dependent, reached a peak at 6 h and returned to baseline at 24 h.The initial transcriptional response involved apoptosis and keratinization, followed by inflammation and immune modulation. In these conditions, microRNA expression was less affected than gene expression., Linked Comment: Hart. Br J Dermatol 2020; 182:1328–1329. Plain language summary available online Respond to this article

Combined measurements of mRNA and protein expression in single cells enables in-depth analysis of cellular states. We present single-cell protein and RNA co-profiling (SPARC), an approach to simultaneously measure global mRNA and large sets of intracellular protein in individual cells. Using SPARC, we show that mRNA expression fails to accurately reflect protein abundance at the time of measurement in human embryonic stem cells, although the direction of changes of mRNA and protein expression are in agreement during cellular differentiation. Moreover, protein levels of transcription factors better predict their downstream effects than do the corresponding transcripts. We further show that changes of the balance between protein and mRNA expression levels can be applied to infer expression kinetic trajectories, revealing future states of individual cells. Finally, we highlight that mRNA expression may be more varied among cells than levels of the corresponding proteins. Overall, our results demonstrate that mRNA and protein measurements in single cells provide different and complementary information regarding cell states. Accordingly, SPARC can offer valuable insights in gene expression programs of single cells.

Yeast Hsp70 nucleotide exchange factor Fes1 is expressed by rare alternative splicing as two isoforms. Fes1L is targeted to the nucleus, and Fes1S localizes to the cytosol and is required for the efficient proteasomal degradation of cytosolic misfolded proteins, as well as of species that are imported into the nucleus for degradation., Cells maintain proteostasis by selectively recognizing and targeting misfolded proteins for degradation. In Saccharomyces cerevisiae, the Hsp70 nucleotide exchange factor Fes1 is essential for the degradation of chaperone-associated misfolded proteins by the ubiquitin-proteasome system. Here we show that the FES1 transcript undergoes unique 3′ alternative splicing that results in two equally active isoforms with alternative C-termini, Fes1L and Fes1S. Fes1L is actively targeted to the nucleus and represents the first identified nuclear Hsp70 nucleotide exchange factor. In contrast, Fes1S localizes to the cytosol and is essential to maintain proteostasis. In the absence of Fes1S, the heat-shock response is constitutively induced at normally nonstressful conditions. Moreover, cells display severe growth defects when elevated temperatures, amino acid analogues, or the ectopic expression of misfolded proteins, induce protein misfolding. Importantly, misfolded proteins are not targeted for degradation by the ubiquitin-proteasome system. These observations support the notion that cytosolic Fes1S maintains proteostasis by supporting the removal of toxic misfolded proteins by proteasomal degradation. This study provides key findings for the understanding of the organization of protein quality control mechanisms in the cytosol and nucleus.

Background MicroRNAs (miRNAs) have been identified as key players in posttranscriptional gene regulation and have a significant impact on basal cell carcinoma (BCC) development. The Sonic hedgehog pathway inhibitor vismodegib has been approved for oral therapy of metastatic or advanced BCC. Here, a high-throughput miRNA sequencing analysis was carried out to identify differentially expressed miRNAs and possible novel miRNA candidates in vismodegib-treated BCC tissue. Additionally, we described our surgical experience with neoadjuvant oral vismodegib therapy. Patients and methods A punch biopsy (4 mm) from a patient with an extensive cranial BCC under oral vismodegib therapy and a corresponding nonlesional epithelial skin biopsy were harvested. Total RNA was isolated, after which a sequencing cDNA library was prepared, and cluster generation was carried out, which was followed by an ultra-high-throughput miRNA sequencing analysis to indicate the read number of miRNA expression based on miRBase 21. In addition to the identification of differentially expressed miRNAs from RNA sequencing data, additional novel miRNA candidates were determined with a tool for identifying new miRNA sequences (miRDeep2). Results We identified 33 up-regulated miRNAs (fold change ≥2) and 39 potentially new miRNA candidates (miRDeep scores 0–43.6). A manual sequence analysis of the miRNA candidates on the genomic locus of chromosome 1 with provisional IDs of chr1_1913 and chr1_421 was further carried out and rated as promising (chr1_1913) and borderline (chr1_421). Histopathology revealed skip lesions in clinically healthy appearing skin at the tumor margins, which were the cause of seven re-excisions by micrographic controlled surgery to achieve tumor-free margins. Conclusion miRNA sequencing revealed novel miRNA candidates that need to be further confirmed in functional Dicer knockout studies. Clinically, on the basis of our surgical experience described here, neoadjuvant vismodegib therapy in BCC appears to impede histopathologic evaluations with effects on surgical therapy. Thus, larger studies are necessary, but are not preferable at this time if other options are available.

SUMMARY We have identified and validated a spaceflight-associated microRNA (miRNA) signature that is shared by rodents and humans in response to simulated, short-duration and long-duration spaceflight. Previous studies have identified miRNAs that regulate rodent responses to spaceflight in low-Earth orbit, and we have confirmed the expression of these proposed spaceflight-associated miRNAs in rodents reacting to simulated spaceflight conditions. Moreover, astronaut samples from the NASA Twins Study confirmed these expression signatures in miRNA sequencing, single-cell RNA sequencing (scRNA-seq), and single-cell assay for transposase accessible chromatin (scATAC-seq) data. Additionally, a subset of these miRNAs (miR-125, miR-16, and let-7a) was found to regulate vascular damage caused by simulated deep space radiation. To demonstrate the physiological relevance of key spaceflight-associated miRNAs, we utilized antagomirs to inhibit their expression and successfully rescue simulated deep-space-radiation-mediated damage in human 3D vascular constructs., Graphical Abstract, In Brief Malkani et al. uncover the role of circulating microRNAs as both a potential biomarker for health risks associated with spaceflight and a countermeasure to mitigate the damage caused to the body by the space environment.

The circulating transcriptome is a valuable source of cancer biomarkers, which, with the exception of microRNAs (miRNAs), remains relatively unexplored. To elucidate which RNAs are present in plasma from melanoma patients and which could be used to distinguish cancer patients from healthy individuals, we used next generation sequencing (NGS), and validation was carried out by qPCR and/or ddPCR. We identified 442 different microRNAs in samples, eleven of which were differentially expressed (p &lt, 0.05). Levels of miR-134-5p and miR-320a-3p were significantly down-regulated (p &lt, 0.001) in melanoma samples (n = 96) compared to healthy controls (n = 28). Differentially expressed protein-encoding mRNA 5&prime, fragments were enriched for the angiopoietin, p21-activated kinase (PAK), and EIF2 pathways. Levels of ATM1, AMFR, SOS1, and CD109 gene fragments were up-regulated (p &lt, 0.001) in melanoma samples (n = 144) compared to healthy controls (n = 41) (AUC = 0.825). Over 40% of mapped reads were YRNAs, a class of non-coding RNAs that to date has been little explored. Expression levels of RNY3P1, RNY4P1, and RNY4P25 were significantly higher in patients with stage 0 disease than either healthy controls or more advanced stage disease (p &lt, 0.001). In conclusion, we have identified a number of novel RNA biomarkers, which, most importantly, we validated in multi-center retrospective and prospective cohorts, suggesting potential diagnostic use of these RNA species.

Motivation MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods. Results To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification. Availability and implementation https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop. Supplementary information Supplementary data are available at Bioinformatics online.

Recent studies suggest that transcription takes place at DNA double-strand breaks (DSBs), that transcripts at DSBs are processed by Drosha and Dicer into damage-induced small RNAs (diRNAs), and that diRNAs are required for DNA repair. However, diRNAs have been mostly detected in reporter constructs or repetitive sequences, and their existence at endogenous loci has been questioned by recent reports. Using the homing endonuclease I-PpoI, we have investigated diRNA production in genetically unperturbed human and mouse cells. I-PpoI is an ideal tool to clarify the requirements for diRNA production because it induces DSBs in different types of loci: the repetitive 28S locus, unique genes and intergenic loci. We show by extensive sequencing that the rDNA locus produces substantial levels of diRNAs, whereas unique genic and intergenic loci do not. Further characterization of diRNAs emerging from the 28S locus reveals the existence of two diRNA subtypes. Surprisingly, Drosha and its partner DGCR8 are dispensable for diRNA production and only one diRNAs subtype depends on Dicer processing. Furthermore, we provide evidence that diRNAs are incorporated into Argonaute. Our findings provide direct evidence for diRNA production at endogenous loci in mammalian cells and give insights into RNA processing at DSBs.

Background To date, numerous nucleic acid species have been detected in the systemic circulation including microRNAs (miRNAs); however, their functional role in this compartment remains unclear. Objective The aim of this study was to determine whether systemic levels of miRNAs abundant in blood, including the neuroendocrine tissue-enriched miR-375, are altered in response to a glucose challenge. Design Twelve healthy males were recruited for an acute crossover study that consisted of two tests each following an 8-hour fasting period. An oral glucose tolerance test (OGTT) was performed, and blood samples were collected over a 3-hour period. Following a period of at least 1 week, the same participants were administered an isoglycemic intravenous glucose infusion (IIGI) with the same blood-collection protocol. Results The glucose response curve following the IIGI mimicked that obtained after the OGTT, but as expected, systemic insulin levels were lower during the IIGI compared with the OGTT (P < 0.05). miR-375 levels in circulation were increased only in response to an OGTT and not during an IIGI. In addition, the response to the OGTT also coincided with the transient increase of circulating glucagon-like peptide (GLP)-1, GLP-2, and glucose-dependent insulinotropic polypeptide. Conclusions The present findings show levels of miR-375 increase following administration of an OGTT and, in light of its enrichment in cells of the gut, suggest that the gastrointestinal tract may play an important role in the abundance and function of this miRNA in the blood.

The molecular basis of many health outcomes attributed to solar ultraviolet radiation (UVR) is unknown. We tested the hypothesis that they may originate from transcriptional changes in blood cells. This was determined by assessing the effect of fluorescent solar simulated radiation (FSSR) on the transcriptional profile of peripheral blood pre- and 6h, 24h and 48h post-exposure in nine healthy volunteers. Expression of 20 genes was down-regulated and one was up-regulated at 6h after FSSR. All recovered to baseline expression at 24h or 48h. These genes have been associated with immune regulation, cancer and blood pressure; health effects attributed to vitamin D via solar UVR exposure. Plasma 25-hydroxyvitamin D

Tumors are composed of an evolving population of cells subjected to tissue-specific selection, which fuels tumor heterogeneity and ultimately complicates cancer driver gene identification. Here, we integrate cancer cell fraction, population recurrence, and functional impact of somatic mutations as signatures of selection into a Bayesian inference model for driver prediction. In an in-depth benchmark, we demonstrate that our model, cDriver, outperforms competing methods when analyzing solid tumors, hematological malignancies, and pan-cancer datasets. Applying cDriver to exome sequencing data of 21 cancer types from 6,870 individuals revealed 98 unreported tumor type-driver gene connections. These novel connections are highly enriched for chromatin-modifying proteins, hinting at a universal role of chromatin regulation in cancer etiology. Although infrequently mutated as single genes, we show that chromatin modifiers are altered in a large fraction of cancer patients. In summary, we demonstrate that integration of evolutionary signatures is key for identifying mutational driver genes, thereby facilitating the discovery of novel therapeutic targets for cancer treatment.

Tumors are composed of an evolving population of cells subjected to tissue-specific selection, which fuels tumor heterogeneity and ultimately complicates cancer driver gene identification. Here, we integrate cancer cell fraction, population recurrence, and functional impact of somatic mutations as signatures of selection into a Bayesian model for driver prediction. We demonstrate that our model, cDriver, outperforms competing methods when analyzing solid tumors, hematological malignancies, and pan-cancer datasets. Applying cDriver to exome sequencing data of 21 cancer types from 6,870 individuals revealed 98 unreported tumor type-driver gene connections. These novel connections are highly enriched for chromatin-modifying proteins, hinting at a universal role of chromatin regulation in cancer etiology. Although infrequently mutated as single genes, we show that chromatin modifiers are altered in a large fraction of cancer patients. In summary, we demonstrate that integration of evolutionary signatures is key for identifying mutational driver genes, thereby facilitating the discovery of novel therapeutic targets for cancer treatment. We acknowledge support of the Spanish Ministry of Economy and Competitiveness, 'Centro de Excelencia Severo Ochoa 2013-2017'. We acknowledge the support of the CERCA Programme/Generalitat de Catalunya. This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 635290. Luis Zapata has been supported by the International PhD scholarship program of La Caixa at CRG.

miRNAs are small 22-nucleotide RNAs that can post-transcriptionally regulate gene expression. It has been proposed that dietary plant miRNAs can enter the human bloodstream and regulate host transcripts; however, these findings have been widely disputed. We here conduct the first comprehensive meta-study in the field, surveying the presence and abundances of cross-species miRNAs (xenomiRs) in 824 sequencing data sets from various human tissues and body fluids. We find that xenomiRs are commonly present in tissues (17%) and body fluids (69%); however, the abundances are low, comprising 0.001% of host human miRNA counts. Further, we do not detect a significant enrichment of xenomiRs in sequencing data originating from tissues and body fluids that are exposed to dietary intake (such as liver). Likewise, there is no significant depletion of xenomiRs in tissues and body fluids that are relatively separated from the main bloodstream (such as brain and cerebro-spinal fluids). Interestingly, the majority (81%) of body fluid xenomiRs stem from rodents, which are a rare human dietary contribution but common laboratory animals. Body fluid samples from the same studies tend to group together when clustered by xenomiR compositions, suggesting technical batch effects. Last, we performed carefully designed and controlled animal feeding studies, in which we detected no transfer of plant miRNAs into rat blood, or bovine milk sequences into piglet blood. In summary, our comprehensive computational and experimental results indicate that xenomiRs originate from technical artifacts rather than dietary intake.

Profilin is vital for actin organisation in eukaryotic cells. It controls actin filament formation by binding monomeric actin and numerous proteins involved in polarised actin assembly. Important for the latter is the interaction surface formed by the N- and C-terminal helices, which pack close to each other on one side of the molecule at a distance from the actin site and mediate binding to poly-proline sequences present in many of the targeted proteins. Via these interactions, profilin contributes to the spatiotemporal control of actin filament growth. Studies of profilin dynamics in living cells by imaging techniques have been hampered by problems to generate fusion constructs with fluorophore proteins without negatively impacting on its poly-proline binding. With the object to circumvent this problem, we have generated an internal fusion of profilin with the green fluorescent variant citrine, here referred to as citrine-profilin. The characterisation of citrine-profilin (CIT-Pfn) demonstrates that it has full capacity to interact with poly-proline and also binds phosphatidylinositol lipids and actin, albeit with 10 times reduced affinity for the latter. Imaging of living cells expressing CIT-Pfn showed a distribution of the fusion protein similar to endogenous profilin. Furthermore, CIT-Pfn rescued the phenotypes observed after the Crispr/Cas9 knockout of the profilin 1 gene, including the lost migratory capacity characterising the knockout cells. Based on this, we conclude that the CIT-Pfn construct will be useful as a tool for displaying profilin localisation in living cells and obtaining information on its dynamic organisation under different conditions and activations of the actin microfilament and microtubule systems.

In response to fasting or hyperglycemia, the pancreatic β-cell alters its output of secreted insulin; however the pathways governing this adaptive response are not entirely established. While the precise role of microRNAs (miRNAs) is also unclear, a recurring theme emphasizes their function in cellular stress responses. We recently showed that miR-184, an abundant miRNA in the β-cell, regulates compensatory proliferation and secretion during insulin resistance. Consistent with previous studies showing miR-184 suppresses insulin release, expression of this miRNA was increased in islets after fasting, demonstrating an active role in the β-cell as glucose levels lower and the insulin demand ceases. Additionally, miR-184 was negatively regulated upon administration of a sucrose-rich diet in Drosophila demonstrating strong conservation of this pathway through evolution. Furthermore, miR-184 and its target Argonaute2 (Ago2) remained inversely correlated as concentrations of extracellular glucose increased, underlining a functional relationship between this miRNA and its targets. Lastly, restoration of Ago2 in the presence of miR-184 rescued suppression of miR-375-targeted genes suggesting these genes act in a coordinated manner during changes in the metabolic context. Together, these results highlight the adaptive role of miR-184 according to glucose metabolism and suggest the regulatory role of this miRNA in energy homeostasis is highly conserved.

Next-generation sequencing now for the first time allows researchers to gage the depth and variation of entire transcriptomes. However, now as rare transcripts can be detected that are present in cells at single copies, more advanced computational tools are needed to accurately annotate and profile them. microRNAs (miRNAs) are 22 nucleotide small RNAs (sRNAs) that post-transcriptionally reduce the output of protein coding genes. They have established roles in numerous biological processes, including cancers and other diseases. During miRNA biogenesis, the sRNAs are sequentially cleaved from precursor molecules that have a characteristic hairpin RNA structure. The vast majority of new miRNA genes that are discovered are mined from small RNA sequencing (sRNA-seq), which can detect more than a billion RNAs in a single run. However, given that many of the detected RNAs are degradation products from all types of transcripts, the accurate identification of miRNAs remain a non-trivial computational problem. Here, we review the tools available to predict animal miRNAs from sRNA sequencing data. We present tools for generalist and specialist use cases, including prediction from massively pooled data or in species without reference genome. We also present wet-lab methods used to validate predicted miRNAs, and approaches to computationally benchmark prediction accuracy. For each tool, we reference validation experiments and benchmarking efforts. Last, we discuss the future of the field.