MicroRNA-mediated regulation of lipid metabolism in virus-infected Emiliania huxleyi

The interactions between Emiliania huxleyi and E. huxleyi virus (EhV) regulate marine carbon and sulfur biogeochemical cycles and play a prominent role in global climate change. As a large DNA virus, EhV has developed a novel "virocell metabolism" model to meet its high metabolic needs. Although it has been widely demonstrated that EhV infection can profoundly rewire lipid metabolism, the epigenetic regulatory mechanisms of lipid metabolism are still obscure. MicroRNAs (miRNAs) can regulate biological pathways by targeting hub genes in the metabolic processes. In this study, the transcriptome, lipidome, and miRNAome were applied to investigate the epigenetic regulation of lipid metabolism in E. huxleyi cells during a detailed time course of viral infection. Combined transcriptomic, lipidomic, and physiological experiments revealed reprogrammed lipid metabolism, along with mitochondrial dysfunction and calcium influx through the cell membrane. A total of 69 host miRNAs (including 1 known miRNA) and 7 viral miRNAs were identified, 27 of which were differentially expressed. Bioinformatic prediction revealed that miRNAs involved in the regulation of lipid metabolism and a dual-luciferase reporter assay suggested that phosphatidylinositol 3-kinase (PI3K) gene might be a target of ehx-miR5. Further qPCR and western blot analysis showed a significant negative correlation between the expression of ehx-miR5 and its target gene PI3K, along with the lower activity of its downstream components (p-Akt, p-TOR, SREBP), indicating that lipid metabolism might be regulated by ehx-miR5 through the PI3K-Akt-TOR signaling pathway. Our findings reveal several novel mechanisms of viral strategies to manipulate host lipid metabolism and provide evidence that ehx-miR5 negatively modulates the expression of PI3K and disturbs lipid metabolism in the interactions between E. huxleyi and EhV.

Contrasting sea ice conditions shape microbial food webs in Hudson Bay (Canadian Arctic)

The transition from ice-covered to open water is a recurring feature of the Arctic and sub-Arctic, but microbial diversity and cascading effects on the microbial food webs is poorly known. Here, we investigated microbial eukaryote, bacterial and archaeal communities in Hudson Bay (sub-Arctic, Canada) under sea-ice cover and open waters conditions. Co-occurrence networks revealed a <3 µm pico‒phytoplankton-based food web under the ice and a >3 µm nano‒microphytoplankton-based food web in the open waters. The ice-edge communities were characteristic of post-bloom conditions with high proportions of the picophytoplankton Micromonas and Bathycoccus. Nano‒ to micro‒phytoplankton and ice associated diatoms were detected throughout the water column, with the sympagic Melosira arctica exclusive to ice-covered central Hudson Bay and Thalassiosira in open northwestern Hudson Bay. Heterotrophic microbial eukaryotes and prokaryotes also differed by ice-state, suggesting a linkage between microbes at depth and surface phytoplankton bloom state. The findings suggest that a longer open water season may favor the establishment of a large phytoplankton-based food web at the subsurface chlorophyll maxima (SCM), increasing carbon export from pelagic diatoms to deeper waters and affect higher trophic levels in the deep Hudson Bay.

A spontaneous thermo-sensitive female sterility mutation in rice enables fully mechanized hybrid breeding

Male sterility enables hybrid crop breeding to increase yields and has been extensively studied. But thermo-sensitive female sterility, which is an ideal property that may enable full mechanization in hybrid rice breeding, has rarely been investigated due to the absence of such germplasm. Here we identify the spontaneous thermo-sensitive female sterility 1 (tfs1) mutation that confers complete sterility under regular/high temperature and partial fertility under low temperature as a point mutation in ARGONAUTE7 (AGO7). AGO7 associates with miR390 to form an RNA-Induced Silencing Complex (RISC), which triggers the biogenesis of small interfering RNAs (siRNAs) from TRANS-ACTING3 (TAS3) loci by recruiting SUPPRESSOR OF GENE SILENCING (SGS3) and RNA-DEPENDENT RNA POLYMERASE6 (RDR6) to TAS3 transcripts. These siRNAs are known as tasiR-ARFs as they act in trans to repress auxin response factor genes. The mutant TFS1 (mTFS1) protein is compromised in its ability to load the miR390/miR390* duplex and eject miR390* during RISC formation. Furthermore, tasiR-ARF levels are reduced in tfs1 due to the deficiency in RDR6 but not SGS3 recruitment by mTFS1 RISC under regular/high temperature, while low temperature partially restores mTFS1 function in RDR6 recruitment and tasiR-ARF biogenesis. A miR390 mutant also exhibits female sterility, suggesting that female fertility is controlled by the miR390-AGO7 module. Notably, the tfs1 allele introduced into various elite rice cultivars endows thermo-sensitive female sterility. Moreover, field trials confirm the utility of tfs1 as a restorer line in fully mechanized hybrid rice breeding.

Exosomal miR-155-5p derived from glioma stem-like cells promotes mesenchymal transition via targeting ACOT12

Tumor-associated exosomes play essential roles in intercellular communication and the foundation of cancer microenvironment in glioma. Many mRNAs, microRNAs (miRNAs) and proteins contained in tumor-associated exosomes can be transferred to recipient cells and contribute to the progression of tumor. Nevertheless, the cellular communication between malignant cells with different heterogeneities or characteristics and resultant tumor progression are still unclear in glioma. Here, we show that exosomes released from glioma stem-like cells (GSCs) contain a significant increasing level of miR-155-5p and could be horizontally transferred to surrounding glioma cells. High expression of miR-155-5p in plasma exosomes from patients was associated with glioma diagnosis and grading. Mechanically, we found that miR-155-5p markedly reduced the expression of acetyl-CoA thioesterase 12 (ACOT12), which played as a tumor suppressor in glioma. Furthermore, mesenchymal transition was significantly promoted in glioma cells treated with GSCs-derived exosomes. In conclusion, GSCs-derived exosomal miR-155-5p play a critical role in glioma progression and facilitating tumor aggressive growth by targeting ACOT12 and promoting mesenchymal transition. Exosomal miR-155-5p is also a potential predictive biomarker for glioma, which may provoke the development of novel diagnostic and therapeutic strategies against glioma.

CircSCAF8 promotes growth and metastasis of prostate cancer through the circSCAF8-miR-140-3p/miR-335-LIF pathway

Circular RNAs (circRNAs) have been increasingly linked to cancer progression. However, the detailed biological functions of circRNAs in prostate cancer (PCa) remain unclear. Using high-throughput circRNA sequencing, we previously identified 18 urine extracellular vesicle circRNAs that were increased in patients with PCa compared with those with benign prostatic hyperplasia. Spearman correlation analysis of the expression levels of the 18 circRNAs between the tumor tissue and matched urine extracellular vesicles in 30 PCa patients showed that circSCAF8 had the highest R² (R² = 0.635, P < 0.001). The Cox proportional hazards regression model was used to estimate the effect of circSCAF8 on progression-free survival. The in vitro and in vivo functional experiments were implemented to investigate the effects of circSCAF8 on the phenotype of PCa. We found that the knockdown of circSCAF8 in PCa cells suppressed the proliferation, migration, and invasion ability, while overexpression of circSCAF8 had the opposite effects. Similar results were observed in vivo. In a cohort of 85 patients who had undergone radical prostatectomy, circSCAF8 expression in PCa tissues was a powerful predictor of progression-free survival (HR = 2.14, P = 0.022). Mechanistically, circSCAF8 can function by binding to both miR-140-3p and miR-335 to regulate LIF expression and activate the LIF-STAT3 pathway that leads to the growth and metastasis of PCa. Collectively, our findings demonstrate that circSCAF8 contributes to PCa progression through the circSCAF8-miR-140-3p/miR-335-LIF pathway.

A novel 3'tRNA-derived fragment tRF-Val promotes proliferation and inhibits apoptosis by targeting EEF1A1 in gastric cancer

At present, it is commonly believed that tRFs and tiRNAs are formed by the specific and selective shear of tRNAs under certain pressure stimulation, rather than by random degradation of tRNA. tRFs and tiRNAs have been reported to contribute to the biological process of a variety of human cancers. However, the evidence for the mechanisms of tRFs and tiRNAs in the occurrence and development of gastric cancer (GC) is still insufficient. Here, we aimed to explore the carcinogenic roles of tRFs and tiRNAs in GC with RNA-sequencing technique, and found a novel 3'tRNA-derived fragment tRF-Val was significantly upregulated in GC tissues and cell lines. tRF-Val expression was positively correlated with tumor size and the depth of tumor invasion in GC tissues. Functionally, tRF-Val promoted proliferation and invasion, and inhibited apoptosis in GC cells. Mechanistically, tRF-Val directly bound to the chaperone molecule EEF1A1, mediated its transport into the nucleus and promoted its interaction with MDM2 (a specific p53 E3 ubiquitin ligase), thus inhibiting the downstream molecular pathway of p53 and promoting GC progression. These findings provided a new potential therapeutic target for GC and a new explanation for the occurrence of GC.

Diversity of bacterial small RNAs drives competitive strategies for a mutual chaperone

Hundreds of bacterial small RNAs (sRNAs) require the Hfq chaperone to regulate mRNA expression. Hfq is limiting, thus competition among sRNAs for binding to Hfq shapes the proteomes of individual cells. To understand how sRNAs compete for a common partner, we present a single-molecule fluorescence platform to simultaneously visualize binding and release of multiple sRNAs with Hfq. We show that RNA residents rarely dissociate on their own. Instead, clashes between residents and challengers on the same face of Hfq cause rapid exchange, whereas RNAs that recognize different surfaces may cohabit Hfq for several minutes before one RNA departs. The prevalence of these pathways depends on the structure of each RNA and how it interacts with Hfq. We propose that sRNA diversity creates many pairwise interactions with Hfq that allow for distinct biological outcomes: active exchange favors fast regulation whereas co-residence of dissimilar RNAs favors target co-recognition or target exclusion.

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

An atypical class of non-coding small RNAs is produced in rice leaves upon bacterial infection

Non-coding small RNAs (sRNA) act as mediators of gene silencing and regulate plant growth, development and stress responses. Early insights into plant sRNAs established a role in antiviral defense and they are now extensively studied across plant-microbe interactions. Here, sRNA sequencing discovered a class of sRNA in rice (Oryza sativa) specifically associated with foliar diseases caused by Xanthomonas oryzae bacteria. Xanthomonas-induced small RNAs (xisRNAs) loci were distinctively upregulated in response to diverse virulent strains at an early stage of infection producing a single duplex of 20-22 nt sRNAs. xisRNAs production was dependent on the Type III secretion system, a major bacterial virulence factor for host colonization. xisRNA loci overlap with annotated transcripts sequences, with about half of them encoding protein kinase domain proteins. A number of the corresponding rice cis-genes have documented functions in immune signaling and xisRNA loci predominantly coincide with the coding sequence of a conserved kinase motif. xisRNAs exhibit features of small interfering RNAs and their biosynthesis depend on canonical components OsDCL1 and OsHEN1. xisRNA induction possibly mediates post-transcriptional gene silencing but they do not broadly suppress cis-genes expression on the basis of mRNA-seq data. Overall, our results identify a group of unusual sRNAs with a potential role in plant-microbe interactions.

First come, first served: superinfection exclusion in Deformed wing virus is dependent upon sequence identity and not the order of virus acquisition

Deformed wing virus (DWV) is the most important globally distributed pathogen of honey bees and, when vectored by the ectoparasite Varroa destructor, is associated with high levels of colony losses. Divergent DWV types may differ in their pathogenicity and are reported to exhibit superinfection exclusion upon sequential infections, an inevitability in a Varroa-infested colony. We used a reverse genetic approach to investigate competition and interactions between genetically distinct or related virus strains, analysing viral load over time, tissue distribution with reporter gene-expressing viruses and recombination between virus variants. Transient competition occurred irrespective of the order of virus acquisition, indicating no directionality or dominance. Over longer periods, the ability to compete with a pre-existing infection correlated with the genetic divergence of the inoculae. Genetic recombination was observed throughout the DWV genome with recombinants accounting for ~2% of the population as determined by deep sequencing. We propose that superinfection exclusion, if it occurs at all, is a consequence of a cross-reactive RNAi response to the viruses involved, explaining the lack of dominance of one virus type over another. A better understanding of the consequences of dual- and superinfection will inform development of cross-protective honey bee vaccines and landscape-scale DWV transmission and evolution.

Angiogenin mediates paternal inflammation-induced metabolic disorders in offspring through sperm tsRNAs

Paternal environmental inputs can influence various phenotypes in offspring, presenting tremendous implications for basic biology and public health and policy. However, which signals function as a nexus to transmit paternal environmental inputs to offspring remains unclear. Here we show that offspring of fathers with inflammation exhibit metabolic disorders including glucose intolerance and obesity. Deletion of a mouse tRNA RNase, Angiogenin (Ang), abolished paternal inflammation-induced metabolic disorders in offspring. Additionally, Ang deletion prevented the inflammation-induced alteration of 5'-tRNA-derived small RNAs (5'-tsRNAs) expression profile in sperm, which might be essential in composing a sperm RNA 'coding signature' that is needed for paternal epigenetic memory. Microinjection of sperm 30-40 nt RNA fractions (predominantly 5'-tsRNAs) from inflammatory Ang^+/+ males but not Ang^-/- males resulted in metabolic disorders in the resultant offspring. Moreover, zygotic injection with synthetic 5'-tsRNAs which increased in inflammatory mouse sperm and decreased by Ang deletion partially resembled paternal inflammation-induced metabolic disorders in offspring. Together, our findings demonstrate that Ang-mediated biogenesis of 5'-tsRNAs in sperm contributes to paternal inflammation-induced metabolic disorders in offspring.

Loss of SNORA73 reprograms cellular metabolism and protects against steatohepatitis

Dyslipidemia and resulting lipotoxicity are pathologic signatures of metabolic syndrome and type 2 diabetes. Excess lipid causes cell dysfunction and induces cell death through pleiotropic mechanisms that link to oxidative stress. However, pathways that regulate the response to metabolic stress are not well understood. Herein, we show that disruption of the box H/ACA SNORA73 small nucleolar RNAs encoded within the small nucleolar RNA hosting gene 3 (Snhg3) causes resistance to lipid-induced cell death and general oxidative stress in cultured cells. This protection from metabolic stress is associated with broad reprogramming of oxidative metabolism that is dependent on the mammalian target of rapamycin signaling axis. Furthermore, we show that knockdown of SNORA73 in vivo protects against hepatic steatosis and lipid-induced oxidative stress and inflammation. Our findings demonstrate a role for SNORA73 in the regulation of metabolism and lipotoxicity.

RNase κ promotes robust piRNA production by generating 2',3'-cyclic phosphate-containing precursors

In animal germlines, PIWI proteins and the associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons. Here we report the extensive sequence and quantitative correlations between 2',3'-cyclic phosphate-containing RNAs (cP-RNAs), identified using cP-RNA-seq, and piRNAs in the Bombyx germ cell line and mouse testes. The cP-RNAs containing 5'-phosphate (P-cP-RNAs) identified by P-cP-RNA-seq harbor highly consistent 5'-end positions as the piRNAs and are loaded onto PIWI protein, suggesting their direct utilization as piRNA precursors. We identified Bombyx RNase Kappa (BmRNase κ) as a mitochondria-associated endoribonuclease which produces cP-RNAs during piRNA biogenesis. BmRNase κ-depletion elevated transposon levels and disrupted a piRNA-mediated sex determination in Bombyx embryos, indicating the crucial roles of BmRNase κ in piRNA biogenesis and embryonic development. Our results reveal a BmRNase κ-engaged piRNA biogenesis pathway, in which the generation of cP-RNAs promotes robust piRNA production.

A small RNA regulates pprM, a modulator of pleiotropic proteins promoting DNA repair, in Deinococcus radiodurans under ionizing radiation

Networks of transcriptional and post-transcriptional regulators are critical for bacterial survival and adaptation to environmental stressors. While transcriptional regulators provide rapid activation and/or repression of a wide-network of genes, post-transcriptional regulators, such as small RNAs (sRNAs), are also important to fine-tune gene expression. However, the mechanisms of sRNAs remain poorly understood, especially in less-studied bacteria. Deinococcus radiodurans is a gram-positive bacterium resistant to extreme levels of ionizing radiation (IR). Although multiple unique regulatory systems (e.g., the Radiation and Desiccation Response (RDR)) have been identified in this organism, the role of post-transcriptional regulators has not been characterized within the IR response. In this study, we have characterized an sRNA, PprS (formerly Dsr2), as a post-transcriptional coordinator of IR recovery in D. radiodurans. PprS showed differential expression specifically under IR and knockdown of PprS resulted in reduced survival and growth under IR, suggesting its importance in regulating post-radiation recovery. We determined a number of potential RNA targets involved in several pathways including translation and DNA repair. Specifically, we confirmed that PprS binds within the coding region to stabilize the pprM (DR_0907) transcript, a RDR modulator. Overall, these results are the first to present an additional layer of sRNA-based control in DNA repair pathways associated with bacterial radioresistance.

Translation and codon usage regulate Argonaute slicer activity to trigger small RNA biogenesis

In the Caenorhabditis elegans germline, thousands of mRNAs are concomitantly expressed with antisense 22G-RNAs, which are loaded into the Argonaute CSR-1. Despite their essential functions for animal fertility and embryonic development, how CSR-1 22G-RNAs are produced remains unknown. Here, we show that CSR-1 slicer activity is primarily involved in triggering the synthesis of small RNAs on the coding sequences of germline mRNAs and post-transcriptionally regulates a fraction of targets. CSR-1-cleaved mRNAs prime the RNA-dependent RNA polymerase, EGO-1, to synthesize 22G-RNAs in phase with translating ribosomes, in contrast to other 22G-RNAs mostly synthesized in germ granules. Moreover, codon optimality and efficient translation antagonize CSR-1 slicing and 22G-RNAs biogenesis. We propose that codon usage differences encoded into mRNA sequences might be a conserved strategy in eukaryotes to regulate small RNA biogenesis and Argonaute targeting.

Unraveling the role of microRNA/isomiR network in multiple primary melanoma pathogenesis

Malignant cutaneous melanoma (CM) is a potentially lethal form of skin cancer whose worldwide incidence has been constantly increasing over the past decades. During their lifetime, about 8% of CM patients will develop multiple primary melanomas (MPMs), usually at a young age and within 3 years from the first tumor/diagnosis. With the aim of improving our knowledge on MPM biology and pathogenesis, we explored the miRNome of 24 single and multiple primary melanomas, including multiple tumors from the same patient, using a small RNA-sequencing approach. From a supervised analysis, 22 miRNAs were differentially expressed in MPM compared to single CM, including key miRNAs involved in epithelial-mesenchymal transition. The first and second melanoma from the same patient presented a different miRNA profile. Ten miRNAs, including miR-25-3p, 149-5p, 92b-3p, 211-5p, 125a-5p, 125b-5p, 205-5p, 200b-3p, 21-5p, and 146a-5p, were further validated in 47 single and multiple melanoma samples. Pathway enrichment analysis of miRNA target genes revealed a more differentiated and less invasive status of MPMs compared to CMs. Bioinformatic analyses at the miRNA isoform (isomiR) level detected a panel of highly expressed isomiRs belonging to miRNA families implicated in human tumorigenesis, including miR-200, miR-30, and miR-10 family. Moreover, we identified hsa-miR-125a-5p|0|-2 isoform as tenfold over-represented in melanoma than the canonical form and differentially expressed in MPMs arising in the same patient. Target prediction analysis revealed that the miRNA shortening could change the pattern of target gene regulation, specifically in genes implicated in cell adhesion and neuronal differentiation. Overall, we provided a putative and comprehensive characterization of the miRNA/isomiR regulatory network of MPMs, highlighting mechanisms of tumor development and molecular features differentiating this subtype from single melanomas.

RNA binding of Hfq monomers promotes RelA-mediated hexamerization in a limiting Hfq environment

The RNA chaperone Hfq, acting as a hexamer, is a known mediator of post-transcriptional regulation, expediting basepairing between small RNAs (sRNAs) and their target mRNAs. However, the intricate details associated with Hfq-RNA biogenesis are still unclear. Previously, we reported that the stringent response regulator, RelA, is a functional partner of Hfq that facilitates Hfq-mediated sRNA-mRNA regulation in vivo and induces Hfq hexamerization in vitro. Here we show that RelA-mediated Hfq hexamerization requires an initial binding of RNA, preferably sRNA to Hfq monomers. By interacting with a Shine-Dalgarno-like sequence (GGAG) in the sRNA, RelA stabilizes the initially unstable complex of RNA bound-Hfq monomer, enabling the attachment of more Hfq subunits to form a functional hexamer. Overall, our study showing that RNA binding to Hfq monomers is at the heart of RelA-mediated Hfq hexamerization, challenges the previous concept that only Hfq hexamers can bind RNA.

Effects of individual base-pairs on in vivo target search and destruction kinetics of bacterial small RNA

Base-pairing interactions mediate many intermolecular target recognition events. Even a single base-pair mismatch can cause a substantial difference in activity but how such changes influence the target search kinetics in vivo is unknown. Here, we use high-throughput sequencing and quantitative super-resolution imaging to probe the mutants of bacterial small RNA, SgrS, and their regulation of ptsG mRNA target. Mutations that disrupt binding of a chaperone protein, Hfq, and are distal to the mRNA annealing region still decrease the rate of target association, k_on, and increase the dissociation rate, k_off, showing that Hfq directly facilitates sRNA-mRNA annealing in vivo. Single base-pair mismatches in the annealing region reduce k_on by 24-31% and increase k_off by 14-25%, extending the time it takes to find and destroy the target by about a third. The effects of disrupting contiguous base-pairing are much more modest than that expected from thermodynamics, suggesting that Hfq buffers base-pair disruptions.

Type I toxin-antitoxin systems contribute to the maintenance of mobile genetic elements in Clostridioides difficile

Toxin-antitoxin (TA) systems are widespread on mobile genetic elements and in bacterial chromosomes. In type I TA, synthesis of the toxin protein is prevented by the transcription of an antitoxin RNA. The first type I TA were recently identified in the human enteropathogen Clostridioides difficile. Here we report the characterization of five additional type I TA within phiCD630-1 (CD0977.1-RCd11, CD0904.1-RCd13 and CD0956.3-RCd14) and phiCD630-2 (CD2889-RCd12 and CD2907.2-RCd15) prophages of C. difficile strain 630. Toxin genes encode 34 to 47 amino acid peptides and their ectopic expression in C. difficile induces growth arrest that is neutralized by antitoxin RNA co-expression. We show that type I TA located within the phiCD630-1 prophage contribute to its stability and heritability. We have made use of a type I TA toxin gene to generate an efficient mutagenesis tool for this bacterium that allowed investigation of the role of these widespread TA in prophage maintenance.

MicroRNA dynamics during hibernation of the Australian central bearded dragon (Pogona vitticeps)

Hibernation is a physiological state employed by many animals that are exposed to limited food and adverse winter conditions. Controlling tissue-specific and organism wide changes in metabolism and cellular function requires precise regulation of gene expression, including by microRNAs (miRNAs). Here we profile miRNA expression in the central bearded dragon (Pogona vitticeps) using small RNA sequencing of brain, heart, and skeletal muscle from individuals in late hibernation and four days post-arousal. A total of 1295 miRNAs were identified in the central bearded dragon genome; 664 of which were novel to central bearded dragon. We identified differentially expressed miRNAs (DEmiRs) in all tissues and correlated mRNA expression with known and predicted target mRNAs. Functional analysis of DEmiR targets revealed an enrichment of differentially expressed mRNA targets involved in metabolic processes. However, we failed to reveal biologically relevant tissue-specific processes subjected to miRNA-mediated regulation in heart and skeletal muscle. In brain, neuroprotective pathways were identified as potential targets regulated by miRNAs. Our data suggests that miRNAs are necessary for modulating the shift in cellular metabolism during hibernation and regulating neuroprotection in the brain. This study is the first of its kind in a hibernating reptile and provides key insight into this ephemeral phenotype.

Small RNA profiling in Pinus pinaster reveals the transcriptome of developing seeds and highlights differences between zygotic and somatic embryos

Regulation of seed development by small non-coding RNAs (sRNAs) is an important mechanism controlling a crucial phase of the life cycle of seed plants. In this work, sRNAs from seed tissues (zygotic embryos and megagametophytes) and from somatic embryos of Pinus pinaster were analysed to identify putative regulators of seed/embryo development in conifers. In total, sixteen sRNA libraries covering several developmental stages were sequenced. We show that embryos and megagametophytes express a large population of 21-nt sRNAs and that substantial amounts of 24-nt sRNAs were also detected, especially in somatic embryos. A total of 215 conserved miRNAs, one third of which are conifer-specific, and 212 high-confidence novel miRNAs were annotated. MIR159, MIR171 and MIR394 families were found in embryos, but were greatly reduced in megagametophytes. Other families, like MIR397 and MIR408, predominated in somatic embryos and megagametophytes, suggesting their expression in somatic embryos is associated with in vitro conditions. Analysis of the predicted miRNA targets suggests that miRNA functions are relevant in several processes including transporter activity at the cotyledon-forming stage, and sulfur metabolism across several developmental stages. An important resource for studying conifer embryogenesis is made available here, which may also provide insightful clues for improving clonal propagation via somatic embryogenesis.

Modification of messenger RNA by 2'-O-methylation regulates gene expression in vivo

Epitranscriptomic modifications of mRNA are important regulators of gene expression. While internal 2'-O-methylation (Nm) has been discovered on mRNA, questions remain about its origin and function in cells and organisms. Here, we show that internal Nm modification can be guided by small nucleolar RNAs (snoRNAs), and that these Nm sites can regulate mRNA and protein expression. Specifically, two box C/D snoRNAs (SNORDs) and the 2'-O-methyltransferase fibrillarin lead to Nm modification in the protein-coding region of peroxidasin (Pxdn). The presence of Nm modification increases Pxdn mRNA expression but inhibits its translation, regulating PXDN protein expression and enzyme activity both in vitro and in vivo. Our findings support a model in which snoRNA-guided Nm modifications of mRNA can regulate physiologic gene expression by altering mRNA levels and tuning protein translation.

Complement C3 activation regulates the production of tRNA-derived fragments Gly-tRFs and promotes alcohol-induced liver injury and steatosis

Complement is known to play a role in alcoholic fatty liver disease (AFLD), but the underlying mechanisms are poorly understood, thereby constraining the development of a rational approach for therapeutic intervention in the complement system. C3 deficiency has been shown to impart protective effects against ethanol-induced hepatic steatosis and inflammation. Here we demonstrate a protection effect in wild-type mice by treatment with CR2-Crry, a specific inhibitor of C3 activation. The expression of glycine transfer (t) RNA-derived fragments (Gly-tRFs) is upregulated in ethanol-fed mice and inhibition of Gly-tRFs in vivo decreases chronic ethanol feeding-induced hepatosteatosis without affecting inflammation. The expression of Gly-tRF was downregulated in C3-deficient or CR2-Crry-treated mice, but not in C5-deficient mice; Gly-tRF expression was restored by the C3 activation products C3a or Asp (C3a-des-Arg) via the regulation of CYP2E1. Transcriptome profiling of hepatic tissues showed that Gly-tRF inhibitors upregulate the expression of sirtuin1 (Sirt1) and subsequently affect downstream lipogenesis and β-oxidation pathways. Mechanistically, Gly-tRF interacts with AGO3 to downregulate Sirt1 expression via sequence complementarity in the 3' UTR. Notably, the expression levels of C3d, CYP2E1 and Gly-tRF are upregulated, whereas Sirt1 is decreased in AFLD patients compared to healthy controls. Collectively, our findings suggest that C3 activation products contribute to hepatosteatosis by regulating the expression of Gly-tRF. Complement inhibition at the C3 activation step and treatment with Gly-tRF inhibitors may be potential and precise therapeutic approaches for AFLD.

The 3' end formation in small RNAs

Small RNAs are a major class of RNAs along with transfer RNAs, ribosomal RNAs, and messenger RNAs. They vary in size from less than 100 nucleotides to several thousand nucleotides and have been identified and characterized both in prokaryotes and eukaryotes. Small RNAs participate in a variety of cellular functions including regulating RNA synthesis, RNA processing, guiding modifications in RNA, and in transport of proteins. Small RNAs are generated by a series of posttranscriptional processing steps following transcription. While RNA 5' end structure, 5' cap formation, and RNA processing mechanisms have been fairly well characterized, the 3' end processing is poorly understood. Recent data point to an emerging theme in small RNAs metabolism in which the 3' end processing is mediated by the exosome, a large multienzyme complex. In addition to removal of nucleotides by the exosome, there is simultaneous rebuilding of the 3' end of some small RNA by adenylation and/or uridylation. This review presents a picture of both degradative and rebuilding reactions operative on the 3' end of some small RNA molecules in prokaryotes and eukaryotes.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.