Endogenous siRNAs promote proteostasis and longevity in germline-less Caenorhabditis elegans.

How lifespan and the rate of aging are set is a key problem in biology. Small RNAs are conserved molecules that impact diverse biological processes through the control of gene expression. However, in contrast to miRNAs, the role of endo-siRNAs in aging remains unexplored. Here, by combining deep sequencing and genomic and genetic approaches in Caenorhabditis elegans, we reveal an unprecedented role for endo-siRNA molecules in the maintenance of proteostasis and lifespan extension in germline-less animals. Furthermore, we identify an endo-siRNA-regulated tyrosine phosphatase, which limits the longevity of germline-less animals by restricting the activity of the heat shock transcription factor HSF-1. Altogether, our findings point to endo-siRNAs as a link between germline removal and the HSF-1 proteostasis and longevity-promoting somatic pathway. This establishes a role for endo siRNAs in the aging process and identifies downstream genes and physiological processes that are regulated by the endo siRNAs to affect longevity.

QsRNA-seq: A protocol for generating libraries for high-throughput sequencing of small RNAs.

Small RNAs (sRNAs) are 20-30 nt long non-coding RNA molecules that regulate essentially all cellular processes. Besides being an intensively studied topic in academic research, sRNAs also hold a promise as clinical biomarkers. While the need for expressional profiling of sRNAs is growing, preparation of sRNA libraries for high-throughput sequencing (HTS) remains technically challenging, due to their small size. The common PAGE-based protocol is time-consuming and inefficient due to material loss, while gel-free protocols generate libraries of insufficient quality. To overcome these shortcomings, we modified the conditions of size-selection by Solid Phase Reversible Immobilization (SPRI) in a way that allows separation of nucleic acids shorter than 100 nt and differing in length by only 20 nt. Implementing the method for preparation of small RNA libraries for HTS resulted in QsRNA-seq, a gel-free, fast and easy-to-perform protocol, amenable to automation, generating very clean libraries that result in high-depth expression data. The protocol also utilizes Unique Molecular Identifiers (UMI) for reduction of library preparation biases and to quantify expression levels. QsRNA-seq provides an excellent solution to the growing needs for small RNA expression profiling for research clinical use.

A-to-I RNA Editing Affects lncRNAs Expression after Heat Shock.

Adenosine to inosine (A-to-I) RNA editing is a highly conserved regulatory process carried out by adenosine-deaminases (ADARs) on double-stranded RNA (dsRNAs). Although a considerable fraction of the transcriptome is edited, the function of most editing sites is unknown. Previous studies indicate changes in A-to-I RNA editing frequencies following exposure to several stress types. However, the overall effect of stress on the expression of ADAR targets is not fully understood. Here, we performed high-throughput RNA sequencing of wild-type and ADAR mutant Caenorhabditis elegans worms after heat-shock to analyze the effect of heat-shock stress on the expression pattern of genes. We found that ADAR regulation following heat-shock does not directly involve heat-shock related genes. Our analysis also revealed that long non-coding RNAs (lncRNAs) and pseudogenes, which have a tendency for secondary RNA structures, are enriched among upregulated genes following heat-shock in ADAR mutant worms. The same group of genes is downregulated in ADAR mutant worms under permissive conditions, which is likely, considering that A-to-I editing protects endogenous dsRNA from RNA-interference (RNAi). Therefore, temperature increases may destabilize dsRNA structures and protect them from RNAi degradation, despite the lack of ADAR function. These findings shed new light on the dynamics of gene expression under heat-shock in relation to ADAR function.

QsRNA-seq: a method for high-throughput profiling and quantifying small RNAs.

The ability to profile and quantify small non-coding RNAs (sRNAs), specifically microRNAs (miRNAs), using high-throughput sequencing is challenging because of their small size. We developed QsRNA-seq, a method for preparation of sRNA libraries for high-throughput sequencing that overcomes this difficulty by enabling a gel-free separation of fragments shorter than 100 nt that differ only by 20 nt in length. The method allows the use of unique molecular identifiers for quantification and is more amenable to automation than gel-based methods. We show that QsRNA-seq gives very accurate, comprehensive, and reproducible results by looking at miRNAs in Caenorhabditis elegans embryos and larvae.

A-to-I RNA editing promotes developmental stage-specific gene and lncRNA expression.

A-to-I RNA editing is a conserved widespread phenomenon in which adenosine (A) is converted to inosine (I) by adenosine deaminases (ADARs) in double-stranded RNA regions, mainly noncoding. Mutations in ADAR enzymes in Caenorhabditis elegans cause defects in normal development but are not lethal as in human and mouse. Previous studies in C. elegans indicated competition between RNA interference (RNAi) and RNA editing mechanisms, based on the observation that worms that lack both mechanisms do not exhibit defects, in contrast to the developmental defects observed when only RNA editing is absent. To study the effects of RNA editing on gene expression and function, we established a novel screen that enabled us to identify thousands of RNA editing sites in nonrepetitive regions in the genome. These include dozens of genes that are edited at their 3' UTR region. We found that these genes are mainly germline and neuronal genes, and that they are down-regulated in the absence of ADAR enzymes. Moreover, we discovered that almost half of these genes are edited in a developmental-specific manner, indicating that RNA editing is a highly regulated process. We found that many pseudogenes and other lncRNAs are also extensively down-regulated in the absence of ADARs in the embryo but not in the fourth larval (L4) stage. This down-regulation is not observed upon additional knockout of RNAi. Furthermore, levels of siRNAs aligned to pseudogenes in ADAR mutants are enhanced. Taken together, our results suggest a role for RNA editing in normal growth and development by regulating silencing via RNAi.

Function of cancer associated genes revealed by modern univariate and multivariate association tests.

Copy number variation (CNV) plays a role in pathogenesis of many human diseases, especially cancer. Several whole genome CNV association studies have been performed for the purpose of identifying cancer associated CNVs. Here we undertook a novel approach to whole genome CNV analysis, with the goal being identification of associations between CNV of different genes (CNV-CNV) across 60 human cancer cell lines. We hypothesize that these associations point to the roles of the associated genes in cancer, and can be indicators of their position in gene networks of cancer-driving processes. Recent studies show that gene associations are often non-linear and non-monotone. In order to obtain a more complete picture of all CNV associations, we performed omnibus univariate analysis by utilizing dCov, MIC, and HHG association tests, which are capable of detecting any type of association, including non-monotone relationships. For comparison we used Spearman and Pearson association tests, which detect only linear or monotone relationships. Application of dCov, MIC and HHG tests resulted in identification of twice as many associations compared to those found by Spearman and Pearson alone. Interestingly, most of the new associations were detected by the HHG test. Next, we utilized dCov's and HHG's ability to perform multivariate analysis. We tested for association between genes of unknown function and known cancer-related pathways. Our results indicate that multivariate analysis is much more effective than univariate analysis for the purpose of ascribing biological roles to genes of unknown function. We conclude that a combination of multivariate and univariate omnibus association tests can reveal significant information about gene networks of disease-driving processes. These methods can be applied to any large gene or pathway dataset, allowing more comprehensive analysis of biological processes.

Genome-wide adaptive complexes to underground stresses in blind mole rats Spalax.

The blind mole rat (BMR), Spalax galili, is an excellent model for studying mammalian adaptation to life underground and medical applications. The BMR spends its entire life underground, protecting itself from predators and climatic fluctuations while challenging it with multiple stressors such as darkness, hypoxia, hypercapnia, energetics and high pathonecity. Here we sequence and analyse the BMR genome and transcriptome, highlighting the possible genomic adaptive responses to the underground stressors. Our results show high rates of RNA/DNA editing, reduced chromosome rearrangements, an over-representation of short interspersed elements (SINEs) probably linked to hypoxia tolerance, degeneration of vision and progression of photoperiodic perception, tolerance to hypercapnia and hypoxia and resistance to cancer. The remarkable traits of the BMR, together with its genomic and transcriptomic information, enhance our understanding of adaptation to extreme environments and will enable the utilization of BMR models for biomedical research in the fight against cancer, stroke and cardiovascular diseases.

Identification of a novel stress-responsive gene Hi95 involved in regulation of cell viability.

cDNA microarray hybridization was used in an attempt to identify novel genes participating in cellular responses to prolonged hypoxia. One of the identified novel genes, designated Hi95 shared significant homology to a p53-regulated GADD family member PA26. In addition to its induction in response to prolonged hypoxia, the increased Hi95 transcription was observed following DNA damage or oxidative stress, but not following hyperthermia or serum starvation. Whereas induction of Hi95 by prolonged hypoxia or by oxidative stress is most likely p53-independent, its induction in response to DNA damaging treatments (gamma- or UV-irradiation, or doxorubicin) occurs in a p53-dependent manner. Overexpression of Hi95 full-length cDNA was found toxic for many types of cultured cells directly leading either to their apoptotic death or to sensitization to serum starvation and DNA damaging treatments. Unexpectedly, conditional overexpression of the Hi95 cDNA in MCF7-tet-off cells resulted in their protection against cell death induced by hypoxia/glucose deprivation or H(2)O(2). Thus, Hi95 gene seems to be involved in complex regulation of cell viability in response to different stress conditions.